svdrx4.m

用matlab中的SVD函数实现线阵或圆阵的music算法的仿真

function p=RX4()
N=8;
K=100;
p=10;
I=eye(8);
SNR=10^(p/10);
s=sqrt(SNR/2)*randn(3,K)+j*sqrt(SNR/2)*rand(3,K);
Xn=sqrt(1/2)*randn(N,K)+j*sqrt(1/2)*rand(N,K);
n=[1:1:N]';
Vs1=exp(j*pi*sin(pi/9)*cos(pi*13/9-2*pi*(n-1)/N));
Vs2=exp(j*pi*sin(2*pi/9)*cos(pi-2*pi*(n-1)/N));
Vs3=exp(j*pi*sin(pi/3)*cos(pi*4/9-2*pi*(n-1)/N));
As=[Vs1,Vs2,Vs3];
Xs=As*s;
X=Xs+Xn;
fi=2*pi/N*[0:1:N-1]';
[U,S,V]=svd(X);
F=U(:,1:3);
xa=[0:2:360]*pi/180;
ya=[0:1:90]*pi/180;
for k=1:length(ya)
    for m=1:length(xa)
        Va=exp(j*pi*sin(ya(k))*cos(xa(m)-fi));
        p(k,m)=1/(Va'*(I-F*F')*Va);
    end
end
mesh(xa,ya,10*log10(p));