CN111220942A - 一种接收换能器阵列幅相一致性近场校准方法 - Google Patents
一种接收换能器阵列幅相一致性近场校准方法 Download PDFInfo
- Publication number
- CN111220942A CN111220942A CN201911257833.5A CN201911257833A CN111220942A CN 111220942 A CN111220942 A CN 111220942A CN 201911257833 A CN201911257833 A CN 201911257833A CN 111220942 A CN111220942 A CN 111220942A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- array
- angle
- phase
- phase consistency
- 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
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
- G01S3/00—Direction-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/80—Direction-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 ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种接收换能器阵列幅相一致性近场校准方法。首先,处于回转台上的接收阵在不同方向采集的入射信号经过聚焦波束形成处理得到一系列的波达方向估计;然后,在频域建立接收信号协方差矩阵,并建立幅相一致性误差估计模型;选取小角度入射角信号,将其邻域角度分别假设为理想入射角,根据幅相一致性误差估计模型估计对应邻域内角度的误差估计结果;最后,根据误差结果补偿阵列,根据试验预设参数选择两已知角度间隔的数据,评估对于角度间隔估计的偏差,根据偏差极小值确定幅相一致性误差,完成校准。本发明能够在不需要辅助信源的情况下提高阵列的幅相一致性及测向性能。
Description
技术领域
本发明涉及的是一种阵列天线的校准和测量方法,具体地说是一种接收换能器阵列幅相一致性近场校准方法。
背景技术
信号波达方向估计中需要采用波束形成方法对来波方位进行扫描,而一些波束形成算法的测向性能极大依赖于阵元间的幅度相位信息,接收换能器阵列由于长时间使用,阵元加工工艺差异,以及环境变化等会不可避免引发水声换能器阵元间幅度相位的不一致性问题,进一步导致测向精度下降、波束图旁瓣级升高。
现有的接收换能器阵列幅相一致性校准方法通常需要借助辅助信源获得方位信息,并且大多采用平面波假设模型。然而借助辅助信源的被动校准方法更易引入误差,并且平面波的假设模型仅适用于远场,一旦被校准对象孔径较大且阵元数较多,若继续采用平面波模型,则校准方法对校准场地尺寸有着较高的要求,而这通常是难以满足的。此外,由于一个周期内时域信号的采样点数目较多,在涉及协方差矩阵运算时,现有校准方法多在时域进行,运算量较大。
近场聚焦波束形成基于球面波假设能够对声源位置聚焦,相对于平面波假设模型,它可以保证距离发射声源较近时依然具有精确的测向结果,为测量条件提供了便利。尤其是在对孔径较大且阵元数较多的换能器阵列进行校准时,相比于传统校准方法具有天然的优势。
发明内容
本发明的目的在于提供一种在校准场地尺寸有限且无辅助信源的情况下能够提高接收换能器阵列幅相一致性的接收换能器阵列幅相一致性近场校准方法。
本发明的目的是这样实现的:
(1)建立近场聚焦波束形成模型,处于回转台上的接收阵在不同方向采集的入射信号经过聚焦波束形成处理得到一系列的波达方向估计;
(2)在频域建立接收信号协方差矩阵,并建立幅相一致性误差估计模型;
(3)选取小角度入射角信号,将其邻域角度分别设为理想入射角,根据幅相一致性误差估计模型估计对应邻域内角度的幅相误差;
(4)根据误差结果补偿阵列,根据预设参数选择两已知角度间隔的数据,评估对于角度间隔估计的偏差,根据偏差极小值确定幅相一致性误差,完成校准。
本发明还可以包括:
1.所述的近场聚焦模型为:
2.步骤(2)具体包括:
建立幅相一致性误差估计模型,Γ和Φ分别作为阵列的幅度相位向量有:
Γ=diag[ρ1,ρ2,…,ρM]
则理想接收信号以及存在误差的接收信号表示为:
X0(t)=AS(t)+N(t)
X(t)=ΓΦ(AS(t)+N(t))
S(t)=[s1(t),s2(t),…,sM(t)]T(s1(t)=s2(t),…,=sM(t))
N(t)=[n1(t),n2(t),…,nM(t)]T
其中,X0(t)和X(t)为以理想接收信号和实际接收信号矩阵;S(t)为信源发射信号组成的发射信号矩阵,N(t)为各阵元接收的噪声组成的噪声矩阵,其中噪声设为高斯白噪声;
将各阵元接收的理想信号X0(t)和实际信号X(t)经过傅里叶变换至频域得到X0 fre(f)和Xfre(f),分别提取距离零频最近频谱峰值X0 P和XP:
Xfre(f)=[x1(f),x2(f),…,xM(f)]T
XP=[P1,P2,…,PM]T
3.步骤(3)中,选择小角度的入射角为θ0的接收信号,在角度范围为ε的邻域中对邻域内K个角度遍历,使得理想角度θs满足θs(k)∈(θ0-ε,θ0+ε),计算对应邻域内各个角度的幅相一致性误差Ω(k)。
4.步骤(4)具体为:
对k组各通道幅度误差ρm(k)和相位误差分别补偿阵列,选择两入射角间隔为Δθ的接收信号,其中参数Δθ0根据回转台获得,分别估计这两组信号的波达方向并计算其角度间隔Δθ(k),考察对角度间隔估计的偏差e(k)=|Δθ0-Δθ(k)|,若在角度邻域(θ0-ε,θ0+ε)内存在唯一极小值,则理想角度为换能器阵列各阵元的幅相误差为和
与现有技术相比,本发明的有益效果是:无需辅助信源而采用回转台提供的角度间隔信息进行校准,试验简便且更易实际操作。相比于平面波假设模型,采用基于球面波假设的近场聚焦波束形成处理方法估计结果更精确,能够克服传统校准方法对于校准场地的尺寸难以满足的问题,通用性强。此外,本发明的校准方法在频域中的一次快拍即可完成协方差矩阵的构建,有效减少了计算量,并且频域处理使得信号在解调频后无需低通滤波环节就能提取有效信息。
附图说明
图1为本发明的流程图;
图2为角度间隔估计偏差示意图;
图3a为阵元间幅度校准前后示意图;
图3b为阵元间相位校准前后示意图;
图4为幅相一致性误差校准前后的波束形成图。
具体实施方式
下面举例对本发明做更详细的描述。
结合图1至图4,本发明的接收换能器阵列幅相一致性近场校准方法主要包括以下几个步骤:
(1)处于回转台上的接收阵在不同方向采集的入射信号经过聚焦波束形成处理得到一系列的波达方向估计;
(2)在频域建立接收信号协方差矩阵,并建立幅相一致性误差估计模型;
(3)选取小角度入射角信号,将其邻域角度分别假设为理想入射角,根据幅相一致性误差估计模型估计对应邻域内角度的幅相误差;
(4)根据误差结果补偿阵列,根据试验预设参数选择两已知角度间隔的数据,评估对于角度间隔估计的偏差,根据偏差极小值确定幅相一致性误差,完成校准。
回转台作为待校准接收阵的载体能够采集不同方向上入射的信号,能够获得任意两组采样数据的入射角度间隔,通过聚焦波束形成估计多组数据的波达方向,选择具有小角度入射的数据遍历其角度邻域,采用协方差矩阵算法计算幅相误差,考察补偿后的阵列对具有两已知入射角间隔的数据的估计偏差,通过偏差极小值存在情况,调节邻域范围,最终确定幅相一致性误差,校准的流程如图1所示。
步骤(1)中涉及的近场聚焦模型为:
其中dm为第m个阵元相对于参考阵元的间距,λ为声波的波长,rn表示在第n个采样时刻下对应的聚焦距离。
步骤(2)建立幅相一致性误差估计模型,Γ和Φ分别作为阵列的幅度相位向量有:
Γ=diag[ρ1,ρ2,…,ρM]
则理想接收信号以及存在误差的接收信号可以表示为:
X0(t)=AS(t)+N(t)
X(t)=ΓΦ(AS(t)+N(t))
S(t)=[s1(t),s2(t),…,sM(t)]T(s1(t)=s2(t),…,=sM(t))
N(t)=[n1(t),n2(t),…,nM(t)]T
其中X0(t)和X(t)为以理想接收信号和实际接收信号矩阵。S(t)为信源发射信号组成的发射信号矩阵,N(t)为各阵元接收的噪声组成的噪声矩阵,其中噪声假设为高斯白噪声;
将各阵元接收的理想信号X0(t)(各阵元幅相一致性误差一致)和实际信号X(t)经过傅里叶变换至频域得到X0 fre(f)和Xfre(f),分别提取距离零频最近频谱峰值X0 P和XP:
Xfre(f)=[x1(f),x2(f),…,xM(f)]T
XP=[P1,P2,…,PM]T
步骤(3)中选择小角度的入射角为θ0的接收信号,在角度范围为ε的邻域中对邻域内K个角度遍历,使得理想角度θs满足θs(k)∈(θ0-ε,θ0+ε),计算对应邻域内各个角度的幅相一致性误差Ω(k)。
步骤(4)对k组各通道幅度误差ρm(k)和相位误差分别补偿阵列,选择两入射角间隔为Δθ的接收信号,其中参数Δθ0根据回转台获得,分别估计这两组信号的波达方向并计算其角度间隔Δθ(k),考察对角度间隔估计的偏差e(k)=|Δθ0-Δθ(k)|,若在角度邻域(θ0-ε,θ0+ε)内存在唯一极小值,则理想角度为换能器阵列各阵元的幅相误差为和
对接收换能器阵列幅相一致性近场校准方法进行了仿真试验分析,待校准的阵列为一阵元数为M=100的均匀线列阵,阵元间距d=3.75mm,由于阵列安装在回转台上使得各数据间入射角度间隔为已知,回转台设定旋转范围为-90°~90°,回转速度为0.6°/s,信号源与接收阵距离为10m处于近场范围,发射高频窄带脉冲信号,其中发射信号中心频率为200kHz,脉冲宽度为0.1ms,信号发射的同时接收阵接收,周期均为1s,则共有300组接收信号,且每相邻两组数据的角度间隔为0.6°,其中以第150组数据的入射角度为理想的参考值,设为0.8°。引入幅度相位误差,各阵元幅度误差满足ρ~N(1,0.32),相位误差满足存在幅相一致性误差情况下使得第150组数据的波束图的旁瓣级升高并导致波达方向估计偏差为0.9°。以0.9°对附近邻域进行搜索,由于设定的理想值为0.8°,则选择邻域为0.7°~0.9°进行遍历,每次遍历增量为0.01°,共21个角度,通过幅相一致性误差估计模型计算得到21组阵元间幅度相位误差,根据21组幅相误差对阵列补偿前,为结果更加明显应选择入射角较大的数据,适当选择两角度间隔为24°的两组数据,选择第100组数据和第140组数据,在每一组幅相误差补偿后,记录对两组数据的波达方向估计结果。经过处理,得到对应不同组数据的角度间隔估计误差。由图2可知,第11组幅相误差对阵列补偿后具有最佳的估计结果,对应的入射角度为0.8°与理想角度参考值一致,并且在邻域内具有唯一的极小值,因此阵列的幅度相位误差可以确定,幅相误差真实值与估计值的对比如图3a和图3b,阵列校准前后对第150组数据进行处理得到的波束图如图4。经过校准,阵列的测向精度得到了提高,波束图的旁瓣得到了抑制。
综上,本发明提供了一种接收换能器阵列幅相一致性近场校准方法。无需辅助信源即可完成校准,尤其是在对孔径较大且阵元数较多的换能器阵列进行校准时,采用的基于近场球面波模型能够克服校准场地尺寸难以满足测量要求的现状。在频域中构造协方差矩阵简化了解调频步骤和计算复杂度,本发明提高了阵列的幅相一致性及测向性能,操作简便,适合工程应用。
Claims (5)
1.一种接收换能器阵列幅相一致性近场校准方法,其特征是包括如下步骤:
(1)建立近场聚焦波束形成模型,处于回转台上的接收阵在不同方向采集的入射信号经过聚焦波束形成处理得到一系列的波达方向估计;
(2)在频域建立接收信号协方差矩阵,并建立幅相一致性误差估计模型;
(3)选取小角度入射角信号,将其邻域角度分别设为理想入射角,根据幅相一致性误差估计模型估计对应邻域内角度的幅相误差;
(4)根据误差结果补偿阵列,根据预设参数选择两已知角度间隔的数据,评估对于角度间隔估计的偏差,根据偏差极小值确定幅相一致性误差,完成校准。
3.根据权利要求2所述的接收换能器阵列幅相一致性近场校准方法,其特征是步骤(2)具体包括:
建立幅相一致性误差估计模型,Γ和Φ分别作为阵列的幅度相位向量有:
Γ=diag[ρ1,ρ2,…,ρM]
则理想接收信号以及存在误差的接收信号表示为:
X0(t)=AS(t)+N(t)
X(t)=ΓΦ(AS(t)+N(t))
S(t)=[s1(t),s2(t),…,sM(t)]T(s1(t)=s2(t),…,=sM(t))
N(t)=[n1(t),n2(t),…,nM(t)]T
其中,X0(t)和X(t)为以理想接收信号和实际接收信号矩阵;S(t)为信源发射信号组成的发射信号矩阵,N(t)为各阵元接收的噪声组成的噪声矩阵,其中噪声设为高斯白噪声;
将各阵元接收的理想信号X0(t)和实际信号X(t)经过傅里叶变换至频域得到X0 fre(f)和Xfre(f),分别提取距离零频最近频谱峰值X0 P和XP:
Xfre(f)=[x1(f),x2(f),…,xM(f)]T
XP=[P1,P2,…,PM]T
4.根据权利要求3所述的接收换能器阵列幅相一致性近场校准方法,其特征是步骤(3)中,选择小角度的入射角为θ0的接收信号,在角度范围为ε的邻域中对邻域内K个角度遍历,使得理想角度θs满足θs(k)∈(θ0-ε,θ0+ε),计算对应邻域内各个角度的幅相一致性误差Ω(k)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911257833.5A CN111220942B (zh) | 2019-12-10 | 2019-12-10 | 一种接收换能器阵列幅相一致性近场校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911257833.5A CN111220942B (zh) | 2019-12-10 | 2019-12-10 | 一种接收换能器阵列幅相一致性近场校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111220942A true CN111220942A (zh) | 2020-06-02 |
CN111220942B CN111220942B (zh) | 2023-01-03 |
Family
ID=70808130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911257833.5A Active CN111220942B (zh) | 2019-12-10 | 2019-12-10 | 一种接收换能器阵列幅相一致性近场校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111220942B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111964757A (zh) * | 2020-08-13 | 2020-11-20 | 浙江苍南仪表集团股份有限公司 | 一种基于回波信号特征参数的换能器一致性评价方法 |
CN113114390A (zh) * | 2021-04-30 | 2021-07-13 | 网络通信与安全紫金山实验室 | 基于宽带信号的阵列ota阵形校准方法及校准*** |
CN117630894A (zh) * | 2024-01-24 | 2024-03-01 | 山东省科学院海洋仪器仪表研究所 | 一种多波束声呐平均阵元间距校准方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644765A (zh) * | 2009-04-23 | 2010-02-10 | 中国科学院声学研究所 | 一种用于水声换能器线阵的幅相误差校正方法 |
EP2396851A1 (en) * | 2009-02-13 | 2011-12-21 | Socowave Technologies Limited | Communication system, apparatus and methods for calibrating an antenna array |
CN102495396A (zh) * | 2011-11-15 | 2012-06-13 | 北京无线电计量测试研究所 | 人体安检***多通道幅度和相位一致性校准方法 |
CN103018739A (zh) * | 2011-09-20 | 2013-04-03 | 中国科学院电子学研究所 | 一种校正多通道幅相误差的微波三维成像方法 |
CN105044660A (zh) * | 2015-06-18 | 2015-11-11 | 许昌学院 | 一种阵列通道不一致性误差快速校正方法 |
CN106371055A (zh) * | 2016-10-26 | 2017-02-01 | 黑龙江大学 | 远近场宽带混合信号超分辨测向阵列幅相误差估计方法 |
CN106788799A (zh) * | 2016-11-22 | 2017-05-31 | 电子科技大学 | 一种处理多种阵列天线信号误差的优化方法 |
EP3229383A1 (fr) * | 2016-04-04 | 2017-10-11 | Thales | Système et procèdè de calibration dynamique d'une ou plusieurs chaines radiofréquence de transmission d'une charge utile de satellite |
CN108872926A (zh) * | 2018-07-11 | 2018-11-23 | 哈尔滨工程大学 | 一种基于凸优化的幅相误差校正及doa估计方法 |
CN109376329A (zh) * | 2018-09-05 | 2019-02-22 | 哈尔滨工程大学 | 一种基于量子鸟群演化机制的阵列幅相误差校正方法 |
-
2019
- 2019-12-10 CN CN201911257833.5A patent/CN111220942B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2396851A1 (en) * | 2009-02-13 | 2011-12-21 | Socowave Technologies Limited | Communication system, apparatus and methods for calibrating an antenna array |
CN101644765A (zh) * | 2009-04-23 | 2010-02-10 | 中国科学院声学研究所 | 一种用于水声换能器线阵的幅相误差校正方法 |
CN103018739A (zh) * | 2011-09-20 | 2013-04-03 | 中国科学院电子学研究所 | 一种校正多通道幅相误差的微波三维成像方法 |
CN102495396A (zh) * | 2011-11-15 | 2012-06-13 | 北京无线电计量测试研究所 | 人体安检***多通道幅度和相位一致性校准方法 |
CN105044660A (zh) * | 2015-06-18 | 2015-11-11 | 许昌学院 | 一种阵列通道不一致性误差快速校正方法 |
EP3229383A1 (fr) * | 2016-04-04 | 2017-10-11 | Thales | Système et procèdè de calibration dynamique d'une ou plusieurs chaines radiofréquence de transmission d'une charge utile de satellite |
CN106371055A (zh) * | 2016-10-26 | 2017-02-01 | 黑龙江大学 | 远近场宽带混合信号超分辨测向阵列幅相误差估计方法 |
CN106788799A (zh) * | 2016-11-22 | 2017-05-31 | 电子科技大学 | 一种处理多种阵列天线信号误差的优化方法 |
CN108872926A (zh) * | 2018-07-11 | 2018-11-23 | 哈尔滨工程大学 | 一种基于凸优化的幅相误差校正及doa估计方法 |
CN109376329A (zh) * | 2018-09-05 | 2019-02-22 | 哈尔滨工程大学 | 一种基于量子鸟群演化机制的阵列幅相误差校正方法 |
Non-Patent Citations (1)
Title |
---|
WEN, ZL等: ""Calibration Method of Amplitude and Phase Consistency of W-Band Transceiver Components"", 《PROCEEDINGS OF 2014 3RD ASIA-PACIFIC CONFERENCE ON ANTENNAS AND PROPAGATION (APCAP 2014)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111964757A (zh) * | 2020-08-13 | 2020-11-20 | 浙江苍南仪表集团股份有限公司 | 一种基于回波信号特征参数的换能器一致性评价方法 |
CN113114390A (zh) * | 2021-04-30 | 2021-07-13 | 网络通信与安全紫金山实验室 | 基于宽带信号的阵列ota阵形校准方法及校准*** |
CN117630894A (zh) * | 2024-01-24 | 2024-03-01 | 山东省科学院海洋仪器仪表研究所 | 一种多波束声呐平均阵元间距校准方法 |
CN117630894B (zh) * | 2024-01-24 | 2024-04-12 | 山东省科学院海洋仪器仪表研究所 | 一种多波束声呐平均阵元间距校准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111220942B (zh) | 2023-01-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111220942B (zh) | 一种接收换能器阵列幅相一致性近场校准方法 | |
CN101644765B (zh) | 一种用于水声换能器线阵的幅相误差校正方法 | |
CN110196414B (zh) | 一种基于补偿天线方向图误差的天线波束指向方法 | |
CN110058193B (zh) | 一种基于单接收通道的数字多波束测角方法及*** | |
CN109856605A (zh) | 一种同时数字多波束形成的二次拟合曲线指向修正方法 | |
CN111025273B (zh) | 一种畸变拖曳阵线谱特征增强方法及*** | |
CN111948598B (zh) | 空域干扰信号检测方法与装置 | |
CN104007413A (zh) | 考虑信源方位误差的阵列位置误差校正方法 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN114137476B (zh) | 一种基于极化修正的短波测向方法 | |
CN112630784A (zh) | 基于凸优化和神经网络的平面阵列幅相误差校正方法 | |
CN110687538A (zh) | 一种基于近场聚焦的超波束形成方法 | |
CN105277932B (zh) | 一种基于下变频波束形成中的多普勒频移校正方法 | |
CN110703219B (zh) | 一种多发多收近场直线阵列获取目标远场rcs的方法 | |
CN116430303A (zh) | 宽带平面阵列多波束形成方法及比幅测角方法 | |
CN111896913A (zh) | 高频雷达单极子/交叉环天线通道增益校准方法及装置 | |
CN111175727B (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN114679227A (zh) | 一种测向误差的空间频域校正方法 | |
Yuan et al. | Correction method for magnitude and phase variations in acoustic arrays based on focused beamforming | |
CN111505590A (zh) | 一种高频地波雷达通道校准方法及*** | |
CN109613474B (zh) | 一种适用于短距离车载雷达的测角补偿方法 | |
CN109471065B (zh) | 一种相干信号的测向方法 | |
CN111007505A (zh) | 一种基于零陷估计的距离通道相位偏差估计方法及*** | |
CN110174639B (zh) | 一种准确提取干涉仪天线阵列单元相位中心的方法 | |
CN107560639B (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 |