一维对流扩散方程matlab代码

以下是一维对流扩散方程的Matlab代码示例：

% 设置网格和时间步长
dx = 0.1; % 空间步长
dt = 0.01; % 时间步长
x = 0:dx:1; % 空间网格
t = 0:dt:1; % 时间网格

% 设置初始条件和边界条件
u0 = exp(-100*(x-0.5).^2); % 初始条件
u = u0;
u(1) = 0; % 左边界条件
u(end) = 0; % 右边界条件

% 设置参数
D = 0.1; % 扩散系数
u_max = 1; % 对流速度
CFL = u_max*dt/dx; % CFL数

% 迭代求解
for n = 2:length(t)
    u_old = u;
    for i = 2:length(x)-1
        % 计算一阶导数
        du_dx = (u_old(i+1)-u_old(i-1))/(2*dx);
        % 计算二阶导数
        d2u_dx2 = (u_old(i+1)-2*u_old(i)+u_old(i-1))/dx^2;
        % 计算更新量
        u(i) = u_old(i) - CFL*du_dx + D*dt*d2u_dx2;
    end
    % 更新边界条件
    u(1) = 0;
    u(end) = 0;
end

% 可视化结果
figure;
plot(x,u0,'--','LineWidth',2);
hold on;
plot(x,u,'LineWidth',2);
xlabel('x');
ylabel('u(x,t)');
legend('初始条件','数值解');

该代码使用有限差分方法对一维对流扩散方程进行求解，并绘制了初始条件和数值解的图像。请根据实际需求进行参数和初始条件的调整。