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