CN1746697A - 一种芯片可实现的多信号分类算法 - Google Patents

一种芯片可实现的多信号分类算法 Download PDF

Info

Publication number
CN1746697A
CN1746697A CN 200510021874 CN200510021874A CN1746697A CN 1746697 A CN1746697 A CN 1746697A CN 200510021874 CN200510021874 CN 200510021874 CN 200510021874 A CN200510021874 A CN 200510021874A CN 1746697 A CN1746697 A CN 1746697A
Authority
CN
China
Prior art keywords
centerdot
prime
matrix
order
algorithm
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.)
Granted
Application number
CN 200510021874
Other languages
English (en)
Other versions
CN100385249C (zh
Inventor
魏平
李立萍
张乐
肖先赐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CNB2005100218746A priority Critical patent/CN100385249C/zh
Publication of CN1746697A publication Critical patent/CN1746697A/zh
Application granted granted Critical
Publication of CN100385249C publication Critical patent/CN100385249C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

一种芯片可实现的多信号分类算法,属于通信技术领域,具体应用于阵列信号处理中的波达方向估计。包括以下步骤:采集数据,构造数据矩阵;计算协方差矩阵;特征值分解;构造噪声子空间;谱峰搜索。其本质是使用了本发明提出的改进坐标旋转算法来实现顺序Jacobi算法中的反正切运算和坐标旋转,而通常使用CORDIC算法需分别实现顺序Jacobi算法中的反正切运算和坐标旋转;同时改进了顺序Jacobi算法的迭代结构,把原来需要进行两次坐标旋转的计算复杂度降低为只需要进行一次坐标旋转。本发明可以提高算法并行性,大大提高运算速度。在芯片中进行EVD运算的时间是原来运行时间的1/3左右,而进行整个多信号分类算法时间降低为1/2左右。

Description

一种芯片可实现的多信号分类算法
技术领域
一种芯片可实现的多信号分类算法,属于通信技术领域,具体应用于阵列信号处理中的波达方向估计(Direction Of Arrival,以下简称DOA)。
背景技术
随着移动通信的蓬勃发展,用户数量迅速增加,频谱资源越来越紧张,如何利用现有频谱资源进一步扩展容量已成为移动通信发展的关键问题。智能天线技术可以产生多个空间定向波束,动态改变覆盖区域形状,使天线主波束对准用户信号到达方向,旁瓣或零陷对准干扰信号到达方向,从而能够降低系统干扰,提高系统容量和频谱效率。
实现智能天线的前提是首先要确定用户信号到达方向,DOA是其中重要的一种方法。现有的DOA方法有很多,如:多信号分类(MUltiple SIgnal Classification,以下简称MUSIC)算法、ESPRIT算法、最大似然估计算法等。MUSIC算法是一种综合考虑性能和复杂度后得到的一种比较优化的算法。现在将MUSIC算法简洁描述如下:
1.采集数据,构造数据矩阵:
当信号从远场(一般认为信号源与阵列天线距离大于信号源信号的波长10倍以上)入射到阵列天线时候,可以近似认为该信号到各天线阵元的方向角一致。常用的阵列天线是均匀的线阵或圆阵,如图1、2所示。
首先在阵列天线的N个阵元并行采集M个数据,形成数据矩阵A:
A = x 1 ( 1 ) x 1 ( 2 ) · · · x 1 ( j ) · · · x 1 ( M ) · · · · · · · · · · · · · · · · · · x i ( 1 ) x i ( 2 ) · · · x i ( j ) · · · x i ( M ) · · · · · · · · · · · · · · · · · · x N ( 1 ) x N ( 2 ) · · · x N ( j ) · · · x N ( M ) - - - ( 1 )
其中xi(j)中的i代表第i(i≤N)个阵元,j代表第j(j≤N)个采集数据。
2.计算协方差矩阵:
在得到数据矩阵A后,则可以计算协方差矩阵R:
R=AAT                                            (2)
从(2)可确定矩阵R是N×N对称实矩阵。
3.特征值分解(EVD):
对协方差矩阵R进行特征值分解,对R的特征值和特征向量矩阵按照R的特征值大小排序,得到:
1,λ2…λj…λN}其中λ1≥λ2…λj…≥λN≥0,
和对应特征向量矩阵
V = v 11 v 12 · · · v 1 j · · · v 1 N · · · · · · · · · · · · · · · · · · V i 1 v i 2 · · · v ij · · · v iN · · · · · · · · · · · · · · · · · · v N 1 v N 2 · · · v Nj · · · v NN - - - ( 3 )
记矩阵V第j列矢量vj为λj对应的特征向量,则V=[V1 V2…Vj…VN]。
4.构造噪声子空间:
对于D个信号入射的情况,R的特征值有如下性质:
λ k = μ k + σ k k = 1,2 · · · D σ k k = D + 1 , · · · N ;
上式中的μk代表入射信号的特征值(μ1≥μ2…≥μD≥0),σk代表噪声方差。一般说来,λD和λD+1在信噪比不恶劣的情况下,都有数量级上的差距,从而能够很好区分出信号个数D;
估计出D大小,从V中构造噪声子空间EN H,其中 E N H = [ V D + 1 · · · V N ] .
5.谱峰搜索:
阵列天线的方向矩阵由信号入射方向θk(k=1,2…D)、天线阵列几何参数(均匀圆阵、均匀线阵)等决定。设方向矩阵的列矢量是α(θ)(当θ一定时候,该值为常数),算法证明,当θ=θk E N H a ( θ k ) = 0 .
于是构造空间谱函数P(θ),在θ∈[0,180]的范围内进行谱峰搜索:
P ( θ ) = 1 | | E N H a ( θ ) | | 2
‖·‖表示矩阵的模2范数,P(θ)中D个峰值对应的D个θ就是待求的信号入射方向θk
MUSIC算法第2步中的特征值分解的计算量占到了整个MUSIC算法的60%。其中,涉及特征值分解的算法有很多,如Jacobi算法、顺序Jacobi算法、QR算法等等。这些算法的特点无一例外的是运算繁杂,都存在非线性运算,不利于芯片实现。相对来说,顺序Jacobi算法是计算MUSIC算法中第2步的特征值分解最理想算法,并行计算时收敛速度最快,非线性计算相对较少,便于芯片实现。
顺序Jacobi算法求解特征值和特征向量的基本步骤如下:
设待特征值分解矩阵是R,特征向量矩阵V(矩阵V的初值是N×N单位矩阵),
1、从R选择消去对象rpq
按照r12,r13,…r1N;r23,r24…r2N;……;r(N-1)N矩阵元素所在的位置的先后顺序从R选择一个消去对象rpq;当选择r(N-1)N后,下次选择r12并重复上述顺序重新进行选择;
2、计算反正切:
θ opt = 1 2 tan - 1 [ 2 r pq r pp - r qq ] - - - ( 4 )
在确定第1步的rpq后,按照rpq从R中提取rpp和rqq,计算反正切得到旋转角θopt
3、计算坐标旋转:
Figure A20051002187400073
R ′ = W pq T RW pq - - - ( 6 )
         V′=Wpq TV                             (7)
在第2步计算得到旋转角θopt后,按照公式(5)构造Wpq T,然后按照公式(6),对R做2次坐标旋转后得到矩阵R′;在计算R′同时,按照公式(7),对V做1次坐标旋转就得到V′,具体计算时候可以用如下表达式
r pp ′ r pq ′ r qp ′ r qq ′ = cos θ opt - sin θ opt sin θ opt cos θ opt r pp r pq r qp r qq cos θ opt sin θ opt - sin θ opt cos θ opt
r pj ′ r qj ′ = cos θ opt - sin θ opt sin θ opt cos θ opt r pj r qj j ≠ p , q ;
v jp ′ v jq ′ = cos θ opt - sin θ opt sin θ opt cos θ opt v jp v jq j任意值
4.判断R′
如果R′不是三角型矩阵,则令R=R′、V=V′后,返回第1步选择下一个消去对象,然后重复步骤1-4;如果R′已经转化为三角型,此时R′的主对角线元素就是特征值,V′就是特征向量矩阵。
5.排序
对R′中的特征值按从大到小顺序,对R′,V′进行对应调整,使得{λ1,λ2…λj…λN}其中λ1≥λ2…λj…≥λN≥0
V={v1 v2…vj…vN}(vj是λj对应的特征向量)
使用上述步骤实现特征值分解时候,运算量主要集中在第2步的反正切和第3步的坐标旋转,因为采用了串行方式来实现反正切和2次坐标旋转,导致整个实现算法并行性很低,运算时间较长,无法满足智能天线在工程应用中的要求。
工程应用中实现反正切和坐标旋转常用的一种算法是可芯片实现的CORDIC算法,它可以在仅使用加法器和移位器的情况下在一个芯片上就实现这两种非线性运算。在使用CORDIC算法实现这两种运算的时候,其运算时间和运算原理基本一致,下面仅就坐标旋转的原理和步骤进行描述,反正切不再赘述。
CORDIC算法实现坐标旋转的原理可以解释如下:
已知
Figure A20051002187400084
和旋转角θ,需要计算:
y ′ x ′ = cos θ - sin θ sin θ cos θ y x - - - ( 8 )
首先定义CORDIC基本的角度集:
1,α2…αk}={arctg2-1,arctg2-2…arctg2-k}(k根据结果精度选择)和CORDIC符号集:
1,ξ2…ξk,ξi=1或-1,i=1,2…k}(k根据结果精度选择),
CORDIC算法实现过程就是确定符号集{ξ1,ξ2…ξk,ξi=1或-1,i=1,2…k},使之满足: θ ≈ Σ i = 1 k ξ i α i , 从而把坐标
Figure A20051002187400092
旋转角θ近似等效为坐标
Figure A20051002187400093
旋转1组CORDIC基本的角{ξ1α1,ξ2α2…ξkαk},所以式(8)可以近似等效
y ′ x ′ = cos ξ k α k - sin ξ k α k sin ξ k α k cos ξ k α k · · · cos ξ 1 α 1 - sin ξ 1 α 1 sin ξ 1 α 1 cos ξ 1 α 1 y x
下面介绍CORDIC实现坐标旋转具体步骤:
把原坐标 和旋转角θ作为计算步骤初值,令i=1,则步骤如下:
1.确定ξi的值
根据θ的正负,判断ξi的值。具体判定标准如下:如果θ大于0,令ξi=1;反之θ小于0,令ξi=-1。
2.计算
Figure A20051002187400096
旋转i个CORDIC基本角之后的坐标
y x = 1 - ξ i 2 - i ξ i 2 - i 1 y x - - - ( 9 )
式(9)使用加法和移位就可以实现。
3.判定是否旋转完毕
如果i=k,说明旋转完毕,转到第4步进行,否则计算θ=θ-ξiαii为固定常数arctg2-i)然后i=i+1,后回到第1步重新计算;
4.模值修正:
把步骤3计算出的坐标
Figure A20051002187400098
乘以固定常数const即可 ( const = Π i = 1 k cos α i ) . 此时输出值就是所求。
发明内容:
本发明的主要任务就是提供一种基于快速、芯片可实现的特征值算法的MUSIC算法来实现波达方向估计。
一种芯片可实现的多信号分类算法,其步骤如下:
步骤一:采集数据,构造数据矩阵
首先在阵列天线的N个阵元并行采集M个数据,形成数据矩阵A:
A = x 1 ( 1 ) x 1 ( 2 ) · · · x 1 ( j ) · · · x 1 ( M ) · · · · · · · · · · · · · · · · · · x i ( 1 ) x i ( 2 ) · · · x i ( j ) · · · x i ( M ) · · · · · · · · · · · · · · · · · · x N ( 1 ) x N ( 2 ) · · · x N ( j ) · · · x N ( M )
其中xi(j)中的i代表第i(i≤N)个阵元,j代表第j(j≤N)个采集数据;数据矩阵A可用存储器RAM存储。
步骤二:计算协方差矩阵:
在得到数据矩阵A后,则可以计算协方差矩阵R:
              R=AAT
协方差矩阵R是一个对称实矩阵,具体表达形式如下:
Figure A20051002187400102
其中(p<q≤N)
该步骤可用
Figure A20051002187400103
个乘法累加器实现。
步骤三:特征值分解(EVD):
已知协方差矩阵是R,设特征向量矩阵V的初值是N×N的单位矩阵。定义CORDIC基本的角度集:
1,α2…αORDER}={arctg2-1,arctg2-2…arctg2-ORDER}和CORDIC符号集:
1,ξ2…ξORDER,ξi=1或-1,i=1,2…ORDER}
其中,ORDER是CORDIC阶数,根据精度选择,是固定常数。
具体步骤如下:
1.从R选择消去对象rpq
按照r12,r13,…r1N;r23,r24…r2N;……;r(N-1)N矩阵元素所在的位置的先后顺序从R选择一个消去对象rpq;当选择r(N-1)N后,下次选择r12并重复上述顺序重新进行选择;本步骤可用基于状态机的数据选择器实现。
2.改进的坐标旋转:
(1)计算γ=2*rpq、β=rpp-rqq;本步骤可用移位器和减法器实现;
(2)根据γ和β确定ξi。如果γ*β≥0,ξi=1,否则ξi=-1;
(3)计算γ′和β′,
γ ′ β ′ = 1 - 2 - 2 i - ξ i 2 - ( i - 1 ) ξ i 2 - ( i - 1 ) 1 - 2 - 2 i γ β ;
本步骤可用移位器和加法器实现;
(4)旋转第i个CORDIC基本旋转角,
r pj ′ r qj ′ = 1 - ξ i 2 - i ξ i 2 - i 1 r pj r qj j ≠ p , q ;
v jp ′ v jq ′ = 1 - ξ i 2 - i ξ i 2 - i 1 v jp v jq j任意值
本步骤可用移位器和加法器实现;
(5)判定是否旋转完所有基本旋转角:
如果i=ORDER,则执行第(6)步;否则i+1,令
γ = γ ′ ; β = β ′ ; r pj r qj = r pj ′ r qj ′ ( j ≠ p , q ) ; v jp v jq = v jp ′ v jq ′ (j是任意值)后,回到第(2)
步重新执行;本步骤可用比较器实现;
(6)模值修正
r pj ′ r qj ′ = const * r pj ′ r qj ′ ( j ≠ p , q ) , v jp v jq = const * v jp ′ v jq ′ (j任意值)
(其中 const = Π i = 1 order cos α i , 在ORDER确定条件下,const近似等于常数);本步骤可用移位器和加法器实现;
(7)在(6)计算的同时,通过查表(由CORDIC符号集{ξ1,ξ2…ξORDER}映射到值tgθopt)和乘法器后求出rpp′、rqq′,具体公式如下:
rpp′=rpp-rpqtgθopt
rqq′=rqq+rpqtgθopt
而对rpq′直接置0:
rpq′=0;
3.判断R′
如果R′不是三角型矩阵,则令R=R′、V=V′后,返回第1步选择下一个消去对象,然后重复步骤1-2;如果R′已经转化为三角型,此时R′的主对角线元素就是特征值,V′就是特征向量矩阵;本步骤可用比较器实现;
4.排序
对R′中的特征值按从大到小顺序,对R′,V′进行对应调整,使得{λ1,λ2…λj…λN}(λ1≥λ2…λj…≥λN)和V={v1 v2…vj…vN},其中Vj是λj对应的特征向量;本步骤可用比较器和数据选择器实现;
步骤四.构造噪声子空间:
根据特征值分解后得到的特征值估计出D大小,然后从V中去构造噪声子空间EN H,其中 E N H = [ v D + 1 · · · v N ] ; 本步骤可用比较器和数据选择器实现;
步骤五.谱峰搜索:
阵列天线的方向矩阵由信号入射方向θk(k=1,2…D)、天线阵列几何位置(均匀圆阵、均匀线阵)等决定。设方向矩阵的列矢量是α(θ)(当θ一定时候,该值为常数),算法证明,当θ=θk E N H a ( θ k ) = 0 .
于是构造空间谱函数P(θ),在θ∈[0,180]的范围内进行谱峰搜索:
P ( θ ) = 1 | | E N H a ( θ ) | | 2
P(θ)中D个峰值对应的D个θ就是待求的信号入射方向θk。本步骤可用加法器、乘法器、移位器、比较器和数据选择器实现。
本发明的本质在于:本发明使用了由本发明提出的改进的坐标旋转算法来实现顺序Jacobi算法中的反正切运算和坐标旋转,而通常使用CORDIC算法需要分别实现顺序Jacobi算法中的反正切运算和坐标旋转(在本发明中不需要使用CORDIC反正切求出坐标旋转的旋转角);同时改进了顺序Jacobi算法的迭代结构,把原来需要进行2次坐标旋转的计算复杂度降低为只需要进行1次坐标旋转。以上两点改进,把顺序Jacobi算法单次迭代中的3次CORDIC运算(2次CORDIC坐标旋转和1次CORDIC反正切运算)大约降低到1.2次CORDIC运算(1次改进的坐标旋转、1次乘法和1次查表)。
本发明正是基于上述两点改进,提出了本发明所述的一种芯片可实现的多信号分类算法,在不增加芯片复杂程度和保证同等精度要求的情况下,可以提高算法并行性,大大提高运算速度,减小运算时间。在芯片中实现EVD的运算时间降低到原来运算时间的1/3左右,而完成整个芯片可实现的多信号分类算法的运算时间降低到1/2左右。比如:对于4个阵元的阵列天线来说,在芯片中实现整个EVD运算时间约为5.76微秒,实现整个MUSIC算法的时间在20微秒以下;而现有的多信号分类算法在芯片实现中整个EVD运算时间约为16.7微秒以上,整个算法运算时间在30微秒以上。
附图说明:
图1:8阶均匀线阵示意图。
图2:8阶均匀圆阵示意图。
图3:MUSIC算法步骤框图。
图4:基于CORDIC算法的顺序Jacobi特征值分解流程图。
图5:本发明所述的基于改进坐标旋转算法的顺序Jacobi特征值分解流程图。

Claims (4)

1、一种芯片可实现的多信号分类算法,其步骤如下:
步骤一:采集数据,构造数据矩阵
首先在阵列天线的N个阵元并行采集M个数据,形成数据矩阵A:
A = x 1 ( 1 ) x 1 ( 2 ) · · · x 1 ( j ) · · · x 1 ( M ) · · · · · · · · · · · · · · · · · · x i ( 1 1 ) x i ( 2 ) · · · x i ( j ) · · · x i ( M ) · · · · · · · · · · · · · · · · · · x N ( 2 ) x N ( 2 ) · · · x N ( j ) · · · x N ( M )
其中xi(j)中的i代表第i(i≤N)个阵元,j代表第j(j≤N)个采集数据;数据矩阵A可用存储器RAM存储;
步骤二:计算协方差矩阵:
在得到数据矩阵A后,则可以计算协方差矩阵R:
                   R=AAT
协方差矩阵R是一个对称实矩阵,具体表达形式如下:
Figure A2005100218740002C2
其中(p<q≤N)
该步骤可用
Figure A2005100218740002C3
个乘法累加器实现;
步骤三:特征值分解(EVD):
已知协方差矩阵是R,设特征向量矩阵V初值是N×N的单位矩阵,定义CORDIC基本的角度集:
1,α2…αORDER}={arctg2-1,arctg2-2…arctg2-ORDER}和CORDIC符号集:
1,ξ2…ξORDER,ξi=1或-1,i=1,2…ORDER}
其中,ORDER是CORDIC阶数,根据精度选择,是固定常数;
具体步骤如下:
1.从R选择消去对象rpq
按照r12,r13,…r1N;r23,r24…r2N;……;r(N-1)N矩阵元素所在的位置的先后顺序从R选择消去对象rpq;当选择r(N-1)N后,下次选择r12并重复上述顺序重新进行选择;本步骤可用基于状态机的数据选择器实现;
2.改进的坐标旋转:
(1)计算γ=2*rpq、β=rpq-rqq,令i=1;本步骤可用移位器和减法器实现;
(2)根据γ和β确定ξi。如果γ*β≥0,ξi=1,否则ξi=-1;
(3)计算γ′和β′, γ ′ β ′ = 1 - 2 - 2 i - ξ i 2 - ( i - 1 ) ξ i 2 - ( i - 1 ) 1 - 2 - 2 i γ β ; 本步骤可用移位器和加法器实现;
(4)旋转第i个CORDIC基本旋转角,
r pj ′ r qj ′ = 1 - ξ i 2 - i ξ i 2 - i 1 r pj r qj - - - j ≠ p , q ;
v jp ′ r jq ′ = 1 - ξ i 2 - i ξ i 2 - i 1 r jp r jq j任意值;
本步骤可用移位器和加法器实现;
(5)判定是否旋转完所有基本旋转角:
如果i=ORDER,则执行第(6)步;否则i+1,令
γ=γ′;β=β′; r pj r qj = r pj ′ r qj ′ ( j ≠ p , q ) ; v jp v jq = v jp ′ v jq ′ (j是任意值)后,回到第(2)
步重新执行;本步骤可用比较器实现;
(6)模值修正
r pj ′ r qj ′ = const * r pj ′ r qj ′ ( j ≠ p , q ) , v jp v jq = const * v jp ′ v jq ′ (j任意值)
(其中 const = Π i = 1 order cos α i , 在ORDER确定条件下,const近似等于常数);本步骤可用移位器和加法器实现;
(7)在(6)计算的同时,通过查表(由确定后的CORDIC符号集{ξ1,ξ2…ξORDER,}映射到tgθopt)和乘法器后求出rpp′、rqq′,具体公式如下:
rpp′=rpp-rpqtgθopt
rqq′=rqq+rpqtgθopt
而对rpq′直接置0:
rpq′=0;
3.判断R′
如果R′不是三角型矩阵,则令R=R′、V=V′后,返回第1步选择下一个消去对象,然后重复步骤1-2;如果R′已经转化为三角型,此时R′的主对角线元素就是特征值,V′就是特征向量矩阵;本步骤可用比较器实现;
4.排序
对R′中的特征值按从大到小顺序,对R′,V′进行对应调整,使得{λ1,λ2…λj…λN}(λ1≥λ2…λj…≥λN)和V={v1 v2…vj…vN},其中vj是λj对应的特征向量;本步骤可用比较器和数据选择器实现;
步骤四.构造噪声子空间:
根据特征值分解后得到的特征值估计出D大小,然后从V中构造噪声子空间EN H,其中 E N H = [ v D + 1 · · · v N ] ; 本步骤可用比较器和数据选择器实现;
步骤五.谱峰搜索:
阵列天线的方向矩阵由信号入射方向θk(k=1,2…D)、天线阵列几何位置(均匀圆阵、均匀线阵)等决定。设方向矩阵的列矢量是a(θ)(当θ一定时候,该值为常数),算法证明,当θ=θk E N H a ( θ k ) = 0 ;
于是构造空间谱函数P(θ),在θ∈[0,180]的范围内进行谱峰搜索:
P ( θ ) = 1 | | E N H a ( θ ) | | 2
P(θ)中D个峰值对应的D个θ就是待求的信号入射方向θk;本步骤可用加法器、乘法器、移位器、比较器和数据选择器实现。
CNB2005100218746A 2005-10-18 2005-10-18 一种利用阵列天线进行波达方向估计的方法 Expired - Fee Related CN100385249C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005100218746A CN100385249C (zh) 2005-10-18 2005-10-18 一种利用阵列天线进行波达方向估计的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005100218746A CN100385249C (zh) 2005-10-18 2005-10-18 一种利用阵列天线进行波达方向估计的方法

Publications (2)

Publication Number Publication Date
CN1746697A true CN1746697A (zh) 2006-03-15
CN100385249C CN100385249C (zh) 2008-04-30

Family

ID=36166317

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005100218746A Expired - Fee Related CN100385249C (zh) 2005-10-18 2005-10-18 一种利用阵列天线进行波达方向估计的方法

Country Status (1)

Country Link
CN (1) CN100385249C (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101051083B (zh) * 2006-04-07 2010-11-03 中国科学院声学研究所 一种二维波达方向估计声纳信号处理方法
CN101325807B (zh) * 2008-07-24 2011-01-12 中国人民解放军理工大学 信号波达方向估计方法
CN101977068A (zh) * 2010-09-28 2011-02-16 电子科技大学 不需目标仰角参数对目标来波方向进行快速测定的方法
CN101149429B (zh) * 2006-09-18 2011-03-30 中国人民解放军空军工程大学导弹学院 适用于均匀圆阵的阵列互耦校正与信源测向方法
CN103294649A (zh) * 2013-05-23 2013-09-11 东南大学 双边cordic运算单元及基于该运算单元的并行雅克比埃尔米特阵特征分解方法和实现电路
CN106169941A (zh) * 2016-08-31 2016-11-30 电子科技大学 一种基于噪声功率的圆阵的doa 估计方法
CN107092006A (zh) * 2017-04-07 2017-08-25 广东精点数据科技股份有限公司 用于doa估计的阵列信号处理方法和系统
CN109426484A (zh) * 2017-08-28 2019-03-05 华为技术有限公司 一种数据排序装置、方法及芯片
CN109633525A (zh) * 2018-11-29 2019-04-16 上海无线电设备研究所 一种快速的均匀圆阵空间谱测向实现方法
CN110110285A (zh) * 2019-04-10 2019-08-09 浙江大学 一种用于FPGA的并行Jacobi计算加速实现方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107402371A (zh) * 2016-05-20 2017-11-28 国家海洋技术中心 一种智能天线阵

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001305202A (ja) * 2000-04-24 2001-10-31 Toyota Central Res & Dev Lab Inc Musicスペクトラム計算方法、その装置及び媒体
CN1297822C (zh) * 2003-02-21 2007-01-31 重庆邮电学院 一种基于td-scdma无线定位来波方向的估计方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101051083B (zh) * 2006-04-07 2010-11-03 中国科学院声学研究所 一种二维波达方向估计声纳信号处理方法
CN101149429B (zh) * 2006-09-18 2011-03-30 中国人民解放军空军工程大学导弹学院 适用于均匀圆阵的阵列互耦校正与信源测向方法
CN101325807B (zh) * 2008-07-24 2011-01-12 中国人民解放军理工大学 信号波达方向估计方法
CN101977068B (zh) * 2010-09-28 2013-12-11 电子科技大学 不需目标仰角参数对目标来波方向进行快速测定的方法
CN101977068A (zh) * 2010-09-28 2011-02-16 电子科技大学 不需目标仰角参数对目标来波方向进行快速测定的方法
CN103294649B (zh) * 2013-05-23 2016-08-10 东南大学 并行双边coridc运算单元、基于该运算单元的并行雅可比运算的埃尔米特阵特征分解实现电路和实现方法
CN103294649A (zh) * 2013-05-23 2013-09-11 东南大学 双边cordic运算单元及基于该运算单元的并行雅克比埃尔米特阵特征分解方法和实现电路
CN106169941A (zh) * 2016-08-31 2016-11-30 电子科技大学 一种基于噪声功率的圆阵的doa 估计方法
CN107092006A (zh) * 2017-04-07 2017-08-25 广东精点数据科技股份有限公司 用于doa估计的阵列信号处理方法和系统
CN109426484A (zh) * 2017-08-28 2019-03-05 华为技术有限公司 一种数据排序装置、方法及芯片
CN109633525A (zh) * 2018-11-29 2019-04-16 上海无线电设备研究所 一种快速的均匀圆阵空间谱测向实现方法
CN110110285A (zh) * 2019-04-10 2019-08-09 浙江大学 一种用于FPGA的并行Jacobi计算加速实现方法
CN110110285B (zh) * 2019-04-10 2020-05-22 浙江大学 一种用于FPGA的并行Jacobi计算加速实现方法

Also Published As

Publication number Publication date
CN100385249C (zh) 2008-04-30

Similar Documents

Publication Publication Date Title
CN1746697A (zh) 一种芯片可实现的多信号分类算法
CN1262130C (zh) 一种估计固定波束空间到达方向的方法
CN1264379C (zh) 时分双工移动通信系统中上下行资源的分配方法和装置
CN1109410C (zh) 在无线通信系统中减少干扰降低噪声影响的信号处理装置及方法
CN1925362A (zh) 一种基于均匀线阵的智能天线的实现方法
CN1738238A (zh) 高速可配置rsa加密算法及协处理器
CN1253044C (zh) 天线阵列系统、控制其方向图的方法以及移动终端
CN1179597C (zh) 无线基站
CN1494773A (zh) 无线基站系统及指向性控制方法
CN1086873C (zh) 合成无线电信号的无线电通信装置
CN101075843A (zh) 一种用于联合检测的同频相邻小区选择的方法及装置
CN1859019A (zh) 一种实现自动频率控制的方法及装置
CN1411189A (zh) 双层加权并行干扰对消算法的简化算法
CN1207859C (zh) 空时解相关码分多址多径信号波达方向估计方法和装置
CN1886956A (zh) 多天线接收机的方法和装置
CN101055467A (zh) 一种基于idp的工业过程动态优化系统及方法
CN1992956A (zh) 一种基于智能天线系统的上下行信号处理方法
CN1860700A (zh) 多址干扰消除
CN1486098A (zh) 双模终端装置和小区搜索方法
CN1520199A (zh) 一种多时隙分组无线资源分配方法
CN1272948C (zh) 一种同步码分多址用户信号二维波达方向估计方法和装置
CN101064581A (zh) 一种解扰解扩的方法与装置
CN1272736C (zh) 基于模块变形和概率局部搜索的集成电路模块级布局方法
CN1160890C (zh) 无线通信系统中数字波束形成方法、模块及其阵列接收机
CN1266861C (zh) 盲自适应波束形成技术

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080430

Termination date: 20101018