 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1970 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
津津愛睡覺金蟲 (小有名氣)
|
[求助]
Matlab小白，急需計(jì)算一個(gè)模型，求懂的大神看看怎么修改。萬分感謝 已有2人參與
|
||
|
Matlab小白，急需計(jì)算一個(gè)模型，解微分方程組，求懂的大神看看怎么修改。萬分感謝 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') 錯(cuò)誤使用 odearguments (line 92) FUNCT 返回的矢量的長(zhǎng)度為 2，但初始條件矢量的長(zhǎng)度為 3。FUNCT 返回的矢量和初始條件矢量的元素?cái)?shù)目必須相同。 出錯(cuò) ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... |
回復(fù)此樓» 猜你喜歡
中科大材料299求調(diào)劑
已經(jīng)有12人回復(fù)
-
274求調(diào)劑
已經(jīng)有14人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑
已經(jīng)有7人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有48人回復(fù)
-
288求調(diào)劑085600材料與化工
已經(jīng)有7人回復(fù)
-
278求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703求調(diào)劑
已經(jīng)有3人回復(fù)
-
322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷
已經(jīng)有3人回復(fù)
-
材料專碩323求調(diào)劑
已經(jīng)有4人回復(fù)
津津愛睡覺
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1100.8
- 散金: 33
- 紅花: 1
- 帖子: 145
- 在線: 60小時(shí)
- 蟲號(hào): 2029058
- 注冊(cè): 2012-09-25
- 性別: GG
- 專業(yè): 信號(hào)理論與信號(hào)處理
|
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') 出錯(cuò) ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 改了程序，現(xiàn)在錯(cuò)誤變成這樣。公式就是想解一個(gè)三個(gè)未知數(shù)的微分方程組，但是里面有未知數(shù) |
3樓2019-06-23 20:57:24
chendequan
鐵蟲 (小有名氣)
- 應(yīng)助: 152 (高中生)
- 金幣: 918.6
- 紅花: 17
- 帖子: 262
- 在線: 208.6小時(shí)
- 蟲號(hào): 994474
- 注冊(cè): 2010-04-11
- 專業(yè): 結(jié)構(gòu)工程
津津愛睡覺2樓2019-06-23 10:06:51
津津愛睡覺
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1100.8
- 散金: 33
- 紅花: 1
- 帖子: 145
- 在線: 60小時(shí)
- 蟲號(hào): 2029058
- 注冊(cè): 2012-09-25
- 性別: GG
- 專業(yè): 信號(hào)理論與信號(hào)處理
|
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') 出錯(cuò) ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 改了程序，現(xiàn)在錯(cuò)誤變成這樣。公式就是想解一個(gè)三個(gè)未知數(shù)的微分方程組，但是里面有未知數(shù) |
4樓2019-06-23 20:57:30
津津愛睡覺
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1100.8
- 散金: 33
- 紅花: 1
- 帖子: 145
- 在線: 60小時(shí)
- 蟲號(hào): 2029058
- 注冊(cè): 2012-09-25
- 性別: GG
- 專業(yè): 信號(hào)理論與信號(hào)處理
5樓2019-06-23 20:58:43
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +13 | 帆船哥 2026-03-04 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生 +10 | another12 2026-03-04 | 15/750 |
|
|
|
[考研] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3 | ipys00zm58 2026-03-06 | 4/200 |
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
|
[考研] 085600材料與化工，一志愿廣州985，求調(diào)劑 +10 | qqyyaill 2026-03-05 | 10/500 |
|
|
|
[考研] 083000，總分284，求調(diào)劑 +3 | 徐yr 2026-03-04 | 3/150 |
|
|
|
[考研] 材料與化工，291，求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
|
[考研] 江西師范大學(xué)王昆-盧章輝課題組招化學(xué)方向碩士研究生兩名 +3 | 哈哈大俠 2026-02-28 | 3/150 |
|
|
|
[考研] 281電子信息求調(diào)劑 +5 | jhtfeybgj 2026-03-02 | 9/450 |
|
|
|
[考研] 材料328求調(diào)劑 +10 | 一個(gè)蘿卜02 2026-03-03 | 10/500 |
|
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 12/600 |
|
|
|
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷 +6 | 熊境喆 2026-03-04 | 6/300 |
|
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
|
[考研]
 0703化學(xué)學(xué)碩 288分
+4
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|