%Appl3_4
%讀入文本數(shù)據(jù)，繪制地震波時間域曲線及頻譜分析圖              
load zdsud.dat;                          %讀取數(shù)據(jù)序列
Xt=zdsud;                                %把數(shù)據(jù)賦值給變量 
Fs=50;                                   %設(shè)定采樣率 單位(Hz)
dt=1/Fs;                                 %求采樣間隔 單位(s)
N=length(Xt);                            %得到序列的長度
Xf=fft(Xt);                              %對信號進行快速Fourier變換(FFT)
figure(1)
subplot(2,1,1),plot([0:N-1]/Fs,zdsud);        %繪制原始值序列
xlabel('時間/s');
title('原始地震波記錄');
grid on;

subplot(2,1,2),plot([0:N-1]/(N*dt),abs(Xf)*2/N);%繪制信號的振幅譜
xlabel(' 頻率/Hz'),
title('幅頻圖');
ylabel('振幅');
xlim([0 Fs/2]);                          %頻率軸只畫出Nyquist頻率之前的部分
grid on;