clear?all;clc;close?all;
%%%%??仿真數據??%%%%
T=0.25;M=100;tc=T/M;fc=1/tc;%%%?一個掃頻周期T?為250ms?；?M為每個周期采樣點100?；?fc采樣頻率
N=256;n=1:N;%%?一個相關積累時間64
%%%
fw=10;%%?雷達工作頻率10MHz
fb=0.102*sqrt（fw）;%%?bragg?頻率
fs=1;%%??信號頻率
a_db=-10;a=10^（a_db/20）;%%?正Bragg峰幅值
b_db=-12;b=10^（b_db/20）;%%?負Bragg峰幅值
c_db=-16;c=10^（c_db/20）;%%?噪聲基底
s_db=-30;s=10^（s_db/20）;%%?信號幅度
bragg=a.*exp（j*2*pi*fb*T.*n）+b.*exp（-j*2*pi*fb*T*n）;%%%?正弦序列產生?Bragg?峰
signal=s.*exp（-j*2*pi*fs*T*n）;
nosic=c.*random（‘Normal‘011N）;
x=bragg+nosic+signal;
%%%
Fx=abs（fftshift（fft（x）））;
freq=（-1/（2*T）:1/（