echo off 
clear all 
close all 
N=1000; 
M=50; 
Rx_av=zeros(1,M+1); 
Sx_av=zeros(1,M+1); 
x=linspace(-25,25,M+1); 
for j=1:10, > 取10次實(shí)現(xiàn)到整體平均 
> 產(chǎn)生N 維均勻分布隨機(jī)變量，在-1/2和1/2之間 
X=rand(1,N)-1/2; 
Rx=Rx_est(X,M); > 本次實(shí)現(xiàn)的自相關(guān)系數(shù) 
Sx=fftshift(abs(fft(Rx))); > 本次實(shí)現(xiàn)到功率譜 
Rx_av=Rx_av+Rx; > 子相關(guān)系數(shù)之和 
Sx_av=Sx_av+Sx; > 功率譜之和 
end; 
Rx_av=Rx_av/10; > 自相關(guān)系數(shù)的整體平均 
Sx_av=Sx_av/10; > 功率譜到整體平均 
> Plotting comments follow 
subplot(2,1,1) 
plot(x,fftshift(Rx_av)) 
grid on 
title('autocorrelation')>畫自相關(guān)函數(shù)圖 
subplot(2,1,2) 
plot(x,fftshift(Sx_av)) >畫功率譜 
grid on 
title('spectrum')