按要求编写matlab代码

为了实现您所提到的要求，即针对一维平流方程的不同情况进行数值模拟并绘制图形，以下是基于MATLAB语言的一个基本示例框架。请注意，由于您的请求没有提供具体的边界条件以及完全明确的时间步进方法，下面的代码将采用一种通用的形式来解决问题的一部分：

% 定义参数
L = 10; % 假设域长度为10个单位
T = 10; % 总时间为10秒
nx = 100; % 空间网格点数
nt = 100; % 时间步长数
dx = L/(nx-1); % 空间步长
dt_values = [0.5, 0.25]; % 不同的时间步长值
c_values = [1, 3, -1]; % 平流速度

% 初始化矩阵存储解
U_central_diff = zeros(nx, nt);
U_forward_backward_diff = U_central_diff;

for dt_index = 1:length(dt_values)
    dt = dt_values(dt_index);
    
    for c_index = 1:length(c_values)
        c = c_values(c_index);
        
        % 初始条件
        u = @(x) (abs(x)<=2)*1;
        U_central_diff(:, 1) = arrayfun(u, linspace(0,L,nx));
        U_forward_backward_diff(:, 1) = U_central_diff(:, 1);

        % 主循环
        for n = 1:nt-1
            U_central_diff(2:nx-1, n+1) = ... 
                % 中央差分格式计算此处...
                
            U_forward_backward_diff(2:nx-1, n+1) = ...
                % 前后差分格式计算此处...
        end
        
        % 绘制结果
        figure;
        subplot(2,1,1);
        plot(linspace(0,L,nx), U_central_diff(:,end));
        title('Central Difference Scheme');
        
        subplot(2,1,2);
        plot(linspace(0,L,nx), U_forward_backward_diff(:,end));
        title('Forward Backward Difference Scheme');
    end
end

注意上述代码中省略号部分需要填充具体的时间离散化公式（如欧拉法等），这取决于实际的问题描述。同时，在`u`函数定义处也需要依据题目给定的具体初值条件来修改。最后，根据题目需求，还需要添加更多的逻辑来处理不同的场景及绘制相应的图像或对比分析。