设计线性相位高通FIR滤波器

调用MATLAB工具箱函数frl设计线性相位高通FIR滤波器。要求通带截止频奉为0.6Πrad,限带截止频率为0.45Πrad,通带最大衰成为0.2 dB.阻带最小衰减为45 dB.显示所设计的单位脉冲响应h(n)的数据，并画出损耗函数曲线。
 

代码

N=46;
Wc=3*pi/10;
wc=Wc/pi;%频率归一化
h=fir1(N,wc, 'high', hanning(N+1));
[H,m]=freqz(h,[1],1024,'whole'); %频率响应                 
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
subplot(2,2,1)
n=0:N;
stem(n,h,'.')
axis([0 N -0.1 0.3])
hold on
n=0:N-1;
x=zeros(N);
plot(n,x,'-')
hold off
xlabel('n')
ylabel('h(n)')
title('h(n)')
subplot(2,2,2)
x=abs(H);
plot(m/pi*200,20*log10(x/max(x)));
grid on;
xlabel('f/Hz');
ylabel('幅度(dB)')
axis([0,500,-80,10]);
title('损耗函数曲线');
 

截图