ip_02_04.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:1k
- % MATLAB script for Illustrative Problem 2.4.
- 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 ten 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
- echo off ;
- end;
- echo on ;
- Rx_av=Rx_av/10; % ensemble average autocorrelation
- Sx_av=Sx_av/10; % ensemble average spectrum
- % Plotting comments follow