SAR 三点目标仿真及RD算法 MATLAB程序

clc;
close all;
clear all;

T=10.e-6;                                          %脉宽10us
B=100.e6;                                          %信号带宽100MHz
fs=120.e6;                                         %采样频率
K=B/T;                                             %调频斜率
r1=9850;                                          %目标1距离
r2=10150;                                          %目标2距离
c=3.e8;                                            %光速
tao1=2*r1/c;
tao2=2*r2/c;
N=2*T*fs;                                          %采样点数
t=linspace（-TTN）;
f=（-N/2:N/2-1）.‘*fs/（N-1）;
tran_signal=exp（pi*j*K*t.^2）.*（abs（t）<=T/2）;       %产生线性调频信号
tran_signal1=（abs（t）<=T/2）; 
t0=linspace（-T+tao1T+tao1N）;
rec_signal1=exp（pi*j*K*（t0-tao1）.^2）.*（abs（t0-tao1）<=T/2）;
rec_signal2=exp（pi*j*K*（t0-tao2）.^2）.*（abs（t0-tao2）<=T/2）;
rec_signal=rec_signal1+rec_signal2;                %产生回波信号  

 figure;
 plot（f/1e6fftshift（abs（fft（rec_signal））））;grid on;xlabel（‘MHZ‘）;title（‘频谱‘）;

pulse_compression1=fftshift（ifft（fft（rec_signal1）.*conj（fft（tran_signal））））;
pulse_compression2=fftshift（ifft（fft（rec_signal2）.*conj（fft（tran_signal））））;
pulse_compression=fftshift（ifft（fft（rec_signal）.*conj（fft（tran_signal））））;
pulse_compression3=ifft（fftshift（fft（rec_signal））.*（exp（pi*j*f.^2/K）.*（abs（f）
figure
plot（t0*c/2rec_signal）;
figure;
plot（t0*c/2abs（pulse_compression））;grid on;xlabel（‘r/m‘）;title（‘脉压后‘）;