appl7_6.m

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

%Appl7_6
load bssud.txt   %加载数据
x=bssud; %输入信号
n=320;    %滤波器阶数
dt=0.02;   %中国数字地震台网的采样间隔为0.02s,采样频率为50Hz
t=[0:length(x)-1]*dt;    %时间序列
f=[0 0.5/25 0.6/25 1.5/25 1.6/25 1];            %频率向量
a=[0 0 1 1 0 0];                 %振幅向量
b=firls(n,f,a);              %采用 firls 设计滤波器
[h,w]=freqz(b);                %计算其频率响应
bb=remez(n,f,a);             %采用 remez 设计滤波器
[hh,w]=freqz(bb);        
figure(1)
subplot (2,1,1),plot (t,x),title ('输入信号')           %绘出输入信号波形
subplot (2,1,2),plot(w/pi,abs(h),w/pi,abs(hh),'r:');         %绘制滤波器的幅频响应
xlabel('归一化频率');ylabel('振幅');
legend('firls','remez');                       %给出图例
%滤波
y1=filtfilt(b,1,x); %产生零相位输出，采用filtfilt 对输入信号滤波
y2=filtfilt(bb,1,x); %产生零相位输出，采用filtfilt 对输入信号滤波
t=[0:length(x)-1]*dt;
figure(2)
subplot (2,1,1),plot (t,y1,'b')                  %绘出输出信号波形
title('输出信号'),xlabel('时间/s')
legend('firls');                                 %给出图例
subplot (2,1,2),plot (t,y2,'r')                  %绘出输出信号波形
title('输出信号'),xlabel('时间/s')
legend('remez');                                 %给出图例