 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 928 | 回復(fù): 3 | ||
鄭美琴琴金蟲 (著名寫手)
|
[求助]
??? Attempted to access C(2); index out of bounds because numel(C)=1. 已有2人參與
|
|
function qt=HSDMJ(t,q) global nr dr drs ... r z Ds C Cs q... v R kf K... C0 p n ncall qr qrr pp m V it y % r網(wǎng)格化 for j=1:nr % 2/r*qr if(j==1) % r=0 qr(j)=2.0*(q(j+1)-q(j))/drs; elseif(j==nr) % r=R qr(j)=2/R*kf/(pp*Ds)*(C-Cs); else qr(j)=2/r(j)*(q(j+1)-q(j-1))/(2*dr); end % qrr if(j==1) qrr(j)=2.0*(q(j+1)-q(j))/drs; elseif(j==nr) qf(j)=q(j-1)+2.0*dr*kf/(pp*Ds)*(C-Cs); qrr(j)=(qf(j)-2.0*q(j)+q(j-1))/drs; else qrr(j)=(q(j+1)-2.0*q(j)+q(j-1))/drs; end % PDEs qt(j)=Ds*(qrr(j)+qr(j)); end % 轉(zhuǎn)置和計數(shù) qt=qt'; ncall=ncall+1; clc; clear all; close all; %% 定義全局變量 global nr dr drs ... r z Ds C q... v R kf K... C0 p n ncall qr qrr F pp V x y m it y Cs %% 模型參數(shù) C0=300; % mg/L F=30; % mL/min m=23; % g d=0.02; % m R=0.0015; % m p=1100; %kg/m^3 pp=718.6; %kg/m^3 Ds=60*4.187e-12; % m^2/min K=32.948; n=2.822; kf=60*2.495e-5; % m/min V=0.025; %L %% r方向網(wǎng)格 % r的范圍為0-R nr=7; r=linspace(0, R, nr); dr=R/(nr-1); for j=1:nr r(j)=(j-1)*dr; end drs=dr^2; tf=180; tout=0:5:tf; nout=37; ncall=0; dt=tf/(nout-1); Cs(1)=0; C(1)=C0; %% 初始條件 for j=1:nr q(j)=0; q0(j)=q(j); end C=zeros(nout,1); for it=1:nout if(it==1) C(it)=C0; else C(it)=(C(it-1)+3*m*kf*dt/(R*pp*V)*Cs(it))/(1+3*m*kf*dt/(R*pp*V)); end C=C(it); Cs=Cs(it); %% ode參數(shù)設(shè)置及方程求解 reltol=1.0e-04; abstol=1.0e-04; options=odeset('RelTol',reltol,'AbsTol',abstol); [t,q]=ode15s(@HSDMJ,tout,q0,options); Cs(it+1)=(q(it,nr)/K)^n; end 報錯： ??? Attempted to access C(2); index out of bounds because numel(C)=1. Error in ==> HSDMJ_1 at 47 C(it)=(C(it-1)+3*m*kf*dt/(R*pp*V)*Cs(it))/(1+3*m*kf*dt/(R*pp*V)); |
回復(fù)此樓» 猜你喜歡
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿華南師大 070300（化學(xué)）304分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有3人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑
已經(jīng)有3人回復(fù)
-
354求調(diào)劑
已經(jīng)有8人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿西南交大，求調(diào)劑
已經(jīng)有5人回復(fù)
-
299求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專業(yè)求調(diào)劑
已經(jīng)有6人回復(fù)
-
311求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 空向量是因為我的自己的式子巧合出現(xiàn)的嗎？
已經(jīng)有6人回復(fù)
- ??? Attempted to access f(31); index out of bounds because numel(f)=30.
已經(jīng)有8人回復(fù)
- MATLAB 定義外部函數(shù)變量
已經(jīng)有14人回復(fù)
- 求助matlab問題
已經(jīng)有5人回復(fù)
材料廖
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
2樓2015-06-14 09:44:33
感謝參與，應(yīng)助指數(shù) +1
|
本帖內(nèi)容被屏蔽 |
3樓2015-06-14 09:50:58
鄭美琴琴
金蟲 (著名寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 1607.8
- 散金: 84
- 紅花: 1
- 帖子: 1002
- 在線: 147.4小時
- 蟲號: 1765560
- 注冊: 2012-04-19
- 性別: MM
- 專業(yè): 化工系統(tǒng)工程
4樓2015-06-14 15:25:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 鄭美琴琴 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +8 | 吃吃吃才有意義 2026-03-19 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人！��！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|