liseven77th

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 49.8
散金: 5
帖子: 36
在線: 21.9小時
蟲號: 3182479
注冊: 2014-05-05
專業(yè): 結(jié)構(gòu)工程

[求助] 求助用matlab解多元指數(shù)型方程組已有1人參與

方程如圖。直接用matlab解方程組解不出。
參考
http://wenku.baidu.com/link?url= ... jMX4Spsljya0IuANYga
這個帖子的方法用牛頓法迭代求解，帶入初值410,450,30迭代得不到結(jié)果。
方程近似解x1=403 x2=446 n=30
有大神能指點一下么

IMG_20150409_230048.jpg

回復(fù)此樓

» 猜你喜歡

0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有103人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab 如何求解多元二階微分方程組已經(jīng)有7人回復(fù)
求助matlab高手指點：大型非線性方程組（牛頓拉佛森法）已經(jīng)有7人回復(fù)
雙指數(shù)方程求解（MATLAB或者直接給公式都可以）已經(jīng)有9人回復(fù)
MATLAB求解五元非線性方程組，在線等已經(jīng)有7人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
matlab如何求解多元多次方程，包括繪圖已經(jīng)有7人回復(fù)
matlab解方程組求助攻已經(jīng)有4人回復(fù)
Matlab多元高次方程組的求解問題已經(jīng)有19人回復(fù)
求助用matlab擬合動力學(xué)方程已經(jīng)有13人回復(fù)
matlab如何求解帶參數(shù)的線性方程組已經(jīng)有6人回復(fù)
求助用matlab或者（1st0pt)編程解方程組已經(jīng)有3人回復(fù)
matlab方程求教已經(jīng)有4人回復(fù)
請教大家一個MATLAB的方程求解問題已經(jīng)有6人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
求助，用matlab解微分方程組，希望高手能給予指點已經(jīng)有17人回復(fù)
求助：matlab解17個未知數(shù)，17個方程。已經(jīng)有4人回復(fù)
matlab解偏微分方程求助已經(jīng)有12人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復(fù)
求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
【求助】matlab 解矩陣方程重金懸賞解決方案已經(jīng)有10人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

1樓 2015-04-09 23:12:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liseven77th: 金幣+10, ★★★★★最佳答案 2015-04-10 22:06:25

引用回帖:

4樓: Originally posted by liseven77th at 2015-04-10 13:34:25
感謝你的回答。
我用1stopt軟件計算了一遍結(jié)果如下
====== 結(jié)果 ======

迭代數(shù): 106
計算用時(時:分:秒:毫秒): 00:00:01:04
計算中止原因: 達到收斂判定標準
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通 ...

欲求得x1=403 x2=446 n=30左右的解，對未知數(shù)限定范圍即可：
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;

算法設(shè)置默認即可

計算結(jié)果：
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
函數(shù)表達式 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)
目標函數(shù)值: 2.65190088855408E-12
x1: 402.479314242841
x2: 444.735331831014
n: 30.0066408207363

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

7樓2015-04-10 17:00:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個回答

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

若限定x1=403 x2=446 n=30左右
僅能找到一般的解：
目標函數(shù)值: 8.94512267343522E-5

x1: 400.003316069477
x2: 421.966847460639
n: 56.0434455220799

不限制未知數(shù)范圍則能找到較好的解：
目標函數(shù)值: 2.71050543121376E-20

x1: 80.2313890183921
x2: 802.313890183922
n: 1.30102999566398

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2015-04-10 08:49:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

抱歉，修正2樓的說法，我在解方程時候，把0.002寫成了0.0002，更正后結(jié)果如下：

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;

計算結(jié)果：
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

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

3樓2015-04-10 09:11:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liseven77th

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 49.8
散金: 5
帖子: 36
在線: 21.9小時
蟲號: 3182479
注冊: 2014-05-05
專業(yè): 結(jié)構(gòu)工程

引用回帖:

3樓: Originally posted by 月只藍 at 2015-04-10 09:11:40
抱歉，修正2樓的說法，我在解方程時候，把0.002寫成了0.0002，更正后結(jié)果如下：

1stopt代碼：
Parameters x1,x2,n;
//Variable ;
Function x1/207000*3/7*(x1/433)^(n-1)-0.0001=0;
x1/207000*3/7* ...

感謝你的回答。
我用1stopt軟件計算了一遍結(jié)果如下
====== 結(jié)果 ======

迭代數(shù): 106
計算用時(時:分:秒:毫秒): 00:00:01:04
計算中止原因: 達到收斂判定標準
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
函數(shù)表達式 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)
目標函數(shù)值: 9.29374614662676E-5
x1: 31445.1779060676
x2: 2220.82781285504
n: -1.13030905578871

====== 計算結(jié)束 ======
差的好多，能截圖給我看一下算法設(shè)置的截面嗎？
PS：方程2中的x1原本是x2，我寫錯了，不過對我的問題不影響。

贊一下

回復(fù)此樓

4樓2015-04-10 13:34:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 北科281學(xué)碩材料求調(diào)劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 306求0703調(diào)劑一志愿華中師范 +6	紙魚ly 2026-03-21	6/300	2026-03-23 12:08 by Iveryant
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 085600材料與化工 +10	安全上岸！ 2026-03-16	10/500	2026-03-23 09:31 by Jy?
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風里
[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 305分求調(diào)劑（食品工程） +4	Sxy112 2026-03-21	6/300	2026-03-22 15:26 by 無懈可擊111
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3	123阿標 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客