本人利用matlab里的fmincon函数求解一个优化问题,主函数:
clc
clear all
lb=[3 0 0 1 0 0 1 0 0 1 0 0 1 0 0 1];
ub=[3 5 5 4 5 5 4 5 5 4 5 5 4 5 5 4];
x0=[3 1 1 2 1 1 2 1 1 2 1 1 2 1 1 2];
options=optimset('display','off','Algorithm','sqp');
[x,fval]=fmincon(@fitness,x0,[],[],[],[],lb,ub,@cony,options);
运行错误提示:
未定义函数或变量 "G"。
出错 extraction (line 35)
q=[G,V3,F,f1,f2];
出错 fitness (line 9)
q=extraction(wx);
出错 finDiffEvalAndChkErr (line 25)
出错 finitedifferences (line 127)
出错 computeFinDiffGradAndJac (line 23)
出错 sqpLineSearch (line 120)
出错 fmincon (line 859)
[X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = sqpLineSearch(funfcn,X,full(A),full(B),full(Aeq),full(Beq), ...
出错 SQP (line 9)
[x,fval]=fmincon(@fitness,x0,[],[],[],[],lb,ub,@cony,options);
fitness函数:
function f = fitness(wx)
a11=900; %浸膏价格G
b11=1100;%挥发油价格V3
c11=120; %原料价格
d11=10; %蒸汽价格
W=120;
q=extraction(wx);
G=q(1);
V3=q(2);
F=q(3);
t=wx(4)+wx(7)+wx(10);
f=(a11*G+b11*V3-d11*F-c11*W)/t;
f=-f;
end
问题就出在q=extraction(wx),这个函数是调用机理模型求解q,难道目标函数fitness里不能调用其他函数吗?并且我确定extraction()函数没有问题,并且程序也能求出q这个值。
cony函数:
function [c,ceq]=cony(x)
q=extraction(x);
f1=q(4);
f2=q(5);
ceq=[];
c=[0.8-f1;f1-1;0.8-f2;f2-1];
end
这个函数就是定义非线性约束的,我也调用了extraction()函数,应该没有问题吧
很着急,谢谢大家了
-5-31 12:27 上传
点击文件名下载附件
13.22 KB, 下载次数: 38
所用程序