clear?all;clc;close?all;
%%%%??仿真數(shù)據(jù)??%%%%
T=0.25;M=100;tc=T/M;fc=1/tc;%%%?一個掃頻周期T?為250ms?；?M為每個周期采樣點(diǎn)100?；?fc采樣頻率
N=256;n=1:N;%%?一個相關(guān)積累時間64
%%%
fw=10;%%?雷達(dá)工作頻率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）;%%?負(fù)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）;%%%?正弦序列產(chǎ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/（