能不能帮忙注解一下下面的程序，再讲解一下每一段作用？

Question

%===============================产生信号===================================
M=1000;  %快拍数
thetas=[5,30]*pi/180;    %信号来波方向
Ns=length(thetas);%信号源个数2
for n=1:Ns
    SS(n,1:M)=randn(1,M)+j*randn(1,M);
    SS(n,1:M)=SS(n,1:M)/sqrt(SS(n,1:M)*SS(n,1:M)'/M);
    A(1:Ne,n)=exp(j*[0:Ne-1].'*beta*d*sin(thetas(n)));
end
%================================Monte-Carlo===============================
SNR=-20:1:20;
monte=100;     %Monte-Carlo模拟的次数
kk=0;
for nsnr=1:length(SNR)
   
    Ns_est=0;
    kk=kk+1;
    correct1=0;
    correct2=0;
    for mk=1:monte
         Renoise=randn(Ne,M);
        %=========产生零均值,方差为1的高斯白噪声==========
        for n=1:Ne
            Renoise(n,:)=Renoise(n,:)./std(Renoise(n,:));
            Renoise(n,:)=Renoise(n,:)-mean(Renoise(n,:));
          
        end
        %================================================
        Ps=10.^(SNR(nsnr)./10);     %信号功率
        S=sqrt(Ps).*SS;     
        Y=A*S+Renoise;    %接收到的信号
       R=Y*Y'./M;   %计算协方差矩阵

Answer 1

%===============================产生信号===================================
M=1000;  %快拍数，即信号样本的数量
thetas=[5,30]*pi/180;    %信号来波方向，以弧度为单位
Ns=length(thetas);%信号源个数2
for n=1:Ns
    SS(n,1:M)=randn(1,M)+j*randn(1,M);  %生成复数高斯随机信号
    SS(n,1:M)=SS(n,1:M)/sqrt(SS(n,1:M)*SS(n,1:M)'/M);  %归一化信号，使其具有单位能量
    A(1:Ne,n)=exp(j*[0:Ne-1].'*beta*d*sin(thetas(n)));  %计算方向向量，这里beta, d, Ne未定义，需要补充
end

这段代码的目的是生成信号源，并为每个信号源生成一个归一化的复数高斯信号。同时，它还计算了与信号来波方向相关的方向向量，但需要额外的参数beta, d, Ne。

%===============================================Monte-Carlo===============================

SNR=-20:1:20;  %信噪比范围，从-20dB到20dB
monte=100;     %Monte-Carlo模拟的次数
kk=0;
for nsnr=1:length(SNR)
    Ns_est=0;  %估计的信号源个数
    kk=kk+1;
    correct1=0;
    correct2=0;
    for mk=1:monte
        Renoise=randn(Ne,M);  %生成零均值，方差为1的高斯白噪声
        for n=1:Ne
            Renoise(n,:)=Renoise(n,:)./std(Renoise(n,:));  %标准化噪声
            Renoise(n,:)=Renoise(n,:)-mean(Renoise(n,:));  %去中心化噪声
        end
        %================================================
        
        Ps=10.^(SNR(nsnr)./10);     %信号功率，根据信噪比计算
        S=sqrt(Ps).*SS;     %生成信号，调整信号功率
        Y=A*S+Renoise;    %接收到的信号，是信号和噪声的和
       R=Y*Y'./M;   %计算协方差矩阵

这段代码执行蒙特卡洛模拟，用于评估在不同信噪比下信号检测的性能。它首先生成高斯白噪声，然后调整信号功率，将信号与噪声相加，最后计算接收到的信号的协方差矩阵。