ideal_lp.m

mb-ofdm-uwb经本人修改

function [hd]=ideal_lp(wc,M)
%理想低通滤波器;hd为0到M-1之间的理想脉冲响应;wc为截止频率;M为理想滤波器的长度
wp=0.10*pi;
ws=0.20*pi;
tr_width=ws-wp;
wc=(ws+wp)/2;
M=ceil(6.6*pi/tr_width)+1;
alpha=(M-1)/2;
n=0:M-1;
m=n-alpha+eps;
hd=sin(wc*m)./(pi*m);
disp(['滤波器的长度为',num2str(M)]);
n=0:M-1;
w_han=(hanning(M))';
h=hd.* w_han;
[db,mag,pha,grd,w]=FREQZ_M(h,[1]);
delta_w=2*pi/1000;
Rp=-(min(db(1:1:wp/delta_w+1)));
disp(['实际通带波动为',num2str(Rp)]);
As=-round(max(db(ws/delta_w+1:1:501)));
disp(['最小阻带衰减为',num2str(As)]);
subplot(2,2,1);stem(n,hd);
title('理想冲激响应');
axis([0 M-1 -0.1 0.3]);
ylabel('hd(n)');
subplot(2,2,2);stem(n,w_han);
title('汉宁窗');
axis([0 M-1 0 1.1]);
ylabel('w(n)');
subplot(2,2,3);stem(n,h);
title('实际冲激响应');
axis([0 M-1 -0.1 0.3]);
xlabel('n');
ylabel('h(n)');
subplot(2,2,4);
stem(w/pi,db);
title('幅度响应(dB)');
axis([0 1 -100 10]);
xlabel('以pi为单位的频率');
ylabel('分贝数');