| 查看: 2683 | 回复: 10 | ||||
[求助]
求助用matlab解多元指数型方程组 已有1人参与
|
||||
|
方程如图。直接用matlab解方程组解不出。 参考 http://wenku.baidu.com/link?url= ... jMX4Spsljya0IuANYga 这个帖子的方法用牛顿法迭代求解，带入初值410,450,30迭代得不到结果。 方程近似解x1=403 x2=446 n=30 有大神能指点一下么  IMG_20150409_230048.jpg |
回复此楼» 猜你喜欢
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有145人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 如何求解多元二阶微分方程组
已经有7人回复
- 求助matlab高手指点：大型非线性方程组（牛顿拉佛森法）
已经有7人回复
- 双指数方程求解（MATLAB或者直接给公式都可以）
已经有9人回复
- MATLAB求解五元非线性方程组，在线等
已经有7人回复
- 利用matlab求一阶变系数微分方程组的数值解。
已经有7人回复
- matlab如何求解多元多次方程，包括绘图
已经有7人回复
- matlab解方程组求助攻
已经有4人回复
- Matlab多元高次方程组的求解问题
已经有19人回复
- 求助用matlab拟合动力学方程
已经有13人回复
- matlab如何求解带参数的线性方程组
已经有6人回复
- 求助用matlab或者（1st0pt)编程解方程组
已经有3人回复
- matlab方程求教
已经有4人回复
- 请教大家一个MATLAB的方程求解问题
已经有6人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 求助：matlab解17个未知数，17个方程。
已经有4人回复
- matlab解偏微分方程求助
已经有12人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- Matlab求解偏微分方程组
已经有22人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-04-10 08:49:55
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
10楼2015-04-12 21:03:08
|
感谢你的回答。 我用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
5楼2015-04-10 13:42:06
6楼2015-04-10 13:44:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liseven77th: 金币+10, ★★★★★最佳答案 2015-04-10 22:06:25
liseven77th: 金币+10, ★★★★★最佳答案 2015-04-10 22:06:25
|
欲求得x1=403 x2=446 n=30左右的解，对未知数限定范围即可： Parameters x1[100,],x2[100,],n[30,]; //Variable ; Function x1/207000*3/7*(x1/433)^(n-1)-0.0001=0; x2/207000*3/7*(x2/433)^(n-1)-0.002=0; Ln(20)/Ln(x2/x1)-n=0; 算法设置默认即可 计算结果： 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 函数表达式 1: x1/207000*3/7*(x1/433)^(n-1)-0.0001-(0) 2: x2/207000*3/7*(x2/433)^(n-1)-0.002-(0) 3: ln(20)/ln(x2/x1)-n-(0) 目标函数值: 2.65190088855408E-12 x1: 402.479314242841 x2: 444.735331831014 n: 30.0066408207363 |
7楼2015-04-10 17:00:24
8楼2015-04-10 22:06:00
9楼2015-04-12 20:23:46
