DSP课程设计---基于DSP的实时FIR滤波器的实现.docVIP

设计题目 IR滤波 设计要求 设计过程 1.设计原理 通过反z变换，数字滤波器的差分方程为： 由上式可以得出如下图所示的直接型结构，这种结构又可以成为卷积型结构。将转置理论应用于图1.1可以得到转置直接型结构。 将式中的系统函数H（z）分解成若干一阶和二阶多项式的连乘积： 则由此式可以构成如图1.2所示的级联型结构。其中为一阶节，为二阶节。每个一阶节、二阶节可用图1.2级联型结构实现。当M1=M2时，即得到图1.3所示的具体结构。这种结构的每一节都便于控制零点，在需要控制传输零点时可以采用。但是它所需要的系数a比直接型的h（n）多。在对滤波器计算时间没有特殊要求的时候可以采用这种形式。若需要严格考虑滤波器的计算时间则需要同时考虑它们的优点及缺点来设计。这在算法设计的时候要使用软件编辑环境来计算运行的时间问题。通常FIR的计算时间都较长。很多时候我们需要牺牲时间来获得想要得到的滤波区功能。 图1.1 FIR滤波器直接型结构图 图1.2 级联型结构图 图1.3 级联型具体结构 FIR滤波器实质上就是一个分节的延迟线，把每一节的输出用滤波器系数进行加权累加，便得到滤波器的输出结果，它总是稳定并且可实现的。在一些工程实际应用（如：图像处理、数据调制解调）中，往往对相位要求较高。FIR滤波器可以实现严格的线性相位，从而得到了广泛应用。它的差分方程数学表达式为： 式中，N是FIR 滤波器的抽头系数，x（n）表示在n时刻输入的信号样值，h（n）表示滤波器的第n级抽头系数。横截型FIR滤波器的结构如图1.4所示： 图1.4 FIR滤波器的横截型结构 2设计方案是无法实现的，FIR的设计就是要寻找一个可以得到的频率特性来逼近，这相当于用一个可实现的单位脉冲响应 h（n）去逼近一个理想单位脉冲响应。可由理想频率特性通过傅氏反变换得到，。一般来说，这样得到的理想单位脉冲响应序列是个无限长序列，因而是非因果的。设有一个截止频率为的理想线性相位低通，延时为，其频率特性是： 得到： 这是一个以n=为中心偶对称的无限长非因果序列，要想用一个有限长的因果序列去逼近它，最简单的方法是截取从0到N-1的一段来表示它，及=;其他N：。 同时，为了保证线性相位，还要满足偶对称。这就像透过一个窗口观看到的一段，因此就表示成和一个窗口函数的乘积，这样对的求解就变为，这里的就被称为窗函数，既然一个频域上的标准的矩形窗口对应于时域是一个无限长的序列，那么在时域上截取一段必然会造成频域的矩形窗口的失真。相应的，截取出的信号也相应失真，从而需要改变原来窗口的形状来修正经过时域截取后的窗口失真。 2.2 方案二: 频率抽样法 窗函数法是从时域出发，用窗函数截取理想的得到，以此有限长的近似，这样得到的频率响应逼近于理想的频率响应。频率抽样法是从频率出发，将给定的理想频率响应加以等间隔抽样 然后以此作为FIR滤波器的频率响应抽样值，再根据DFT（离散傅氏变换）定义由频域这N个抽样值来唯一定义一个有限长序列，同样也可以算出FIR滤波器的系统抽样函数及频率响应，可以推出频率响应是频率抽样值与线性相位因子及如下内插函数的线性组合。 所以，在各频率取样点上，实际滤波器的频率响应是严格地和所要求的滤波器的频率响应一致的，逼近误差为零，但在抽样点之间的频率响应是各取样点的内插函数的延伸叠加而成，有一定的逼近误差，误差大小取决于频率响应曲线的圆滑程度和抽样点的密度。为了减少误差就要增加抽样点数目即增大采样频率，抽样点之间的理想频率特性变化越陡，则逼近误差越大，在理想频率特性的不连续点附近会产生肩峰和纹波。频率抽样法的优点是可以在频域设计，适于利用最优化方法，而且这种方法特别适用于窄带选频滤波器，但频率抽样法的抽样频率只能是2π/N的整数倍或2π/N的整数倍加上π/N。不能保证截止频率的准确取值，要实现精确的就必须取N足够大，相应的计算量也很大。此外，它的阻带最大衰减一般，也只有30-50dB左右，很难满足频域特性要求较高的场合。 2.3 方案三: Chebyshev逼近法 窗函数法和频率采样法设计的滤波器的频率特性都是在不同意义上对所给理想频率特性的逼近。由数值逼近理论可知，对某个函数f（x）的逼近一般有以下三种方法： ①插值法（Interpolating Way） ②最小平方逼近法（Least Square Approaching Way） ③一致逼近法（Consistent Approaching Way） 切比雪夫最佳一致逼近的基本思想是，对于给定区间[a,b]上的连续函数f(x)，在所有n次多项式的集合中，寻找一个多项式p(x)，使它在