points.m

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

function [rho,theta]=points(t);
% function POINTS
% Calculates RHO and THETA when called by samples.m

total=t(1);disc=t(8);
freq=t(2);bw=t(3);
vel=t(4);acc=t(5);
bursts=t(6);pul_burst=t(7);
theta0=-.5*(vel*total+.5*acc*total^2);

if disc,
  theta0=-vel/2;
  dtheta=vel/(bursts-1);
  theta=zeros(bursts,pul_burst);
  for a=1:bursts,
    theta(a,:)=theta0+(a-1)*dtheta*ones(1,pul_burst);end
else
  dt=total/(bursts*pul_burst-1);
  t=zeros(bursts,pul_burst);
  for a=1:bursts, for b=1:pul_burst,
    t(a,b)=(a-1)*pul_burst+b;end;end
  t=dt.*t;theta=vel.*t+.5.*acc.*t.^2;
  theta=theta-(max(max(theta))+min(min(theta)))/2;
end

drho=bw/(pul_burst-1);
rho=zeros(bursts,pul_burst);
for a=1:pul_burst,
rho(:,a)=freq+(a-1)*drho*ones(bursts,1);end