求汉宁窗设计FIR低通滤波器 程序

wp=0.2*piws=0.3*pi

wd=ws-wp %计算过渡带宽度

N=ceil(8*pi/wd)%根据汉宁窗计算所需h(n)长度去大于等于x的最小整数

wn=(0.2+0.3)*pi/2 %计算理想理想低通滤波器通带截止频率

b=fir1(N,wn/pi,hanning(N+1)) %调用fir1计算低通FIR数字滤波器

freqz(b,0.25,512)

程序的错误有两点

设计滤波器的采样频率与信号采样频率不同。在用数字信号处理时，要求参与运算的序列采样频率相同。

对于f1 f2 f3三个频率来讲，采样频率fs太低了，fs至少是f3的4倍以上，因为正弦信号按照2f3采样肯定失真

因此可以考虑将dt=0.02删掉，在程序开头写fs=500(是按照10f3算的)，dt=1/fs，并将[H,f]=freqz(b,1,512,100)中的100改写为fs。而后试试，应该是你要的结果

---