%%%%%%參數設置%%%%%%
clc
clear?all
Tp=4e-6;%脈沖寬度
B=300e6;%信號帶寬
k=B/Tp;%調頻斜率
Fs=30e6;%采樣頻率
prf=50;%脈沖周期
c=3.0e8;
lamda=0.05;%波長
w=0.05;%轉速
Rx=[-3?2.5?0?0];%散射點橫坐標
%?Ry=[80?80?78.5?16];%散射點縱坐標
?Ry=[17?13.5?7?-10];%散射點縱坐標
R=sqrt（Rx.^2+Ry.^2）;%散射點到中心點距離
theta0=atan2（RyRx）;%散射點角度
Rmax=10100;%距離門最大距離
Rmin=9900;%距離門最小距離
Rref=10000;%參考距離
pr=c/2/B;%距離分辨率
pa=pr/2;%要求的方位分辨率
delta_theta=lamda/2/pa?%要求的轉角
t_n=delta_theta/w%要求的觀測時間
npulse=ceil（t_n*prf）;%要求的積累脈沖數
nfast=round（（2*Rmax/c-2*Rmin/c+Tp）*Fs）;%距離門采樣數
ndistance=round（Tp*Fs）;%目標距離向采樣數
t=linspace（2*Rmin/c-Tp/22*Rmax/c+Tp/2nfast）;%距離門內時間
tn=linspace（0Tpndistance）;%脈內時間
sr=zeros（npulsenfast）;
ss=zeros（npulsenfast）;
s0=zeros（npulsenfast）;
for?n=1:npulse