| 查看: 2651 | 回復(fù): 10 | ||||
liseven77th新蟲 (初入文壇)
|
[求助]
求助用matlab解多元指數(shù)型方程組 已有1人參與
|
|||
|
方程如圖。直接用matlab解方程組解不出。 參考 http://wenku.baidu.com/link?url= ... jMX4Spsljya0IuANYga 這個帖子的方法用牛頓法迭代求解,帶入初值410,450,30迭代得不到結(jié)果。 方程近似解x1=403 x2=446 n=30 有大神能指點(diǎn)一下么 IMG_20150409_230048.jpg |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
若限定x1=403 x2=446 n=30左右 僅能找到一般的解: 目標(biāo)函數(shù)值: 8.94512267343522E-5 x1: 400.003316069477 x2: 421.966847460639 n: 56.0434455220799 不限制未知數(shù)范圍則能找到較好的解: 目標(biāo)函數(shù)值: 2.71050543121376E-20 x1: 80.2313890183921 x2: 802.313890183922 n: 1.30102999566398 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
抱歉,修正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 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

新蟲 (初入文壇)
|
感謝你的回答。 我用1stopt軟件計算了一遍結(jié)果如下 ====== 結(jié)果 ====== 迭代數(shù): 106 計算用時(時:分:秒:毫秒): 00:00:01:04 計算中止原因: 達(dá)到收斂判定標(biāo)準(zhǔn) 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 函數(shù)表達(dá)式 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) 目標(biāo)函數(shù)值: 9.29374614662676E-5 x1: 31445.1779060676 x2: 2220.82781285504 n: -1.13030905578871 ====== 計算結(jié)束 ====== 差的好多,能截圖給我看一下算法設(shè)置的截面嗎? PS:方程2中的x1原本是x2,我寫錯了,不過對我的問題不影響。 |
新蟲 (初入文壇)
新蟲 (初入文壇)
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
欲求得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è)置默認(rèn)即可 計算結(jié)果: 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 函數(shù)表達(dá)式 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) 目標(biāo)函數(shù)值: 2.65190088855408E-12 x1: 402.479314242841 x2: 444.735331831014 n: 30.0066408207363 |

新蟲 (初入文壇)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +3 | z1z2z3879 2026-03-18 | 3/150 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評估A類院;瘜W(xué)工程與工藝,誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|