isar_code_new ISAR点目标RD成像算法

ISAR点目标RD成像算法,距离向包络对齐，多普勒中心相位补偿

% domain scatter compensation
% 单特显点法相位校正
% gj 2007 6 1
clear
%% phase compasation
load R_evp.mat;

[NM]=size（R_abs）;
R_sum=sum（abs（R_abs）1）/N;%每一列均值
R_sqr=sum（R_abs.*conj（R_abs）1）/N;
Sigma=ones（1M）;
Sigma=1-（R_sum.*conj（R_sum））./R_sqr;%相关系数

[Yk]=min（Sigma）;
if（Y>0.12）
    error（‘none mainpoint‘）;
end

k
phase=zeros（1N）;
for l=1:N
    phase（l）=R_abs（lk）/abs（R_abs（lk））;
    R_abs（l:）=R_abs（l:）.*conj（phase（l））;
end


plot（R_sum）;
imwrite（mat2gray（angle（R_abs））‘Phasecompastion.jpg‘）;
img=fft（R_abs[]1）;
img = [img（N/2+1:N:）; img（1:N/2:）];  % equivalent to fftshift
imwrite（mat2gray（abs（img））‘image_mainpoint.jpg‘）;

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

     文件      12186  2012-07-30 21:23  isar_code_new\absctract.jpg

     文件        663  2012-07-29 21:56  isar_code_new\dmnsct.asv

     文件        681  2012-07-30 11:09  isar_code_new\dmnsct.m

     文件        691  2012-07-24 20:20  isar_code_new\dpcmp.asv

     文件        913  2012-07-25 09:23  isar_code_new\dpcmp.m

     文件       3100  2012-07-30 21:22  isar_code_new\dstcpmss.asv

     文件       3133  2012-07-30 21:26  isar_code_new\dstcpmss.m

     文件      19148  2012-07-30 21:23  isar_code_new\envolop.jpg

     文件       3082  2012-07-30 21:23  isar_code_new\image_Doppler.jpg

     文件       2658  2012-07-24 20:33  isar_code_new\image_mainpoint.jpg

     文件      56382  2012-07-30 21:23  isar_code_new\Phasecompastion.jpg

     文件     230976  2012-07-30 20:39  isar_code_new\rawdata.jpg

     文件   20171755  2012-07-30 20:39  isar_code_new\rawdata.mat

     文件      30767  2012-07-30 21:23  isar_code_new\R_compress.jpg

     文件    2023539  2012-07-30 21:23  isar_code_new\R_evp.mat

     文件       2378  2012-07-24 18:55  isar_code_new\simdata.asv

     文件       2412  2012-07-30 20:39  isar_code_new\simdata.m

     目录          0  2012-07-24 21:42  isar_code_new

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

             22564464                    18