| 查看: 2591 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
求助用matlab解多元指数型方程组 已有1人参与
|
||||
|
方程如图。直接用matlab解方程组解不出。 参考 http://wenku.baidu.com/link?url= ... jMX4Spsljya0IuANYga 这个帖子的方法用牛顿法迭代求解，带入初值410,450,30迭代得不到结果。 方程近似解x1=403 x2=446 n=30 有大神能指点一下么  IMG_20150409_230048.jpg |
回复此楼» 猜你喜欢
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有22人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
物理学I论文润色/翻译怎么收费?
已经有297人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有17人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有11人回复
11楼2015-04-12 21:52:42
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-04-10 08:49:55
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
|
抱歉，修正2楼的说法，我在解方程时候，把0.002写成了0.0002，更正后结果如下： 1stopt代码： Parameters x1,x2,n; //Variable ; Function x1/207000*3/7*(x1/433)^(n-1)-0.0001=0; x1/207000*3/7*(x2/433)^(n-1)-0.002=0; Ln(20)/Ln(x2/x1)-n=0; 计算结果： Algorithms: Universal Global Optimation(UGO1) Function 1: x1/207000*3/7*(x1/433)^(n-1)-0.0001-(0) = 9.341412192E-6 2: x1/207000*3/7*(x2/433)^(n-1)-0.002-(0) = -2.610202075E-7 3: ln(20)/ln(x2/x1)-n-(0) = -1.573571353E-7 Objective Function (Min.): 8.72718617816277E-11 x1: 406.638344029709 x2: 444.682149647678 n: 33.4959878584276 |
3楼2015-04-10 09:11:40
|
感谢你的回答。 我用1stopt软件计算了一遍结果如下 ====== 结果 ====== 迭代数: 106 计算用时(时:分:秒:毫秒): 00:00:01:04 计算中止原因: 达到收敛判定标准 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 函数表达式 1: x1/207000*3/7*(x1/433)^(n-1)-0.0001-(0) 2: x1/207000*3/7*(x2/433)^(n-1)-0.002-(0) 3: ln(20)/ln(x2/x1)-n-(0) 目标函数值: 9.29374614662676E-5 x1: 31445.1779060676 x2: 2220.82781285504 n: -1.13030905578871 ====== 计算结束 ====== 差的好多，能截图给我看一下算法设置的截面吗？ PS：方程2中的x1原本是x2，我写错了，不过对我的问题不影响。 |
4楼2015-04-10 13:34:25