clear?all;
close?all;
clc;
%%設置信號為窄帶信號，幅度為1
S=[111]
imag1=sqrt（-1）;
element_num=8;%陣元數為8
d_lamda=1/2;%陣元間距d與波長lamda的關系
theta=linspace（-pi/2pi/22000）;
theta0=0/180*pi;%來波方向?（我覺得應該是天線陣的指向）
%theta0=0;%來波方向
theta1=12/180*pi;
theta2=30/180*pi;
L=20;
kk_sum=0;
SNR=8;
w_l0=exp（-imag1*2*pi*d_lamda*sin（theta0）*[0:element_num-1]‘）;
w_l1=exp（-imag1*2*pi*d_lamda*sin（theta1）*[0:element_num-1]‘）;
w_l2=exp（-imag1*2*pi*d_lamda*sin（theta2）*[0:element_num-1]‘）;
A=[w_l0w_l1w_l2];
S=randn（length（S）L）
x=A*S;
X=awgn（xSNR‘measured‘）;
R=X*X‘/L;


[VD]=eig（R）;
Ci=inv（R）;?

for??j=1:length（theta）??%（我認為是入射角度，即來波方向計算陣列流形矩陣A）
?????a=exp（-imag1*2*pi*d_lamda*sin（theta（j））*[0:element_num-1].‘）;
?????
?????p（j）=1/（a‘*Ci*a）;??%（matlab中的‘默認為共軛轉置，如果要計算轉置為w.‘*a）
end

figure（1）
plot（theta*180/pidb（abs（p）/max（abs（p））））;
xlabel（‘角度/°‘）;
ylabel（‘歸一化增益/dB‘）;
title（‘0°，12°，30°時capon算法天線方向圖‘）
figure（2）?????????????????????????????%繪出應用MUSIC法估計的波達方向圖????????????
polar（thetaabs（p））;
%axis（[-90?90?-60?0]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1152??2017-12-15?20:43??capon.m
?????文件????????1703??2017-12-15?21:19??capon_L.m
?????文件????????1420??2017-12-15?22:21??capon_MUSIC.m
?????文件????????1634??2017-12-15?21:05??capon_SNR.m