CA-CFAR代码实现

置顶

qq_36722544

于 2018-08-21 19:50:28 发布

阅读量1.2w

点赞数 10

文章标签： CFAR 电子恒虚警平均算法

本文链接：https://blog.csdn.net/qq_36722544/article/details/81913411

版权

本文档介绍了雷达信号处理中的恒虚警（CFAR）算法，特别是连续平均（CA）CFAR的MATLAB实现。首先，通过设定参数创建了包含Bragg峰、信号和噪声的仿真数据。接着，对数据进行傅里叶变换，并应用CA-CFAR算法进行处理。最后，对处理结果进行归一化幅值的绘制，展示了CA-CFAR算法的效果。

摘要由CSDN通过智能技术生成

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=-35;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’,0,1,1,N);
x=bragg+nosic+signal;
%%%
Fx=abs(fftshift(fft(x)));
freq=(-1/(2*T):1/(T*N):(1/(2*T)-1/(T*N)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CA-CFAR %%%%%%%%%%%%%%%%%%%%%%%%%
M=N;P_fa=10.^(-6);
R=16;%参考单元数
n=R/2;
L_slipper=R+1;%