应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求指点，1stopt如何是所求值为正值
51/1返回列表
查看: 893  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

houchuyue

银虫 (初入文坛)

[求助] 求指点，1stopt如何是所求值为正值 已有2人参与

1.我的目标是求X1,X2,X3,且是正数的解，如何修改，应该如何添加约束条件。2。为何每次运行，结果都不同，如何处理？
Constant a=0.01195,b=0.0256,k=0.0105,B3=0.0027;
Parameters x1, x2, x3;
Function
a+b=x1+x2;
x3+2*b=x1+2*x2;
x2*x3*y/x1=k;
I=0.5*(x1+4*x2+x3+4*b);
f=-2*(1.303*ln(1+1.2*I^(0.5))+1.303*I*(1/(1+1.2*I^(0.5))*0.5*1.2*I^(-0.5)));
B1=0.2103+(2*0.4711/(1.4^(2)*I))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))));
B2=-2*0.4711/(1.4^(2)*I^(2))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))))-2*0.4711/(1.4^(2)*I)*(0.5*1.4*I^(-0.5)*exp(-1.4*I^(0.5))+(1+1.4*I^(0.5))*exp(-1.4*I^(0.5))*(-0.5*1.4*I^(-0.5)));
B4=0.2+(2*2.7/(1.4^(2)*I))*(1-(1+1.4*I^(0.5))*exp(-1.4*I^(0.5)))+(2*(-30.7)/(12^(2)*I))*(1-(1+12*I^(0.5))*exp(-12*I^(0.5)));
B5=-2*2.7/(1.4^(2)*I^(2))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))))-2*2.7/(1.4^(2)*I)*((0.5*1.4*I^(-0.5))*exp(-1.4*I^(0.5))+(1+1.4*I^(0.5))*exp(-1.4*I^(0.5))*(-0.5*1.4*I^(-0.5)))-2*(-30.7)/(12^(2)*I^(2))*(1-(1+12*I^(0.5)*exp(-12*I^(0.5))))-2*(-30.7)/(12^(2)*I)*((0.5*12*I^(-0.5))*exp(-12*I^(0.5))+(1+12*I^(0.5))*exp(-12*I^(0.5))*(-0.5*12*I^(-0.5)));
B6=0.375+(2*1.5/(1.4^(2)*I))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))));
B7=-2*1.5/(1.4^(2)*I^(2))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))))-2*1.5/(1.4^(2)*I)*(0.5*1.4*I^(-0.5)*exp(-1.4*I^(0.5))+(1+1.4*I^(0.5))*exp(-1.4*I^(0.5))*(-0.5*1.4*I^(-0.5)));
ln(Y)=f+2*(x1-x3)*B1+4*x3*x1*B2+2*(x2+x3)*B3+2*((x2+x3)*(x3+b)+x3*x2)*0.0416+2*b*(2*x2*B4+B5+2*x1*B6-B7);

迭代数: 72
计算用时(时:分:秒:毫秒): 00:00:31:104
计算中止原因: 达到收敛判定标准
优化算法: 准牛顿法(BFGS) + 通用全局优化法
函数表达式 1: 0.01195+0.0256-(x1+x2)
         2: x3+2*0.0256-(x1+2*x2)
         3: x2*x3*y/x1-(0.0105)
         4: i-(0.5*(x1+4*x2+x3+4*0.0256))
         5: f-(-2*(1.303*ln(1+1.2*i^(0.5))+1.303*i*(1/(1+1.2*i^(0.5))*0.5*1.2*i^(-0.5))))
         6: b1-(0.2103+(2*0.4711/(1.4^(2)*i))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5)))))
         7: b2-(-2*0.4711/(1.4^(2)*i^(2))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5))))-2*0.4711/(1.4^(2)*i)*(0.5*1.4*i^(
            -0.5)*exp(-1.4*i^(0.5))+(1+1.4*i^(0.5))*exp(-1.4*i^(0.5))*(-0.5*1.4*i^(-0.5))))
         8: b4-(0.2+(2*2.7/(1.4^(2)*i))*(1-(1+1.4*i^(0.5))*exp(-1.4*i^(0.5)))+(2*(-30.7)/(12^(2)*i))*(1-(1+12*i^(0.5))
            *exp(-12*i^(0.5))))
         9: b5-(-2*2.7/(1.4^(2)*i^(2))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5))))-2*2.7/(1.4^(2)*i)*((0.5*1.4*i^(-0.5))
            *exp(-1.4*i^(0.5))+(1+1.4*i^(0.5))*exp(-1.4*i^(0.5))*(-0.5*1.4*i^(-0.5)))-2*(-30.7)/(12^(2)*i^(2))*(1
            -(1+12*i^(0.5)*exp(-12*i^(0.5))))-2*(-30.7)/(12^(2)*i)*((0.5*12*i^(-0.5))*exp(-12*i^(0.5))+(1+12*i^(0.5))
            *exp(-12*i^(0.5))*(-0.5*12*i^(-0.5))))
         10: b6-(0.375+(2*1.5/(1.4^(2)*i))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5)))))
         11: b7-(-2*1.5/(1.4^(2)*i^(2))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5))))-2*1.5/(1.4^(2)*i)*(0.5*1.4*i^(-0.5)
            *exp(-1.4*i^(0.5))+(1+1.4*i^(0.5))*exp(-1.4*i^(0.5))*(-0.5*1.4*i^(-0.5))))
         12: ln(y)-(f+2*(x1-x3)*b1+4*x3*x1*b2+2*(x2+x3)*0.0027+2*((x2+x3)*(x3+0.0256)+x3*x2)*0.0416+2*0.0256*(2*x2
            *b4+b5+2*x1*b6-b7))
目标函数值: 0.0346057342149256
x1: -0.613628828746281
x2: 0.651181234844569
x3: 0.637530987966978
y: 0.0356211237987828
i: 1.36550919200013
f: -3.04451434010513
b1: 0.0981303646848507
b2: 0.149337201828133
b4: 0.869550118362684
b5: 0.855878503182631
b6: 0.0178474428783668
b7: 0.475498544863951

====== 计算结束 ======

[ Last edited by houchuyue on 2014-3-12 at 16:18 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-12 16:11:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
houchuyue: 金币+2, ★★★很有帮助 2014-03-12 17:08:20
不加大于0的范围约束求得的x1、x2、x3也都是大于0的啊，且稳定唯一：

目标函数值(最小): 2.44738091577421E-22
x1: 0.0174024208040297
x2: 0.0201476131995002
x3: 0.00649756064023006
y: 1.39581797264244
i: 0.103445344352546
f: -1.21342501037417
b1: -1.12354414258258
b2: 15.7971962194402
b4: -1.4883351519721
b5: 81.0412352437473
b6: -3.87200958150146
b7: 50.2988629360697
赞一下 回复此楼
3楼2014-03-12 16:46:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
houchuyue: 金币+5, ★★★很有帮助 2014-03-12 16:33:27
Parameters x1[0,], x2[0,], x3[0,];
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-03-12 16:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houchuyue

银虫 (初入文坛)
引用回帖:
3楼: Originally posted by dingd at 2014-03-12 16:46:46
不加大于0的范围约束求得的x1、x2、x3也都是大于0的啊，且稳定唯一：

目标函数值(最小): 2.44738091577421E-22
x1: 0.0174024208040297
x2: 0.0201476131995002
x3: 0.00649756064023006
y: 1.3958179726424 ...

我用的是1.5破解版。算不出满意的结果，一直变动。请问是版本原因？还是我需要在算法设置上要进行调整？
赞一下 回复此楼
4楼2014-03-12 17:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

1.5还是破解版，原因就说不清了。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

houchuyue
5楼2014-03-12 17:11:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭