fft1.m

高频地波雷达回波信号二维傅里叶变换仿真处理

%构造一个线性调频信号作为差频信号进行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);