%構造一個線性調頻信號作為差頻信號進行2D FFT處理
clear all;
clc;
%調頻帶寬：B=5MHz;
B=5e6;
%調頻起始頻率f0=5MHz;
f0=5e6;
%調頻周期：T=1ms;
T=1e-3;
%AD采樣頻率：fad=1MHz;
fad=1e6;
t=0:1/fad:T-1/fad;
%目標距離:R0=3500m;
R0=3.5e3;
%目標速度:v=2.6m/s;
v=2.6;
%光速c
c=3e8;
%相參積累周期數N2=256;
N2=256;
n=0:N2-1;
k=B/T;%調頻速率
Rn=R0+v*T.*n;%Rn為一個256維的行向量
fdv=2*f0*v/c; %fdv為速度多普勒頻率
fRn=2*k.*Rn/c; %fRn為距離多普勒頻率，為一256維行向量,相參積累
fdvRn=4*k*v.*Rn/(c^2); % fdvRn為距離速度多普勒頻率，為一256維行向量，相參積累
phase0=4*pi*f0.*Rn/c-4*pi*k.*(Rn.*Rn)/(c^2);%phase0為差頻信號的初始相位，為一256維行向量，相參積累
phase1000=zeros(1000,N2);
for count=1:1000
    phase1000(count,:)=phase0;
end
phasen=2*pi.*(t'*(fdv+fRn+fdvRn))+phase1000;
si=cos(phasen);%差頻信號,1000*256,1000代表采樣點數，256為相參積累數
fftre1=abs(fft(si,1024));% 1024*256
fftre11=fftre1(1:512,:);%512*256
size(fftre11)
x1=(0:255);
y1=1e6/(2*511).*(0:511);
figure(1);
%mesh(x1,y1,fftre11);
fftre2=abs(fft(fftre11',N2));%256*512
fftre22=fftre2(1:N2/2,:);%128*512
[maxcol1,i1]=max(fftre11)
[maxcol2,i2]=max(fftre22')

x2=1e6/(2*511).*(0:511);
y2=1/(T*2*127).*(0:127);
figure(2);
%mesh(x2,y2,fftre22);