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

登錄注冊

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

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

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

返回列表

天才超人

新蟲 (初入文壇)

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

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

各位大神走過路過拜托幫幫忙 🙏

我需要用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換來換去，實(shí)在是解不出來。所以來拜托各位大神了 🙏🙏  謝謝！

回復(fù)此樓

» 猜你喜歡

材料考研調(diào)劑已經(jīng)有7人回復(fù)
270調(diào)劑已經(jīng)有7人回復(fù)
275學(xué)碩081000服從調(diào)劑到其他專業(yè)，保不住本專業(yè)了已經(jīng)有6人回復(fù)
0710生物學(xué)，325求調(diào)劑已經(jīng)有3人回復(fù)
296材料專碩求調(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ù)

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的回帖

陳時(shí)通

鐵蟲 (小有名氣)

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

不會(huì) 我的問題和你一樣

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

贊一下

回復(fù)此樓

無

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

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

天才超人

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105
帖子: 22
在線: 9.3小時(shí)
蟲號(hào): 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小時(shí)
蟲號(hào): 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小時(shí)
蟲號(hào): 1945549
注冊: 2012-08-19
專業(yè): 環(huán)境工程

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

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

收到了錯(cuò)誤信息：
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小時(shí)
蟲號(hào): 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小時(shí)
蟲號(hào): 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ā)表

[考研] 311求調(diào)劑 +5	勇敢的小吳 2026-04-02	5/250	2026-04-02 09:27 by syjjj0321
[考研] 314求調(diào)劑 +9	1xiaojun23 2026-03-31	9/450	2026-04-02 09:18 by olim
[考研] 材料調(diào)劑 +12	一樣YWY 2026-04-01	12/600	2026-04-02 09:15 by olim
[考研] 310求調(diào)劑 +11	爭取九點(diǎn)睡 2026-03-30	11/550	2026-04-02 09:09 by sanrepian
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +14	睿08 2026-03-27	15/750	2026-04-02 08:46 by 花??
[考研] 336求調(diào)劑 +5	kiyy 2026-04-01	5/250	2026-04-02 07:14 by chixmc
[考研] 279求調(diào)劑 +7	莫xiao 2026-04-01	7/350	2026-04-01 22:05 by 客爾美德
[考研] 調(diào)劑 +3	好好讀書。 2026-04-01	5/250	2026-04-01 15:39 by 好好讀書。
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7	hanamiko 2026-03-30	7/350	2026-04-01 13:22 by JourneyLucky
[考研] 0856材料化工調(diào)劑總分330 +18	zhubinhao 2026-03-27	18/900	2026-04-01 09:37 by oooqiao
[考研] 254材料與化工求調(diào)劑 +3	翰冬林楠 2026-03-30	4/200	2026-03-31 17:53 by yishunmin
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考研] 353求調(diào)劑 +3	江上楓_26 2026-03-28	3/150	2026-03-31 15:53 by jp9609
[考研] 一志愿大連理工大學(xué)材料求調(diào)劑 +6	Gymno 2026-03-30	6/300	2026-03-31 07:26 by 無際的草原
[考研] 293求調(diào)劑 +3	末未mm 2026-03-30	5/250	2026-03-30 17:23 by 王保杰33
[考研] 283求調(diào)劑（080500） +14	A child 2026-03-27	14/700	2026-03-30 12:06 by 探123
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 283求調(diào)劑 +3	A child 2026-03-28	3/150	2026-03-28 15:41 by ms629
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar

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

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

» 猜你喜歡

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

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