appl3_4.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 21 行

%Appl3_4
%读入文本数据，绘制地震波时间域曲线及频谱分析图              
load zdsud.dat;                          %读取数据序列
Xt=zdsud;                                %把数据赋值给变量 
Fs=50;                                   %设定采样率 单位(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;