appl7_1.m

MATLAB7.x数字信号处理 光盘内容

%Appl7_1
clear all;
wp=[0.032 0.2];N=320;    %通带边界频率（归一化频率）和滤波器阶数
dt=0.02;   %中国数字地震台网的采样间隔为0.02s,采样频率为50Hz
b=fir1(N,wp,hanning(N+1));           %设计FIR带通滤波器
figure(1)
[H,f]=freqz(b,1,512,1/dt);        %求出频率特性
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('频率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('频率/Hz');ylabel('相位/^o');grid on;
load ChangChun.txt   %加载地震波形记录
y=filtfilt(b,1,ChangChun); %产生零相位输出，采用filtfilt 对输入信号滤波
t=[0:length(ChangChun)-1]*dt;
figure(2)
subplot (2,1,1), plot (t,ChangChun),title ('输入信号')%绘出输入信号波形
subplot (2,1,2),plot (t,y)                  %绘出输出信号波形
title('输出信号'),xlabel('时间/s')