哈尔滨工业大学威海校区《数字信号处理》实验二.docVIP

数 字 信 号 处 理 实验报告 实验名称： 实验二 IIR滤波器设计 实验日期： 2011．11．19 姓 名： 尤伟 学 号： 090240328 哈尔滨工业大学（威海） 实验二 IIR滤波器设计 实验目的 掌握冲激响应法和双线性变换法设计IIR滤波器的原理及具体设计方法，熟悉用双线性设计法设计低通、带通和高通IIR数字滤波器的计算机程序； 熟悉模拟Butterworth滤波器和Chebyshev滤波器的设计，掌握冲激响应法和双线性变换法设计数字IIR滤波器的方法。 实验原理和方法 连续系统系统函数和离散时间系统函数之关系用模拟信号的采样信号作为两者之间联系的桥梁。设某连续时间系统，其频率响应，冲击响应；对采样得到即 上式进行拉氏变换，得到 令离散序列h[n]=h[nT]，则 实验内容 不同阶次模拟巴特沃兹滤波器的频率响应 程序： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 数字信号处理实验 - IIR滤波器设计 % 1 利用Matlab系统函数计算不同阶次的巴特沃兹滤波器频率响应 % 2 根据提出模拟IIR的指标，计算IIR的阶次、系数等 % 姓名：尤伟 % 学号：090240328 % 时间：2011.11.19 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 1 利用Matlab系统函数计算不同阶次的巴特沃兹滤波器频率响应 clc,clear all,close all; freq_axis = [ 0: 0.001 :2 ]; %归一化的频率 N_array = [ 3 5 10 12 ]; %滤波器的阶次 figure(1); for I=1:4 N = N_array(I); [z p k] = buttap(N); %获得5阶巴特沃兹滤波器的参数 [b a] = zp2tf( z, p, k ); %得到传递函数 tf:trans-function disp( [ N= num2str(N) 极点: ] ); disp(p); disp( [ N= num2str(N),零点: ] ); disp(z); [H ] = freqs(b, a, freq_axis); subplot( 2, 2, I ); plot( freq_axis, abs(H)); axis( [0 2 0 1] ); %限制坐标轴 xlabel( w/wc -- 归一化频率 ); ylabel( 幅度谱 ); title( [巴特沃兹滤波器幅度谱 N=,num2str(N)] ); end %2 根据提出模拟IIR的指标，计算IIR的阶次、系数等 J = sqrt(-1); fp = 5e3; %5KHz 通带截止频率***模拟频率 alfa_p = 2; %最大通带衰减2dB fs = 12e3; %12Kz 阻带截止频率***模拟频率 alfa_s = 30;%最小阻带衰减30dB k_sp=sqrt((10^(0.1*alfa_p)-1)/(10^(0.1*alfa_s)-1));%求ksp lamda_sp = 2*pi*fs/(2*pi*fp); %求lanmda sp N = ceil( -log(k_sp)/log(lamda_sp) ); %求N，ceil为取就进整数 p = zeros(1,N); for I=1:N p(I)=exp(J*pi/2+J*pi*(2*I-1)/2/N ) ; %求极点Sk end z = []; %表示无零点 [b1 a1] = zp2tf(z,p,1); %得到传递函数 tf:trans-function fc=fp*(10^(0.1*alfa_p)-1)^(-1/2/N) ; %3dB截至频率 figure(2); freq_axis1 =[0:50:15e3]; [H1 ] = freqs(b1, a1, freq_axis1/fc);%对Wc作频率归一化 %Frequency response of analog filters plot( freq_axis1 , 20*log10(abs(H1))