clc;
clear all;
g=9.8;                       %重力加速度
c=3e8;                       %光速
f0=6.75e6;                   %雷達工作頻率
f_Bragg=sqrt(g*f0/(pi*c));   %Bragg頻率
fb1=f_Bragg-0.02;            %第一個窄帶信號的中心頻率
fb2=-f_Bragg-0.01;           %第二個窄帶信號的中心頻率，兩者與＋、－Bragg頻率基本吻合
fs=0.5;                      %模擬目標在多普勒譜上所在頻率
B1=0.625;                    %調頻參數
B2=0.775;                    %調頻參數
f_m1=0.08;                   %調頻參數
f_m2=0.065;                  %調頻參數
sita1=2*pi*rand(1,1);        %隨機初相
sita2=2*pi*rand(1,1);        %隨機初相
T=0.5;                       %采樣間隔（掃頻周期）
Num_T=256;                   %采樣點數（掃頻周期數）
t=T*(0:Num_T-1);
doppler=1/T/Num_T*(-Num_T/2:Num_T/2-1);
dop_Bragg=doppler/f_Bragg;
v_doppler=c*doppler/(2*f0);
s1=(1.1-0.2*rand(1,1))*exp(j*(2*pi*fb1*t+B1*cos(2*pi*f_m1*t+sita1)));   %模擬海雜波
s2=(1.1-0.2*rand(1,1))*exp(j*(2*pi*fb2*t+sita2+B2*cos(2*pi*f_m1*t)));   %模擬海雜波
s3=exp(j*2*pi*fs*t);          %模擬目標
SCR=-8;                       %信雜比
a_s=10^(SCR/20);              %目標信號相對海雜波的幅度
CNR=20;                       %雜噪比
a_c=10^(CNR/20);              %海雜波相對噪聲的幅度
s=a_c*(a_s*s3+s1+s2)+randn(1,Num_T);        %模擬的窄帶信號
FT_s=fftshift(fft(s));
figure(1)
plot(doppler,20*log10(abs(FT_s)/max(abs(FT_s))));
xlabel('f (Hz)');
ylabel('Am (dB)');
title('before accumulation');
grid
Num_C=1;                       %循環對消次數
[FT_s_ac,s_ac]=Func_OCC_root(s.',T,Num_C);     %循環對消
figure(2)
plot(doppler,20*log10(abs(FT_s_ac)/max(abs(FT_s_ac))));
xlabel('f (Hz)');
ylabel('Am (dB)');
title('before accumulation');
grid