基于capon算法的MATLAB代码

基于capon算法的MATLAB代码，针对雷达DOA估计的MATLAB一维capon的MATLAB代码

%%Capon算法DOA估计
%MUSIC算法
%%错误：相位s_phase没有与j相乘，这样就不是改变相位，而是改变幅度了。
%%问题：频率相同，是不是相干信号？
%%仿真现像：同频能估计出所有的入射方向，但是有偏差；不同频，估计出所有的入射方向，偏差较小；
%%U（1:D）为信号子空间，U（D+1:end）为噪声子空间。
%%如果信号子空间维数小于信号源个数，就是有信号子空间扩散到了噪声子空间，这样会导致某些导...
%%向矢量与噪声子空间不完全正交，而无法正确估计信号源方向。
%%如果信号子空间维数大于信号源个数，噪声子空间扩散到了信号子空间，扩散的部分与剩下的噪声...
%%子空间正交，而出现假峰。
%%由上可知，准确获知信号源的个数对正确估计信号的方向起着重要作用。此外，相干信号会引起信...
%%号子空间的扩散

clear all;clc
c=3e8;%光速
s_angle=pi*[-30 10 50]/180;%入射角度
M=8;%阵元数
D=3;%信号源数
s_freq=[3.02e3 3.01e3 3e3];%信号频率为30MHz
s_wavelength=c./s_freq;
s_sample=7e3;%采样频率为70MHzs
d=s_wavelength（2）/2;%阵元间距为半波长
s_phase=[0.1 0 -0.3]*（2*pi）;%初相
s_snapshoot=100;%快拍数
N=normrnd（01Ms_snapshoot）;%不相关噪声，Mxs_snapshoot
SNR=[10 10 10];%信噪比
%Amp=sqrt（10.^（SNR/10））;%各复信号幅度，其中噪声功率为1
Amp=sqrt（2*10.^（SNR/10））;%各实正弦信号幅度，其中噪声功率为1

%入射信号S（Dxs_snapshoot）每一行，为同一信号不同时间的采样
%%（1）
n=[1:s_snapshoot];
for k=1:D
    S（k:）=Amp（k）*sin（2*pi*n*s_freq（k）/s_sample+s_phase（k））;
    %S（k:）=Amp（k）*exp（j*（2*pi*n*s_freq（k）/s_sample+s_phase

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2490  2010-07-12 20:11  capon.m

----------- ---------  ---------- -----  ----

                 2490                    1