myfir.m

fft滤波器

function myfir(wc,N)
%加矩形窗
%wc为理想低通的截止频率，N为窗函数的窗口大小
N1=100;
N2=N1*wc/pi;
%线性相位低通滤波器
Hd1=zeros(1,N1+1);
stepk=pi/N1;
w=0:stepk:pi;
for k=1:N2
    Hd1(k)=exp(-j*N*w(k)/2);            
end
%低通滤波器幅频特性
aHd1=abs(Hd1);    
%低通滤波器冲激响应
hd1=zeros(1,N1+1);
for k=1:(N1+1)
    y=@(w) exp(-j*N*w/2+j*w*k)/(2*pi);
    hd1(k)=quad(y,-wc,wc);             
end
%作图
n=1:(N1+1);
subplot(2,2,1)
stem(n,aHd1,'.');
title('低通滤波器幅频特性');
subplot(2,2,2)
stem(n,hd1,'.');
title('低通滤波器冲激响应');

%矩形窗函数
x=zeros(1,N1+1);
for k=1:N
    x(k)=1;                            
end
%加窗
h=hd1.*x;  
%求fir滤波器的频率响应
a=1;
[H,w]=freqz(h,a,256,'whole');  
%求fir滤波器的幅频特性
aH=abs(H); 
%作图
n1=length(w);
m=0:0.5*2/n1:0.5-0.5*2/n1;
m1=length(m);
subplot(2,2,3)
stem(n,h,'.');
title('fir滤波器的幅频特性');
subplot(2,2,4)
stem(m,aH(1:n1/2),'.');
title('fir滤波器的幅频特性');