ip_02_04.m

来自「Communition MATLAB (Proakis)书中的源程序」· M 代码 · 共 17 行

% MATLAB script for Illustrative Problem 4, Chapter 2.
echo on
N=1000;
M=50;
Rx_av=zeros(1,M+1);
Sx_av=zeros(1,M+1);
for j=1:10,		    		% take the ensemble average over 10 realizations
   X=rand(1,N)-1/2;	      		% N i.i.d. uniformly distributed random variables
			      		% between -1/2 and 1/2
   Rx=Rx_est(X,M); 	      		% Autocorrelation of the realization
   Sx=fftshift(abs(fft(Rx))); 		% Power spectrum of the realization
   Rx_av=Rx_av+Rx;	      		% sum of the autocorrelations
   Sx_av=Sx_av+Sx;	      		% sum of the spectrums
end;
Rx_av=Rx_av/10;	              		% ensemble average autocorrelation
Sx_av=Sx_av/10;		      		% ensemble average spectrum
% Plotting comments follow