%?Capon?algorithm?for?DOA?estimation
clc
clear?all
close?all
sensor_number=10;?%陣元數(shù)
source_number=2;?%信源數(shù)
T=200;??%陣列采樣快拍數(shù)
sig1=exp（j*（0.1*pi*（0:T-1）））.‘;
sig2=exp（j*（0.2*pi*（0:T-1）））.‘;??%窄帶復正弦信號
S0=[sig1.‘;sig2.‘];??%信號矩陣形式
snr=10;?%信噪比dB
source_doa1=-45*pi/180;
source_doa2=60*pi/180;
source_doa=[source_doa1?source_doa2];?%信號入射角度
A=zeros（sensor_numbersource_number）;
A=[（exp（-j*pi*（0:sensor_number-1）*sin（source_doa（1））））.‘?（exp（-j*pi*（0:sensor_number-1）*sin（source_doa（2））））.‘];??%方向矩陣
%%%高斯白噪聲*********************************************************
n=zeros（sensor_numberT）;?
real_noise0=randn（sensor_numberT）;
imag_noise0=randn（sensor_numberT）;
mean_real_noise=mean（mean（real_noise0））;
mean_imag_noise=mean（mean（imag_noise0））;
real_noise=real_noise0-mean_real_