tic
%-------將數(shù)據(jù)共軛重排，后進(jìn)行正交運(yùn)算。
clc;
clear;
close?all;
f=4*10^9?;???????%1?.8750e+009;%?信源載波頻率
c=3e+8;
w=c/f;%?信源波長(zhǎng)
d=w*0.4;%?
md=10;%?陣元數(shù)
th=[10?20];%?來波方向角
nd=length?（th）;
snr=0;%:2:30;%?信噪比
snr1=10^（snr/20）;
kp=200;
%?b2c=0.1;
rad=pi/180;

%--------------------------------------------------------------------------
for?m=1:md
????for?n=1:nd
????????a?（mn）=exp?（-j*2*pi*d*（m-1）*sin?（th?（n）*rad）/w）;
????end
end
%--------------------------------------------------------------------------

?for?p=1:nd
????for?q=1:kp
????????p1=rand?（11）;
????????p2=rand?（11）;?%?生成兩個(gè)隨機(jī)數(shù)?（在0-1之間）
????????sr?（pq）=sqrt?（-2*snr1*snr1*log?（p1））*cos?（2*pi*p2）;????%?實(shí)部
????????si?（pq）=sqrt?（-2*snr1*snr1*log?（p1））*sin?（2*pi*p2）;????%?虛部
????????s?（pq）=sr?（pq）+j*si?（pq）;?????????????????????????%?模擬生入射信號(hào)矩陣S?（t）
????end
end
%-------------------------------------------------------------------------
n=wgn?（mdkp0‘complex‘）;%?生成復(fù)信號(hào)高斯白噪聲，md*kp階，功率0dBW即功率為1，幅度為1
%--------------------------------------------------------------------------
x=a*s+n;%?生成模擬入射信號(hào)矩陣
%--------------------------------------------------------------------------
%?X1=x（1:nd1:md）;
%?X2=x（nd+1:end1:md）;
%%%%%%共軛PM
nm=eye（mdmd）;
jiaohuan=rot90（nm）;
ll=conj（x）;
bn=jiaohuan*ll;
z=[xbn];%為md*2kp維的
r=z*z‘/2*kp;%?構(gòu)成協(xié)方差矩陣為md維的

%?K=bn（1:nd1:md）;
%?L=bn（nd+1:end1:md）;
%?rr=size（X1*X1‘+K*K‘）
%?tt=size（X2*X1‘+L*K‘）
%?G=[X1*X1‘+K*K‘;X2*X1‘+L*K‘];
%?H=[X1*X2‘+K*L‘;X2*X2‘+L*L‘];
%?P=inv（G‘*G）*G‘*H;
%?Q=[P‘-diag（ones（1md-nd））];
G=r（:1:nd）;
H=r（:nd+1:end）;
P=inv（G‘*G）*G‘*H;
Q=[P‘-diag（ones（1md-nd））];
%?D=（Q*Q‘）^（-0.5）;
F=Q‘*（（Q*Q‘）^（-0.5））;
theta=-90:0.1:90;
for?i=1:length?（theta）
????AA=exp?（-j*2*pi*d*（0:md-1）‘*sin?（theta?（i）*rad）/w）;
????WW=AA‘*F*F‘*AA;
????pu（i）=abs（1/WW）;

end
pu?=10*log?（pu?）;
theta=linspace?（-90901801）;
figure（1）
plot（thetapu‘k‘‘linewidth‘2）;%?thetapuu‘b‘
grid
xlabel?（‘角度/\[Degree]‘）?
ylabel?（‘空間譜/dB‘）;
title?（‘PM算法均勻線陣譜‘）;

grid?on

%-----------------------------------------
%-----------------------------------------
n=0;%?掃描出譜峰值并顯示出估計(jì)角
for?i=2:1799
????while?（pu?（i）>pu?（i-1））&&（pu?（i）>pu?（i+1））&&（pu?（i）>-16）
????????n=n+1;
????????dd?（n）=pu?（i）;
????????angle?（n）=（i-901）/10??%?將采樣點(diǎn)數(shù)轉(zhuǎn)化為度，一度為10個(gè)步長(zhǎng)
????????break
????end
end

%-----在譜峰輸入搜索數(shù)值----------------------------------------------------
[maxdkd]=max?（dd）;
for?i=1:n
????doa=num2str?（angle?（i））;
????if?i????????text?（angle?（i）dd?（i）doa）;%?在橫坐標(biāo)angle，縱坐標(biāo)dd顯示數(shù)值doa
????else
????????text?（angle?（i）dd?（i）doa）;
????end
end
toc

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2757??2012-06-08?11:03??gong-e-pm\gongepm.asv

?????文件???????2775??2012-06-08?08:28??gong-e-pm\gongepm.m

?????文件???????2893??2012-06-01?20:20??gong-e-pm\gongepmshuju.asv

?????文件???????2893??2012-06-01?20:21??gong-e-pm\gongepmshuju.m

?????目錄??????????0??2014-06-29?09:55??gong-e-pm

-----------?---------??----------?-----??----

????????????????11318????????????????????5