问题标题：
用窗函数设计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...
廖晓峰回答：
???Undefinedfunctionorvariable'wd'.多谢帮助，可是还有小问题啊
宋炯回答：
呵呵，不好意思，我根据你给的数据就改了前面，忘了改后面（wd应全改为wc），试试这个看，一样的（注释和title跟上面一样，可以自己添加，得出的图像你觉得效果不佳可以改axis([-4554-0.92]);对应坐标参数）wp=0.6*pi;ws=0.4*pi;wd=(wp+ws)/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');axis([-4554-0.92]);subplot(222),stem(nn,w,'k');axis([-4554-0.52]);xlabel('n');subplot(223),stem(nn,h,'k');axis([-4554-0.92]);xlabel('n');w=(-512:511)/511;subplot(224),plot(w,HH,'k');axis([-34-250250]);xlabel('omega/pi');
查看更多
