文章目录
前言
雷达算法仿真学习,主要流程包括生成目标回波,将基带信号加在载波上其次是加上目标信息和干扰信号后发射,解调回波信号得到基带信号然后MTI,MTD,最后再检测。
一、基础理解
1.慢时间维度和快时间维度
-
连续M个脉冲回波经过相干解调之后的基带数据形成的一个二维数据矩阵,上图只是雷达三维数据块的一个切面,其中数据的每一行代表的是对同一个距离单元的一连串脉冲测量,每一列对应的是一个脉冲回波的连续采样,也叫连续的距离门,其中每个小方块的元素都是一个复数,代表每个个距离单元的I,Q分量。
-
快时间维或距离维的采样率至少为发射脉冲的带宽,数量级在MHz左右,慢时间维,按雷达的脉冲重复周期进行采样,故采样率为PRF=1/PRT,数量级一般在KHz,所以慢时间维和快时间维其中的慢和快从这里我们也能知道一些原因。
2.脉冲压缩
- 脉冲压缩滤波器实质上就是匹配滤波器,匹配滤波器是以输出最大信噪比为准则设计出来的最佳线性滤波器。假设噪声为高斯白噪声,写出信噪比的表达式,根据schwarz不等式可以得到匹配滤波器的响应,对应的冲激函数形式
- 根据冲击函数形式可以知道匹配滤波之后的输出信号就是输入信号的自相关函数,其傅立叶变换结果就是信号功率谱,则信号带宽越大,输出信号越窄,距离分辨力越好。所以,当宽脉冲的脉内频率或相位经过调制后,信号带宽增大,经过匹配滤波器后就会被压缩为窄脉冲。当然脉冲被压缩的特点也可以从群时延的角度去理解。滤波器的群时延恰好与信号相反,当信号经过滤波器就会呈现多个信号同相位输出,在某个时刻输出压缩成单一载频的窄脉冲。
3.MTI
- 在MTI处理中,沿慢时间维对快时间/慢时间数据矩阵进行高通滤波,得到一个新的数据矩阵,在新的数据矩阵中的杂波分量已经被衰减,此操作可以判断目标的存在,但不能确定目标的运动。
- 典型的MTI滤波器为延迟对消器,进行脉冲间对消。
4. 脉冲多普勒处理
- 直接对慢时间数据序列进行一个谱分析,从而代替滤波处理,处理结果仍然是一个数据矩阵,但是坐标变换为快时间和多普勒频率,此时根据多普勒频率可以求的目标运动的速度。最简单的是对数据矩阵中每一行慢时间序列的离散傅里叶变换(DFT)
5.cacfar检测
- 主要为了的到RD平面的阈值,然后根据阈值判断的出找到多普勒频率,从而得到想要目标的速度。
代码如下:
function [ index, XT ] = ca_cfar( xc, N, pro_N, PAD)
%
alpha=N.*(PAD.^(-1./N)-1);
index=1+N/2+pro_N/2:length(xc)-N/2-pro_N/2;%这句话提前可以提高效率
for j=1:256
XT(j,:)=zeros(1,length(index)+20);
for i=index
cell_left=xc(j,i-N/2-pro_N/2:i-pro_N/2-1);
cell_right=xc(j,i+pro_N/2+1:i+N/2+pro_N/2);
Z=(sum(cell_left)+sum(cell_right))./N;
XT(j,i-N/2-pro_N/2)=Z.*alpha;
end
end
end
二、具体流程
1.产生回波信息加入目标时延以及噪声
A=1;
u=(2*pi*B)/tao;% u的具体含义k*2pi
t=(