基于blackman窗的fir數(shù)字帶通濾波器設計(已改無錯字)
2023-07-23 01:33:06 本頁面
【正文】 le(H)。grd=grpdelay(b,a,w)。（2）定義濾波器系數(shù)向量函數(shù)function hd=ideallp(wc,M)。al=(M1)/2。n=[0:(M1)]。m=nal+eps。hd=sin(wc*m)./(pi*m)。（3）主程序：% 數(shù)字濾波器指標ws1=*pi。 wp1=*pi。ws2=*pi。 wp2=*pi。As=60。tr_width=min((wp1ws1),(wp2ws2))。M=ceil(11*pi/tr_width)+1。n=[0:1:M1]。wc1=(ws1+wp1)/2。wc2=(ws2+wp2)/2。hd=ideal_lp(wc2,M)ideal_lp(wc1,M)。% 生成blackman窗w_bla=(blackman(M))39。h=hd.*w_bla% 頻域圖像的繪制freqz(h,[1])figure(2)。subplot(2,2,1),stem(n,hd)。title(39。idael impulse response39。)axis([0 M1 ])。xlabel(39。n39。)。ylabel(39。hd(n)39。)xa=0.*n。hold onplot(n,xa,39。k39。)。hold offsubplot(2,2,2),stem(n,w_bla)。title(39。blackman window39。)axis([0 M1 0 ])。xlabel(39。n39。)。ylabel(39。w(n)39。)subplot(2,2,3),stem(n,h)。title(39。actual impulse response39。)axis([0 M1 ])。xlabel(39。n39。)。ylabel(39。h(n)39。)hold onplot(n,xa,39。k39。)。hold offn=0:600。k=n/2000。f1=2*pi*450。f2=2*pi*600。x=sin(f1*k)+sin(f2*k)。y=filter(h,[1],x)。x1=sin(f1*k)。x2=sin(f2*k)。figure(2)subplot(2,2,1)。plot(x1)。title(39。x1的圖形39。)。grid。axis([0 100*pi 5 5])。ylabel(39。x139。)。subplot(2,2,2)。plot(x2)。title(39。x2的圖形39。)。grid。axis([0 100*pi 5 5])。ylabel(39。x2)39。)。subplot(2,2,3)。plot(x)。title(39。x的圖形39。)。grid。axis([0 100*pi 5 5])。ylabel(39。x39。)。subplot(2,2,4)。plot(y)。title(39。y的圖形39。)。grid。axis([0
環(huán)評公示相關推薦
鄂ICP備17016276號-1