24小时热门版块排行榜    

查看: 2655  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

抱歉,修正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的回帖

liseven77th

新虫 (初入文坛)

附上我的算法设置界面
因为对数值计算不懂,所以直接用了默认的算法设置及参数
5楼2015-04-10 13:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271材料工程求调剂 +5 .6lL 2026-03-18 5/250 2026-03-19 03:07 by 无懈可击111
[考研] 一志愿天大材料与化工(085600)总分338 +4 蔡大美女 2026-03-13 4/200 2026-03-18 22:47 by walc
[考研] 267一志愿南京工业大学0817化工求调剂 +8 SUICHILD 2026-03-12 8/400 2026-03-18 20:55 by winsuccess
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 307求调剂 +3 冷笙123 2026-03-17 3/150 2026-03-18 09:55 by macy2011
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见