MATLAB 用凯撒窗函数设计线性相位低通FIR数字滤波器

MATLAB 用凯撒窗函数设计线性相位低通FIR数字滤波器

clear all;
close all;
fp=1500;fs=2500;rs=40;Fs=10000;
wp=2pifp/Fs;ws=2pifs/Fs;
Bt=ws-wp; %计算过渡带
alph=0.5842*(rs-21)^0.4+0.07886*(rs-21);
M=ceil(rs-8/2.285/Bt); %求阶数M
wc=(wp+ws)/2/pi; %计算理想高通滤波器通带截止频率
hn=fir1(M,wc,kaiser(M+1,alph)); %求hamming窗hn
n=0:M;subplot(211);
stem(n,hn,’.’,‘r’);
axis([0,30,-0.1,0.4]); title(‘h(n)波形’);
xlabel(‘n’);ylabel(‘h(n)’);
[Hw,w]=freqz(hn);
subplot(212);
%plot(wFs/20pi,20log10(abs(Hw)));
% 绘制损耗函数的曲线
w=w/pi;
plot(w,20log10(abs(Hw)));
axis([0,1,-80,0]); title(‘损耗函数曲线’);
xlabel(‘w/pi’);ylabel(‘20lg|Hg(w)|’);