版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

sunjiahuan

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1798.4
散金: 43
紅花: 2
帖子: 221
在線: 45小時
蟲號: 916576
注冊: 2009-11-30
性別: GG
專業(yè): 可再生與替代能源利用中的

[求助] 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個方程，I1=f(x1)，I2=f(x2)，I3=f(x3)，I4=f(x4),I5=f(x5); 關(guān)于方程里的lambertw函數(shù)，可以在附件（太陽電池I-V方程顯示求解原理研究及應(yīng)用中的第二章中找到）。

而要求解的方程組是由它們5個組成的,如下：
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 中是方程組里的每一個方程分別對x1,x2 x3 x4 x5 求一階導(dǎo)數(shù)。其實(shí)就是求方程組的雅克比矩陣；

在newton.m里是求解的編程，是我在這個網(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)該沒問題。因為我的matlab只學(xué)了個皮毛，不知道這里面哪出錯了，還請各位高手給指點(diǎn)一下。

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : newton.m

2012-02-22 16:56:51, 775 bytes

附件 2 : dfun.m

2012-02-22 16:56:51, 120 bytes

附件 3 : fun.m

2012-02-22 16:56:52, 969 bytes

附件 4 : 太陽電池I_V方程顯式求解原理研究及應(yīng)用.nh

2012-02-22 16:57:53, 4.29 M

1樓 2012-02-22 17:39:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

★ ★
xiegangmai: 金幣+2, 專家考核, 謝謝參與！ 2013-04-17 23:53:43

猜初值很頭疼的，下面是正解：

x1: 2.60751007153198
x2: 2.56027779100342
x3: 2.43893790346484
x4: 2.29898071559569
x5: 20.094293518415

贊一下

回復(fù)此樓

9樓2013-04-17 10:28:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 17 個回答

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(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

錯誤提示信息是矩陣奇異，
可能初值不太合適
我用自帶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]

贊一下(2人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

sunjiahuan

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2012-02-22 19:55:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunjiahuan

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1798.4
散金: 43
紅花: 2
帖子: 221
在線: 45小時
蟲號: 916576
注冊: 2009-11-30
性別: GG
專業(yè): 可再生與替代能源利用中的

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2012-02-22 19:55:48:
錯誤提示信息是矩陣奇異，
可能初值不太合適
我用自帶fsolve解了下
結(jié)果不是很好
[/code]
function aa
[X,FVAL,EXITFLAG,OUTPUT] =fsolve(@fun,[18.7783 0.01 18.1315 18.0062 -2],optimset('MaxFun ...

您好，非常感謝你的幫助，您能不能把你的matlab代碼發(fā)到我的郵箱里么sunjiahuan0370@sina.com？我把你貼出來的代碼移到我的matlab里居然沒成功。謝謝了。

贊一下

回復(fù)此樓

3樓2012-02-22 20:29:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunjiahuan

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1798.4
散金: 43
紅花: 2
帖子: 221
在線: 45小時
蟲號: 916576
注冊: 2009-11-30
性別: GG
專業(yè): 可再生與替代能源利用中的

sunjiahuan: 回帖置頂 2012-02-22 20:46:39

不好意思，剛看到I1，I2,I3,I4方程中最后一項的里的x1，x2,x3,x4, 應(yīng)該為-x1,-x2,-x3,-x4,而x5不變。所以我附件里fun.m的方程組里的也要相應(yīng)改變，因無法編輯，在此說明一下。比如I1中的最后一項1.019*10^(-6)*(exp(x1/1.3403)-1);應(yīng)改為1.019*10^(-6)*(exp(-x1/1.3403)-1);

贊一下

回復(fù)此樓

4樓2012-02-22 20:46:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 17 個回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 284求調(diào)劑 +12	junqihahaha 2026-03-26	13/650	2026-03-30 13:48 by Iveryant
[考研] 311求調(diào)劑 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6	niko- 2026-03-27	6/300	2026-03-30 10:25 by herarysara
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4	12.27 2026-03-27	4/200	2026-03-29 14:40 by zhshch
[考研] 22408 359分調(diào)劑 +4	Qshers 2026-03-27	6/300	2026-03-29 14:16 by Qshers
[考研] 求收留 +5	1943443204 2026-03-28	5/250	2026-03-29 12:11 by 無際的草原
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落睿可思
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 296求調(diào)劑 +4	汪��？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪�。浚�
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[考研] 一志愿南航材料專317分求調(diào)劑 +5	炸呀炸呀炸薯條 2026-03-23	5/250	2026-03-24 16:52 by 星空星月
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

sunjiahuan

[求助] matlab求解非線性方程組 已有1人參與

» 本帖附件資源列表

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

dingd

【答案】應(yīng)助回帖

dbb627

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

sunjiahuan

sunjiahuan

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] matlab求解非線性方程組已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: