box.m

matla窗口制作

%___________________________常见窗函数的实现FIR滤波器_____________
%_______________________    条件：N=40，wc=0.2*pi----------------
N=14;
n=[0:1:N-1];
wc=0.2*pi;
hd=ideal_lp1(wc,N);
%____________________________矩形窗函数___________________________
w_box=(boxcar(N))';
h=hd.*w_box;                                %理想滤波器与窗函数相乘
[dp,mag,pha,grd,w]=freqz_m(h,[1]);          %求幅值响应
delta_w=2*pi/1000;
Ap=-(min(db(1:1:wc/delta_w+1)));
subplot(221);
stem(n,hd);
title('理想单位脉冲响应hd(n)');
subplot(222)
stem(n,w_box);
%stepz(n);
title('矩形窗w(n)');
axis([-5,50,0,1.2]);
subplot(223);
stem(n,h);
title('实际单位脉冲响应');
subplot(224);
plot(w/pi,dp);
title('幅度响应/db');
axis([0,1,-200,10]);
grid on;