fir滤波器设计.txt

压缩包中共有7个文件

设计带通滤波器
[f1,f2]=freqspace(21,'meshgrid');
Hd=ones(21);
r=sqrt(f1.^2+f2.^2);
Hd((r<0.1)|(r>0.5))=0;
colormap(jet(64));
mesh(f1,f2,Hd);

h1=fwind1(Hd,hamming(21));%一维汉明窗
figure,freqz2(h1)
h2=fwind1(Hd,bartlett(21));
figure,freqz2(h2)
h3=fwind1(Hd,hanning(21));
figure,freqz2(h3)
h4=fwind1(Hd,blackman(21));
figure,freqz2(h4)


二维高斯窗设计近似圆对称的带通滤波器
[f1,f2]=freqspace(21,'meshgrid');
Hd=ones(21);
r=sqrt(f1.^2+f2.^2);
Hd((r<0.1)|(r>0.5))=0;%通频带为0.5－0.1
colormap(jet(64));
mesh(f1,f2,Hd);
%生成旋转对称的二维高斯窗函数
win=fspecial('gaussian',21,2);
win=win./max(win(:));%将窗口的最大值归1
mesh(win)
h1=fwind2(Hd,win);%用前一步生成的窗函数设计滤波器
freqz2(h1)


二维拉普拉斯－高斯窗设计近似圆对称的高通滤波器
[f1,f2]=freqspace(21,'meshgrid');
Hd=ones(21);
r=sqrt(f1.^2+f2.^2);
Hd(r<0.5)=0;
colormap(jet(64));
mesh(f1,f2,Hd);
win=fspecial('log',21,2);
win=win./max(win(:));
mesh(win);
h2=fwind2(Hd,win);
freqz2(h2)