二维边值问题（II）波动方程、动图处理

波动方程

①源代码：

clear

t=linspace(0,100,5000);

x=linspace(0,10,500);alpha=1.2;x0=0;c=1;  %初始化

dx=x(2)-x(1);

dt=t(2)-t(1);

r=c*(dt/dx);

for j=1:2

    for i=1:length(x)

        y(i,j)=sin(i*dx-x0)*exp(-1.2*i*dx);

    end

end