雷达信号处理技术与系统Word格式.doc
- 文档编号:8444232
- 上传时间:2023-05-11
- 格式:DOC
- 页数:27
- 大小:989.50KB
雷达信号处理技术与系统Word格式.doc
《雷达信号处理技术与系统Word格式.doc》由会员分享,可在线阅读,更多相关《雷达信号处理技术与系统Word格式.doc(27页珍藏版)》请在冰点文库上搜索。
雷达发射机产生线性调频信号,通过天线辐射出去,如果传播过程中遇到目标,就会反射回一部分电磁波,由雷达接收机接收。
这就是回波信号,回波信号中包含有目标的距离,速度,角度等给方面的信息。
目标信号包括期望目标和非期望目标,如图4所示。
图4目标回波产生
由于目标和雷达之间的距离和相对速度的影响,回波信号会产生一定的延迟,以及多普勒频移。
使用传播响应函数来描述目标回波产生的过程,相对发射电磁波传播响应函数如式3-1:
(3-1)
式4-2中,表示目标个数,为噪声信号,和分别为雷达天线发射幅度增益和接收幅度增益,为传播衰减,为各目标的延时时间。
各目标的延时时间满足式3-2:
(3-2)
其中,表示各个目标相对于雷达的速度。
对于雷达天线发射幅度增益和接收幅度增益,采用低旁瓣天线功率方向图进行仿真,使用函数描述天线方向图,即
(3-3)
其中,为天线最大幅度增益,为目标目标偏离雷达发射方向的角度,表示天线3dB带宽。
l传播衰减
回波信号在空气中传播,会发生一定的损耗,称之为传播衰减。
主要包括两个方面,一是大气损耗,二是功率稀释。
大气损耗是雷达工作频率、目标距离和仰角的函数。
雷达工作频率越高,大气损耗越大。
所以,在频率较低的频段(3GHz以下),大气损耗可以不予考虑,在频率较高的频段,进行选择性的考虑。
在本次实验中,雷达工作在中重频下,并没有考虑大气损耗的影响。
功率稀释是由于天线是向所有方向均匀发射能量的,也就是说天线具有球形辐射方向图,所以目标处接收到的电磁波能量的功率密度为,其中为雷达发射功率,为雷达和目标之间的距离,即单程的功率稀释为,类似的,双程的功率稀释为。
最后,传播衰减要同时考虑大气损耗和功率稀释两方面的影响,即。
l调制到中频
雷达是利用物体反射电磁波的特性来发现并确定目标参数的,雷达发射的信号应该是一个载波受到调制的大功率射频信号。
雷达工作频率是按照雷达的用途来确定的,为了调高雷达系统的工作性能和抗干扰能力,有时要求它能在几个频率上跳变工作或者同时工作。
调制到中频,就是对发射信号乘上一个载频信号,即
(3-4)
l加入噪声
在雷达接收机中,除了目标回波信号之外的任何其他信号都成为噪声。
它包括雷达系统之外的干扰信号和雷达接收机内部产生的热噪声。
热噪声(电子的热骚动)和散射效应噪声(半导体的载流子密度的变化)是雷达接收机中两种主要的内部噪声源。
在本仿真实验中,使用了简化模型,即假设加入的噪声是服从高斯分布的,模型如式(3-5):
,(3-5)
其中,,
四、信号处理
雷达信号处理的流程如下:
接收机
正交双通道采样
匹配滤波
MTI
MTD
测距
测速
测角
下面具体介绍如下。
l正交双通道采样
正交双通道处理就是中频回波信号经过两个相似的支路分别处理,其差别仅是其基准的相参电压相位差90°
,这两路称为:
同相支路(InphaseChannel)——I支路
正交支路(QuadratureChannel)——Q支路
传统方法使用的是模拟正交双通道处理,正交I、Q通道处理是将接收机输出的中频回波信号分别与正交的两路相参信号混频(采用模拟乘法器),然后进行低通滤波,从而得到I、Q两路基带信号,再通过A/D变换给出同相分量和正交分量的数字量,如图5所示:
低通滤波
A/D
中频带通信号
图5正交双通道采样结构图
正交双通道处理的优点(相对于单通道处理):
Ø
可区分,以确定目标相对运动方向。
能消除盲相(单通道MTI时目标多普勒信号的相位取样对消导致零输出)。
l匹配滤波
脉冲压缩的目的是集中单个雷达发射信号的所有能量,获的最大输出信噪比。
方法是进行匹配滤波,在接收机中设置一个与发射信号频率相匹配的压缩网络,使经过调制的宽脉冲的回波信号变成窄脉冲,保持良好的距离分辨力。
脉冲压缩网络实际上就是一个匹配滤波器网络。
匹配滤波器是指输出信噪比最大准则下的最佳线性滤波器。
根据匹配理论,匹配滤波器的传输特性:
(4-1)
时域表示(冲激响应)为:
(4-2)
其中,K为幅度归一化常数,S(ω)是发射信号,x(t)是回波信号。
线性调频信号的匹配滤波有两种方法:
时域匹配滤波、频域匹配滤波。
时域匹配滤波:
滑动滤波器(FIR)。
运算量大,难以满足实时处理的要求。
频域匹配滤波:
傅立叶变换后频谱相乘。
可采用FFT算法大幅度降低运算量,满足实时处理的要求。
傅立叶变换后频谱相乘,具有简单的频域解析表达式。
通过加窗能够获得很低的旁瓣,如图6所示。
,
图6频域匹配滤波算法
l动目标显示MTI
动目标显示(MTI)即MovingTargetIndication,是利用MTI滤波器滤除相应杂波,从而提高目标检测性能。
固定目标频谱的谱线位于脉冲重复频率的整数倍点处,而运动目标回波信号存在多普勒频移,动目标显示滤波器利用运动目标回波和杂波在频谱上的区别,有效地抑制杂波而提取信号。
最直接的方法是将相邻重复周期的回波信号相减,则固定目标回波由于振幅不变而互相抵消,运动目标回波相减后剩下相邻重复周期振幅变化的部分。
实验中用到的就是这种传统的非递归型一次对消器,即二脉冲对消。
结构如图7:
图7二脉冲对消结构图
时域方程为:
,传输函数为:
,它是一个单零点系统,零点的位置在,频率响应为:
频率响应如图8,在脉冲重复频率的整数倍点处有凹口,所以固定目标回波在通过MTI滤波器后将受到很大的抑制,理想状态下,输出为零。
图8MTI滤波器频率响应
l动目标检测MTD
MTD也就是一种相参积累和多普勒滤波的结合,相干积累的目的为:
1、集中多个脉冲重复周期/调频周期内雷达发射的所有信号所有能量,获取最大输出信噪比。
2、减小目标RCS起伏对目标检测的影响。
动目标检测(MTD)即MovingTargetDetection,根据最佳线性滤波理论,在杂波背景下检测运动目标回波,除了杂波抑制滤波器外,还应串接有对脉冲串信号匹配的滤波器。
MTD利用了回波脉冲串的相参性进行相参积累。
实际工作中,采用一组相邻且部分重叠的滤波器组覆盖整个多普勒频率范围,这就是窄带多普勒滤波器组。
N个相邻的多普勒滤波器组的实现是由N个输出的横向滤波器(N个脉冲和N-1根迟延线)经过各脉冲不同的加权并求和后形成的。
结构如图9:
图9MTD滤波器组成框图图10MTD滤波器频率响应N=8
设加在第k个滤波器的第i个输出端头的加权值为:
k表示标号从0到N-1的滤波器,每一个k值对应一组不同的加权值,相应地对应一个不同的多普勒滤波器响应。
图10中所示滤波器响应是N=8时加权所得各标记k的滤波器频率响应,k取0~7。
该滤波器的频率覆盖范围为0到fr。
在仿真实验中,通常是通过快速傅里叶变换FFT来实现的。
l恒虚警检测CFAR
在对雷达回波信号作了脉冲压缩、MTI、MTD滤波后,接着就要对目标的存在进行判决:
过门限检测。
这里采用的恒虚警概率检测,即CFAR,将检测门限计算成使雷达接收机能保持恒定的预定虚警率。
检测的原则如下:
(4-3)
图11CFAR检测原理图
恒虚警检测的重点就是确定恒虚警检测的门限。
式4-4给出了门限值VT和虚警概率Pfa之间的关系:
(4-4)
其中,为噪声的功率,由于噪声的功率是一直变化的,为了保持恒定的虚警概率,必须依据噪声方差的估计连续更新门限值。
连续改变门限值以保持恒定虚警概率的过程叫做恒虚警概率(CFAR)。
在仿真实验中,假设,要计算出门限值,就要先估算出噪声的功率。
由恒定,可以推导出:
但是噪声的功率是不知道的,必须要对它进行估计。
因为目标信号的幅度值是比较大的,为了消除目标信号对噪声估计的影响,在被检测目标左右,都有保护单元,不参与对噪声功率的估计,原理图如图12所示,估计的噪声功率如式4-5:
(4-5)
图12噪声估计
l测距
雷达工作时,发射机经天线向空间发射一串重复周期一定的高频脉冲。
如果在电磁波传播的过程上有目标的存在,那么雷达就可以接受到由目标反射回来的回波。
由于回波信号往返于雷达和目标之间,它将滞后于发射脉冲一段时间,如图13所示,电磁波的能量是以光速传播的,设目标的距离为,则传播的距离为光速乘上时间间隔,即
(4-6)
式中,为目标到雷达站的单程距离,单位为m(米),的单位为s(秒),因子是考虑到往返的时间延迟。
图13雷达测距原理
与往返延迟时间(脉冲周期)对应的距离称为雷达的非模糊距离。
因为接收到的回波可以看成使离它最近的发射脉冲的回波,也可以是前一个发射脉冲的回波,在这种情况下,
或
可以看到,该回波就有了距离模糊。
所以,最大无模糊距离必须对应于脉冲周期的一半,即
(4-7)
其中,为脉冲周期,为脉冲重复频率。
l测速
有些雷达除了要确定目标的位置外,也需要确定运动目标的相对速度。
目标运动的速度可以从测量确定时间间隔的距离变化量来定,即。
这种办法称为目标距离微分法测速,这种测速方法需要较长的时间,且不能测定其瞬时速度。
一般来说,测量的准确度也差,其数据只能作为粗侧用。
测速常使用的方法是多普勒测频法。
我们知道,当目标与雷达站之间存在相对速度时,接收到回波信号的载频相对于发射信号的载频会产生一个频移,这个频移在物理学上称为多普勒频移,它的数值为
(4-8)
式中,为多普勒频移,单位为Hz;
为雷达和目标之间的径向速度,单位为m/s;
为载波波长,单位为m。
当目标向着雷达站运动时,,回波载频提高;
反之,回波载频降低,如图14。
雷达只需测量出目标回波的多普勒频移,就可以求得目标相对于雷达的相对速度。
(a)靠近目标(b)后退目标
图14显示多普勒频移的雷达接收信号频谱
多普勒测频法测量速度可靠性比较高,但是容易出现多普勒模糊,脉冲串线谱的包络为形,谱线的间隔为脉冲重复频率,如图15所示。
只要预期的多普勒频移小于各个滤波器带宽的二分之一(即一个FFT门的宽度的二分之一),则多普勒滤波器组就能够解出目标的多普勒频移。
所以,在无模糊的情况下,目标最大多普勒频移为,即脉冲雷达的最大无模糊预期目标的相对径向速度为:
(4-15)
(a)多普勒得到了分辨(b)频谱移到了下一个多普勒滤波器中,发生模糊
图15发射和接收波形的频谱和多普勒组
l振幅和差式测角
雷达测角的物理基础是电波在均匀介质中传播的直线性和雷达天线的方向性。
测角的方法可以分为振幅法和相位法两大类。
在本次实验中,用到的是属于振幅法的比幅单脉冲测角法,这里只对这种测角方法进行简单介绍。
比幅单脉冲测角法采用两个相同且彼此部分重叠的波束,其方向图如图16所示。
如果目标处在两波束的交叠轴OA方向,则由两波束收到的信号强度相等,否则一个波束收到的信号强度高于另一个。
所以常常称OA为等信号轴,这种测角方法也叫做等信号测角法。
当两个波束收到的回波信号相等时,等信号轴所指的方向即为目标方向。
如果目标处在OB方向,波束2收到的回波比波束1的强,处在OC方向时,波束2的回波较波束1的弱,因此,比较两个波束的强弱就可以判断目标偏离等信号轴的方向,并可用查表的办法估计出偏离等信号轴角度的大小。
图16比幅单脉冲测角原理
波束1接收到的回波信号,波束2接收到的回波信号的电压值为,其中,为天线电压方向性函数,为目标方向偏离等信号轴的角度,为与波束最大值方向的偏角。
由和可以求得其差值以及和值,在等信号轴附近可以进行近似,即
归一化的和差值为:
因为正比于目标偏离的角度,故可以用它来判断的方向和大小。
五、仿真结果
l参数
从PD雷达系统构成图中,可以看到该雷达接收到回波信号分为和、差信号两路,雷达发射信号是线性调频脉冲信号,和波束和差波束进行信号处理的方式是完全相同的,在进行CFAR后,进行振幅和差式测角,使用和波束处理后的结果进行测距和测速。
实验中,PD雷达系统及目标的仿真参数设置如下:
表5-1雷达系统仿真参数表
雷达系统仿真参数
雷达发射功率Pt
50W
中心频率Fc
1GHz
脉冲宽度Tp
100us
脉冲重复频率Fr
2kHz
带宽B
5MHz
采样率Fs
10MHz
天线孔径D
5m
天线3dB波束宽度
0.88*λ/D
A、B波束与天线轴向的夹角
0.8*B/2
波束主瓣初始指向
20度
CPI周期数
64
目标参数
目标距离
50Km
目标平均后向散射截面积
1m2
目标方位角
20度
目标速度
100m/s
l结果
发射信号(LFM):
回波信号:
匹配滤波:
MTI:
MTD:
CFAR:
测距:
49.9890m
测速:
100.3125m/s
测角:
20.1815度
六、心得
通过对雷达系统仿真与性能评估这门课的学习,对雷达系统有了一定深度的了解,简单学习了雷达系统各个部分的原理,包括雷达发射机,雷达接收机,天线,目标回波,信号处理机以及数据处理机部分。
对脉压,MTI,MTD,CFAR,测距,测速,测角的原理都有了一定的理解。
在此基础上,学习了雷达系统仿真的方法,包括功能级仿真,信号级仿真,分布式交互仿真和半实物仿真。
课程中,主要是使用Matlab编程实现雷达系统各部分的功能。
七、代码
lLFM信号
c=3e8;
%光速
Fc=1e9;
%中心频率【Hz】
Tp=50e-6;
%脉冲宽度【微秒】
Fr=2000;
%脉冲重复频率【Hz】
B=1e6;
%带宽【Hz】
Fs=10e6;
%采样率【Hz】
K=B/Tp;
%调频率【Hz】
Tr=1/Fr;
%脉冲重复周期【秒】
CPI=10*Tr;
%仿真持续时间【秒】
Delta_t=1/Fs;
%时域采样点时间间隔【秒】
s_lfm_I=zeros(length(round(CPI*Fs)),1);
%定义信号实部数组
s_lfm_Q=zeros(length(round(CPI*Fs)),1);
%定义信号虚部数组
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%数据流仿真方式%%%%%
t_set_Tp=0:
Delta_t:
Tp;
%一个脉冲内的时间采样点数组
s_lfm_Tp_I=cos(pi*K*(t_set_Tp-Tp/2).^2);
%脉冲信号实部
s_lfm_Tp_Q=sin(pi*K*(t_set_Tp-Tp/2).^2);
%脉冲信号虚部
N=length(s_lfm_Tp_I);
forn=1:
10
pulse_start=round((n-1)*Tr*Fs)+1;
%第n个脉冲起始采样点
pulse_stop=round((n-1)*Tr*Fs)+N;
%第n个脉冲结束采样点
s_lfm_I(pulse_start:
pulse_stop)=s_lfm_Tp_I;
%将第n个脉冲信号的实部添加到其在信号数组中的对应位置;
s_lfm_Q(pulse_start:
pulse_stop)=s_lfm_Tp_Q;
%将第n个脉冲信号的虚部添加到其在信号数组中的对应位置;
end
figure
plot(0:
Delta_t*(length(s_lfm_I)-1),s_lfm_I)
axis([0CPI-1.251.25])
title('
Realpart'
)
xlabel('
time[sec]'
ylabel('
magnitude'
l信号处理
clear
%%%%%雷达系统仿真参数%%%%%
%光速
k=1.38e-23;
%玻尔兹曼常数
Pt=50;
%发射功率【W】
%中心频率【Hz】
Wavelength=c/Fc;
%工作波长【m】
Tp=100e-6;
%脉冲宽度【微秒】
%脉冲重复频率【Hz】
B=5e6;
%带宽【Hz】
%采样率【Hz】
F=10^(6.99/10);
%噪声系数
%调频率【Hz】
%脉冲重复周期【秒】
%时域采样点时间间隔【秒】
D=5;
%天线孔径【m】
Ae=1*pi*(D/2)^2;
%天线有效面积【m^2】
G=4*pi*Ae/Wavelength^2;
%天线增益
BeamWidth=0.88*Wavelength/D;
%天线3dB波束宽度【deg】
BeamShift=0.8*BeamWidth/2;
%A、B波束与天线轴向的夹角【deg】
Theta0=20*pi/180;
%波束主瓣初始指向【度】
Wa=0;
2*pi/1;
%天线波束转速【rad/sec】
Num_Tr_CPI=64+1;
%CPI周期数
%%%%%目标仿真参数%%%%%
R_set=[50e3];
%目标距离【m】
RCS=[1];
%目标平均后向散射截面积【m^2】
Theta_target_set=[20]*pi/180;
%目标方位角【deg】
V_set=[100];
%目标速度【m/s】
RCS_Ground_0=10^(-30/10);
%地面单位面积后向散射截面积【m^2】
%%%%%%%%%%%%%%%%%%%%%%%%%%
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 雷达 信号 处理 技术 系统