压缩感知off grid代码论文见2013

压缩感知off grid代码论文见2013 Sparse Frequency diverse MIMO radar imaging for Off-Grid target based on adaptive iterative MAP

function [ target traget_index position_x position_y ] = SACR_iMAP（ signal RefMat Hx Hy x_bound y_bound）
%SAC-iMAP算法
% siganl  为回波信号
% RefMat 为参考矩阵
% Hx 为参考矩阵x方向倒数
% Hy 为参考矩阵y方向倒数

figure（1000）;
[m n] = size（RefMat）;
Mx = sqrt（n）;

target_n = zeros（n1）;
for i = 1:n
    target_n（i） = RefMat（:i）‘*signal/（RefMat（:i）‘*RefMat（:i））;
end

delta_x = zeros（1 n）;
delta_y = zeros（1n）;

target_n1 = target_n;
p = 0.4;

lower_bound = [-ones（1 n）*x_bound -ones（1 n）*y_bound].‘;
upper_bound = - lower_bound;
% delta = [delta_x delta_y];

% ksi为噪声功率的估计值
ksi = norm（signal - RefMat*target_n）/m;

while（1）
    target_n = target_n1;
    fprintf（‘ksi =%f ......................................................................\n‘ksi）;
    alpha = p*ksi/2;
%     alpha = 0.2;
    A = RefMat + [Hx Hy]*[diag（delta_x）;diag（delta_y）];
    %内循环使用了FOCUSS算法
    while（1）
        temp_target_n1 = target_n1;
        diagMat               = diag（abs（target_n1）.^（1-p/2））.^2;
        tran                       = A*diagMat*A‘+alpha*eye（m）;
        transform            = pinv（tran）;
        target_n1             = diagMat*A‘*transform*signal;
        if（norm（target_n1 - temp_target_n1） / norm（temp_target_n1） < 1e-4）
            break;
        end
    end
    
%     imagesc（abs（reshape（