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