设理想带阻滤波器频率响应为 { 1 0<=w<=π/3 He(e(jw))={ 0 π<=|w|<=2π/3 { -1 2π/3<=|w|<=π 利用凯泽窗函数设计长度为45的带通滤波器,阻带衰减为60Hz。 要求:实验... 设理想带阻滤波器频率响应为 { 1 0<=w<=π/3
He(e(jw))={ 0 π<=|w|<=2π/3
{ -1 2π/3<=|w|<=π
利用凯泽窗函数设计长度为45的带通滤波器,阻带衰减为60Hz。
要求:实验基本原理,源程序及其标注。
%凯泽窗函数设计的带通滤波器幅频特性
%利用kaiser窗函数设计一个长度为奇数的带通滤波器(长度为奇数对应滤波器是偶数阶的,
% 由此能满足带通滤波器的要求),通带范围是15hertz到40hertz,阻带范围是0hertz到10
% hertz、50hertz到250hertz,阻带波纹为0.01decibel,通带波纹为0.05decibel,信号
%采样率为500hertz
clear all;
fs=500;%采样频率
fcuts=[10,15,40,50];%频率参数
mags=[0 1 0];%用于制定过渡带(通带或阻带)频率段的理想滤波器幅度值,在通带内其元
%素值取1,阻带内取0
devs=[0.01,0.05,0.01];%设置各通带和阻带内允许的幅度的Z大误差
[n,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fs);
n=n+rem(n,2);%求偶数阶,rem求余数
hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');%线性相位通带标准滤波器的设计
[H,f]=freqz(hh,1,1024,fs);%频率响应
plot(f,abs(H),'r','LineWidth',2);
grid on;
title('凯泽窗函数设计的带通滤波器幅频特性');
xlabel('频率(f)');
ylabel('幅度');
%%%%%%%%%%%%%%%%%%%%%%%设置背景色为白色%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
set(gcf,'color','w');