这段时间要做下keystone算法 编程后有点问题 不知道出在哪里 希望有人指点一下：
主要问题1.我的仿真数据能不能用keystone算法 2.程序中进行插值时频率f怎么求？
下面是我编的一段程序
fc=2e8;
c=3e8;
v=400;
theat=pi/3;
fs=500;
im1=zeros(100,960);
r=zeros(1,960);
l=zeros(1,100);
for i=1:960
    r(i)=sqrt((100000-v*cos(theat)*(i-1)/fs)^2+(v*sin(theat)*(i-1)/fs)^2);
    l(i)=round((r(i)-99000)/20);
    im1(l(i),i)=exp(j*4*pi*fc*r(i)/c);
end
x=fft(im);
y=zeros(100,960);
for i=1:100
    for m=1:960
        n=1:960;
        y(i,m)=sum(x(i,:).*sinc(（fc/fc+f）*m-n));
    end
end
F=ifft(y);
figure;
imagesc(abs(F))
求指点[em2][em2][em2][em2]