MATLAB 代码优化：模拟波的传播

本文提供一个使用 MATLAB 代码模拟波传播的示例，并对代码进行优化，使其结构和风格与原始版本有更大的区别。

代码优化后的版本

% 清除命令行窗口和工作区
clc;
clear;

% 创建一个初始矩阵并设定初始条件
field = zeros(151);
field(1, 100:120) = 3;

% 设置参数
v = -0.5;
dt = 0.1;
dx = 0.1;

% 进行数值求解
for n = 1:150
    for i = 1:150
        field(n+1, i) = field(n, i) - v * dt / dx * (field(n, i+1) - field(n, i));
    end
end

% 设置边界条件
field(:, 1) = 0;
field(:, end) = 0;

% 可视化结果
C = 0:0.05:15;
Z = zeros(150, 302);
Z(1, 200:240) = 3;
for n = 1:150
    Z(n+1, 1:300) = Z(n, 2:301);
end

figure('Name', '马天行-12334108');
writerObj = VideoWriter('MTX+HW.avi', 'MPEG-4');
open(writerObj);

x = 0:0.1:15;
for n = 1:151
    % 画出传播的波
    plot(x, field(n, :), 'LineWidth', 2);
    set(gca, 'FontSize', 15, 'FontName', 'Times New Roman');
    xlabel('x', 'FontSize', 15);
    ylabel('u(x,t)', 'FontSize', 15);
    title('马天行-12334108');
    axis([1, 15, -0.5, 5]);
    hold on;
    
    % 画出理论波形
    D = Z(n, :);
    D(:, 301) = [];
    plot(C, D, 'LineWidth', 2);
    hold off;
    
    legend('传播的波', '理论波形');
    
    frame = getframe(gcf);
    writeVideo(writerObj, frame);
end

close(writerObj);

在上述代码中，我进一步修改了变量名称、调整了代码的注释风格、修改了图形窗口的标题，以及修改了一些格式和细节。希望这次修改能够满足您的要求，使代码在结构和风格上与原始版本有更大的区别。