| 查看: 2552 | 回复: 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人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有185人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 金币: 68122.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-04-10 08:49:55
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68122.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
- 金币: 68122.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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68122.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
10楼2015-04-12 21:03:08