| 查看: 1716 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
1stopt 计算问题
|
|||
|
命令如下： Title "Type your title here"; constant k=1863; constant m1=14.86; constant m2=66.80; constant m3=1.46; constant m4=17.66; Parameters x(1:20)[0,1]; Function 0=exp((13836+25.44*k)/8.314/k)*x1^3*x2-x5; 0=exp((35288+204.59*k)/8.314/k)*x1^12*x2^7-x6; 0=exp((18946+25.44*k)/8.314/k)*x1*x2-x7; 0=exp((18566+23.68*k)/8.314/k)*x1*x2^2-x8; 0=exp((-4752.67-75.14*k)/8.314/k)*x1*x2^6-x9; 0=exp((107800+13.44*k)/8.314/k)*x1^3*x3^2-x10; 0=exp((122731+8.85*k)/8.314/k)*x1^2*x3-x11; 0=exp((302587+66.24*k)/8.314/k)*x1^3*x3^2-x12; 0=exp((86231+0.995*k)/8.314/k)*x1*x3-x13; 0=exp((-19512+28.19*k)/8.314/k)*x2^3*x3^2-x14; 0=exp((354932-417.63*k)/8.314/k)*x1^2*x2*x3-x15; 0=exp((1252543-277.06*k)/8.314/k)*x1*x2*x3^3-x16; 0=exp((203172+21.38*k)/8.314/k)*x1^3*x4^2-x17; 0=exp((75158+21.38*k)/8.314/k)*x1*x4-x18; 0=exp((313417-282.66*k)/8.314/k)*x1*x2*x3-x19; 0=exp((2279-13.54*k)/8.314/k)*x2*x4-x20; 0=(0.5*x1+3*x5+12*x6+x7+x8+x9+3*x10+2*x11+3*x12+x13+2*x15+x16+3*x17+x18+x19)*(m1/40+m2/102+m3/64+m4/84)-m1/40; 0=(x2+x5+7*x6+x7+2*x8+6*x9+3*x14+x15+x16+x20)*(m1/40+m2/102+m3/64+m4/84)-m2/102; 0=(x3+x10+x11+2*x12+x13+2*x14+x15+2*x16+x19)*(m1/40+m2/102+m3/64+m4/84)-m3/64; 0=(x4+2*x17+x18+x19+x20)*(m1/40+m2/102+m3/64+m4/84)-m4/84; 0=(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20)-1; 这个程序计算出的结果： 迭代数: 10740 计算用时(时:分:秒:毫秒): 00:00:22:238 计算中止原因: 达到收敛判定标准 优化算法: 标准简面体爬山法 + 通用全局优化法 目标函数值: 0.337300657640182 x1: 0.00280473183833806 x2: 0.389556948473223 x3: 2.47955133843895E-9 x4: 0.028798627469645 x5: 4.47794033080389E-7 x6: 7.19960278061123E-12 x7: 0.079174458763718 x8: 0.0243537011949764 x9: 7.21467656804983E-12 x10: 0.0181130086301488 x11: 2.49572142641952E-10 x12: 3.93442470085121E-11 x13: 2.20909868564788E-9 x14: 2.29566803469628E-11 x15: 4.50226794231978E-11 x16: 7.36016488772663E-9 x17: 0.000119112732787977 x18: 0.13534118437126 x19: 7.8257263248991E-12 x20: 0.00255012081464932 请问其中的目标函数值是什么意思？ 另外 为什么计算出的数值相加不符合第21个方程等于1 的条件？ 能否用其他方法进一步计算其数值？ 能否用其他方法估算计算的初值使其可以应用于matlab或其他需要赋初值的软件中？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有82人回复
-
感觉是下周放榜了
已经有8人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt 拟合计算
已经有6人回复
- 毕不了业了，发帖求助Origin8.0多元非线性拟合的问题
已经有9人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 1stOPT1.5版每次仿真结果不一样，求助！急！
已经有3人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 动力学方程参数估计方法
已经有14人回复
- 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- matlab数据拟合问题求助
已经有7人回复
- origin或1stopt计算积分面积
已经有4人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- Gulp 拟合势函数 初值设置的问题
已经有12人回复
- origin进行非线性拟合如何选择函数或者自定义合适的函数
已经有19人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 为什么同样的数据，origin 和excel做成的拟合公式不同？
已经有13人回复
- 请教下已知曲线以及曲线上的各个散点如何求出这条曲线的方程表达式
已经有10人回复
- 用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手
已经有7人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2013-09-03 20:59:38
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-09-03 18:41:55
3楼2013-09-03 20:56:59
5楼2013-09-04 08:26:23
