数字信号处理课程设计

Fs=8192; rs=53; fp=Fs/4; fs=fp+0.04*Fs; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; Bt=ws-wp; alph=0.5842*(rs-21)^0.4+0.07886*(rs-21); N=ceil((rs-8)/2.285/Bt); b=fir1(N,Fs/4,alph); wc=(wp+ws)/2/pi; hn=fir1(N,wc,kaiser(N+1,alph)); subplot(211); stem(hn,'.');grid; subplot(212); plot(b,1,2048);grid;