| 查看: 2235 | 回复: 1 | |||
[求助]
MATLAB非线性规划结果出现exitflag=5和No active inequalities.急！急！
|
|
源程序如下：（本来是等式约束，效果不理想，后改为不等式约束），请各位高手指教！！！ fun=' 0.6*abs((exp(x(4)) *77467+exp(x(5))*59866+exp(x(6))*137334)/100000-exp(4.338)* exp(x(13))/100000)+ 0.1*abs (exp(x(7))/10-exp(x(8))* exp(4.9352)/10000)+ 0.1*abs((exp(x(9))*77467+exp(x(10))*59866+exp(0.65196)*137334)/10000-exp(10.2418)* exp(2.71)/10000)+ 0.2*abs((exp(x(11)) *77467+exp(x(12)) *59866+exp(1.6212)*137334)/100000-exp(14.043)* exp(0.4112)/100000) '; x0=[0.6428,-2.389, -2.9392,3.0302 , 2.8285,3.1759, 6.5327,8.517,0.1495,0.1,2.4425 ,1.5792 ,11.238]; A=[ 0.2843 -0.1147 0 1 0 0 0 0 0 0 0 0 0 0.3662 0 0 0 1 0 0 0 0 0 0 0 0 -0.2069 0 0 0 0 1 0 0 0 0 0 0 0 -0.0829 0 0 0 0 0 0 0 0 0 0 0 1 0 0 -0.0989 0 0 0 1 0 0 0 0 0 0 0 -0.05405 0 0 0 0 0 1 0 0 0 0 0 0 0.3834 -0.1281 0 0 0 0 0 1 0 0 0 0 0 0.617 0 0 0 0 0 0 0 1 0 0 0 0 0 0.169 0 0 0 0 0 0 0 1 0 0 0 0 0.0918 0 0 0 0 0 0 0 0 1 0 ]; %线性等式约束 b =[3.5 3.07 3.06 11.18 6.87 8.68 -0.39 -1.33 1.99 1.35]; Aeq=[ ]; beq=[ ]; lb=[0.64,-2.6, -3.2,2.7302 , 2.7,3.0, 6.5,8.5,0.12,0.05,2.2 ,1.5 ,11.23]; ub=[0.68,-2.1, -2.6,3.2302 , 3.2,3.3, 6.6,8.6,0.16,0.2,2.6 ,1.7 ,11.26 ]; [x,fval,exitflag,output]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub) |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有291人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
|
各位大神帮忙看看呀！最初的程序是下面的（因为效果不理想改为上面程序）： fun=' abs((exp(x(5)) *77467+exp(x(6))*59866+exp(x(7))*137334)/100000-exp(4.338)* exp(x(14))/100000)+ abs (exp(x(8))/10-exp(x(9))* exp(4.9352)/10000)+ abs((exp(x(10))*77467+exp(x(11))*59866+exp(0.65196)*137334)/10000-exp(10.2418)* exp(2.71)/10000)+ abs((exp(x(12)) *77467+exp(x(13)) *59866+exp(1.6212)*137334)/100000-exp(14.043)* exp(0.4112)/100000) '; x0=[0.6428,-2.389,-2.617,-2.9392,3.0302 , 2.8285,3.1759, 6.5327,8.517,0.1495,0.1,2.4425 ,1.5792 ,11.238, 10.26]; Aeq =[ 0.2843 -0.1147 0 0 1 0 0 0 0 0 0 0 0 0 0 0.3662 0 0 0 0 1 0 0 0 0 0 0 0 0 0 -0.2069 0 0 0 0 0 1 0 0 0 0 0 0 0 0 -0.0829 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 -0.0989 0 0 0 1 0 0 0 0 0 0 0 0 -0.05405 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0.3834 0 -0.1281 0 0 0 0 0 1 0 0 0 0 0 0 0.617 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0.05045 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0.169 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0.0918 0 0 0 0 0 0 0 0 1 0 0 ]; %线性等式约束 beq =[3.487 3.06386 3.04288 11.17834 6.858 8.64885 -0.3822 -1.3254 -10.12 1.975 1.3404]; A=[ ]; b=[ ]; lb=[ ]; ub=[ ]; [x,fval,exitflag,output]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub) |
2楼2015-03-27 00:28:05
5