嗯。我有一个不同的错误消息从你的代码!我纠正它,如下所示。
使用直线法解决PDE
98(30天)
显示旧的评论
大家好！
我试图解决PDE通过直线法,使用ODE15s。
我得到错误消息
不足够的输入参数。
错误在TracerFlow(3号线)
新西兰dCdt = 0 (1);
脚本的代码粘贴。
谢谢
% RunTracerFlow
关闭所有
清晰的所有
clc
%的数据
tspan = linspace (0, 61);
L = 1;
新西兰= 100;
CA0init = 0.1;
dz = l /新西兰;
Da = 2 e 1;
U = 2 e 1;
k = 1;
IC = 0(1、新西兰);
%的能手
[t c] = ode15s (@TracerFlow tspan, IC,[],新西兰CA0init, dz,哒,U, k)
%重新计算
C (: 1) = CA0init + 1./900。* (60。* t t ^ 2);
C: N + 1) = 1. / 3。* (4。* C(:,新西兰)- C (:, Nz-1));% dCdt = 0
%绘制
新西兰xaxis = linspace (0 L + 1);
桠溪= tspan;
图像(xaxis桠溪C)
包含(的轴向位置)
ylabel (“时间间隔”)
colormap飞机
colorbar
函数dCdt = TracerFlow (t、C、新西兰、CA0init dz,哒,U, k)
%预
新西兰dCdt = 0 (1);
%定义边界条件
C (1) = CA0init + 1./900。* (60。* t t ^ 2);
C(新西兰+ 1)= 1. / 3。* (4。* C(新西兰)- C (Nz-1));% dCdt = 0
为我= 2:新西兰
dCdz (i) = 1. / (2。* dz)。* (C (i + 1) - C(张));%为中心
d2Cdz2 (i) = 1. / (dz ^ 2)。* (C (i + 1) 2。* C(我)+ C(张));
dCdt (i) = Da。* d2Cdz2 (i) - u . * dCdz (i) - k . * C (i)。^ 2;
结束
结束
接受的答案
艾伦·史蒂文斯
2020年11月3日
关键线路需要更改
[t c] = ode15s (@TracerFlow tspan, IC,[],新西兰CA0init, dz,哒,U, k)
%重新计算
C (: 1) = CA0init + 1./900。* (60。* t t ^ 2);
C: N + 1) = 1. / 3。* (4。* C(:,新西兰)- C (:, Nz-1));% dCdt = 0
他们应该
[t、C] = ode15s (@TracerFlow tspan, IC,[],新西兰CA0init, dz,哒,U, k);
%重新计算
C (: 1) = CA0init + 1./900。* (60。* t t ^ 2);
C(:,新西兰+ 1)= 1. / 3。* (4。* C(:,新西兰)- C (:, Nz-1));% dCdt = 0
同时,我更换图片浏览。
