Matlab编程实现数字滤波器怎么设计

你好好查点资料，这个有太多东西可以写，我给你个简单的例子：
Wc=0.5*pi;
N=21;
window=hanning(N+1); %选择要用的窗函数
b=fir1(N,Wc/pi,window); %加窗的线性相位设计函数
[H,W]=freqz(b,1,512); %求幅频响应
subplot(3,1,1)
plot(W/pi,20*log10(abs(H))); %画图增益
xlabel('Frequency(piHz)');
ylabel('Hanning');
window1=hamming(N+1);
b1=fir1(N,Wc/pi,window1);
[H,W]=freqz(b1,1,512); %求幅频响应
subplot(3,1,2)
plot(W/pi,20*log10(abs(H))); %画图增益
xlabel('Frequency(piHz)');
ylabel('Hamming');
window2=blackman(N+1);
b2=fir1(N,Wc/pi,window2);
[H,W]=freqz(b2,1,512); %求幅频响应
subplot(3,1,3)
plot(W/pi,20*log10(abs(H))); %画图增益
xlabel('Frequency(piHz)');
ylabe