appl7_4.m

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

%Appl7_4
figure(1)
load fjdos0351.txt;                   %读取数据序列
Xt= fjdos0351;                   %把数据赋值给变量
Fs=50; dt=0.02;              %设定采样率(单位Hz)和采样间隔(单位s)
N=length(Xt); t=[0:N-1]*dt;      %得到序列的长度和时间序列
subplot(2,1,1),plot(t,Xt);        %绘制原始值序列
title('滤波前');
ylabel('振幅');
Fcp=1;  Fcs=2;   %设置的通带边界频率，
Wp=Fcp/(Fs/2);  Ws=Fcs/(Fs/2);              %转换为标准频率
Wn=(Wp+Ws)/2;                    %给出设计时用的边界频率
wdelta=Ws-Wp;               %过渡带宽
N=ceil(8*pi/wdelta);           %按Hamming窗求最小阶数
pa=(N-1)/2/Fs;                 %计算相位延迟
b=fir1(N,Wn);           % 设计FIR带阻滤波器
Yt=filtfilt(b,1,Xt);                  %对信号进行滤波
t=t-pa;                            %计算相位延迟
subplot(2,1,2),plot(t,Yt);           %绘制滤波后波形序列
xlim([0 300])
xlabel('时间/s');ylabel('振幅');title('滤波后信号');
figure(2)
[H,f]=freqz(b,1,128,50);
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;