版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、 应用统计与随机课程 课程实验报告 题目: 相关正态随机过程的仿真 实验1 相关正态随机过程的仿真实验目的: 以正态随机过程为例,掌握离散时间随机过程的仿真方法,理解正态分布随机过程与均匀分布随机过程之间的相互关系,理解随机过程的相关函数等数值特征;培养计算机编程能力。实验内容:程序代码: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)
2、bar(n2)实验结果:结果分析:因为两个独立序列是随机产生,且在0,1均匀分布,故将0,1分为十个等宽区间时,落在每个区间的数目应该大致相等。实验内容:程序代码: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)实验结果:结果分析:绘制出的图形符合白色正态分布实验内容: 程序代码:clc;%第(3)题u1=rand(1,100000);u2=
3、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函数绘制分布直方图%第(4)题sum=0;for i=1:100000-1 sum=sum+x(i);%-表示x(n)的1到100000项的累加和endmx=sum/1
4、00000%-算出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);%-表示集合统计的方法计算出来的相关函数 endr%-算出r的值%第(5)问%先算出在各区间上数据出现的比例num1=0;num2=0;num3=0;num4=0;for i=1:1:100000 if (x(i)
5、<-2) num1=num1+1; else if (x(i)>=-2)&(x(i)<=0) num2=num2+1; else if (x(i)>0)&(x(i)<=2) num3=num3+1; else num4=num4+1; end end endend disp('实验值为')p1=num1/100000p2=num2/100000p3=num3/100000p4=num4/100000%再算出理论值p2=0;for i=1:200000 p2=p2+1/(sqrt(2*pi)*2)*exp(-(i*0.00001)*(i*0.00001)/(2*2*2)*0.00001;endp3=p2;p1=(1-2*p2)/2;p4=p1;disp('理想值为')p1,p2,p3,p4实验结果: 结果分析: 随着k值增大,相关性越弱,故r值越小,(4)(5)问理论值与理想值基本一致。 实验体会: 在本次试验的开始,因为对MATLAB的不熟练,不知道如何编程,经过查阅资料,了解了相关函数的使用方法之后,成功做完了本次实验。 通过本次试验,我掌握了
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年南京客运考试口诀是什么意思
- 2024年郴州考从业资格证客运试题
- 2024年安徽烟台客运上岗证考试题
- 2024年西藏驾校客运从业资格证模拟考试题
- 抗病毒药物行业竞争格局与投资战略研究咨询报告
- 2024年南京客运资格证专业知识试题及答案
- 2024年阜阳客运资格证仿真试题
- 在线AI教育行业投资机会分析与策略研究报告
- 企业风险保险行业发展趋势及前景展望分析报告
- 抗感染药物行业消费者群体特征分析
- 2020教科版新教材小学科学三年级上册《天气》单元教材解读
- 各类地质特征描述重点
- (完整版)眼部药物代谢动力学的研究
- 50303《建筑电气工程施工质量验收规范》常用表格
- 公路、电力线路交叉跨越一般施工方法
- 矩形灯杆强度挠度计算书3-以10米灯杆为例
- (完整版)牛津树单词卡(可放大彩打)
- 企业共享用工模式分析
- 浅谈锅炉烟气及风量的测量
- 《田家四季歌》优质课件
- 《用人单位内部劳动保障规章制度》
评论
0/150
提交评论