 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (3540)
- >考研 (721)
- >導(dǎo)師招生 (595)
- >蟲(chóng)友互識(shí) (224)
- >文獻(xiàn)求助 (89)
- >休閑灌水 (75)
- >碩博家園 (48)
- >招聘信息布告欄 (33)
- >考博 (33)
- >博后之家 (18)
- >綠色求助(高懸賞) (17)
- >論文投稿 (16)
- >基金申請(qǐng) (12)
- >攻關(guān)文獻(xiàn)(高獎(jiǎng)勵(lì)) (8)
- >教師之家 (8)
- >外文書(shū)籍求助 (7)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1499 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
天才超人新蟲(chóng) (初入文壇)
|
[求助]
Matlab解微分方程 求助 已有1人參與
|
|
|
各位大神走過(guò)路過(guò) 拜托幫幫忙 🙏 我需要用Matlab解3個(gè)ODE方程，用的是bvp4c solinit=bvpinit(linspace(0,LA,1000),[E_A U_A C_OH_A]); options=bvpset('RelTol',1e-3,'AbsTol',[1e-12 1e-12 1e-12]); sol=bvp4c(@odefun,@odebc,solinit,options); 方程的file如下： function [dydx]=odefun(y) global F R T Z_Na Z_Cl Z_OH Z_H C_Na_A C_Cl_A i U_A FixedCharge D_OH e D_H C_Na=C_Na_A*exp(-Z_Na*F/R/T*(y(2)-U_A)); % conc of Na+ C_Cl=C_Cl_A*exp(-Z_Cl*F/R/T*(y(2)-U_A)); % conc of Cl- dydx(1)=F/e*(C_Na-C_Cl+10^(-8)/y(3)-y(3)+FixedCharge); % Electric field dydx(2)=-y(1); % Electric potential J_H=-D_H*(-10^(-8)/y(3)^2*dydx(3)+Z_H*F/R/T*10^(-8)/y(3)*dydx(2)); %Flux of H+ dydx(3)=-(i/F-Z_H*J_H)/Z_OH/D_OH-Z_OH*F/RT*y(3)*dydx(2); % Conc of OH- end 邊界條件如下： function [res]=odebc(ya,~) global C_OH_A U_A E_A res(1)=ya(1) - E_A; res(2)=ya(2) - U_A; res(3)=ya(3) - C_OH_A; res=res'; end 運(yùn)行后一直收到如下的錯(cuò)誤信息： Error using BPM_ode_Dec1417>odefun Too many input arguments. Error in bvparguments (line 105) testODE = ode(x1,y1,odeExtras{:}); Error in bvp4c (line 130) bvparguments(solver_name,ode,bc,solinit,options,varargin); Error in BPM_ode_Dec1417 (line 40) sol=bvp4c(@odefun,@odebc,solinit,options); 盯著程序兩天了，各種variable換來(lái)?yè)Q去，實(shí)在是解不出來(lái)。所以來(lái)拜托各位大神了 🙏🙏 謝謝！ |
回復(fù)此樓» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有7人回復(fù)
-
270調(diào)劑
已經(jīng)有7人回復(fù)
-
275學(xué)碩081000服從調(diào)劑到其他專(zhuān)業(yè)，保不住本專(zhuān)業(yè)了
已經(jīng)有6人回復(fù)
-
0710生物學(xué)，325求調(diào)劑
已經(jīng)有3人回復(fù)
-
296材料專(zhuān)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
385分 生物學(xué)（071000）求調(diào)劑
已經(jīng)有5人回復(fù)
-
377求調(diào)劑
已經(jīng)有3人回復(fù)
-
【求調(diào)劑】新能源材料本科，一志愿211，初試321
已經(jīng)有6人回復(fù)
-
316求調(diào)劑
已經(jīng)有11人回復(fù)
-
321求調(diào)劑
已經(jīng)有9人回復(fù)
3樓2017-12-16 13:36:47
2樓2017-12-16 07:09:13
天才超人
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 105
- 帖子: 22
- 在線: 9.3小時(shí)
- 蟲(chóng)號(hào): 1945549
- 注冊(cè): 2012-08-19
- 專(zhuān)業(yè): 環(huán)境工程
4樓2017-12-17 03:39:10
天才超人
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 105
- 帖子: 22
- 在線: 9.3小時(shí)
- 蟲(chóng)號(hào): 1945549
- 注冊(cè): 2012-08-19
- 專(zhuān)業(yè): 環(huán)境工程
5樓2017-12-18 06:20:18
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請(qǐng)]
請(qǐng)問(wèn)共同通訊和共同一作的認(rèn)可度問(wèn)題
 30+3
|
psa1234 2026-04-01 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 288求調(diào)劑 一志愿哈工大 材料與化工 +30 | 洛神哥哥 2026-03-31 | 30/1500 |
|
|
[考研] 材料工程322分 +7 | 哈哈哈吼吼吼哈 2026-04-01 | 7/350 |
|
|
|
[考研] 0710生物學(xué)求調(diào)劑 +7 | manman511 2026-04-01 | 7/350 |
|
|
|
[考研] 085602化學(xué)工程268分蹲調(diào)劑 +8 | 月照花林。 2026-04-01 | 8/400 |
|
|
|
[考研] 290求調(diào)劑 +5 | dfffsar 2026-03-29 | 5/250 |
|
|
|
[考研] 362求調(diào)劑 +13 | 西南交材料專(zhuān)碩3 2026-03-31 | 13/650 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +7 | @taotao 2026-03-30 | 7/350 |
|
|
|
[考研] 一志愿中農(nóng)0710生物學(xué)，微生物方向總分338求調(diào)劑 +3 | 柒xxxx. 2026-03-26 | 3/150 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +12 | 26要上岸 2026-03-28 | 12/600 |
|
|
|
[考研] 0805求調(diào)劑 +6 | 是水分 2026-03-31 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
|
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4 | 6266jl 2026-03-30 | 8/400 |
|
|
|
[考研] 315求調(diào)劑 +6 | akie... 2026-03-28 | 7/350 |
|
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
|
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專(zhuān)碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
|
[考研] 材料專(zhuān)碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|