哈工大随机信号实验报告..docVIP

Harbin Institute of Technology 实验报告 课程名称： 随机信号分析 院 系： 电信学院 班 级： 姓 名：哈尔滨工业大学 实验一 各种分布随机数的产生 实验目的 在很多系统仿真的过程中，需要产生不同分布的随机变量。利用计算机可以很方便地产生不同分布的随机变量，各种分布的随机变量的基础是均匀分布的随机变量。有了均匀分布的随机变量，就可以用函数变换等方法得到其他分布的随机变量。 实验内容 产生均匀分布的随机数、高斯分布的随机数和其它分布的随机数。 实验原理 均匀分布随机数的产生原理 产生伪随机数的一种实用方法是同余法，它利用同余运算递推产生伪随机数序列。最简单的方法是加同余法 为了保证产生的伪随机数能在[0,1]内均匀分布，需要M为正整数，此外常数c和初值y0亦为正整数。加同余法虽然简单，但产生的伪随机数效果不好。另一种同余法为乘同余法，它需要两次乘法才能产生一个[0,1]上均匀分布的随机数 式中，a为正整数。用加法和乘法完成递推运算的称为混合同余法，即 用混合同余法产生的伪随机数具有较好的特性，一些程序库中都有成熟的程序供选择。 常用的计算语言如Basic、C和Matlab都有产生均匀分布随机数的函数可以调用，只是用各种编程语言对应的函数产生的均匀分布随机数的范围不同，有的函数可能还需要提供种子或初始化。 Matlab提供的函数rand()可以产生一个在[0,1]区间分布的随机数，rand(2,4)则可以产生一个在[0,1]区间分布的随机数矩阵，矩阵为2行4列。Matlab提供的另一个产生随机数的函数是random(unif,a,b,N,M)，unif 表示均匀分布，a和b是均匀分布区间的上下界，N和M分别是矩阵的行和列。 随机变量的仿真 根据随机变量函数变换的原理，如果能将两个分布之间的函数关系用显式表达，那么就可以利用一种分布的随机变量通过变换得到另一种分布的随机变量。 若X是分布函数为F(x)的随机变量，且分布函数F(x)为严格单调升函数，令Y=F(X)，则Y必为在[0,1]上均匀分布的随机变量。反之，若Y是在[0,1]上均匀分布的随机变量，那么 即是分布函数为FX(x)的随机变量。式中为的反函数。这样，欲求某个分布的随机变量，先产生在[0,1]区间上的均匀分布随机数，再经上式变换，便可求得所需分布的随机数。 高斯分布随机数的仿真 广泛应用的有两种产生高斯随机数的方法，一种是变换法，一种是近似法。 如果X1,X2是两个互相独立的均匀分布随机数，那么下式给出的Y1,Y2 便是数学期望为m，方差为??的高斯分布随机数，且互相独立，这就是变换法。 另外一种产生高斯随机数的方法是近似法。在学习中心极限定理时，曾提到n个在[0,1]区间上均匀分布的互相独立随机变量Xi (i=1,2…,n)，当n足够大时，其和的分布接近高斯分布。当然，只要n不是无穷大，这个高斯分布是近似的。由于近似法避免了开方和三角函数运算，计算量大大降低。当精度要求不太高时，近似法还是具有很大应用价值的。 各种分布随机数的仿真 有了高斯随机变量的仿真方法，就可以构成与高斯变量有关的其他分布随机变量，如瑞利分布、指数分布和?分布随机变量。 实验过程和结果分析 1、均匀分布、高斯分布随机数的产生与仿真 思路：利用已知matlab函数直接产生随机数。 程序 x=random(unif,5,10,1,1000);%产生1000个服从于（5,10）的随机数 y=random(normal,0,1,1,3000);%产生3000个服从于N（0,1）的随机数 subplot(211),plot(x);title(均匀分布随机数) subplot(212),plot(y);title(高斯分布随机数) 分析：产生的随机数呈现中间多，两头少的趋势，普遍集中于期望附近。 2、瑞利分布、指数分布及分布随机数的产生与仿真 ①　思路：利用已知matlab函数的变换加和产生随机数。 ②　程序 N=5000; G1=random(Normal,0,1,1,N); G2=random(Normal,0,1,1,N); G3=random(Normal,0,1,1,N); G4=random(Normal,0,1,1,N); R=sqrt(G1.*G1+G2.*G2); E=G1.*G1+G2.*G2; X=G1.*G1+G2.*G2+G3.*G3+G4.*G4; subplot(311);plot(R);title(瑞利分布随机数) subplot(312)