CN1664611A - 一种基于非直线天线阵列的无源通道校正方法 - Google Patents
一种基于非直线天线阵列的无源通道校正方法 Download PDFInfo
- Publication number
- CN1664611A CN1664611A CN2005100184383A CN200510018438A CN1664611A CN 1664611 A CN1664611 A CN 1664611A CN 2005100184383 A CN2005100184383 A CN 2005100184383A CN 200510018438 A CN200510018438 A CN 200510018438A CN 1664611 A CN1664611 A CN 1664611A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- array
- msup
- mtr
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000012937 correction Methods 0.000 claims description 48
- 238000002592 echocardiography Methods 0.000 claims description 47
- 238000005457 optimization Methods 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 9
- 238000013519 translation Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 abstract 2
- 230000009977 dual effect Effects 0.000 abstract 2
- 239000011159 matrix material Substances 0.000 abstract 2
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 238000001228 spectrum Methods 0.000 description 20
- 238000001514 detection method Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000004220 aggregation Methods 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 241000283973 Oryctolagus cuniculus Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/0218—Very long range radars, e.g. surface wave radar, over-the-horizon or ionospheric propagation systems
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/536—Discriminating between fixed and moving objects or between objects moving at different speeds using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves
-
- 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
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于非直线天线阵列的无源通道校正方法,将天线阵列设置为含有平移不变阵元偶组的非直线形式;通过平移不变阵元偶组检测出单方位海洋回波;利用单方位海洋回波估计通道幅度失配系数,实现幅度校正;利用已经过幅度校正的单方位海洋回波和已知的阵列位置信息估计通道相位失配系数,实现相位校正。本发明完全不需要任何辅助信号源,是一种真正的无源通道校正方法,其优势在于:不存在有源通道校正方法面临的船只回波干扰、多径效应等棘手问题;利用了大量高强度单方位海洋回波,具有良好的精度和稳健性;运算量较小;能长期不间断稳定工作;大大改善了雷达的应用灵活性;在提高探测性能的同时,大幅降低了雷达的研制成本和维护费用。
Description
技术领域
本发明涉及一种基于非直线天线阵列对高频地波雷达进行无源通道校正的方法。
背景技术
高频地波雷达利用高频电磁波在导电海洋表面绕射传播衰减小的特点,采用垂直极化天线辐射电波,能超视距探测海平面视线以下出现的舰船、飞机和导弹等运动目标。另外,高频地波雷达利用海洋表面对高频电磁波的一阶散射和二阶散射机制,从雷达回波中提取风场、浪场、流场等海态信息,可实现对海洋环境大范围、高精度和全天候的实时监测。
由于硬件本身的差异、接收通道的非理想特性、周边环境影响等多种因素的作用,实际中天线阵列各个通道的幅相特性是有差异的,导致回波信号经过不同通道后的幅度和相位变化(复增益)不一致,一般称之为通道失配。通道失配引起波束扫描和方位估计的误差增大,甚至完全失效,是影响高频地波雷达探测性能的关键问题之一。为了保证雷达能有效工作,必须采取措施使通道失配限制在一定的范围内:一方面,通过适当的措施(如元器件筛选)使各通道在制作时尽量保证其一致性;另一方面,可以对通道失配系数进行测量或估计,然后通过校正进一步缩小通道特性的差异。
现有的通道校正方法可分为有源和无源两类。在有源校正方法中,将辅助信号源置于天线阵前方足够远的开阔场地发射校正信号,再测量各接收通道的输出,根据已知信号源方位扣除阵列空间位置引起的相位差,即可得到通道失配信息。在无源校正方法中,无需方位准确已知的辅助信号源,直接利用接收的实测数据和一些先验知识(如阵列形式)估计通道失配系数,然后进行补偿校正。有些无源校正方法还可以实现信号方位和通道失配的联合估计。在刘德树、罗景青等编著的《空间谱估计及其应用》(中国科技大学出版社1997年)一书中对一般无源校正方法有详细阐述。
受地形条件、工作波长、电波传播、雷达体制、天线阵列、(硬)目标回波、海洋杂波、噪声干扰等多种因素的影响,高频地波雷达的通道校正实现起来比较困难,现有方法都只能解决部分问题,且费时费力,成本高昂,代价很大。雷达天线阵前方是海面,如采用有源校正方法则辅助信号源只能放在舰船或岛屿上,其维护十分麻烦而昂贵,难以长期稳定工作。现有的无源校正方法一般需要多次复杂的迭代运算,计算量很大,不一定能满足实时性要求,且有可能收敛于局部最优,而不是全局最优,以致出现完全错误的估计值。实际雷达系统与理想模型的差异,也使得一般无源校正方法的适用条件不能得到满足,难以实用。通道校正已成为制约高频地波雷达探测性能并阻碍其实际应用的重大技术难题,必须加以妥善解决。
武汉大学电波传播实验室曾考虑过将海上已知天然或人工物体对雷达电波的反射信号作为校正信号。只要反射源的距离、速度已知,就可以从回波中检测出校正信号,然后根据已知反射源方位估计出各通道的失配系数。其具体实施细节可参考03128238.5号中国发明专利申请“一种利用海洋回波进行阵列通道校正的方法”。该发明可利用探测海域内已知的岛屿、灯塔和钻井平台等固定反射物,不存在辅助信号源的放置和维护问题,也无需额外的硬件开销,实现了在线实时自动校正,具有一定的实用价值。但该发明实际上是一种特殊的有源校正方法,应用范围和实际效果有限,不适用于没有已知固定反射物的海域,且仍然受到噪声干扰、舰船回波、多径效应等不利因素的影响。由于该发明提出了一种对频谱不重叠的单方位海洋回波的分离与检测技术,满足了本发明提出的无源通道校正方法的基本要求,因此下面将对其进行重点介绍。
高频地波雷达一般采用调频中断连续波(frequency modulated interrupted continuouswave,缩写FMICW)体制。Rafaat Khan等人发表的题为“高频地波雷达目标探测与跟踪”(Target Detection and Tracking With a High Frequency Ground Wave Radar,IEEE Journal ofOceanic Engineering,1994,19(4):540~548)的论文中对此有详细描述。在该波形体制下,海洋回波(包括海杂波和硬目标回波)进入接收机后,经混频、低通滤波、A/D转换和二维FFT(如图1所示)可得距离-多普勒(速度)二维回波谱(如图2所示)。在二维回波谱中,雷达接收的大量海洋回波按距离和速度进行了分离,分散在很多谱点上。当第二次FFT(多普勒变换)相干积累时间较长(10分钟左右)时,雷达可获得很高的速度分辨力,二维回波谱中与海洋回波对应的谱点可达1000个以上,很适合用统计方法将其中频谱不重叠的单方位回波检测出来。
单方位回波的检测是通过对特定形式阵列(如图3所示)二维回波谱输出的统计分析实现的。特定形式阵列由阵元1~4构成,其位置坐标为(xi,yi),对应的某个二维回波谱点输出为Yi,i=1,2,3,4。阵元1和2组成阵元偶A1,3和4组成阵元偶A2,A1与A2之间具有平移不变性,则有
令
,容易证明,在没有噪声的理想情况下,二维回波谱中单方位谱点对应的η1是一个只与通道失配有关的不变量,不妨记为η1′。实际系统中噪声是不可避免的,单方位谱点对应的η1集中分布在η1′附近。另一方面,通过简单分析和数值模拟可知,多方位谱点对应的η1是一个与目标距离、(径向)速度、方位和回波信号振幅均有关的变量,由于这些目标参数的随机性,η1将呈随机分散状态。综合以上分析,将二维回波谱中所有超过一定信噪比门限的谱点对应的η1标在复平面上,则有且仅有一个区域(η1′附近)出现高度聚集现象,其中大多数η1值对应于单方位谱点。令
,由与前面类似的分析可知,η2在复平面上也会出现聚集区,其中大多数η2值对应于单方位谱点。令
,η3在复平面上同样会出现聚集区,其中大多数η3值对应于单方位谱点。
经过理论分析和数值模拟发现,多方位谱点对应的η1、η2和η3同时落入各自聚集区概率极小,因此可以用η1、η2和η3能否同时落入各自聚集区作为检测单方位谱点的判据。A1和A2作为一个平移不变阵元偶组,构成了检测单方位回波(谱点)的特定形式阵列。若阵列中平移不变阵元偶组超过一个,则将是否同时被多个阵元偶组检测到作为判据,可进一步筛选出单方位谱点。含有平移不变阵元偶组的阵列形式是很常见的,如均匀线阵(或均匀平面阵)。
发明内容
针对现有方法的局限性,本发明的目的是利用高频地波雷达接收的单方位海洋回波,提供一种基于非直线天线阵列的无源通道校正方法,以减少通道幅相失配,提高雷达系统性能。这里的非直线天线阵列是指并非所有天线阵元都位于同一直线上的阵列。
为了实现上述目的,本发明提供的技术方案是:一种基于非直线天线阵列的无源通道校正方法,将天线阵列设置为含有平移不变阵元偶组的非直线形式;通过平移不变阵元偶组检测出单方位海洋回波;利用单方位海洋回波估计通道幅度失配系数,实现幅度校正;利用已经过幅度校正的单方位海洋回波和已知的阵列位置信息估计通道相位失配系数,实现相位校正。
在上述方法中,可利用单方位海洋回波通过
估计通道幅度失配系数,实现幅度校正,其中
为阵元i的通道幅度失配系数估计值,i=1,2,…,M,M为阵元个数,Yi(l)为阵元i接收的第l个单方位回波输出,l=1,2,…,L,L为单方位回波个数。
本发明还可以利用已经过幅度校正的单方位海洋回波和已知的阵列位置信息通过 估计通道相位失配系数,实现相位校正,其中
Ψ=[θ1,θ2,…,θL,φ2,φ3,…,φM]T
Yi=[Yi(1),Yi(2),…,Yi(L)]T
θl为第l个单方位回波的到达角,φi为阵元i的通道相位失配系数,(xi,yi)为阵元位置坐标,阵元1为坐标原点,即(x1,y1)=(0,0),λ为回波信号波长,
为Ψ的估计值。
本发明利用对三角阵的处理降低全局最优化的维数,并通过初值预估尽可能采用局部最优化方法,以降低多维参数估计的运算量。
1)先从所有阵元中选出3个阵元构成一个三角阵,再选出3个单方位回波用于参数估计;2)以该三角阵中的任一阵元为基准通道,采用全局最优化方法估计另两个通道的相位失配系数和3个单方位回波的到达角;3)增加一个单方位回波用于三角阵的参数估计,得出该单方位回波的到达角;4)按步骤3)得出其它单方位回波的到达角;5)将三角阵与另一个阵元组合成4元阵,再把所有单方位回波用于该4元阵的参数估计,得到新加入阵元的通道相位失配系数;6)按步骤5)得出其它阵元的通道相位失配系数。
本发明也可在经过步骤4)后,将所有单方位回波用于三角阵的参数估计,以已求得的单方位回波到达角和通道相位失配系数的估计值为初值,采用局部最优化方法求得这些参数更精确的估计值,然后进行后面的步骤5)和6)。
在经过步骤6后,将所有单方位回波用于整个阵列的参数估计,以已求得的单方位回波到达角和通道相位失配系数的估计值为初值,采用局部最优化方法求得这些参数更精确的估计值。
对于M元L形阵,将阵元1、2、M作为步骤1中所选的三角阵用于参数估计,实现通道相位校正。
对于4元T形阵,将阵元1、2、4或阵元2、3、4作为步骤1中所选的三角阵用于参数估计,实现通道相位校正。
本发明的优势在于其出色的实用性能:完全不需要任何辅助信号源,是一种真正的无源通道校正方法。本发明只利用了单方位回波,不存在有源通道校正方法面临的船只回波干扰、多径效应等棘手问题;利用了大量(可达100个以上)高强度单方位海洋回波,信息利用率较高,具有良好的精度和稳健性;采用一些特殊处理时运算量较小,能满足实时性要求;由于海洋回波总是大量存在,该方法能长期不间断稳定工作;大大改善了雷达的应用灵活性,天线系统可随意更换、增减、移动,这在以前是难以想象的;在提高探测性能的同时,大幅降低了雷达的研制成本和维护费用。
附图说明
图1为高频地波雷达工作原理图;
图2为高频地波雷达距离-多普勒(速度)二维回波谱图;
图3为用于检测单方位回波的特定形式阵列示意图;
图4为本发明M元任意非直线阵列示意图;
图5为三角阵示意图;
图6为M元L形阵示意图;
图7为4元矩形阵示意图;
图8为4元T形阵示意图;
具体实施方式
本发明的关键在于通过建立非直线天线阵列接收的单方位回波信号模型,将通道校正问题转化成了参数估计问题,并由此得到了较为精确的通道失配估计。
首先考虑图4所示的M(M≥3)元任意非直线天线阵列的情况,以下分步骤对本发明在该情况下的具体实施方式进行描述。
(一)信号模型
设图4所示的非直线天线阵列的阵元位置坐标分别为(xi,yi)(i=1,2,…,M),其中阵元1为坐标原点,即(x1,y1)=(0,0)。海洋回波可看作平面波,假设从距离-多普勒(速度)二维回波谱中检测出的单方位回波有L(L≥3)个,则阵元i接收的第l(l=1,2,…,L)个单方位回波输出为
其中,θl和A(l)分别为第l个单方位回波的到达角和复振幅,gi和φi分别为阵元i的通道幅度失配系数和相位失配系数,λ为回波信号波长,Wi(l)为加性噪声。以阵元1为基准通道,即
,则由(1)式有
Y1(l)=A(l)+W1(l) (2)
对加性噪声Wi(l)采用以下假设:
1)对应于不同i或l的Wi(l)相互独立;
2)Wi(l)为具有相同方差σ2的高斯白噪声。
则(1)和(2)式构成了通道失配估计的信号模型。
(二)通道幅度校正
阵元i的通道幅度失配估计可直接通过
(三)通道相位校正
阵列接收通道经过幅度校正后,阵元i接收的第l个单方位回波输出为
定义
Yi=[Yi(1),Yi(2),…,Yi(L)]T i=2,3,…,M
Ψ=[θ1,θ2,…,θL,φ2,φ3,…,φM]T
以Ψ为待估参数矢量,根据(4)式和噪声模型假设,采用最大似然法(请参考张贤达所著《现代信号处理》,清华大学出版社,1994年)可得Ψ的估计值为
其中,
表示后面表达式极小化时Ψ的取值,‖X‖2表示矢量X的2-范数。显然,通道相位失配系数和单方位回波到达角实现了联合估计。
A(l)包含在有噪声的Yi(l)中,是不可能直接得到的,fi(Ψ)也是不可能直接构造的,实际中用于通道相位失配估计的不可能是(5)式,需要进行改进。在一般信噪比条件(≥20dB)下,由(2)式可得Y1(l)≈A(l),用Y1(l)替换(5)式中的A(l),得到Ψ的估计值实际表达式为
其中
根据
即可对通道相位失配进行校正。
理论分析和仿真实验表明,上述通道相位校正方法的适用条件是采用非直线阵列且存在至少两个角度相差不为0°或180°的单方位回波。通道相位失配估计由于涉及所有阵元,实际上是一个多维参数估计问题,要通过多维搜索求得,这就涉及最优化方法的选择。由于存在局部极小,必须采用全局最优化方法(请参考唐芳、王凌发表的论文“从局部极小到全局最优”,计算机工程与应用,2002.6:56~58)对通道相位失配进行估计,在参数较多(可达100个以上)的情况下运算量十分惊人,达不到实时性要求。本发明采用了特殊处理以降低运算量,下面将对此作详细说明。
图5所示为最简单的非直线阵列:由不共直线的3个阵元构成的三角阵。若只利用3个单方位回波联合估计该三角阵的通道相位失配系数和单方位回波到达角,则根据(6)式可判断为5维搜索(有一个阵元作为基准通道)。由于维数较少,即使采用模拟退火、进化计算、混沌搜索、随机抽样等全局最优化方法,该5维搜索的运算量也不大,可满足实时性要求。
对于图4所示的M元任意非直线阵列,为了降低多维参数估计的运算量,可利用其中包含的某个三角阵先对待估参数进行初值预估,再采用局部最优化方法(如最速下降法)求得更精确估计,具体步骤如下:
1)先从所有阵元中选出3个阵元构成一个三角阵,再从大量单方位回波中选出3个用于参数估计。
2)以该三角阵中的某个阵元为基准通道,采用全局最优化方法估计另两个通道的相位失配系数和3个单方位回波的到达角,这是5维搜索。
3)增加一个单方位回波用于三角阵的参数估计,则(6)式中待估参数也增加了一个(即新增单方位回波的到达角),由5维搜索变成了6维搜索;将步骤2)求得的5个参数估计值代入含有6个待估参数的(6)式,再采用只剩一个待估参数的(6)式估计新增单方位回波的到达角,这是1维搜索。
4)除了步骤1)中选出的3个单方位回波,其它L-3个单方位回波的到达角估计值都可采用步骤3)中的方法求得。
5)将所有L个单方位回波用于三角阵的参数估计,以已求得的L+2个参数(L个单方位回波到达角和2个通道相位失配系数)估计值为初值,采用局部最优化方法求得这些参数更精确的估计值。
6)将三角阵与另一个阵元组合成4元阵,再把所有L个单方位回波用于该4元阵的参数估计,则(6)式中待估参数为L+3个;将步骤5)已求得的L+2个参数估计值代入(6)式,则(6)式只剩一个待估参数,即新加入阵元的通道相位失配系数,可采用1维搜索求得其估计值。
7)除了步骤1)中选出的三角阵,其它M-3个阵元的通道相位失配估计值都可采用步骤6)中的方法求得。
8)将所有L个单方位回波用于整个M元阵列的参数估计,以已求得的L+M-1个参数(L个单方位回波到达角和M-1个通道相位失配系数)估计值为初值,采用局部最优化方法求得这些参数更精确的估计值。
在参数初值预估误差不大的情况下,局部最优化所得结果也是全局最优,而计量量却比全局最优化小得多。以上特殊处理中的初值预估实际上占用了大部分计算时间,但其中最多只涉及5维搜索的全局最优化,因而能满足实时性要求。以上8个步骤是特殊处理的典型方式,在实际应用中可根据具体情况进行简化、充实、调整或改进,其基本思想是利用对三角阵的处理降低全局最优化的维数,并通过初值预估尽可能采用局部最优化方法,以达到降低多维参数估计运算量的目的。
图6所示为M元L形阵,是本发明较常见的一个实施例。阵元1~M-1构成一个均匀线阵,其中可划分出多个用于检测单方位回波的平移不变阵元偶组。通道幅度校正可通过(3)式实现,关键是通道相位校正。将阵元1、2、M构成的三角阵用于参数初值预估,通过(6)式和特殊处理即可实现通道相位校正。
图7所示为4元矩形阵。在该实施例中,4个阵元仅构成了一个可用于检测单方位回波的平移不变阵元偶组,其中任意3个阵元构成的三角阵都可用于参数初值预估。
图8所示为4元T形阵。在该实施例中,阵元1~3构成3元均匀线阵,只能划分出一个用于检测单方位回波的平移不变阵元偶组。阵元1、2、4和阵元2、3、4构成的两个三角阵都可用于参数初值预估。
虽然本发明描述的通道校正方法首先在高频地波雷达上获得了成功,但从本质上讲,该方法也有可能应用于其它接收了大量单方位信号的探测系统或通信系统。
Claims (9)
1.一种基于非直线天线阵列的无源通道校正方法,其特征在于:将天线阵列设置为含有平移不变阵元偶组的非直线形式;通过平移不变阵元偶组检测出单方位海洋回波;利用单方位海洋回波估计通道幅度失配系数,实现幅度校正;利用已经过幅度校正的单方位海洋回波和已知的阵列位置信息估计通道相位失配系数,实现相位校正。
2.根据权利要求1所述的方法,其特征在于:利用单方位海洋回波通过 估计通道幅度失配系数,实现幅度校正,其中
为阵元i的通道幅度失配系数估计值,i=1,2,…,M,M为阵元个数,Yi(l)为阵元i接收的第l个单方位回波输出,l=1,2,…,L,L为单方位回波个数。
4.根据权利要求3所述的方法,其特征在于:利用对三角阵的处理降低全局最优化的维数,并通过初值预估尽可能采用局部最优化方法,以降低多维参数估计的运算量。
5.根据权利要求4所述的方法,其特征在于:1)先从所有阵元中选出3个阵元构成一个三角阵,再选出3个单方位回波用于参数估计;2)以该三角阵中的任一阵元为基准通道,采用全局最优化方法估计另两个通道的相位失配系数和3个单方位回波的到达角;3)增加一个单方位回波用于三角阵的参数估计,得出该单方位回波的到达角;4)按步骤3)得出其它单方位回波的到达角;5)将三角阵与另一个阵元组合成4元阵,再把所有单方位回波用于该4元阵的参数估计,得到新加入阵元的通道相位失配系数;6)按步骤5)得出其它阵元的通道相位失配系数。
6.根据权利要求5所述的方法,其特征在于:在经过步骤4)后,将所有单方位回波用于三角阵的参数估计,以已求得的单方位回波到达角和通道相位失配系数的估计值为初值,采用局部最优化方法求得这些参数更精确的估计值,然后进行后面步骤。
7.根据权利要求5所述的方法,其特征在于:在经过步骤6)后,将所有单方位回波用于整个阵列的参数估计,以已求得的单方位回波到达角和通道相位失配系数的估计值为初值,采用局部最优化方法求得这些参数更精确的估计值。
8.根据权利要求5所述的方法,其特征在于:对于M元L形阵,将阵元1、2、M作为步骤1)中所选的三角阵用于参数估计,实现通道相位校正。
9.根据权利要求5所述的方法,其特征在于:对于4元T形阵,将阵元1、2、4或阵元2、3、4作为步骤1)中所选的三角阵用于参数估计,实现通道相位校正。
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2005100184383A CN1664611A (zh) | 2005-03-24 | 2005-03-24 | 一种基于非直线天线阵列的无源通道校正方法 |
AU2006227098A AU2006227098B2 (en) | 2005-03-24 | 2006-03-21 | A passive channel adjustment method based on a non-linear antenna array |
PCT/CN2006/000453 WO2006099800A1 (fr) | 2005-03-24 | 2006-03-21 | Procede de reglage de canal passif base sur un reseau d’antennes |
EP06722105A EP1879043A4 (en) | 2005-03-24 | 2006-03-21 | METHOD FOR PASSIVE CHANNEL ADJUSTMENT BASED ON A NONLINEAR ANTENNA GROUP |
CN200610071360A CN100585429C (zh) | 2005-03-24 | 2006-03-22 | 一种基于非直线天线阵列的无源通道校正方法 |
US11/859,746 US7667639B2 (en) | 2005-03-24 | 2007-09-22 | Passive channel calibration method based on non-linear antenna array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2005100184383A CN1664611A (zh) | 2005-03-24 | 2005-03-24 | 一种基于非直线天线阵列的无源通道校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1664611A true CN1664611A (zh) | 2005-09-07 |
Family
ID=35035795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2005100184383A Pending CN1664611A (zh) | 2005-03-24 | 2005-03-24 | 一种基于非直线天线阵列的无源通道校正方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7667639B2 (zh) |
EP (1) | EP1879043A4 (zh) |
CN (1) | CN1664611A (zh) |
AU (1) | AU2006227098B2 (zh) |
WO (1) | WO2006099800A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1804656B (zh) * | 2006-01-20 | 2010-04-14 | 武汉大学 | 一种利用电离层回波进行高频雷达天线阵列通道校正的方法 |
CN101923157A (zh) * | 2010-07-29 | 2010-12-22 | 西安空间无线电技术研究所 | 一种星载双通道角跟踪校准系统及方法 |
CN104267386A (zh) * | 2014-10-17 | 2015-01-07 | 武汉大学 | 一种外辐射源雷达旋转阵列的无源通道校正方法及系统 |
CN104698444A (zh) * | 2013-12-23 | 2015-06-10 | 中国科学院电子学研究所 | 确定阵元位置误差对阵列天线三维成像质量影响的方法 |
CN107783086A (zh) * | 2017-10-17 | 2018-03-09 | 中国电子科技集团公司第三十八研究所 | 用于诊断天线阵口径幅相场的畸变位置的方法 |
CN109856627A (zh) * | 2019-01-11 | 2019-06-07 | 中国船舶重工集团公司第七二四研究所 | 基于空间补零插值的三角阵相控阵2d-fft测向方法 |
CN111505391A (zh) * | 2020-05-09 | 2020-08-07 | 北京航空航天大学云南创新研究院 | 一种检测阵列天线中的失配天线单元的方法 |
CN112986935A (zh) * | 2021-02-24 | 2021-06-18 | 湖北中南鹏力海洋探测系统工程有限公司 | 一种基于地貌共形布阵的多元阵列通道无源校准方法 |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102012004308B4 (de) * | 2012-03-01 | 2018-09-20 | Atlas Elektronik Gmbh | Annäherungssensor und Verfahren zum Bestimmen der Annäherung an einen elektrisch leitfähigen Körper |
CN104854473B (zh) * | 2012-11-05 | 2017-05-24 | 泰克年研究发展基金会公司 | 使用多普勒聚焦的亚奈奎斯特雷达处理 |
JP5788860B2 (ja) * | 2012-12-28 | 2015-10-07 | 日本電信電話株式会社 | 無線通信システムおよび無線通信方法 |
RU2540982C1 (ru) * | 2013-09-10 | 2015-02-10 | ОТКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО "НИИ измерительных приборов-Новосибирский завод имени Коминтерна" /ОАО "НПО НИИИП-НЗиК"/ | Способ определения координат целей (варианты) и комплекс для его реализации (варианты) |
RU2567114C1 (ru) * | 2014-05-05 | 2015-11-10 | Открытое акционерное общество "Азовский оптико-механический завод" (ОАО "АОМЗ") | Система для измерения координат объекта навигации |
US10175203B2 (en) * | 2014-09-11 | 2019-01-08 | Cpg Technologies, Llc | Subsurface sensing using guided surface wave modes on lossy media |
US10101444B2 (en) | 2014-09-11 | 2018-10-16 | Cpg Technologies, Llc | Remote surface sensing using guided surface wave modes on lossy media |
CN104378170B (zh) * | 2014-10-27 | 2017-02-08 | 西安电子科技大学 | 一种临近空间动态等离子鞘套信道建模及模拟的方法 |
CN105554780B (zh) * | 2015-12-23 | 2019-03-12 | 哈尔滨工业大学 | 毫米波下Massive MIMO多小区协作波束分配方法 |
US20170307755A1 (en) * | 2016-04-20 | 2017-10-26 | YoR Labs | Method and System for Determining Signal Direction |
DE102017204496A1 (de) * | 2017-03-17 | 2018-09-20 | Robert Bosch Gmbh | Verfahren und Radarvorrichtung zum Ermitteln von radialer relativer Beschleunigung mindestens eines Zieles |
CN110007267B (zh) * | 2019-01-29 | 2020-08-18 | 杭州电子科技大学 | 一种基于混合基线的均匀圆阵干涉仪测向解模糊方法 |
CN110286370B (zh) * | 2019-05-08 | 2022-03-22 | 自然资源部第一海洋研究所 | 航行状态下的船载地波雷达运动补偿方法 |
US11998391B1 (en) | 2020-04-02 | 2024-06-04 | yoR Labs, Inc. | Method and apparatus for composition of ultrasound images with integration of “thick-slice” 3-dimensional ultrasound imaging zone(s) and 2-dimensional ultrasound zone(s) utilizing a multi-zone, multi-frequency ultrasound image reconstruction scheme with sub-zone blending |
US11547386B1 (en) | 2020-04-02 | 2023-01-10 | yoR Labs, Inc. | Method and apparatus for multi-zone, multi-frequency ultrasound image reconstruction with sub-zone blending |
US11832991B2 (en) | 2020-08-25 | 2023-12-05 | yoR Labs, Inc. | Automatic ultrasound feature detection |
US11344281B2 (en) | 2020-08-25 | 2022-05-31 | yoR Labs, Inc. | Ultrasound visual protocols |
US11704142B2 (en) | 2020-11-19 | 2023-07-18 | yoR Labs, Inc. | Computer application with built in training capability |
US11751850B2 (en) | 2020-11-19 | 2023-09-12 | yoR Labs, Inc. | Ultrasound unified contrast and time gain compensation control |
WO2023033835A1 (en) * | 2021-09-03 | 2023-03-09 | Zeku, Inc. | Apparatus and method of beam parameter estimation for line-of-sight determination |
CN115130362B (zh) * | 2022-07-08 | 2023-04-07 | 中国科学院声学研究所 | 一种基于信道匹配的阵列阵元位置校准方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5361072A (en) * | 1992-02-28 | 1994-11-01 | Codar Ocean Sensors, Ltd. | Gated FMCW DF radar and signal processing for range/doppler/angle determination |
FR2762397B1 (fr) * | 1997-04-18 | 1999-07-09 | Thomson Csf | Dispositif d'auto test de la chaine d'emission et de reception d'un radar, notamment pour automobile |
CN1321331C (zh) * | 2003-06-30 | 2007-06-13 | 武汉大学 | 一种利用海洋回波进行阵列通道校正的方法 |
CN1328594C (zh) * | 2004-05-17 | 2007-07-25 | 中国人民解放军空军雷达学院 | 用高精度辅助阵元进行阵列校正与信源测向的方法 |
-
2005
- 2005-03-24 CN CN2005100184383A patent/CN1664611A/zh active Pending
-
2006
- 2006-03-21 EP EP06722105A patent/EP1879043A4/en not_active Withdrawn
- 2006-03-21 WO PCT/CN2006/000453 patent/WO2006099800A1/zh not_active Application Discontinuation
- 2006-03-21 AU AU2006227098A patent/AU2006227098B2/en not_active Ceased
-
2007
- 2007-09-22 US US11/859,746 patent/US7667639B2/en not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1804656B (zh) * | 2006-01-20 | 2010-04-14 | 武汉大学 | 一种利用电离层回波进行高频雷达天线阵列通道校正的方法 |
CN101923157A (zh) * | 2010-07-29 | 2010-12-22 | 西安空间无线电技术研究所 | 一种星载双通道角跟踪校准系统及方法 |
CN101923157B (zh) * | 2010-07-29 | 2013-05-01 | 西安空间无线电技术研究所 | 一种星载双通道角跟踪校准系统 |
CN104698444A (zh) * | 2013-12-23 | 2015-06-10 | 中国科学院电子学研究所 | 确定阵元位置误差对阵列天线三维成像质量影响的方法 |
CN104698444B (zh) * | 2013-12-23 | 2017-04-26 | 中国科学院电子学研究所 | 确定阵元位置误差对阵列天线三维成像质量影响的方法 |
CN104267386A (zh) * | 2014-10-17 | 2015-01-07 | 武汉大学 | 一种外辐射源雷达旋转阵列的无源通道校正方法及系统 |
CN107783086A (zh) * | 2017-10-17 | 2018-03-09 | 中国电子科技集团公司第三十八研究所 | 用于诊断天线阵口径幅相场的畸变位置的方法 |
CN107783086B (zh) * | 2017-10-17 | 2021-01-05 | 中国电子科技集团公司第三十八研究所 | 用于诊断天线阵口径幅相场的畸变位置的方法 |
CN109856627A (zh) * | 2019-01-11 | 2019-06-07 | 中国船舶重工集团公司第七二四研究所 | 基于空间补零插值的三角阵相控阵2d-fft测向方法 |
CN111505391A (zh) * | 2020-05-09 | 2020-08-07 | 北京航空航天大学云南创新研究院 | 一种检测阵列天线中的失配天线单元的方法 |
CN111505391B (zh) * | 2020-05-09 | 2022-05-31 | 北京航空航天大学云南创新研究院 | 一种检测阵列天线中的失配天线单元的方法 |
CN112986935A (zh) * | 2021-02-24 | 2021-06-18 | 湖北中南鹏力海洋探测系统工程有限公司 | 一种基于地貌共形布阵的多元阵列通道无源校准方法 |
Also Published As
Publication number | Publication date |
---|---|
EP1879043A4 (en) | 2013-01-30 |
EP1879043A1 (en) | 2008-01-16 |
AU2006227098A1 (en) | 2006-09-28 |
AU2006227098B2 (en) | 2010-07-29 |
US7667639B2 (en) | 2010-02-23 |
US20080012753A1 (en) | 2008-01-17 |
WO2006099800A1 (fr) | 2006-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1664611A (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
CN1740812A (zh) | 高频表面波雷达均匀直线阵接收通道的近场校正方法 | |
CN110456334B (zh) | 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法 | |
CN108549059B (zh) | 一种复杂地形条件下的低空目标仰角估计方法 | |
CN103901417B (zh) | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 | |
CN103529437B (zh) | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 | |
CN105137410B (zh) | 基于ofdm的高分辨率雷达通信一体化的波形优化方法 | |
CN109765562A (zh) | 一种三维前视声像声纳系统和方法 | |
CN1847877A (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
CN110412587B (zh) | 一种基于解卷积的下视合成孔径三维成像方法及系统 | |
CN110346798B (zh) | 一种双基合成孔径雷达波数域高效成像处理方法 | |
CN109188386B (zh) | 基于改进二维esprit算法的mimo雷达高分辨率参数估计方法 | |
CN104076337B (zh) | 基于阵元幅相误差校正的机载雷达杂波抑制方法 | |
CN105807275A (zh) | 基于部分杂波先验知识的mimo-ofdm-stap稳健波形设计 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN109765521A (zh) | 一种基于子阵划分的波束域成像方法 | |
CN110412570A (zh) | 基于空间脉冲相位编码的hrws-sar成像方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN114488148B (zh) | 一种基于稀疏tops-sar成像模式实现方法 | |
CN111896913A (zh) | 高频雷达单极子/交叉环天线通道增益校准方法及装置 | |
CN109946655B (zh) | Mimo雷达lfm信号置零波形设计方法 | |
CN108872947B (zh) | 一种基于子空间技术的海杂波抑制方法 | |
CN114114188A (zh) | 一种低旁瓣的fda雷达通信一体化波形设计方法 | |
CN104122548B (zh) | 机扫米波阵列雷达测角方法 | |
CN111929682B (zh) | 基于距离向分块的机载dbf-tops sar成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |