山大数字信号处理课件13第八章复习总结.pdfVIP

第8章复习 FIR滤波器MATLAB 函数设计 JMY © Copyrights Reserved, SDU 山东大学数字信号处理精品课程1 / 43 • MATLAB 提供的窗函数 – W=boxcar(M) returns the M-point rectangular window function in array W. – W=triang(M) returns the M-point Bartlett window function in array W. – W=hanning(M) returns the M-point Hanning window function in array W. – W=hamming(M) returns the M-point Hamming window function in array W. – W=blackman(M) returns the M-point Blackman window function in array W. – W=kaiser(M,beta) returns the beta-valued the M-point rectangular window function in array W. JMY © Copyrights Reserved, SDU 山东大学数字信号处理精品课程2 / 43 • 用这些窗函数设计FIR 滤波器还需理想低通冲激响应 function hd = ideal_lp(wc,M); % ideal lowpass filter computation % [hd] = ideal_lp(wc,M) % hd = ideal impulse response between 0 to M-1 % wc = cutoff frequency in radians % M = length of the ideal filter alpha = (M-1)/2; n = [0:1:(M-1)]; m = n - alpha + eps; hd = sin(wc*m) ./ (pi*m); JMY © Copyrights Reserved, SDU 山东大学数字信号处理精品课程3 / 43 freqz_m 函数给出幅度绝对值响应, dB 函数, 相位响应,群延迟响应. • function [db,mag,pha,w] = freqz_m(b,a); • % modified version of freqz subroutine • % [db,mag,pha,w] = freqz_m(b,a); • % db = relative magnitude in dB computed over 0 to pi radians • % mag = absolute magnitude computed over 0 to pi radians • % pha = phase response in radians over 0 to pi radians • % w = 501frequency samples between 0 to pi radians • % b = numerator polynomial of H(z) (for FIR:b=h) • % a =denominator polynomial of H(z) (for FIR:b=h) [H,w] = freqz(b,a,1000,whole); H = (H(1:1:501)); w = (w(1:1:501)); mag = abs(H); db = 20*log10((mag+eps)/max(mag)); pha = angle(H); grd = grpdelay(b,a,w); JMY © Copyright