| 查看: 4351 | 回復(fù): 16 | ||||
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ù),可以在附件(太陽(yáng)電池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é)者用過(guò)的,得到的結(jié)果都挺好。所以方法應(yīng)該沒(méi)問(wèn)題。因?yàn)槲业膍atlab只學(xué)了個(gè)皮毛,不知道這里面哪出錯(cuò)了,還請(qǐng)各位高手給指點(diǎn)一下。 |
材料+新能源+催化+外語(yǔ) |
銅蟲 (小有名氣)
榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +4 |
|
存為.m文件 命令窗口 >> syms x1 x2 x3 x4 x5 >> F=[x1+x2+x3+x4+x5-30; 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); 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); 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); 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)]; [r,n]=mulNewton(F,[18.7783 0.01 18.1315 18.0062 -2]) Warning: Matrix is singular to working precision. > In mulNewton at 17 r = NaN NaN NaN NaN NaN n = 2 [r,n]=mulNewton(F,[18.7783 20 18.1315 18.0062 -20]) r = 20.0741 20.0739 20.0728 20.0718 -50.2925 n = 5 |

榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +4 |

榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +4 |
|
錯(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] |

榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +4 |
|
將上述代碼拷到editor里面,存為.m文件,點(diǎn)擊editor中run按鈕 結(jié)果如下 Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. X = 20.0741 20.0739 20.0728 20.0718 -50.2925 FVAL = 1.0e-007 * 0.0000 -0.8502 -0.8544 -0.8544 -0.8544 EXITFLAG = 1 OUTPUT = iterations: 13 funcCount: 74 algorithm: 'trust-region dogleg' firstorderopt: 4.6500e-007 message: [1x695 char] |

銅蟲 (小有名氣)
送鮮花一朵 |
您好,非常感謝你的幫助,您能不能把你的matlab代碼發(fā)到我的郵箱里么sunjiahuan0370@sina.com?我把你貼出來(lái)的代碼移到我的matlab里居然沒(méi)成功。謝謝了。 |
銅蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
鐵桿木蟲 (著名寫手)
ส็็็

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 11408總分309,一志愿東南大學(xué)求調(diào)劑,不挑專業(yè) +4 | 天賦帶到THU 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +14 | jhybd 2026-03-23 | 19/950 |
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|