| 查看: 1009 | 回复: 5 | |||
[求助]
求指点，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 ] |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有282人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有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人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68141.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.6小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-03-12 16:25:04
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
houchuyue: 金币+2, ★★★很有帮助 2014-03-12 17:08:20
感谢参与，应助指数 +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
4楼2014-03-12 17:07:29
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
houchuyue5楼2014-03-12 17:11:57
6楼2014-03-12 17:13:15
精华III:
20