fsoum2.m

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

function [s,V] = fsoum1(z,rho,theta)
% Soumekh Interpolation Method in Frequency
%
% Usage:  [s,V] = fsoum1(z,rho,theta)
% Where:  z = M by N matrix with samples in frequency domain
%         rho, theta = radius and angles of data points
%         s = interpolated matrix
%         V = vector that includes [X0 dkx dky fl]

close all

% Variables
bursts = size(z,1);
pul_burst = size(z,2);
rho = rho .* (4 * pi / 3e8);	%% frequency matrix in rad/mt
rhov = rho(1,:);
ac = max(max(theta)) - min(min(theta));
ca = (max(max(theta)) + min(min(theta))) / 2;

% Points of Interpolation
x = linspace(min(rhov), max(rhov) * cos(ac / 2), bursts);
y = min(rhov) * sin(max(theta(:,1))) * linspace(-1,1,pul_burst);
[yy,xx]=meshgrid(y,x);
fxi=xx.*cos(ca)-yy.*sin(ca);
fyi=yy.*cos(ca)+xx.*sin(ca);
clear xx yy
fxd=rho.*cos(theta);
fyd=rho.*sin(theta);
dkx=fxi(1,2)-fxi(1,1);
dky=fyi(2,1)-fyi(1,1);
X0=pi/(rhov(2)-rhov(1));
figure;plot(fxd,fyd,'.',fxi,fyi,'.');pause;close

% Interpolation
flops(0)
s=zeros(size(z));
al=rhov(2)-rhov(1);
be=theta(1,2)-theta(1,1);
for a=1:bursts,
for b=1:pul_burst,
a,b
arg=X0.*sqrt((fxi(a,b)-fxd).^2+(fyi(a,b)-fyd).^2);
I=z.*(X0^2).*sqrt(fxd.^2+fyd.^2).*bessel(1,arg)./arg;
ss=sum(sum(I));
s(a,b)=ss;
end;end
fl=flops;

V=[X0 dkx dky fl];