CN104991121A - 一种扫频阵列校准方法 - Google Patents
一种扫频阵列校准方法 Download PDFInfo
- Publication number
- CN104991121A CN104991121A CN201510334031.5A CN201510334031A CN104991121A CN 104991121 A CN104991121 A CN 104991121A CN 201510334031 A CN201510334031 A CN 201510334031A CN 104991121 A CN104991121 A CN 104991121A
- Authority
- CN
- China
- Prior art keywords
- array
- frequency
- phase
- calibration method
- frequency sweep
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种扫频阵列校准方法,其改进之处在于,包括以下步骤:步骤1,校准源发射扫频信号;步骤2,阵列接收设备依次采集单频信号;步骤3,对步骤2采集的单频信号进行积累处理;步骤4,对步骤3所得的相位矢量Pn按阵列相位模型进行拟合,使得实际测量数据与理论模型的偏差最小,得到拟合相位矢量Qn;步骤5,利用步骤4所得的拟合相位矢量Qn,得到阵列的通道相位差异矢量Q0。本发明所公开的扫频阵列校准方法,通过多频点的整体相位模型对阵列进行校准,无需预知校准源位置,利用统计减小了到达角预计等因素带来的偏差,消除了校准值随频率的波动。
Description
技术领域
本发明涉及阵列校准领域,尤其涉及一种扫频阵列校准方法。
背景技术
常规的阵列校准方式通常逐个频点进行,存在以下问题:1、需要预知校准源的方位;2、单频校准结果中存在难以消除的偏差;3、受各种干扰因素影响,校准值随频率波动较大。这些问题在高频段更加突出。由于天线阵列需要的面积较大,阵列往往受周围环境影响,使得来波的波前面改变,不能准确反映到达角。因此,即使知道校准源方位,也不能正确计算其到达角,给校准结果带来较大偏差。
发明内容
本发明所要解决的技术问题就是提供一种无需预知校准源位置并可有效抑制随机干扰影响的扫频阵列校准方法。
本发明采用如下技术方案:
一种扫频阵列校准方法,其改进之处在于,包括以下步骤:
步骤1,校准源发射扫频信号,扫频信号包括频率依次等间隔递增的N个单频信号,扫频的起始频率为f0;
步骤2,阵列接收设备依次采集单频信号,阵列接收设备包括L个通道,L>1,第l通道采集的第n个单频信号为sl,n;
步骤3,对步骤2采集的单频信号进行积累处理,由sl,n积累出一个标量,标量的相位表征了sl,n的相位,获取该相位;sl,n的相位为pl,n,所有通道的相位构成阵列相位矢量Pn=[p1,n,p2,n,…,pL,n];
步骤4,对步骤3所得的相位矢量Pn按阵列相位模型进行拟合,使得实际测量数据与理论模型的偏差最小,得到拟合相位矢量Qn;
步骤5,利用步骤4所得的拟合相位矢量Qn,得到阵列的通道相位差异矢量Q0。
进一步的,步骤1中所述的扫频信号恰好覆盖阵列接收设备的工作频段,N>30,所述单频信号采用可积累的信号形式。
进一步的,步骤2中所述的通道具有稳定的相位特性。
进一步的,步骤4中所述的阵列相位模型表达式为ql,n=ql,1+l·(n-1)·α·fΔ,其中ql,n为所述pl,n的拟合量,ql,n=pl,n+el,n,所述el,n为pl,n对所述阵列相位模型的偏差,所述el,n为α的函数,所述fΔ为所述扫频信号的频率步进,所述α为常数,与所述扫频信号的到达角相关。
进一步的,步骤4中所述的拟合过程准则为所述el,n的统计特征最小,其执行步骤为:
步骤41,构筑所述el,n的统计量E,所述E同样为α的函数,可记为E(α);
步骤42,搜索所述α,使得所述E最小,记为α|E→min,简记为αm;
步骤43,将所述αm代入所述阵列相位模型,可得到所述拟合相位矢量Qn=[q1,n,q2,n,…,qL,n]。
进一步的,步骤5中所述通道相位差异矢量Q0,其表达式为Q0=[q1,0,q2,0,…,qL,0],ql,0= ql,1- l·α·f1。
本发明的有益效果在于:
本发明所公开的扫频阵列校准方法,通过多频点的整体相位模型对阵列进行校准,无需预知校准源位置,利用统计减小了到达角预计等因素带来的偏差,消除了校准值随频率的波动。具体的说,校准源发射扫频信号,阵列接收设备采集和处理每个频率的信号,获取其相位。所有通道的相位构成阵列相位矢量,该矢量为工作频率的函数,符合特定的模型。将实测数据与模型进行拟合,可以得到通道间的固有差异。相对传统技术,该方法无需预知校准源位置,可有效抑制随机干扰的影响。
附图说明
图1为本发明实施例1所使用的阵列布局图;
图2为本发明实施例1所公开的校准流程图;
图3为使用本发明实施例1所公开的校准方法校准前后的比较图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明所公开的扫频阵列校准方法,校准源发射的扫频信号,恰好覆盖阵列接收设备的工作频段,一边反映全频段的特性。N尽量多,使得统计特性更加准确,一般取N>30。所述单频信号采用可积累的信号形式,如单频信号和线性调频信号。单频信号经过傅里叶变换可得到一个谱峰;线性调频信号通过脉冲压缩也可获得一个窄峰。该峰值点可作为表征信号的标量,其相位代表信号的相位。阵列接收设备的通道具有稳定的相位特性,不随时间大幅变化。阵列相位模型表达式为ql,n=ql,1+l·(n-1)·α·fΔ。其中,ql,n为pl,n的拟合量,ql,n=pl,n+el,n;el,n为pl,n对阵列相位模型的偏差,α的函数; fΔ为扫频信号的频率步进;α为常数,与所述扫频信号的到达角相关。拟合的准则为所述el,n的统计特征最小。典型的统计特征为均方根。其执行步骤可描述如下:步骤4-1,构筑所述el,n的统计量E,所述E同样为α的函数,可记为E(α);步骤4-2,搜索所述α,使得所述E最小,记为α|E→min,简记为αm;步骤4-3,将所述αm代入所述阵列相位模型,可得到所述拟合相位矢量Qn=[q1,n,q2,n,…,qL,n]。通道相位差异矢量Q0=[q1,0,q2,0,…,qL,0],其中,ql,0= ql,1- l·α·f1。
实施例1,将本发明所公开的扫频阵列校准方法用于如图1所示阵列,阵列为8元均匀线阵线,阵元1至阵元8从右至左依次排列,阵列间距8m,校准源位于偏离法线30°、距离约10km的位置,校准源发射功率在1W左右。
校准过程如图2所示,包括:发射扫频信号、采集信号、获取相位信息、拟合、计算通道相位差异,共5个步骤。
步骤1,校准源发射扫频信号,扫频信号频率为4MHz~23.9MHz,频率步进为0.1MHz,共N=100个频率,信号形式为线性调频,每个频率发射512个连续的线性调频脉冲。
步骤2,阵列接收设备依次采集单频信号,采集过程使用IQ分离技术。阵列接收设备包括L=8个通道,与阵元一一对应,记为通道1至通道8。通道具有稳定的相位特性,稳定度优于3°/24h。第l通道采集的第n个单频信号为sl,n。sl,n为一段采样序列,每个采样点表示为一个复数,具有幅度和相位信息。
步骤3,对步骤2采集的单频信号进行积累处理。首先对单个线性调频脉冲进行脉压,然后对脉压结果进行累加。累加结果的峰值相位为pl,n,可表征通道相位,所有通道的相位构成阵列相位矢量Pn=[p1,n,p2,n,…,pL,n]。Pn随频率变化的曲线如图3左侧所示。
步骤4,对步骤3所得的相位矢量Pn按阵列相位模型进行拟合,使得实际测量数据与理论模型的偏差最小,得到拟合相位矢量Qn。根据阵列信号处理的基本理论,阵列相位模型表达式为ql,n=ql,1+l·(n-1)·α·fΔ。其中,ql,n为pl,n的拟合量,ql,n=pl,n+el,n;el,n为pl,n对阵列相位模型的偏差,是α的函数; fΔ为扫频信号的频率步进;α为常数,与所述扫频信号的到达角相关。本例中α≈2π·8m·sin(-30°)/(3e8m/s)≈8.38e-8s,该值从实际测量曲线中也可以大致看出。拟合采用最小均方根准则,选择适当的α,使el,n的均方根最小。具体步骤如下:
步骤41,构筑el,n的均方根E(α)=∑l∑n(el,n)2/(L·N);步骤42,搜索所述α,使得所述E最小,记为α|E→min,简记为αm,本例中可将搜索范围确定在8.38e-8s附近;步骤43,将所述αm代入所述阵列相位模型,可得到所述拟合相位矢量Qn=[q1,n,q2,n,…,qL,n]。Qn随频率变化的曲线如图3右侧所示。
步骤5,利用步骤4所得的拟合相位矢量Qn,得到阵列的通道相位差异矢量Q0;Q0为阵列的固有特性,与频率无关。通道相位差异矢量Q0=[q1,0,q2,0,…,qL,0],ql,0= ql,1- l·α·4MHz。
Claims (6)
1.一种扫频阵列校准方法,其特征在于,包括以下步骤:
步骤1,校准源发射扫频信号,扫频信号包括频率依次等间隔递增的N个单频信号,扫频的起始频率为f0;
步骤2,阵列接收设备依次采集单频信号,阵列接收设备包括L个通道,L>1,第l通道采集的第n个单频信号为sl,n;
步骤3,对步骤2采集的单频信号进行积累处理,由sl,n积累出一个标量,标量的相位表征了sl,n的相位,获取该相位;sl,n的相位为pl,n,所有通道的相位构成阵列相位矢量Pn=[p1,n,p2,n,…,pL,n];
步骤4,对步骤3所得的相位矢量Pn按阵列相位模型进行拟合,使得实际测量数据与理论模型的偏差最小,得到拟合相位矢量Qn;
步骤5,利用步骤4所得的拟合相位矢量Qn,得到阵列的通道相位差异矢量Q0。
2.根据权利要求1所述的扫频阵列校准方法,其特征在于:步骤1中所述的扫频信号恰好覆盖阵列接收设备的工作频段,N>30,所述单频信号采用可积累的信号形式。
3.根据权利要求1所述的扫频阵列校准方法,其特征在于:步骤2中所述的通道具有稳定的相位特性。
4.根据权利要求1所述的扫频阵列校准方法,其特征在于:步骤4中所述的阵列相位模型表达式为ql,n=ql,1+l·(n-1)·α·fΔ,其中ql,n为所述pl,n的拟合量,ql,n=pl,n+el,n,所述el,n为pl,n对所述阵列相位模型的偏差,所述el,n为α的函数,所述fΔ为所述扫频信号的频率步进,所述α为常数,与所述扫频信号的到达角相关。
5.根据权利要求4所述的扫频阵列校准方法,其特征在于:步骤4中所述的拟合过程准则为所述el,n的统计特征最小,其执行步骤为:
步骤41,构筑所述el,n的统计量E,所述E同样为α的函数,可记为E(α);
步骤42,搜索所述α,使得所述E最小,记为α|E→min,简记为αm;
步骤43,将所述αm代入所述阵列相位模型,可得到所述拟合相位矢量Qn=[q1,n,q2,n,…,qL,n]。
6.根据权利要求1所述的扫频阵列校准方法,其特征在于:步骤5中所述通道相位差异矢量Q0,其表达式为Q0=[q1,0,q2,0,…,qL,0],ql,0= ql,1- l·α·f1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510334031.5A CN104991121B (zh) | 2015-06-16 | 2015-06-16 | 一种扫频阵列校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510334031.5A CN104991121B (zh) | 2015-06-16 | 2015-06-16 | 一种扫频阵列校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104991121A true CN104991121A (zh) | 2015-10-21 |
CN104991121B CN104991121B (zh) | 2018-01-16 |
Family
ID=54302958
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510334031.5A Active CN104991121B (zh) | 2015-06-16 | 2015-06-16 | 一种扫频阵列校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104991121B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109214088A (zh) * | 2018-09-07 | 2019-01-15 | 厦门大学深圳研究院 | 一种最小间距可控的大规模超稀疏平面阵列快速布局方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3062880B2 (ja) * | 1996-09-10 | 2000-07-12 | テクトロニクス・インコーポレイテッド | 位相測定装置及び方法 |
EP2402772A1 (en) * | 2009-02-27 | 2012-01-04 | Furuno Electric Co., Ltd. | Phase determining device and frequency determining device |
CN102404033A (zh) * | 2011-11-24 | 2012-04-04 | 北京交通大学 | 一种ofdm系统中的天线阵列校准方法和装置 |
CN102468805A (zh) * | 2010-11-03 | 2012-05-23 | 北京普源精电科技有限公司 | 一种扫频信号发生器及扫频信号的产生方法 |
CN104037499A (zh) * | 2014-05-30 | 2014-09-10 | 中国电子科技集团公司第二十二研究所 | 一种有源短波天线及其组成的有源短波天线阵 |
-
2015
- 2015-06-16 CN CN201510334031.5A patent/CN104991121B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3062880B2 (ja) * | 1996-09-10 | 2000-07-12 | テクトロニクス・インコーポレイテッド | 位相測定装置及び方法 |
EP2402772A1 (en) * | 2009-02-27 | 2012-01-04 | Furuno Electric Co., Ltd. | Phase determining device and frequency determining device |
CN102468805A (zh) * | 2010-11-03 | 2012-05-23 | 北京普源精电科技有限公司 | 一种扫频信号发生器及扫频信号的产生方法 |
CN102404033A (zh) * | 2011-11-24 | 2012-04-04 | 北京交通大学 | 一种ofdm系统中的天线阵列校准方法和装置 |
CN104037499A (zh) * | 2014-05-30 | 2014-09-10 | 中国电子科技集团公司第二十二研究所 | 一种有源短波天线及其组成的有源短波天线阵 |
Non-Patent Citations (1)
Title |
---|
宋孝果: "天地波高频雷达阵列校准和直达波抑制", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109214088A (zh) * | 2018-09-07 | 2019-01-15 | 厦门大学深圳研究院 | 一种最小间距可控的大规模超稀疏平面阵列快速布局方法 |
CN109214088B (zh) * | 2018-09-07 | 2023-05-02 | 厦门大学 | 一种最小间距可控的大规模超稀疏平面阵列快速布局方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104991121B (zh) | 2018-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109471095B (zh) | 一种基于快速迭代插值的fmcw雷达距离估计方法 | |
CN104515909B (zh) | 一种基于相关法的大天线方向图测量方法 | |
CN113075461B (zh) | 基于宽带甚高频辐射信号探测的超短基线闪电三维定位方法 | |
CN105510766A (zh) | 射频电缆故障定位检测装置及检测方法 | |
CN109212527B (zh) | 用于高频地波雷达的大孔径分布式多站目标定位方法 | |
CN103472463A (zh) | 一种卫星导航接收机设备时延标定方法 | |
CN107390192A (zh) | 相控阵天气雷达快速幅度、相位一致性测量方法 | |
CN104535959A (zh) | 时空欠采样下信号频率及doa联合测量方法及装置 | |
CN104270208B (zh) | 一种远端射频单元rru驻波比检测的方法及装置 | |
US9960862B2 (en) | Method and device for detecting standing-wave ratio | |
CN111580137B (zh) | 一种高精度导航接收机射频通道群时延特性的拟合方法 | |
CN104914393A (zh) | 一种用于梳状谱发生器相位谱校准的装置及方法 | |
CN102788909B (zh) | 采用全息法测量天线面精度的方法 | |
CN103968864B (zh) | 用于准确测量布里渊谱的频移的最大相似匹配分析方法 | |
CN111190197A (zh) | 一种导航卫星信号质量在轨优化与保持方法 | |
CN104991121A (zh) | 一种扫频阵列校准方法 | |
CN102540153A (zh) | 一种基于站间直达波干扰的阵列幅相误差校准方法 | |
CN109752705B (zh) | 高频水声阵列性能参数测量方法及系统、设备及存储介质 | |
CN103869334A (zh) | Gnss空间信号畸变的自动识别与处理方法 | |
US10649018B2 (en) | Time alignment of lightning emissions at LF-MF using waveform feature comparison | |
CN105785225A (zh) | 接收端信号时延估算误差的补偿方法 | |
CN110456317B (zh) | 一种基于流星尾迹回波的相控阵雷达系统定标方法 | |
CN104880697A (zh) | 基于稀疏约束的线性调频信号参数估计方法 | |
CN104579516A (zh) | 一种驻波比检测方法和设备 | |
CN110531298A (zh) | 特高频局部放电传感器性能改进检测方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |