CN104020454B - 利用发射泄露对太赫兹fmcw成像雷达实时校正的方法 - Google Patents

利用发射泄露对太赫兹fmcw成像雷达实时校正的方法 Download PDF

Info

Publication number
CN104020454B
CN104020454B CN201410095062.5A CN201410095062A CN104020454B CN 104020454 B CN104020454 B CN 104020454B CN 201410095062 A CN201410095062 A CN 201410095062A CN 104020454 B CN104020454 B CN 104020454B
Authority
CN
China
Prior art keywords
signal
frequency
target
spectrum
radar
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201410095062.5A
Other languages
English (en)
Other versions
CN104020454A (zh
Inventor
成彬彬
江舸
邓贤进
张健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Electronic Engineering of CAEP
Original Assignee
Institute of Electronic Engineering of CAEP
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Electronic Engineering of CAEP filed Critical Institute of Electronic Engineering of CAEP
Priority to CN201410095062.5A priority Critical patent/CN104020454B/zh
Publication of CN104020454A publication Critical patent/CN104020454A/zh
Application granted granted Critical
Publication of CN104020454B publication Critical patent/CN104020454B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提出了利用发射泄露对太赫兹FMCW成像雷达实时校正的方法,利用对FMCW成像雷达系统不利但又难以消除的发射泄露信号进行收发信道进行校正,变不利为有利,校正可以对每个脉冲实时完成,能做到实时校正,不需要多长时间预热,并且校正率较高。

Description

利用发射泄露对太赫兹FMCW成像雷达实时校正的方法
技术领域
本发明涉及太赫兹雷达与成像技术,具体是一种利用发射泄露对太赫兹FMCW成像雷达实时校正的方法。
背景技术
太赫兹FMCW雷达由于其超高的成像分辨力而得到广泛的关注,然而,由于经过多次倍频和放大,这些器件的非线性会使雷达发射支路的调频连续波信号及接收支路混频器的本振信号都发生畸变,形成幅度和相位调制,分别用dj(t)和A(t)来表示,他们的存在会影响信号带内平坦度和与频率有关的相位延迟。
发射和本振信号的相位和幅度调制也会被引入到中频信号,因此,一个单一目标的检测信号不再是理想的正弦波,而是变为:
式中幅度和相位的调制项来自雷达信号和接收本振的乘积:
AIF(t,R)=ALO(t)·AT(t-2R/c)
与距离相关的幅度和相位调制来自发射信号到达和目标回波信号的传输时间。
对于浅度调制,中频信号包含一个频率为fIF的单频信号,其频谱会被一个表达式为的频谱宽度展宽。该展宽会使雷达的距离分辨力恶化,并使太赫兹雷达图像变坏。需要对此非线性调制进行校正。
目前常用的校正方法为事先校正法,在对实际的目标进行测试前,先在目标位置附近R0处放置一个强反射体,利用太赫兹雷达对该目标进行照射,并存储采集到的数字中频信号,
接下来雷达所有的测量在进行谱分析之前都要扣除这一校正波形,因此,检测到的中频信号的幅度和相位调制表达式需要做如下修正:
A IF ( t , R ) → A IF ( t , R ) A 0 ( t , R 0 )
信号的频率也将偏移一个已知量DfIF=2KR0/c,该值可以加回去从而得到绝对距离。只要幅度调制的补偿达到统一并且相位调制补偿到0,距离的展宽将会不明显。
该方法具有如下缺陷:
1.使用不方便,需随机器携带一个强反射体,并且校正前需将机器预热到实际工作时的状态,一般为半小时。
2.校正无法做到实时,工作条件(如温度、目标距离等)一旦改变,需重新校正。
由于太赫兹FMCW成像雷达多应用于对近程目标进行高分辨力成像(几米到几十米),且一般采用收发分离的架构,发射信号不可避免的会泄露到接收机中,在一般系统中,需尽量进行隔离,以较小对接收机的影响。本设计中,正是利用该泄露信号完成对太赫兹雷达收发信道的实时校正。
发明内容
本发明提出了一种利用发射泄露信号对太赫兹FMCW成像雷达进行实时校正的方法,利用对FMCW成像雷达系统不利但又难以消除的发射泄露信号进行收发信道进行校正,变不利为有利,校正可以对每个脉冲实时完成,能做到实时校正,不需要多长时间预热,并且校正率较高。
本发明的技术方案如下:
利用发射泄露对太赫兹FMCW成像雷达实时校正的方法,其特征在于:发射信号经发射天线的旁瓣耦合到接收天线的旁瓣上,信号被接收天线收集,假定发射天线到接收天线电磁波的传播路径长度为2R0,根据FMCW雷达理论,耦合到接收机的发射信号会在FMCW雷达的中频上产生一个频率为DfIF0=2KR0/c的频移,其中,DfIF0为耦合产生的中频偏移频率,K为调频斜率,c为电磁波传播速度;同时,发射信号照射到距离为R的被测目标上经目标散射的回波信号也被接收天线收集,同样在会中频上产生一个频率为DfIF=2KR/c的频移,其中DfIF为目标产生的中频偏移频率,R为被测成像目标与雷达的中心距离。由于目标可能不是点目标,故DfIF可能不是一个单频信号,而是具有一定的带宽。由于在一般的应用条件下,满足R>>R0,从而有DfIF>>DfIF0,故在频域上,耦合产生的中频偏移频率DfIF0和目标产生的中频偏移频率DfIF可通过滤波实现分离,从而分别获得耦合信号产生的中频信号S0(t,R0)(I,Q两路,复信号)和由于目标散射产生的回波中频信号S(t,R)(I,Q两路,复信号),计算S(t,R)/S0(t,R0),获得校正后的目标回波信号SIF(t,R')=S(t,R)/S0(t,R0),从而完成对收发信道非线性的实时校正。其中,t为时间,R'为经过校正后计算获得的目标距离,目标实际距离R=R'+R0
其中,发射天线到接收天线的传播路径长度为2R0(该值由收发天线间距决定,一般为几十厘米),实际成像目标的中心距离为R(一般为几米到几十米)。
该方法的具体处理流程如下:
A.FMCW雷达开机并发射太赫兹信号,照射被成像目标;
B.对信号处理单元中A/D采集每个脉冲得到的中频信号,进行N点复数的快速傅里叶变换(Fast Fourier Transform,FFT)得到中频信号频谱H(f);
C.取出中频信号频谱H(f)中频率范围在[0,f0]∪[(fs-f0),fs](DfIF0<f0<DfIF)范围内的信号,获得泄露产生的频谱HR0(f),即将H(f)频谱中f0~(fs-f0)部分的频谱值置0得到HR0(f),其中,fs为信号处理系统采样率,f0为分离低通滤波器的截止频率,DfIF0<f0<DfIF,其大小可根据收发天线的距离和目标的位置进行选取,具体方法为:根据雷达调频斜率参数K和收发天线距离2R0估算DfIF0,根据雷达调频斜率参数K和目标位置R估算DfIF,f0取两者之间的值均可,简便起见,可取f0=(DfIF+DfIF0)/2。
D.取出H(f)频谱中频率范围在区间[f0,(fs-f0)]内的频谱值,H(f)频谱中其余部分置0得到目标回波产生的频谱HR(f);
E.分别对泄露产生的频谱HR0(f)和目标回波产生的频谱HR(f)进行N点逆快速傅里叶变换(Inversed FFT,IFFT)得到耦合信号产生的中频信号S0(t,R0)和目标散射产生的回波中频信号S(t,R);
F.计算S(t,R)/S0(t,R0),获得校正后的目标回波信号SIF(t,R')=S(t,R)/S0(t,R0);
G.利用校正后的数据SIF(t,R')完成成像,目标上散射点的相对位置不变,距离雷达的位置通过R=R'+R0完成补偿。
本发明的有益效果如下:
1)利用对FMCW成像雷达系统不利但又难以消除的发射泄露信号对收发信道的非线性进行校正,变不利为有利,这也是本设计的创新点。
2)校正可以对每个脉冲实时完成,环境变化对信道和回波造成的漂移可被校正在内,不需再单独考虑或者重新校正。
3)除了增加部分的处理代码,无需额外增加或者改动硬件即可实现实时校正,也不需自带校正用的强散射体。
附图说明
图1为采用本发明进行的实时校正示意图
具体实施方式
如图1所示,驱动发射链的点频源和驱动接收链的点频源共用参考信号,以保证收发系统的相参性,基带信号一般为线性调频信号,通过与发射链驱动点频源和接收链驱动点频源混频分别调制到发射和接收链路,带通滤波后作为发射倍频放大链和接收倍频放大链的驱动信号。发射倍频放大链的输出经发射喇叭天线辐射出去形成太赫兹发射信号,发射信号经发射天线的旁瓣耦合到接收天线的旁瓣上,同时发射信号经发射天线的主瓣照射到目标上,经目标散射的回波信号也会进入接收天线主瓣,旁瓣耦合信号和目标回波信号均通过接收天线收集,假定发射天线到接收天线的传播路径长度为2R0(该值由收发天线间距决定,一般为几十厘米),则从太赫兹波传播路径上来讲,相当于在R0处放置一反射目标,对于FMCW雷达来说,会在中频上产生一个频率为DfIF0=2KR0/c的频移;假设实际成像目标的中心距离为R(一般为几米到几十米),则目标产生的频率偏移在DfIF=2KR/c附近,由于R>>R0,故在频域上,DfIF0和DfIF不会产生混叠,可以通过滤波的方法将其分离。从而分别获得耦合信号产生的中频信号S0(t,R0)(I,Q两路,复信号)和由于目标散射产生的回波中频信号S(t,R)(I,Q两路,复信号),计算S(t,R)/S0(t,R0),获得校正后的目标回波信号SIF(t,R')=S(t,R)/S0(t,R0),从而完成对收发信道非线性的实时校正。其中,t为时间,R'为经过校正后计算获得的目标距离,目标实际距离R=R'+R0
该方法的具体处理流程如下:
A.FMCW雷达开机并发射太赫兹信号,照射被成像目标;
B.对信号处理单元中A/D采集每个脉冲得到的中频信号,进行N点复数的快速傅里叶变换(Fast Fourier Transform,FFT)得到中频信号频谱H(f);
C.取出中频信号频谱H(f)中频率范围在[0,f0]∪[(fs-f0),fs](DfIF0<f0<DfIF)范围内的信号,获得泄露产生的频谱HR0(f),即将H(f)频谱中f0~(fs-f0)部分的频谱值置0得到HR0(f),其中,fs为信号处理系统采样率,f0为分离低通滤波器的截止频率,DfIF0<f0<DfIF,其大小可根据收发天线的距离和目标的位置进行选取,具体方法为:根据雷达调频斜率参数K和收发天线距离2R0估算DfIF0,根据雷达调频斜率参数K和目标位置R估算DfIF,f0取两者之间的值均可,简便起见,可取f0=(DfIF+DfIF0)/2。
D.取出H(f)频谱中频率范围在区间[f0,(fs-f0)]内的频谱值,H(f)频谱中其余部分置0得到目标回波产生的频谱HR(f);
E.分别对泄露产生的频谱HR0(f)和目标回波产生的频谱HR(f)进行N点逆快速傅里叶变换(Inversed FFT,IFFT)得到耦合信号产生的中频信号S0(t,R0)和目标散射产生的回波中频信号S(t,R);
F.计算S(t,R)/S0(t,R0),获得校正后的目标回波信号SIF(t,R')=S(t,R)/S0(t,R0);
G.利用校正后的数据SIF(t,R')完成成像,目标上散射点的相对位置不变,距离雷达的位置通过R=R'+R0完成补偿。

Claims (2)

1.利用发射泄露对太赫兹FMCW成像雷达实时校正的方法,其特征在于:发射信号经发射天线的旁瓣耦合到接收天线的旁瓣上,发射信号被接收天线收集,耦合到接收机的发射信号会在FMCW雷达的中频上产生一个频率为ΔfIF0=2KR0/c的频移,其中,2R0为设定的发射天线到接收天线电磁波的传播路径长度,ΔfIF0为耦合产生的中频偏移频率,K为调频斜率,c为电磁波传播速度;同时,发射信号照射到距离为R的被测目标上经目标散射的回波信号也被接收天线收集,同样在会中频上产生一个频率为ΔfIF=2KR/c的频移,其中ΔfIF为目标产生的中频偏移频率,R为被测成像目标与雷达的中心距离;
当R>>R0时,ΔfIF>>ΔfIF0,故在频域上,耦合产生的中频偏移频率ΔfIF0和目标产生的中频偏移频率ΔfIF通过滤波实现分离,从而分别获得耦合信号产生的中频信号S0(t,R0)和由于目标散射产生的回波中频信号S(t,R),中频信号S0(t,R0)和回波中频信号S(t,R)均为I、Q两路复信号,计算S(t,R)/S0(t,R0),获得校正后的目标回波信号SIF(t,R')=S(t,R)/S0(t,R0),从而完成对收发信道非线性的实时校正;其中,t为时间,R'为经过校正后计算获得的目标距离,目标实际距离R=R'+R0
2.根据权利要求1所述的利用发射泄露对太赫兹FMCW成像雷达实时校正的方法,其特征在于具体处理流程如下:
A.FMCW雷达开机并发射太赫兹信号,照射被成像目标;
B.对信号处理单元中A/D采集每个脉冲得到的中频信号,进行N点复数的快速傅里叶变换得到中频信号频谱H(f);
C.取出中频信号频谱H(f)中频率范围在[0,f0]∪[(fs-f0),fs]范围内的信号,其中ΔfIF0<f0<ΔfIF,获得泄露产生的频谱HR0(f),即将H(f)频谱中f0~(fs-f0)部分的频谱值置0得到HR0(f),其中,fs为信号处理系统采样率,f0为分离低通滤波器的截止频率,ΔfIF0<f0<ΔfIF,其大小可根据收发天线的距离和目标的位置进行选取,具体方法为:根据雷达调频斜率参数K和收发天线距离2R0估算ΔfIF0,根据雷达调频斜率参数K和目标位置R估算ΔfIF,f0可取为f0=(ΔfIF+ΔfIF0)/2;
D.取出H(f)频谱中频率范围在区间[f0,(fs-f0)]内的频谱值,H(f)频谱中其余部分置0得到目标回波产生的频谱HR(f);
E.分别对泄露产生的频谱HR0(f)和目标回波产生的频谱HR(f)进行N点逆快速傅里叶变换得到耦合信号产生的中频信号S0(t,R0)和目标散射产生的回波中频信号S(t,R);
F.计算S(t,R)/S0(t,R0),获得校正后的目标回波信号SIF(t,R')=S(t,R)/S0(t,R0);
G.利用校正后的数据SIF(t,R')完成成像,目标上散射点的相对位置不变,距离雷达的位置通过R=R'+R0完成补偿。
CN201410095062.5A 2014-03-14 2014-03-14 利用发射泄露对太赫兹fmcw成像雷达实时校正的方法 Expired - Fee Related CN104020454B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410095062.5A CN104020454B (zh) 2014-03-14 2014-03-14 利用发射泄露对太赫兹fmcw成像雷达实时校正的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410095062.5A CN104020454B (zh) 2014-03-14 2014-03-14 利用发射泄露对太赫兹fmcw成像雷达实时校正的方法

Publications (2)

Publication Number Publication Date
CN104020454A CN104020454A (zh) 2014-09-03
CN104020454B true CN104020454B (zh) 2016-08-17

Family

ID=51437311

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410095062.5A Expired - Fee Related CN104020454B (zh) 2014-03-14 2014-03-14 利用发射泄露对太赫兹fmcw成像雷达实时校正的方法

Country Status (1)

Country Link
CN (1) CN104020454B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104330154B (zh) * 2014-10-16 2016-08-24 中国电子科技集团公司第五十研究所 窄线宽的太赫兹探测器
CN105842685B (zh) * 2016-03-18 2019-02-12 浙江大华技术股份有限公司 一种多目标雷达探测方法
CN106019246B (zh) * 2016-04-08 2018-08-17 中国工程物理研究院电子工程研究所 一种基于内部目标的主动太赫兹快速安检仪实时校正方法
CN106199188B (zh) * 2016-07-20 2019-01-18 中国科学院紫金山天文台 一种利用环形器校准矢量场测量中移动电缆相位变化的装置及方法
CN108459305A (zh) * 2017-02-20 2018-08-28 北京雷测科技有限公司 太赫兹雷达收发系统和单发多收太赫兹相参雷达
CN108459326A (zh) * 2017-02-20 2018-08-28 北京雷测科技有限公司 太赫兹雷达收发系统和单发多收太赫兹相参雷达
CN109765533A (zh) * 2018-12-29 2019-05-17 成都聚利中宇科技有限公司 一种通道实时校准方法及系统及成像设备、检测设备
CN110907902B (zh) * 2019-09-23 2023-12-05 成都锦江电子系统工程有限公司 天气雷达定标方法
CN112462336B (zh) * 2020-10-10 2022-11-29 东华大学 一种fmcw异物检测雷达泄漏信号的自适应消除方法
CN113671495B (zh) * 2021-08-18 2024-02-23 上海无线电设备研究所 基于Zynq平台的太赫兹雷达探测系统以及方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7773205B2 (en) * 2007-06-06 2010-08-10 California Institute Of Technology High-resolution three-dimensional imaging radar
US8144052B2 (en) * 2008-10-15 2012-03-27 California Institute Of Technology Multi-pixel high-resolution three-dimensional imaging radar

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
0.14THz高分辨力成像雷达信号处理;成彬彬 等;《强激光与粒子束》;20130630;第25卷(第6期);全文 *
Real-Time Imaging With a 140 GHz Inverse Synthetic Aperture Radar;BinBin Cheng 等;《IEEE TRANSACTIONS ON TERAHERTZ SCIENCE AND TECHNOLOGY》;20130930;第3卷(第5期);全文 *
THz Imaging Radar for Standoff Personnel Screening;Ken B.Cooper 等;《IEEE TRANSACTIONS ON TERAHERTZ SCIENCE AND TECHNOLOGY》;20110930;第1卷(第1期);正文第2节A、C部分 *

Also Published As

Publication number Publication date
CN104020454A (zh) 2014-09-03

Similar Documents

Publication Publication Date Title
CN104020454B (zh) 利用发射泄露对太赫兹fmcw成像雷达实时校正的方法
EP3153876B1 (en) Method of system compensation to reduce the effects of self interference in frequency modulated continuous wave altimeter systems
EP3220162B1 (en) Radar device and position-determination method
CN107607923A (zh) 基于lfmcw雷达的振动监测系统与信号处理方法
JP5197138B2 (ja) マルチスタティックレーダ装置
CN106932765A (zh) 具有相位噪声估计的雷达设备
CN103529444A (zh) 一种车载毫米波雷达动目标识别器及识别方法
CN104076352A (zh) 低截获测速方法及雷达装置
CN106353748A (zh) 用于fmcw雷达测距系统的信号处理装置及方法
US11774551B2 (en) Method and device for compensating for interfering influences
US20160245913A1 (en) Apparatus and method for measuring precipitation in the atmosphere using k-band frequency-modulated continuous wave (fmcw) weather radar system
CN101464514B (zh) 步进频率雷达系统的校准方法与校准处理器
CN103064065B (zh) 双周期夹断式伪随机码的波形设计和对回波的信号处理方法
CN109001729B (zh) 太赫兹成像中调频连续波线性度实时校准方法及其系统
CN103954937B (zh) 一种宽范围高精度微波测距雷达系统设计方法
CN101788671A (zh) 应用于外差探测啁啾调幅激光测距装置的多周期调制方法
KR100940918B1 (ko) 음영지역 탐지를 위한 펄스 압축 레이더의 펄스 파형 송신방법, 이를 이용한 펄스 압축 레이더 및 레이더 네트워크
Ma et al. High-resolution microwave photonic radar with sparse stepped frequency chirp signals
Suksmono et al. Signal processing of range detection for SFCW radars using Matlab and GNU radio
JP2016156620A (ja) レーダ装置
JP5104425B2 (ja) 距離測定方法及び距離測定装置
US20070063888A1 (en) Single side band radar
CN104007438B (zh) 一种基于回波数据的步进频率合成孔径雷达频差估计方法
CN102998659B (zh) 基于脉间调制的多普勒频谱赋形方法及系统
KR101235057B1 (ko) 디지털고주파기억장치

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20170314