24小时热门版块排行榜    

查看: 2652  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +5 司空. 2026-03-18 5/250 2026-03-18 21:52 by 无懈可击111
[考研] 274求调剂 +5 S.H1 2026-03-18 5/250 2026-03-18 21:27 by guosr9609
[考研] 材料专硕英一数二306 +4 z1z2z3879 2026-03-18 4/200 2026-03-18 20:00 by 楤哥
[考研] 085600材料与化工 +5 安全上岸! 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0703化学调剂 ,六级已过,有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
信息提示
请填处理意见