capon算法的DOA估计MATLAB代码

% Capon algorithm for DOA estimation
clc
clear all
close all
sensor_number=10; %阵元数
source_number=2; %信源数
T=200;  %阵列采样快拍数
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_