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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

天才超人

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105
帖子: 22
在線: 9.3小時
蟲號: 1945549
注冊: 2012-08-19
專業(yè): 環(huán)境工程

[求助] 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

運(yùn)行后一直收到如下的錯誤信息：
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換來換去，實在是解不出來。所以來拜托各位大神了 🙏🙏  謝謝！

回復(fù)此樓

» 猜你喜歡

805600專碩材料與化工348分求調(diào)劑已經(jīng)有6人回復(fù)
重慶大學(xué)材料與化工085600，初試370+，求求調(diào)劑建議已經(jīng)有5人回復(fù)
309求調(diào)劑已經(jīng)有3人回復(fù)
一志愿北交大材料工程，總分358 已經(jīng)有4人回復(fù)
302求調(diào)劑一志愿北航070300，本科鄭大化學(xué) 已經(jīng)有11人回復(fù)
化學(xué)工程專碩324分，一志愿中國礦業(yè)大學(xué)求調(diào)劑已經(jīng)有7人回復(fù)
材料求調(diào)劑已經(jīng)有11人回復(fù)
【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有4人回復(fù)
336求調(diào)劑已經(jīng)有5人回復(fù)
材料求調(diào)劑已經(jīng)有8人回復(fù)

1樓 2017-12-16 04:57:32

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

chengle80

祝福

發(fā)自小木蟲Android客戶端

2樓2017-12-16 07:09:13

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

陳時通

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 163.5
紅花: 1
帖子: 122
在線: 10.4小時
蟲號: 3871989
注冊: 2015-05-15
性別: GG
專業(yè): 流體力學(xué)

不會我的問題和你一樣

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

無

3樓2017-12-16 13:36:47

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

天才超人

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105
帖子: 22
在線: 9.3小時
蟲號: 1945549
注冊: 2012-08-19
專業(yè): 環(huán)境工程

拜托各位高手幫忙看一下。。

贊一下

回復(fù)此樓

4樓2017-12-17 03:39:10

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

天才超人

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105
帖子: 22
在線: 9.3小時
蟲號: 1945549
注冊: 2012-08-19
專業(yè): 環(huán)境工程

拜托拜托！！感謝感謝��！

贊一下

回復(fù)此樓

5樓2017-12-18 06:20:18

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

天才超人

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105
帖子: 22
在線: 9.3小時
蟲號: 1945549
注冊: 2012-08-19
專業(yè): 環(huán)境工程

我把bvp4c換成了ode45，并簡化了最后一個方程，如下：

y0=[E_A;U_A;C3A];
xspan=[0 10^(-4)];
sol=ode45(@odefun,xspan,y0);

function [dydx]= odefun(x,y)
global F R T Z1 Z2 Z3 C1A C2A i U_A FixedCharge D3 e
C1=C1A.*exp(-Z1.*F/R/T.*(y(2)-U_A));
C2=C2A.*exp(-Z2.*F/R/T.*(y(2)-U_A));
dydx(1)=F/e*(C1-C2+10^(-8)/y(3)-y(3)+FixedCharge);
dydx(2)=-y(1);
dydx(3)=-i/Z3/D3-Z3*F/R/T*y(3)*dydx(2);
end

收到了錯誤信息：
Error using bvp4c (line 251)
Unable to solve the collocation equations -- a singular Jacobian encountered.
Error in ODEtest (line 37)
sol=bvp4c(@odefun,@odebc,solinit);

可以幫忙解釋一下這是怎么回事嗎？謝謝！

贊一下

回復(fù)此樓

6樓2017-12-18 16:08:21

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

chendequan

鐵蟲 (小有名氣)

應(yīng)助: 152 (高中生)
金幣: 918.6
紅花: 17
帖子: 262
在線: 208.6小時
蟲號: 994474
注冊: 2010-04-11
專業(yè): 結(jié)構(gòu)工程

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

感謝參與，應(yīng)助指數(shù) +1

內(nèi)容已刪除

贊一下

回復(fù)此樓

QQ：516477448，真心幫助解決MATLAB相關(guān)問題，提供詳細(xì)資料，Word文檔明確具體問題及要求，盡力而為！

7樓2017-12-18 18:22:17

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

天才超人

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105
帖子: 22
在線: 9.3小時
蟲號: 1945549
注冊: 2012-08-19
專業(yè): 環(huán)境工程

已解決

回復(fù)此樓

8樓2017-12-19 13:31:42

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主天才超人的主題更新

返回列表

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

[考研] 生物學(xué)308分求調(diào)劑（一志愿華東師大） +4	相信必會光芒萬?/a> 2026-03-31	4/200	2026-04-02 06:52 by ilovexiaobin
[考研] 385分生物學(xué)（071000）求調(diào)劑 +4	qf626 2026-04-01	4/200	2026-04-02 06:44 by ilovexiaobin
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +12	hanamiko 2026-03-27	12/600	2026-04-01 22:34 by 運(yùn)氣yunqi
[考研] 一志愿鄭大材料工程290求調(diào)劑 +18	Youth_ 2026-03-30	18/900	2026-04-01 19:23 by 無際的草原
[考研] 309求調(diào)劑 +19	誰不是少年 2026-03-29	19/950	2026-04-01 15:47 by jp9609
[考研] 0856初試324分求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 15:28 by 考研上岸呀一定?/a>
[考研] 348求調(diào)劑 +9	zzzzyk123 2026-04-01	9/450	2026-04-01 14:24 by w蟲蟲123
[考研] 一志愿同濟(jì)大學(xué)323分（080500）求調(diào)劑 +4	yikeniu 2026-04-01	4/200	2026-04-01 14:06 by asdfzly
[考研] 考研調(diào)劑 +11	Amber00 2026-03-31	11/550	2026-04-01 11:32 by wangjy2002
[考研] 314求調(diào)劑 +6	1xiaojun23 2026-03-31	6/300	2026-04-01 11:28 by 今天星期八LQ
[考研] 0856材料與化工調(diào)劑，339 +10	10213207 2026-03-31	10/500	2026-04-01 08:28 by Linda Hu
[考研] 322求調(diào)劑 +4	熹僖XX 2026-03-31	4/200	2026-04-01 08:21 by JourneyLucky
[考研] 一志愿西交大080500材料學(xué)碩349 +6	jqx1258 2026-03-31	7/350	2026-03-31 21:08 by yuq
[考研] 282求調(diào)劑不挑專業(yè) 求收留 +4	Yam. 2026-03-30	5/250	2026-03-31 14:41 by 王亮_大連醫(yī)科大
[考研] 327求調(diào)劑 +5	小卡不卡. 2026-03-29	5/250	2026-03-30 19:30 by Wang200018
[考研] 332求調(diào)劑 +6	@MZB382400 2026-03-28	6/300	2026-03-30 16:57 by 無際的草原
[考研] 310求調(diào)劑 +10	爭取九點睡 2026-03-30	10/500	2026-03-30 16:45 by ztnimte
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 調(diào)劑考研 +3	王杰一 2026-03-29	3/150	2026-03-29 08:09 by fmesaito

24小時熱門版塊排行榜

[求助] Matlab解微分方程 求助 已有1人參與

» 猜你喜歡

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

[求助] Matlab解微分方程求助已有1人參與