CN111708023B - 毫米波多发多收线阵的相位不一致性校正方法及成像*** - Google Patents
毫米波多发多收线阵的相位不一致性校正方法及成像*** Download PDFInfo
- Publication number
- CN111708023B CN111708023B CN202010840929.0A CN202010840929A CN111708023B CN 111708023 B CN111708023 B CN 111708023B CN 202010840929 A CN202010840929 A CN 202010840929A CN 111708023 B CN111708023 B CN 111708023B
- Authority
- CN
- China
- Prior art keywords
- millimeter wave
- transmitting
- array
- echo signal
- receiving
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
- G01S7/4056—Means for monitoring or calibrating by simulation of echoes specially adapted to FMCW
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种毫米波多发多收线阵的相位不一致性校正方法及成像***,利用毫米波多发多收线阵录取目标回波信号;采取金属平板作为参考目标,录取参考目标回波信号;结合金属平板的尺寸、位置和毫米波多发多收线阵的阵元分布方式,采用电磁计算软件得到相同场景下金属平板的仿真回波信号;利用参考目标回波信号和仿真回波信号对目标回波信号进行补偿,得到相位不一致性校正后的目标回波信号,校正后的目标回波信号直接用于目标的成像与重构。本发明采用电磁计算软件仿真的金属平板回波信号作为距离恢复的补偿数据,无需计算需要补偿的真实距离,可直接恢复多发多收稀疏阵列成像场景下真实的成像几何关系,保证成像质量。
Description
技术领域
本发明属于雷达信号处理技术领域,具体地涉及一种毫米波多发多收线阵的相位不一致性校正方法及成像***。
背景技术
毫米波通常指频率在30GHz到300GHz(对应波长1mm-10mm)之间的电磁波,是近年来研究最热门的电磁频谱。相比于微波以及更低频段的电磁波,毫米波雷达载频更高,波长更短,更容易实现大带宽的雷达信号,因此可以实现更高距离和方位分辨率的成像,在安全检查、无损探测等领域具有重大的应用前景。
近年来,随着毫米波源和检测器件的重大突破,已经有多个国内外研究机构开发了毫米波安检成像***,其按照阵列结构主要可以分为两大类:准单站阵列成像***和稀疏多站阵列成像***。准单站阵列成像***的发射阵元间隔和接收阵元间隔都是均匀的,收发线阵左右错开半个阵元间隔,形成准单站分布。这种阵列***的优势是可以实现对目标的全方位波束覆盖,***结构、成像算法简单,但是阵元通道数多,成本较高。稀疏多站阵列成像***的发射阵元和接收阵元是按照一定的设计规则稀疏分布的,通过等效阵元中心分布或者采用相应的稀疏阵列成像算法,可以以较少的阵元数目实现与准单站阵列成像***相同的功能,极大的降低了硬件成本和***功耗,是目前毫米波成像***的主流发展趋势。
对于毫米波阵列成像***,由于器件制造差异和传输通道电尺寸长度的影响,相同的源信号通过不同发射通道和接收通道,对应的射频信号和本振信号的初相往往是不一致的。这个不一致的相位会破坏阵列录取目标回波的横向多普勒信息,进而导致成像结果散焦甚至模糊。对于传统的准单站毫米波阵列***,通常采用一个金属平板作为参考目标校正相位的不一致性,再补偿金属平板到线阵之间的直线距离从而恢复真实的成像几何关系。对于多发多收的稀疏阵列,收发阵元之间存在一定的距离,导致准单站近似失效,在采取金属平板作为参考目标校正相位的不一致性后,需要补偿的距离大于金属平板到线阵之间的距离,此时准确的计算收发阵元之间需要补偿的距离存在较大的难度。以一种四发八收线阵为例,其阵元分布如图1所示。在图1中,圆圈为接收阵元,叉号为发射阵元,线阵总长度为0.17米。毫米波***工作在超外差模式,发射信号为宽带线性调频信号,脉冲宽度为T P 。四个发射阵元依次发射该信号,阵元之间发射信号的时间间隔为T A (T A >T P ),每个阵元发射信号的同时,八个接收通道同步开启同时接收目标的回波信号,以此保证信号之间的相参性。由于天线等器件水平差异和传输通道电尺寸长度不一致的影响,相同的源信号通过不同发射通道和接收通道,对应的射频信号和本振信号的初相也存在一定的差异。不一致的相位会破坏阵列录取目标回波的真实横向多普勒信息,导致成像结果散焦甚至模糊。
目前,大部分文献的毫米波稀疏线阵由一个发射通道和一个接收通道通过扫描等效得到,此时不存在通道间的相位不一致性问题。对于实孔径多发多收线阵的相位不一致性,可以通过测量的手段测出通道间的相位差异,在实际应用中对该相位差进行预补偿处理,但是测量不可避免会存在一定误差,而且当收发阵元非常多的时候,工作量巨大,精确的测量不易实现。现有技术中也有采用一个金属小球作为参考目标对***初始相位进行校正,这种校正方法只适用于阵列长度远小于参考目标和成像目标距离的成像场景,不需要再进一步补偿阵元到参考目标之间的距离。当阵列尺寸与目标距离相当时,这种方法并不适用。此外,金属小球属于滑动型散射中心,其双站散射强度很弱,会导致参考数据的信噪比较低,进而引入新的误差相位。
发明内容
针对毫米波多发多收线阵成像***在实际应用中存在的阵元通道相位不一致的问题,本发明提出了一种毫米波多发多收线阵的相位不一致性校正方法及成像***,保证毫米波多发多收线阵的成像性能。
为实现上述技术目的,本发明采用的具体技术方案如下:
毫米波多发多收线阵的相位不一致性校正方法,包括:
(1)利用毫米波多发多收线阵录取目标回波信号S r 。
(2)利用金属平板作为参考目标,金属平板与毫米波多发多收线阵的扫描平面平行,录取金属平板的回波信号,将其作为相位不一致性的参考目标回波信号S ref ;
(3)利用电磁计算软件仿真得到与步骤(2)相同场景下的金属平板的仿真回波信号S fek ,其中相同场景是指仿真中的毫米波多发多收线阵的阵元分布、***的工作频率、金属平板的尺寸以及设置位置均与步骤(2)相同;
作为优选方案,本发明所述毫米波多发多收线阵为四发八收的毫米波多发多收线阵。
作为优选方案,本发明所述毫米波多发多收线阵顺序发射宽带调频连续波。
作为优选方案,本发明步骤(2)中所设置的金属平板需要满足毫米波多发多收线阵中的发射阵元的主瓣波束照射不到金属平板的边缘。
作为优选方案,本发明步骤(3)中利用电磁计算软件FEKO仿真得到与步骤(2)相同场景下的金属平板的电磁散射回波信号。
作为优选方案,本发明步骤(1)中,利用毫米波多发多收线阵录取目标回波信号,其中毫米波多发多收线阵中任意收发阵元组合的目标回波信号表达式为:
其中坐标为(x p ,y p ,z p )点目标与毫米波多发多收线阵中发射阵元天线和接收阵元天线的距离分别为R Tm和R Rn ,其中m为发射阵元天线序号,n为接收阵元天线序号,目标散射系数为σ(x p ,y p ,z p ),位置为(x Tm ,z)的发射阵元天线到目标的距离为,位置为(x Rn ,z)的发射阵元天线到目标的距离为;,t为距离快时间,f c 为雷达载频,γ为调频率,c为光速,exp表示e为底的指数函数,θ m 为第m个发射阵元天线的附加相位,为第n个接收阵元天线的附加相位。
本发明步骤(2)中,不考虑电磁波传输路线的距离,设毫米波多发多收线阵中任意收发阵元之间的传输距离为R ref_mn ,毫米波多发多收线阵中任意收发阵元组合的参考目标回波信号的表达式为:
其中σ mn 对应第m个发射阵元天线发射信号和第n个接收阵元天线接收信号时金属平板的散射强度。
本发明步骤(3)中,电磁计算软件仿真得到的毫米波多发多收线阵中任意收发阵元组合的金属平板的电磁散射回波信号的表达式为:
本发明步骤(4)中,相位不一致性校正后的目标回波信号S t_comp ,其中校正后的毫米波多发多收线阵中任意收发阵元组合的目标回波信号表达式为:
本发明还提供一种毫米波多发多收线阵成像***,采用上述毫米波多发多收线阵的相位不一致性校正方法获取相位不一致性校正后的目标回波信号用于目标的成像与重构。
本发明的有益效果如下:
本发明采取金属平板作为参考目标录取参考目标回波信号,在录取参考目标回波信号后,测量金属平板到毫米波多发多收线阵的直线距离,相比于大量不同收发阵元位置与金属平板复杂的几何关系,金属平板到毫米波多发多收线阵的直线距离可以很容易的直接测量得到。获得直线距离后,结合金属平板的尺寸和毫米波多发多收线阵的阵元分布方式,可以采用电磁计算软件仿真得到相同场景下金属平板的电磁散射回波信号,电磁计算数据结果的相位信息包含了需要补偿的参考目标距离。此外,电磁计算数据结果的幅度信息对应着金属平板真实的电磁散射强度,利用该电磁计算数据结果可以同时实现参考目标引入的偏移距离和散射强度不一致性的补偿。本发明提供的方法是最准确的毫米波多发多收线阵的相位不一致性校正方法。
本发明方法不需要繁琐的测量毫米波多发多收线阵的相位不一致的数值。在采取光滑金属平板作为参考目标时,对参考目标引入的偏移距离通过精确的相同场景下的电磁计算数据来补偿恢复,无需计算毫米波多发多收线阵其阵元分布与参考目标复杂几何关系下的补偿距离,本发明的补偿精度更高,更易于实现,极大的提升了工作效率。另一方面,本发明提出采用电磁计算软件仿真的金属平板回波信号作为距离恢复的补偿数据,无需计算需要补偿的真实距离,可以直接恢复多发多收稀疏阵列成像场景下真实的成像几何关系,进而保证***的成像质量。
附图说明
图1是四发八收线阵阵元分布示意图;
图2是四发八收阵列的成像几何关系示意图;
图3是金属平板回波距离示意图;
图4是FEKO电磁计算场景示意图;
图5是本发明流程图。
具体实施方式
为了使本发明的技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用于解释本发明,并不用于限定本发明。
实施例1:
本实施例提供一种毫米波多发多收线阵的相位不一致性校正方法,包括:
S1: 目标宽带回波获取。
本实施例针对的是宽带雷达***,毫米波阵列雷达顺序发射宽带调频连续波。
以四发八收的毫米波多发多收线阵为例,其成像的几何关系如图2所示。毫米波多发多收线阵任意位置的发射信号的都具有如下式所示的模式。
其中,t为距离快时间,T P 为调频连续波雷达的脉冲宽度,f c 为雷达载频,γ为调频率,rect( )表示方脉冲。假设坐标为(x p ,y p ,z p )的点目标与发射阵元天线和接收阵元天线的距离分别为R Tm和R Rn ,其中m为发射阵元序号,n为接收阵元序号。目标散射系数为σ(x p ,y p ,z p ),任意收发阵元组合的回波可以写为:
以发射信号作为本振信号,***通过解线频调的方式接收得到的目标回波信号,在消除视频相位项并且忽略信号的包络后,此时任意收发阵元组合的回波信号表达式为:
也可以表示为
S2: 参考目标宽带回波获取。
参考目标为一个金属平板,与图2中稀疏线阵的扫描平面平行放置,距离向坐标为y ref 。金属平板的尺寸需要保证毫米波多发多收线阵中任意发射天线的波束照射不到金属平板的边缘。根据镜面反射定理,每一组收发阵元得到的金属平板的回波信号只存在一个固定的距离。以任意一组收发阵元为例,其距离示意图如图3所示。其中入射波与水平面的夹角与反射波与水平面的夹角相同,此时接收阵元才能够接收到对应发射阵元经金属平板散射后的回波信号。
不考虑电磁波传输路线的距离,设任意收发阵元之间的传输距离为R ref_mn ,毫米波多发多收线阵中任意收发阵元组合的参考目标回波信号可表达为
其中σ mn 对应第m个发射阵元天线发射信号和第n个接收阵元天线接收信号时金属平板的散射强度。
S3: 利用电磁计算软件FEKO仿真和图3完全相同的场景下的金属平板的散射回波,其仿真场景如图4所示。箭头为发射阵元位置,叉号为接收阵元位置,目标为纯金属平板,采用物理光学法计算目标的电磁散射回波,单次计算可得到一个发射阵元发射信号,8个接收阵元同时接收的回波信号。通过MATLAB软件控制FEKO软件的.PRE文件可以改变发射阵元的位置,通过多次计算即可得到完整的稀疏阵列的回波信号。
由于电磁计算软件计算目标散射数据时,仿真条件均处于理想条件,不存在通道间的相位不一致的问题。此时电磁计算的毫米波多发多收线阵中任意收发阵元组合的金属平板的电磁散射回波信号可以表达为:
S4: 利用参考数据和电磁计算数据对目标的回波信号进行补偿,可得
此时完成了毫米波多发多收线阵通道间相位不一致性的校正。
实施例2:
本实施例提供一种毫米波多发多收线阵成像***,采用实施例1所述毫米波多发多收线阵的相位不一致性校正方法获取相位不一致性校正后的目标回波信号,校正后的目标回波信号直接用于目标的成像与重构。
综上所述,虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明,任何本领域普通技术人员,在不脱离本发明的精神和范围内,当可作各种更动与润饰,因此本发明的保护范围当视权利要求书界定的范围为准。
Claims (10)
1.毫米波多发多收线阵的相位不一致性校正方法,其特征在于,包括:
(2)利用金属平板作为参考目标,金属平板与毫米波多发多收线阵的扫描平面平行,录取金属平板的回波信号,将其作为相位不一致性的参考目标回波信号S ref ;
(3)利用电磁计算软件仿真得到与步骤(2)相同场景下的金属平板的仿真回波信号S fek ,其中相同场景是指仿真中的毫米波多发多收线阵的阵元分布、***的工作频率、金属平板的尺寸以及设置位置均与步骤(2)相同;
2.根据权利要求1所述的毫米波多发多收线阵的相位不一致性校正方法,其特征在于,所述毫米波多发多收线阵为四发八收的毫米波多发多收线阵。
3.根据权利要求1所述的毫米波多发多收线阵的相位不一致性校正方法,其特征在于,所述毫米波多发多收线阵顺序发射宽带调频连续波。
4.根据权利要求1、2或3所述的毫米波多发多收线阵的相位不一致性校正方法,其特征在于,步骤(2)中所设置的金属平板的尺寸需要满足毫米波多发多收线阵中的发射阵元的主瓣波束照射不到金属平板的边缘。
5.根据权利要求4所述的毫米波多发多收线阵的相位不一致性校正方法,其特征在于,步骤(3)中利用电磁计算软件FEKO仿真得到与步骤(2)相同场景下的金属平板的电磁散射回波信号。
6.根据权利要求3所述的毫米波多发多收线阵的相位不一致性校正方法,其特征在于,步骤(1)中,利用毫米波多发多收线阵录取目标回波信号,其中毫米波多发多收线阵中任意收发阵元组合的目标回波信号表达式为:
10.一种毫米波多发多收线阵成像***,其特征在于,采用权利要求1至9中任一权利要求所述毫米波多发多收线阵的相位不一致性校正方法获取相位不一致性校正后的目标回波信号用于目标的成像与重构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010840929.0A CN111708023B (zh) | 2020-08-20 | 2020-08-20 | 毫米波多发多收线阵的相位不一致性校正方法及成像*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010840929.0A CN111708023B (zh) | 2020-08-20 | 2020-08-20 | 毫米波多发多收线阵的相位不一致性校正方法及成像*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111708023A CN111708023A (zh) | 2020-09-25 |
CN111708023B true CN111708023B (zh) | 2020-11-06 |
Family
ID=72547303
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010840929.0A Active CN111708023B (zh) | 2020-08-20 | 2020-08-20 | 毫米波多发多收线阵的相位不一致性校正方法及成像*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111708023B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116840940B (zh) * | 2023-06-07 | 2024-04-05 | 珠海微度芯创科技有限责任公司 | 基于毫米波圆柱校准体算法***的通道校准方法及*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6352502B1 (en) * | 1998-12-03 | 2002-03-05 | Lightouch Medical, Inc. | Methods for obtaining enhanced spectroscopic information from living tissue, noninvasive assessment of skin condition and detection of skin abnormalities |
CN109478974B (zh) * | 2016-11-09 | 2022-03-29 | Lg 电子株式会社 | 在无线通信***中用于无线通信的方法及其装置 |
US10270455B2 (en) * | 2017-02-20 | 2019-04-23 | Qualcomm Incorporated | Multi-phase clock generation employing phase error detection in a controlled delay line |
CN107528643B (zh) * | 2017-09-04 | 2020-11-24 | 中国计量大学 | 一种多通道发射机及其通道一致性误差自动补偿方法 |
CN107809292A (zh) * | 2017-10-09 | 2018-03-16 | 武汉滨湖电子有限责任公司 | 一种发射通道多通道相位校正装置和校正方法 |
CN109283546A (zh) * | 2018-12-12 | 2019-01-29 | 唐之初 | 一种新型相位式激光测距装置及方法 |
CN110632595B (zh) * | 2019-09-24 | 2022-11-01 | 西南交通大学 | 一种主动毫米波成像方法及***、储存介质、成像设备 |
-
2020
- 2020-08-20 CN CN202010840929.0A patent/CN111708023B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111708023A (zh) | 2020-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10663563B2 (en) | On-site calibration of array antenna systems | |
CN109239684B (zh) | 基于电磁散射模型的雷达目标回波模拟***及其模拟方法 | |
CN106226761B (zh) | 一种高性能相干高频雷达多频探测方法 | |
CN110764068B (zh) | 一种多探头准远场电磁散射截面(rcs)外推测试*** | |
CN105974389B (zh) | 基于迭代处理的机扫米波雷达单脉冲测角方法 | |
CN108459312B (zh) | 基于复合多径因子估计的加权多频最大似然仰角估计方法 | |
CN106054165B (zh) | 一种实现分布式无源雷达目标检测的方法 | |
CN107229036B (zh) | 基于信号处理的多通道阵列雷达幅相误差在线检测方法 | |
CN109375213A (zh) | 基于子阵划分的频率分集阵列信号处理方法 | |
CN111708023B (zh) | 毫米波多发多收线阵的相位不一致性校正方法及成像*** | |
CN108152836A (zh) | 一种基于fda天线抗三星频差定位体制的欺骗方法 | |
CN115542243A (zh) | 基于阵列天线的干涉仪测向方法及*** | |
CN110109077B (zh) | 基于时间反演的mimo雷达相干角度估计方法 | |
CN109521418B (zh) | 基于干涉场的地基雷达测角方法 | |
CN113156437B (zh) | 一种高轨sar对低轨sar成像的射频干扰影响评估方法 | |
KR20150076755A (ko) | 위상배열 레이더의 배열면 정렬 방법 | |
CN112147593B (zh) | 一种高速密集***物破片目标的四维参数估计方法 | |
Xinghua et al. | Performance gain bounds of coherently combining multiple radars in a target-based calibration manner | |
CN110308445A (zh) | 一种基于车载数字阵调频连续波雷达的成像方法 | |
CN115469286A (zh) | 基于毫米波汽车雷达最小冗余mimo阵列的超分辨测角方法 | |
Hyder et al. | A new target localization method for bistatic FDA radar | |
CN110392387A (zh) | 无线信号的角度测量方法和设备 | |
CN111983575A (zh) | 一种有源和无源融合定标方法及装置 | |
Sharma et al. | Mutual coupling compensation in circular dipole array for FM based passive coherent location radar | |
Tan et al. | UWB-AOA estimation method based on a spare antenna array with virtual element |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |