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