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）?????????????break;
????????end
????end
????
%?????imagesc（abs（reshape（