本发明属于阵列信号处理领域，尤其涉及一种阵列信号的接收通道的幅相误差的自校正方法。本发明方法不使用外部辅助信号源，通过在不同的采样时间段，切换阵列的阵元与阵元信号接收通道之间的对应关系，就可接收到不同对应关系下的时域信号，选取两相邻时间段内的信号，利用本论文中的自校正方法，对每个通道的幅度和相位误差进行估计，并将误差的估计值代入最后推导的测向公式中，得到信号的来波方向估计。
1.一种阵列信号的接收通道的幅相误差的自校正方法，其特征在于，具体步骤如下：
S1、由M个阵元组成均匀线阵，作为接收阵列，接收阵列中任意两相邻阵元间的距离为d，令接收阵列中的每个阵元对应一个接收通道，对接收阵列中的阵元进行编号，阵元依次编号为c1,c2,...,cM，在t0时刻，接收阵列的阵元与接收通道的对应连接关系变为：编号为c1的阵元连接到的接收通道记为d1，编号为c2的阵元连接到的接收通道记为d2，一直到编号为cM的阵元连接到的接收通道记为dM，在此连接关系下，对每个接收通道的输出采集L个离散时刻的信号，得到一组采样信号Y1＝[y1(t1) y1(t2) … y1(tL)]＝[y11 y12 … y1L]，其中，y1i是阵列接收到的第i个采样时刻的信号向量，y1(tn)为接收阵列在tn时刻的采样信号，s(tn)为tn时刻的入射信号向量，s(tn)＝[s1(tn),s2(tn),...,sK(tn)]T，si(tn)为第i个入射信号，K为入射信号个数，[·]T表示转置操作， 是在该采样时刻，存在接收通道的幅相误差时的信号的方向矩阵，A为理想情况下的信号的方向矩阵，A＝[a(θ1) a(θ2) ... a(θK)]，a(θi)为理想情况下第i个入射信号的方向向量， θi为第i个入射信号的入射角度， 为接收阵列对第i个入射信号进行接收时，由波程差引起的任意两相邻阵元的接收信号之间相位差， λi为第i个入射信号的波长，G1为接收通道的幅相误差矩阵，G1是对角矩阵，G1＝diag(g1,g2,...gM)，diag(·)为对角矩阵符号，g1,g2,...gM为G1对角线上的元素，其余元素为0，gi是编号为ci的阵元连接的接收通道的幅相误差参数，gi＝αiexp(jφi),i＝1,2,...M，αi为该接收通道的幅度误差，φi为该接收通道的相位误差，v(tn)是tn时刻的白噪声向量，其噪声方差为σ2，所述理想情况指不存在接收通道的幅相误差的情况；
S2、在t0+ΔT时刻，切换各阵元对应的接收通道，接收阵列的阵元与接收通道的对应连接关系变为：编号为c1的阵元连接到的接收通道记为d'1，编号为c2的阵元连接到的接收通道记为d'2，一直到编号为cM的阵元连接到的接收通道记为d'M，在此连接关系下，同样对每个接收通道的输出采集L个离散时刻的信号，得到另一组采样信号Y2＝[y2(t1+ΔT) y2(t2+ΔT) … y2(tL+ΔT)]＝[y21 y22 … y2L]，其中，y2i是阵列接收到的第i个采样时刻的信号向量，是在该采样时刻，存在接收通道的幅相误差时的信号的方向矩阵，
G2＝diag(g′1,g'2,g'3,...,g'M)，G2是对角矩阵，g′1,g'2,g'3,...,g'M为G2对角线上的元素，其余元素为0，g′i为编号为ci的阵元连接的接收通道的幅相误差参数，g′i＝α′iexp(jφ′i),i＝1,2,...M，α′i为该接收通道的幅度误差，φ′i为该接收通道的相位误差；
S3、由S1和S2中得到的两组采样信号Y1＝[y11 y12 … y1L]和Y2＝[y21 y22 … y2L]，分别计算两组信号的自相关矩阵的估计值，计算公式为： 和S4、对两个自相关矩阵的估计值R1和R2分别进行特征值分解，可得R1的M个特征值，按降序排列为：λ1≥λ2≥…≥λK＞λK+1＝…＝λM＝σ2，各特征值对应的归一化特征向量依次为：u1,u2,…,uM，同理，可得R2的M个特征值，按降序排列为：λ′1≥λ′2≥…≥λ'K＞λ'K+1＝…＝λ'M＝σ2，各特征值对应的归一化特征向量依次为：u′1,u′2,…,u′M；
S5、构造两组采样信号的Y1和Y2的噪声子空间EN1和EN2，其中，EN1是由λK+1,…,λM对应的归一化特征向量uK+1,…,uM张成的噪声子空间，EN1＝span{uK+1,…,uM}；EN2是由λ'K+1,…,λ'M对应的归一化特征向量u'K+1,…,u'M张成的噪声子空间，EN2＝span{u'K+1,…,u'M}；
S6、令S1中的 a'1(θi)是在t0时刻，存在接收通道的幅相
同理，令S2中的 a'2(θi)是在t0+ΔT时刻，存在接收通道的幅
根据MUSIC算法原理，S5所述噪声子空间EN1的归一化特征向量uK+1,…,uM与存在幅相误差时的信号方向向量正交，即 令U1＝[u K +1 ,… ,u M] ，U2 ＝ [u 'K + 1 ,… ,u 'M] ，可得 正 交关 系 方程 组 ：其中，O是M-K行1列的全零矩阵；
S7、令G1＝F1g G2＝F2g，其中g＝[g1,g2,g3,...,gM]T，F1和F2是表征阵元与接收通道之间的连接关系的M阶满秩矩阵，每一行每一列都仅有1个元素为1，其余元素均为0，F1的对角线上元素均为1，其余元素均为0，F2中的元素，以第一行的元素为例，在t0+ΔT时刻，若编号为c1的阵元连接的接收通道是编号为ck的阵元在t0时刻连接的接收通道，则接收通道的误差g′1＝gk，相应的F2中的第一行的第k列元素为1，第一行的其余元素为0，F2中其余的第2到M行元素中元素1所在的列号，按照同样的方法由编号为c2到cM的阵元连接的接收通道确定；
S8、将S7中得到的关系式G1＝F1g，G2＝F2g代入S6中的正交关系方程组中，可得如下正交关系方程组的等价变换： 令g中的第一个元素g(1)＝1，继续改写方程： 令 b＝[O' 1]T，
且Q为列满秩矩阵，使得方程组Qg＝b存在最小二乘解gLS，所述最小二乘解的表达式为：
其中， Φ(·)表示将括号中的向量变换成对角矩阵的操作，对
角矩阵的对角线上的元素为括号中向量的元素，[·]为共轭转置操作， 为矩阵的伪逆操作；
S9、根据S8中的结果，可得如下的MUSIC谱函数：
2.根据权利要求1所述的一种阵列信号的接收通道的幅相误差的自校正方法，其特征在于：S2所述切换通道的条件为：阵元所在矩阵为列满秩矩阵。
一种阵列信号的接收通道的幅相误差的自校正方法
技术领域
[0001] 本发明属于阵列信号处理领域，尤其涉及一种阵列信号的接收通道的幅相误差的自校正方法。
背景技术
[0002] 阵列信号处理是信号处理领域内的一个重要分支，在近30年来得到迅速发展，其应用涉及雷达、通信、声呐、地震、勘探、射电天文以及生物医学工程等众多军事及国民经济领域。与传统的单个定向传感器相比，阵列信号处理具有灵活的波束控制、高的信号增益、较强的抗干扰能力及高的空间超分辨能力等优点，因而受到了极大关注，与此相关的研究工作不断发展与深入，其应用范围也不断扩大。同时随着微电子技术、数字信号处理技术、并行处理技术的迅猛发展，阵列信号处理的理论和实际应用也得到迅速发展。
[0003] 阵列信号处理最主要的两个研究方向是自适应空域滤波(自适应阵列处理)和空间谱估计。空间谱估计是近30年来发展起来的一门新兴的空域信号处理技术，它是在波束形成技术以及时域谱估计技术基础上发展起来的，其主要目标是研究提高在处理带宽内空间信号角度的估计精度、角度分辨力和提高运算速度的各种算法。传统的波达方向估计方法，阵列的角度分辨力受瑞利限的约束，即瑞利限内的空间目标是不可分辨的。空间谱估计技术把超分辨估计技术应用到DOA估计中，可以突破瑞利限的约束，得到很高的角度分辨率，它在雷达、通信、声呐等众多领域有极为广阔的应用前景。
[0004] 经过数十年的密集研究，空间谱估计测向算法得到了很大的发展，出现了许多较为成熟的算法以及这些算法的改进算法，因此在实际中得到了广泛的应用。
[0005] 在理想条件下，空间谱估计测向算法可获得非常优良的测向性能，但在实际应用中，由阵列信号的各接收通道失配等非理想因素引起的通道的幅相误差往往难以避免，其对测向算法的影响在设计过程中是必须考虑的问题，因而在实际系统中必须对通道的幅相误差进行校正。
发明内容
[0006] 本发明的目的是针对阵列测向方法中，阵列中各阵元信号的接收通道之间存在幅相不一致性，导致阵列信号通过不同的通道时，会有不同的幅度和相位偏差，此即为通道的幅相误差。本专利对此种误差提出了一种自校正方法。该方法不使用外部辅助信号源，通过在不同的采样时间段，切换阵列的阵元与阵元信号接收通道之间的对应关系，就可接收到不同对应关系下的时域信号，选取两相邻时间段内的信号，利用本论文中的自校正方法，对每个通道的幅度和相位误差进行估计，并将误差的估计值代入最后推导的测向公式中，得到信号的来波方向估计。
[0007] 一种阵列信号的接收通道的幅相误差的自校正方法的具体步骤如下：
[0008] S1、由M个阵元组成均匀线阵，作为接收阵列，接收阵列中任意两相邻阵元间的距离为d，令接收阵列中的每个阵元对应一个接收通道，对接收阵列中的阵元进行编号,阵元依次编号为c1,c2,...,cM，在t0时刻，接收阵列的阵元与接收通道的对应连接关系变为：编号为c1的阵元连接到的接收通道记为d1，编号为c2的阵元连接到的接收通道记为d2，一直到编号为cM的阵元连接到的接收通道记为dM，在此连接关系下，对每个接收通道的输出采集L个离散时刻的信号，得到一组采样信号Y1＝[y1(t1) y1(t2) … y1(tL)]＝[y11 y12 … y1L]，其中，y1i是阵列接收到的第i个采样时刻的信号向量，y1(tn)为接收阵列在tn时刻的采样信号， s(tn)为tn时刻的入射信号向量，s(tn)＝[s1(tn),s2(tn),...,sK(tn)]T，si(tn)为第i个入射信号，K为入射信号个数，[·]T表示转置操作， 是在该采样时刻，存在接收通道的幅相误差时的信号的方向矩阵，A为理想情况下的信号的方向矩阵，A＝[a(θ1) a(θ2) ... a(θK)]，a(θi)为理想情况下第i个入射信号的方向向量， θi为第i个入
射信号的入射角度， 为接收阵列对第i个入射信号进行接收时，由波程差引起的任意两相邻阵元的接收信号之间相位差， λi为第i个入射信号的波长，G1为接收通道
的幅相误差矩阵，G1是对角矩阵，G1＝diag(g1,g2,...gM)，diag(·)为对角矩阵符号，g1,g2,...gM为G1对角线上的元素，其余元素为0，gi是编号为ci的阵元连接的接收通道的幅相误差参数，gi＝αiexp(jφi),i＝1,2,...M，αi为该接收通道的幅度误差，φi为该接收通道的相位误差，v(tn)是tn时刻的白噪声向量，其噪声方差为σ2，所述理想情况指不存在接收通道的幅相误差的情况；
[0009] S2、在t0+ΔT时刻，切换各阵元对应的接收通道，接收阵列的阵元与接收通道的对应连接关系变为：编号为c1的阵元连接到的接收通道记为d'1，编号为c2的阵元连接到的接收通道记为d'2，一直到编号为cM的阵元连接到的接收通道记为d'M，在此连接关系下，同样对每个接收通道的输出采集L个离散时刻的信号，得到另一组采样信号Y2＝[y2(t1+ΔT) y2(t2+ΔT) … y2(tL+ΔT)]＝[y21 y22 … y2L]，其中，y2i是阵列接收到的第i个采样时刻的信号向量，是在该采样时刻，存在接收通道的幅相误差时的信号的方向矩阵，G2＝diag
(g′1,g'2,g'3,...,g'M)，G2是对角矩阵，g′1,g'2,g'3,...,g'M为G2对角线上的元素，其余元素为0，g′i为编号为ci的阵元连接的接收通道的幅相误差参数，g′i＝α′iexp(jφ′i),i＝1,
2,...M，α′i为该接收通道的幅度误差，φ′i为该接收通道的相位误差；
[0010] S3、由S1和S2中得到的两组采样信号Y1＝[y11 y12 … y1L]和Y2＝[y21 y22 … y2L]，分别计算两组信号的自相关矩阵的估计值，计算公式为： 和
[0011] S4、对两个自相关矩阵的估计值R1和R2分别进行特征值分解，可得R1的M个特征值，按降序排列为：λ1≥λ2≥…≥λK＞λK+1＝…＝λM＝σ2，各特征值对应的归一化特征向量依次为：u1,u2,…,uM，同理，可得R2的M个特征值，按降序排列为：λ′1≥λ′2≥…≥λ'K＞λ'K+1＝…＝λ'M＝σ2，各特征值对应的归一化特征向量依次为：u′1,u′2,…,u′M；
[0012] S5、构造两组采样信号的Y1和Y2的噪声子空间EN1和EN2，其中，EN1是由λK+1,…,λM对应的归一化特征向量uK+1,…,uM张成的噪声子空间，EN1＝span{uK+1,…,uM}。EN2是由λ'K+1,…,λ'M对应的归一化特征向量u'K+1,…,u'M张成的噪声子空间，EN2＝span{u'K+1,…,u'M}；
[0013] S6、令S1中的 a'1(θi)是在t0时刻，存在接收通道的幅相误差时，第i个入射信号的方向向量，
同理，令S2中的 a'2(θi)是在t0+ΔT时刻，存在接收通道的幅相
误差时，第i个入射信号的方向向量，
根据MUSIC算法原理，S5所述噪声子空间EN1的归一化特征向量uK+1,…,uM与存在幅相误差时的信号方向向量正交，即 令U1＝
[uK+1,…,uM]，U2＝[u'K+1,…,u'M]，可得正交关系方程组：
其中，O是M-K行1列的全零矩阵
[0014] S7、令G1＝F1g G2＝F2g，其中g＝[g1,g2,g3,...,gM]T，F1和F2是表征阵元与接收通道之间的连接关系的M阶满秩矩阵，每一行每一列都仅有1个元素为1，其余元素均为0，F1的对角线上元素均为1，其余元素均为0，F2中的元素，以第一行的元素为例，在t0+ΔT时刻，若编号为c1的阵元连接的接收通道是编号为ck的阵元在t0时刻连接的接收通道，则接收通道的误差g′1＝gk，相应的F2中的第一行的第k列元素为1，第一行的其余元素为0，F2中其余的第2到M行元素中元素1所在的列号，按照同样的方法由编号为c2到cM的阵元连接的接收通道确定；
[0015] S8、将S7中得到的关系式G1＝F1g，G2＝F2g代入S6中的正交关系方程组中，可得如下正交关系方程组的等价变换： 令g中的第一个元素g(1)＝1，继续改写方程： 令 b＝
[O' 1]T，且Q为列满秩矩阵，使得方程组Qg＝b存在最小二乘解gLS，所述最小二乘解的表达式为： 其中， Φ(·)表示将括号中的向量变换成对角矩阵的操
作，对角矩阵的对角线上的元素为括号中向量的元素，[·]H为共轭转置操作， 为矩阵的伪逆操作；
[0016] S9、根据S8中的结果，可得如下的MUSIC谱函数：
[0017] 进一步地，S2所述切换通道的条件为：阵元所在矩阵为列满秩矩阵。
[0018] 本发明的有益效果是：
[0019] 本发明不使用辅助信号源，就可完成对通道幅相误差的校正，得到有效的来波方向估计角。因此可有效降低设备的投入成本，促进空间谱估计方法在实际工程中的推广应用。本发明方法可广泛应用于通信、雷达、声纳、遥测、无线电监测与管理等领域。
附图说明
[0020] 图1本校正方法的MUSIC谱估计结果图
[0021] 图2本校正方法的角度估计结果的对比图
[0022] 图3本校正方法的角度估计结果的均方根误差的对比图
具体实施方式
[0023] 下面结合实施例和附图，详细说明本发明的技术方案。
[0024] 本实施方式采用9元均匀线阵，且有9个接收通道，每个阵元依次从左到右编号为c1,c2,...,c9，接收通道从左到右依次编号为d1,d2,...,d9，阵元间距为λ4，同时存在9个阵元接收通道，两个远场入射信号，入射角分别为4.3°和14.4°，信噪比SNR＝9dB,阵元接收信2
号的噪声是零均值的高斯白噪声，设噪声功率为σ＝1。
[0025] 两种连接关系，分别为：
[0026] 第一种连接关系：阵元编号c1,c2,...,c9分别连接到通道d1,d2,...,d9；
[0027] 第二种连接关系：阵元编号c1,c2,...,c9分别连接到通道d8,d6,d4,d2,d9,d1,d3,d5,d7
[0028] 每种连接关系采样一组序列，每组序列含有64个采样点。第一种连接关系的采样序列为第一组采样序列，第二种连接关系的采样序列为第二组采样信号。
[0029] 此种切换满足S8所述切换条件，矩阵Q的秩为9，为列满秩矩阵。
[0030] 本发明的实施方式的具体流程如下：
[0031] 步骤1.由上述两种连接关系，可获得的两组采样序列分别为：
[0032] 第一组：采样时刻n＝1,2,...,64，其中前3个时刻及后3个时刻的采样序列[0033]
[0034]
[0035] 第二组：采样时刻n＝65,66,...,128，其中前3个时刻及后3个时刻的采样序列[0036]
[0037]
[0038] 步骤2.由步骤1中的两组采样序列，分别计算其自相关矩阵的估计，第一组采样序列的自相关矩阵的估计R1为：矩阵为9行9列，其中前3列及后3列为：
[0039]
[0040]
[0041] 由第一组采样序列的自相关矩阵的估计R1，求得其噪声子空间矩阵U1＝[uK+1,…,uM]为：矩阵为9行7列，其中前3列及最后一列为：
[0042]
[0043]
[0044] 第一种连接关系下，存在接收通道的幅相误差时，信号的方向矩阵为：
[0045]
[0046] 第二组采样序列的自相关矩阵的估计R2为：矩阵为9行9列，其中前3列及后3列为：
[0047]
[0048]
[0049] 由第二组采样序列的自相关矩阵的估计R2，求得其噪声子空间矩阵U2＝[u'K+1,…,u'M]为：矩阵为9行7列，其中前3列及最后1列为：
[0050]
[0051]
[0052] 第二种连接关系下，存在接收通道的幅相误差时，信号的方向矩阵为：
[0053]
[0054] 步骤3.根据技术方案S9中得到的MUSIC谱估计公式：
[0055]
[0056] 最后的角度估计结果为4.5000度和14.000度，与真实入射角度4.3度和14.4度相比，误差分别为0.2度和0.4度，仿真结果见图1，由仿真图可以看出，本专利提出的自校正算法可实现对幅相误差的有效校正，得到有效的角度估计值。
[0057] 步骤4.改变信噪比，其余仿真条件不变，对比信噪比SNR在5～50dB范围内的自校正算法估计结果以及理想情况下的角度估计结果，仿真结果见图2。由仿真图可以看出，随着信噪比的增大，自校正算法得到的角度估计结果逐渐趋近理想情况下的角度估计结果。
[0058] 步骤5.改变信噪比，其余仿真条件不变，对比信噪比SNR在5～50dB范围内的自校正算法估计结果的均方根误差以及理想情况下的角度估计结果的均方根误差，仿真结果见图3。由仿真图可以看出，随着信噪比的增大，自校正算法的角度估计结果的均方根误差逐渐趋于0，证明了本专利提出的自校正算法的有效性。
[0059] 当接收阵列的接收通道存在幅相不一致性时，MUSIC谱估计法无法得到有效的角度估计结果，必须对接收通道的幅相不一致进行校正。本发明的自校正方法不借助辅助信号源，既能保证对接收通道的幅相误差的校正，而且还能有效降低设备的复杂性以及成本，促进了空间谱估计方法在实际工程中的应用。
