04实验四白化滤波器的设计实验报告

《04实验四白化滤波器的设计实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验五白化滤波器的设计【实验目的】了解白化滤波器的用途，掌握白化滤波器的设计方法。【实验器材】1.设备：通用计算机2.软件：MATLAB6.5版木以上【实验原理】虽然理论中我们常常以白噪声最为研宄对象，但是在统计信号处理中，往往会遇到等待处理的随机信号是非白色的，这样会给问题的解决带来困难。克服这一网难的措施之一是，对色噪声进行白化处理。主要内容是设计一个稳定的线性滤波器，将输入的色噪声变成输出的白噪声。在这里，我们就对一般的具有功率谱Gv（的的平稳随机过程X（t）白化处理问题进行讨论。为了具体的进行分析和计算，假设G、.伸）可以表达成有理数的

2、形式，即其中分子、分母为多项式。这个假设对于通常见到的功率谱是很近似的，而II有可行的方法用有理数去逼近任意的功率谱密度。由Tga（的是功率谱，它的平稳随机过程相关函数的傅里叶变换具有非负的实函数和偶函数的性质。这些性质必然在其有理函数的表示式中体现出来，特别是，Gv（的的零、极点的分布和数量会具冇若干个特点。由于Gx（的是实函数，因此有：Ga（CO）=G；（CO）,tz2是实数，G，加）的零、极点是共扼成对的。从而也可以把GA加）的表示式写成如下形式：(>694-6Z,)...…（加+乂）aO+A)……(加+A)(-jco+a^(-jeo^a

3、k)把仍开拓到复平面s中去，另=+用S代替就可以把函数GA.（⑻扩大到整个复平面。G/的的零、极点必将对称于轴，如下图所示：由于Gv（fd）是偶函数，因此不难判断，GA（fzJ）的零、极点是象限对成的，从而对于轴也是对称的。由于GA加）20,因此分子的虚根必然是偶数倍数个，否则GA（叫会出现负值。这就是说y⑺轴上的零、极点必将成对的出现。由于gv（的是可积的，因此分子的阶数不能大于分母的阶数，这就是说零点令：G:⑽=G;⑽=则有总数不会大于极点总数，而且分母不可能有虚根，这意味着轴上没有极点。综合上述情况，在S平面的零、极点的可能为置如上图所示

4、：(./必+汉1)……(y•仍+乂)^(>69+^,)……(j>y+A)(->69+or,)……(-jco^ak)a(-jo)D……(-加+A)Gv⑽=G:⑽［G:⑽］'=G:⑽G;⑽其中G了（⑼代表零、极点均在s左平面的部分，代表零、极点均在s右平面的部分。若在jTy轴上有零点的话，必是成对的。则将一个放在6^（⑼内，将另一个放在⑽内。实质上，G;⑽对应的时域函数在负时间域为零，而G:（⑼对应的时域函数在正时间域为零。根据上述的讨论，可以求得闩化滤波器的的解析式为:

5、/w)l2=由于

6、H,⑻)

7、2=什⑽H;⑽=H'⑽仏(-co)G»=G:(a

8、))G；(co)=G：⑽G:(-的故得：G:⑽若运用傅里叶变换进行分析计算，以s代替jTy，可得白化滤波器公式H,（5）=—-—其中+G:⑴我们知道，仏咏）的傅里叶反变换是白化律波器在时域的单位冲击响应⑽，⑽零、极点在S左半平面’因此志的零、极点也是在S左半平面。故它对应的时域函数^（0在负时域时为零，也就是说，上述白化滤波器是物理可实现的。白化滤波器的设计方法是：首先计算色噪声自相关函数，根据色噪声的自相关函数，计算出色噪声的功率谱（色噪声的自和关函数和功率谱构成一对傅里叶变换对），然后将色噪声的功率谱分解成在S的左半平面和右半平面。取功率谱

9、在S的左半平面的那些值，然后找出它们的零、极点。则根据公式G:⑽二……(jco+ak、a……(7^+A)和⑽二忐计算出白化滤波器的频谱。6Z为零点，A为极点。【实验任务3（1）编写白化滤波程序用Matlab语言编写白化滤波器程序，设计的白化滤波器能够将有色噪声转变为白噪声。白化滤波器的设计步骤：①产生色噪声②计算色噪声的功率谱（先求自相关，在求功率谱）③计算白化滤波器的传输函数（取色噪声功率谱在S平面的左半部分，求倒数）④色噪声通过白滤波器（使用快速卷积的实现方法）⑤计算白化滤波器的统计特性。（2）用软件方法产生一个有色噪声信号，色噪声的数学表

10、达式为：2*exp（l）*rand（）,使该信号通过白化滤波器。(3)引用实验一中编制的测试均值、方差、和关函数(包括自相关函数)及实验三中功率谱密度的程序计算有色噪声的这些参数，并将计算结果存入文件中O【实验报告要求】(1)简述实验目的及实验原理。(2)按实验的记录数据画岀有色噪声信号和其通过白化滤波器后的均值方差、相关函数及功率谱密度图形。(3)按实验步骤纪荥观察到的均值方差、相关函数、频谱及功率谱密度，并对所得结果进行分析和解释。%**林林林*色噪声的读取林本林fsM4000;[xl,Fs]=audioread(，C:UsersAdm

11、inistratorMusicfe噪声.wav’)；sound(xl,fs);x=xl(1:44000);t=(0:length(xl)-l)/f