CN111257823A - 一种基于耦合阵列的测向方法 - Google Patents

一种基于耦合阵列的测向方法 Download PDF

Info

Publication number
CN111257823A
CN111257823A CN201911369799.0A CN201911369799A CN111257823A CN 111257823 A CN111257823 A CN 111257823A CN 201911369799 A CN201911369799 A CN 201911369799A CN 111257823 A CN111257823 A CN 111257823A
Authority
CN
China
Prior art keywords
signal source
array
array antenna
coupling
phase detection
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.)
Granted
Application number
CN201911369799.0A
Other languages
English (en)
Other versions
CN111257823B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201911369799.0A priority Critical patent/CN111257823B/zh
Publication of CN111257823A publication Critical patent/CN111257823A/zh
Application granted granted Critical
Publication of CN111257823B publication Critical patent/CN111257823B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction
    • G01S3/46Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
    • G01S3/48Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • 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
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/04Details
    • G01S3/12Means for determining sense of direction, e.g. by combining signals from directional antenna or goniometer search coil with those from non-directional antenna
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q21/00Antenna arrays or systems
    • H01Q21/06Arrays of individually energised antenna units similarly polarised and spaced apart
    • H01Q21/061Two dimensional planar arrays

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Multimedia (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种基于耦合阵列的测向方法。方法采用测向系统针对耦合阵列建立混合等效电路模型,通过CMOS图像传感器拍摄信号源的图像并处理后获得入射角信息,同时利用基于耦合阵列的信号相位检测系统实时检测信号源发出的电磁入射波被接收时的接收电流相位差信息;最后利用接收电流信息和入射角信息相结合通过优化算法求解混合等效电路模型,采用混合等效电路模型对待测信号源采集获得的接收电流相位差信息进行处理,获得信号源的位置,实现了信号源的定位测量。本发明能够实现基于耦合阵列的信号源的测向,测向方法可延伸应用到波束扫描、赋形和测向等多个领域中。

Description

一种基于耦合阵列的测向方法
技术领域
本发明公开了一种耦合阵列的应用方法,其中涉及了一种采用耦合阵列的 相关理论推导实现测向应用的方法。
背景技术
随着射频和基带集成电路的发展,基于阵列天线的微波系统小型化需求逐 渐增强。虽然射频前端和基带系统小型化会改善工作性能,但是系统中小型化 的阵列天线会导致阵列天线单元间产生强耦合,使得各射频通路电流中具有非 独立成分并难以分离,使得系统工作性能极大降低。
阵列天线单元间的耦合关系可利用一个电路阻抗有效表征,但是实际上很 难测量该电路阻抗,尤其是对于大尺寸耦合阵列或者与接收机模块相接的耦合 阵列。很少有一个有效的模型和方法可以描述耦合阵列的耦合关系以及测量耦 合矩阵,因此在未知电路阻抗的情况下利用耦合阵列进行无线应用是十分困难 的。
无线电定位测向,目的是为了估计非协作信号源发出的射频或微波信号的 入射角。传统的定位方法包括旋转或切换有向天线、Watson-Watt测向法、到达 时差法(TDOA)、和差分析、多普勒分析(Doppler)以及应用在更为复杂的 相控阵系统和多输入多输出(MIMO)系统中的多信号分类(MUSIC)、基于 旋转不变性的信号参数估计(ESPRIT)、最大似然估计(MLE)等波达方向角 估计(DOA)方法。随着射频、微波和毫米波集成电路的发展,基于阵列天线 的小型化测向系统设计需求越来越大,但是阵列天线的小型化会导致强耦合的 产生,导致传统的测向方法失效。
发明内容
为了有效解决上述问题,本发明提出了一种耦合阵列的测向应用方法,解 决了采用耦合阵列对非协作信号源进行无线微波定位的技术问题。
本发明针对具有强耦合特点的测向应用,本发明中耦合是指阵列天线单元 之间的接收电流会相互影响。耦合阵列通常是指多个间距小于二分之一波长的 阵列天线单元排布而成的阵列天线结构。本发明中强耦合通常是指耦合阵列上 的阵列天线单元接收端口回波损耗由于耦合的影响小于15dB的情况。
本发明采用的技术方案是:
方法采用测向系统,所述的测向系统包括一个CMOS图像传感器、一个基 于耦合阵列的信号相位检测系统、一个信号源和电脑,CMOS图像传感器朝向 信号源拍摄,信号源位于CMOS图像传感器视野范围内,CMOS图像传感器输 出连接到电脑;基于耦合阵列的信号相位检测系统包括耦合阵列、接收机模块、 相位检测模块和单片机;耦合阵列是由L个阵列天线单元(13-18)规则排布构 成,接收机模块是由L个接收机构成,相位检测模块(25-29)是由N-1个相位 检测芯片电路构成,相位检测芯片电路有2个输入端,阵列天线单元经接收机连接到其余L-1个相位检测芯片电路的输入端,其余L-1个阵列天线单元依次 经各自的接收机连接到相位检测芯片电路的另一个输入端,L-1个相位检测芯 片电路的输出连接到单片机的输入端,单片机输出端连接到电脑。
针对耦合阵列建立混合等效电路模型,通过CMOS图像传感器拍摄信号源 的图像并处理后获得入射角信息,同时利用基于耦合阵列的信号相位检测系统 实时检测信号源发出的电磁入射波被接收时的接收电流相位差信息;最后利用 接收电流信息和入射角信息相结合通过优化算法求解混合等效电路模型,采用 混合等效电路模型对待测信号源采集获得的接收电流相位差信息进行处理,获 得信号源的位置,实现了信号源的定位测量。
本发明中,每个阵列天线单元包括接收信号、耦合参数、接收电流的多方 面因素;耦合参数是指阵列天线单元之间存在的固有的耦合阻抗参数。
所述方法具体为:
1)建立混合等效电路模型
所述的混合等效电路模型中作为入射到阵列天线单元表面时信号源所发射 的电磁入射波等效的电压向量V,作为阵列天线单元的耦合参数等效的电路阻 抗矩阵Z以及作为阵列天线单元接收到信号源所发射电磁入射波的接收电流等 效的电路电流向量I,表示为:
V=(r+Z)I
上述公式两边均除以V中任意等效电压的模值|V|后表示为:
Figure BDA0002339370290000021
Figure BDA0002339370290000022
r=η/cosθ
其中,r为由r构成的对角矩阵,r为入射到耦合阵列表面时电磁入射波的 匹配参数,η为波阻抗,θ和
Figure BDA0002339370290000031
分别对应入射角中的俯仰角和方位角;
2)将信号源布置在一个空间位置上,CMOS图像传感器拍摄信号源的图像 并进行图像分析处理,获得信号源相对于CMOS图像传感器图像成像面的视角 坐标
Figure BDA0002339370290000032
3)将信号源布置在一个空间位置上进行测量,通过基于耦合阵列的信号相位 检测系统同步获得相位差信息,信号源发出的电磁入射波被阵列天线单元接收 时的接收电流相位为β12,...,βL,L表示基于耦合阵列的信号相位检测系统中 的阵列天线单元的总数,采用基于耦合阵列的信号相位检测系统实时检测接收 电流间的相位差信息,以第1个阵列天线单元为参考时,其余L-1个阵列天线 单元上的接收电流相位均与第1个阵列天线单元上的接收电流相位进行相减归 一化,得到的相位差分别为α1=0,α2=β21,...,αi=βi1,...,αL=βL1,可作 为一组接收电流相位差信息。
4)重复步骤2)和3)多次,将信号源布置在不同空间位置上分别进行测量, 获得视角坐标
Figure BDA0002339370290000033
及其对应的接收电流信息,将视角坐标
Figure BDA0002339370290000034
通过三维空间 坐标变换计算得到电磁入射波的入射角
Figure BDA0002339370290000035
θ和
Figure BDA0002339370290000036
分别对应入射角中的俯仰 角和方位角,根据阵列天线远场计算公式和阵列天线单元具体排布方式构造成 为正向等式的复数形式输出Vn,将接收电流相位差信息进行构造,成为正向等 式的复数形式输入I'n,在得到多组Vn和I'n后,输入到以下公式的正向等式中, 采用优化算法求解获得正向等式中的电路阻抗、系数矩阵、权值矩阵和偏置常 数;
Vn=(r+Z)In=(r+Z)W[f(σnI'n+bn)]
其中,σn表示第n组系数矩阵,σn=[σn1n2,...,σnL],σn1为第n组系数矩阵 中的第1个元素,σn2为第n组系数矩阵中的第2个元素,σnL为第n组系数矩阵 中的第L个元素;W为权值矩阵,bn为第n个偏置常数,f(·)为非线性函数矢量;
Figure BDA0002339370290000037
归一化接收电流向量为
Figure BDA0002339370290000038
A1n为第1个阵列天线单元上的接收电流归一化幅度,A2n为第2个阵列单元上的 接收电流归一化幅度,ALn为第L个阵列天线单元上的接收电流归一化幅度;Z 表示电路阻抗矩阵。
本发明中,采用复数形式的批梯度下降法对(r+Z)W、σn以及bn进行优化, 当均方误差最小时即可实现求解。
本发明中,阵列天线单元的数目为L时,需要优化的电路阻抗、系数矩阵、 权值和偏置常数的数目为L(3L+1)个。
5)对于待测信号源,重复实施步骤3)获得接收电流相位差信息,将接收电流 相位差信息进行构造得到I'n,I'n输入正向等式后在0°至90°范围内以固定步长扫 描θ输出多个Vn,对每个Vn进行计算处理得到多个入射角,将得到的多个入射 角中的θ与正向等式计算过程中对应的扫描角进行匹配,入射角中的θ与扫描 角相等时的入射角为待测信号源的位置。
所述的优化算法采用包括但不限于全复数反向传播(FCBP)方法、蒙特卡 洛算法。
所述的应用方法均可以用硬件描述语言写入FPGA或CPLD等可编程逻辑 器件中实现波束扫描、赋形和测向等正问题应用。
所述的CMOS图像传感器和基于耦合阵列的信号相位检测系统相位位置保 持固定。
所述的CMOS图像传感器与电脑通过USB端口直接相连,单片机通过串口 转USB接口与电脑直接相连。
所述的信号源处于阵列天线远场区域范围内。
本发明建立了一种耦合阵列的混合等效电路模型,不同于传统的等效电路 模型,该混合等效电路模型包括对阵列天线单元接收到信号源所发射电磁入射 波的接收信号、耦合参数和接收电流的建模。在考虑各种因素影响下,根据该 混合等效电路模型并求解了关于阵列天线单元间耦合参数的正向等式,基于该 正向等式可实现耦合阵列的测向。
本发明的测向方法可延伸应用到波束扫描、赋形和测向等多个领域中。
在多径效应较小的不同场合下,已知混合电路等效模型中的电压向量Vn和 归一化接收电流向量I'n组成信息,经过方法求解优化均得到模型。
本发明分析了耦合阵列在工作时的物理过程并解决了现存采用耦合阵列进 行测向的问题,所具有的有益效果是:
1.该方法可使得耦合阵列在波束扫描、赋形及测向等领域进行应用。
2.相同频率下,实现了使用尺寸更小的耦合阵列进行测向。对于由6个阵列 天线单元组成的耦合阵列,每次进行测向的时间仅需50ms。
附图说明
附图1是混合等效电路模型的示意图。
图中:1、2是电磁入射波,3、4是阵列天线单元。
附图2是实施例1中采用的阵列天线单元以及由其构成的耦合阵列。
图中:(a)是阵列天线单元结构,(b)是由阵列天线单元构成的耦合阵列。5 表示接收端口,6、7、8、9、10、11表示阵列天线单元。
附图3是实施例1中FCBP算法优化的收敛曲线图。
附图4是实施例1中的测向结果计算值和实际值对比。
附图5是基于耦合阵列的测向软硬件系统及模块示意图。
图中:12表示信号源,13、14、15、16、17、18表示阵列天线单元,19、 20、21、22、23、24表示接收机模块,25、26、27、28、29表示相位检测模块, 30表示单片机,31表示电脑,32表示CMOS图像摄像头
附图6是实施例2中采用的倒F圆极化耦合阵列结构。
图中:33、34、35、36、37、38表示倒F圆极化阵列天线单元。
附图7是实施例2中FCBP算法优化的收敛曲线图。
附图8是实施例2中的测向结果计算值和实际值对比图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
如图5所示,方法采用测向系统,测向系统包括一个CMOS图像传感器32、 一个基于耦合阵列的信号相位检测系统、一个信号源12和电脑31,CMOS图像 传感器32朝向信号源12拍摄,信号源12位于CMOS图像传感器32视野范围 内,同时信号源12处于阵列天线远场区域范围内。CMOS图像传感器32输出 连接到电脑31;基于耦合阵列的信号相位检测系统包括耦合阵列、接收机模块 19-24、相位检测模块25-29和单片机30;耦合阵列是由L个阵列天线单元13-18 规则排布构成,接收机模块19-24是由L个接收机构成,相位检测模块25-29是 由N-1个相位检测芯片电路构成,相位检测芯片电路有2个输入端,经由2个 输入端输入到电路中的信号可以经过鉴相得到相位差。阵列天线单元13-18和接 收机的数量相同,相位检测芯片电路的数量比阵列天线单元13-18和接收机的数 量少一个,阵列天线单元17经接收机连接到其余L-1个相位检测芯片电路的输 入端,其余L-1个阵列天线单元13-16、18依次经各自的接收机连接到相位检 测芯片电路的另一个输入端,L-1个相位检测芯片电路的输出连接到单片机30 的输入端,单片机30输出端连接到电脑31。
CMOS图像传感器32和基于耦合阵列的信号相位检测系统中的耦合阵列相 位位置保持固定。CMOS图像传感器32与电脑31通过USB端口直接相连,单 片机30通过串口转USB接口与电脑31直接相连。
信号源12发出电磁波信号,被信号相位检测系统的耦合阵列的阵列天线单 元13-18接收,进而经接收机模块19-24的L个接收机和L-1个相位检测模块 25-29的相位检测仪器后将信号发送到单片机30,经单片机30发送到电脑31 中进行分析处理。
针对耦合阵列建立混合等效电路模型,通过CMOS图像传感器32拍摄信号 源12的图像并处理后获得入射角信息,同时利用基于耦合阵列的信号相位检测 系统实时检测信号源12发出的电磁入射波被接收时的接收电流相位差信息;最 后利用接收电流信息和入射角信息相结合通过优化算法求解混合等效电路模型, 进而用接收电流信息和入射角信息对信号源12所需的系统内部进行校准,采用 混合等效电路模型对待测信号源采集获得的接收电流相位差信息进行处理,获 得信号源的位置,实现了信号源的定位测量。
图1所示为两个阵列天线单元3、4接收到电磁入射波1、2时的混合等效 电路模型。作为入射到阵列天线单元表面时信号源8所发射的电磁入射波等效 的电压向量V=[V1=|V|,V2=|V|e-jkdsinθ]T,作为阵列天线单元的耦合参数等效的 电路阻抗矩阵:
Figure BDA0002339370290000061
以及作为阵列天线单元接收到信号源8所发射电磁入射波的接收电流等效 的电路电流向量
Figure BDA0002339370290000062
基于V=(r+Z)I,可以得到Vn=[V1= 1,V2=e-jkdsinθ]T,
Figure BDA0002339370290000063
Figure BDA0002339370290000064
k为波矢,d为两个阵列天线单元的间距。
本发明方法实施过程可以为:
1)建立混合等效电路模型
混合等效电路模型中作为入射到阵列天线单元表面时信号源12所发射的电 磁入射波等效的电压向量V,作为阵列天线单元的耦合参数等效的电路阻抗矩 阵Z以及作为阵列天线单元接收到信号源12所发射电磁入射波的接收电流等效 的电路电流向量I,表示为:
V=(r+Z)I
上述公式两边均除以V中任意等效电压的模值|V|后表示为:
Figure BDA0002339370290000065
Figure BDA0002339370290000066
r=η/cosθ
其中,r为由r构成的对角矩阵,r表征入射到耦合阵列表面时电磁入射波 的匹配情况,η为波阻抗,θ和
Figure BDA0002339370290000067
分别对应入射角中的俯仰角和方位角;
2)将信号源12布置在一个空间位置上,CMOS图像传感器32拍摄信号源 12的图像并进行图像分析处理,获得信号源12相对于CMOS图像传感器32图 像成像面的视角坐标
Figure BDA0002339370290000071
视角坐标
Figure BDA0002339370290000072
是在信号源12在视角坐标系下的相对于CMOS图像传感器 32光心的图像位置坐标,视角坐标系是以CMOS图像传感器32的光心为原点、 以过成像面中心的水平方向作为x轴,以过成像面中心的竖直方向作为y轴、以 CMOS图像传感器32的光心和成像面中心连线作为z轴的三维坐标,成像面即 为图像两者重合,θc为信号源12在图像中位置向x轴作投影的投影点到CMOS 图像传感器32光心的连线与成像面中心到CMOS图像传感器32光心的连线之 间的夹角,
Figure BDA0002339370290000073
为信号源12在图像中位置向y轴作投影的投影点到CMOS图像传 感器32光心的连线与成像面中心到CMOS图像传感器32光心的连线之间的夹 角。
3)将信号源12布置在一个空间位置上进行测量,通过基于耦合阵列的信号 相位检测系统同步获得相位差信息,信号源12发出的电磁入射波被阵列天线单 元接收时的接收电流相位为β12,...,βL,L表示基于耦合阵列的信号相位检测 系统中的阵列天线单元的总数,信号源12在每一个位置上时,采用基于耦合阵 列的信号相位检测系统实时检测接收电流间的相位差信息,以第1个阵列天线 单元为参考时,其余L-1个阵列天线单元上的接收电流相位均与第1个阵列天 线单元上的接收电流相位进行相减归一化,得到的相位差分别为α1=0, α2=β21,...,αi=βi1,...,αL=βL1,可作为一组接收电流相位差信息;
4)重复步骤2)和3)多次,将信号源12布置在不同空间位置上分别进行测量, 获得视角坐标
Figure BDA0002339370290000074
及其对应的接收电流信息,将这样已知的视角坐标
Figure BDA0002339370290000075
通 过三维空间坐标变换计算得到电磁入射波的入射角
Figure BDA0002339370290000076
θ和
Figure BDA0002339370290000077
分别对应入射 角中的俯仰角和方位角,根据阵列天线远场计算公式和阵列天线单元具体排布 方式构造成为正向等式的复数形式输出Vn,阵列天线单元的具体排布方式指的 是阵列天线单元在空间中的不同分布方式,根据阵列天线单元在空间中的相对 位置关系,求得电磁波信号入射到阵列天线单元上时个阵列天线单元间的相对 相位值,输出Vn被构造为以e为底,幂数为j乘以相对相位值的指数。j为虚数 单位,e为自然常数2.71828。然后将这样已知的接收电流相位差信息进行构造, 即构造为以e为底,幂数为j乘以接收电流相位差的指数。j为虚数单位,e为 自然常数2.71828;成为正向等式的复数形式输入I'n,在得到多组Vn和I'n后,输 入到以下公式的正向等式中,采用优化算法求解获得正向等式中的电路阻抗(即 耦合参数)、系数矩阵、权值矩阵和偏置常数;
Vn=(r+Z)In=(r+Z)W[f(σnI′n+bn)]
其中,σn表示第n组系数矩阵,σn=[σn1n2,...,σnL],σn1为第n组系数矩阵 中的第1个元素,σn2为第n组系数矩阵中的第2个元素,σnL为第n组系数矩阵 中的第L个元素;W为权值矩阵,bn为第n个偏置常数,f(·)为非线性函数矢量;
Figure BDA0002339370290000081
归一化接收电流向量为
Figure BDA0002339370290000082
A1n为第1个阵列天线单元上的接收电流归一化幅度,A2n为第2个阵列单元上的 接收电流归一化幅度,ALn为第L个阵列天线单元上的接收电流归一化幅度;Z 表示电路阻抗矩阵;
5)确定正向等式后,对于待测信号源,重复实施步骤3)获得接收电流相位差 信息,将接收电流相位差信息进行构造得到I'n,I'n输入正向等式后在0°至90°范 围内以固定步长扫描θ输出多个Vn,对每个Vn进行计算处理得到多个入射角, 将得到的多个入射角中的θ与正向等式计算过程中对应的扫描角进行匹配,入 射角中的θ与扫描角相等时的入射角为待测信号源的位置。是指在一定范围内, 如0至90°的区间内按照一定步长逐个代入到正向等式中进行计算,求得的角度 位置θ与此时用于计算的扫描角θ进行比较,如果相等此时求得的角度位置即 为待测信号源的位置。
本发明的实施例如下:
实施例1:
实施例1采用图2所示的阵列天线单元以及由其构成的耦合阵列来证明上 述方法。阵列天线单元间距为0.3倍的工作波长,基板介电常数为6.15,基板厚 度为1.6mm。
首先根据实际情况自行设定电磁入射波的入射角
Figure BDA0002339370290000083
取值范围,将已知的 入射角进行构造可以得到Vn,通过CSTTM Microwave Studio对图2(b)所示的耦 合阵列进行仿真得到电路阻抗矩阵Z,基于仿真得到的电路阻抗矩阵Z、入射角 和正向等式可以计算得到阵列天线单元6、7、8、9、10、11对应的In,然后以 阵列天线单元6的接收电流为基准进行归一化,可以得到阵列天线单元6、7、8、 9、10、11的归一化接收电流向量I'n。Vn与I'n可以共同构成用于FCBP算法优化 的先验信息。
FCBP算法优化的收敛曲线如图3所示,FCBP算法在经过5000余次的迭代 后达到设定的均方误差MSE阈值。
图4将测向结果的计算值和实际值进行了比较,θ的平均绝对误差为0.17°, 均方根误差为0.71°;
Figure BDA0002339370290000091
的平均绝对误差为0.24°,均方根误差为1.04°。结果验 证了本发明中提出的耦合阵列应用方法的有效性。
实施例2:
实施例2采用图5所示的测向软硬件系统对所述方法进行证明。
图5所示,硬件系统包括由阵列天线单元13、14、15、16、17、18构成的 耦合阵列,接收机模块19、20、21、22、23、24,相位检测模块25、26、27、 28、29、单片机30、电脑31以及CMOS图像摄像头32组成。信号源12发出 的电磁入射波被阵列天线单元13、14、15、16、17、18接收后,在阵列天线单 元上的接收电流通过接收端口输入到接收机模块19、20、21、22、23、24中。 以接收机模块23输出的中频信号为参考信号,接收机模块19、20、21、22、24 输出的中频信号输入到相位检测模块后,通过鉴相得到接收机模块19、20、21、 22、24与接收机模块23输出的中频信号之间的信号相位差信息,相位检测模块 的输出信号输入到单片机30中经过AD采样后输出,然后利用串口通信的方式 输入到电脑31中,在电脑中可收到相位差的具体数值。CMOS图像摄像头32 与电脑31连接可以显示信号源12的图像位置,采用坐标变换等方式可以计算 得到非协作信号源在耦合阵列坐标系下的角度位置,所得角度位置也是信号源 发出的电磁入射波的入射角。
图6所示为实施例中所采用的倒F圆极化耦合阵列结构,阵列天线单元间 距为0.3倍波长,其中阵列天线单元33、34、35、36、37、38分别对应图5中 的阵列天线单元13、14、15、16、17、18。
首先根据实际情况自行设定信号源12发出的电磁入射波的入射角范围,将 已知的入射角进行构造可以得到Vn。将信号源放置在CMOS图像摄像头视野范 围内的多个位置,可以在电脑上测得相位检测模块25、26、27、28、29输出的 5组相位差α23456,然后基于5组相位差可以构造得到接收到的归一 化电流向量I'n。Vn与I'n可以共同构成用于FCBP算法优化的先验信息。
FCBP算法优化的收敛曲线如图7所示,FCBP算法在经过17000余次的迭 代后达到设定的均方误差MSE阈值。
图8将测向结果的计算值和实际值进行了比较,θ的平均绝对误差为0.86°, 均方根误差为1.00°;
Figure BDA0002339370290000092
的平均绝对误差为1.11°,均方根误差为1.26°。结果验 证了本发明中提出的耦合阵列应用方法的有效性。

Claims (6)

1.一种基于耦合阵列的测向方法,其特征在于:方法采用测向系统,所述的测向系统包括一个CMOS图像传感器(32)、一个基于耦合阵列的信号相位检测系统、一个信号源(12)和电脑(31),CMOS图像传感器(32)朝向信号源(12)拍摄,信号源(12)位于CMOS图像传感器(32)视野范围内,CMOS图像传感器(32)输出连接到电脑(31);基于耦合阵列的信号相位检测系统包括耦合阵列、接收机模块(19-24)、相位检测模块(25-29)和单片机(30);耦合阵列是由L个阵列天线单元(13-18)规则排布构成,接收机模块(19-24)是由L个接收机构成,相位检测模块(25-29)是由N-1个相位检测芯片电路构成,相位检测芯片电路有2个输入端,阵列天线单元(17)经接收机连接到其余L-1个相位检测芯片电路的输入端,其余L-1个阵列天线单元(13-16、18)依次经各自的接收机连接到相位检测芯片电路的另一个输入端,L-1个相位检测芯片电路的输出连接到单片机(30)的输入端,单片机(30)输出端连接到电脑(31);针对耦合阵列建立混合等效电路模型,通过CMOS图像传感器(32)拍摄信号源(12)的图像并处理后获得入射角信息,同时利用基于耦合阵列的信号相位检测系统实时检测信号源(12)发出的电磁入射波被接收时的接收电流相位差信息;最后利用接收电流信息和入射角信息相结合通过优化算法求解混合等效电路模型,采用混合等效电路模型对待测信号源采集获得的接收电流相位差信息进行处理,获得信号源的位置,实现了信号源的定位测量。
2.根据权利要求1所述的一种基于耦合阵列的测向方法,其特征在于:所述方法具体为:
1)建立混合等效电路模型
所述的混合等效电路模型中作为入射到阵列天线单元表面时信号源(12)所发射的电磁入射波等效的电压向量V,作为阵列天线单元的耦合参数等效的电路阻抗矩阵Z以及作为阵列天线单元接收到信号源(12)所发射电磁入射波的接收电流等效的电路电流向量I,表示为:
V=(r+Z)I
上述公式两边均除以V中任意等效电压的模值|V|后表示为:
Figure FDA0002339370280000011
Figure FDA0002339370280000021
r=η/cosθ
其中,r为由r构成的对角矩阵,r为入射到耦合阵列表面时电磁入射波的匹配参数,η为波阻抗,θ和
Figure FDA0002339370280000024
分别对应入射角中的俯仰角和方位角;
2)将信号源(12)布置在一个空间位置上,CMOS图像传感器(32)拍摄信号源(12)的图像并进行图像分析处理,获得信号源(12)相对于CMOS图像传感器(32)图像成像面的视角坐标
Figure FDA0002339370280000027
3)将信号源(12)布置在一个空间位置上进行测量,通过基于耦合阵列的信号相位检测系统同步获得相位差信息,信号源(12)发出的电磁入射波被阵列天线单元接收时的接收电流相位为β12,...,βL,L表示基于耦合阵列的信号相位检测系统中的阵列天线单元的总数,采用基于耦合阵列的信号相位检测系统实时检测接收电流间的相位差信息,以第1个阵列天线单元为参考时,其余L-1个阵列天线单元上的接收电流相位均与第1个阵列天线单元上的接收电流相位进行相减归一化,得到的相位差分别为α1=0,α2=β21,...,αi=βi1,...,αL=βL1,可作为一组接收电流相位差信息;
4)重复步骤2)和3)多次,将信号源(12)布置在不同空间位置上分别进行测量,获得视角坐标
Figure FDA0002339370280000028
及其对应的接收电流信息,将视角坐标
Figure FDA0002339370280000029
通过三维空间坐标变换计算得到电磁入射波的入射角
Figure FDA0002339370280000026
θ和
Figure FDA0002339370280000025
分别对应入射角中的俯仰角和方位角,根据阵列天线远场计算公式和阵列天线单元具体排布方式构造成为正向等式的复数形式输出Vn,将接收电流相位差信息进行构造,成为正向等式的复数形式输入I'n,在得到多组Vn和I'n后,输入到以下公式的正向等式中,采用优化算法求解获得正向等式中的电路阻抗、系数矩阵、权值矩阵和偏置常数;
Vn=(r+Z)In=(r+Z)W[f(σnI'n+bn)]
其中,σn表示第n组系数矩阵,σn=[σn1n2,...,σnL],σn1为第n组系数矩阵中的第1个元素,σn2为第n组系数矩阵中的第2个元素,σnL为第n组系数矩阵中的第L个元素;W为权值矩阵,bn为第n个偏置常数,f(·)为非线性函数矢量;
Figure FDA0002339370280000022
归一化接收电流向量为
Figure FDA0002339370280000023
A1n为第1个阵列天线单元上的接收电流归一化幅度,A2n为第2个阵列单元上的接收电流归一化幅度,ALn为第L个阵列天线单元上的接收电流归一化幅度;Z表示电路阻抗矩阵;
5)对于待测信号源,重复实施步骤3)获得接收电流相位差信息,将接收电流相位差信息进行构造得到I'n,I'n输入正向等式后在0°至90°范围内以固定步长扫描θ输出多个Vn,对每个Vn进行计算处理得到多个入射角,将得到的多个入射角中的θ与正向等式计算过程中对应的扫描角进行匹配,入射角中的θ与扫描角相等时的入射角为待测信号源的位置。
3.根据权利要求2所述的一种基于耦合阵列的测向方法,其特征在于:
所述的优化算法采用包括但不限于全复数反向传播(FCBP)方法、蒙特卡洛算法。
4.根据权利要求1所述的一种基于耦合阵列的测向方法,其特征在于:
所述的CMOS图像传感器(32)和基于耦合阵列的信号相位检测系统相位位置保持固定。
5.根据权利要求1所述的一种基于耦合阵列的测向方法,其特征在于:
所述的CMOS图像传感器(32)与电脑(31)通过USB端口直接相连,单片机(30)通过串口转USB接口与电脑(31)直接相连。
6.根据权利要求1所述的一种基于耦合阵列的测向方法,其特征在于:
所述的信号源(12)处于阵列天线远场区域范围内。
CN201911369799.0A 2019-12-26 2019-12-26 一种基于耦合阵列的测向方法 Active CN111257823B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911369799.0A CN111257823B (zh) 2019-12-26 2019-12-26 一种基于耦合阵列的测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911369799.0A CN111257823B (zh) 2019-12-26 2019-12-26 一种基于耦合阵列的测向方法

Publications (2)

Publication Number Publication Date
CN111257823A true CN111257823A (zh) 2020-06-09
CN111257823B CN111257823B (zh) 2022-03-22

Family

ID=70945316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911369799.0A Active CN111257823B (zh) 2019-12-26 2019-12-26 一种基于耦合阵列的测向方法

Country Status (1)

Country Link
CN (1) CN111257823B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111914506A (zh) * 2020-08-05 2020-11-10 中国电子科技集团公司第三十六研究所 天线阵测向方法、装置及电子设备
CN113468844A (zh) * 2021-06-17 2021-10-01 浙江大学 一种耦合阵列波束综合的解析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009168498A (ja) * 2008-01-11 2009-07-30 Sony Corp 方向検出システム並びに方向検出装置
CN103633442A (zh) * 2012-08-24 2014-03-12 香港城市大学 相控阵列、相干源阵列、天线阵列以及用于控制它们的系统
CN107085198A (zh) * 2017-06-23 2017-08-22 中国电子科技集团公司第三十六研究所 一种构建四阵元立体阵列的方法和装置
CN109932683A (zh) * 2019-01-24 2019-06-25 中国电子科技集团公司第二十九研究所 一种基于仿生技术的小尺度阵列测向方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009168498A (ja) * 2008-01-11 2009-07-30 Sony Corp 方向検出システム並びに方向検出装置
CN103633442A (zh) * 2012-08-24 2014-03-12 香港城市大学 相控阵列、相干源阵列、天线阵列以及用于控制它们的系统
CN107085198A (zh) * 2017-06-23 2017-08-22 中国电子科技集团公司第三十六研究所 一种构建四阵元立体阵列的方法和装置
CN109932683A (zh) * 2019-01-24 2019-06-25 中国电子科技集团公司第二十九研究所 一种基于仿生技术的小尺度阵列测向方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王丽轩: "基于接收信号强度的可见光室内定位算法研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111914506A (zh) * 2020-08-05 2020-11-10 中国电子科技集团公司第三十六研究所 天线阵测向方法、装置及电子设备
CN111914506B (zh) * 2020-08-05 2024-05-14 中国电子科技集团公司第三十六研究所 天线阵测向方法、装置及电子设备
CN113468844A (zh) * 2021-06-17 2021-10-01 浙江大学 一种耦合阵列波束综合的解析方法
CN113468844B (zh) * 2021-06-17 2023-08-04 浙江大学 一种耦合阵列波束综合的解析方法

Also Published As

Publication number Publication date
CN111257823B (zh) 2022-03-22

Similar Documents

Publication Publication Date Title
US10663563B2 (en) On-site calibration of array antenna systems
CN111273224B (zh) 一种基于可视化阵列天线的测量方法
CN111257823B (zh) 一种基于耦合阵列的测向方法
CN111693978B (zh) 基于mimo毫米波雷达的散点检测方法
CN106707255B (zh) 一种相控阵雷达仿真系统及方法
CN111711500B (zh) 一种仿真天线阵列标校与射频信号监测系统
CN105388442B (zh) 一种基于可移动短路板的自由空间法校准方法
CN108168559B (zh) 一种基于分布式天线的室内定位系统及方法
Jam et al. A submillimeter-wave near-field measurement setup for on-wafer pattern and gain characterization of antennas and arrays
CN105445555B (zh) 一种相控阵天线单元有源驻波比计算方法
CN109839543B (zh) 一种天线的幅相一致性测试系统及测试方法
CN110208736B (zh) 基于四阶累量的非圆信号均匀阵列波达方向角估计方法
CN107942145A (zh) 一种角变化球面三元组天线辐射电磁场测量系统和方法
CN108717188B (zh) 适用于mimo穿墙雷达成像的快速墙体补偿方法
WO2023016571A1 (zh) 射频指标测量方法、装置、系统、电子设备及存储介质
CN115002653A (zh) 基于一体式uwb基站的轻量级二维高精度定位方法
Mueller et al. Broadband biquad UHF antenna array for DOA
CN111505391B (zh) 一种检测阵列天线中的失配天线单元的方法
CN110174639B (zh) 一种准确提取干涉仪天线阵列单元相位中心的方法
CN110444886B (zh) 一种降低由极化变化引起的天线间相位差误差的方法
Gu et al. Analysis of amplitude-phase error of phased array calibration in mid-field
Tan et al. UWB-AOA estimation method based on a spare antenna array with virtual element
Luchin et al. Algorithm of radio direction finding in the HF range in polarization fading conditions with collocated antennas
Dang et al. A simple method for calculating the output voltage at receive antenna terminals using the complex antenna factor
CN110320489B (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