24小时热门版块排行榜    

查看: 2654  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liseven77th

新虫 (初入文坛)

[求助] 求助用matlab解多元指数型方程组 已有1人参与

方程如图。直接用matlab解方程组解不出。
参考
http://wenku.baidu.com/link?url= ... jMX4Spsljya0IuANYga
这个帖子的方法用牛顿法迭代求解,带入初值410,450,30迭代得不到结果。
方程近似解x1=403 x2=446 n=30
有大神能指点一下么
求助用matlab解多元指数型方程组
IMG_20150409_230048.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

抱歉,修正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
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2015-04-10 09:11:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
若限定x1=403 x2=446 n=30左右
仅能找到一般的解:
目标函数值: 8.94512267343522E-5

x1: 400.003316069477
x2: 421.966847460639
n: 56.0434455220799

不限制未知数范围则能找到较好的解:
目标函数值: 2.71050543121376E-20

x1: 80.2313890183921
x2: 802.313890183922
n: 1.30102999566398
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2015-04-10 08:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liseven77th

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 月只蓝 at 2015-04-10 09:11:40
抱歉,修正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* ...

感谢你的回答。
我用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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liseven77th

新虫 (初入文坛)

附上我的算法设置界面
因为对数值计算不懂,所以直接用了默认的算法设置及参数
5楼2015-04-10 13:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 求材料调剂 +8 隔壁陈先生 2026-03-12 8/400 2026-03-18 22:19 by li123456789.
[考研] 一志愿吉林大学材料学硕321求调剂 +4 Ymlll 2026-03-18 6/300 2026-03-18 22:15 by li123456789.
[考研] 321求调剂 +3 何润采123 2026-03-18 3/150 2026-03-18 21:27 by li123456789.
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 271材料工程求调剂 +4 .6lL 2026-03-18 4/200 2026-03-18 20:05 by 楤哥
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见