| 查看: 1049 | 回复: 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)); |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 空向量是因为我的自己的式子巧合出现的吗？
已经有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
- 专业: 金属材料的合金相、相变及
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
- 专业: 化工系统工程
4楼2015-06-14 15:25:05
20