随机信号分析试验报告

1、HarbinInstituteofTechnology实验报告课程名称：随机信号分析院系：电子与信息工程学院班级：姓名：学号：指导教师：实验时间：实验一、各种分布随机数的产生1 .均匀分布随机数的产生原理产生伪随机数的一种实用方法是同余法，它利用同余运算递推产生伪随机数序列。最简单的方法是加同余法yn1=ync(modM)yn1_M为了保证产生的伪随机数能在0,1内均匀分布，需要M为正整数，此外常数c和初值y0亦为正整数。加同余法虽然简单，但产生的伪随机数效果不好。另一种同余法为乘同余法，它需要两次乘法才能产生一个0,1上均匀分布的随机数yn1yn1=ayn(modM)Xn1二M式中，a为正整

2、数。用加法和乘法完成递推运算的称为混合同余法，即yn1yn1=aync(modM)Xn.1=M用混合同余法产生的伪随机数具有较好的特性，一些程序库中都有成熟的程序供选择。常用的计算语言如Basic、C和Matlab都有产生均匀分布随机数的函数可以调用，只是用各种编程语言对应的函数产生的均匀分布随机数的范围不同，有的函数可能还需要提供种子或初始化。Matlab提供的函数rand()可以产生一个在0,1区间分布的随机数，rand(2,4)则可以产生一个在0,1区间分布的随机数矩阵，矩阵为2行4歹限Matlab提供的另一个产生随机数的函数是random('unif,a,b,N,M),unif

3、表示均匀分布，a和b是均匀分布区间的上下界，N和M分别是矩阵的行和列。2 .随机变量的仿真根据随机变量函数变换的原理，如果能将两个分布之间的函数关系用显式表达，那么就可以利用一种分布的随机变量通过变换得到另一种分布的随机变量。若X是分布函数为F(x)的随机变量，且分布函数F(x)为严格单调升函数，令Y=F(X),则Y必为在0,1上均匀分布的随机变量。反之，若Y是在0,1上均匀分布的随机变量，那么X=Fx(Y)即是分布函数为FX(x)的随机变量。式中FX1()为Fx()的反函数。这样，欲求某个分布的随机变量，先产生在0,1区间上的均匀分布随机数，再经上式变换，便可求得所需分布的随机数。3 .高斯

4、分布随机数的仿真广泛应用的有两种产生高斯随机数的方法，一种是变换法，一种是近似法。如果X1,X2是两个互相独立的均匀分布随机数，那么下式给出的Y1,Y2jYi=o；-2lnXicos(2ttX2)+m丫2=-J-2lnXsin(271X2)+m便是数学期望为m方差为1的高斯分布随机数，且互相独立，这就是变换法另外一种产生高斯随机数的方法是近似法。在学习中心极限定理时，曾提到n个在0,1区间上均匀分布的互相独立随机变量Xi（i=1,21）,当n足够大时，具和的分布接近高斯分布。当然，只要n不是无穷大，这个高斯分布是近似的。由于近似法避免了开方和三角函数运算，计算量大大降低。当精度要求不太高时，近

5、似法还是具有很大应用价值的。4 .各种分布随机数的仿真有了高斯随机变量的仿真方法，就可以构成与高斯变量有关的其他分布随机变量，如瑞利分布、指数分布和？2分布随机变量。（二）实验目的在很多系统仿真的过程中，需要产生不同分布的随机变量。利用计算机可以很方便地产生不同分布的随机变量，各种分布的随机变量的基础是均匀分布的随机变量。有了均匀分布的随机变量，就可以用函数变换等方法得到其他分布的随机变量。（三）实验结果高斯分布随机数-50052x10"/分布随机数附：源程序subplot(2,2,1);x=random('unif,2,5,1,1024);plot(x);t田e('

6、均匀分布随机数')subplot(2,2,2);G1=random('Normal',0,1,1,20000);plot(G1);title('高斯分布随机数')subplot(2,2,3);G2=random('Normal',0,1,1,20000);R=sqrt(G1.*G1+G2.*G2);plot(R);title('瑞利分布随机数')subplot(2,2,4);G3=random('Normal',0,1,1,20000);G4=random('Normal',0,1,1,20

7、000);X=G1.*G1+G2.*G2+G3.*G3+G4.*G4;plot(X);title('xA2分布随机数')实验二、随机变量检验1、均值的计算在实际计算时，如果平稳随机序列满足各态历经性，则统计均值可用时间均值代替。这样，在计算统计均值时，并不需要大量样本函数的集合，只需对一个样本函数求时间平均即可。甚至有时也不需要计算Ntg时的极限，况且也不可能。通常的做法是取一个有限的、计算系统能够承受的N求时间均值和时间方差。根据强调计算速度或精度的不同，可选择不同的算法。n设随机数序列xi,x2，,xn,一种计算均值的方法是直接计算下m=上£咒Nn式中，xn为随机

8、数序列中的第n个随机数。另一种方法是利用递推算法，第n次迭代的均值也亦即前n个随机数的均值为mnMnmnJxn=mn+1(xn-mn)迭代结束后，便得到随机数序列的均值nnnm=mN递推算法的优点是可以实时计算均值，这种方法常用在实时获取数据的场合。n当数据量较大时，为防止计算误差的积累，也可米用m=mi+£仅-mJ式Nn：中，mi是取一小部分随机数计算的均值。2、方差的计算计算方差也分为直接法和递推法。仿照均值的做法21,N221”22二-二不'(xn-m);u7y'xn-mNnziNnd方差的递推算法需要同时递推均值和方差1/、_2n-1p_21,、2“mn=mn

9、一(xn-mn二)；-n=；-n1一一mn/)nnn迭代结束后，得到随机数序列的方差为二-2=；n其它矩函数也可用类似的方法得到。3、统计随机数的概率密度直方图假定被统计的序列x(n)的最大值和最小值分别为a和b。将(a,b)区间等分M(M应与被统计的序列x(n)的个数N相适应，否则统计效果不好。)份后的区间小,(b-a)、,(b-a)2*(b-a)、为(a,a+-),(a+-,a+-),，MMM(b-a)(i-1)2*(b-a)*i(b-a)(M-1)(a.,a.)，(a.,b)。用f，表小MMM序列x(n)的值落在(a+(ba)(iNa+2*(ba)*i)区间里的个数，统计序列MMx(n)

10、的值在各个区间的个数f(i),i=0,2,，M-1,则f就粗略地反映了随机序列的概率密度的情况。用图形方式显示出来就是随机数的概率密度直方图。(二)实验目的随机数产生之后，必须对它的统计特性做严格的检验。一般来讲，统计特性的检验包括参数检验、均匀性检验和独立性检验等。事实上，我们如果在二阶矩范围内讨论随机信号，那么参数检验只对产生的随机数一、二阶矩进行检验。我们可以把产生的随机数序列作为一个随机变量，也可以看成随机过程中的一个样本函数。不论是随机变量还是随机过程的样本函数，都会遇到求其数字特征的情况，有时需要计算随机变量的概率密度直方图等。（三）实验结果3000502460均匀分布随机数直方图

11、1001500100050005高斯分布随机数直方图2000瑞利分布随机数直方图15002000100010000*20246/分布随机翻直方图2000r500°202040士Meanl土Mean2土MeanSFhMean4田田田田VariancelVariance2Variance3.5197-0.00711.25453.99370.77191.0056043158*0490附：源程序subplot(2,2,1);x=random('unif,2,5,1,1024);hist(x,2:0.2:5);title('均匀分布随机数直方图');s1=0forn1=

12、1:1024s1=x(n1)+s1;endMean1=s”1024;t1=0forn1=1:1024t1=(x(n1)-Mean1)A2+t1;endVariance1=t11024;subplot(2,2,2);G1=random('Normal',0,1,1,20000);hist(G1,-4:0.2:4);title('高斯分布随机数直方图');s2=0forn2=1:20000s2=G1(n2)+s2;endMean2=s220000;t2=0forn2=1:20000t2=(G1(n2)-Mean2)A2+t2;endVariance2=t220000

13、;subplot(2,2,3);G2=random('Normal',0,1,1,20000);R=sqrt(G1.*G1+G2.*G2);hist(R,0:0.2:5);title('瑞利分布随机数直方图');s3=0forn3=1:20000s3=R(n3)+s3;endMean3=s320000;t3=0forn3=1:20000t3=(R(n3)-Mean3)A2+t3;endVariance3=t320000;subplot(2,2,4);G3=random('Normal',0,1,1,20000);G4=random('No

14、rmal',0,1,1,20000);X=G1.*G1+G2.*G2+G3.*G3+G4.*G4;hist(X,0:0.5:30);title('xA2分布随机数直方图')s4=0forn4=1:20000s4=X(n4)+s4;endMean4=s420000;t4=0forn4=1:20000t4=(X(n4)-Mean4)A2+t4;end实验三、中心极限定理的验证（一）实验原理如果n个独立随机变量的分布是相同的，并且具有有限的数学期望和方差，当n无穷大时，它们之和的分布趋近于高斯分布。这就是中心极限定理中的一个定理。我们以均匀分布为例，来解释这个定理。若n个随机

15、变量Xi（i=1,2,n）都为0,1区间上的均匀分布的随机变量，且互相独立，当n足够大时，其和Y=gXi的分布接近高斯分布。i1（二）实验目的利用计算机产生均匀分布的随机数。对相互独立的均匀分布的随机变量做和，可以很直观看到均匀分布的随机变量的和，随着做和次数的增加分布情况的变化，通过实验对中心极限定理的进行验证。（三）实验结果40123均匀分布随机数直方图300200100五个均匀分布之和随机数直方图十个均匀分布之和随机数直方图高斯分布随机数直方图0分析：随n取值的增大，均匀分布随机序列求和的图形越发接近于高斯分布附：源程序X0=random('unif,0,1,1,1024);X1

16、=random('unif,0,1,1,1024);X2=random('unif,0,1,1,1024);X3=random('unif,0,1,1,1024);X4=random('unif,0,1,1,1024);X5=random('unif,0,1,1,1024);X6=random('unif,0,1,1,1024);X7=random('unif,0,1,1,1024);X8=random('unif,0,1,1,1024);X9=random('unif,0,1,1,1024);G=random('

17、normal',0,1,1,1024);Y1=X0+X1+X2+X3+X4;Y2=X0+X1+X2+X3+X4+X5+X6+X7+X8+X9;subplot(2,2,1);hist(X0,0:0.2:2);title('均匀分布随机数直方图')subplot(2,2,2);hist(Y1,0:0.2:6);title('五个均匀分布之和随机数直方图')subplot(2,2,3);hist(Y2,0:0.2:8);title('十个均匀分布之和随机数直方图')subplot(2,2,4);hist(G,-4:0.2:4);title(&#

18、39;高斯分布随机数直方图')(一)实验原理在实际应用中，我们可以把产生的随机数序列看成随机过程中的一个样本函数。如果平稳随机序列满足各态历经性，则统计自相关序列可用时间自相关序列代替。当数据的样本数有限时，也只能用有限个数据来估计时间自相关序列，统计自相关序列的估值。若各态历经序列X(n)的一个样本有N个数据x(0),x(1)，x(N-1),由于实序列自相关序列是对称的，自相关函数的估值为1N和|R(m)=N0刍X(n)X(+|m|)(二)实验目的在随机信号理论中，自相关函数是非常重要的概念。在实际系统仿真中也会经常计算自相关函数。通过本试验学生可以亲自动手计算自相关函数，加深对概念

19、的理解，并增强实际动手能力。(三)实验结果分析：分别生成均值为0和1,方差为1的高斯随机数，由图形可以明显看出两者自相关函数的差异。附：源程序N=256;xn=random('norm',0,1,1,N);Rx=xcorr(xn,'biased');m=-N+1:N-1;subplot(2,1,1);plot(m,Rx);title('均值为0,方差为1的高斯分布的自相关函数');axis(-NN-1-0.51.5);N=256;xn=random('norm',1,1,1,N);Xk=fft(xn,2*N);Rx=ifft(ab

20、s(Xk)A2)/N);m=-N:N-1;subplot(2,1,2);plot(m,fftshift(Rx);title('均值为1,方差为1的高斯分布的自相关函数');axis(-NN-1-0.51.5);实验五、功率谱密度(一)实验原理一般把平稳随机序列的功率谱定义为自相关序列的傅里叶变换。如果自相关X(n)的功率谱与自相关序列的关系为Sx()八Rx(m)emm-_：C,、1，C，、Rx(m)=Sx()e2*与实平稳过程一样，实平稳序列的功率谱也是非负偶函数，即Sx()-0Sx()=Sx(-)可以证明，功率谱还可表示为N2£x(n)e"|)n-_N当x

21、(n)为各态历经序列时，可去掉上式中的统计均值计算，将随机序列x(n)用它的一个样本序列x(n)代替。在实际应用中，由于一个样本序列的可用数据个数N有限，功率谱密度也只能是估计1N1Sxm唱x(n)i式中，x()是必0的傅里叶变换。这是比较简单的一种估计方法，这种功率谱密度的估计方法称为周期图方法。如果直接利用数据样本做离散傅里叶变换，可得到x()的离散值。由于这种方法可借助FFT算法实现，所以得到了广泛的应用。(二)实验目的在随机信号理论中，功率谱密度和自相关函数一样都是非常重要的概念。在实际系统仿真中也会经常计算。通过本试验学生可以亲自动手，加深对概念的理解，并增强实际动手能力。（三）实验

22、结果oo-4均值为0,方差为1的高斯分布的功率谱密度50100200250150f/Hz300均值为1,方差为1的高斯分布的功率清密度50100200250150(Hz300附：源程序N=256;x1=random('normal',0,1,1,N);Sx1=abs(fft(x1)A2)/N;subplot(2,1,1);plot(10*log10(Sx1);t田e('均值为0,方差为1的高斯分布的功率谱密度');xlabel('f/Hz')ylabel('SxldB')x2=random('normal',1,1

23、,1,N);Sx2=abs(fft(x2).A2)/N;subplot(2,1,2);plot(10*log10(Sx2);t田e('均值为1,方差为1的高斯分布的功率谱密度');xlabel('f/Hz')ylabel('Sx2dB')实验六、随机信号经过线性系统前后信号仿真（一）实验原理需要先仿真一个指定系统，再根据需要仿真输入的随机信号，然后使这个随机信号通过指定的系统。通过对实际系统建模，计算机可以对很多系统进行仿真。在信号处理中，一般将线性系统分解为一个全通放大器（或衰减器）和一个特定频率响应的滤波器。由于全通放大器可以用一个常数代替，

24、因此线性系统的仿真往往只需设计一个数字滤波器。滤波器设计可采用MATLA蜒供的函数，也可利用相应的方法自行设计。MATLA班供了多个设计滤波器的函数，可以很方便地设计低通、带通、高通、多带通、带阻滤波器。（二）实验目的系统仿真是信号仿真处理的一个重要部分，通过该实验要求学生掌握系统仿真的基本概念，并学会系统的仿真方法。（三）实验结果1、低通滤波器低通滤'波喜业E00.51mxn的功率谱密度4010.201初经低通滤波器的自相关函数0.3-0120001000010002000m乂n餐低通滤波富的功率谱密度20刈的自相关图数m乂n的功率谱密度50rO3、高通滤波器带通滤渡器刈经帚通通流波

25、器的自相关函数0.202-2000-1000010002000m2020高通滤波器刈经高通通滤'波器的自相关函数条带逋滤波器-0500.5111由经多带通通滤波器的目相关由教-2000-1000010002000xn经多帚通滤波器的功率诺密度5、带阻滤波器15带逋滤液器10.5-051乂n的功率谓密度刈经帚逋通滤波器的自相关函数0.202-2000-10000m100020002050omp35ano汾20Hz附：源程序1 .X(n)N=2000;fs=400;Nn=random('normal',0,1,1,N);t=(0:N-1)/fs;fi=random('

26、;unif,0,1,1,2)*2*pi;xn=sin(2*pi*50*t+fi(1)+Nn;Rx=xcorr(xn,'biased');m=-N+1:N-1;Sx=abs(fft(xn).A2)/N;f=(-N/2:N/2-1)*fs/N;subplot(211),plot(m,Rx);xlabel('m')ylabel('Rx(m)')title('xn的自相关函数');subplot(212),plot(f,fftshift(10*log10(Sx(1:N);xlabel('f/Hz')ylabel('

27、Sx/dB')title('xn的功率谱密度');2 .低通滤波器h=fir1(100,0.4);H=fft(h,2*N);HW=abs(H).A2;Rx=xcorr(xn,'biased');Sx=abs(fftshift(fft(xn,2*N).A2)/(2*N);Sy=Sx.*HW;Ry=fftshift(ifft(Sy);f=(-N:N-1)*fs/(2*N);m=(-N:N-1);subplot(311);plot(-N:N-1)/N,fftshift(abs(HW(1:2*N);title('低通滤波器');subplot(3

28、12),plot(m,Ry);xlabel('m')ylabel('Ry(m)')title('xn经低通滤波器的自相关函数');subplot(313),plot(f,fftshift(10*log10(Sy(1:2*N);axis(-200200-2020);xlabel('f/Hz')ylabel('Sy/dB')title('xn经低通滤波器的功率谱密度');3 .带通滤波器h=fir1(100,0.10.5);H=fft(h,2*N);HW=abs(H)A2;Rx=xcorr(xn,

29、9;biased');Sx=abs(fftshift(fft(xn,2*N).A2)/(2*N);Sy=Sx.*HW;Ry=fftshift(ifft(Sy);f=(-N:N-1)*fs/(2*N);m=(-N:N-1);subplot(311);plot(-N:N-1)/N,fftshift(abs(HW(1:2*N);title('带通滤波器');subplot(312),plot(m,Ry);xlabel('m')ylabel('Ry(m)')title('xn经带通通滤波器的自相关函数');subplot(313)

30、,plot(f,fftshift(10*log10(Sy(1:2*N);axis(-200200-2020);xlabel('f/Hz')ylabel('Sy/dB')title('xn经带通滤波器的功率谱密度');4 .高通滤波器h=fir1(100,0.6,'high');H=fft(h,2*N);HW=abs(H).A2;Rx=xcorr(xn,'biased');Sx=abs(fftshift(fft(xn,2*N).A2)/(2*N);Sy=Sx.*HW;Ry=fftshift(ifft(Sy);f=(-

31、N:N-1)*fs/(2*N);m=(-N:N-1);subplot(311);plot(-N:N-1)/N,fftshift(abs(HW(1:2*N);title('高通滤波器');subplot(312),plot(m,Ry);xlabel('m')ylabel('Ry(m)')title('xn经高通通滤波器的自相关函数');subplot(313),plot(f,fftshift(10*log10(Sy(1:2*N);axis(-200200-2020);xlabel('f/Hz')ylabel('Sy/dB')title('xn经高通滤波器的功率