final.m

进行ISAR运动补偿的程序 kylsdone 方法来消除越距离单元走动,注释完整

% FINAL.M interpolation program
% Modify it to select the interpolation function
% to be used.  Saves results in SDATA.MAT
% Loads gdata(conttaining RHO and THETA) and zdata (containing Z)
% Example using threest.m

clear all; close all; clc; tic;
load gdata; ; load zdata;

%%% Interpolation in Frequency domain
disp('Interpolation')
[s,V] = threest(z,rho,theta); % Select the interpolation method in this line
%[s,V]=euclide(z,rho,theta);
%[s,V]=nene(z,rho,theta);
%[s,V]=fsoum1(z,rho,theta);
%[s,V]=fsoum2(z,rho,theta);

%%% Hamming window
disp('Hamming window')
wz1 = hamming(size(z,1)); wz2=hamming(size(z,2)); wz = wz1 * wz2.';
zk=z; sk=s;
z = z .* wz;
s = s .* wz;
clear wz

%%% Transformation of data to space domain
disp('Transformation of data')
hs = abs(fftshift(fft2(flipud(z),size(z,1),size(z,2))));
hs = hs ./ max(max(hs)); 
hc = abs(fftshift(fft2(flipud(s),size(z,1),size(z,2))));
hc = hc ./ max(max(hc)); clear z s
hsk = abs(fftshift(fft2(flipud(zk),size(zk,1),size(zk,2))));
hck = abs(fftshift(fft2(flipud(sk),size(zk,1),size(zk,2))));

noisefloor=0.01;    % -40db. Change as desired.
for a = 1:size(hs,1), for b = 1:size(hs,2),
if hs(a,b) < noisefloor, hs(a,b) = noisefloor; end;end;end
for a = 1:size(hc,1), for b = 1:size(hc,2),
if hc(a,b) < noisefloor, hc(a,b) = noisefloor; end
end; end

X0=V(1); dkx=V(2); dky=V(3); fl=V(4); time=toc

save sdata hc hs noisefloor dkx dky X0 time fl
soplot
%ents=entropy(hsk)*100;
%entc=entropy(hck)*100;