对于一组地震波数据的fft分析

clear all
clc

%地震图形绘制
load practice.dat
x=practice;
fs=50;                                 %采样频率50Hz
dt=1/fs;
N=length（x）;
t=[0:N-1]*dt;
figure（1）;
plot（tx‘k‘）;
xlabel（‘时间/s‘）;ylabel（‘振幅‘）;
title（‘地震波形信号‘）;
grid on

%fft计算
fs1=0.002;N1=4096;                     %采样平率2Hz，点数4096
n1=0:N1-1;t1=n1/fs1;
Y1=fft（xN1）;
mag1=abs（Y1）;
f1=n1*fs1/N1;
figure（2）;
subplot（211）plot（f1（1:N1/2）mag1（1:N1/2）*2/N1‘r‘）;
xlabel（‘频率/Hz‘）;ylabel（‘幅值‘）;
title（‘fs1=2Hz N1=4096 fft       幅值曲线‘）;
grid on;
b1=angle（Y1）;
subplot（212）;plot（f1（1:N1/2）b1（1:N1/2）*2/N1）;
xlabel（‘频率/Hz‘）;ylabel（‘相位‘）;
title（‘相位曲线‘）;
grid on;

fs2=10000;N1=4096;                     %采样平率10Hz，点数4096
n2=0:N1-1;t1=n2/fs2;
Y2=fft（xN1）;
mag2=abs（Y2）;
f2=n2*fs2/N1;
figure（3）;
subplot（211）plot（f2（1:N1/2）mag2（1:N1/2）*2/N1‘r‘）;
xlabel（‘频率/Hz‘）;ylabel（‘幅值‘）;
title（‘fs1=10Hz N1=4096 fft       幅值曲线‘）;
grid on;
b2=angle（Y2）;
subplot（212