%capon和CBF波束形成
clear?all;close?all;clc
M=16;?????????%均勻線陣列數目
P=2;??????????%信號源數目
f0=2500;f1=1000;%信號頻率Hz
nn=2000;??????%快拍數
angle1=-10;angle2=30;%信號入射角度
th=[angle1;angle2]‘;
SN1=-40;SN2=-40;%信噪比
sn=[SN1;SN2];
degrad=pi/180;
tt=0:1/2000:（1-1/2000）;
x0=exp（-j*2*pi*f0*tt）;?%2個信號x0、x1
x1=exp（-j*2*pi*f1*tt）;?
t=1:nn;
S=[x0（t）;x1（t）];
nr=randn（Mnn）;
ni=randn（Mnn）;
u=nr+j*ni;??????%復高斯白噪聲
pn=2?;??????????%復高斯噪聲功率
tmp=sqrt（pn*10.^（sn/10））;%信號功率
S2=diag（tmp）*S+u（[1?2]:）;?%計算信號幅值，滿足所設信噪比
AA1=[exp（-pi*j*[0:M-1]*sin（angle1*pi/180））]‘;%陣列流型;系數的設置對結果影響非常大%
AA2=[exp（-pi*j*[0:M-1]*sin（angle2*pi/180））]‘;%估計是涉及到周期的問題????????????%
A=[AA1AA2];
X=A*S2;??%接收到的信號
Rxx=X*X‘./nn;??%相關矩陣
invRxx=inv（Rxx）;
%信號搜索
th2=[-pi/2:pi/180:pi/2];
for?n=1:length（th2）
????a（n）=th2（n）;
????A2=[exp（-1.1*pi*j*[0:M-1]*sin（a（n）））]‘;%搜索空間的方向向量表示??????????????????????????????
????den=A2‘*invRxx*A2;
????P_capon（n）=abs（1./den）;%capon波束形成
????P_cbf（n）=A2‘*Rxx*A2;%CBF波束形成
end;
semilogy（th2*180/piP_capon‘r‘）;
hold?on;
semilogy（th2*180/piP_cbf‘b‘）;
axis（[-90?90?1e-4?1e5]）;
grid?on;
xlabel（‘角度‘）;ylabel（‘P功率譜‘）;
title（‘SNR1=-40?SNR2=20?th1=-10?th2=30‘）;
legend（‘capon‘‘CBF‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1325??2002-01-27?21:17??capon_MVM.m
?????文件?????????756??2015-05-03?23:10??esprit.m
?????文件????????1010??2015-04-21?20:02??music.m
?????文件????????2096??2015-03-31?21:02??music_capon_practice.m
?????文件?????????919??2015-04-21?22:38??music_pinghua.m