ml.m

％极大似然方法 gcn为高斯有色噪声！！！比较不同点

%极大似然方法  gcn为高斯有色噪声！！！比较不同点！！！！！
clear all;
M=5;
I=3;
N=1024;
n=2;
for run=1:n
A=manifoldmusic(M,I);
%w=normrnd(0,0.1,M,N);
for m=1:M
    w(m,:)=gcn(N);  %
end
S=exprnd(0.9,I,N);
sigma2s=cum2x(S,S,0,0,0,'biased')
sigma2w=cum2x(w,w,0,0,0,'biased')
SNR=10*log10(sigma2s/sigma2w)
y=A*S+w;
R=zeros(M,M);
R0=zeros(M,M);
for k1=1:M
    for k2=1:M
          R(k1,k2)=mean(y(k1,:).*conj(y(k2,:)));
    end
end
R0=R0+R;
end
R0=R0/n;
R1=inv(R0);
for m=1:M
  for n=1:180
     A1(m,n)=exp(-1*j*(m-1)*pi*cos(n*pi/180));
   end
end

%P=10 * log(abs(A1'*R1*A1));  %这个为何与下一个有很大差别，MUSIC为何相同！？？？
P=1./abs(A1'*R1*A1);
n=1:180;
figure(1)
plot(n,abs(P(91,1:180)));
zoom on