 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4347 | 回復(fù): 16 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
sunjiahuan銅蟲 (小有名氣)
|
[求助]
matlab求解非線性方程組 已有1人參與
|
|||
|
I1=(1.3403/0.323)*(354.3091/443.6719-lambertw((1.0885e-004/443.6719)*exp((330.7*x1+354.3091)/443.6719)))-x1/(0.323+330.7)+1.019*10^(-6)*(exp(x1/1.3403)-1); I2=(1.3403/0.323)*(349.7334/437.9477-lambertw((1.0744e-004/437.9477)*exp((326.4292*x2+349.7334)/437.9477)))-x2/(0.323+326.4292)+1.019*10^(-6)*(exp(x2/1.3403)-1); I3=(1.3403/0.323)*(337.973/423.2354-lambertw((1.0383e-004/423.2354)*exp((315.4524*x3+337.973)/423.2354)))-x3/(0.323+315.4524)+1.019*10^(-6)*(exp(x3/1.3403)-1); I4=(1.3403/0.323)*(324.4644/406.3363-lambertw((9.9677e-005/406.3363)*exp((302.844*x4+324.4644)/406.3363)))-x4/(0.323+302.844)+1.019*10^(-6)*(exp(x4/1.3403)-1); I5=1.019*10^(-6)*(exp(x5/1.3403)-1)； 上面是5個(gè)方程，I1=f(x1)，I2=f(x2)，I3=f(x3)，I4=f(x4),I5=f(x5); 關(guān)于方程里的lambertw函數(shù)，可以在附件（太陽電池I-V方程顯示求解原理研究及應(yīng)用中的第二章中找到）。 而要求解的方程組是由它們5個(gè)組成的,如下： x1+x2+x3+x4+x5-30=0 I1-I2=0; I1-I3=0; I1-I4=0; I1-I5=0; 上面的方程組已經(jīng)寫好，在fun.m里。 而在dfun.m 中是方程組里的每一個(gè)方程分別對(duì)x1,x2 x3 x4 x5 求一階導(dǎo)數(shù)。其實(shí)就是求方程組的雅克比矩陣； 在newton.m里是求解的編程，是我在這個(gè)網(wǎng)上（http://hi.baidu.com/aillieo/blog ... c9a59647106493.html）找到的，只是把里面的方程及相關(guān)部分改了一下。 我在matlab里輸入newton([0.01 0.01 0.01 0.01 0.01],0.00001,200) 最后得到的結(jié)果是 In newton at 6 Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN. > In newton at 6 ans = NaN NaN NaN NaN NaN。 這種方法是好多學(xué)者用過的，得到的結(jié)果都挺好。所以方法應(yīng)該沒問題。因?yàn)槲业膍atlab只學(xué)了個(gè)皮毛，不知道這里面哪出錯(cuò)了，還請(qǐng)各位高手給指點(diǎn)一下。 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : newton.m
- 附件 2 : dfun.m
- 附件 3 : fun.m
- 附件 4 : 太陽電池I_V方程顯式求解原理研究及應(yīng)用.nh
2012-02-22 16:56:51, 775 bytes
2012-02-22 16:56:51, 120 bytes
2012-02-22 16:56:52, 969 bytes
2012-02-22 16:57:53, 4.29 M
» 收錄本帖的淘帖專輯推薦
 材料+新能源+催化+外語 |
» 猜你喜歡
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有11人回復(fù)
-
332求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料專碩 085600求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有14人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有14人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有11人回復(fù)
-
環(huán)境科學(xué)與工程334分求調(diào)劑
已經(jīng)有8人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助解復(fù)雜非線性方程組的好的方法
已經(jīng)有24人回復(fù)
- matlab如何求解帶參數(shù)的線性方程組
已經(jīng)有6人回復(fù)
- 跪求求解線性方程組的c語言程序
已經(jīng)有3人回復(fù)
- 用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？
已經(jīng)有11人回復(fù)
- matlab如何求解一個(gè)非線性微分方程組
已經(jīng)有8人回復(fù)
- matlab求解高階方程組
已經(jīng)有3人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 求Matlab解非線性方程代碼
已經(jīng)有12人回復(fù)
- 非線性方程組的解法
已經(jīng)有5人回復(fù)
- matlab求解非線性方程組，錯(cuò)誤提示怎么解決
已經(jīng)有5人回復(fù)
- MATLAB求解代數(shù)方程組問題，請(qǐng)幫忙
已經(jīng)有9人回復(fù)
- 求助 MATLAB解方程組-fslove
已經(jīng)有7人回復(fù)
- 求助matlab解非線性方程。求代碼
已經(jīng)有5人回復(fù)
- 求助matlab---fsolve解非線性方程組
已經(jīng)有6人回復(fù)
- matlab求解非線性方程組，求助！
已經(jīng)有6人回復(fù)
- 求matlab求解一個(gè)方程組的問題
已經(jīng)有6人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- matlab解方程組
已經(jīng)有8人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實(shí)數(shù)解！
已經(jīng)有18人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
sunjiahuan
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1798.4
- 散金: 43
- 紅花: 2
- 帖子: 221
- 在線: 45小時(shí)
- 蟲號(hào): 916576
- 注冊(cè): 2009-11-30
- 性別: GG
- 專業(yè): 可再生與替代能源利用中的
4樓2012-02-22 20:46:34
dbb627
榮譽(yù)版主 (著名寫手)
-
專家經(jīng)驗(yàn): +4 - 仿真EPI: 6
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時(shí)
- 蟲號(hào): 149791
- 注冊(cè): 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計(jì)算模擬
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
ben_ladeng(金幣+2): 謝謝應(yīng)助 2012-02-22 21:59:15
ben_ladeng(專家考核): 2012-02-22 21:59:24
caemechanics(金幣+5): 感謝應(yīng)助 2012-02-23 16:18:30
ben_ladeng(金幣+2): 謝謝應(yīng)助 2012-02-22 21:59:15
ben_ladeng(專家考核): 2012-02-22 21:59:24
caemechanics(金幣+5): 感謝應(yīng)助 2012-02-23 16:18:30
|
錯(cuò)誤提示信息是矩陣奇異， 可能初值不太合適 我用自帶fsolve解了下 結(jié)果不是很好 [/code] function aa [X,FVAL,EXITFLAG,OUTPUT] =fsolve(@fun,[18.7783 0.01 18.1315 18.0062 -2],optimset('MaxFunEvals',1000)) function f=fun(x) x1=x(1); x2=x(2); x3=x(3); x4=x(4); x5=x(5); f1=x1+x2+x3+x4+x5-30; f2=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-4.1495*(0.7986-lambertw(2.4533e-007*exp((326.4292*x2+349.7334)/437.9477)))+x2/(0.323+326.4292)-1.019*10^(-6)*(exp(x2/1.3403)-1); f3=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-4.1495*(0.7985-lambertw(2.4532e-007*exp((315.4524*x3+337.973)/423.2354)))+x3/(0.323+315.4524)-1.019*10^(-6)*(exp(x3/1.3403)-1); f4=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-4.1495*(0.7985-lambertw( 2.4531e-007*exp((302.844*x4+324.4644)/406.3363)))+x4/(0.323+302.844)-1.019*10^(-6)*(exp(x4/1.3403)-1); f5=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-1.019*10^(-6)*(exp(x5/1.3403)-1); f=[f1 f2 f3 f4 f5]; [/code] Solver stopped prematurely. fsolve stopped because it exceeded the function evaluation limit, options.MaxFunEvals = 1000 (the selected value). X = 18.7895 19.0894 19.0963 19.1094 -46.0868 FVAL = -0.0022 -0.0180 -0.0172 -0.0175 2.7252 EXITFLAG = 0 OUTPUT = iterations: 180 funcCount: 1001 algorithm: 'trust-region dogleg' firstorderopt: 0.0022 message: [1x143 char] |
» 本帖已獲得的紅花（最新10朵）
sunjiahuan
2樓2012-02-22 19:55:48
sunjiahuan
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1798.4
- 散金: 43
- 紅花: 2
- 帖子: 221
- 在線: 45小時(shí)
- 蟲號(hào): 916576
- 注冊(cè): 2009-11-30
- 性別: GG
- 專業(yè): 可再生與替代能源利用中的
3樓2012-02-22 20:29:00
sunjiahuan
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1798.4
- 散金: 43
- 紅花: 2
- 帖子: 221
- 在線: 45小時(shí)
- 蟲號(hào): 916576
- 注冊(cè): 2009-11-30
- 性別: GG
- 專業(yè): 可再生與替代能源利用中的
5樓2012-02-22 20:53:37
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 292求調(diào)劑 +8 | 是妍子也是研子 2026-03-30 | 8/400 |
|
|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 7/350 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25 | RichLi_ 2026-03-25 | 25/1250 |
|
|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪�。�！ 2026-03-25 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|