| 5 | 1/1 | 返回列表 |
| 查看: 2688 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
liseven77th新蟲 (初入文壇)
|
[求助]
求助用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 |
新蟲 (初入文壇)
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +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)驗: +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 |

新蟲 (初入文壇)
|
感謝你的回答。 我用1stopt軟件計算了一遍結(jié)果如下 ====== 結(jié)果 ====== 迭代數(shù): 106 計算用時(時:分:秒:毫秒): 00:00:01:04 計算中止原因: 達到收斂判定標(biāo)準(zhǔn) 優(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) 目標(biāo)函數(shù)值: 9.29374614662676E-5 x1: 31445.1779060676 x2: 2220.82781285504 n: -1.13030905578871 ====== 計算結(jié)束 ====== 差的好多,能截圖給我看一下算法設(shè)置的截面嗎? PS:方程2中的x1原本是x2,我寫錯了,不過對我的問題不影響。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +16 | 新時代材料 2026-03-27 | 16/800 |
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 研研,接電話 2026-03-24 | 8/400 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考博] 26申博 +3 | 加油沖。 2026-03-26 | 3/150 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|