加了高斯白噪声后的阵列接收信号
xBsawg
xs
r加了高斯白噪声后的阵列接收信号
v11111111ra
d
Qdiagv噪声协方差矩阵对角线值相同
vra
d
1M
Qdiagv噪声协方差矩阵对角线值不同
Qra
d
M噪声协方差矩阵不为对角阵
RxxQ数据协方差矩阵
以下是MUSIC的程序
UVeigR
UUU1MP
theta900590
谱峰搜索
forii1le
gththeta
AAzeros1le
gthM
forjj0M1
AA1jjexpj2jjpidsi
thetaii180pil
e
d
Pmusiciiabs1AAUUUUAA
e
d
Pmusic10log10PmusicmaxPmusiceps
plotthetaPmusic
xlabelA
glethetadegree
ylabelPthetadB
titleMUSIC非相关源测向仿真fo
tsize13fo
tweightboldfo
t
ame隶书colorred
grido
MUSIC程序结束
附录3ESPRIT算法程序
N1024快拍数
doa102345180pi信号到达角
w100030002000信号频率
M8阵元数
14
fPle
gthw信号个数，也可以用特征分解的大特征值数来决定
l150波长
dl2阵元间距
s
r15信噪比
阵列流形矩阵
BzerosPM
fork1P
Bkexpj2pidsi
doakl0M1
e
d
BB
s10s
r201sqrt2ra
d
3Njra
d
3N仿真信号随机信号
s10s
r20si
w0N1仿真信号正弦信号
s10s
r20expjw0N1仿真信号（指数信号）
s10s
r20sawtoothw0N1仿真信号（锯齿波信号）
s2expjw1N生成信号
xBs
xBs1sqrt2ra
d
MNjra
d
MN
xBsawg
xs
r加了高斯白噪声后的阵列接收信号
v11111111ra
d
Qdiagv噪声协方差矩阵对角线值相同
vra
d
1M
Qdiagv噪声协方差矩阵对角线值不同
Qra
d
M噪声协方差矩阵不为对角阵
RxxQ数据协方差矩阵
以下是ESPRIT程序
RxxR1M11M1M1维的自相关函数
RxyR1M12MM1维的互相关函数
bzeros1M2eyeM2
bbzerosM11
CxxRxxmi
eigRxxeyeM1
CxyRxymi
eigRxxb
aeigCxxCxy
找出最接近1的a值其对应的角度即为φ
a1absabsa1
fori1P
cdmi
a1
a1d1000
bbiad
ad1000
e
d
ifP1
15
fdispThea
glesofsig
alsare
else
dispThea
gleofsig
alis
e
d
Adoa180pi
doaA
DOAasi
a
glebbpipi180
DOAsortDOA
stemdoaDOAr
xlabelDOA实际角度
ylabelDOA估计角度
titleESPRIT非相关源测向仿真fo
tsize13fo
tweightboldfo
t
amecolorred
grido
ESPRIT程序结束
附录4GEESE算法程序
N1024快拍数
doa30030180pi信号到达角
w100030002000信号频率
M8阵元数
Ple
gthw信号个数，也可以用特征分解的大特征值数来决定
Lambda150波长
dLambda2阵元间距
s
r15信噪比
阵列流形矩阵
BzerosPM
fork1P
Bkexpj2pidsi
doakLambda0M1
e
r