oceanmodel.m

海杂波模拟和循环对消法实现杂波抑制的matlab源程序,包含不同信噪比下的目标回波以及最小二乘实现的滤波器完成杂波抑制

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