[ÇóÖú] ¹ØÓÚÈý´ÎÏîÏàλÐźŵÄÒ»´Î²î·ÖÔËËãµÄÎÊÌâÇë´ó¼Ò°ïæ¿´¿´ ÒÑÓÐ1È˲ÎÓë

clear;
clc;
Tr=80e-6;   %Âö³åÖظ´Ê±¼ä80us
T=0.5e-6;     %×ÓÂö³å¿í¶È0.5us
fp=2e6;     %ƵÂʲ½½øÁ¿2MHz
f0=10e9;    %·¢ÉäÔØƵ10G
M=1024;      %µäÐÍ»ýÀÛÂö³åÊý128,¾àÀë·Ö±æÂÊ0.236m
c=3e8;      %¹âËÙ
v=1500;     %ËÙ¶È
a=500;     %¼ÓËÙ¶È

R1=10003;

t=0:Tr:M*Tr-Tr;
y311=exp(-j*4*pi/c*((f0+fp/Tr*t).*(R1+v*(t+T/2+2*R1/c)+1/2*a*(t+T/2+2*R1/c).*(t+T/2+2*R1/c))));
y312=exp(-j*4*pi/c*((f0+fp/Tr*(t+Tr)).*(R1+v*((t+Tr)+T/2+2*R1/c)+1/2*a*((t+Tr)+T/2+2*R1/c).*((t+Tr)+T/2+2*R1/c))));
r1=y311.*conj(y312);
z1=fftshift(fft(r1));