相關(guān)正態(tài)隨機(jī)過(guò)程的仿真

1、 應(yīng)用統(tǒng)計(jì)與隨機(jī)課程 課程實(shí)驗(yàn)報(bào)告 題目： 相關(guān)正態(tài)隨機(jī)過(guò)程的仿真 實(shí)驗(yàn)1 相關(guān)正態(tài)隨機(jī)過(guò)程的仿真實(shí)驗(yàn)?zāi)康模?以正態(tài)隨機(jī)過(guò)程為例，掌握離散時(shí)間隨機(jī)過(guò)程的仿真方法，理解正態(tài)分布隨機(jī)過(guò)程與均勻分布隨機(jī)過(guò)程之間的相互關(guān)系，理解隨機(jī)過(guò)程的相關(guān)函數(shù)等數(shù)值特征;培養(yǎng)計(jì)算機(jī)編程能力。實(shí)驗(yàn)內(nèi)容：程序代碼：u1=rand(1,100000);u2=rand(1,100000);%-在0,1區(qū)間用rand函數(shù)生成兩個(gè)相互獨(dú)立的是隨機(jī)序列n1=hist(u1,10)%-用hist函數(shù)繪制分布直方圖subplot(121)%-將兩幅分布圖顯示在一個(gè)窗口bar(n1)n2=hist(u2,10)subplot(122)

2、bar(n2)實(shí)驗(yàn)結(jié)果：結(jié)果分析：因?yàn)閮蓚€(gè)獨(dú)立序列是隨機(jī)產(chǎn)生，且在0,1均勻分布，故將0,1分為十個(gè)等寬區(qū)間時(shí)，落在每個(gè)區(qū)間的數(shù)目應(yīng)該大致相等。實(shí)驗(yàn)內(nèi)容：程序代碼：clc;u1=rand(1,100000);u2=rand(1,100000);%-在0,1區(qū)間用rand函數(shù)生成兩個(gè)相互獨(dú)立的是隨機(jī)序列en=sqrt(-2*log(u1).*cos(2*pi*u2);-定義白色正態(tài)分布e(n)n=hist(en,100);%-用hist函數(shù)繪制分布直方圖bar(n)實(shí)驗(yàn)結(jié)果:結(jié)果分析：繪制出的圖形符合白色正態(tài)分布實(shí)驗(yàn)內(nèi)容： 程序代碼：clc;%第（3）題u1=rand(1,100000);u2=

3、rand(1,100000);%-在0,1區(qū)間用rand函數(shù)生成兩個(gè)相互獨(dú)立的是隨機(jī)序列en=sqrt(-2*log(u1).*cos(2*pi*u2);%-定義白色正態(tài)分布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 %-生成隨機(jī)過(guò)程x(n)hist(x,100)%-用hist函數(shù)繪制分布直方圖%第（4）題sum=0;for i=1:100000-1 sum=sum+x(i);%-表示x(n)的1到100000項(xiàng)的累加和endmx=sum/1

4、00000%-算出mx的值for i=1:100000-1 sum=sum+x(i)*x(i);%-表示x(n)*x(n)的1到100000項(xiàng)的累加和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);%-表示集合統(tǒng)計(jì)的方法計(jì)算出來(lái)的相關(guān)函數(shù) endr%-算出r的值%第(5)問(wèn)%先算出在各區(qū)間上數(shù)據(jù)出現(xiàn)的比例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('實(shí)驗(yàn)值為')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實(shí)驗(yàn)結(jié)果： 結(jié)果分析： 隨著k值增大，相關(guān)性越弱，故r值越小，（4）（5）問(wèn)理論值與理想值基本一致。 實(shí)驗(yàn)體會(huì)： 在本次試驗(yàn)的開(kāi)始，因?yàn)閷?duì)MATLAB的不熟練，不知道如何編程，經(jīng)過(guò)查閱資料，了解了相關(guān)函數(shù)的使用方法之后，成功做完了本次實(shí)驗(yàn)。 通過(guò)本次試驗(yàn)，我掌握了