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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求指点，1stopt如何是所求值为正值
51/1返回列表
查看: 939  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +6 z1z2z3879 2026-03-16 8/400 2026-03-17 14:03 by 勇敢太监王公公
[考研] 332求调剂 +5 Zz版 2026-03-13 5/250 2026-03-17 11:01 by JineShine
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 290求调剂 +6 @将就将就看 2026-03-10 11/550 2026-03-17 10:22 by @将就将就看
[文学芳草园] 伙伴们，祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 283求调剂，材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
关闭
请填处理意见
关闭