数字信号处理（II）.docVIP

数字信号处理(II) 实验一 维纳滤波器的计算机实现 班级：硕0711班 姓名： 学号： 三．实验内容与步骤 1. 仔细阅读维纳滤波原理，根据图1.1(附录)给出的框图编制维纳滤波程序。 维纳滤波源程序见附录1 2． 运行维纳滤波程序，选择L=500,N=10，观察并记录实验结果，分析比较下列三 个问题： ① 与s(n)比较，信号x(n)在维纳滤波前后有何差别？滤波效果如何？ ② 估计出的h~ (n)和理想的h(n)比较，近似程度如何？ ③ 理想的维纳滤波和FIR 维纳滤波效果有何差异？ 选择L=500，N=10，实验结果如下： 记录三个均方误差的打印结果如下： Ex =0.9946 Ei =0.2518 Er =0.3160 其中Ei表示s(n)和理想估计的s(n)的均方误差，Er表示s(n)和实验估计的s(n)的均方误差，Ex表示x(n)和s(n)的均方误差。 (1).由图1、图3和图4可以看出，与s(n)比较，信号x(n)再滤波后信号明显好的多。滤波后的x(n)在大体上与s(n)十分近似，只是细节上稍有不同。理想滤波的均方误差是0.2597，实际滤波后的均方误差是0.3347，滤波效果一般。 (2).由图2可以看出，估计出的h~ (n)和理想的h(n)基本相同，只是在局部上下有些小小的浮动。 (3).图3和图4告诉我们，理想的维纳滤波和FIR 维纳滤波效果基本是很相似的，只有些微小的不同，但由于理想的维纳滤波是用阶数为无穷大的滤波器滤波（很难实现），所以效果还是要好一些，从两个的均方误差也可以看出来。 3．固定L=500, 分别取N=3、20, 根据实验结果，观察N 的大小对h~(n)的估计和 滤波效果的影响。记录实验结果。 a.L=500,N=3, 实验结果如下： 打印的三个均方误差分别为 Ex =1.0189 Ei =0.2168 Er =0.3866 b. L=500,N=20，实验结果如下： 打印的三个均方误差分别为 Ex = 1.0176 Ei =0.2026 Er =0.2427 由上面两次的实验可以看出来，固定L，随着N的增大，h~(n)越接近理想滤波的h(n)，精度也越高。由图和打印出的均方误差可以知道，N越大，滤波的效果也越好，和s(n)也越相似。 4． 固定N=10, 改变L=200,1000, 根据实验结果，观察并记录L 的大小对h~(n)的精度和滤波效果的影响。 a.L=200,N=10，实验结果如下： 打印的三个均方误差分别为 Ex =1.0854 Ei =0.3105 Er =0.2576 b.L=1000,N=10，实验结果如下： 打印的三个均方误差分别为 Ex =1.0044 Ei =0.2657 Er =0.3388 由上面两次的实验可以看出来，固定N，随着L的增大，h~(n)越接近理想滤波的h(n)，精度也越高。滤波效果也会更好。 5． 仔细阅读有关维纳预测原理，弄清AR 模型参数估计与维纳预测的关系，根据框图1.2(附录) 编制AR 模型参数估计程序。 AR 模型参数估计源程序见附录2 6． 运行AR 模型参数估计程序，选择p=1, ＝－0.6，L=100, ，观察并纪录与的最佳估值，与理论值进行比较。 实验结果如下： 选择p=1, =-0.6,L=100, =1，运行AR参数程序，得到估计的= -0.5536,噪声方差=1.3503。实验所得到的估计值接近于理论值。 7． 固定p=1, ＝－0.6，L=100, =1改变L=50,500, 观L 的大小对AR 模型参数估计的精度的影响。 选择p=1, =-0.6,L=50,噪声方差为1，运行AR参数程序，得到估计的=-0.6302,噪声方差=1.1760。 选择p=1, =-0.6,L=500,噪声方差为1，运行AR参数程序，得到估计的=-0.6333,噪声方差= 1.0805。 由实验可以看出，固定p, ,L, ，随着L的次数的增加，AR模型参数估计的精度也就越高。 四．思考题 1.推导出公式(1-13)，验证(1-16)的结果。 (1)推导公式(1-13) 由题知 s(n)与v(n)互不相关， 又 作Z变换 又 其中由单位圆内的零极点组成，由单位圆外的零极点组成。 物理实现 对于, 只取 (2)验证(1-16)的结果 又 取单位圆为积分围线，上式等于单位圆内的极点(Z=0.95及Z=0.724)的留数之和。 2.由公式1-21，怎样得到, ，和？ Levinson-D