算法实现如下：

        parameter：雷达参数设置。

        antVet：目标点 对应的非相参积累天线数据。

function [angle,doa_abs] = caponMethod(parameter,antVec)txAntenna = parameter.txAntenna;  % 发射天线 [1  1]rxAntenna = parameter.rxAntenna;  % 接收天线 [1  1  1  1]txNum = length(txAntenna);        % 发射天线个数  2rxNum = length(rxAntenna);        % 接收天线个数  4lambda = parameter.lambda;        % 波长dx = parameter.dx;   % 天线间隔virtualAntenna = parameter.virtualAntenna;  % 虚拟天线个数Rx = antVec * antVec' / virtualAntenna; deg = -90:0.1:90;a = zeros(virtualAntenna,1);kk = 1;for degscan = degfor txId = 1:txNumfor rxId = 1:rxNumvirtualAntennaId = (txId-1) * rxNum + rxId - 1;dphi = 2 * pi / lambda * dx * virtualAntennaId * sind(degscan);a((txId-1) * rxNum + rxId) = exp(-1i * dphi);endendRxInv = inv(Rx);  % inv矩阵的逆P_out(kk) = 1/(a'*RxInv*a);kk = kk + 1;enddoa_abs = abs(P_out);[pk,loc]=max(doa_abs);angle = deg(loc);end