三角调频连续波雷达 matlab仿真

clc;
close all;
fa=5e6;
fb=10e6;%上下频率
T=0.001;   %调制时间
B=fb-fa;  %带宽
fs=4*B; %采样频率
Ts=1/fs;%采样间隔
n=8;%周期数
num=round（T/Ts）;%周期采样点数
k=2*B/T;  %斜率
tao=0.01*T;%延迟时间
num_late=round（tao/Ts）;%延迟点数

 t=（0:num/2-1）/fs;%上变频时间序列
 t1=（num/2-1:-1:0）/fs;%下变频时间
 
 
%%%%%%%%%%%%%%%%%%%%%%%上下变频信号
freq_up=fa+k*t;
freq_down=fb-k*t;

f_up=cos（2*pi*freq_up.*t）;
f_down=cos（2*pi*freq_down.*t1）;

% figure（）;
% subplot 411;
% plot（freq_up）;
% subplot 412;
% plot（freq_down）;
% subplot 413;
% f_up=f_up（1:2^12）;
% y1=fftshift（abs（fft（f_up）））;
% y1=y1（1:length（y1）/2）;
% plot（y1）;
% subplot 414;
% f_down=f_down（1:2^12）;
% y2=fftshift（abs（fft（f_down）））;
% y2=y2（1:length（y2）/2）;
% plot（y2）;


%%%%%%%%%%%%%%%%%%
ft_period=[f_upf_down];
for i=1:n-1