appl7_5.m

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

%Appl7_5
clear all;
wp=2/25;N=400;    %通带边界频率（归一化频率）和滤波器阶数
dt=0.02;   %中国数字地震台网的采样间隔为0.02s,采样频率为50Hz
b=fir1(N,wp,hamming(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 cz02.txt   %加载数据
x=cz02'; %输入信号
y=filtfilt(b,1,x); %产生零相位输出，采用filtfilt 对输入信号滤波
t=[0:length(x)-1]*dt;
figure(2)
subplot (2,1,1), plot (t,x),title ('输入信号')%绘出输入信号波形
subplot (2,1,2),plot (t,y)                  %绘出输出信号波形
title('输出信号'),xlabel('时间/s')