一种极化敏感阵列参数估计的降维music算法
一种极化敏感阵列参数估计的降维music算法
【技术领域】
[0001] 本发明属于阵列信号处理技术领域，特别涉及极化敏感阵列信号参数估计算法的 改进，具体的讲是减小算法运算量，提高运算效率的降维MUSIC算法。
【背景技术】
[0002] 极化敏感阵列不仅能够通过阵元的空间位置的差异获取信号的空间到达角信息， 又能够通过阵元正交相位的差异获取信号本身所具有的极化信息，其阵元输出的是含有极 化信息的矢量信息而不是通常的标量信息，这一特性奠定了进一步提高极化敏感阵列的性 能的基础，同时也使得我们有可能进一步提高阵列信号的处理能力。相比于普通标量阵列， 极化敏感阵列存在的性能优势，可以从以下四个方面表现出来：（1)较强的抗干扰能力；
[2] 稳健的检测能力；（3)较高的分辨能力；(4)极化多址的能力。
[0003] 信号参数估计问题是阵列信号处理的一个基本问题，自20世纪60年代以来，D0A 估计技术取得飞速的发展，涌现了大量的优良算法。首先是Burg在1967年提出的最大熵谱 估计方法和Schmidt R 0在1969年提出的最小方差谱估计方法，这两种算法标志着阵列 天线在分辨能力方面取得改善。多重信号分类（Multiple Signal Classification, MUSIC) 算法的提出，使得高分辨测向技术发展有了飞跃性地突破。MUSIC算法实质上是利用噪声子 空间和信号子空间的正交原理而得到的。MUSIC算法自从被提出来，就得到了广泛的研宄 和应用，是一种可产生渐进无偏估计的高分辨率特征结构法。
[0004] 将多重信号分类算法用于极化敏感阵列的信号源到达角和极化参数估计的思想 主要是通过对接收数据进行特征分解或奇异值分解等数学分解，将接收数据划分为两个相 互正交的子空间，即信号子空间与噪声子空间，然后利用阵列流型矢量和子空间的关系来 构造联合谱，从而根据谱峰得到信号源到达角和极化参数估计。传统MUSIC联合谱估计方 法如下：
[0005] 1)先求出阵列输出协方差矩阵；
[0006] 2)对协方差矩阵进行特征值分解，得到噪声子空间Un;
[0007] 3)由噪声子空间与阵列信号矢量>!(久叭;V7)张成的子空间与噪声子空间Un正交 的原理，构造MUSIC联合谱估计：
[0009] 由式（1)可以看出，要估计信号源的到达角和极化参数需要进行四维谱峰搜索， 需要付出巨大的运算量和储存量，不易实现。
【发明内容】
[0010] 本发明提供了一种极化敏感阵列参数估计的降维MUSIC算法，在没有降低算法估 计误差性能的同时，将四维谱峰搜索降低为二维谱峰搜索，使得算法的计算量级从〇{n4}降 为〇 In2} (n为搜索范围内点数），大大降低了算法的运算复杂度，提高了系统的运算效率。 [0011] 本发明解决上述问题采用的技术方案为：
[0012] 一种极化敏感阵列参数估计的降维MUSIC算法，其特征在于包括以下步骤：
[0013] 步骤1.建立阵列信号的数学模型
[0014] 空间有K个不相关的信号入射到阵列中，根据入射信号的到达角的供)以及极化参 数（Y，n)和噪声得到阵列接收信号模型的表达式
[0016]其中，0 G[-31/2, JT/2)为入射信号俯仰角，沪e[0,2;r)为入射信号方位角， y G [0, 31/2)为极化辅角，n e [-31，31)为极化相位差。整个阵列的流形矩阵A = [a" a2, ? ? ?，aK]，s(n) = h (n), s2(n), ? ? ?，sK(n)] T为信号矢量，v(n)为噪声。阵列扫描矢 量ak是极化导向矢量a kp和空间导向矢量a kE^ Kronecker积，即
[0018] 其极化导向矢量表示为：
[0020] 其中，匕@，0表示球坐标系与直角坐标系单位矢量之间的变换关系矩阵， E k(y, n)表示完全极化波的Jones矢量。
[0021] 需注意的是，步骤1中根据需求的不同的六个分量可以不用全部使 用。
[0022] 步骤2.利用步骤1得到的阵列接收的N次快拍数据，得到估计信号相关矩阵
[0024]步骤3.通过对相关矩阵R进行特征值分解，得到噪声子空间Un
[0026] 式中，P为取~(沒承；^7)的分量个数，M为阵元个数。A =iag{>1，人2，…，入PM}， 并且最小的PM-K个特征值相等，均等于阵列接收噪声强度。K个大的特征值对应的特征矢 量构成信号子空间<S> = span{Us}，PM-K个小特征值对应的特征矢量构成噪声子空间<N> =span{U n}。阵列信号矢量张成子空间和信号子空间相同，并且与噪声子空间正交，表不 为：
[0027] span {A} = span {Uj (7)
[0028] span {A}丄 span {Un} (8)
[0029] 步骤4.计算信号源的俯仰角和方位角过程
[0030] 降维MUSIC算法将MUSIC谱估计中波达方向角和极化参数进行剥离，根据传统的 MUSIC谱估计
，定义函数
[0034] 由子空间原理，阵列流形矢量张成的子空间与噪声子空间正交，即：
[0035] span {A}丄 span {Un} (11)
[0036] 将式（11)代入式（9)得到
[0038]当y G (〇,jt/2)时，由式⑷可知EH(y，q)是列满秩的，当为非满秩， 即del0(久4 = 0，式（12)成立，。则俯仰角和方位角的估计
[0040]步骤5.计算信号源的极化参数过程
[0041] 将所求俯仰角和方位角代入中间变量，得到
[0043] 则极化参数估计：
[0044]
[0045] 由式（13)和式（15)可以看出，本发明将信号源到达角和极化参数分离开来搜索， 将一个四维谱峰搜索变为两个二维谱峰搜索，通过本发明的方法，估计入射信号角度和极 化参数的计算量大大降低，且具有良好的正确性和精确度。
【附图说明】
[0046] 图1为由双正交电偶极子对构成的均匀线阵。
【具体实施方式】
[0047] 如图1所示由M个双正交电偶极子对沿y轴排列构成均匀线阵，双正交电偶极子 对分别沿x轴和y轴放置，阵元间距为d，本例中降维MUSIC算法的具体实现如下：
[0048] 步骤1?建立阵列信号数学模型：
[0049] 已知空间有K个信号入射到极化敏感阵列中，入射窄带信号互不相关， 0 e [-31/2, Ji/2)为入射信号俯仰角，pe[CU；r)为入射信号方位角，y e [0,jt/2)为极 化辅角，n e [-31，31)为极化相位差，阵列的接收信号模 型x(n)可以表示为：
[0051]式中A = [a^ a2, . . . aK]为整个阵列的流形矩阵，其中ak(k = 1，2. . . K)为信号极 化域-空域联合导向矢量；信号矢量s(n) = [8101)，82(11)，...，％(11)]\奴11)为均值为零， 方差为〇2的复高斯白噪声。ak是极化导向矢量（akp)和空间导向矢量（aks)的Kronecker 积，即：
[0053]其中，
[0055] 式中，
A表示波长。
[0056] 对于沿x轴和y轴放置的电偶极子对，其极化导向矢量akp只使用两个分量，表不 为：
[0058]将式（19)代入式（17)得到
[0060]步骤2.求信号的自相关矩阵[0061]阵列接收信号的自相关矩阵为：
[0062] Rx= E [x (n) x H (n) ] = ARSAH+ 〇 2I (21)
[0063] 其中，A为整个阵列的流形矩阵，Rs为入射信号的自相关函数，I为单位阵。
[0064] 实际工程中，通常采用N个观测样本值x (0)，x (1)，. . .，x (N-1)得到相关矩阵的估 计；，用义代替艮进行计算，其中：
[0066]步骤3?求解噪声子空间：
[0067]将步骤3求得的相关矩阵的估计^进行特征值分解，求得噪声子空间：
[0069] 其中A 3和A n为特征值矩阵，分别对应于K个大特征值和2M-K个小特征值。U s 为K个大特征值所对应的特征矢量，其张成的空间为信号子空间，队为2M-K个小特征值所 对应的特征矢量，其张成的空间为噪声子空间。
[0070] 步骤4.计算信号源的俯仰角和方位角
[0071] 首先定义检测量为：
[0073]将步骤1中公式（20)代入式（24)得到
[0075]由于
诚立，则式（25)可以表示为：
[0077]由式（26)可以看出入射信号的波达方向角信息和极化参数信息被分离开来，此 时定义中间变量：
[0079] 由子空间原理，信号导向矢量张成的子空间与噪声子空间的特征向量正交
[0080] span {A}丄 span {Un} (28)
[0081] 将式（28)代入式（26)可以得到
[0083]当yG (〇,jt/2)时EH(y，q)是列满秩的，要使得
非 满秩，即
。根据最优化搜索方法得到估计零谱的K个极小值，这些极小值的 位置即为信号源到达角的估计。为了方便起见，通常用零谱的倒数来搜索其K个极 大值。定义降维MUSIC谱：
[0085] 其中K个谱峰对应的坐标即为信号源的到达角估计
[0086] 步骤5?计算信号源的极化参数：
[0087]将步骤4得到的信号源到达角代入式（27)中，得到中间变量的估计
[0089] 对于极化参数的估计，可以由如下的谱估计搜索得到：
[0091] 对于极化敏感阵列的信号源到达角和极化参数估计，传统的联合MUSIC算法需要 对四个参数同时进行谱峰搜索，其计算量为〇 (n4)，其中n为搜索范围内点数，而本发明的算 法将信号源的到达角和极化参数进行分离搜索，其计算量为〇(2n2)，由此可见大大降低了 运算量，提高了系统的运算效率。
【主权项】
1. 一种极化敏感阵列参数估计的降维MUSIC算法，其特征在于包括以下步骤： 步骤1.建立阵列信号的数学模型 空间有K个不相关的信号入射到阵列中，根据入射信号的到达角以及极化参数 (γ，η)和噪声得到阵列接收信号模型的表达式：(2) 其中，Θ e [- π /2, π /2)为入射信号俯仰角，沪e [0,2;r)为入射信号方位角， γ e [〇, jt/2)为极化辅角，n e [_π, π)为极化相位差；整个阵列的流形矩阵A = Iia1, a2, · · ·，aK]，s(n) = [S1Oi), s2(n), · · ·，sK(n)] τ为信号矢量，ν(η)为噪声；阵列扫描矢 量ak是极化导向矢量a kp和空间导向矢量a kE^ Kronecker积，即： α,(θ,φ,γ,η) = α^θ.,φ)? α!φ(β..φ.ν^η) (3) 其极化导向矢量akp表示为：其中，6(0#)表示球坐标系与直角坐标系单位矢量之间的变换关系矩阵，Ek(Y，η) 表示完全极化波的Jones矢量； 步骤2.利用步骤1得到的阵列接收的N次快拍数据，得到估计信号相关矩阵：(5) 步骤3.通过对相关矩阵R进行特征值分解，得到噪声子空间Un(6) 式中，P为取iV:艮叭;"）的分量个数，M为阵元个数。Λ =CliagU1, λ 2，...，λ PM}，并 且最小的PM-K个特征值相等，均等于阵列接收噪声强度；K个大的特征值对应的特征矢量 构成信号子空间<S> = span {Us}，PM-K个小特征值对应的特征矢量构成噪声子空间<N> = span{Un};阵列信号矢量张成子空间和信号子空间相同，并且与噪声子空间正交，表不为： span {A} = span {Uj (7) span {A}丄 span {Un} (8) 步骤4.计算信号源的俯仰角和方位角过程 将MUSIC谱估计中波达方向角和极化参数进行剥离，根据传统的MUSIC谱估计定义函数：由子空间原理，阵列流形矢量张成的子空间与噪声子空间正交，即： span {A} _L span {Un} (11) 将式（11)代入式（9)得到(12) 当γ e (〇, π/2)时，由式⑷知ΕΗ(γ，τι)是列满秩的；当为非满秩，即 式（12)成立； 则俯仰角和方位角的估计：(13) 步骤5.计算信号源的极化参数过程 将所求俯仰角和方位角(^，^|代入中间变量，得到(14)则极化 参数估计：(卜）
【专利摘要】本发明公开了一种极化敏感阵列参数估计的降维MUSIC算法，属于阵列信号处理技术领域。本发明将阵列接收数据得到的信号相关矩阵进行特征值分解得到噪声子空间；将信号的导向矢量中波达方向角和极化信息分离开来，得到中间变量；由中间变量构造波达方向角的搜索函数，谱峰搜索得到波达方向角的估计；将波达方向角的估计值代入降维MUSIC谱估计，由谱峰对应的坐标得到极化参数估计。本发明的算法将四维谱峰搜索降低为二维谱峰搜索，使得算法的计算量级从o{n4}降为o{n2}，大大降低了算法的运算复杂度，提高了系统的运算效率。
【IPC分类】G01S3/14
【公开号】CN104898085
【申请号】CN201510245865
【发明人】李会勇, 张远芳, 周正, 李洋, 谢菊兰
【申请人】电子科技大学
【公开日】2015年9月9日
【申请日】2015年5月14日
【技术领域】
[0001] 本发明属于阵列信号处理技术领域，特别涉及极化敏感阵列信号参数估计算法的 改进，具体的讲是减小算法运算量，提高运算效率的降维MUSIC算法。
【背景技术】
[0002] 极化敏感阵列不仅能够通过阵元的空间位置的差异获取信号的空间到达角信息， 又能够通过阵元正交相位的差异获取信号本身所具有的极化信息，其阵元输出的是含有极 化信息的矢量信息而不是通常的标量信息，这一特性奠定了进一步提高极化敏感阵列的性 能的基础，同时也使得我们有可能进一步提高阵列信号的处理能力。相比于普通标量阵列， 极化敏感阵列存在的性能优势，可以从以下四个方面表现出来：（1)较强的抗干扰能力；
[2] 稳健的检测能力；（3)较高的分辨能力；(4)极化多址的能力。
[0003] 信号参数估计问题是阵列信号处理的一个基本问题，自20世纪60年代以来，D0A 估计技术取得飞速的发展，涌现了大量的优良算法。首先是Burg在1967年提出的最大熵谱 估计方法和Schmidt R 0在1969年提出的最小方差谱估计方法，这两种算法标志着阵列 天线在分辨能力方面取得改善。多重信号分类（Multiple Signal Classification, MUSIC) 算法的提出，使得高分辨测向技术发展有了飞跃性地突破。MUSIC算法实质上是利用噪声子 空间和信号子空间的正交原理而得到的。MUSIC算法自从被提出来，就得到了广泛的研宄 和应用，是一种可产生渐进无偏估计的高分辨率特征结构法。
[0004] 将多重信号分类算法用于极化敏感阵列的信号源到达角和极化参数估计的思想 主要是通过对接收数据进行特征分解或奇异值分解等数学分解，将接收数据划分为两个相 互正交的子空间，即信号子空间与噪声子空间，然后利用阵列流型矢量和子空间的关系来 构造联合谱，从而根据谱峰得到信号源到达角和极化参数估计。传统MUSIC联合谱估计方 法如下：
[0005] 1)先求出阵列输出协方差矩阵；
[0006] 2)对协方差矩阵进行特征值分解，得到噪声子空间Un;
[0007] 3)由噪声子空间与阵列信号矢量>!(久叭;V7)张成的子空间与噪声子空间Un正交 的原理，构造MUSIC联合谱估计：
[0009] 由式（1)可以看出，要估计信号源的到达角和极化参数需要进行四维谱峰搜索， 需要付出巨大的运算量和储存量，不易实现。
【发明内容】
[0010] 本发明提供了一种极化敏感阵列参数估计的降维MUSIC算法，在没有降低算法估 计误差性能的同时，将四维谱峰搜索降低为二维谱峰搜索，使得算法的计算量级从〇{n4}降 为〇 In2} (n为搜索范围内点数），大大降低了算法的运算复杂度，提高了系统的运算效率。 [0011] 本发明解决上述问题采用的技术方案为：
[0012] 一种极化敏感阵列参数估计的降维MUSIC算法，其特征在于包括以下步骤：
[0013] 步骤1.建立阵列信号的数学模型
[0014] 空间有K个不相关的信号入射到阵列中，根据入射信号的到达角的供)以及极化参 数（Y，n)和噪声得到阵列接收信号模型的表达式
[0016]其中，0 G[-31/2, JT/2)为入射信号俯仰角，沪e[0,2;r)为入射信号方位角， y G [0, 31/2)为极化辅角，n e [-31，31)为极化相位差。整个阵列的流形矩阵A = [a" a2, ? ? ?，aK]，s(n) = h (n), s2(n), ? ? ?，sK(n)] T为信号矢量，v(n)为噪声。阵列扫描矢 量ak是极化导向矢量a kp和空间导向矢量a kE^ Kronecker积，即
[0018] 其极化导向矢量表示为：
[0020] 其中，匕@，0表示球坐标系与直角坐标系单位矢量之间的变换关系矩阵， E k(y, n)表示完全极化波的Jones矢量。
[0021] 需注意的是，步骤1中根据需求的不同的六个分量可以不用全部使 用。
[0022] 步骤2.利用步骤1得到的阵列接收的N次快拍数据，得到估计信号相关矩阵
[0024]步骤3.通过对相关矩阵R进行特征值分解，得到噪声子空间Un
[0026] 式中，P为取~(沒承；^7)的分量个数，M为阵元个数。A =iag{>1，人2，…，入PM}， 并且最小的PM-K个特征值相等，均等于阵列接收噪声强度。K个大的特征值对应的特征矢 量构成信号子空间<S> = span{Us}，PM-K个小特征值对应的特征矢量构成噪声子空间<N> =span{U n}。阵列信号矢量张成子空间和信号子空间相同，并且与噪声子空间正交，表不 为：
[0027] span {A} = span {Uj (7)
[0028] span {A}丄 span {Un} (8)
[0029] 步骤4.计算信号源的俯仰角和方位角过程
[0030] 降维MUSIC算法将MUSIC谱估计中波达方向角和极化参数进行剥离，根据传统的 MUSIC谱估计
，定义函数
[0034] 由子空间原理，阵列流形矢量张成的子空间与噪声子空间正交，即：
[0035] span {A}丄 span {Un} (11)
[0036] 将式（11)代入式（9)得到
[0038]当y G (〇,jt/2)时，由式⑷可知EH(y，q)是列满秩的，当为非满秩， 即del0(久4 = 0，式（12)成立，。则俯仰角和方位角的估计
[0040]步骤5.计算信号源的极化参数过程
[0041] 将所求俯仰角和方位角代入中间变量，得到
[0043] 则极化参数估计：
[0044]
[0045] 由式（13)和式（15)可以看出，本发明将信号源到达角和极化参数分离开来搜索， 将一个四维谱峰搜索变为两个二维谱峰搜索，通过本发明的方法，估计入射信号角度和极 化参数的计算量大大降低，且具有良好的正确性和精确度。
【附图说明】
[0046] 图1为由双正交电偶极子对构成的均匀线阵。
【具体实施方式】
[0047] 如图1所示由M个双正交电偶极子对沿y轴排列构成均匀线阵，双正交电偶极子 对分别沿x轴和y轴放置，阵元间距为d，本例中降维MUSIC算法的具体实现如下：
[0048] 步骤1?建立阵列信号数学模型：
[0049] 已知空间有K个信号入射到极化敏感阵列中，入射窄带信号互不相关， 0 e [-31/2, Ji/2)为入射信号俯仰角，pe[CU；r)为入射信号方位角，y e [0,jt/2)为极 化辅角，n e [-31，31)为极化相位差，阵列的接收信号模 型x(n)可以表示为：
[0051]式中A = [a^ a2, . . . aK]为整个阵列的流形矩阵，其中ak(k = 1，2. . . K)为信号极 化域-空域联合导向矢量；信号矢量s(n) = [8101)，82(11)，...，％(11)]\奴11)为均值为零， 方差为〇2的复高斯白噪声。ak是极化导向矢量（akp)和空间导向矢量（aks)的Kronecker 积，即：
[0053]其中，
[0055] 式中，
A表示波长。
[0056] 对于沿x轴和y轴放置的电偶极子对，其极化导向矢量akp只使用两个分量，表不 为：
[0058]将式（19)代入式（17)得到
[0060]步骤2.求信号的自相关矩阵[0061]阵列接收信号的自相关矩阵为：
[0062] Rx= E [x (n) x H (n) ] = ARSAH+ 〇 2I (21)
[0063] 其中，A为整个阵列的流形矩阵，Rs为入射信号的自相关函数，I为单位阵。
[0064] 实际工程中，通常采用N个观测样本值x (0)，x (1)，. . .，x (N-1)得到相关矩阵的估 计；，用义代替艮进行计算，其中：
[0066]步骤3?求解噪声子空间：
[0067]将步骤3求得的相关矩阵的估计^进行特征值分解，求得噪声子空间：
[0069] 其中A 3和A n为特征值矩阵，分别对应于K个大特征值和2M-K个小特征值。U s 为K个大特征值所对应的特征矢量，其张成的空间为信号子空间，队为2M-K个小特征值所 对应的特征矢量，其张成的空间为噪声子空间。
[0070] 步骤4.计算信号源的俯仰角和方位角
[0071] 首先定义检测量为：
[0073]将步骤1中公式（20)代入式（24)得到
[0075]由于
诚立，则式（25)可以表示为：
[0077]由式（26)可以看出入射信号的波达方向角信息和极化参数信息被分离开来，此 时定义中间变量：
[0079] 由子空间原理，信号导向矢量张成的子空间与噪声子空间的特征向量正交
[0080] span {A}丄 span {Un} (28)
[0081] 将式（28)代入式（26)可以得到
[0083]当yG (〇,jt/2)时EH(y，q)是列满秩的，要使得
非 满秩，即
。根据最优化搜索方法得到估计零谱的K个极小值，这些极小值的 位置即为信号源到达角的估计。为了方便起见，通常用零谱的倒数来搜索其K个极 大值。定义降维MUSIC谱：
[0085] 其中K个谱峰对应的坐标即为信号源的到达角估计
[0086] 步骤5?计算信号源的极化参数：
[0087]将步骤4得到的信号源到达角代入式（27)中，得到中间变量的估计
[0089] 对于极化参数的估计，可以由如下的谱估计搜索得到：
[0091] 对于极化敏感阵列的信号源到达角和极化参数估计，传统的联合MUSIC算法需要 对四个参数同时进行谱峰搜索，其计算量为〇 (n4)，其中n为搜索范围内点数，而本发明的算 法将信号源的到达角和极化参数进行分离搜索，其计算量为〇(2n2)，由此可见大大降低了 运算量，提高了系统的运算效率。
【主权项】
1. 一种极化敏感阵列参数估计的降维MUSIC算法，其特征在于包括以下步骤： 步骤1.建立阵列信号的数学模型 空间有K个不相关的信号入射到阵列中，根据入射信号的到达角以及极化参数 (γ，η)和噪声得到阵列接收信号模型的表达式：(2) 其中，Θ e [- π /2, π /2)为入射信号俯仰角，沪e [0,2;r)为入射信号方位角， γ e [〇, jt/2)为极化辅角，n e [_π, π)为极化相位差；整个阵列的流形矩阵A = Iia1, a2, · · ·，aK]，s(n) = [S1Oi), s2(n), · · ·，sK(n)] τ为信号矢量，ν(η)为噪声；阵列扫描矢 量ak是极化导向矢量a kp和空间导向矢量a kE^ Kronecker积，即： α,(θ,φ,γ,η) = α^θ.,φ)? α!φ(β..φ.ν^η) (3) 其极化导向矢量akp表示为：其中，6(0#)表示球坐标系与直角坐标系单位矢量之间的变换关系矩阵，Ek(Y，η) 表示完全极化波的Jones矢量； 步骤2.利用步骤1得到的阵列接收的N次快拍数据，得到估计信号相关矩阵：(5) 步骤3.通过对相关矩阵R进行特征值分解，得到噪声子空间Un(6) 式中，P为取iV:艮叭;"）的分量个数，M为阵元个数。Λ =CliagU1, λ 2，...，λ PM}，并 且最小的PM-K个特征值相等，均等于阵列接收噪声强度；K个大的特征值对应的特征矢量 构成信号子空间<S> = span {Us}，PM-K个小特征值对应的特征矢量构成噪声子空间<N> = span{Un};阵列信号矢量张成子空间和信号子空间相同，并且与噪声子空间正交，表不为： span {A} = span {Uj (7) span {A}丄 span {Un} (8) 步骤4.计算信号源的俯仰角和方位角过程 将MUSIC谱估计中波达方向角和极化参数进行剥离，根据传统的MUSIC谱估计定义函数：由子空间原理，阵列流形矢量张成的子空间与噪声子空间正交，即： span {A} _L span {Un} (11) 将式（11)代入式（9)得到(12) 当γ e (〇, π/2)时，由式⑷知ΕΗ(γ，τι)是列满秩的；当为非满秩，即 式（12)成立； 则俯仰角和方位角的估计：(13) 步骤5.计算信号源的极化参数过程 将所求俯仰角和方位角(^，^|代入中间变量，得到(14)则极化 参数估计：(卜）
【专利摘要】本发明公开了一种极化敏感阵列参数估计的降维MUSIC算法，属于阵列信号处理技术领域。本发明将阵列接收数据得到的信号相关矩阵进行特征值分解得到噪声子空间；将信号的导向矢量中波达方向角和极化信息分离开来，得到中间变量；由中间变量构造波达方向角的搜索函数，谱峰搜索得到波达方向角的估计；将波达方向角的估计值代入降维MUSIC谱估计，由谱峰对应的坐标得到极化参数估计。本发明的算法将四维谱峰搜索降低为二维谱峰搜索，使得算法的计算量级从o{n4}降为o{n2}，大大降低了算法的运算复杂度，提高了系统的运算效率。
【IPC分类】G01S3/14
【公开号】CN104898085
【申请号】CN201510245865
【发明人】李会勇, 张远芳, 周正, 李洋, 谢菊兰
【申请人】电子科技大学
【公开日】2015年9月9日
【申请日】2015年5月14日
最新回复(0)