版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

houchuyue

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 248
帖子: 20
在线: 12.6小时
虫号: 2690455
注册: 2013-09-29
专业: 冶金物理化学与冶金原理

[求助] 求指点，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 ]

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有217人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助1stOpt作图已经有6人回复
求助1stopt含约束非线性拟合已经有5人回复
1stopt求最小值已经有37人回复
谁用1stopt这个软件啊，离散数据拟合方程，方程未知已经有5人回复
求大神帮忙用1stOpt 拟合复数数据已经有13人回复
求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙已经有6人回复
求高版本1stopt拟合，已经有11人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有9人回复
求大侠帮助跑一下1stopt 已经有12人回复
哪位朋友有注册版1stOpt软件？急需已经有5人回复
1stOPT1.5版每次仿真结果不一样，求助！急！已经有3人回复
运算过程有复数出现，求高版本1stOpt 已经有6人回复
1stopt运行程序，总是出现字符函数错误的提示，请高手帮助已经有5人回复
请问matlab拟合二元方程的问题已经有3人回复
麻烦高人指点一下1stopt和spss的非线性回归中weibull的两参数的公式应该怎么表达啊？已经有9人回复
MATLAB拟合工具如何引入误差棒？已经有11人回复
学1stopt要多久？已经有12人回复
用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手已经有7人回复

1楼 2014-03-12 16:11:48

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68121.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★
感谢参与，应助指数 +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的回帖

houchuyue

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 248
帖子: 20
在线: 12.6小时
虫号: 2690455
注册: 2013-09-29
专业: 冶金物理化学与冶金原理

引用回帖:

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的回帖