用窗函数设计FIR滤波器时怎么定义理想高通滤波器响应函数wp=0.6*pi;wr=0.4*pi;tr_width=wp-wr;N=ceil(6.2*pi/tr_width);n=0:1:N-1;wc=(wr+wp)/2;hd=ideal_hp(wc,N);
来源:学生作业帮助网 编辑:作业帮 时间:2024/05/14 05:39:25
用窗函数设计FIR滤波器时怎么定义理想高通滤波器响应函数wp=0.6*pi;wr=0.4*pi;tr_width=wp-wr;N=ceil(6.2*pi/tr_width);n=0:1:N-1;wc=(wr+wp)/2;hd=ideal_hp(wc,N);
用窗函数设计FIR滤波器时怎么定义理想高通滤波器响应函数
wp=0.6*pi;wr=0.4*pi;
tr_width=wp-wr;
N=ceil(6.2*pi/tr_width);
n=0:1:N-1;
wc=(wr+wp)/2;
hd=ideal_hp(wc,N);
用窗函数设计FIR滤波器时怎么定义理想高通滤波器响应函数wp=0.6*pi;wr=0.4*pi;tr_width=wp-wr;N=ceil(6.2*pi/tr_width);n=0:1:N-1;wc=(wr+wp)/2;hd=ideal_hp(wc,N);
给你个完整的,不知道是不是你想要的!
wp=0.6*pi;wr=0.4*pi;
wc=(wr+wp)/2;
N=33;M=(N-1)/2;
nn=-M:M;
n=nn+eps;
hd=2*((-1).^n).*sin(wd*n)./(pi*n); %理想冲击响应
w=blackman(N)'; %海明窗
h=hd.*w; %实际冲击响应
H=20*log10(abs(fft(h,1024))); %实际滤波器的分贝幅度特性
HH=[H(513:1024) H(1:512)];
subplot(221),stem(nn,hd,'k');
xlabel('n'); title('理想冲击响应');axis([-45 54 -0.9 2]);
subplot(222),stem(nn,w,'k');axis([-45 54 -0.5 2]);
title('海明窗');xlabel('n');
subplot(223),stem(nn,h,'k');
axis([-45 54 -0.9 2]);xlabel('n');title('实际冲击响应');
w=(-512:511)/511;
subplot(224),plot(w,HH,'k');
axis([-3 4 -250 250]);xlabel('\omega/\pi');title('滤波器分贝幅度特性');