boxcar15.m

是一些关于matlab的实验fir滤波器代码

wp1 = 0.35*pi; wp2 = 0.45*pi; 
ws1 = 0.3*pi; ws2 = 0.5*pi;
As =20; Ap=3;
tr_width = min((wp1-ws1),(ws2-wp2)) % 1)过渡带宽
n=[0:15]; 
wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2; %2) 截止频率2个
hd= ideal_lp(wc2,16) -ideal_lp(wc1,16); % 3)hd(n)
w_bla= (boxcar(16))';                % 4) w(n)
h = hd.* w_bla;                       % 5) h(n)
[db,mag,pha,grd,w] = freqz_m(h,[1]);  % 6) h(n)频谱计算
delta_w = 2*pi/1000;
Rp= -min(db(wp1/delta_w+1:1:wp2/delta_w))         %7) MAX Rp
As = -round(max(db(ws2/delta_w+1:1:501)))         %7) Min As
subplot(2,2,1); 
stem(n,hd); title('理想冲击响应') 
axis([0 15 -0.4 0.5]); xlabel('n'); ylabel('hd(n)')
subplot(2,2,2);
stem(n,w_bla);title('15点矩形窗') 
axis([0 15 0 1.1]);  ylabel('w(n)')
subplot(2,2,3);
stem(n,h);title('实际冲击响应')
axis([0 15 -0.4 0.5]); ylabel('h(n)')
subplot(2,2,4);
plot(w/pi,db);title('衰减幅度');
axis([0 1 -150 50]); ylabel('Decibels')
grid;