24小时热门版块排行榜

返回列表

nmsz

新虫 (正式写手)

应助: 1 (幼儿园)
金币: 47.6
散金: 673
红花: 3
帖子: 741
在线: 283小时
虫号: 964979
注册: 2010-03-08
专业: 数量经济学

[求助] 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)

回复此楼

» 猜你喜欢

求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有177人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

1楼2015-03-27 00:24:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

nmsz

新虫 (正式写手)

应助: 1 (幼儿园)
金币: 47.6
散金: 673
红花: 3
帖子: 741
在线: 283小时
虫号: 964979
注册: 2010-03-08
专业: 数量经济学

各位大神帮忙看看呀！最初的程序是下面的（因为效果不理想改为上面程序）：
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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 nmsz 的主题更新

返回列表