 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1977 | 回復: 9 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
津津愛睡覺金蟲 (小有名氣)
|
[求助]
Matlab小白，急需計算一個模型，求懂的大神看看怎么修改。萬分感謝 已有2人參與
|
|
|
Matlab小白，急需計算一個模型，解微分方程組，求懂的大神看看怎么修改。萬分感謝 function dy = funct(t,y) I = 50 % mA F = 96487 % C/mol c0 = 10 % mM Vol = 3.344 % ml Q = 6.5 % ml/min A = 33.75 % cm2 syms x eqn = 0.86*x+0.4*(x*exp(-asinh(theta/(-2*x)))+x*exp(asinh(theta/(-2*x))))==y(3) cma = double(solve(eqn,x)) cmion = (y(3)-0.86*cma)/0.4 phi_sp = asinh(-3/y(2)) phi_ele = asinh(-3/cma) phi_mem = (I*170)/(-1.12*F*(cma.*cosh(phi_ele)+y(2).*cosh(phi_sp))) J = -(1.12*1e-3/170)*(2*cma.*cosh(phi_ele)-2*y(2).*cosh(phi_sp))+3*phi_mem dy=zeros(3,1); dy(1) = (1e4*I/A)/(116*F) dy(2) = -J/(250e-6*0.5)+(c0-y(2))/(0.5*(Vol/Q)*60) dy(3) = J/(290*1e-6) end clear all; c0=10 tspan=[0 600] [t,y] = ode45(@funct,tspan,[0 c0 1.66*c0]) plot(t, y(:,2)) xlabel('t'),ylabel('c') 錯誤使用 odearguments (line 92) FUNCT 返回的矢量的長度為 2，但初始條件矢量的長度為 3。FUNCT 返回的矢量和初始條件矢量的元素數目必須相同。 出錯 ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... |
回復此樓» 猜你喜歡
新疆大學地質與礦業(yè)工程學院招生
已經有32人回復
-
0856材料與化工353分求調劑
已經有4人回復
-
282分材料專業(yè)求調劑院校
已經有13人回復
-
315食品工程求調劑
已經有5人回復
-
環(huán)境調劑
已經有3人回復
-
接受26屆調劑生
已經有20人回復
-
考研材料與化工，求調劑
已經有3人回復
-
歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組
已經有35人回復
-
304求調劑（085602一志愿985）
已經有4人回復
-
294 英二數二物化 求調劑
已經有6人回復
hzlhm
至尊木蟲 (著名寫手)
- 應助: 387 (碩士)
- 金幣: 17808.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時
- 蟲號: 1517335
- 注冊: 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動力系統(tǒng)
7樓2019-06-24 13:20:05
chendequan
鐵蟲 (小有名氣)
- 應助: 152 (高中生)
- 金幣: 918.6
- 紅花: 17
- 帖子: 262
- 在線: 208.6小時
- 蟲號: 994474
- 注冊: 2010-04-11
- 專業(yè): 結構工程
津津愛睡覺2樓2019-06-23 10:06:51
津津愛睡覺
金蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1100.8
- 散金: 33
- 紅花: 1
- 帖子: 145
- 在線: 60小時
- 蟲號: 2029058
- 注冊: 2012-09-25
- 性別: GG
- 專業(yè): 信號理論與信號處理
|
function dydt = funct(t,x, y) I = 50; % mA F = 96487; % C/mol c0 = 10; % mM Vol = 3.344; % ml Q = 6.5; % ml/min A = 33.75; % cm2 y1=y(1); y2=y(2); y3=y(3); y3 = 0.86*x+0.4*(x*exp(-asinh(y1/(-2*x)))+x*exp(asinh(y1/(-2*x)))) phi_sp = asinh(-3/y2) phi_ele = asinh(-3/x) phi_mem=(I*170)/(-1.12*F*(x*cosh(phi_ele)+y2*cosh(phi_sp))) J =-(1.12*1e-3/170)*(2*x*cosh(phi_ele)-2*y2*cosh(phi_sp))+3*phi_mem dy1dt = (1e4*I/A)/(116*F) dy2dt = -J/(250e-6*0.5)+(c0-y2)/(0.5*(Vol/Q)*60) dy3dt = J/(290*1e-6) end clear all; I = 50; % mA F = 96487; % C/mol c0 = 10; % mM Vol = 3.344; % ml Q = 6.5; % ml/min A = 33.75; % cm2 tspan=[0 600] x0=c0 [t,x,y] = ode45(@funct,tspan, c0, [0 c0 1.66*c0]) plot(t, theta) xlabel('t'),ylabel('c') 出錯 ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 改了程序，現(xiàn)在錯誤變成這樣。公式就是想解一個三個未知數的微分方程組，但是里面有未知數 |
3樓2019-06-23 20:57:24
津津愛睡覺
金蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1100.8
- 散金: 33
- 紅花: 1
- 帖子: 145
- 在線: 60小時
- 蟲號: 2029058
- 注冊: 2012-09-25
- 性別: GG
- 專業(yè): 信號理論與信號處理
|
function dydt = funct(t,x, y) I = 50; % mA F = 96487; % C/mol c0 = 10; % mM Vol = 3.344; % ml Q = 6.5; % ml/min A = 33.75; % cm2 y1=y(1); y2=y(2); y3=y(3); y3 = 0.86*x+0.4*(x*exp(-asinh(y1/(-2*x)))+x*exp(asinh(y1/(-2*x)))) phi_sp = asinh(-3/y2) phi_ele = asinh(-3/x) phi_mem=(I*170)/(-1.12*F*(x*cosh(phi_ele)+y2*cosh(phi_sp))) J =-(1.12*1e-3/170)*(2*x*cosh(phi_ele)-2*y2*cosh(phi_sp))+3*phi_mem dy1dt = (1e4*I/A)/(116*F) dy2dt = -J/(250e-6*0.5)+(c0-y2)/(0.5*(Vol/Q)*60) dy3dt = J/(290*1e-6) end clear all; I = 50; % mA F = 96487; % C/mol c0 = 10; % mM Vol = 3.344; % ml Q = 6.5; % ml/min A = 33.75; % cm2 tspan=[0 600] x0=c0 [t,x,y] = ode45(@funct,tspan, c0, [0 c0 1.66*c0]) plot(t, theta) xlabel('t'),ylabel('c') 出錯 ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 改了程序，現(xiàn)在錯誤變成這樣。公式就是想解一個三個未知數的微分方程組，但是里面有未知數 |
4樓2019-06-23 20:57:30
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 標題:撿漏預警|08工科/09農學調劑!英語要求低，過線即有機會! +4 | 馬超放煙花 2026-03-07 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 復試調劑 +5 | 呼呼？~+123456 2026-03-08 | 6/300 |
|
|
|
[考研] 一志愿蘇大材料工程專碩293求調劑 +5 | 鋼鐵大炮 2026-03-04 | 5/250 |
|
|
[考研] 293一志愿華東理工 0817化學工程與技術 調劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
|
[考研] 求0856調劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
|
[考研] 083000環(huán)境科學與工程調劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
|
[考研] 材料化工求調劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 第一志愿上海大學，專業(yè)化學工程與技術，總分288，求調劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 材料專碩290求調劑 +8 | 杰尼龜aaa 2026-03-04 | 8/400 |
|
|
|
[考研] 289求調劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
|
[考研] 080500材料科學與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調劑 +6 | 好運好運接接 2026-03-04 | 6/300 |
|
|
|
[考研] 0856材料專碩274能調劑去哪里？ +3 | 22735 2026-03-04 | 4/200 |
|
|
|
[考研] 材料學碩080500復試調劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
|
[考研] 286 +6 | ksncj 2026-03-04 | 6/300 |
|
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
|
[考研] 281求調劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 0703化學306調劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
|
[考研] 化學0703求調劑 學碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|