基于传递矩阵的近远场数据变换方法及装置制造方法及图纸

本发明专利技术涉及一种基于传递矩阵的近远场数据变换方法及装置，方法包括如下步骤：计算进行近远场数据变换所需要的最小合成孔径角Φ；对目标区进行网格化处理；构造近场观测矩阵A

【技术实现步骤摘要】
基于传递矩阵的近远场数据变换方法及装置


[0001]本专利技术涉及RCS测试与数据处理
，尤其涉及将近场条件下电磁散射测试数据外推至远场的数据变换方法及装置。

技术介绍

[0002]目标RCS(雷达散射截面积)的定义是一个远场的概念，要求照射目标的入射波和到达接收天线的散射波均为平面波。但真正理想的平面波是不存在的；在工程上，室外场RCS测试，一般采取距离远场的方式，即目标与测量系统之间的距离足够远，入射波和散射波近似为平面波；在室内则是利用抛物反射面等设备，在一定区域内，将馈源辐射的球面波转换成平面波。当测试条件受限，无法实现近似平面波时，研究人员采用近远场变换的测量方法，从近场散射回波中提取目标的远场散射信息，进而得到RCS。
[0003]现有的近远场变换方法主要包括：基于天线耦合公式的近远场变换方法(近场散射与远场RCS链条传递关系式)、平面波综合法、基于目标散射分布函数不变法、以及基于合成孔径的近远场变换方法。其中基于天线耦合公式的近远场变换方法需要目标的完全双站信息，工程实现起来较为困难；平面波综合法用一个天线单元先后或一组单元同时在各个单元位置上作加权发射，空间某点的场也依次或同时叠加，在一定的区域范围内球面波加权求和形成近似的平面波的方法，该方法目前主要用于天线测试，对于大型目标散射特性测试效率较低；基于目标散射分布函数不变的方法，利用ISAR成像基本原理，由近场测试值反演出目标散射分布函数，该方法缺点在于需要首先对目标进行成像再提取散射中心重构，流程较为复杂，同时对于点频数据也无法处理。
[0004]基于合成孔径的近远场变换方法CNFFFT(Circular Near
‑
Field to Far
‑
Field Transformation)，是对基于散射分布函数不变法的一种改进，该方法无需成像，利用汉克尔函数对相位项进行修正从而实现从目标近场散射特性至远场RCS的变换，是目前近场测试工程应用较广的一种方法。CNFFT方法的问题在于，该方法在推导时基于柱面波积分和展开时使用了Bessel函数的叠加定理与大宗量Hankel函数近似，因此主要应用于圆周均匀采样模式下的近远场变换，但实际近场测试多采用直线扫描模式，因此CNFFT方法难以直接适用。
[0005]针对以上不足，需要提供一种形式简单、易于实现且通用性强的改进的目标近场/远场数据转换方法，解决实际工程测试中目标近场散射特性变换至远场RCS的问题。

技术实现思路

[0006]本专利技术要解决的技术问题在于近远场变换方法形式复杂、通用性较差，针对现有技术中的缺陷，基于散射分布函数不变原理，提供一种形式简单、易于实现且通用性强的改进的目标近场/远场数据转换方法，解决实际工程测试中目标近场散射特性变换至远场RCS的问题。
[0007]为了解决上述技术问题，本专利技术提供了一种基于传递矩阵的近远场数据变换方
法，包括如下步骤：计算进行近远场数据变换所需要的最小合成孔径角Φ；对目标区进行网格化处理；构造近场观测矩阵A
N
，近场观测矩阵A
N
实现近场散射特性数据至目标散射中心的映射；构造远场重构矩阵A
F
，远场重构矩阵A
F
实现目标散射中心至远场散射特性数据的映射；基于近场观测矩阵A
N
与远场重构矩阵A
F
构造传递矩阵A；通过传递矩阵A实现目标的近场散射特性数据至远场散射特性数据的映射。
[0008]优选地，最小合成孔径角Φ的含义为，对于特定的方位角Φ0，近场测试的方位角度区域的最小值应该小于或等于Φ0
‑
Φ/2，近场测试的方位角度区域的最大值应该大于或等于Φ0+Φ/2。
[0009]优选地，计算进行近远场数据变换所需要的最小合成孔径角Φ的步骤中，最小合成孔径角Φ表示为：
[0010][0011]其中，R为测试距离、D为目标水平最大尺寸。
[0012]优选地，对目标区进行网格化处理的步骤中，网格化区域以目标区的中心为中心，构建横向长度、纵向长度均为D的区域，其中横向网格化点数Nx和纵向网格化点数Ny应满足：
[0013][0014]其中，Θ为近场散射特性数据相对于目标的孔径张角，f为测试数据最高频率，B为测试数据对应的带宽，c0表示真空光速。
[0015]优选地，构造近场观测矩阵A
N
的步骤中，近场观测矩阵A
N
表示为：
[0016][0017]其中，表示近场观测矩阵A
N
的第(m,n)个元素；R
m,n
表示以测试天线位置为坐标系原点，第n个测试状态下目标区第m个点距离测试天线的距离；x
m,n
与y
m,n
表示第m个点的坐标；f
n
是第n个测试状态对应的频率。
[0018]优选地，构造远场重构矩阵A
F
的步骤中，远场重构矩阵A
F
表示为：
[0019][0020]其中，表示远场重构矩阵A
F
的第(k,j)个元素；R
k,,j
表示以目标区中心为坐标系原点，第k个测试状态下目标区第j个点距离坐标系原点的距离，x
j
与y
j
表示第j个点的坐标；f
k
与θ
k
是第k测试状态对应的频率和角度。
[0021]优选地，基于近场观测矩阵A
N
与远场重构矩阵A
F
构造传递矩阵A的步骤中，传递矩
阵A表示为：
[0022]A＝A
F
·
A
N
。
[0023]优选地，通过传递矩阵A实现目标的近场散射特性数据至远场散射特性数据的映射的步骤中，S
F
＝A
·
S
N
；其中，S
N
为目标的近场散射特性数据，S
F
为目标的远场散射特性数据。
[0024]本专利技术还提供了一种基于传递矩阵的近远场数据变换装置，包括：第一计算模块，第一计算模块用于计算进行近远场数据变换所需要的最小合成孔径角Φ；网格处理模块，网格处理模块用于对目标区进行网格化处理；近场观测矩阵构造模块，近场观测矩阵构造模块用于构造近场观测矩阵A
N
，近场观测矩阵A
N
实现近场散射特性数据至目标散射中心的映射；远场重构矩阵构造模块，远场重构矩阵构造模块用于构造远场重构矩阵A
F
，远场重构矩阵A
F
实现目标散射中心至远场散射特性数据的映射；传递矩阵构造模块，传递矩阵构造模块用于基于近场观测矩阵A
N
与远场重构矩阵A
F
构造传递矩阵A；第二计算模块，第二计算模块用于通过传递矩阵A实现目标的近场散射特性数据至远场散射特性本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于传递矩阵的近远场数据变换方法，其特征在于，包括如下步骤：计算进行近远场数据变换所需要的最小合成孔径角Φ；对目标区进行网格化处理；构造近场观测矩阵A
N
，所述近场观测矩阵A
N
实现近场散射特性数据至目标散射中心的映射；构造远场重构矩阵A
F
，所述远场重构矩阵A
F
实现目标散射中心至远场散射特性数据的映射；基于所述近场观测矩阵A
N
与所述远场重构矩阵A
F
构造传递矩阵A；通过所述传递矩阵A实现目标的近场散射特性数据至远场散射特性数据的映射。2.根据权利要求1所述的基于传递矩阵的近远场数据变换方法，其特征在于：所述最小合成孔径角Φ的含义为，对于特定的方位角Φ0,近场测试的方位角度区域的最小值应该小于或等于Φ0
‑
Φ/2,近场测试的方位角度区域的最大值应该大于或等于Φ0+Φ/2。3.根据权利要求2所述的基于传递矩阵的近远场数据变换方法，其特征在于：所述计算进行近远场数据变换所需要的最小合成孔径角Φ的步骤中，所述最小合成孔径角Φ表示为：其中，R为测试距离、D为目标水平最大尺寸。4.根据权利要求1所述的基于传递矩阵的近远场数据变换方法，其特征在于：所述对目标区进行网格化处理的步骤中，网格化区域以所述目标区的中心为中心，构建横向长度、纵向长度均为D的区域，其中横向网格化点数Nx和纵向网格化点数Ny应满足：其中，Θ为近场散射特性数据相对于目标的孔径张角，f为测试数据最高频率，B为测试数据对应的带宽，c0表示真空光速。5.根据权利要求1所述的基于传递矩阵的近远场数据变换方法，其特征在于：所述构造近场观测矩阵A
N
的步骤中，所述近场观测矩阵A
N
表示为：其中，表示所述近场观测矩阵A
N
的第(m,n)个元素；R
m,n
表示以测试天线位置为坐标系原点，第n个测试状态下目标区第m个点距离测试天线的距离；x
m,n
与y
m,n
表示第m个点的坐标；f
n
是第n个测试状态对应的频率。6.根据权利要求1所述的基于传递矩阵的近远场数据变换方法，其特征在于：所述构造远场重构矩阵A
F
的步骤中，所述远场重构矩阵A
F
表示为：
其中，表示所述远场重构矩阵A
F
的第(k,j)个元素；R
k,,j
表示以目标区中心为坐标系原点，第k个测试状态下目标区第j个点距离坐标系原点的距离，x
j
与y
j
表示第j个点的坐标；f
k
与θ
k
是第k测试状态对应的频率和角度。7.根据权利要求1所述的基于传递矩阵的近远场数据变换方法，其特征在于...

【专利技术属性】
技术研发人员：吕鸣，
申请(专利权)人：北京环境特性研究所，
类型：发明
国别省市：