阵列信号处理中改进型的盖氏圆信源数估计.zip

阵列信号处理中改进型的盖氏圆信源数估计详细程序，MATLAB代码

clc
clear all
close all
f=3*10^9;
c=3*10^8;
wl=c/f;d=wl/2;
th=[20 40];
N=length（th）;
M=8;
L=100;
rad=pi/180;
SNR=10;
snr=10.^（SNR/20）;
for m=1:M
    for n=1:N
        a（mn）=exp（-j*2*pi*（m-1）*d*sin（th（n）*rad）/wl）;
    end
end
for p=1:N
    for q=1:L
        p1=rand（11）;
        p2=rand（11）;
        si（pq）=sqrt（-2*snr*snr*log（p1））*cos（2*pi*p2）;
        sr（pq）=sqrt（-2*snr*snr*log（p1））*sin（2*pi*p2）;
        s（pq）=si（pq）+j*sr（pq）;
    end
end
RS=s*s‘/L;
for p=1:M
    for q=1:M
        xiangguan=0.3;gonglv=1;
        n（pq）=1*xiangguan^abs（p-q）*exp（j*（p-q）*0.77*pi）;
    end
end
xx=a*RS*a‘+n;
for p1=1:M-1
    for q1=1:M-1
        xx1（p1q1）=xx（p1q1）;
    end
end
[VD]=eig（xx1）;
t=zeros（M-11）;
T=[V t;t‘ 1];
RT=T‘*xx*T;
for i=1:M-1
    r（i）=abs（RT（iM））;%盖氏圆半径
end
for k=1:M-2
    gde（k）=r（k）-sum（r）/（L*（M-1）/80）;
end
kk=1;
while gde（kk）>0
    kk=kk+1;
end
KK1=kk-1;%n=0
%%%%%% 
E=[teye（M-1）;1t‘];
for n=1:M-1
   TT=E*T;
   
RT=TT‘*xx*TT;
for i=1:M-1
    r（i）=abs（RT（iM））;%盖氏圆半径
end
for k=1:M-2
    gde（k）=r（k）-sum（r）/（L*（M-1）/80）;
end
kk=1;
while gde（kk）>0
    kk=kk+1;
end
KK（n）=kk-1;
T=TT;
end
NNN=（sum（KK）+KK1）/M





        
        
        
        
        
        
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1302  2012-11-15 20:11  MGDE.m
     文件          46  2018-10-13 19:14  readme.txt