版塊導(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í)熱門(mén)版塊排行榜

北京石油化工學(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
專(zhuān)業(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
專(zhuān)業(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 ...

您好，非常感謝你的幫助，您能不能把你的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的回帖

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

dbb627

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

專(zhuān)家經(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
專(zhuān)業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

感謝參與，應(yīng)助指數(shù) +1
ben_ladeng(金幣+2): 謝謝應(yīng)助 2012-02-22 21:59:15
ben_ladeng(專(zhuān)家考核): 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
專(zhuān)業(yè): 可再生與替代能源利用中的

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

不好意思，剛看到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ù)此樓

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

已閱回復(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
專(zhuān)業(yè): 可再生與替代能源利用中的

引用回帖:

贊一下

回復(fù)此樓

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

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

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

普通表情龍兔虎貓

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

[考研] 生物考研337分求調(diào)劑 +3	cgxin 2026-03-30	5/250	2026-03-31 11:47 by cgxin
[考研] 285求調(diào)劑 +6	AZMK 2026-03-30	10/500	2026-03-31 11:41 by AZMK
[考研] 262求調(diào)劑 +7	ZZ..000 2026-03-30	8/400	2026-03-31 10:05 by cal0306
[考研] 335求調(diào)劑 +3	321* 2026-03-31	3/150	2026-03-31 09:00 by 功夫瘋狂
[考研] 環(huán)境工程 085701，267求調(diào)劑 +11	minht 2026-03-29	11/550	2026-03-31 08:21 by sunshine0013
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4	6266jl 2026-03-30	4/200	2026-03-30 21:48 by zhuangyan123
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +7	蘿卜燉青菜 2026-03-28	8/400	2026-03-30 19:59 by 無(wú)際的草原
[考研] 材料化工340求調(diào)劑 +3	jhx777 2026-03-30	3/150	2026-03-30 17:54 by JourneyLucky
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +8	26要上岸 2026-03-28	8/400	2026-03-30 16:39 by herarysara
[考研] 329求調(diào)劑 +8	星野? 2026-03-26	8/400	2026-03-30 13:41 by chemdavid
[考研] 一志愿南航 335分 \| 0856 \| GPA 4.07 \| 有科研經(jīng)歷 +8	cccchenso 2026-03-29	8/400	2026-03-29 23:53 by 我是小康
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +7	AZMK 2026-03-28	11/550	2026-03-28 21:01 by xxxsssccc
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣(mài)魚(yú) 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671
[考研] 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