| 5 | 1/1 | 返回列表 |
| 查看: 2666 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
liseven77th新蟲 (初入文壇)
|
[求助]
求助用matlab解多元指數(shù)型方程組 已有1人參與
|
|||
|
方程如圖。直接用matlab解方程組解不出。 參考 http://wenku.baidu.com/link?url= ... jMX4Spsljya0IuANYga 這個(gè)帖子的方法用牛頓法迭代求解,帶入初值410,450,30迭代得不到結(jié)果。 方程近似解x1=403 x2=446 n=30 有大神能指點(diǎn)一下么 IMG_20150409_230048.jpg |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
抱歉,修正2樓的說法,我在解方程時(shí)候,把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; 計(jì)算結(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 |
|
若限定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 |

新蟲 (初入文壇)
|
感謝你的回答。 我用1stopt軟件計(jì)算了一遍結(jié)果如下 ====== 結(jié)果 ====== 迭代數(shù): 106 計(jì)算用時(shí)(時(shí):分:秒:毫秒): 00:00:01:04 計(jì)算中止原因: 達(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 ====== 計(jì)算結(jié)束 ====== 差的好多,能截圖給我看一下算法設(shè)置的截面嗎? PS:方程2中的x1原本是x2,我寫錯(cuò)了,不過對(duì)我的問題不影響。 |
新蟲 (初入文壇)
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 269專碩求調(diào)劑 +5 | 金恩貝 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿211,0703化學(xué)310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 |  ̄^ ̄゜汗 2026-03-19 | 4/200 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|