通信系统计算机模拟第七讲.ppt

通信系统的计算机模拟 第七讲 Announcement Submit your proposal today 第六讲 回顾 IIR and FIR Filter IRR滤波器的实现 产生转置直接Ⅱ型滤波器的规则 矩阵形式 Maltab生成 FIR滤波器的实现 IIR滤波器：综合方法与滤波器特性 冲激不变滤波器 阶跃不变滤波器 双线性z变换 数字滤波器的计算机辅助设计 5.5 FIR滤波器：综合技术与滤波器特性 如果一个滤波器的冲激响应是有限长的，or 冲激响应截短为有限长的，则离散时域的滤波器输出信号为： FIR滤波器的吸引力 1．通信中有的滤波器不能用拉普拉斯变换表示，即不能直接用IIR，如平方根升余弦（SQRC）脉冲整形滤波器和Jakes多普勒滤波器（Jakes Doppler filter）。FIR设计方法则可以轻易地实现这些滤波器的仿真。 2．在许多仿真应用中，滤波器数据可能是以测量的频率响应或脉冲响应的形式经验性地给出， FIR仿真容易。（虽然有些方法可以拟合成AMMA自回归滑动平均，但频率响应数据不平滑时误差大） 3．采用FIR滤波器，可以给定任意幅度响应和相位响应。 4．FIR滤波器没有反馈，因此总是稳定的。 FIR仿真模型的不足 计算效率不如IIR滤波器实现的高。 （5-60）定义的卷积运算，产生一个输出采样点就需要N次（复数）乘法和加法运算，比如说N1024 点的话，FIR算法运行起来比典型的IIR滤波器算法 慢得多。 IIR的复杂度由滤波器的阶数决定的，而不同于FIR滤波器，由冲激响应的长度决定。 克服FIR仿真模型不足的措施 时域卷积与基于DFT运算的卷积的计算效率 时域卷积与基于DFT运算的卷积的计算效率之比值跟log2N/N成比例。 如果冲激响应长度少于128个采样点的话，时域卷积与基于DFT/FFT的卷积在计算量上没有太大的区别。 如果对冲激响应的长度进行截短，FIR滤波器方法计算效率可以得到改进。截短等价于使用一个长度为N个采样的窗函数w[k]与冲激响应h[k](k=1,2,….)相乘。 N值的选择？该窗函数内至少包含有总能量的98%。 窗函数 时域的截短（加窗）等价于频域的卷积 矩形窗 窗函数 窗函数 5.5.1 根据幅度响应进行的FIR设计 依据：滤波器的频率响应（包括幅度和相位）和单位冲激响应是一个傅里叶变换对。 推导单位冲激响应的一般做法：给定期望的幅度响应A(f)，再计算其傅里叶反变换 通常目标幅度响应为实偶函数，单位冲激响应也是实偶函数型的。因为冲激响应是偶函数，是非因果的. 如何变成时域上可实现的系统？ 必须对所得的冲激响应进行截短使其为有限长 问题及解决办法？截短冲激响应时，必须小心以避免引入显著的误差，适当的窗函数可以减少误差的影响 还要在时间轴上进行平移，使之成为因果函数 在时间轴上平移冲激响应，相当于使滤波器获得一个大小等于该时间平移长度的群时延的线性相位响应，可以得到具有任意幅度响应和线性相位响应的滤波器。 想要一个同时具有目标幅度响应和相位响应的滤波器，就要给定一个复传递函数。假定截断误差不显著，那么除了定常群时延外，最终得到的滤波器可以同时满足目标滤波器对幅度响应和相位响应的要求。 根据幅度响应进行的FIR设计方法 数字滤波器的幅度响应，是频率f的连续函数，是以采样频率为周期的周期函数，可表示为傅里叶级数。 这些离散的傅里叶级数的系数，给出期望的数字滤波器的冲激响应。 期望频率响应的IFT or期望频率响应采样点的IFFT?冲激响应 这些方法很基本，小心使用会得到有用的滤波器设计 根据幅度响应进行的FIR设计方法 由于数字滤波器的幅度响应具有以采样频率为周期的周期性，其幅度响应可以展开为傅里叶级数，形式如下： 进一步化简 设M=2L+1，利用这一替换并改变序数k=n-L，可得： 分析 h1[k]=h[k+L]这是由式（5-64）所定义的因果滤波器的冲激响应平移L个采样点后得到的。 （5-68）定义滤波器不是一个因果滤波器。 但是，它便于采用基本的傅里叶变换法进行设计，并且只要简单地对冲激响应进行合适的平移就可以变为因果的。 很显然，由h[n]和h1[n]所确定的滤波器的幅度响应是相同的，而传递函数只是相差一个由式（5-67）所定义的线性相位平移。 ?因此，设计过程基于式（5-68）。 设计过程 假设一个滤波器将具有给定幅度响应 求得基本设计方程 这是我们的基本设计方程 例5-8 理想数字低通滤波器 我们设计了一个数字低通滤波器，来近似一个带宽为 的理想数字滤波器，这里fN为奈奎斯特频率，fs为采样频率，?是介于0和1之间的一个参数。因此，期望的幅度响应为 解： 从（5－73）