应用统计与随机过程实验报告实验一.docx

《应用统计与随机过程实验报告实验一.docx》由会员分享，可在线阅读，更多相关《应用统计与随机过程实验报告实验一.docx（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、HUNAN UNIVERSITY应用统计学与随机过程实验报告学生姓名 阳教授 学生学号 201308030511 专业班级 通信工程1305班 实验题目 相关正态分布离散随机过程的产生实验一 相关正态分布离散随机过程的产生一、 实验目的 以正态随机过程为例，掌握离散时间随机过程的仿真方法，理解正态分布随机过程与均匀分布随机过程之间的相互关系，理解随机过程的相关函数等数值特征;培养计算机编程能力。二、 实验要求1) 利用计算机语言的0,1区间均匀分布随机数产生函数生成两个相互独立的序列 U1(n)|n=1,2,100000，U2(n)|n=1,2,1000002) 生成均值为m=0，根方差=1的

2、白色正态分布序列 e(n)|n=1,2,，100000 3) 假设离散随机过程x(n)服从均值为=0、根方差为、相关函数为 功率谱函数为 随机过程x(n)的生成方法为 (n=1,2,100000) 给定初始条件x(0)=04) 采用集合统计的方法计算 验证计算出来的统计参数与理论值是否一致，差异大小。5) 采用计算机程序计算正态分布的区间积分 根据已生成的序列x(n)，在100000个数据中，分别计算（-， -2），-2，0，(0，2，2，)区间上数据出现的比例P1，P2，P3，P4。比较P1，P2，P3，P4与理想值（0.5-P）,P,P,（0.5-P）的一致性。三、 实验代码及结果1. 利

3、用计算机语言的0,1区间均匀分布随机数产生函数生成两个相互独立的序列U1(n)|n=1,2,100000，U2(n)|n=1,2,100000代码：u1=rand(1,100000); u2=rand(1,100000);%-在0,1区间用rand函数生成两个相互独立的是随机序列n1=hist(u1,10)%-用hist函数绘制分布直方图subplot(121)%-将两幅分布图显示在一个窗口bar(n1)n2=hist(u2,10)subplot(122)bar(n2)实验结果：分析:利用随机函数产生了两个随机序列,区间为0,1。hist函数默认将区间划分为10等份。2. 生成均值为m=0，根

4、方差=1的白色正态分布序列 e(n)|n=1,2,，100000 代码：clc; u1=rand(1,100000); u2=rand(1,100000);%-在0,1区间用rand函数生成两个相互独立的是随机序列 en=sqrt(-2*log(u1).*cos(2*pi*u2);%-定义白色正态分布e(n) n=hist(en,100);%-用hist函数绘制分布直方图 bar(n)实验结果：3. 假设离散随机过程x(n)服从均值为=0、根方差为、相关函数为 功率谱函数为 随机过程x(n)的生成方法为 (n=1,2,100000) 给定初始条件x(0)=0代码：clc; %第（3）题u1=r

5、and(1,100000);u2=rand(1,100000);%-在0,1区间用rand函数生成两个相互独立的是随机序列en=sqrt(-2*log(u1).*cos(2*pi*u2);%-定义白色正态分布e(n)a=0.6; x(1)=2*sqrt(1-a*a)*en(1);%-初始化for n=1:100000-1; x(n+1)=a*x(n)+2*sqrt(1-a*a).*en(n+1);end %-生成随机过程x(n)hist(x,100)%-用hist函数绘制分布直方图实验结果：分析：生成服从均值为=0、根方差为、相关函数为 的离散随机过程x(n)。4. 采用集合统计的方法计算 验

6、证计算出来的统计参数与理论值是否一致，差异大小。代码：clc; %第（4）题sum=0;for i=1:100000-1 sum=sum+x(i);%-表示x(n)的1到100000项的累加和endmx=sum/100000%-算出mx的值for i=1:100000-1 sum=sum+x(i)*x(i); %-表示x(n)*x(n)的1到100000项的累加和endax=sqrt(sum/100000);%-算出ax的值for k=1:4 sum=0;%-始化sum的值 for j=1:100000-k sum=sum+x(j)*x(j+k); end r(k)=sum/(100000-k

7、);%-表示集合统计的方法计算出来的相关函数endr%-算出r的值输出结果：5. 采用计算机程序计算正态分布的区间积分 根据已生成的序列x(n)，在100000个数据中，分别计算（-，-2），-2，0，(0，2，2，)区间上数据出现的比例P1，P2，P3，P4。比较P1，P2，P3，P4与理想值（0.5-P）,P,P,（0.5-P）的一致性。代码：clc;%第(5)问%先算出在各区间上数据出现的比例num1=0;num2=0;num3=0;num4=0;for i=1:1:100000 if (x(i)=-2)&(x(i)0)&(x(i)分析：通过将积分运算转化为小区间内的值的相加，可以得到p1,p2,p3,p4的实验值，与理想值进行对比相差不大。四、实验体会： 通过本次试验，我掌握了离散时间随机过程的仿真方法，进一步理解了课堂上学习正态分布随机过程与均匀分布随机过程之间的相互关系，理解随机过程的相关函数等数值特征，并且学会了相关函数的使用，培养计算机编程能力。