參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  glmhlhyx 2017-12-15 13:53:33

  Matlab中，函數(shù)fir1（）和fir2（）利用加窗傅里葉級數(shù)法設(shè)計FIR濾波器。函數(shù)fir1（）用來設(shè)計傳統(tǒng)的LP（低通）、HP（高通）、BP（帶通）、BS（帶阻）和多頻帶FIR濾波器；而函數(shù)fir2（）用來設(shè)計具有任意幅度響應(yīng)的的FIR濾波器。 一、函數(shù)fir1的各種形式如下： b = fir1(N，Wn) b = fir1(N，Wn，'ftype') b = fir1(N，Wn，window) b = fir1(N，Wn，'ftype'，window) b = fir1(...，'normalization') 其中，1、N為濾波器節(jié)點個數(shù)； 2、Wn（0<Wn<1）為歸一化截止頻率； 3、ftype為濾波器類型（默認(rèn)（缺省時）是low）： 'high' for a highpass filter with cutoff frequency Wn. 'stop' for a bandstop filter， if Wn = [w1 w2]. The stopband frequency range is specified by this interval. 'DC-1' to make the first band of a multiband filter a passband. 'DC-0' to make the first band of a multiband filter a stopband. 4、window為所加窗的類型（默認(rèn)的是hamming窗），e.g hamming、chebwin、blackman、hanning、kaiser等。 5、normalization為是否將濾波器的幅度進(jìn)行歸一化： f is a vector of frequency points in the range from 0 to 1， where 1 corresponds to the Nyquist frequency.The frequency points must be in increasing order. m is a vector containing the desired magnitude response at the points specified in f. Duplicate frequency points are allowed， corresponding to steps in the frequency response. 'scale' (default): Normalize the filter so that the magnitude response of the filter at the center frequency of the passband is 0 dB. 'noscale': Do not normalize the filter. fir1（）舉例：設(shè)計一個低通濾波器 t = 0:0.01:2; f =2*sin(2*pi*20*t)+4*sin(2*pi*60*t); N = 11; %濾波器節(jié)點個數(shù) wc = 0.5; %歸一化截止頻率 hd = fir1(N，wc，'low'); % 基于加窗函數(shù)的FIR濾波器設(shè)計 ft = conv(f，hd); figure(1) plot(abs(fft(f))); title('原始信號f'); figure(2) plot(abs(fft(ft))); title('濾波后信號ft'); 二、函數(shù)fir2的各種形式如下： b = fir2(n，f，m) b = fir2(n，f，m，window) b = fir2(n，f，m，npt) b = fir2(n，f，m，npt，window) b = fir2(n，f，m，npt，lap) b = fir2(n，f，m，npt，lap，window) 其中，向量f是指定頻率點的幅度響應(yīng)樣本，與m定義的幅度響應(yīng)樣本對應(yīng)；f和m具有相同的長度，并且f的diyi個和Z后一個分量分別是0和1；可以對f中的頻點進(jìn)行復(fù)制，從而跳變地逼近幅度響應(yīng)指標(biāo)。 npt指定了函數(shù)fir2（）進(jìn)行內(nèi)插得頻率響應(yīng)的柵格點數(shù)目，默認(rèn)值為512。 lap指定了在f中重復(fù)頻率點間插入的區(qū)域大小。 舉例：設(shè)計一個30階的低通濾波器 f = [0 0.6 0.6 1]; m = [1 1 0 0]; b = fir2(30，f，m); [h，w] = freqz(b，1，128); plot(f，m，w/pi，abs(h)) legend('Ideal'，'fir2 Designed') title('Comparison of Frequency Response Magnitudes')

  贊(18)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論