nynw.net
当前位置:首页 >> 一个关于用mAtlAB求函数最优解的问题 >>

一个关于用mAtlAB求函数最优解的问题

你写的f(x,y)式子很矛盾,不应该有min。设题目是去除min后,求f(x,y)的最小值。 xi = [0 40 20 0 20 40 0 40];yi = [30 30 0 0 60 0 60 60];li = [16.9 25.4 27.8 31.9 35 36.6 38.3 42];f = @(x,y)sum(abs(sqrt((x-xi).^2+(y-yi).^2)-li));fmin...

function f=myfun(x) f=x(1)+x(2)+81; end %%函数,其中x=x(1)y=x(2) function [c,ceq]=mycon(x) c(1)=100000-x(1)*x(2); c(2)=1-x(1); c(3)=1-x(2); ceq=[]; end %%条件 clc,clear all x0=[400 ;400]; [x,fval]=fmincon(@myfun,x0,[],[],[],[...

function f=myfun(x) f=x(1)+x(2)+81; end %%函数,其中x=x(1)y=x(2) function [c,ceq]=mycon(x) c(1)=100000-x(1)*x(2); c(2)=1-x(1); c(3)=1-x(2); ceq=[]; end %%条件 clc,clear all x0=[400 ;400]; [x,fval]=fmincon(@myfun,x0,[],[],[],[...

这种题一般只要写出它的目标函数,以及约束条件就可以通过MATLAB求解。 设购买甲x(1)Kg, 购买乙x(2)Kg 这里的目标函数很明显是:费用最低 min 5*x(1)+4*x(2) 约束条件: (1)A元素大于等于12: 2*x(1)+4*x(2)>=12 (2)B元素大于等于36: 2*x(1)+x(2...

参考代码(里面有些系数请再仔细检查,我虽已尽量小心,但不敢保证录入无误): z = @(x)(1.8-x(1))/400*x(2)^2 + (-1.8*x(1)^2+146.24*x(1)-1270)/50*x(2) ... +(-3.24*x(1)^3+262.872*x(1)^2-3771*x(1)-82375)/100;x = fmincon(@(x)-z(x),[0 0...

Wij数据提供不完整,无法求f(x)最优解。 求f(x)最优解可以用matlab的fmincon最小值函数或ga遗传算法函数等来求解。

看楼主的意思,每次求得的解应该是相同的大小的一维数组,因为不知道具体的解的维度,大概写一下,应该像下面这样吧: for i=1:20 x(:,i)=fsolve(@energe,x0,[],Qxi(i),Qyi(i),cx(i),cy(i)); end 注意要把i=i+1去掉,for循环会自动加的。。。 ...

function main() clc; clear all; close all; options = optimset('Algorithm', 'interior-point', 'Display', 'off'); A = []; b = []; Aeq = []; beq = []; lb = [1; -inf]; ub = [inf; inf]; x0 = [1 2]; [x,fval,exitflag] = fmincon(@net_f...

您好:unction main() clc; clear all; close all; options = optimset('Algorithm', 'interior-point', 'Display', 'off'); A = []; b = []; Aeq = []; beq = []; lb = [1; -inf]; ub = [inf; inf]; x0 = [1 2]; [x,fval,exitflag] = fmincon(@...

function f=myfun(x) f=x(1)+x(2)+81; end %%函数,其中x=x(1)y=x(2) function [c,ceq]=mycon(x) c(1)=100000-x(1)*x(2); c(2)=1-x(1); c(3)=1-x(2); ceq=[]; end %%条件 clc,clear all x0=[400 ;400]; [x,fval]=fmincon(@myfun,x0,[],[],[],[...

网站首页 | 网站地图
All rights reserved Powered by www.nynw.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com