waveshow.m

来自「Lagrange插值、Newton插值函数等插值函数」· M 代码 · 共 39 行

function waveshow(n,velocity)

%首先定义一个点随时间的变化；

t=[0:.05:n];%设定时间为一秒；
T=1;%周期定为一秒;
omige=2*pi/T;
%om=4*pi/T;
A=100;
u=A*sin(omige*t);
%u1=200*sin(om*t);
%u2=800*sin(om*4*t);  
%定义波在空间的传播;
%u=u+u1+u2;
v=velocity;
x=v*t;
nu=zeros(1,size(t,2));
for i=1:size(t,2)
    for j=1:i
        nu(j,i)=u(i-j+1);
    end
end


%动态显示波形；

for i=1:size(t,2)
    for j=1:i
        plot(x(1:j),nu(1:j,i),'*b')
        title(['the waveshape at time',num2str(x(1,j))])
        grid on
        axis on
        axis([0 max(x) min(u) max(u)])
        xlabel('time')
        ylabel('wave-shape')
    end
    pause(0)
end