P_S_waves.zip_s波_waves_动画 振动_地震波_振动图像

clc,clear %纵波与横波质点运动方向和波的传播方向示意图 f = 2; T = 1/f; omega = 2*pi*f; vel = 3000; x1 = 0:50:5000; x2 = x1; nx = length(x1); k = 2*pi*f/vel; n = 40; M = moviein(n); for it=1:n t = T*it/n; for xi = 1:nx A1(xi) = 0.2*cos(k*x1(xi) - omega*t); A2(xi) = 0.0; x2(xi) = x1(xi) + 100*cos(k*x1(xi) - omega*t); end %绘制横波 plot(x1*0.001,A1 + 1.0,'o') hold on plot(x1*0.001,A1 + 1.15,'o') plot(x1*0.001,A1 + 1.85,'o') plot(x1*0.001,A1 + 0.7,'o') plot(x1(20)*0.001,A1(20) + 1.3,'ro') plot(x1(20)*0.001,A1(20) + 1.15,'ro') plot(x1(20)*0.001,A1(20) + 1,'ro') plot(x1(20)*0.001,A1(20) + 0.85,'ro') plot(x1(20)*0.001,A1(20) + 0.7,'ro') plot(x1(1)*0.001,A1(1) + 1.3,'ro') plot(x1(1)*0.001,A1(1) + 1.15,'ro') plot(x1(1)*0.001,A1(1) + 1,'ro') plot(x1(1)*0.001,A1(1) + 0.85,'ro') plot(x1(1)*0.001,A1(1) + 0.7,'ro') text(0,1.8,'S-wave') %绘制纵波 plot(x2*0.001,A2 - 1.0,'o') hold on plot(x2*0.001,A2 - 1.15,'o') plot(x2*0.001,A2 - 1.85,'o') plot(x2*0.001,A2 - 0.7,'o') plot(x2(20)*0.001,A2(20) - 1.3,'ro') plot(x2(20)*0.001,A2(20) - 1.15,'ro') plot(x2(20)*0.001,A2(20) - 1,'ro') plot(x2(20)*0.001,A2(20) - 0.85,'ro') plot(x2(20)*0.001,A2(20) - 0.7,'ro') plot(x2(1)*0.001,A2(1) - 1.3,'ro') plot(x2(1)*0.001,A2(1) - 1.15,'ro') plot(x2(1)*0.001,A2(1) - 1,'ro') plot(x2(1)*0.001,A2(1) - 0.85,'ro') plot(x2(1)*0.001,A2(1) - 0.7,'ro') text(0,-0.4,'P-wave') hold off axis([-0.5,5,-2.2,2.2]) ylabel('Ampitude') xlabel('Distance(km)') M(:,it) = getframe; end movie(M,10)