资源描述:
《实验一白噪声测试》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实用标准文档白噪声测试一、实验目的⑴了解白噪声信号的特性,包括均值(数学期望)、均方值、方差、相关函数、概率密度、频谱及功率谱密度等。⑵掌握白噪声信号的分析方法。二、实验原理所谓白噪声是指它的概率统计特性服从某种分布而它的功率谱密度又是均匀的。确切的说,白噪声只是一种理想化的模型,因为实际的噪声功率谱密度不可能具有无限宽的带宽,否则它的平均功率将是无限大,是物理上不可实现的。然而白噪声在数学处理上比较方便,所以它在通信及电子工程系统的分析中有十分重要的作用。一般地说,只要噪声的功率谱密度的宽度远大于它所作用的系统的带宽,并且
2、在系统的带内,它的功率谱密度基本上是常数,就可以作为白噪声处理了。白噪声的功率谱密度为:其中为单边功率谱密度。白噪声的自相关函数为:白噪声的自相关函数是位于τ=0处、强度为的冲击函数。这表明白噪声在任何两个不同的瞬间的取值是不相关的。同时也意味着白噪声能随时间无限快的变化,因为它的带宽是无限宽的。下面我们给出几种分布的白噪声。随机过程的几种分布前人已证明,要产生一个服从某种分布的随机数,可以先求出其分布函数的反函数的解析式,再将一个在[0,1]区间内的均匀分布的随机数的值代入其中,就可以计算出服从某种分布的随机数。下面我们就
3、求解这些随机数。[0,1]区间均匀分布随机信号的产生:采用混合同余法产生[0,1]区间的均匀分布随机数。混合同余法产生随机数的递推公式为:n=0,1,2……文案大全实用标准文档n=1,2,3……由上式的出如下实用算法:其中:,其中k为计算几种数字尾部的字长,t为任意选定的正整数,为任意非负整数,为奇数Matlab语言中的rand()函数是服从[0,1]均匀分布的,所以在以后的实验中如果用到均匀分布的随机数,我们统一使用rand()函数。正态分布(高斯分布)随机信号的产生:高斯分布的密度函数为:采用变换法产生正态分布随机数,若
4、、示[0,1]均匀分布随机数,则有正态分布随机数:指数分布随机信号的产生:指数分布的密度函数为:当x>0时,当x≤0时f(x)=0,其中α>0它的反函数(指数分布随机数)为:其中r为[0,1]区间均匀分布的随机数。三、实验内容与结果1.产生五种概率分布的信号Matlab程序:文案大全实用标准文档%生成各种分布的随机数x1=unifrnd(-1,1,1,1024);%生成长度为1024的均匀分布x2=normrnd(0,1,1,1024);%生成长度为1024的正态分布x3=exprnd(1,1,1024);%生成长度为102
5、4的指数分布均值为零x4=raylrnd(1,1,1024);%生成长度为1024的瑞利分布x5=chi2rnd(1,1,1024);%生成长度为1024的卡方分布%时域特性曲线:figure;subplot(3,2,1),plot(1:1024,x1);gridon;title('均匀分布');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);subplot(3,2,2),plot(1:1024,x2);gridon;title('正态分布');xlabel('时间(t)');yl
6、abel('幅度');axis([01024-22]);subplot(3,2,3),plot(1:1024,x3);gridon;title('指数分布');xlabel('时间(t)');ylabel('幅度');axis([01024-15]);subplot(3,2,4),plot(1:1024,x4);gridon;title('瑞利分布');xlabel('时间(t)');ylabel('幅度');axis([01024-14]);subplot(3,2,5),plot(1:1024,x5);gridon;tit
7、le('卡方分布');xlabel('时间(t)');ylabel('幅度');axis([01024-15]);文案大全实用标准文档2.均值:均值E[x(t)]表示集合平均值或数学期望值。基于随机过程的各态历经性,可用时间间隔t内的幅值平均值表示:均值表达了信号变化的中心趋势,或称之为直流分量。在MATLAB中,可以用mean()函数来计算。文案大全实用标准文档%求各种分布的均值figure;m1=mean(x1);m2=mean(x2);m3=mean(x3);m4=mean(x4);m5=mean(x5);subplo
8、t(3,2,1),plot(1:1024,m1);title('均匀分布均值');xlabel('时间(t)');ylabel('幅度');axis([01024-22]);subplot(3,2,2),plot(1:1024,m2);title('高斯分布均值');xlabel('时间(