杭电_数字信号处理课程设计_实验4.rar

clear all; N=21; wp=0.578*pi; ws=0.5*pi; wn=(wp+ws)/2; b=fir1(N-1,wn/pi,hamming(N)); freqz(b,1,512); disp(b);