26年CH通信工程 基于MATLAB的IIR数字滤波器设计与仿真终稿全(2588字).doc

附录第3章关键代码3.3.1 脉冲响应不变法设计低通/带通滤波器%低通Fs = 8000; fp = 1500; fs = 2000; Rp = 1; Rs = 40;Wp_lp = fp / (Fs/2); Ws_lp = fs / (Fs/2);[n_lp, ~] = buttord(Wp_lp, Ws_lp, Rp, Rs);[z, p, k] = buttap(n_lp);[B, A] = zp2tf(z, p, k);[Bt, At] = lp2lp(B, A, 2*pi*fp);[b_lp_imp, a_lp_imp] = impinvar(Bt, At, Fs);%带通fp1 = 1200; fp2 = 1800; fs1 = 800; fs2 = 2200;Wp_bp = [fp1, fp2] / (Fs/2); Ws_bp = [fs1, fs2] / (Fs/2);[n_bp, ~] = buttord(Wp_bp, Ws...