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

实验一 相关正态分布离散随机过程的产生 实验目的 以正态随机过程为例，掌握离散时间随机过程的仿真方法，理解正态分布随机过程与均匀分布随机过程之间的相互关系，理解随机过程的相关函数等数值特征;培养计算机编程能力。 实验要求 利用计算机语言的[0,1]区间均匀分布随机数产生函数生成两个相互独立的序列 {U1(n)|n=1,2,…100000}，{U2(n)|n=1,2,…100000} 生成均值为m=0，根方差=1的白色正态分布序列 {e(n)|n=1,2,…，100000} 假设离散随机过程x(n)服从均值为=0、根方差为、相关函数为 功率谱函数为 随机过程x(n)的生成方法为 (n=1,2,…100000) 给定初始条件x(0)=0 采用集合统计的方法计算 验证计算出来的统计参数与理论值是否一致，差异大小。 采用计算机程序计算正态分布的区间积分 根据已生成的序列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）的一致性。 实验代码及结果 利用计算机语言的[0,1]区间均匀分布随机数产生函数生成两个相互独立的序列 {U1(n)|n=1,2,…100000}，{U2(n)|n=1,2,…100000} 代码： u1=rand(1,100000); u2=rand(1,100000); subplot(1,2,1); hist(u1); subplot(1,2,2); hist(u2); 实验结果： 分析:利用随机函数产生了两个随机序列,区间为[0,1]。hist函数默认将区间划分为10等份。 生成均值为m=0，根方差=1的白色正态分布序列 {e(n)|n=1,2,…，100000} 代码： clc; u1=rand(1,100000); u2=rand(1,100000); e=sqrt(-2*log(u1)).*cos(2*pi*u2); hist(e,100); 实验结果： 假设离散随机过程x(n)服从均值为=0、根方差为、相关函数为 功率谱函数为 随机过程x(n)的生成方法为 (n=1,2,…100000) 给定初始条件x(0)=0 代码： clc; u1=rand(1,100000); u2=rand(1,100000); en=sqrt(-2*log(u1)).*cos(2*pi*u2); 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).*e(n+1); end hist(x,100); 实验结果： 分析：生成服从均值为=0、根方差为、相关函数为 的离散随机过程x(n)。 采用集合统计的方法计算 验证计算出来的统计参数与理论值是否一致，差异大小。 代码： u1=rand(1,100000); u2=rand(1,100000); en=sqrt(-2*log(u1)).*cos(2*pi*u2); 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 sum=0; for i=1:100000-1 sum=sum+x(i); end mx=sum/100000 for i=1:100000-1 sum=sum+x(i)*x(i); end ax=sqrt(sum/100000) for k=1:4 sum=0; for j=1:100000-k sum=sum+x(j)*x(j+k); end r(k)=sum/(100000-k); end r 输出结果： 采用计算机程序计算正态分布的区间积分 根据已生成的序列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）的一致性。 代码： num1=0;num2=0;num3=0;num4=0; for i=1:1