求助!谁能帮我看看这程序问题出在哪里?

daisy041216

我要设计一个低通fir数字滤波器,要求用汉宁窗,通带截止频率Wp=0.75л,阻带截止频率Ws=1л，通带衰减不大于2dB，阻带衰减不小于40dB，N=42,采样频率f=16000Hz.只要观察频谱图就可以了.
          这是我写的程序,但是出来的图不满足要求的设计指标,这是为什么呢?

wp=0.75*pi;ws=1*pi;fs=16000 %技术指标
deltaw=ws-wp; %计算过渡带△ω
N=42; %滤波器长度
n=[0:1:N-1]; %滤波器阶数
wc=(ws+wp)/2; %截止频率
hd=ideallp(wc,N); %求理想脉衝响应
w_han=(hanning(N))'; %求窗函数
h=hd.*w_han;
[H,w]=freqz(h,[1]); %求频率特性
mag=abs(H); %取幅频特性
db=20*log10((mag+eps)/max(mag)); %取对数并归一
plot(w/pi,db); %画幅频图
title('幅度响应');
grid;
axis([0 1,-100 10]);
ylabel('db');xlabel('频率')


function hd=ideallp(wc,N); %定义ideallp
alpha=(N-1)/2;
n=[0:N-1];
m=[n-alpha]+eps; %加小数避免除零操作
hd=sin(wc*m)./(pi*m); %hd=0到N-1之间的理想脉衝响应

     这是运行后的图,但是就是不满足要求!

还有阻带截止频率要从哪里开始看?

怎么计算滤波器的实际阻带波动和通带波动?

滤波器的系数要怎么求并显示出来呢?