资源描述：
《系统建模与仿真实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、系统建模与仿真实验报告报告一：产生10中独立分布的随机数，并检验其一、二阶距的性质。1、[0,1]区间的均匀分布采用乘同余法产生均匀分布在（0,1）之间的随机数。乘同余法的递推公式为：一般情况下，，a为整数，M于计算机的字长有关，,m为16或32，x的初值为。利用MATLAB实现，代码如下：functionu=undistribution(a,b,m);%乘同余法lam=8*a-3;M=pow2(m);x(1)=pow2(b)+1;fori=2:10000;y=lam*x(i-1);x(i)=mod(y,M);endu=x/M;end调用函数，并检验产生随机数的数字特性。y=undist
2、ribution(3,2,32);hist(y,50);E=mean(y);D=var(y);title('0-1均匀分布直方图');text(0,-20,strcat('均值为',num2str(E)));text(0.77,-20,strcat('均值为',num2str(D)));2、标准正态分布高斯分布的概率密度函数：；首先利用前面产生均匀分布随机数的方法生成两组均匀分布的随机数u1,u2；利用公式：，Z服从高斯分布。MATLAB实现代码如下：u1=undistribution(3,2,32);u2=undistribution(2,3,32);z=sqrt(-2*log(u1)
3、).*cos(2*pi*u2);hist(z,100);E=mean(z);D=var(z);title('标准正态分布直方图')text(-6,-40,strcat('均值为',num2str(E)));text(6,-40,strcat('方差为',num2str(D)));3、指数分布指数分布的概率密度函数如下：；首先利用前面产生均匀分布随机数的方法生成一组均匀分布的随机数u；则数列，为均值为，方差为的指数分布随机数列。MATLAB实现代码如下：u=undistribution(3,2,32);lam=8;y=-log(u)/lam;hist(y,150);E=mean(y);D=
4、var(y);title('指数分布直方图')text(0,-130,strcat('均值为',num2str(E)));text(2.3,-130,strcat('方差为',num2str(D)));4、广义指数分布广义指数分布的概率密度函数如下：，设有两组高斯分布的数列x、y，另，其中s为信噪比，则服从广义指数分布。MATLAB实现代码入下：u1=undistribution(3,2,32);u2=undistribution(2,3,32);x=sqrt(-2*log(u1)).*cos(2*pi*u2);y=sqrt(-2*log(u1)).*sin(2*pi*u2);s=2;x
5、1=x+sqrt(2*s);z=x1.^2+y.^2;hist(z,100)E=mean(z);D=var(z);title('广义指数分布直方图');text(0,-60,strcat('均值为',num2str(E)));text(63,-60,strcat('方差为',num2str(D)));5、瑞利分布瑞利分布的概率密度函数为：；瑞利分布数组可以利用高斯分布数列来产生，设有高斯分布的数列x、y，则服从瑞利分布。利用MATLAB实现代码如下：u1=undistribution(3,2,32);u2=undistribution(2,3,32);x=sqrt(-2*log(u1))
6、.*cos(2*pi*u2);y=sqrt(-2*log(u1)).*sin(2*pi*u2);z=sqrt(x.^2+y.^2);hist(z,50)E=mean(z);D=var(z);title('瑞利分布直方图');text(0,-65,strcat('均值为',num2str(E)));text(5.5,-65,strcat('方差为',num2str(D)));6、广义瑞利分布广义瑞利分布的概率密度函数如下：；同样，可以利用产生瑞利分布的方法来生产广义瑞利分布数列。设有高斯分布数列x、y，则服从广义瑞利分布，其中a为常数。MATLAB实现代码如下：u1=undistribut
7、ion(3,2,32);u2=undistribution(2,3,32);x=sqrt(-2*log(u1)).*cos(2*pi*u2);y=sqrt(-2*log(u1)).*sin(2*pi*u2);a=2;x1=x+a;z=sqrt(x1.^2+y.^2);hist(z,100)E=mean(z);D=var(z);title('广义瑞利分布直方图');text(0,-30,strcat('均值为',num2str(E))
