数字滤波器设计实验报告.doc

文档介绍：-
. z
数字滤波器设计实验报告
实验目的
研究数字滤波器的设计思想，理解数字频域，模拟频域的关系，掌握数字系统el('频率/Hz');ylabel('振幅|H(e^jw)|');
%% *_2=cos(2*pi*150*t)
Ap=1;As=60;% 定义通带及阻带衰减
dev=[10^(-As/20),(10^(Ap/20)-1)/(10^(Ap/20)+1),10^(-As/20)];% 计算偏移量
mags=[0,1,0];% 带通
fcuts=[80,120,180,220];% 边界频率
[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数
hh2=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计
*_2=filter(hh2,1,*);% 滤波
*_2(1:ceil(N/2))=[];% 群延时N/2，删除无用信号局部
L=length(*_2);N=2^(ne*tpow2(L));Hw_2=fft(*_2,N);
figure(3);subplot(2,1,1);plot(t(1:L),*_2);
grid on;title('*_2=cos(2*pi*150*t)');*label('时间/s');
subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_2));% 查看信号频谱
grid on;title('滤波后信号*_2频谱图');*label('频率/Hz');ylabel('振幅|H(e^jw)|');
%% *_3=5*cos(2*pi*600*t)
Ap=1;As=60;% 定义通带及阻带衰减
dev=[10^(-As/20),(10^(Ap/20)-1)/(10^(Ap/20)+1)];% 计算偏移量
mags=[0,1];% 高通
fcuts=[500,550];% 边界频率
[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数
hh2=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计
-
. z
*_3=filter(hh2,1,*);% 滤波
*_3(1:ceil(N/2))=[];% 群延时N/2，删除无用信号局部
L=length(*_3);N=2^(ne*tpow2(L));Hw_3=fft(*_3,N);
figure(4);subplot(2,1,1);plot(t(1:L),*_3);
grid on;title('*_3=5*cos(2*pi*600*t)');*label('时间/s');
subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_3));% 查看信号频谱
grid on;title('滤波后信号*_3频谱图');*label