fir6.m

matlab离散信号时域分析

%给定通带截止频率和阻带截止频率以及阻带最小衰减， 如何用窗函数法设计线性相位低通滤波器

%数据的初始化

Wp=0.3*pi;Ws=0.5*pi;As=40;

%利用Kaiser窗函数设计
Rs=10^(-40/20);
f=[0.3 0.5];
a=[1 0];
dev=Rs*ones(1,length(a));

[M,Wc,beta,ftype]=kaiserord(f,a,dev);

%利用一型滤波器设计
M=mod( M, 2) + M;

%h=fir1(M,Wc,ftype,kaiser( M + 1 ,beta));
h = fir1(M,Wc,kaiser(M+1,beta),ftype);
omega=linspace(0,pi,512);
mag=freqz(h,[1],omega);
plot(omega/pi,20*log10(abs(mag)));
xlabel('Normalized frequency');
ylabel('Gain ,dB');
grid;