相位干涉仪测向算法及其在tms320c6711上的实现.docVIP

相位干涉仪测向算法及其在TMS320C6711上的实现 ; 相位干涉仪测向技术广泛应用于天文、雷达、声纳等领域。将干涉仪原理用于无线电测向始于上世纪五十年代和六十年代,随着数字信号处理器的出现,通过数字信号处理器来实现高精度实时测向成为可能。; 本文在对一维和二维相位干涉仪进行研究的基础上给出了五通道相位干涉仪的基本关系式,分析了测向精度,并对解相位模糊问题和信道校正问题进行了探讨。采用多基线五元圆形天线阵列为模型,由天线阵列接收到的信号求解出五元天线阵列的互相关信号,并由此提取测向所需的方位。本文以五通道相位干涉仪硬件实现为目标,采用高速浮点数字信号处理芯片TMS320C6711进行测向处理。; 1 相位干涉仪测向原理; 1.1 一维相位干涉仪测向原理; 图1所示为一个最简单的一维双阵元干涉仪模型。图中,间隔为d(d称为基线)的两根天线A1和A2所接收的远场辐射信号之间的相位差为:; 鸡眨剑ǎ处校洌λ)cosθ; (1); 式(1)中,λ为接收电磁波的波长。因此,只要测量出φ,就能算出辐射源的到达方向θ:; θ=arccos(φλ/4πd)nbsp;nbsp; (2); 1.2 测向误差的分析; 在实际系统中,两根天线A1和A2接收的信号为:; xi(t)=s(t)exp[(-1)jj2πd/λcosθ]+ni(t),i=1,2 (3); 其中,ni代表对应阵元i接收的噪声,两阵元的噪声统计相互独立,且与信号统计独立。; 两个阵元接收信号的互相关为:; r=E{x1(t)x2*(t)}=Psexp(j4πd/λcosθ)nbsp;nbsp; (4); 式中,E代表期望运算,“*”代表复共轭运算,Ps代表信号功率,相关以后噪声得到抑制。; 由(4)式有:; θ=arccos[(λ/4πd)arg(r21)+kλ/2d; (5); 式中,arccos表示反余弦函数,arg代表复数取幅角运算,区间为[-π,π]。k为整数,且满足:; -2d/λ-arg(r21)/2π≤k≤2d/λ-arg(r21)/2π; (6); 在(6)式中,当d/λgt;0.5时,k的取值不唯一,θ有多个解,由此产生测向模糊。; 对(5)式求导,有:; |Δθ|=λ/4πd|sinθ|Δarg(r21); (7); 由(7)式可以得出以下结论:sinθ越大,即方位角与干涉仪法线方向的夹角越小,测向精度越高;反之,测向精度降低,直至测向无效。当θ=±90°(即信号从干涉仪法线方向入射)时,精度最高;θ=0°或180°(即信号从干涉仪基线方向入射)时, 接收信号互相关的幅角arg(r21)反映不出方位角的变化,测向无效。但单基线干涉仪不能同时测量俯仰角和方位角,此时至少需要另一条独立基线的干涉仪对测得的数据联合求解。; 1.3 二维干涉仪测向原理及去模糊处理; 1.3.1 多基线五元圆形天线模型; 五通道相位干涉仪采用宽口径、多基线的五元圆形天线阵,五边形的五个阵元均匀分布在半径为R的圆上,五个阵源分别为1、2、3、4、5,如图2所示。天线阵平面与地面平行,测得的方位角θ为以天线到地面的垂足为原点,目标在地面上的方位角。测得的俯仰角φ对应于目标到原点的距离(俯仰角0°对应原点)。; 两个阵元接收信号之间的互相关为:; ri,j+1=E{xi(t)x*i+1(t)}=GiGi+1Psexp{j2π(r/λ)sinφ?[cos(θ+54°-72°i)-cos(θ-18°-72°i)]}; i=1~5,定义r56=r51; 方位角θ和俯仰角φ的具体计算如下:; Q ri,i+1的幅角为αi,i+1=arg(ri,i+1)+2k2π=4π(r/λ)cos54°sinφcos(θ+108°-72°i); ri+3,i+4的幅角为αi+3,i+4=arg(ri+3,i+4)+2k1π=4π(r/λ)cos54°sinφcos(θ-108°-72°i); ∴θ=atan2[αi+3,i+4-αi,i+1)csc108°,(αi+3,i+4+αi,i+1)sec108°]+72°i; (8); 式中,i=1~5,令r56=r51、r67=r12、r78=r23、r89=r34;atan2(y,x)代表四象限求反正切函数;arcsin代表反正弦函数。k1、k2为整数,且满足:; (4r/λ)sinφcos54°sin108°-[arg(ri+3,i+4)-arg(ri,i+1)]/2π≤k1-k2≤(4r/λ)sinφcos54°sin108°-[arg(ri+3,i+4)-arg(ri,i+1)]/2π; (10); (4r/λ)sinφcos54°cos108°-[arg(ri+3,i+4)-arg(ri,i+1)]/2π≤k1+k2≤(4r/λ)sinφcos54