一种利用非圆信号测定天线阵列接收机幅相响应的方法技术

本发明专利技术针对实际应用中天线阵列信号处理算法所使用的天线阵列接收信号模型中存在幅相响应误差的问题，在两个或两个以上的方向已知的非圆信号与方向未知的信号同时存在的情况下，使用非圆信号作为校正源，利用天线阵列的接收信号向量的扩展向量的样本自相关矩阵的噪声子空间与非圆信号的扩展方向向量之间的正交关系来估计天线阵列的幅相响应，从而实现为波达方向估计、波束形成等天线阵列信号处理所使用的天线阵列接收信号模型提供准确的接收机幅相响应估计的目的。

【技术实现步骤摘要】
一种利用非圆信号测定天线阵列接收机幅相响应的方法
本专利技术属于电子信息
中的天线阵列接收机幅相响应的测定方法，特别是在两个或两个以上的方向已知的非圆信号与方向未知的信号同时存在的情况下利用信号的非圆特性测定天线阵列接收机幅相响应的方法。
技术介绍
利用天线阵列接收信号进行信息获取与探测的技术已广泛应用于现代电子侦察、雷达、通信、声纳、地震、射电天文等诸多领域。天线阵列信号处理算法一般都假设所使用的天线阵列接收信号模型中各天线位置、各接收机的幅度和相位响应(简称幅相响应)等模型参数是精确己知的。然而，在当今的加工工艺水平和实际的工程应用中各种误差不可避免，环境的温度、湿度、天线阵列平台的震动、有源器件的老化等都会引起天线阵列信号处理算法所使用的天线阵列接收信号模型中出现阵元位置误差、接收机幅相响应误差。若在天线阵列信号处理算法所使用的天线阵列接收信号模型中不考虑这些误差因素，天线阵列信号处理算法将出现性能严重恶化甚至失效的情况。由于天线阵列信号处理算法所使用的天线阵列接收信号模型中的误差是高精度、高分辨率的天线阵列信号处理技术走向实用化的一个瓶颈，因此，估计天线阵列接收信号模型中的误差具有重要的实用价值，是天线阵列信号处理技术工程实用化的关键之一。天线阵列接收信号模型中的误差估计技术是伴随着天线阵列信号处理技术同时发展的，常见的天线阵列接收信号模型中的误差估计方法是通过对特定方向的天线阵列方向向量直接进行测量来实现的，需要准确已知辅助信号的方向。由于应用环境的物理变化或天线阵列各个天线、接收机的维护更换等原因，天线阵列接收信号模型中的实际幅相响应也会出现相应变化，若不进行重新估计，则天线阵列信号处理算法所使用的天线阵列接收信号模型中始终存在幅相响应误差，仍然无法避免天线阵列信号处理算法出现性能严重恶化甚至失效的情况。若在进行天线阵列信号处理的同时能实时的估计天线阵列接收信号模型中的实际幅相响应，则可用于对消天线阵列信号处理算法所使用的天线阵列接收信号模型中存在的幅相响应误差，避免天线阵列信号处理算法出现性能严重恶化甚至失效的情况。
技术实现思路
为了在进行天线阵列信号处理的同时能实时的估计天线阵列接收信号模型中的实际幅相响应，本专利技术首先利用信号的非圆特性建立天线阵列的接收信号向量的扩展向量和扩展方向向量，再利用天线阵列接收信号的扩展向量的样本自相关矩阵的噪声子空间与非圆信号的扩展方向向量之间的正交关系估计天线阵列的接收机幅相响应，从而实现为波达方向估计、波束形成等天线阵列信号处理所使用的天线阵列接收信号模型提供准确的接收机幅相响应估计的目的。本专利技术采用的天线阵列的接收信号向量通常表示为：其中x(t)为天线阵列的接收信号向量，向量维数等于天线阵列的天线个数M，t为采样时刻，sk(t)、φk、θk和a(θk)分别表示第k个非圆信号的实数发射信号、相位角、相对于天线阵列的方向和方向θk对应的天线阵列的方向向量，sk(t)是实数而不是复数就是本专利技术要利用的信号的非圆特性，k＝1,2,…,K，K为非圆信号的个数，v(t)为天线阵列的接收机噪声向量，∑表示求和，G为一个对角矩阵，其第m个对角元素G(m,m)表示第m个阵元接收机的幅相响应，在无误差的情况下G等于M阶的单位矩阵。本专利技术方法建立的天线阵列的接收信号向量的扩展向量为其中，表示向量或矩阵的共轭转置，表示向量或矩阵的转置。相应的，本专利技术方法建立的非圆信号的扩展方向向量为其中，φk、θk和a(θk)分别表示第k个非圆信号的相位角、相对于天线阵列的方向和方向θk对应的天线阵列的方向向量，k＝1,2,...,K，K为非圆信号的个数。天线阵列的接收信号向量的扩展向量的样本自相关矩阵为其中t为采样时刻，t＝1,2,…,P，P表示与采样时刻数对应的天线阵列的接收信号向量的个数。本专利技术方法利用的样本自相关矩阵的特征值分解为：其中矩阵Λ是对角矩阵，对角线上的元素为样本自相关的特征值，按降序排列即λ1≥λ2＞λ3≥…≥λ2M，矩阵U是由自相关矩阵的特征向量u1,u2,u3,…,u2M构成的矩阵，与特征值一一对应，UH表示矩阵U的共轭转置。记样本自相关矩阵的噪声子空间为：Q＝[uK+1uK+2…u2M]其中K为信号的个数，可以采用
技术介绍
中常用的大特征值判定方法确定非圆信号的个数K，M为天线阵列的天线个数。若不考虑接收机噪声，由式(1)可知样本自相关矩阵的噪声子空间与非圆信号的扩展方向向量b(θk)之间存在正交关系：QHb(θk)＝0，k＝1,2,…,K即使存在噪声，上述正交关系也近似成立。假设第1和第2个信号为方向已知的非圆信号，则有：其中Q1和Q2分别表示由矩阵Q的上面M行向量和下面M行向量组成的矩阵。注意，这里只是以两个非圆信号的情况为例进行说明，本专利技术方法同样可以用于非圆信号个数大于2的情况。由式(3)和式(4)式可得其中g是M阶向量，第m个元素等于第m个阵元接收机的幅相响应，表示把一个向量转换为对角矩阵的操作。将上述两式联合在一起，得因此，向量g为矩阵D的等于1的特征值所对应的特征向量。因此，可以通过计算矩阵D的等于1的特征值所对应的特征向量的方式确定向量g，进而得到天线阵列幅相响应的估计。本专利技术首先利用信号的非圆特性建立天线阵列的接收信号向量的扩展向量，再利用天线阵列接收信号的扩展向量的样本自相关矩阵的噪声子空间与非圆信号的扩展方向向量之间的正交关系估计天线阵列的接收机幅相响应，从而实现为波达方向估计、波束形成等天线阵列信号处理所使用的天线阵列接收信号模型提供准确的接收机幅相响应估计的目的。本专利技术方法包括以下步骤：步骤1.初始化处理：将天线阵列的天线数(记为M)，天线阵列的接收信号向量的个数(记为P)初始化存入内存；步骤2.确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵：首先经过常规方法处理以确定天线阵列的接收信号向量，其次由天线阵列的接收信号向量生成天线阵列的接收信号向量的扩展向量，然后确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵；步骤3.确定样本自相关矩阵的噪声子空间：对样本自相关矩阵进行特征值分解，确定样本自相关矩阵的噪声子空间；步骤4：确定天线阵列幅相响应的估计：利用样本自相关矩阵的噪声子空间与非圆信号的扩展方向向量之间的正交关系，确定天线阵列幅相响应估计。在步骤2中所述经过常规方法处理以获取天线阵列的接收信号向量，其处理方法为I/Q双通道接收方法或希尔伯特变换处理方法。在步骤2中所述天线阵列的接收信号向量的样本通常表示为：x(t)＝[x1(t)x2(t)…xM(t)]T其中x(t)为天线阵列接收的信号向量，向量维数等于天线阵列的天线个数M，t为采样时刻，t＝1,2,…,P，P表示与采样时刻数对应的天线阵列的接收信号向量的个数，xm(t)表示天线阵列的接收信号向量x(t)的第m个元素，m＝1,2,…,M，表示矩阵或向量的转置。在步骤2中所述由天线阵列的接收信号向量生成扩展向量为：其中，向量x*(t)表示向量x(t)的共轭。在步骤2中所述确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵为：其中，表示扩展向量的样本自相关矩阵，t为采样时刻，每一采样时刻对一个接收信号向量采样，t＝1,2,…,P，P表示与采样时刻数对应的天线阵列的接收本文档来自技高网...

【技术保护点】
一种利用非圆信号测定天线阵列接收机幅相响应的方法，包括： 步骤1.初始化处理：将天线阵列的天线数（记为M），天线阵列的接收信号向量的个数（记为P）初始化存入内存； 步骤2.确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵：首先经过常规方法处理以确定天线阵列的接收信号向量，其次由天线阵列的接收信号向量生成天线阵列的接收信号向量的扩展向量，然后确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵； 步骤3.确定样本自相关矩阵的噪声子空间：对样本自相关矩阵进行特征值分解，确定样本自相关矩阵的噪声子空间； 步骤4.确定天线阵列幅相响应的估计：利用样本自相关矩阵的噪声子空间与非圆信号的扩展方向向量之间的正交关系，确定天线阵列幅相响应估计。

【技术特征摘要】
1.一种利用非圆信号测定天线阵列接收机幅相响应的方法，包括：步骤1.初始化处理：将天线阵列的天线数，天线阵列的接收信号向量的个数初始化存入内存；步骤2.确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵：首先经过常规方法处理以确定天线阵列的接收信号向量，其次由天线阵列的接收信号向量生成天线阵列的接收信号向量的扩展向量，然后确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵；经过常规方法处理以确定天线阵列的接收信号向量，其处理方法为I/Q双通道接收方法或希尔伯特变换处理方法，确定的天线阵列的接收信号向量为：x(t)＝[x1(t)x2(t)…xM(t)]T其中x(t)为天线阵列的接收信号向量，向量维数等于天线阵列的天线个数M，t为采样时刻，t＝1,2,…,P，P表示与采样时刻数对应的天线阵列的接收信号向量的个数，x1(t),x2(t),…,xM(t)分别表示天线阵列的接收信号向量x(t)的第1,2,...,M个元素，表示矩阵或向量的转置；由天线阵列的接收信号向量生成天线阵列的接收信号向量的扩展向量为：其中，向量x*(t)表示向量x(t)的共轭，x(t)为天线阵列的接收信号向量；确定天线阵列的接收信号向量的扩展向量的样本自相关矩阵为：其中，表示扩展向量的样本自相关矩阵，每一采样时刻对一个接收信号向量采样，表示向量或矩阵的共轭转置；步骤3.确定扩展向量的样本自相关矩阵的噪声子空间：对扩展向量的样本自相关矩阵进行特征值分解，确定扩展向量的样本自相关矩阵的噪声子空间；扩展向量的样本自相关矩阵的特征值分解为：其中矩阵Λ是对角矩阵，对角向上的元素分别对应扩展向量的样本自相关矩阵的特征值，按降序排列即λ1≥λ2＞λ3≥…≥λ2M，矩阵U是由扩展向量的样本自相关矩阵的特征向量u1,u2,u3,…,u2M构成的矩阵，与特征值一一对应，表示向量或矩阵的共轭转置；...

【专利技术属性】
技术研发人员：徐保根，万群，万义和，汤四龙，文飞，邹继锋，龚辉，丁学科，周志平，
申请(专利权)人：同方电子科技有限公司，电子科技大学，
类型：发明
国别省市：江西;36