| 查看: 264 | 回复: 0 | |||
dage16wo木虫 (小有名气)
|
[求助]
fmincon优化的问题
|
|
function f=objtctfun(x) global G1 L1 G4 L4 Gg Gd Gw G05 L05 G04 L04 G03 L03 G02 L02 G01 L01 k l a c f=(G1*L1+G4*L4+(Gg+Gd)*L1+Gw*L4+G05*L05+G04*L04+G03*L03+G02*L02+G01*L01)*cos(k*pi/180)*(l+x(3))/(sqrt(a^2*c^2+x(2)^2*x(1)^2-2*a*c*x(1)*x(2))); 以上为目标函数 function [c,ceq]=conobject(x) global a c1 k l h1 k1 k2 pi c=[(l+x(3))-sqrt(a^2+x(1)^2)-sqrt(x(2)^2+c1^2); sqrt(x(2)^2+c1^2)-sqrt(a^2+x(1)^2)-(l+x(3)); sqrt(x(2)^2+c1^2)-sqrt(a^2+x(1)^2)-(l+x(3)); (l+2*x(3))-(sqrt(a^2+x(1)^2)+sqrt(x(2)^2+c1^2)); (sqrt(x(2)^2+c1^2))-(sqrt(a^2+x(1)^2)+(l+2*x(3))); sqrt(x(2)^2+c1^2)-(sqrt(a^2+x(1)^2)+(l+2*x(3))); x(1)-h1; c1-x(1); x(2)-k2; k1-x(2); a^2+x(1)^2+c1^2+x(2)^2-2*(a*x(2)+x(1)*c1)*cos(k*pi/180)+2*sqrt(a^2*c1^2+x(2)^2*x(1)^2-2*a*c1*x(1)*x(2))*sin(k*pi/180) -(l+2*x(3))^2; ]; ceq=a^2+x(1)^2+c1^2+x(2)^2-2*(a*x(2)+x(1)*c1)-(l+x(3))^2; 以上为约束函数 G1=410 L1=21550 G4=200 L4=20730 Gw=49.79 Gg=199.16 Gd=100 G05=0 L05=0 G04=0 L04=0 G03=830.72 L03=17707 G02=1214.5 L02=11517 G01= 1659.95 L01=4445 a=100 c1=390 k=30 l=500 h1=4000 k1=1000 k2=2500 初始值x0=[500:500:500] 非线性优化为题 出现问题 ??? Error using ==> vertcat CAT arguments dimensions are not consistent. Error in ==> conobject at 3 c=[(l+x(3))-sqrt(a^2+x(1)^2)-sqrt(x(2)^2+c1^2); Error in ==> fmincon at 654 [ctmp,ceqtmp] = feval(confcn{3},X,varargin{:}); Caused by: Failure in initial user-supplied nonlinear constraint function evaluation. FMINCON cannot continue. |
回复此楼» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
职称评审没过，求安慰
已经有22人回复
-
垃圾破二本职称评审标准
已经有15人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
