DSP课程设计概要.ppt

IIR：冲激响应不变法 双线性变换法 FIR：Hamming窗函数设计 * * * * 王庆玲 王鹏 秦绯雨 冲激响应不变法设计IIR数字滤波器 用数字滤波器的单位冲激响应 模仿模拟滤波器的单位冲激响应 基本原理 ZT[h(n)]=H(z) LT[ha(t)]=Ha(s) S平面到z平面的映射关系 冲激响应不变法特点 h(n)完全模仿模拟滤波器的单位抽样响应 相位保持线性关系 频率响应混迭 参数设定 ： 阻带截止频率850Hz ：阻带最大衰减40dB ：通带截止频率650Hz ：通带最大衰减1dB [N, OmegaC] = buttord(OmegaP, OmegaS, Rp, As, s); %通过buttord函数，由通带阻带截止频率计算阶数和3dB截止频率 [b,a] = butter(N, OmegaC, ‘s’); %由计算的阶数，通过butter计算频率响应的分子分母系数 ，完成模拟滤波器设计 [bz,az] = impinvar(b,a,Fs);%通过impinvar函数将模拟滤波器转换为数字滤波器 dbHx = 0.8017 41.9396 双线性变换实现IIR数字滤波器 一、变换原理： 使数字滤波器的频率响应 与模拟滤波器的频率响应相似。 避免多值映射 优缺点 优点：避免了频率响应的混迭现象 s 平面与 z 平面为单值变换 缺点：除了零频率附近，? 与w之间有严重的非线性。 模拟滤波器的数字化方法 参数选择 Rp=3; As=40; wp=0.132*2*pi; ws=0.15*2*pi; 加的是高斯白噪声 实现步骤 1)频率预畸 将数字域各临界频率预畸为各个模拟域临界频率 Wp=2*fs*tan(wp/2);Ws=2*fs*tan(ws/2); 2）利用buttord函数，由通带阻带截止频率计算阶数和3dB截止频率 [N,wc]=buttord(Wp,Ws,3,40,s); 3）由计算的阶数，通过butter计算频率响应的分子分母系数 ，完成模拟滤波器设计 [b,a]=butter(N,wc,s); 4）bilinear把模拟滤波器转化为数字滤波器 [bz,az] = bilinear(b,a,fs); Fir df-窗函数法 窗函数设计法 设计思想： 时域逼近理想数字滤波器 一、设计方法： 要选择合适的形状和长度 w(n)：窗函数序列 线性相位FIR低通滤波器设计 1）求数字频率 2）确定理想低通频响并求hd(n) 3）选择窗函数：由As=50db确定hamming窗 4）确定N 值 5）确定FIR滤波器的h(n) 6）求 H及As ，验证 根据自己的音频设定模拟域的截至频率，以及抽样频率 OmegaP=2*pi*5000; 通带截止频率 OmegaS=2*pi*7000; 阻带截止频率 Fs=44100; 抽样频率 设定阻带衰减，选择窗函数 As=50; 阻带衰减 select Hamming Window 窗的点数 验证阻带衰减是否满足要求 满足 加噪 加入两个高频余弦信号 和信噪比为25db的高斯白噪声 * *