firlp.asv

滤波器的MATLAB实现程序 内附多个源代码

clc
clear all
close all
T1=0;
T2=0.5925;
N=60;
alpha=(N-1)/2;
l=0:N-1;wl=(2*pi/N)*l;
Hrs=[ones(1,7),T2,T1,zeros(1,43),T1,T2,ones(1,6)];
Hdr=[1 1 0 0];
wdl=[0 0.2 0.3 1];
k1=0:floor((N-1)/2);k2=floor((N-1)/2)+1:N-1;
angH=[-alpha*(2*pi)/N*(k1),alpha*(2*pi)/N*(N-k2)];
Hdk=Hrs.*exp(j*angH);
h=ifft(Hdk,N);
h2=fft(h,1024);
figure
plot(l,Hrs,'.',wdl,Hdr);
title('lowpass')
grid
hold on
plot([0:1023]/1024,abs(h2),'k');
title('Magnitude response')
grid
axis([0 0.5 -0.1 1.1]);