 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 929 | 回復: 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網格化 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 % 轉置和計數 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 %% 模型參數 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方向網格 % 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參數設置及方程求解 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)); |
回復此樓» 猜你喜歡
化學求調劑
已經有5人回復
-
一志愿中國石油大學（華東） 本科齊魯工業(yè)大學
已經有3人回復
-
332求調劑
已經有4人回復
-
265求調劑
已經有9人回復
-
328求調劑，英語六級551，有科研經歷
已經有8人回復
-
求調劑
已經有3人回復
-
278求調劑
已經有6人回復
-
307求調劑
已經有10人回復
-
279分求調劑 一志愿211
已經有12人回復
-
297求調劑
已經有9人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 空向量是因為我的自己的式子巧合出現的嗎？
已經有6人回復
- ??? Attempted to access f(31); index out of bounds because numel(f)=30.
已經有8人回復
- MATLAB 定義外部函數變量
已經有14人回復
- 求助matlab問題
已經有5人回復
材料廖
木蟲 (正式寫手)
- 應助: 15 (小學生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
2樓2015-06-14 09:44:33
感謝參與，應助指數 +1
|
本帖內容被屏蔽 |
3樓2015-06-14 09:50:58
鄭美琴琴
金蟲 (著名寫手)
- 應助: 3 (幼兒園)
- 金幣: 1607.8
- 散金: 84
- 紅花: 1
- 帖子: 1002
- 在線: 147.4小時
- 蟲號: 1765560
- 注冊: 2012-04-19
- 性別: MM
- 專業(yè): 化工系統(tǒng)工程
4樓2015-06-14 15:25:05
相關版塊跳轉
我要訂閱樓主 鄭美琴琴 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 324分 085600材料化工求調劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩，總分308求調劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 295求調劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 北科281學碩材料求調劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 招收調劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 材料學碩318求調劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 材料與化工求調劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
|
[考研] 0817調劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 311求調劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|