untitled604.m

来自「这部分是用MATLAB产生随机过程的一系列源代码.含产生高斯分布随机变量,多变量」· M 代码 · 共 21 行

N=1000;
M=50;
Rx_av=zeros(1,M+1);
Sx_av=zeros(1,M+1);
for j=1:10   % 取10次实现的整体平均
    X=rand(1,N)-1/2;   % N i.i.d.均匀分布随机变量在-1/2和1/2间
    Rx=Rx_est(X,M);    % 本次实现的自相关函数
    Sx=fftshift(abs(fft(Rx))); % 本次实现的功率谱
    Rx_av=Rx_av+Rx;  %  自相关系数之和
    Sx_av=Sx_av+Sx;   % 功率谱之和
end;
Rx_av=Rx_av/10;     %  自相关系数的整体平均
Sx_av=Sx_av/10;     %  功率谱的整体平均
n=[0:1:50];
f=[0:1:50];
subplot(2,1,1);
plot(n,Rx_av);
title('自相关函数')
subplot(2,1,2);
plot(f,Sx_av);
title('功率谱')