| 查看: 1316 | 回复: 12 | ||
郑美琴琴金虫 (著名写手)
|
[求助]
求助，程序中的错误时什么原因 已有1人参与
|
|
这是Command window的程序： global nr nz dr dz drs dzs... r z Dc Dt ca Tk... cae Tke h k E R... rk0 v rho Cp Tw dH... ncall %模型参数 ca0=0; cae=0.01; Tk0=305.0; Tke=305.0; Tw=355.0; r0=2.0; z1=100.0; v=1; Dc=0.1; Dt=0.1; k=0.01; h=0.01; rho=1.0; Cp=0.5; rk0=1.5e+09; dH=-10000.0; E=15000.0; R=1.987; %x轴网格 nz=20; dz=z1/nz; for i=1:nz; z(i)=i*dz; end %半径网格化 nr=7; dr=r0/(nr-1); for j=1:nr; r(j)=(j-1)*dr; end drs=dr^2; %自变量 tf=200.0; tout=[0:50:tf]'; nout=5; ncall=0; %初始条件 for i=1:nz; for j=1:nr ca(i,j)=ca0; Tk(i,j)=Tk0; y0((i-1)*nr+j)=ca(i,j); y0((i-1)*nr+j*nz*nr)=Tk(i,j); end end %ODE集成 reltol=1.0e-04;abstol=1.0e-04; options=odeset('RelTol',reltol,'AbsTol',abstol); [t,y]=ode15s(@pde_13,tout,y0,options); 这是function pde_的程序： function yt=pde_13(t,y) global nr nz dr dz drs dzs... r z Dc Dt ca Tk... cae Tke h k E R... rk0 v rho Cp Tw dH... ncall for i=1:nz for j=1:nr ij=(i-1)*nr+j; ca(i,j)=y(ij); Tk(i,j)=y(ij+nr*nz); end end for i=1:nz for j=1:nr if(j==1) car(i,j)=2*(ca(i,j+1)-ca(i,j))/drs; Tkr(i,j)=2*(Tk(i,j+1)-Tk(i,j))/drs; elseif(j==nr) car(i,j)=0.0; Tkr(i,j)=(1/r(j))*(h/k)*(Tw-Tk(i,j)); else car(i,j)=(1/r(j))*(ca(i,j+1)-ca(i,j-1))/(2*dr); Tkr(i,j)=(1/r(j))*(Tk(i,j+1)-Tk(i,j-1))/(2*dr); end if(j==1) carr(i,j)=2*(ca(i,j+1)-ca(i,j))/drs; Tkrr(i,j)=2*(Tk(i,j+1)-Tk(i,j))/drs; elseif(j==nr) carr(i,j)=2*(ca(i,j-1)-ca(i,j))/drs; Tkf(i,j)=Tk(i,j-1)+2*dr*h/k*(Tw-Tk(i,j)); Tkrr(i,j)=(Tkf-2.0*Tk(i,j)+Tk(i,j-1))/drs; else carr(i,j)=(ca(i,j+1)-2.0*ca(i,j)+ca(i,j-1))/drs; Tkrr(i,j)=(Tk(i,j+1)-2.0*Tk(i,j)+Tk(i,j-1))/drs; end if(i==1) caz(i,j)=(ca(i,j)-cae)/dz; Tkz(i,j)=(Tk(i,j)-Tke)/dz; else caz(i,j)=(ca(i,j)-ca(i-1,j))/dz; Tkz(i,j)=(Tk(i,j)-Tk(i-1,j))/dz; end rk=rk0*exp(-E/(R*Tk(i,j)))*ca(i,j)^2; cat(i,j)=Dc*(carr(i,j)+car(i,j))-v*caz(i,j)-rk; Tkt(i,j)=Dt*(Tkrr(i,j)+Tkr(i,j))-v*Tkz(i,j)-dH/(rho*Cp)*rk; end end for i=1:nz for j=1:nr ij=(i-1)*nr+j; yt(ij)=cat(i,j); yt(ij+nr*nz)=Tkt(i,h); end end yt=yt'; ncall=ncall+1; 运行结果： Warning: Divide by zero. > In pde_13 at 44 In funfun\private\odearguments at 110 In ode15s at 227 Warning: Divide by zero. > In pde_13 at 44 In funfun\private\odearguments at 110 In ode15s at 227 Warning: Divide by zero. > In pde_13 at 44 In funfun\private\odearguments at 110 In ode15s at 227 Warning: Divide by zero. > In pde_13 at 44 In funfun\private\odearguments at 110 In ode15s at 227 Warning: Divide by zero. > In pde_13 at 44 In funfun\private\odearguments at 110 In ode15s at 227 Warning: Divide by zero. > In pde_13 at 44 In funfun\private\odearguments at 110 In ode15s at 227 ??? Subscripted assignment dimension mismatch. Error in ==> pde_13 at 32 Tkrr(i,j)=(Tkf-2.0*Tk(i,j)+Tk(i,j-1))/drs; Error in ==> funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode15s at 227 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... 各位大神，帮帮忙。感激不尽！ |
回复此楼» 猜你喜欢
285求调剂
已经有12人回复
-
求调剂
已经有6人回复
-
085600材料与化工301分求调剂院校
已经有19人回复
-
277工科求调剂
已经有11人回复
-
277求调剂 数一104分
已经有13人回复
-
304求调剂
已经有10人回复
-
336求调剂，一志愿中科大
已经有6人回复
-
071000生物学，一志愿深圳大学296分，求调剂
已经有3人回复
-
一志愿北京化工085600 310分求调剂
已经有19人回复
-
274求调剂求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp 编译主程序时候出错，求好心人看一下
已经有17人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- matlab为啥出现以下这种错误
已经有6人回复
- 很简单的Fortran赋值出错啦，求大神们帮忙找bug...
已经有11人回复
- 求助：fortran程序中换行及return报错
已经有12人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求助：intel visual fortran的格式问题
已经有14人回复
- RVA测黏度中遇到问题
已经有12人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 求助！！！！安装MS5.0程序过程中出现的错误提示
已经有12人回复
- ｃ＋＋中的负数求模问题
已经有4人回复
- 在网上求助计算化学问题时的注意事项-Sobereva
已经有6人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【求助】Error termination request processed by link 9999是什么错误？程序在下面
已经有13人回复
- 【求助】运行ms的时候跳出“由于另一个程序正在运行中，此操作无法完成。请选择切换…
已经有28人回复
- 【求助】此程序也是没有问题了，但是依然运行不出结果
已经有7人回复
- 【求助】fortran程序未出错，但运行不出结果，中途中断！
已经有6人回复
- 【求助】有很多错误的程序，找不出来了
已经有18人回复
- 【求助】程序运行后出现的错误 access violation
已经有15人回复
- 【求助】HTTP500内部程序错误
已经有8人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
郑美琴琴: 金币+100, ★★★★★最佳答案 2015-04-19 16:14:21
感谢参与，应助指数 +1
郑美琴琴: 金币+100, ★★★★★最佳答案 2015-04-19 16:14:21
|
本帖内容被屏蔽 |
2楼2015-04-19 11:23:19
郑美琴琴
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 1607.8
- 散金: 84
- 红花: 1
- 帖子: 1002
- 在线: 147.4小时
- 虫号: 1765560
- 注册: 2012-04-19
- 性别: MM
- 专业: 化工系统工程
3楼2015-04-19 16:15:14
|
本帖内容被屏蔽 |
4楼2015-04-19 16:39:45
郑美琴琴
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 1607.8
- 散金: 84
- 红花: 1
- 帖子: 1002
- 在线: 147.4小时
- 虫号: 1765560
- 注册: 2012-04-19
- 性别: MM
- 专业: 化工系统工程
5楼2015-04-20 09:19:08
郑美琴琴
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 1607.8
- 散金: 84
- 红花: 1
- 帖子: 1002
- 在线: 147.4小时
- 虫号: 1765560
- 注册: 2012-04-19
- 性别: MM
- 专业: 化工系统工程
|
大神，能否再咨询个问题，上面的编程我是按照书上别人的程序编程的，有个地方看不懂，能否帮忙解答下？实在感激不尽！ function yt=pde_13(t,y)中， 一开始有个： for i=1:nz for j=1:nr ij=(i-1)*nr+j; ca(i,j)=y(ij); Tk(i,j)=y(ij+nr*nz); end 后面有： for i=1:nz for j=1:nr ij=(i-1)*nr+j; yt(ij)=cat(i,j); yt(ij+nr*nz)=Tkt(i,j); end 书上的解释分别是 % 1D to 2D matrices % 2D to 1D matrices. 请问这是什么意思？为什么要这样设置呢？  |
6楼2015-04-21 17:21:19
|
本帖内容被屏蔽 |
7楼2015-04-21 18:27:17
郑美琴琴
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 1607.8
- 散金: 84
- 红花: 1
- 帖子: 1002
- 在线: 147.4小时
- 虫号: 1765560
- 注册: 2012-04-19
- 性别: MM
- 专业: 化工系统工程
太谢谢了！ 8楼2015-04-21 19:59:57
郑美琴琴
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 1607.8
- 散金: 84
- 红花: 1
- 帖子: 1002
- 在线: 147.4小时
- 虫号: 1765560
- 注册: 2012-04-19
- 性别: MM
- 专业: 化工系统工程
9楼2015-04-22 15:34:09
|
本帖内容被屏蔽 |
10楼2015-04-22 16:29:47
