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

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

24小時(shí)熱門版塊排行榜

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

返回列表

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

sunjiahuan

銅蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1798.4
散金: 43
紅花: 2
帖子: 221
在線: 45小時(shí)
蟲(chóng)號(hào): 916576
注冊(cè): 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個(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)寫(xiě)好，在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)一下。

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(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 : 太陽(yáng)電池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的回帖

sunjiahuan

銅蟲(chóng) (小有名氣)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

5樓2012-02-22 20:53:37

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

查看全部 17 個(gè)回答

dbb627

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

專家經(jīng)驗(yàn): +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲(chóng)號(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

錯(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]

贊一下(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

銅蟲(chóng) (小有名氣)

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

送鮮花一朵

引用回帖:

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

贊一下

回復(fù)此樓

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

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

sunjiahuan

銅蟲(chóng) (小有名氣)

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

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

贊一下

回復(fù)此樓

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

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

查看全部 17 個(gè)回答

普通表情龍兔虎貓

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

[考研] 省雙一流重點(diǎn)一本大學(xué)招收調(diào)劑 +3	wwwwffffff 2026-03-31	5/250	2026-03-31 19:49 by 曾曾曾！
[考研] 311求調(diào)劑 +7	李芷新1 2026-03-31	7/350	2026-03-31 18:10 by 無(wú)際的草原
[考研] 085602 307分求調(diào)劑 +9	不知道叫什么！ 2026-03-26	9/450	2026-03-31 18:10 by 544594351
[考研] 環(huán)境工程 085701，267求調(diào)劑 +13	minht 2026-03-29	13/650	2026-03-31 17:47 by 544594351
[考研] 調(diào)劑 +4	好好讀書(shū)。 2026-03-28	5/250	2026-03-31 16:29 by 我就是我ok
[考研] 286求調(diào)劑 +6	Faune 2026-03-30	6/300	2026-03-31 14:37 by jp9609
[考研] 266分，求材料相關(guān)專業(yè)調(diào)劑 +10	哇呼哼呼哼 2026-03-30	12/600	2026-03-31 11:00 by 熊一刀
[考研] 310求調(diào)劑 +10	爭(zhēng)取九點(diǎn)睡 2026-03-30	10/500	2026-03-30 16:45 by ztnimte
[考研] 324求調(diào)劑 +9	hanamiko 2026-03-26	11/550	2026-03-30 14:27 by JourneyLucky
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 材料與化工304求B區(qū)調(diào)劑 +4	邱gl 2026-03-26	7/350	2026-03-30 08:39 by 探123
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 296求調(diào)劑 +4	汪�。�！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4	為學(xué)666 2026-03-24	4/200	2026-03-25 11:27 by BruceLiu320