白噪声的产生与测试实验汇总

1、白噪声的产生与测试第二组m噪声的产生与测试、实验目的了解白噪声信号自身的特性，包括均值、均方值、方差、相关函数、概率密 度、频谱及功率谱密度等。掌握白噪声的分析方法。熟悉常用的信号处理仿真软 件平台matlab软件仿真。了解估计功率谱密度的几种方法，掌握功率谱密度估 计在随机信号处理中的作用、实验原理所谓白噪声是指它的概率统计特性服从某种分布而它的功率谱密度又是均 匀的。确切的说，白噪声只是一种理想化的模型，因为实际的噪声功率谱密度不 可能具有无限宽的带宽，否则它的平均功率将是无限大，是物理上不可实现的。然而白噪声在数学处理上比较方便，所以它在通信系统的分析中有十分重要的作 用。一般地说，只要

2、噪声的功率谱密度的宽度远大于它所作用的系统的带宽，并 且在系统的带内，它的功率谱密度基本上是常数，就可以作为白噪声处理了。白 噪声的功率谱密度为:Sn（f）N其中N0/2就是白噪声的均方值。白噪声的自相关函数位:R（）牛（）白噪声的自相关函数是位于T =0处、强度为严的冲击函数。这表明白噪声在任何两个不同的瞬间的取值是不相关的。同时也意味着白噪声能随时间无限快的变化，因为它的带宽是无限宽的。下面我们给出几种分布的白噪 声。随机过程的几种分布均匀分布随机信号、正态分布（高斯分布）随机信号、指数分布随机信号等。实验任务与要求以及从随机噪 通过实验要求掌握几种分布的随机噪声共同点和不同点, 声的相关

3、和功率普中得到白噪声的特征，重点在于系统测试与分析。实验系统框图如图2-1、图2-2所示:图2-1 各种分布随机信号测试5介RS机a程V信号特性测试及绘制曲线图2-2随机信号叠加后的特性测试 自选matlab或C/C+软件之一产生几种概率分布的仿真随机信号：随机数的 长度N=1024这些随机数包括均匀分布、正态分布、指数分布、瑞利分布、 分布。并计算这些随机数的均值、均方值、方差，自相关函数、概率密度、频谱 及功率谱密度并绘图。分析实验结果，搞清楚均值、均方值、方差，自相关函数、 频谱及功率谱密度的物理意义。验证当N增大时，白噪声的功率普密度逼近号。设产生N=2048a长度的（2, 3）正态随

4、机随机数，从中取1024、10240、20480个点的功率普密度，做比较, 观察这些随机数的功率谱密度随长度的变化。实际的白噪声功率普密度不是常 数。 根据白噪声的特性，确定哪些随机信号属于白噪声范畴。根据分析确定白噪声与概率分布有关系吗?通过编程分别确定当5个均匀分布过程、5个指数分布分别叠加时，结果 是否是高斯分布。叠加次数对结果的影响?四、实验步骤与结果1.产生五种概率分布的信号，如下图:瑞利分布时域特性曲线3.53D18卡方分布时域特性曲线16141264200.806040.2-0.4-0.G心82004008001000600 时间t庖0 ly-0.2> 10蝕 loB0 8

5、12001200W2.均匀分布、正态分布、指数分布、瑞利分布、2方分布均值、均方值、方差，自相关函数、概率密度、频谱及功率谱密度等参数图像：，如下图:均值：均值 Ex(t)表示集合平均值或数学期望值。基于随机过程的各态历经性, 可用时间间隔t内的幅值平均值表示：JT-lEk©=篙0阿c=4均值表达了信号变化的中心趋势，或称之为直流分量。在MATLAB，可以用mean()函数来计算。高斯分布均值0502001000400600800时间(t)-0,512002指数分布均值1 50,502001000400600800时间(t)1200瑞利分布均值卡方分布均值方差：随机过程的方差函数描述

6、了随机过程所有样本函数在t时刻的函数值相对于其数学期望的偏离程度。定义：M-1臣=CO - ERG)尸fNc=a其中(T (t)是随机过程的标准差。当随即过程表征的是接收机输出端的噪声电压时，T 2(t)表示小号在单位电阻上的瞬时交流功率统计平均值，而T (t)表示噪声电压相对于电压统计平均值的交流分量。在MATLAB，可以用std()函数计算出标准差t (t)，再平方就可以得到方 差。re均匀分布方差0.50-0.50200WOO400600800时间1200高斯分布方差2O5Q瑞利分布方差2001000400600800时间5 Q541200卡方分布方差(2)_unuL nu5i=i对于平

7、稳随机过程x(t) 可以用相关函数表示。自相关：信号的相关性是指客观事物变化量之间的相依关系。 和y(t)在两个不同时刻t和t+ T的起伏值的关联程度， 在离散情况下，信号x(n)和y(n)的相关函数定义为：jr-iM-iT=0 t随机信号的自相关函数表示波形自身不同时刻的相似程度。与波形分析、频 谱分析相比，它具有能够在强噪声干扰情况下准确地识别信号周期的特点。卡方分布自相关函数图指数分布自相关函数图瑞利分布e相关®数團概率密度函数：一维分布函数为:若Fx（x i；t 1）对xi的一阶偏导存在，则一维概率密度为:0F愛（尤1；如）&（X；r訂=典'丿在MATLAB，

8、可以用ksdensity（）函数来计算一维概率密度。频谱：信号频谱分析是采用傅立叶变换将时域信号x(t)变换为频域信号x(f)，从另一个角度来了解信号的特征。时域信号x(t )的傅氏变换为：在MATLAB中，对信号进行快速傅立叶变换fft()就可以得到频谱函数。60卡方分布幅频特性频率也均匀分布幅频特性4035302520151050 W0 200300400500600700 RD0 900 WOO频率也ly频率H2频率H2功率谱密度：随机信号的功率谱密度是随机信号的各个样本在单位频带内的频谱分量消耗在一欧姆电阻上的平均功率的统计均值，是从频域描述随机信号的平均统计参量，表示x(t)的平均功

9、率在频域上的分布。它只反映随机信号的振幅信息，而没 有反映相位信息。在MATLAB中，可由下式得到功率谱密度：S3 = lim半服ST)鬥2.5°0正态分布功率谱密度2150.520010001200400600800频率Hz频率Hz卡方分布功率谱密度454as2001000400600800频率Hz120035«利分布功率谱密度2.51502001000400600800频率Hz12003.5指数分布功率谱密度32.5202D0400>WIOQ630频率WBOO WOD 1200取20480个点时的功率谱密度和自相关函数，如下图(1)功率谱密度:频率Hz频率Hz14

10、5|W>g加4和点正态分布功率谙S5 X山却随机信号叠加:频率出250均匀分布盏加2001501006033.544.551.522.50.52502001501005040455015202530350510指数分布S加4.随机信号检验:指数分布正态检验0 9990 9970 950,96 950.90n- n- 口 uacziADOo o1 Zr I51015202530354UDala459 9 Ao no o0.51.522.533.544.5Dala五、实验总结这次试验让我们对白噪声有了很大的理解，最主要是在实验过程中用到了好久不用的matlab软件，

11、由于好长时间不用好多的函数的功能都忘记了，而且实验 过程中用到的好多函数以前都没接触过，所以还得花好长时间去查阅相关资料。这次试验的目的其实让我们学会是利用 matlab软件对信号分析，同时加深我们 对信号和噪声参数处理的理解，锻炼我们的实践动手能力。参考文献:随机信号分析教程 高等教育出版社 主编：李兵兵matlab7.x程序设计语言 西安电子科技大学出版社 主编：楼顺天现代通信原理与技术西安电子科技大学出版社主编：张辉高频电子线路高等教育出版社主编：曾兴文%生成各种分布的随机数x1=unifrnd(-1,1,1,1024)%; 生成长度为 1024的均匀分布 x2=normrnd(0,1,

12、1,1024)%; 生成长度为 1024的正态分布 x3=exprnd(1,1,1024)%; 生成长度为 1024的指数分布均值为零 x4=raylrnd(1,1,1024)%; 生成长度为 1024的瑞利分布 x5=chi2rnd(1,1,1024)%; 生成长度为 1024的卡方分布 %时域特性曲线： figure;subplot(3,2,1),plot(1:1024,x1);gridn;titleC均匀分布 时域特性曲线);xlabel(时间(t)');ylabel幅度(v)');axis(0 1000 -2 2 );subplot(3,2,2),plot(1:1024

13、,x2);gr0ln;title('正态分布');xlabel(时间(t) ');ylabel(幅度(v) ');axis(0 1000 -2 5 );subplot(3,2,3),plot(1:1024,x3);gridn;title('指数分布');xlabel(时间(t) ');ylabel(幅度(v) ');axis(0 1000 -2 5 );subplot(3,2,4),plot(1:1024,x4);gridn;title('瑞利分布');xlabel(时间(t) ');ylabel(幅度(v

14、) ');axis(0 1000 -2 5 );subplot(3,2,5),plot(1:1024,x5);gridn;title('卡方分布');xlabel(时间(t) ');ylabel(幅度(v) ');axis(0 1000 -2 5 );%求各种分布的均值 figure;m1=mean(x1),m2=mean(x2),m3=mean(x3),m4=mean(x4),m5=mean(x5);sub plot(3,2,1), plot(1:1024,m1);grioln;title(' 均匀分布均值);xlabel(时间(t) '

15、;);ylabel幅 度(v)');sub plot(3,2,2), plot(1:1024,m2);gridn;title(' 高斯分布均值);xlabel(时间(t) ');ylabel幅 度(v)');subplot(3,2,3),plot(1:1024,m3);gridn;titleC指数分布均值);xlabel时间t)');ylabel幅度(v)');subplot(3,2,4),卩10妇:10244)妙加;世哄瑞利分布均值);xlabel时间t)');ylabel幅度(v)');subplot(3,2,5),卩10妇:

16、10245)妙01门;世哄卡方分布均值);xlabel时间t)');ylabel幅度(v)');%求各种分布的方差 figure;v1=var(x1),v2=var(x2),v3=var(x3),v4=var(x4),v5=var(x5);subplot(3,2,1),plot(1:1024,v1);groin;title('均匀分布方差);xlabel(时间t)');ylabel幅度');subplot(3,2,2),plot(1:1024,v2);groin;title('高斯分布方差);xlabel(时间t)');ylabel幅度&

17、#39;);subplot(3,2,3),plot(1:1024,v3);groin;title('指数分布方差);xlabel时间t)');ylabel幅度');subplot(3,2,4),plot(1:1024,v4);groin;title('瑞利分布方差);xlabel时间t)');ylabel幅度');subplot(3,2,5),plot(1:1024,v5);groin;title('卡方分布方差);xlabel(时间t)');ylabel幅度');%求各种分布的自相关函数 figure;title(

18、9;自相关函数图');x_c1 ,lags=xcorr(x1 ,200'u, nbiased)'x_c2,lags=xcorr(x2,200'u,nbiased)'x_c3,lags=xcorr(x3,200'u,nbia sed)'x_c4,lags=xcorr(x4,200'u, nbiased)'x_c5,lags=xcorr(x5,200'u,nbiased)'subplot(3,2,1),plotQagS'XCDQriDhtitleC均匀分布自相关函数图');subplot(3,2

19、,2),plot(lags,x_c2);groin;title('正态分布');subplot(3,2,3),plot(lags,x_c3);groin;title('指数分布');subplot(3,2,4),plotQagS'XCgrohtitleC瑞利分布');subplot(3,2,5),plot(lags,x_c5);groin;title('卡方分布');%求各种分布的概率密度函数 y1=unifpdf(x1,-1,1);y2=normpdf(x2,0,1);y3=exppdf(x3,1);y4=raylpdf(x4,

20、1);y5=chi2pdf(x5,1);%各种分布的概率密度估计 figure;k1,n1=ksdensity(x1);k2,n2=ksdensity(x2);k3,n3=ksdensity(x3);k4,n4=ksdensity(x4);k5,n5=ksdensity(x5);subplot(3,2,1),plot(n1,k1);gri(0n;title('均匀分布概率密度');xlabel(时间');ylabel(幅度') subplot(3,2,2),plot(n2,k2);gri(0n;title('正态分布');xlabel(时间

21、9;);ylabel幅度') subplot(3,2,3),plot(n3,k3);gri(on;title('指数分布');xlabel(时间');ylabel幅度(') subplot(3,2,4),plot(n4,k4);gri(on;title('瑞利分布');xlabel(时间');ylabel幅度') subplot(3,2,5),plot(n5,k5);gri(on;title('卡方分布');xlabel(时间');ylabel幅度') %功率谱密度 figure;f1=ff

22、t(x1,1024);f2=fft(x2,1024);f3=fft(x3,1024);f4=fft(x4,1024);迺豊)_eqe牙feo報罂)呂菸a二OOL 0 SOL oDs一xe-(电sqe)10_d-(eze)10_dqns迺豊_eqe_yc(-ZH 慝)_eqe_x( feo彩3L)呂話B二09 0 esL oDs一X2(dsqe)lo_d-(0ze)lo_dqns(ZH 慝)_eqe_x(峯期慝® 异)呂話B二09 0 SOL oDs一xe-(l4)sqe)10_d-(Lze)10_dqns豊)_eqe_y (-ZHi?慝)_eqe_XT豊)_eqe_yc(-ZHi?)

23、_eqe_XTgn64&&&&&&&&&&&&&&& 慕ffi峯期慝®&二OL 0 LLgoDS_xe&(mfeogep _岂耄 q(gd)sqe-寸 00L)10_d-(gze)10_dqns二OL 0 LLgoDS_xe&(mfeo 巨晖)邑耄 q (寸 d)sqe-寸sL)10_d-(寸 ze)10_dqns二OL 0 LLgoDS_xe&(m定0報罂互岂耄 q(ed)sqe-t70oL)10_d-(eze)10_dqns二OL

24、 0 LLgoDS_xe&(mq(0d)sqe-t70oL)lo_d-(0ze)lo_dqns二gVOL*OL gVOL9 lsods一xe&F迺冒互岂耄 q(Ld)sqe-寸sL)10_d-(Lze)10_dqns&s U(£)U8* £)USEHgd &s 17(寸汪U8* 寸)useh寸d-寸SIV(二U8*目ueEHed &00 v(du8*a)ueeEH0d &00 V(I4)U8* EueeEHLd于 0OL-gx)EH£V') subplot(3,2,4),plot(abs(f4),axis(0

25、 1023 0 50);goidtitle('瑞利分布');Xlabel(频率 Hz'); Ylabel(幅值V') subplot(3,2,5),plot(abs(f5),axis(0 1023 0 100);goidtitle('卡方分布');Xlabel(频率Hz'); Ylabel(幅值V') %随机数长度为 20480时 figure;x11=normrnd(2,3,1,20480%); 生成长度为 20480的正态分布 x_c11,lags=xcorr(x11,20480'u,nbiased)'f_c1

26、=fft(x11,20480);f_c2=fft(x11,10240);f_c3=fft(x11,1024);p11=mean(f_c1.*conj(f_c1)/20480;p22=mean(f_c2.*conj(f_c2)/10240;p33=mean(f_c3.*conj(f_c3)/1024;subplot(2,2,1);plot(1:20480,abs(p11);gion；title('N为20840正态分布功率谱密度');Xlabel(频率Hz');Ylabel(幅值') subplot(2,2,2);plot(1:10240,abs(p22);gion；ti