clear;close?all;clc;
%=========參數(shù)初始化=========%
num=8;%陣元數(shù)
f=10e9;%信號(hào)載波頻率
fs1=20;%信號(hào)頻率
finF1=10;%干擾1頻率
finF2=70;%干擾2頻率
a1=90;%信號(hào)方向
inF1=30;%干擾1方向
inF2=150;%干擾2方向
c=3e8;%光速
lmda=c/f;%波長
d=lmda/2;%陣元間距
thita=0:1:180;
fs=500;
Ts=[0:1/fs:1];
snapshot=length（Ts）;
SNR=10;%信噪比
INR1=30;%干噪比
INR2=0;%干噪比
s=10^（SNR/20）*sin（2*pi*fs1*Ts）;%信號(hào)包絡(luò)
s_inF1=10^（INR1/20）*sin（2*pi*finF1*Ts）;%干擾1包絡(luò)
s_inF2=10^（INR2/20）*sin（2*pi*finF2*Ts）;%干擾1包絡(luò)
athita=exp（j*2*pi*d/lmda*cos（a1*pi/180）*[0:num-1]）.‘;%信號(hào)導(dǎo)向矢量
athita_inF1=exp（j*2*pi*d/lmda*cos（inF1*pi/180）*[0:num-1]）.‘;%干擾1導(dǎo)向矢量
athita_inF2=exp（j*2*pi*d/lmda*cos（inF2*pi/180）*[0:num-1]）.‘;%干擾2導(dǎo)向矢量
x=athita*s+athita_inF1*s_inF1+athita_inF2*s_inF2+sqrt（0.5）*（randn（numsnapshot）+j*randn（numsnapshot））;%輸入
%===========Capon=============%
R=x*x‘/snapshot;
Wopt=inv（R）*athita/（athita‘*inv（R）*athita）;
for?i=1:length（thita）
????y（i）=Wopt‘*exp（j*2*pi*d/lmda*cos（thita（i）*pi/180）.*[0:num-1]）.‘;
end
figure（1）
plot（thita20*log10（abs（y）/max（abs（y））））;grid;hold?on;
plot（a1-50:0‘.-‘）
plot（inF1-50:0‘.-‘）
plot（inF2-50:0‘.-‘）
title（‘capon方向圖‘）
xlabel（‘方位角（度）‘）
ylabel（‘輸出功率（dB）‘）
%===========Capon_DOA============%
for?i=1:length（thita）
????A=exp（j*2*pi*d/lmda*cos（thita（i）*pi/180）.*[0:num-1]）.‘;
????Scapon（i）=1/（A‘*inv（R）*A）;
end
figure（2）
plot（thita20*log10（abs（Scapon）））;grid;hold?on;
title（‘capon測角‘）
xlabel（‘方位角（度）‘）
ylabel（‘輸出功率（dB）‘）

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

?????目錄??????????0??2007-08-12?22:12??DOA

?????文件???????1063??2007-06-03?22:58??DOA\face_direction_pic.m

?????文件???????1531??2007-06-06?22:10??DOA\capon_example.m

?????文件???????1660??2007-06-03?20:44??DOA\ESPRINT.m

?????文件???????1197??2007-06-05?22:39??DOA\MN_music.m

?????文件???????1027??2007-06-04?21:23??DOA\music.m

?????文件????????606??2007-05-29?18:49??DOA\ULA1.m

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

?????????????????7302????????????????????8