close?all
clear?all
clc
c=3.0e8;
RF=77e9;
B=350e6;
Tchirp=38.9e-6;%調制周期
N=256;%采樣點數
k=128;%積累周期
S=B/Tchirp;%調制斜率
la=c/RF;%載波波長
BL=la/2;%天線間距
TS=Tchirp/N;%采樣間隔
FS=1/TS;%采樣頻率
%%%%目標參數%%%
targetn=4;
R=[5?20?30?50];
V=[0?5??10?15];
jiao=[pi/4?pi/6?pi/9?pi/18];
r=BL*sin（jiao）;%波程差
xc=2*pi*BL*sin（jiao）/la;%相位差

%%%%%航跡追蹤顯示%%%%%%%%%%
Mx=zeros（2564）;
My=zeros（2564）;
for?i=1:256;
???for?m=1:4;
d=Tchirp*128*V（m）*i+R（m）;
x=d.*sin（jiao）;
y=d.*cos（jiao）;
Mx（im）=x（:m）;
My（im）=y（:m）;
k=（i-1）*4+m;
Wx（k）=Mx（i