 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1498 | 回復: 7 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
天才超人新蟲 (初入文壇)
|
[求助]
Matlab解微分方程 求助 已有1人參與
|
|
|
各位大神走過路過 拜托幫幫忙 🙏 我需要用Matlab解3個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 運行后一直收到如下的錯誤信息： 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換來換去，實在是解不出來。所以來拜托各位大神了 🙏🙏 謝謝！ |
回復此樓» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有7人回復
-
270調(diào)劑
已經(jīng)有7人回復
-
275學碩081000服從調(diào)劑到其他專業(yè)，保不住本專業(yè)了
已經(jīng)有6人回復
-
0710生物學，325求調(diào)劑
已經(jīng)有3人回復
-
296材料專碩求調(diào)劑
已經(jīng)有5人回復
-
385分 生物學（071000）求調(diào)劑
已經(jīng)有5人回復
-
377求調(diào)劑
已經(jīng)有3人回復
-
【求調(diào)劑】新能源材料本科，一志愿211，初試321
已經(jīng)有6人回復
-
316求調(diào)劑
已經(jīng)有11人回復
-
321求調(diào)劑
已經(jīng)有9人回復
天才超人
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 105
- 帖子: 22
- 在線: 9.3小時
- 蟲號: 1945549
- 注冊: 2012-08-19
- 專業(yè): 環(huán)境工程
4樓2017-12-17 03:39:10
2樓2017-12-16 07:09:13
3樓2017-12-16 13:36:47
天才超人
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 105
- 帖子: 22
- 在線: 9.3小時
- 蟲號: 1945549
- 注冊: 2012-08-19
- 專業(yè): 環(huán)境工程
5樓2017-12-18 06:20:18
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導師 +14 | 傾____蕭 2026-03-27 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 一志愿9初試366 本雙非求調(diào)劑 +3 | 運氣來得若有似?/a> 2026-04-02 | 3/150 |
|
|
|
[考研] 314求調(diào)劑 +9 | 1xiaojun23 2026-03-31 | 9/450 |
|
|
|
[考研] 282求調(diào)劑 +7 | 呼吸都是減肥 2026-04-01 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +10 | 陽陽哇塞 2026-04-01 | 10/500 |
|
|
|
[考研] 材料工程322分 +7 | 哈哈哈吼吼吼哈 2026-04-01 | 7/350 |
|
|
|
[考研] 286求調(diào)劑 +16 | PolarBear11 2026-03-26 | 16/800 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +13 | 大饞小子 2026-03-27 | 14/700 |
|
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學 ，080500材料科學與工程學碩，總分289分 +10 | @taotao 2026-03-29 | 10/500 |
|
|
|
[考研] 349求調(diào)劑 +6 | zwjjjjjj 2026-03-31 | 6/300 |
|
|
|
[考研] 一志愿a區(qū)211，085601-307分求調(diào)劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|
|
[考研] 求調(diào)劑 生物學 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
|
[考研] 江蘇蘇北高校誠邀調(diào)劑同學 +3 | zzll406 2026-03-31 | 3/150 |
|
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
|
[考研] 323分 食品與營養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
|
[考研] 復試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|