实验3 第6、7章.docVIP

编程练习 1、基于巴特沃斯模拟滤波器原型使用冲激不变转换方法设计数字滤波器，滤波器采样率为1000Hz要求具有下面的参数指标： 通带截止频率： 通带波动值： 阻带截止频率： 阻带波动值： 2、要求采用双线性变换方法设计一个数字滤波器，滤波器采样率为1000Hz其他指标参数不变： 通带截止频率： 通带波动值： 阻带截止频率： 阻带波动值： 用汉宁窗和布莱克曼窗设计FIR低通滤波器。设N=11，ωc=0.2π。 在MATLAB中，先编写直接型到并联型的转换子程序tf2par.m文件和理想低通滤波器的脉冲响应子程序ideallp.m文件以及cplxcomp.m文件以备调用 调用的tf2par.m文件： function [C,B,A] = tf2par(b,a); %直接型到并联型的转换子程序 M = length(b); N = length(a); [r1,p1,C] = residuez(b,a); p = cplxpair(p1,1e-9); I = cplxcomp(p1,p); r = r1(I); K = floor(N/2); B = zeros(K,2); A = zeros(K,3); if K*2 == N; for i=1:2:N-2 pi = p(i:i+1,:); ri = r(i:i+1,:); [Bi,Ai] = residuez(ri,pi,[]); B(fix((i+1)/2),:) = real(Bi); A(fix((i+1)/2),:) = real(Ai); end [Bi,Ai] = residuez(r(N-1),p(N-1),[]); B(K,:) = [real(Bi) 0]; A(K,:) = [real(Ai) 0]; else for i=1:2:N-1 pi = p(i:i+1,:); ri = r(i:i+1,:); [Bi,Ai] = residuez(ri,pi,[]); B(fix((i+1)/2),:) = real(Bi); A(fix((i+1)/2),:) = real(Ai); end end ... 调用的ideallp.m文件： function hd = ideallp(wc,N); % 理想低通滤波器的脉冲响应子程序 % hd = ideallp(wc,N) % ---------------------------------------------- % hd = 点0到N-1之间的理想脉冲响应 % wc = 截止频率 % N = 理想滤波器的长度 tao = (N-1)/2; % 理想脉冲响应的对称中心 n = [0: (N-1)]; % 设定脉冲响应长度 m = n - tao + eps; % 加一个小数以避免零作除数 hd = sin(wc*m) ./ (pi*m); % 理想脉冲响应 调用的cplxcomp.m文件: function I = cplxcomp(p1,p2) % 比较两个包含同样标量元素，但（可能）有不同下标的复数对 I=[]; for j=1:1:length(p2) for i=1:1:length(p1) if (abs(p1(i)-p2(j)) 0.0001) I=[I,i]; end end end I=I; ... 解： MATLAB的源程序为： wp = 0.2*pi; % ws = 0.3*pi; % Rp = 1; % As = 15; % % ?￡?a?-Dí??±ê???μ?êμ???ó3é? T = 0.001; Fs=1/T; % OmegaP = wp*Fs;