Matlab实验7.pptVIP

* 实验七：滤波器的设计 1无限冲激响应滤波器阶数估计和滤波器设计 巴特沃兹带阻滤波器 Ws = [0.4 0.6]; Wp = [0.2 0.8]; Rp = 0.4; Rs = 50; % Estimate the Filter Order [N1, Wn1] = buttord(Wp, Ws, Rp, Rs); % Design the Filter [num,den] = butter(N1,Wn1,stop); % Display the transfer function disp(Numerator Coefficients are );disp(num); disp(Denominator Coefficients are );disp(den); % Compute the gain response [g, w] = gain(num,den); % Plot the gain response plot(w/pi,g);grid axis([0 1 -60 5]); xlabel(\omega /\pi); ylabel(Gain in dB); title(Gain Response of a Butterworth Bandstop Filter) 2吉布斯现象； 在傅立叶分布、谐波数列叠加合成时，会产生尖锐不连续的区间，好比方波的下降及上升，这就是吉布斯现象，这种现象产生的不仅与谐波的叠加合成的数量有关，而且其变化的宽度随合成的个数增加而变窄。 单位方波： 傅立叶级数表示： function y = gibbs(x,M,duty) % GIBBS % x： 待估价的值 (0x0.5) % M： 区间 % duty：占空比(0-1) if (nargin3) duty = 0.5; end p = (-M:M); for (j=1:length(x)) y(j) = sum(duty*sinc(duty*p).*exp(-i*2*pi*p*abs(x(j)))); end 正规重构：覆盖整个单位方波范围 s = linspace(-0.5,0.5,200); plot(s,gibbs(s,M)); M=5 M=25 结论：随着M的增大，摆动的振幅不会变成0。 3有限冲激响应滤波器阶数估计和滤波器设计； 用凯塞窗设计一FIR低通滤波器， 低通边界频率 阻带边界频率 阻带衰减不小于50dB。 解： 首先由过渡带宽和阻带衰减来决定凯塞窗的N和 ?凯塞窗设计对应的MATLAB程序为： wn=kaiser(30,4.55); nn=[0:1:29]; alfa=(30-1)/2; hd=sin(0.4*pi*(nn-alfa))./(pi*(nn-alfa)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-80,10]); grid; xlabel(归一化频率/?) ylabel(幅度/dB) 4设计滤波器，对于给定的一个CEG和弦，分别提取出C、E、G和弦。 提取C和弦的程序如下： clear all; close all; % Read in the sound data [d,r] = wavread(CEG.wav); % r is the sampling rate % d is the data d=d(0.3*r:0.6*r); d=d(:); % make sure its a column d=d/max(d); % normalize *