CN112485765B - 一种电磁特性测试无源极化校准方法 - Google Patents

一种电磁特性测试无源极化校准方法 Download PDF

Info

Publication number
CN112485765B
CN112485765B CN202011221815.4A CN202011221815A CN112485765B CN 112485765 B CN112485765 B CN 112485765B CN 202011221815 A CN202011221815 A CN 202011221815A CN 112485765 B CN112485765 B CN 112485765B
Authority
CN
China
Prior art keywords
polarization
pss
calibration
measurement
scattering matrix
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
Application number
CN202011221815.4A
Other languages
English (en)
Other versions
CN112485765A (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.)
Beijing Electromechanical Engineering Research Institute
Original Assignee
Beijing Electromechanical Engineering Research Institute
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 Beijing Electromechanical Engineering Research Institute filed Critical Beijing Electromechanical Engineering Research Institute
Priority to CN202011221815.4A priority Critical patent/CN112485765B/zh
Publication of CN112485765A publication Critical patent/CN112485765A/zh
Application granted granted Critical
Publication of CN112485765B publication Critical patent/CN112485765B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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

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

本发明涉及一种电磁特性测试无源极化校准方法，属于雷达领域。本发明利用极化选择表面(Polarization Selective Surface,PSS)通过两种方案实现对目标的极化校准，PSS具有的极化散射特性兼具了金属二面角反射器和金属平板的极化散射特征，避免了测量过程中先后两次更换定标体而可能引入的摆放误差，从而保证了极化校准的精度，另一方面，作为定标体，PSS的复RCS因子只需一次测量便可获得，因此，PSS作为定标体用于极化校准简化了定标测量过程。根据PSS的极化散射矩阵特征，通过巧妙地设计PSS的摆放姿态，可获得简单矩阵形式，大大简化了极化校准的运算过程。PSS具有成本低廉的优点，且其在很宽频带内具有稳定的极化散射特性，因此，适用于室内宽带极化雷达测量系统的极化校准。

Description

一种电磁特性测试无源极化校准方法
技术领域
本发明属于雷达技术领域，具体涉及一种电磁特性测试无源极化校准方法。
背景技术
在雷达目标极化散射矩阵的测量中，由于受到雷达各极化通道间存在增益 失衡、交叉串扰及测量环境等因素影响，导致目标极化散射矩阵的测量值与其 理论值存在偏差，极化散射矩阵测量值与理论值的关系可表示为图1。
图1中，雷达发射通道和接收通道的传输特性分别用传输矩阵T和R表示(其 中，第一位下标表示接收极化，第二位下标表 示发出极化)，背景的极化散射矩阵用I表示，目标极化散射矩阵测量值M与 其真值S可由如下关系式表征
M＝R·(S+I)·T (1)
式中，M、S、I、R、T均为2×2矩阵。
上式可整理为
M＝R·S·T+Im (2)
式中，Im为背景极化散射矩阵的测量值。
极化校准的意义在于，从目标极化散射矩阵测量值中尽可能地还原出其真 实值，从而准确地获得被测目标的真实极化散射特性。由式(2)，可整理得到 如下表达式
S＝R-1·(M-Im)·T-1 (3)
可见，若已知雷达发射通道传输矩阵T、雷达接收通道传输矩阵R和背景极化散 射矩阵测量值Im，通过式(3)计算，便可实现对任意目标的极化校准。Im可通 过直接测量背景获得，因此，极化校准的关键是如何准确获得矩阵T和R。
求解矩阵T和R的通常做法为：通过设计一个或多个极化散射矩阵已知的定 标体，通过式(2)建立定标体极化散射矩阵测量值与其理论值关系，通过构建 多个方程组求得矩阵T和R。
专利申请201510097312.3中提出了一种可旋转双天线有源极化校准器 (Dual-antenna Polarimetric Active Radar Calibrators,DPARC)用于极化 校准，其原理框图如图2所示。接收天线和发射天线分别可绕各自电轴作0～360° 旋转，实现以不同极化接收和转发，从而获得多种不同形式的极化散射矩阵。 为提高天线极化纯度，申请人在收发天线口面处安装了极化滤波装置。
该专利中提出3中极化校准方案，其中与本发明相关的两种校准方案如下：
方案一的校准方案如下：
(1)转发天线与接收天线极化方式相互正交作同步旋转，此时，DPARC的 极化散射矩阵转发天线极化角θ分别取0°、90°、45° 三种姿态，可分别获得三种姿态下的极化散射矩阵/> α为DPARC的复RCS因子，可通过相对定标法测量获得；
(2)分别测量θ为0°、90°、45°状态下DPARC的极化散射矩阵M1、M2、 M3，在测量中，雷达信号经延时线及软件门处理后，目标区背景杂波Im被消除；
(3)将M1、M2、M3及S1、S2、S3通过式(3)建立等式，建立各极化分量测 量值与理论值的关系模型，通过联立8个方程组求解收雷达系统接受和发射通 道的传输矩阵R和T。
方案二的校准方案如下：
(1)转发天线与接收天线极化方式相同作同步旋转，此时，DPARC的极化 散射矩阵
(2)测量收发天线旋转状态下的极化散射矩阵M，在测量中，雷达信号经 延时线及软件门处理后，目标区背景杂波Im被消除。将M中各元素应用傅里叶 级数分别进行展开，有
其中wv表示所有的极化状态；
(3)将M及S通过式(3)建立等式，建立各极化分量测量值与理论值的 关系模型，提取式(4)常数项系数和2阶项系数，令对应项系数相等，求解收 雷达系统接收和发射通道的传输矩阵R和T。
该专利提出的可旋转双天线有源极化校准装置可提供多种极化散射特性， 据此可发展多种极化校准算法用于极化校准。但仍存在以下不足：
(1)可旋转双天线有源极化校准装置研发成本相对较高，因此，应用此种 方案进行极化校准测量的试验成本一般会比较高。
(2)单台DPARC工作带宽有限，通常适用于某一个波段的极化校准，无法 适用于超宽带极化校准。
另一文献“Fujita M.,Masuda T.,Satake M.SIR-C Polarimetric CalibrationExperiment Using Polarization Selective Dihedral[J].IEEE Transactions onGeoscience and Remote Sensing,Vol.36,No.6, pp:1872-1878,1998”提出了利用有源极化校准器(Polarimetric Active Radar Calibrators,PARC)和具有极化选择特性的二面角反射器(Polarization Selective Dihedrals，PSD)对SIR-C进行极化校准。该文献中提出极化校准 方案如下：
(1)将PSD分别绕其视线轴以0°、90°姿态放置，可分别获得两种姿态 下的极化散射矩阵γ为PSD的复RCS因子，可通过 测量获得。将PARC以某一种姿态放置，可获得第三组极化散射矩阵 />α为PSD的复RCS因子，可通过测量获得。
(2)分别测量0°、90°状态下PSD的极化散射矩阵M1、M2，测量某姿态 下的PARC的极化散射矩阵M3。
(3)将M1、M2、M3及S1、S2、S3通过式(3)建立等式，建立各极化分量测 量值与理论值的关系模型，通过联立方程组求解收雷达系统接收和发射通道的 传输矩阵R和T。
该文献提出的应用PARC和PSD的极化校准方案具有校准方法简单且有效的 优点，但仍存在以下不足：
(1)应用PARC进行极化校准测量的试验成本相对较高，且单台PARC的工 作带宽有限，不适用于超宽带校准；
(2)针对静态雷达系统的校准需求，在测量极化定标体PSM的过程中，需 要先后更换定标体，在此过程中很容易使两个定标体的位置摆放不一致，从而 造成系统的校准精度下降。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是如何提供一种电磁特性测试无源极化校准方法， 以将极化选择表面(Polarization Selective Surface，PSS)用于极化校准， 利用极化选择表面对不同极化电磁波的所呈现出其独特的极化散射特征，结合 背景技术中的专利申请中提出的校准算法，可实现低成本、高精度、超宽带的 极化校准工作。
(二)技术方案
为了解决上述技术问题，本发明提出一种电磁特性测试无源极化校准方法， 该方法包括如下步骤：
步骤S1：雷达开机测量金属球定标体，记录金属球的测量数据；
步骤S2：雷达开机测量背景，记录背景的极化散射矩阵Im；
步骤S3：将极化选择表面PSS以0°姿态放置于目标支架上，雷达开机测 量，记录此姿态下的极化散射矩阵测量值将PSS以90°姿态放置于目标支 架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值/>将PSS以45° 姿态放置于目标支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值 />
步骤S4：将各姿态下的极化散射矩阵测量值分别与背景的 极化散射矩阵Im作差运算，得到PSS的极化散射矩阵测量值M0、M90、M45，以 去除背景的影响；
步骤S5：利用金属球测量数据进行定标处理，得到PSS的复RCS因子β；
步骤S6：将被测目标放置于目标支架上，雷达开机测量，记录此姿态下的 极化散射矩阵测量值
步骤S7：将与背景的极化散射矩阵Im作差运算，得到被测目标的极化散 射矩阵测量值Mt，以去除背景的影响；
步骤S8：将PSS的极化散射矩阵测量值M0、M90、M45，应用式(16)～式 (25)可求得测量系统的及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参数；其中，各计 算公式如下：
将Mt和各参数值代入式(9)即可完成该目标的极化校准，得到目标的估计值 St，其中，公式(9)如下：
本发明还提供一种电磁特性测试无源极化校准方法，该方法包括如下步骤：
步骤-1：雷达开机测量金属球定标体，记录金属球的测量数据；
步骤-2：雷达开机测量背景，记录背景的极化散射矩阵Im；
步骤-3：使极化选择表面PSS绕雷达视线方向均匀慢速顺时针旋转N圈， N>＝1，雷达测量并存储PSS转过N圈的全部极化散射矩阵测量值PSM数据，记 为
步骤-4：将极化散射矩阵测量值与背景的极化散射矩阵Im作差运算，得 到PSS的极化散射矩阵测量值MR，去除背景的影响；
步骤-5：利用金属球测量数据进行定标处理，得到PSS的复RCS因子β；
步骤-6：将被测目标放置于目标支架上，雷达开机测量，记录此姿态下的 极化散射矩阵测量值
步骤-7：将与背景的极化散射矩阵Im作差运算，得到被测目标的极化散 射矩阵测量值Mt，去除背景的影响；
步骤-8：将MR各极化分量分别应用傅里叶级数进行展开，有
其中wv表示所有的极化状态，结合式(26)至式(29)，令对应项系数相 等，可求得测量系统的及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参 数；其中，结合式(26)至式(29)以及公式(9)如下：
将Mt和各参数值代入式(9)即可完成该目标的极化校准，得到目标的估计值 St，其中，公式(9)如下：
(三)有益效果
本发明提出一种电磁特性测试无源极化校准方法，本项发明中提出将PSS 作为新型定标体用于极化校准，存在以下3个优点：
(1)它所具有的极化散射特性兼具了金属二面角反射器和金属平板的极化 散射特征，为此，只需一个PSS定标体便可替代传统上必须应用二面角反射器 和金属平板组合才能实现的极化校准方案，避免了测量过程中先后两次更换定 标体而可能引入的摆放误差，从而保证了极化校准的精度。另一方面，作为定 标体，二面角反射器和金属平板的复RCS因子需要先后两次测量才能获得，而 PSS的复RCS因子只需一次测量便可获得，因此，PSS作为定标体用于极化校准 简化了定标测量过程。
(2)根据PSS的极化散射矩阵特征，通过巧妙地设计PSS的摆放姿态，可 获得等简单矩阵形式，大大简化了极化校准的运算过程。
(3)PSS具有成本低廉的优点，且其在很宽频带内具有稳定的极化散射特 性，因此，适用于室内宽带极化雷达测量系统的极化校准。
附图说明
图1为现有技术中极化散射矩阵测量误差示意图；
图2为现有技术中可旋转双天线有源转发器；
图3为电磁波入射PSS示意图；
图4为正入射电磁波极化与PSS位置关系示意图。
具体实施方式
为使本发明的目的、内容和优点更加清楚，下面结合附图和实施例，对本 发明的具体实施方式作进一步详细描述。
极化选择表面(Polarization Selective Surface，PSS)为由若干根平行 金属导线按照一定间隔均匀排布而成的周期性结构，如图3所示。图中灰色部 分为金属导线，要求其应具有良好的电导率。绿色平板为介质基板，要求其应 具有良好的透波特性。
当一束平面波照射到一无PSS上时(如图3所示)，随着入射电磁波极化 方式不同，PSS对入射电磁波表现为阻抗或容抗特征，从而实现其通态或阻态特 性。PSS的理想极化选择特性与雷达工作频率及其自身参数有关。
对于理想的PSS，当正入射PSS的电磁波的极化与金属导线平行时，PSS呈 现为阻态特性，此时电磁波能量全部被反射；当正入射PSS的电磁波的极化垂 直于金属导线时，PSS呈现通态特性，此时电磁波能量将全部通过，无电磁波能 量被反射。PSS对极化电磁波的响应特性适用于很宽的频带范围。
将PSS如图4(a)所示水平放置，使PSS金属导线与水平面平行，此时PSS 的极化散射矩阵可表示为β为PSS的复RCS因子，可通过测量获 得，为简化书写，后续推导过程省略β。以0°姿态的PSS为基准状态，使PSS 绕雷达视线方向顺时针旋转/>如图4(b)所示，此时，PSS的极化散射矩阵可 表示为
S＝A-1·S0·A (5)
式中，为旋转变换矩阵，对式(5)进行整理，有
可见，PSS的极化散射矩阵各分量是其旋转角的二阶三角函数，其极化散射 特性相当于具有相同复RCS因子的金属平板与金属二面角反射器叠加，因此， 只通过一个PSS便替代传统的金属二面角反射器、金属平板组合应用方案，且 无需更换定标体，避免先后两次定标测量，传统的极化校准方法可推广应用于 PSS。
由于背景对目标极化散射矩阵测量的影响可通过背景对消技术予以消除， 因此，为便于表示，后续均不考虑背景Im的影响，式(2)可表示为
对(7)式进行分解，有
式中，为雷达测量系统的交叉极化因子，表征了各极化通道间相互串扰的程度。
对式(8)两边作矩阵运算，可得任意被测目标的估计矩阵
可见，只要获得及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参 数，通过对(9)式进行运算，便可对任意被测目标实现极化校准。
基于上述极化校准原理，本发明提出将PSS作为定标体用于极化校准，并 通过巧妙地设计极化校准方案，大大简化了极化校准算法的求解过程。本发明 设计以下两种极化校准方案。
方案1：将PSS分别置于0°、90°、45°，各对应姿态下的极化散射矩阵 为
将式(10)、(11)、(12)分别代入式(8)进行展开,有
由式(13)、(14)、(15)中HH分量，有
联立(16)、(17)可得
将(18)代入(17)，可求出据此，再结合式(13)、(14)，可求出
至此，将上述各求解参数代入至式(9)便实现了对任意目标的测量极化散 射矩阵的极化校准。
方案2：具体校准过程如下。
使PSS绕雷达视线方向顺时针旋转，PSS的极化散射矩阵形式如式(6)所 示，将式(6)代入式(8)进行展开，省略PSS的复RCS因子β，有
将PSS的极化散射矩阵测量值M的各极化分量分别进行傅里叶级数展开， 并仅保留常数项和二阶展开项，有
按照(26)和(30)、式(27)和(31)、式(28)和(32)、式(29) 和(33)中各阶系数的对应关系，可建立若干等式关系求解及 RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参数。以HH极化分量为例，可建立如下关 系式
由式(34)、(35)、(36)可求得
RHH·THH＝c0,HH-a2,HH (37)
同理，对于VV极化分量，可分别求得
RVV·TVV＝c0,VV+a0,VV (40)
将式(38)、(39)、(41)、(42)分别代入式(27)、(28)，并与 式(31)、(32)中各阶项系数建立等式关系，便可求出RHH·TVV、RVV·THH。
至此，将上述各求解参数代入至式(9)便实现了对任意目标的测量极化散 射矩阵的极化校准。
除了上述列举出的两种极化校准方案以外，还可考虑应用PSS和二面角反 射器结合用于极化校准，0°、90°姿态的PSS可提供等简单形式 极化散射矩阵，45°姿态的二面角反射器可提供/>利用方案1的方法及过 程同样可完成极化校准，在此不作展开。
实施应用举例：
为进一步说明采用本项发明在极化校准中如何具体应用，现分别对采用方 案1和方案2的情况加以说明。
采用方案1的测量和极化校准处理步骤如下：
步骤-1：雷达开机测量金属球定标体，记录金属球的测量数据。
步骤-2：雷达开机测量背景，记录背景的极化散射矩阵Im。
步骤-3：将PSS以0°姿态放置于目标支架上，雷达开机测量，记录此姿态 下的极化散射矩阵测量值将PSS以90°姿态放置于目标支架上，雷达开机 测量，记录此姿态下的极化散射矩阵测量值/>将PSS以45°姿态放置于目标 支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值/>
步骤-4：将各姿态下的极化散射矩阵测量值分别与背景的极 化散射矩阵Im作差运算，得到PSS的极化散射矩阵测量值M0、M90、M45，去除 背景的影响。
步骤-5：利用金属球测量数据进行定标处理，得到PSS的复RCS因子β。
步骤-6：将被测目标放置于目标支架上，雷达开机测量，记录此姿态下的 极化散射矩阵测量值
步骤-7：将与背景的极化散射矩阵Im作差运算，得到被测目标的极化散 射矩阵测量值Mt，去除背景的影响。
步骤-8：将PSS的极化散射矩阵测量值M0、M90、M45，应用式(16)-式 (25)可求得测量系统的及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参数。其中，将Mt和 各参数代入式(9)即可完成该目标的极化校准，得到目标的估计值St。其中， />
采用方案2的测量和极化校准处理步骤如下：
步骤-1：雷达开机测量金属球定标体，记录金属球的测量数据。
步骤-2：雷达开机测量背景，记录背景的极化散射矩阵Im。
步骤-3：使PSS绕雷达视线方向均匀慢速顺时针旋转N圈(N应不小于1)， 雷达测量并存储PSS转过N圈的全部PSM数据，记为/>
步骤-4：将极化散射矩阵测量值与背景的极化散射矩阵Im作差运算，得 到PSS的极化散射矩阵测量值MR，去除背景的影响。
步骤-5：利用金属球测量数据进行定标处理，得到PSS的复RCS因子β。
步骤-6：将被测目标放置于目标支架上，雷达开机测量，记录此姿态下的 极化散射矩阵测量值
步骤-7：将与背景的极化散射矩阵Im作差运算，得到被测目标的极化散 射矩阵测量值Mt，去除背景的影响。
步骤-8：将MR各极化分量分别应用傅里叶级数进行展开，有
其中wv表示所有的极化状态，结合式(26)至式(29)，令对应项系数相 等，可求得测量系统的及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参 数。将Mt代入式(9)即可完成该目标的极化校准，得到目标的估计值St。
需要注意的是：
(1)除了采用本发明选择的0°、90°、45°摆放姿态外，可选择其他任 意三种不同PSS的摆放姿态；
(2)应用二面角反射器和金属平板组合进行极化校准的方案均适用于PSS。
本发明的技术关键点和欲保护点为：
(1)本发明提出的利用PSS进行极化校准的技术方案。
由于测量背景对目标极化散射矩阵测量的影响可通过背景对消技术予以消 除，因此，测量背景的影响可暂时不予考虑。目标极化散射矩阵的测量值与理 论值之间关系可表示为
对(A.1)式进行分解，有
式中，为雷达测量系统的交叉极化因子，表征了各极化通道间相互串扰的程度。
对式(A.2)两边作矩阵运算，可得任意被测目标的估计矩阵
可见，只要获得及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参 数，通过对(A.3)式进行运算，便可对任意被测目标实现极化校准。/>
基于上述极化校准原理，本发明提出将PSS作为定标体用于极化校准，并 通过巧妙地设计极化校准方案，大大简化了极化校准算法的求解过程。
记0°姿态的PSS对应为其金属导线与水平面平行时的姿态，以0°姿态PSS 为基准状态，使PSS绕雷达视线方向顺时针旋转此时，PSS的极化散射矩阵 可表示为
式中，β为PSS的复RCS因子，可通过测量获得，为简化书写，后续推导过 程省略β。由式(A.4)可见，PSS的极化散射矩阵各分量是其旋转角的二阶三 角函数，其极化散射特性相当于具有相同复RCS因子的金属平板与金属二面角 反射器叠加，因此，只通过一个PSS便替代传统的金属二面角反射器、金属平 板组合应用方案，且无需更换定标体，避免先后两次定标测量，传统的极化校 准方法可推广应用于PSS。
基于PSS极化散射特性特殊性，本发明设计以下两种极化校准方案。
方案1：将PSS分别置于0°、90°、45°，各对应姿态下的极化散射矩阵 为
将式(A.5)、(A.6)、(A.7)分别代入式(A.2)进行展开,有
由式(A.8)、(A.9)、(A.10)中HH分量，有
联立(A.11)、(A.12)可得
/>
将(A.13)代入式(A.12)，可求出
据此，再结合式(A.8)、(A.9)，可求出
至此，将上述各求解参数代入至式(A.3)便实现了对任意目标的测量极化 散射矩阵的极化校准。
方案2：参考前人提出将傅里叶级数用于极化校准的思想，该方法可推广应 用于本发明，具体校准过程如下。
使PSS绕雷达视线方向顺时针旋转，PSS的极化散射矩阵形式如式(A.4) 所示，将式(A.4)代入式(A.2)进行展开
将PSS的极化散射矩阵测量值M的各极化分量分别进行傅里叶级数展开， 并仅保留常数项和二阶展开项，有
按照(A.21)和(A.25)、式(A.22)和(A.26)、式(A.23)和(A.27)、 式(A.24)和(A.28)中各阶系数的对应关系，可建立若干等式关系求解 及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参数。以HH极化分量为例，可建 立如下关系式
由式(A.29)、(A.30)、(A.31)可求得
RHH·THH＝c0，HH-a2，HH (A.32)
同理，对于VV极化分量，可分别求得
RVV·TVV＝c0，VV+a0，VV (A.35)
将式(A.33)、(A.34)、(A.36)、(A.37)分别代入式(A.22)、(A.23)， 并与式(A.26)、(A.27)中各阶项系数建立等式关系，便可求出RHH·TVV、RVV·THH。
至此，将上述各求解参数代入至式(A.3)便实现了对任意目标的测量极化 散射矩阵的极化校准。
(2)应用PSS及PSS结合二面角反射器进行极化校准的有关方案。
应用PSS和二面角反射器结合用于极化校准，0°、90°姿态的PSS可提供 等简单形式极化散射矩阵，45°姿态的二面角反射器可提供/>采用与方案1的相同的方法及过程同样可完成极化校准，在此不作展开。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通 技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变 形，这些改进和变形也应视为本发明的保护范围。

Claims (7)

1.一种电磁特性测试无源极化校准方法，其特征在于，该方法包括如下步骤：
步骤S1：雷达开机测量金属球定标体，记录金属球的测量数据；
步骤S2：雷达开机测量背景，记录背景的极化散射矩阵Im；
步骤S3：将极化选择表面PSS以0°姿态放置于目标支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值将PSS以90°姿态放置于目标支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值/>将PSS以45°姿态放置于目标支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值/>
步骤S4：将各姿态下的极化散射矩阵测量值分别与背景的极化散射矩阵Im作差运算，得到PSS的极化散射矩阵测量值M0、M90、M45，以去除背景的影响；
步骤S5：利用金属球测量数据进行定标处理，得到PSS的复RCS因子β；
步骤S6：将被测目标放置于目标支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值
步骤S7：将与背景的极化散射矩阵Im作差运算，得到被测目标的极化散射矩阵测量值Mt，以去除背景的影响；
步骤S8：将PSS的极化散射矩阵测量值M0、M90、M45，应用式(16)～式(25)可求得测量系统的及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参数；其中，各计算公式如下：
将Mt和各参数值代入式(9)即可完成该目标的极化校准，得到目标的估计值St，其中，公式(9)如下：
2.如权利要求1所述的电磁特性测试无源极化校准方法，其特征在于，所述步骤S3中的0°、90°、45°三种摆放姿态可替换为其他任意三种不同PSS的摆放姿态。
3.如权利要求1所述的电磁特性测试无源极化校准方法，其特征在于，所述极化选择表面PSS为由若干根平行金属导线按照一定间隔均匀排布而成的周期性结构，其金属导线应具有良好的电导率，其介质基板应具有良好的透波特性。
4.如权利要求3所述的电磁特性测试无源极化校准方法，其特征在于，随着入射电磁波极化方式不同，PSS对入射电磁波表现为阻抗或容抗特征，从而实现其通态或阻态特性，PSS的理想极化选择特性与雷达工作频率及其自身参数有关；当入射电磁波极化与金属导线平行时，极化选择表面呈现为阻态特性，此时电磁波能量全部被反射；当入射电磁波极化垂直于金属导线时，极化选择表面呈现通态特性，此时电磁波能量将全部通过，无电磁波能量被反射。
5.一种电磁特性测试无源极化校准方法，其特征在于，该方法包括如下步骤：
步骤-1：雷达开机测量金属球定标体，记录金属球的测量数据；
步骤-2：雷达开机测量背景，记录背景的极化散射矩阵Im；
步骤-3：使极化选择表面PSS绕雷达视线方向均匀慢速顺时针旋转N圈，N>＝1，雷达测量并存储PSS转过N圈的全部极化散射矩阵测量值PSM数据，记为
步骤-4：将极化散射矩阵测量值与背景的极化散射矩阵Im作差运算，得到PSS的极化散射矩阵测量值MR，去除背景的影响；
步骤-5：利用金属球测量数据进行定标处理，得到PSS的复RCS因子β；
步骤-6：将被测目标放置于目标支架上，雷达开机测量，记录此姿态下的极化散射矩阵测量值
步骤-7：将与背景的极化散射矩阵Im作差运算，得到被测目标的极化散射矩阵测量值Mt，去除背景的影响；
步骤-8：将MR各极化分量分别应用傅里叶级数进行展开，有
其中wv表示所有的极化状态，结合式(26)至式(29)，令对应项系数相等，可求得测量系统的及RHH·THH、RHH·TVV、RVV·THH、RVV·TVV参数；其中，结合式(26)至式(29)以及公式(9)如下：
将Mt和各参数值代入式(9)即可完成该目标的极化校准，得到目标的估计值St，其中，公式(9)如下：
6.如权利要求5所述的电磁特性测试无源极化校准方法，其特征在于，所述极化选择表面PSS为由若干根平行金属导线按照一定间隔均匀排布而成的周期性结构，其金属导线应具有良好的电导率，其介质基板应具有良好的透波特性。
7.如权利要求5所述的电磁特性测试无源极化校准方法，其特征在于，随着入射电磁波极化方式不同，PSS对入射电磁波表现为阻抗或容抗特征，从而实现其通态或阻态特性，PSS的理想极化选择特性与雷达工作频率及其自身参数有关；当入射电磁波极化与金属导线平行时，极化选择表面呈现为阻态特性，此时电磁波能量全部被反射；当入射电磁波极化垂直于金属导线时，极化选择表面呈现通态特性，此时电磁波能量将全部通过，无电磁波能量被反射。
CN202011221815.4A 2020-11-05 2020-11-05 一种电磁特性测试无源极化校准方法 Active CN112485765B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011221815.4A CN112485765B (zh) 2020-11-05 2020-11-05 一种电磁特性测试无源极化校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011221815.4A CN112485765B (zh) 2020-11-05 2020-11-05 一种电磁特性测试无源极化校准方法

Publications (2)

Publication Number Publication Date
CN112485765A CN112485765A (zh) 2021-03-12
CN112485765B true CN112485765B (zh) 2023-11-03

Family

ID=74928276

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011221815.4A Active CN112485765B (zh) 2020-11-05 2020-11-05 一种电磁特性测试无源极化校准方法

Country Status (1)

Country Link
CN (1) CN112485765B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7345625B1 (en) * 2005-09-28 2008-03-18 Lockheed Martin Corporation Radar polarization calibration and correction
CN104614715A (zh) * 2015-03-05 2015-05-13 北京航空航天大学 一种可用于目标双站雷达散射截面测量定标与极化校准装置及其测量校准方法
CN104659498A (zh) * 2015-03-05 2015-05-27 北京航空航天大学 一种可旋转双天线有源极化校准装置及其极化校准方法
CN106154240A (zh) * 2016-06-20 2016-11-23 西北工业大学 一种极化散射矩阵快速测量系统及方法
CN108646226A (zh) * 2018-05-10 2018-10-12 北京航空航天大学 用于极化散射矩阵测量的无源极化校准体

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8416127B2 (en) * 2011-03-31 2013-04-09 Raytheon Company Dynamic calibration radar system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7345625B1 (en) * 2005-09-28 2008-03-18 Lockheed Martin Corporation Radar polarization calibration and correction
CN104614715A (zh) * 2015-03-05 2015-05-13 北京航空航天大学 一种可用于目标双站雷达散射截面测量定标与极化校准装置及其测量校准方法
CN104659498A (zh) * 2015-03-05 2015-05-27 北京航空航天大学 一种可旋转双天线有源极化校准装置及其极化校准方法
CN106154240A (zh) * 2016-06-20 2016-11-23 西北工业大学 一种极化散射矩阵快速测量系统及方法
CN108646226A (zh) * 2018-05-10 2018-10-12 北京航空航天大学 用于极化散射矩阵测量的无源极化校准体

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
全极化探地雷达采集系统及校准技术初探;李丽丽;《中国优秀硕士学位论文全文数据库基础科学辑》(第8期);全文 *

Also Published As

Publication number Publication date
CN112485765A (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
US10663563B2 (en) On-site calibration of array antenna systems
Van Zyl Calibration of polarimetric radar images using only image parameters and trihedral corner reflector responses
CN104659498B (zh) 一种可旋转双天线有源极化校准装置及其极化校准方法
CN104914408B (zh) 基于中国余数定理的频率、doa联合测量方法以及装置
CN108037374B (zh) 一种阵列天线近场标定方法
CN101149429A (zh) 适用于均匀圆阵的阵列互耦校正与信源测向方法
Fulton et al. Calibration of a digital phased array for polarimetric radar
CN108802669B (zh) 二维波达方向估计方法、二维波达方向估计装置及终端
CN110275139B (zh) 一种基于旋转式基元复用的超短基线定位系统及方法
CN111487478B (zh) 基于深度神经网络的角度依赖型复杂阵列误差校准方法
Jiang et al. Two new estimation algorithms for sensor gain and phase errors based on different data models
US20230155649A1 (en) Array Shape Reconstruction for Distributed Systems
CN113189592A (zh) 考虑幅相互耦误差的车载毫米波mimo雷达测角方法
CN113835063B (zh) 一种无人机阵列幅相误差与信号doa联合估计方法
CN112485765B (zh) 一种电磁特性测试无源极化校准方法
CN108646226A (zh) 用于极化散射矩阵测量的无源极化校准体
CN112556618B (zh) 安装误差测量方法、装置、电子设备和可读存储介质
CN108663577A (zh) 一种多探头球面近场交叉极化校正方法
CN108387868B (zh) 一种基于信号到达角度差测量的扫描式辐射源定位方法
CN116070082A (zh) 一种基于空间谱估计的极化敏感阵列幅相校正方法
CN113109807B (zh) 基于压缩感知的频率分集阵列雷达地下目标三维成像方法
Friedlander et al. Performance of direction-finding systems with sensor gain and phase uncertainties
JP3710409B2 (ja) 受信アレーアンテナキャリブレーション装置
Henault et al. Effects of mutual coupling on the accuracy of adcock direction finding systems
CN107515841A (zh) 计算电磁波波达角的方法

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