利用MATLAB繪制隨機相位正弦波

實驗二利用MATLAB繪制隨機相位正弦波X(t)=asin(wt+0)的均值，方差和自相關(guān)函數(shù)的圖像［實驗?zāi)康模萃ㄟ^繪制圖像，深入理解隨機相位正弦波的均值，方差和自相關(guān)函數(shù)。［實驗學(xué)時］課外完成［實驗準(zhǔn)備］1．熟悉隨機相位正弦波的均值，方差和自相關(guān)函數(shù)含義，掌握其計算方法；2．編寫隨機相位正弦波的均值，方差和自相關(guān)函數(shù)的書面程序代碼。［實驗內(nèi)容及步驟］繪制隨機相位正弦波X(t)=asin(wt+0)的均值，方差和自相關(guān)函數(shù)的圖像說明：根據(jù)課本中計算隨機過程數(shù)字特征的方法來確定X(t)=asin(wt+0)的均值，方差和自相關(guān)函數(shù)的取值；確定上述均值，方差和自相關(guān)函數(shù)繪圖的坐標(biāo)范圍；利用MATLAB語言編寫繪圖程序，并進行調(diào)試；(4)運行實驗結(jié)果，編寫本次實驗的實驗報告。設(shè)定該正弦函數(shù)的隨機相位。在(-n,n)服從均勻:

X(t)=asin(wt+0)已知條件。的概率密度為:f(x)=,—n<%<n,0,其他其均值為Ex二J*asin(wt+0)—^d。_n 2n=0其方差為Dx=E(x2)-E(x)2二aiE(l—cos(23t+2e))2二a22Matlab的操作源代碼如下：均值函數(shù)clear;swsita;t=[-10:0.1:10];a=3;s=[-5:0.1:5];Ex=a/(2*pi)*int(sin(w*t+sita),sita,-pi,pi);subplot(2,2,1);plot(t,double(s(Ex)));xlabel('t:-10to10','FontSize',16);ylabel('y=Ex','FontSize',16);title('相位隨機變量的均值函數(shù)圖');圖像：相位隨機變量的均值函數(shù)圖方差函數(shù)clear;swsita;t=[-10:0.1:10];a=3;s=[-5:0.1:5];Dx=aA2-aA2/(2*pi)*int(cos(2*w*t+2*sita),sita,-pi,pi);subplot(2,2,2);plot(t,double(s(Dx)));xlabel('t:-10to10','FontSize',16);ylabel('y=Dx','FontSize',16);title('相位隨機變量的方差函數(shù)圖');圖像：相位隨機變量的方差函數(shù)圖

109JI0.5x95自相關(guān)函數(shù)109JI0.5x95自相關(guān)函數(shù)8 ■ ■ ■ -10 -5 0 5 10t:-10to10Symsstsita;a=3;w=pi/20;Rx=aA2/(2*pi)*int(sin(w*t+sita).*sin(w*s+sita),sita,-pi,pi);ezsurf(Rx);ezsurf(t,s,Rx,[-