【求助】论文上的程序怎么是错的请高手指导下 3Q - MATLAB爱好者论坛-LabFans.com

wr171903566 · 2007-06-30, 02:36

function f=wolunfun(x)
f=x(1)*(x(2)+x(3)*20)/2
function[c,ceq]=woluncon(x)
c(1)=15-atand(x(3)/x(2));
c(2)=atand(x(3)/x(2))-30;
c(3)=1.949*10^6*exp(-1.7762*2*x(2)/(x(2)+20*x(3)))*(1/(x(1)^3*(x(2)+20*x(3))^3))^0.5-244;
c(4)=114800/(x(1)^3*x(3)*(0.5+(x(2)+1)^0.5))-45.08;
c(5)=((118552/(x(1)*x(2)))^2+(34519/(x(1)*x(3)))^2)^0.5*(1.179*10^(-8)*x(1)^3*x(3)^3)/(x(1)*x(2)-2*x(1)-0.4)^4-0.001*x(1)*x(2);
ceq(1)=x(1)-6.3;
ceq(3)=x(3)-2;
x0=[8,10,2];
lb=[2,7,1];
ub=[10,25,3];
options=optimset('Display','iter','LargeScale','off');;
[x,fval,exitflag,output,lambda]=fmincon(@wolunfun,x0,[],[],[],[],lb,ub,@woluncon,options)

2007-06-30, 02:36	#1
wr171903566 初级会员注册日期: 2007-06-30 帖子: 1 声望力: 0	【求助】论文上的程序怎么是错的请高手指导下 3Q function f=wolunfun(x) f=x(1)(x(2)+x(3)20)/2 function[c,ceq]=woluncon(x) c(1)=15-atand(x(3)/x(2)); c(2)=atand(x(3)/x(2))-30; c(3)=1.94910^6exp(-1.77622x(2)/(x(2)+20x(3)))(1/(x(1)^3(x(2)+20x(3))^3))^0.5-244; c(4)=114800/(x(1)^3x(3)(0.5+(x(2)+1)^0.5))-45.08; c(5)=((118552/(x(1)x(2)))^2+(34519/(x(1)x(3)))^2)^0.5(1.17910^(-8)x(1)^3x(3)^3)/(x(1)x(2)-2x(1)-0.4)^4-0.001x(1)x(2); ceq(1)=x(1)-6.3; ceq(3)=x(3)-2; x0=[8,10,2]; lb=[2,7,1]; ub=[10,25,3]; options=optimset('Display','iter','LargeScale','off');; [x,fval,exitflag,output,lambda]=fmincon(@wolunfun,x0,[],[],[],[],lb,ub,@woluncon,options)