第七章有限冲激响应滤波器的设计.pptVIP

2021/3/23 * 其他FIR数字滤波器的文件 remez.m:设计采用切比雪夫最佳一致逼近FIR数字 滤波器 remexord.m:采用切比雪夫一致逼近设计FIR数字滤 波器所需要的滤波器阶数 sgolay.m设计savitzky-golay平滑滤波器 firls.m:用最小平方法设计线性相位FIR数字滤波器 fircls.m:用带约束的最小平方法设计线性相位FIR数 字滤波器 Fircls1.m:用带约束的最小平方法设计线性相位FIR 低通和高通滤波器 Firrcos.m:用来设计低通线性相位FIR数字滤波器 2021/3/23 * Example of a filter % Example of a length 31 lowpass filter: h=firpm(30,[0 .1 .2 .5]*2,[1 1 0 0]); fvtool(h) % Example of a low-pass differentiator: h=firpm(44,[0 .3 .4 1],[0 .2 0 0],differentiator); % Example of a type 4 highpass filter: h=firpm(25,[0 .4 .5 1],[0 0 1 1],h); 2021/3/23 * Example--Design a Lowpass Filter With pmin = 4 and pmax = 12 b=firlpnorm(30,[0 0.4 0.45 1],[0 0.4 0.45 1],[1 1 0 0],[1 1 10 10],[4 12]); [H,W,S]=freqz(b,1,1024); S.plot = mag; fvtool(b); 2021/3/23 * Example--Return a Minimum Phase Bandstop Filter b=firlpnorm(21, [0 .25 .35 .7 .8 1], [0 .25 .35 .7 .8 1],[1 1 0 0 1 1], [1 1 5 5 1 1], [2 4], minphase); fvtool(b) 2021/3/23 * sgolay.m用来设计Savitzky-Golay平滑滤波器。其调用格式是 b=sgolay（k,f） 式中k是多项式的阶次,f是拟合的双边点数。要求 k＜f ,且f为奇数。 N = 4;F = 21;b=sgolay(N,F); （4） Savitzky-Golay平滑滤波器 2021/3/23 * fircls.m用带约束的最小平方法设计线性相位FIR数字滤波器。可设计任意给定的理想幅频特性。 n = 51; f = [0 0.4 0.8 1]; a = [0 1 0]; up = [ 0.02 1.02 0.01]; lo = [-0.02 0.98 -0.01]; b = fircls(n,f,a,up,lo); 线性相位FIR数字滤波器 2021/3/23 * 例7.3 令N=10,分别用矩形窗和海明窗重复例7.1。 解 根据要求编制MATLAB程序如下: clear all; N=10; b1=fir1(N,0.25,boxcar(N+1)); b2=fir1(N,0.25,hamming(N+1)); M=128; 2021/3/23 * h1=freqz(b1,1,M); h2=freqz(b2,1,M); % 分别求两个滤波器的频率响应; t=0:10; subplot(221) stem(t,b2,.);hold on; plot(t,zeros(1,11));grid; f=0:0.5/M:0.5-0.5/M; M1=M/4; 2021/3/23 * 程序 for k=1:M1 hd(k)=1; hd(k+M1)=0; hd(k+2*M1)=0; hd(k+3*M1)=0; end subplot(222) plot(f,abs(h1),b-,f,abs(h2),g-,f,hd,-);grid; ? 2021/3/23 * 运行结果 2021/3/23 * 例7.4设计一多带滤波器,要求理想幅频响应在归一化频率0.2~0.3,0.6~0.8之间