脉冲压缩雷达与匹配滤波【定稿材料】
- 格式:doc
- 大小:229.58 KB
- 文档页数:11
雷达系统中的信号处理技术摘要本文介绍了雷达系统及雷达系统信号处理的主要内容,着重介绍与分析了雷达系统信号处理的正交采样、脉冲压缩、MTD和恒虚警检测几种现代雷达技术,雷达系统通过脉冲压缩解决解决雷达作用距离和距离分辨力之间的矛盾,通过MTD来探测动目标,通过恒虚警〔CFAR〕来实现整个系统对目标的检测。
关键词雷达系统正交采样脉冲压缩MTD 恒虚警检测1雷达系统概述雷达是Radar〔Radio Detection And Ranging〕的音译词,意为“无线电检测和测距”,即利用无线电波来检测目标并测定目标的位置,这也是雷达设备在最初阶段的功能。
雷达的任务就是测量目标的距离、方位和仰角,还包括目标的速度,以及从目标回波中获取更多有关目标的信息。
典型的雷达系统如图1,它主要由雷达发射机、天线、雷达接收机、收发转换开关、信号处理机、数据处理机、终端显示等设备组成。
图1雷达系统框图随着现代电子技术的不断发展,特别是数字信号处理技术、超大规模集成数字电路技术、电脑技术和通信技术的告诉发展,现代雷达信号处理技术正在向着算法更先进、更快速、处理容量更大和算法硬件化方向飞速发展,可以对目标回波与各种干扰、噪声的混叠信号进行有效的加工处理,最大程度低剔除无用信号,而且在一定的条件下,保证以最大发现概率发现目标和提取目标的有用信息。
雷达发射机产生符合要求的雷达波形,然后经馈线和收发开关由发射天线辐射出去,遇到目标后,电磁波一部分反射,经接收天线和收发开关由雷达接收机接收,然后对雷达回波信号依次进行信号处理、数据处理,就可以获知目标的相关信息。
雷达信号处理的流程如下:图 2 雷达信号处理流程2雷达信号处理的主要内容雷达信号处理是雷达系统的主要组成部分。
信号处理消除不需要的杂波,通过所需要的目标信号,并提取目标信息。
内容包括雷达信号处理的几个主要部分:正交采样、脉冲压缩、MTD和恒虚警检测。
正交采样是信号处理的第一步,担负着为后续处理提供高质量数据的任务。
雷达脉冲压缩算法研究雷达脉冲压缩算法是一种通过对短脉冲信号进行加窗和相关运算,从而实现高分辨率雷达成像的算法。
这种算法在目标探测、识别以及跟踪等领域中有着广泛的应用。
在本文中,我们将深入地探讨雷达脉冲压缩算法的基本原理、发展历史以及未来的研究方向。
一、基本原理脉冲雷达技术中,发射的信号被目标反射后接收到信号会被传回雷达接收机。
然而,目标信号在传输过程中会遭受多径效应的干扰,这导致接收到的信号在时间域上发生扩展,时间分辨率会降低。
为了解决这个问题,雷达脉冲压缩技术应运而生。
雷达脉冲压缩算法主要基于短脉冲信号的性质,即其具有宽带性和瞬时功率很大。
算法的基本步骤为:先对短脉冲信号进行加窗,使其具有良好的频谱特性;然后进行相关运算,使反射信号会在一段极短的时间内被压缩,从而提高时间分辨率。
加窗操作的目的是消除反射信号的频率偏移,使其具有宽带性。
常用的窗函数有海明窗、布莱克曼窗、汉宁窗等。
这些窗函数在保留谱线的同时,在频域上也可以压缩主瓣宽度。
相关运算的基本原理是将原始信号与一个匹配滤波器进行卷积,从而使信号被在一小段时间内压缩。
匹配滤波器通常是原始信号的逆时域复共轭,其功率频谱密度与信号的功率频谱密度接近,但是带宽更宽。
二、历史发展雷达脉冲压缩算法的诞生最早可以追溯到20世纪50年代初。
当时,人们意识到脉冲雷达系统的时间分辨率受到多径效应的限制,无法满足目标识别和跟踪的需求。
为解决这个问题,一些科学家开始研究如何对反射信号进行压缩,并尝试应用于实际应用中。
在此后的数十年中,雷达脉冲压缩算法经历了一个逐步发展的过程。
20世纪70年代末,复合式高分辨雷达(SAR)系统的出现使得脉冲压缩技术得到了广泛的应用。
90年代初,人们开始对逆问题进行研究,从而进一步提高了脉冲压缩算法的效率和精度。
三、未来研究方向在当今的信息技术快速发展的时代,雷达脉冲压缩算法如何更好地适应未来的发展成为了一个重要的问题。
未来研究方向主要包括以下三个方面:1. 面向多异步输入的实时压缩算法。
从匹配滤波的角度说明脉冲压缩的基本原理
脉冲压缩是雷达信号处理中的一种重要技术,其基本原理是将宽脉冲信号通过匹配滤波器进行压缩,以提高雷达的分辨率和探测能力。
下面从匹配滤波的角度解释脉冲压缩的基本原理。
匹配滤波器是一种特殊的线性滤波器,其输出信号的频谱与输入信号的频谱成共轭对称。
在雷达信号处理中,匹配滤波器被用于接收和处理回波信号,以获得最大的信噪比。
在脉冲压缩雷达中,发射信号通常是一个宽脉冲,其频谱具有较大的带宽。
当这个宽脉冲信号照射到目标后,目标的散射回波信号被接收。
由于目标的距离不同,回波信号的延迟时间也不同。
为了实现高分辨率,需要对回波信号进行匹配滤波处理。
匹配滤波器对输入信号的频谱具有选择性,只有与滤波器频谱相匹配的信号才能通过滤波器。
在脉冲压缩雷达中,匹配滤波器的频谱与发射信号的频谱相匹配,因此,当回波信号通过匹配滤波器时,其频谱会被压缩。
由于匹配滤波器输出的信号与输入信号的频谱成共轭对称,因此输出信号的带宽与输入信号相同。
通过匹配滤波器的压缩处理,回波信号的脉冲宽度被显著压缩,从而提高了雷达的分辨率。
同时,由于匹配滤波器能够使回波信号的信噪比最大化,因此提高了雷达的探测能力。
总之,脉冲压缩的基本原理是通过匹配滤波器的压缩处理,将宽脉冲信号转换为窄脉冲信号,同时保持其带宽不变,从而实现高分辨率和探测能力的提升。
匹配滤波技术在雷达信号处理中的应用一 、匹配滤波原理在输入为确知加白噪声的情况下,所得输出信噪比最大的线性滤波器就是匹配滤波器,设一线性滤波器的输入信号为 :)(t x)()()(t n t s t x +=其中)(t s 为确知信号,)(t n 为均值为零的平稳白噪声,功率谱密度为2/No 。
设线性滤波器系统的冲击响应为)(t h ,其频率响应为)(ωH ,其输出响应:)()()(t n t s t y o o +=输入信号能量:∞<=⎰∞∞-dt t s s E )()(2输入、输出信号频谱函数:dte t s S t j ⎰∞∞--=ωω)()()()()(ωωωS H S o =ωωωπωωd eS H t s tj o ⎰∞-=)()(21)(输出噪声的平均功率:ωωωπωωπd P H d P t n E n n o o⎰⎰∞∞-∞∞-==)()(21)(21)]([22)()()(21)()(2122ωωωπωωπωωd P H d eS H SNR n t j o o⎰⎰∞∞-∞∞-=利用Schwarz 不等式得:ωωωπd P S SNR n o ⎰∞∞-≤)()(212式取等号时,滤波器输出功率信噪比o SNR 最大取等号条件:otj n e P S H ωωωαω-=)()()(*当滤波器输入功率谱密度是2/)(o n N P =ω的白噪声时,MF 的系统函数为:,)()(*o t j e kS H ωωω-=oN k α2=k 为常数1,)(*ωS 为输入函数频谱的复共轭,)()(*ωω-=S S ,也是滤波器的传输函数)(ωH 。
oso N E SNR 2=Es 为输入信号)(t s 的能量,白噪声)(t n 的功率谱为2/o No SNR 只输入信号)(t s 的能量Es 和白噪声功率谱密度有关。
白噪声条件下,匹配滤波器的脉冲响应:)()(*t t ks t h o -=如果输入信号为实函数,则与)(t s 匹配的匹配滤波器的脉冲响应为:)()(t t ks t h o -=k 为滤波器的相对放大量,一般1=k 。
标准文档脉冲压缩雷达的仿真脉冲压缩雷达与匹配滤波的MATLAB仿真姓名:--------学号:----------2014-10-28西安电子科技大学一、 雷达工作原理雷达,是英文Radar 的音译,源于radio detection and ranging 的缩写,原意为"无线电探测和测距",即用无线电的方法发现目标并测定它们的空间位置。
因此,雷达也被称为“无线电定位”。
利用电磁波探测目标的电子设备。
发射电磁波对目标进行照射并接收其回波,由此获得目标至电磁波发射点的距离、距离变化率(径向速度)、方位、高度等信息。
雷达发射机的任务是产生符合要求的雷达波形(Radar Waveform ),然后经馈线和收发开关由发射天线辐射出去,遇到目标后,电磁波一部分反射,经接收天线和收发开关由接收机接收,对雷达回波信号做适当的处理就可以获知目标的相关信息。
但是因为普通脉冲在雷达作用距离与距离分辨率上存在自我矛盾,为了解决这个矛盾,我们采用脉冲压缩技术,即使用线性调频信号。
二、 线性调频(LFM )信号脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。
这种体制采用宽脉冲发射以提高发射的平均功率,保证足够大的作用距离;而接受时采用相应的脉冲压缩算法获得窄脉冲,以提高距离分辨率,较好的解决雷达作用距离与距离分辨率之间的矛盾。
脉冲压缩雷达最常见的调制信号是线性调频(Linear Frequency Modulation )信号,接收时采用匹配滤波器(Matched Filter )压缩脉冲。
LFM 信号的数学表达式:(2.1)其中c f 为载波频率,()t rect T为矩形信号:(2.2)其中BKT=是调频斜率,信号的瞬时频率为()22cT Tf Kt t+ -≤≤,如图(图2.1.典型的LFM信号(a)up-LFM(K>0)(b)down-LFM(K<0))将式1改写为:(2.3)其中(2.4)是信号s(t)的复包络。
目录一.匹配滤波器————————————————P2 一.脉冲压缩雷达———————————————P2 三.线性调频脉冲压缩—————————————P3 二.线性调频脉冲信号匹配滤波—————————P5 五.附录———————————————————P8一.匹配滤波器匹配滤波器是在白噪声背景中检测信号的最佳线性滤波器,其输出信噪比在某个时刻达到最大。
如果一直输入信号u(t),其频谱为U(ω),则可以证明匹配滤波器在频率域的特性为)ex p()()(0*jwt kU H -=ωω式中,U*(ω)为频谱U(ω)的共轭值;k 为滤波器的增益常数;t 0是使滤波器实际上能够实现所必需的延迟时间,在t 0时刻将有信号的最大输出。
同样可以证明,匹配滤波器在时间域的函数,即其脉冲响应为)()(h 0*t t ku t -=式中,u*(t 0-t)为输入信号的镜像,它与输入信号u(t)的波形相同,但从时间t 0开始反转过来。
在对匹配滤波器作理论研究时,延时t 0和增益常数k 可以不予考虑,因此匹配滤波器的上述方程式特性可以简化为)()()()(**t u t h U H -==ωω从上式可以看出:匹配滤波器的传输函数是输入信号频谱的复共轭值,匹配滤波器的脉冲响应是输入信号的镜像函数。
还可以进一步证明,匹配滤波器在输出端给出的最大瞬时信噪比为max 2)(N E N S =式中,N 0是输入噪声的谱密度,它是匹配滤波器输入端单位频带内的噪声功率;E 是输入信号能量:⎰⎰+∞∞-+∞∞-==dt |u(t)|df |U(f)|E 22二.脉冲压缩雷达我们知道,雷达的距离分辨力取决于信号带宽。
在普通脉冲雷达中,雷达信号的时宽带宽积为一敞亮(约为1),因此不能兼顾距离分辨力和速度分辨力两项指标。
近年来,从改进雷达体制方面来矿大作用距离和提高距离分辨力方面已有很大进展。
这种体制就是脉冲压缩(PC )雷达体制,它采用宽脉冲发射以提高发射的平均功率,保证足够的最大作用距离,而在接收时则采用相应的脉冲压缩法获得窄脉冲,以提高距离分辨力,因而能较好地解决作用距离和分辨能力之间的矛盾。
目录一.匹配滤波器————————————————P2 一.脉冲压缩雷达———————————————P2 三.线性调频脉冲压缩—————————————P3 二.线性调频脉冲信号匹配滤波—————————P5 五.附录———————————————————P8一.匹配滤波器匹配滤波器是在白噪声背景中检测信号的最佳线性滤波器,其输出信噪比在某个时刻达到最大。
如果一直输入信号u(t),其频谱为U(ω),则可以证明匹配滤波器在频率域的特性为)ex p()()(0*jwt kU H -=ωω式中,U*(ω)为频谱U(ω)的共轭值;k 为滤波器的增益常数;t 0是使滤波器实际上能够实现所必需的延迟时间,在t 0时刻将有信号的最大输出。
同样可以证明,匹配滤波器在时间域的函数,即其脉冲响应为)()(h 0*t t ku t -=式中,u*(t 0-t)为输入信号的镜像,它与输入信号u(t)的波形相同,但从时间t 0开始反转过来。
在对匹配滤波器作理论研究时,延时t 0和增益常数k 可以不予考虑,因此匹配滤波器的上述方程式特性可以简化为)()()()(**t u t h U H -==ωω从上式可以看出:匹配滤波器的传输函数是输入信号频谱的复共轭值,匹配滤波器的脉冲响应是输入信号的镜像函数。
还可以进一步证明,匹配滤波器在输出端给出的最大瞬时信噪比为max 2)(N E N S =式中,N 0是输入噪声的谱密度,它是匹配滤波器输入端单位频带内的噪声功率;E 是输入信号能量:⎰⎰+∞∞-+∞∞-==dt |u(t)|df |U(f)|E 22二.脉冲压缩雷达我们知道,雷达的距离分辨力取决于信号带宽。
在普通脉冲雷达中,雷达信号的时宽带宽积为一敞亮(约为1),因此不能兼顾距离分辨力和速度分辨力两项指标。
近年来,从改进雷达体制方面来矿大作用距离和提高距离分辨力方面已有很大进展。
这种体制就是脉冲压缩(PC )雷达体制,它采用宽脉冲发射以提高发射的平均功率,保证足够的最大作用距离,而在接收时则采用相应的脉冲压缩法获得窄脉冲,以提高距离分辨力,因而能较好地解决作用距离和分辨能力之间的矛盾。
脉冲压缩雷达的仿真脉冲压缩雷达与匹配滤波的MATLAB仿真姓名:--------学号:----------2014-10-28西安电子科技大学信息对抗技术一、 雷达工作原理雷达,是英文Radar 的音译,源于radio detection and ranging 的缩写,原意为"无线电探测和测距",即用无线电的方法发现目标并测定它们的空间位置。
因此,雷达也被称为“无线电定位”。
利用电磁波探测目标的电子设备。
发射电磁波对目标进行照射并接收其回波,由此获得目标至电磁波发射点的距离、距离变化率(径向速度)、方位、高度等信息。
雷达发射机的任务是产生符合要求的雷达波形(Radar Waveform ),然后经馈线和收发开关由发射天线辐射出去,遇到目标后,电磁波一部分反射,经接收天线和收发开关由接收机接收,对雷达回波信号做适当的处理就可以获知目标的相关信息。
但是因为普通脉冲在雷达作用距离与距离分辨率上存在自我矛盾,为了解决这个矛盾,我们采用脉冲压缩技术,即使用线性调频信号。
二、 线性调频(LFM )信号脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。
这种体制采用宽脉冲发射以提高发射的平均功率,保证足够大的作用距离;而接受时采用相应的脉冲压缩算法获得窄脉冲,以提高距离分辨率,较好的解决雷达作用距离与距离分辨率之间的矛盾。
脉冲压缩雷达最常见的调制信号是线性调频(Linear Frequency Modulation )信号,接收时采用匹配滤波器(Matched Filter )压缩脉冲。
LFM 信号的数学表达式:(2.1)其中c f 为载波频率,()t rect T为矩形信号:(2.2)其中BKT=是调频斜率,信号的瞬时频率为()22cT Tf Kt t+ -≤≤,如图(图2.1.典型的LFM信号(a)up-LFM(K>0)(b)down-LFM(K<0))将式1改写为:(2.3)其中(2.4)是信号s(t)的复包络。
由傅立叶变换性质,S(t)与s(t)具有相同的幅频特性,只是中心频率不同而以,因此,Matlab仿真时,只需考虑S(t)。
以下Matlab程序产生(2.4)的LFM 信号,并作出其时域波形和幅频特性。
%%线性调频信号的产生T=10e-6; %持续时间是10usB=30e6; %调频调制带宽为30MHzK=B/T; %调频斜率Fs=2*B;Ts=1/Fs; %采样频率和采样间隔N=T/Ts;N=T/Ts;t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.^2); %产生线性调频信号subplot(211)plot(t*1e6,real(St));xlabel('时间/us');title('LFM的时域波形');grid on;axis tight;subplot(212)freq=linspace(-Fs/2,Fs/2,N);plot(freq*1e-6,fftshift(abs(fft(St)))); xlabel('频率/MHz'); title('LFM 的频域特性'); grid on;axis tight;(图2.2:LFM 信号的时域波形和频域特性)三、 压缩脉冲的匹配滤波信号()s t 的匹配滤波器的时域脉冲响应为:(3.1)0t 是使滤波器物理可实现所附加的时延。
理论分析时,可令0t =0,重写3.1式,(3.2)将2.1式代入3.2式得:(3.3)图3.1:LFM 信号的匹配滤波如图3.1,()s t 经过系统()h t 得输出信号()o s t ,(3.4)当0t T ≤≤时,(3.5)当0T t -≤≤时,(3.6)合并3.5和3.6两式:(3.7)3.7式即为LFM 脉冲信号经匹配滤波器得输出,它是一固定载频c f 的信号。
当t T ≤时,包络近似为辛克(sinc )函数。
(3.8)图3.2:匹配滤波的输出信号如图3.2,当πB t=±π时,t=±1/B为其第一零点坐标;当πB t=±π/2时,t=±1/2B,习惯上,将此时的脉冲宽度定义为压缩脉冲宽度。
(3.9)LFM信号的压缩前脉冲宽度T和压缩后的脉冲宽度 之比通常称为压缩比D,(3.10)式3.10表明,压缩比也就是LFM信号的时宽频宽积。
由 2.1,3.3,3.7式,s(t),h(t),so(t)均为复信号形式,Matab仿真时,只需考虑它们的复包络S(t),H(t),So(t)。
以下Matlab程序段仿真了图3.1所示的过程,并将仿真结果和理论进行对照。
%%线性调频信号的匹配滤波T=10e-6;B=30e6;K=B/T;Fs=10*B;Ts=1/Fs;N=T/Ts;t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.^2); %产生线性调频信号Ht=exp(-j*pi*K*t.^2); %匹配滤波器Sot=conv(St,Ht); %线性调频信号经过匹配滤波器subplot(211)L=2*N-1;t1=linspace(-T,T,L);Z=abs(Sot);Z=Z/max(Z); %归一化Z=20*log10(Z+1e-6);Z1=abs(sinc(B.*t1)); %sinc函数Z1=20*log10(Z1+1e-6);t1=t1*B;plot(t1,Z,t1,Z1,'r.');axis([-15,15,-50,inf]);grid on;legend('仿真','sinc');xlabel('时间sec \times\itB');ylabel('振幅,dB');title('线性调频信号经过匹配滤波器');subplot(212) %放大N0=3*Fs/B;t2=-N0*Ts:Ts:N0*Ts;t2=B*t2;plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');axis([-inf,inf,-50,inf]);grid on;set(gca,'Ytick',[-13.4,-4,0],'Xtick',[-3,-2,-1,-0.5,0,0.5,1,2,3]);xlabel('时间sec \times\itB');ylabel('振幅,dB');title('线性调频信号经过匹配滤波器(放大)');结果:图3.3:线性调频信号的匹配滤波上图中,时间轴进行了归一化,(t/(1/B)=t x B)。
图中反映出理论与仿真结果吻合良好。
第一零点出现在±1(即±1/B)处,此时相对幅度-13.4dB。
压缩后的脉冲宽度近似为1/B(±1/2B),此时相对幅度-4dB,这理论分析(图3.2)一致。
四、Matlab仿真1.任务:对以下雷达系统仿真。
雷达发射信号参数:幅度:1.0信号波形:线性调频信号频带宽度:30MHz脉冲宽度:10us中心频率:1GHzHz雷达接收方式:正交解调接收距离门:10Km~15Km目标:Tar1:10.5KmTar2:11KmTar3:12KmTar4:12Km+5mTar5:13KmTar6:13Km+2m2.系统模型:结合以上分析,用Matlab仿真雷达发射信号,回波信号,和压缩后的信号的复包络特性,其载频不予考虑(实际中需加调制和正交解调环节),仿真信号与系统模型如下图。
图4.1:雷达仿真等效信号与系统模型3.线性调频脉冲压缩雷达仿真程序LFM_radar仿真程序模拟产生理想点目标的回波,并采用频域相关方法(以便利用FFT)实现脉冲压缩。
函数LFM_radar的参数意义如下:T:LFM信号的持续脉宽;B:LFM信号的频带宽度;Rmin:观测目标距雷达的最近位置;Rmax:观测目标距雷达的最远位置;R:一维数组,数组值表示每个目标相对雷达的距离;RCS:一维数组,数组值表示每个目标的雷达散射截面。
在Matlab指令窗中输入:LFM_radar(10e-6,30e6,10000,15000,[10500,11000,12000,12005,13000,13002],[1,1,1,1,1,1])得到的仿真结果如下图。
五、心得通过这次使用Matlab对脉冲压缩雷达的仿真,让我充分理解到了脉冲压缩雷达的工作原理,以及脉冲压缩雷达与普通脉冲雷达的差异,这让我对与雷达原理这门课有了更加深入的理解,对于匹配滤波的深入了解,使得在课堂中没有充分理解的地方清晰的展现在眼前。
这次实验不仅仅会促进我雷达原理课程的学习,也为我以后学习雷达专业提供了一种可靠的方法。
六、附录Matlab代码(LFM_radar.m)%%脉冲压缩雷达仿真function LFM_radar(T,B,Rmin,Rmax,R,RCS)if nargin==0T=10e-6; %脉冲持续时间10usB=30e6; %频带宽度30MHzRmin=10000;Rmax=15000; %作用范围R=[10500,11000,12000,12008,13000,13002]; %目标位置RCS=[1 1 1 1 1 1]; %雷达散射面end%%参数设定C=3e8; %设定速度为光速K=B/T; %调频斜率Rwid=Rmax-Rmin; %距离Twid=2*Rwid/C; %时间Fs=5*B;Ts=1/Fs; %采样频率和采样间隔Nwid=ceil(Twid/Ts);%%回波t=linspace(2*Rmin/C,2*Rmax/C,Nwid); %接收范围(2*Rmin/C < t < 2*Rmax/C)M=length(R); %目标数量td=ones(M,1)*t-2*R'/C*ones(1,Nwid);Srt=RCS*(exp(j*pi*K*td.^2).*(abs(td)<T/2)); %雷达回波%%利用FFT和IFFT进行数字信号处理Nchirp=ceil(T/Ts); %多脉冲持续时间Nfft=2^nextpow2(Nwid+Nwid-1);Srw=fft(Srt,Nfft); %雷达回波的fft计算t0=linspace(-T/2,T/2,Nchirp);St=exp(j*pi*K*t0.^2); %线性调频信号Sw=fft(St,Nfft); %线性调频信号的fft计算Sot=fftshift(ifft(Srw.*conj(Sw))); %脉冲压缩后的信号N0=Nfft/2-Nchirp/2;Z=abs(Sot(N0:N0+Nwid-1));Z=Z/max(Z);Z=20*log10(Z+1e-6);%产生图像subplot(211)。