CN105137454A - 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 - Google Patents
一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 Download PDFInfo
- Publication number
- CN105137454A CN105137454A CN201510433748.5A CN201510433748A CN105137454A CN 105137454 A CN105137454 A CN 105137454A CN 201510433748 A CN201510433748 A CN 201510433748A CN 105137454 A CN105137454 A CN 105137454A
- Authority
- CN
- China
- Prior art keywords
- covariance matrix
- passage
- matrix
- signal
- iteration
- 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
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 118
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000013598 vector Substances 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000012544 monitoring process Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 230000017105 transposition Effects 0.000 claims description 5
- 230000021615 conjugation Effects 0.000 claims description 3
- 235000002595 Solanum tuberosum Nutrition 0.000 claims description 2
- 244000061456 Solanum tuberosum Species 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000001629 suppression Effects 0.000 abstract 1
- 238000005303 weighing Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法及实现装置,属于抗干扰技术领域。本发明采用数字下变频模块将天线阵列输入的中频信号转换为基带信号;通道一致性监测与补偿模块对N个通道下变频后的基带信号进行补偿;协方差矩阵计算模块是计算补偿后的信号的协方差矩阵;协方差矩阵特征分解模块对协方差矩阵进行特征分解,求出特征值和特征向量;加权输出模块根据求出的特征向量进行加权得到最优权值,最终输出抗干扰之后的信号。本发明在FPGA实现时资源占用少,权值更新速度快,在输入信号之前加入了通道一致性检测与补偿模块,减少了由于通道不一致性带来的干扰抑制能力的损失,提高了抗干扰算法的鲁棒性。
Description
技术领域
本发明属于抗干扰技术领域,具体涉及一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法及实现装置。
背景技术
在***战争中,美军前所未有的大量使用了GPS(GlobalPositioningSystem)精确制导技术,而经过GPS技术制导的武器,命中率急剧提高,成为左右战争形势的关键因素之一。随着战争的结束,人们对卫星导航的认识逐步加深,经过几十年的发展,许多国家有了自己的卫星导航***(包括中国的北斗***),这些***统称全球导航卫星***(GlobalNavigationSatelliteSystem,GNSS),GNSS在军用、民用领域都发挥了巨大的优势,已成为信息体系的重要基础设施。
由于导航卫星固有的脆弱性,在信号传播过程中受到复杂的电磁环境影响,以及对接收机的抗干扰研究不够深入等原因,使得卫星导航***很容易受到各种有意无意的干扰而影响其精度甚至失去导航定位的能力。因此,研究抗干扰技术,对提高GNSS的抗干扰能力、准确无误的利用其提供的信息具有重要意义。
目前各类抗干扰技术理论主要针对GNSS接收机,其中基于输入信号协方差矩阵特征分解的抗干扰技术发展较早、理论成熟,该技术是通过对输入信号协方差矩阵进行特征分解,分类出相互正交的干扰和噪声子空间,最后利用干扰空间和噪声空间的正交性得出最优权值加权到天线阵列实现抗干扰;但是,由于协方差矩阵是复数矩阵,现有的矩阵特征分解理论在复数域并不成立,并且矩阵分解计算量巨大,这些因素是FPGA实现上述抗干扰算法的瓶颈。2009年,申正义等在《均匀圆阵DOA估计的抗干扰预处理》中提出了一种将协方差矩阵由复数域转换为实数域的方法,从而避免对复数矩阵进行分解,但该方法仅适用于基于谱峰搜索的波达方向(DOA)估计,在基于协方差矩阵特征分解的抗干扰算法中需要求出复数权值,该方法显然不适合。
发明内容
本发明针对现有技术中提到的基于特征分解的抗干扰算法FPGA实现困难的问题,提出了一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法及实现装置,所述实现方法计算量小,权值更新速度快,有效提高了抗干扰能力。
所述的实现装置包括在FPGA上设置的数字下变频(DDC)模块、通道一致性监测与补偿模块、协方差矩阵计算模块、协方差矩阵特征分解模块和加权输出模块。
所述的数字下变频(DDC)模块将天线阵列输入的中频信号转换为基带信号。包括:DDS(DirectDigitalSynthesizer)产生相互正交的正弦、余弦信号;将输入的中频信号分别和正、余弦信号进行混频,利用低通滤波器对混频的信号进行滤波,得到基带信号的同相分量和正交分量。
所述的通道一致性监测与补偿模块,包括:求出N(N为抗干扰天线阵元个数)个通道的瞬时相位和瞬时幅度并做平均;将除第一通道外的其它N-1个通道和第一通道的平均瞬时幅度做除法,平均瞬时相位做减法,求出除第一通道外的其它N-1个通道相对于第一通道的幅度不一致性和相位不一致性;对N个通道下变频后的基带信号进行补偿。
所述的协方差矩阵计算模块是计算补偿后的信号的协方差矩阵。
所述的协方差矩阵特征分解模块是对协方差矩阵进行特征分解,求出特征值和特征向量,该模块运用迭代的方法对协方差矩阵进行分解。包括:初始化、设置当前迭代状态、进迭代运算、条件判断。
所述的加权输出模块是根据求出的特征向量进行加权得到最优权值wopt,最终输出抗干扰之后的信号。加权公式为:y=wopt HX,其中y为输出,X为补偿后的信号组成的矩阵。
本发明具有的优点和积极效果在于:
(1)FPGA可实现性:本发明提出的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,克服了传统抗干扰理论在FPGA实现时特征分解模块计算量巨大,实现困难的问题,通过直接在复数域利用迭代的方法对协方差矩阵进行特征分解,实现抗干扰,FPGA实现简单,资源占用少。
(2)权值更新速度快:本发明提出的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,协方差矩阵特征分解模块在利用迭代的方法进行特征分解,收敛速度较快,对于一个复矩阵迭代开始到结束耗时不到1μs,最优权值由协方差矩阵的特征向量构成,因此最优权值的更新速度也比较快。
(3)鲁棒性:本发明提出的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,在输入信号之前加入了通道一致性检测与补偿模块,减少了由于通道不一致性带来的干扰抑制能力的损失,提高了抗干扰算法的鲁棒性。
附图说明
图1是本发明的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法的工作流程图。
具体实施方式
下面将结合附图对本发明作进一步的详细说明。
本发明首先提出的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现装置,如图1所示,在FPGA上实现了五个模块:通道一致性监测与补偿模块、数字下变频模块、协方差矩阵计算模块、协方差矩阵特征分解模块和加权输出模块。
所述的数字下变频模块接收天线阵列的中频信号,并将中频信号转换为基带信号发送给通道一致性监测与补偿模块,在所述的通道一致性监测与补偿模块中对基带信号进行补偿;所述的协方差矩阵计算模块根据补偿后的信号计算协方差矩阵;协方差矩阵特征分解模块利用迭代的方法对协方差矩阵进行特征分解,求出特征值和特征向量,并在加权输出模块中根据所述的特征值和特征向量进行加权输出。
基于上述的实现装置,本发明还提供一种所述实现装置的实现方法,具体步骤如下:
步骤一、数字下变频(DDC)模块接收天线阵列的中频信号;
数字下变频模块将天线阵列输入的中频信号转换为基带信号,并输出基带信号的同相分量和正交分量,具体为:
(1.1)产生DDS(DirectDigitalSynthesizer,直接数字式频率合成器),DDS时钟设置为采样率,DDS以采样率为数据率输出相互正交的正弦信号sin(w0t)和余弦信号cos(w0t)。
(1.2)将DDS产生的正余弦信号分别和输入的中频信号进行混频,实现方式是有符号定点乘法器,输入的中频信号和正弦信号相乘得到的是包含高频分量和基带信号ddcq,输入的中频信号和余弦信号相乘得到的是包含高频分量和基带信号的ddci。
(1.3)对ddci和ddcq分别进行滤波,滤波器选择两个低通滤波器,截止频率分别为ddci和ddcq信号带宽的一半,滤除高频分量,输出分别是:基带信号的同相分量xi,基带信号的正交分量xq,对N个通道做相同的处理,则第n个通道下变频模块输出的基带信号的同相分量为xni,基带信号的正交分量xnq。
步骤二、通道一致性监测与补偿模块实现通道不一致性的监测与补偿功能,具体为:
(2.1)对N(N为天线阵列中抗干扰天线阵元个数)个通道分别采集下变频后M(M>1024)个采样点的数据,求出每个采样点的瞬时幅度和瞬时相位,计算公式为:
其中An(k)是第n个通道的第k个采样点的瞬时幅度,θn(k)是第n个通道第k个采样点的瞬时相位,xnq(k)、xni(k)分别是步骤一中数字下变频模块输出的对应于第n个通道第k个采样点的正交分量和同相分量。
(2.2)对M个采样点的瞬时幅度和瞬时相位做平均,求出N个通道的平均瞬时幅度和平均瞬时相位计算公式为:
其中,M为采样点的个数;求平均是为了保证每个采样点都能得到较为理想的补偿从而使整体补偿误差减小。
(2.3)将除第一通道外的其他N-1个通道和第一通道的平均瞬时幅度做除法,平均瞬时相位做减法,求出除第一通道外的其它N-1个通道相对于第一通道的幅度不一致性(平均瞬时幅度的比值)和相位不一致性(平均瞬时相位的差值);设除第一通道外的其它N-1个通道相对于第一通道的幅度不一致性和相位不一致性分别是:
则补偿就是将各通道下变频输出与由幅度不一致性、相位不一致性构成的复数做复数乘法,复数乘法公式为:
Xnin=xni+jxnq为步骤一中第n通道的数字下变频输出的基带信号,xni,xnq分别为数字下变频模块输出的对应于第n个通道的正交分量和同相分量。
Andiffer为第n通道相对于第一通道的幅度不一致性,θndiffer为第n通道相对于第一通道的相位不一致性。
步骤三、协方差矩阵计算模块是计算补偿后的信号Xnout的协方差矩阵,具体为:
(3.1)利用复数乘法器计算协方差矩阵,计算公式为:Rx=XXH,H表示共轭转置;其中X是由步骤二中的各通道输出补偿后的信号组成的矩阵,定义为:
XH=[Xlout *…XNout *],*表示共轭。
(3.2)对上一步中的Rx做累加:M为采样点的个数,累加用加法器实现,为保证累加准确,在FPGA中该加法器的工作时钟必须不小于输入信号数据率的两倍。
(3.3)对Rxx做平均,则协方差矩阵M为采样点个数,R是一个N行、N列的矩阵,N是抗干扰天线阵元的个数。
步骤四、协方差矩阵特征分解模块利用迭代的方法对协方差矩阵R进行特征分解,求出特征值和特征向量,该模块是FPGA实现抗干扰算法的关键,具体为:
(4.1)初始化。初始化迭代过程中的一些参数,这些参数和数值分别是:矩阵E,E初始化为N*N元的单位矩阵;当前迭代轮次l=0,迭代轮次上限L,5≤L≤50;非对角线元素2范数门限ε,0≤ε≤10-5。
(4.2)设置当前迭代状态。就是设置当前迭代要消去的元素所在的行、列,按照矩阵元素从左到右、从上到下的顺序设置非对角线元素所在的行、列。以4*4矩阵B为例,要消去的元素分别为:B(1,2),B(1,3),B(1,4),B(2,3),B(2,4),B(3,4),其中B(1,2)表示矩阵B的第一行第二列的元素,以此类推。
(4.3)迭代运算。迭代运算是构造一个正交矩阵,对协方差矩阵和单位矩阵进行一系列正交运算,使得协方差矩阵经过这些变换后化为近似的对角矩阵,由于正交运算后矩阵特征值不变,因此对角矩阵对角线上的元素就是协方差矩阵的特征值。设第i次迭代后协方差矩阵变为R(i),单位矩阵变为E(i),第i+1次迭代的正交矩阵为P,则,
R(i+1)=PR(i)PH
E(i+1)=E(i)PH
其中,
S*为S的共轭,PH为P的共轭转置,R(i)(m,n)表示第i次迭代中矩阵R的第m行第n列元素。
迭代运算的具体步骤为:
(4.3.1):求出C,S。假设步骤4.2设置的迭代状态为第m行第n列,第i次迭代后,第m行m列和第n行m列元素的值分别xmm=R(i)(m,m),xnm=R(i)(n,m)。则:
|·|表示取模。
(4.3.2):构造矩阵P和PH。
矩阵左侧m,n表示第m行,第n行,下侧m,n表示第m列,第n列。
(4.3.3):迭代运算。第i+1次迭代的运算为:
R(i+1)=PR(i)PH
E(i+1)=E(i)PH
其中,R(i)是第i次迭代后的协方差矩阵,E(i)表示第i次迭代后的单位矩阵。
(4.3.4):条件判断。条件判断控制整个迭代运算是继续还是终止,条件判断的判断逻辑为:
(a)计算当前迭代轮次l。在步骤(4.2)中,若状态集里的所有迭代状态都经历过一次迭代,当最后一个迭代状态结束后,称一轮迭代完成,用迭代轮次l表示经历了多少轮的迭代。若在步骤(4.2)中设置的状态是状态集里的最后一个状态,则l=l+1,否则返回步骤(4.2)按照状态集里的状态设置下一个状态。
(b)计算每一轮迭代完毕后,协方差矩阵非对角元素的2范数。计算公式为:
其中,||r||2表示矩阵非对角元素的二范数,R(m,n)表示矩阵R的第m行第n列元素,|.|表示取绝对值。
(c)判断。判断逻辑为:若l>L或者||r||2<ε(其中l是步骤(a)中计算出的当前迭代轮次,L是步骤(4.1)初始化时设置的迭代轮次上限;||r||2是步骤(b)中计算出的矩阵非对角线元素的2范数,ε是步骤(4.1)初始化时设置的非对角线元素2范数门限),则终止迭代,输出此结果,结果为:此时R的对角线的值就是协方差矩阵的特征值λ1,…λN,E的列向量就是协方差矩阵特征值所对应的特征向量η1,…ηN,否则返回步骤(4.2)。
步骤五、加权输出。根据步骤四中求出的特征值和特征向量,对特征值λ1,…λN按照从小到大的顺序进行排序,找出和最小特征值λ1比值小于4dB的其它特征值,这些特征值为λ1,…λD,D<N,则这些特征值对应的特征向量η1,…ηD之和就是最优权值wopt,即:
wopt=η1+…+ηD
得到最优权值后对步骤二中的输出Xnout进行加权,加权公式为:
其中y为输出,X为由信号Xnout组成的矩阵,wopt H表示最优权值wopt的共轭转置。
Claims (5)
1.一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,其特征在于:具体步骤如下,
步骤一、数字下变频模块接收天线阵列的中频信号;
数字下变频模块将天线阵列输入的中频信号转换为基带信号,并输出基带信号的同相分量和正交分量;
步骤二、通道一致性监测与补偿模块实现通道不一致性的监测与补偿功能,具体为:
(2.1)对天线阵列中N个通道分别采集数字下变频后M个采样点的数据,求出每个采样点的瞬时幅度和瞬时相位,计算公式为:
其中An(k)是第n个通道的第k个采样点的瞬时幅度,θn(k)是第n个通道第k个采样点的瞬时相位,xnq(k)、xni(k)分别是步骤一中数字下变频模块输出的对应于第n个通道第k个采样点的正交分量和同相分量,N为天线阵列中抗干扰天线阵元个数;
(2.2)对M个采样点的瞬时幅度和瞬时相位做平均,求出N个通道的平均瞬时幅度和平均瞬时相位计算公式为:
其中,M为采样点的个数;
(2.3)将除第一通道外的其他N-1个通道和第一通道的平均瞬时幅度做除法,平均瞬时相位做减法,求出除第一通道外的其它N-1个通道相对于第一通道的幅度不一致性和相位不一致性;设除第一通道外的其它N-1个通道相对于第一通道的幅度不一致性和相位不一致性分别是:
则补偿就是将各通道数字下变频模块输出与由幅度不一致性、相位不一致性构成的复数做复数乘法,复数乘法公式为:
Xnin=xni+jxnq为步骤一中第n通道的数字下变频输出的基带信号,xni,xnq分别为数字下变频模块输出的对应于第n个通道的正交分量和同相分量;
步骤三、协方差矩阵计算模块是计算补偿后的信号Xnout的协方差矩阵,具体为:
(3.1)利用复数乘法器计算协方差矩阵,计算公式为:Rx=XXH,H表示共轭转置;其中X是由步骤二中的各通道输出补偿后的信号组成的矩阵,定义为:
(3.2)对上一步中的Rx做累加:累加用加法器实现;
(3.3)对Rxx做平均,则协方差矩阵M为采样点个数,R是一个N行、N列的矩阵,N是抗干扰天线阵元的个数;
步骤四、协方差矩阵特征分解模块利用迭代的方法对协方差矩阵R进行特征分解,求出特征值和特征向量,具体为:
(4.1)初始化;
初始化迭代过程中的参数,这些参数和数值分别是:矩阵E,E初始化为N*N元的单位矩阵;当前迭代轮次l=0,迭代轮次上限L,5≤L≤50;非对角线元素2范数门限ε,0≤ε≤10-5;
(4.2)设置当前迭代状态;
设置当前迭代要消去的元素所在的行、列,按照矩阵元素从左到右、从上到下的顺序设置非对角线元素所在的行、列;
(4.3)迭代运算;
迭代运算是构造一个正交矩阵,对协方差矩阵和单位矩阵进行一系列正交运算,使得协方差矩阵经过这些变换后化为近似的对角矩阵,由于正交运算后矩阵特征值不变,因此对角矩阵对角线上的元素就是协方差矩阵的特征值;迭代后单位矩阵的列向量就是协方差矩阵的特征值所对应的特征向量;
步骤五、加权输出;
根据步骤四中求出的特征值和特征向量,对特征值λ1,…λN按照从小到大的顺序进行排序,找出和最小特征值λ1比值小于4dB的其它特征值,这些特征值为λ1,…λD,D<N,则这些特征值对应的特征向量η1,…ηD之和就是最优权值wopt,即:
wopt=η1+…+ηD
得到最优权值后对步骤二中的输出Xnout进行加权,加权公式为:
其中y为输出,X为由信号Xnout组成的矩阵,wopt H表示最优权值wopt的共轭转置。
2.根据权利要求1所述的所述的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,其特征在于:步骤一具体实现方法如下,
(1.1)产生DDS,DDS时钟设置为采样率,DDS以采样率为数据率输出相互正交的正弦信号sin(w0t)和余弦信号cos(w0t);
(1.2)将DDS产生的正信号、余弦信号分别和输入的中频信号进行混频,实现方式是符号定点乘法器,输入的中频信号和正弦信号相乘得到的是包含高频分量和基带信号ddcq,输入的中频信号和余弦信号相乘得到的是包含高频分量和基带信号的ddci;
(1.3)对ddci和ddcq分别进行滤波,滤波器选择两个低通滤波器,截止频率分别为ddci和ddcq信号带宽的一半,滤除高频分量,输出分别是:基带信号的同相分量xi,基带信号的正交分量xq,对N个通道做相同的处理,则第n个通道经过数字下变频模块输出的基带信号的同相分量为xni,基带信号的正交分量xnq,n=1,2,…,N。
3.根据权利要求1所述的所述的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,其特征在于:所述加法器的工作时钟不小于输入信号数据率的两倍。
4.根据权利要求1所述的所述的一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现方法,其特征在于:所述的迭代运算的具体步骤为,
(4.3.1):求出C,S;
假设步骤4.2设置的迭代状态为第m行第n列,第i次迭代后,第m行m列和第n行m列元素的值分别xmm=R(i)(m,m),xnm=R(i)(n,m),则:
|·|表示取模;
(4.3.2):构造正交矩阵P和PH:
矩阵左侧m,n表示第m行,第n行,下侧m,n表示第m列,第n列;
(4.3.3):迭代运算:第i+1次迭代的运算为:
R(i+1)=PR(i)PH
E(i+1)=E(i)PH
其中,R(i)是第i次迭代后的协方差矩阵,E(i)表示第i次迭代后的单位矩阵;
(4.3.4):条件判断;
条件判断的判断逻辑为:
(a)计算当前迭代轮次l;在步骤(4.2)中,若状态集里的所有迭代状态都经历过一次迭代,当最后一个迭代状态结束后,称一轮迭代完成,用迭代轮次l表示经历了多少轮的迭代;若在步骤(4.2)中设置的状态是状态集里的最后一个状态,则l=l+1,否则返回步骤(4.2)按照状态集里的状态设置下一个状态;
(b)计算每一轮迭代完毕后,协方差矩阵非对角元素的2范数,计算公式为:
其中,||r||2表示矩阵非对角元素的二范数,R(m,n)表示矩阵R的第m行第n列元素,|.|表示取绝对值;
(c)判断;
判断逻辑为:若l>L或者||r||2<ε,则终止迭代,输出此结果,结果为:此时R的对角线的值就是协方差矩阵的特征值λ1,…λN,E的列向量就是协方差矩阵特征值所对应的特征向量η1,…ηN,否则返回步骤(4.2)。
5.一种基于协方差矩阵特征分解的抗干扰算法的FPGA实现装置,其特征在于:所述的实现装置包括在FPGA上实现的五个模块,分别为:通道一致性监测与补偿模块、数字下变频模块、协方差矩阵计算模块、协方差矩阵特征分解模块和加权输出模块;
所述的数字下变频模块接收天线阵列的中频信号,并将中频信号转换为基带信号发送给通道一致性监测与补偿模块,在所述的通道一致性监测与补偿模块中对基带信号进行补偿;
所述的协方差矩阵计算模块根据补偿后的信号计算协方差矩阵;协方差矩阵特征分解模块利用迭代的方法对协方差矩阵进行特征分解,求出特征值和特征向量,并在加权输出模块中根据所述的特征值和特征向量进行加权输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510433748.5A CN105137454B (zh) | 2015-07-22 | 2015-07-22 | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510433748.5A CN105137454B (zh) | 2015-07-22 | 2015-07-22 | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105137454A true CN105137454A (zh) | 2015-12-09 |
CN105137454B CN105137454B (zh) | 2017-08-08 |
Family
ID=54722860
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510433748.5A Expired - Fee Related CN105137454B (zh) | 2015-07-22 | 2015-07-22 | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137454B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106324625A (zh) * | 2016-08-30 | 2017-01-11 | 北京航空航天大学 | 一种基于二范数多目标优化的卫星导航***自适应抗干扰方法 |
CN106680837A (zh) * | 2016-12-14 | 2017-05-17 | 中国电子科技集团公司第二十研究所 | 一种卫星导航干扰抑制算法 |
CN108829903A (zh) * | 2017-11-09 | 2018-11-16 | 北京广利核***工程有限公司 | 判定fpga冗余设计的代码与综合后电路一致性的方法和*** |
CN111591248A (zh) * | 2020-04-16 | 2020-08-28 | 杭州勒格智能设备有限公司 | 车门控制通信***及方法 |
CN113075701A (zh) * | 2021-03-25 | 2021-07-06 | 陕西凌云电器集团有限公司 | 一种具有差分定位和抗干扰功能的卫星信号接收方法及模块 |
CN116028770A (zh) * | 2023-01-18 | 2023-04-28 | 珠海微度芯创科技有限责任公司 | 适用于实、复协方差矩阵的特征值分解硬件实现方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7508875B1 (en) * | 2005-09-15 | 2009-03-24 | Bae Systems Information And Electronic Systems Integration Inc. | Correlation matrix estimation via adaptive array perturbational processing to enable jammer location |
CN102170302A (zh) * | 2011-03-23 | 2011-08-31 | 西安电子科技大学 | 基于fpga的智能天线抗干扰片上***及方法 |
CN103630911A (zh) * | 2013-12-11 | 2014-03-12 | 北京北斗星通导航技术股份有限公司 | 一种导航信号的处理方法及装置 |
CN104345321A (zh) * | 2014-11-03 | 2015-02-11 | 西安烽火电子科技有限责任公司 | 一种卫星导航接收机用空时自适应抗干扰方法 |
CN104360354A (zh) * | 2014-12-05 | 2015-02-18 | 北京北斗星通导航技术股份有限公司 | 空频抗干扰处理方法和装置 |
-
2015
- 2015-07-22 CN CN201510433748.5A patent/CN105137454B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7508875B1 (en) * | 2005-09-15 | 2009-03-24 | Bae Systems Information And Electronic Systems Integration Inc. | Correlation matrix estimation via adaptive array perturbational processing to enable jammer location |
CN102170302A (zh) * | 2011-03-23 | 2011-08-31 | 西安电子科技大学 | 基于fpga的智能天线抗干扰片上***及方法 |
CN103630911A (zh) * | 2013-12-11 | 2014-03-12 | 北京北斗星通导航技术股份有限公司 | 一种导航信号的处理方法及装置 |
CN104345321A (zh) * | 2014-11-03 | 2015-02-11 | 西安烽火电子科技有限责任公司 | 一种卫星导航接收机用空时自适应抗干扰方法 |
CN104360354A (zh) * | 2014-12-05 | 2015-02-18 | 北京北斗星通导航技术股份有限公司 | 空频抗干扰处理方法和装置 |
Non-Patent Citations (2)
Title |
---|
何訸: "GNSS接收机抗干扰技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
倪丁华: "分块并行DBF算法及其实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106324625A (zh) * | 2016-08-30 | 2017-01-11 | 北京航空航天大学 | 一种基于二范数多目标优化的卫星导航***自适应抗干扰方法 |
CN106324625B (zh) * | 2016-08-30 | 2019-02-19 | 北京航空航天大学 | 一种基于二范数多目标优化的卫星导航***自适应抗干扰方法 |
CN106680837A (zh) * | 2016-12-14 | 2017-05-17 | 中国电子科技集团公司第二十研究所 | 一种卫星导航干扰抑制算法 |
CN106680837B (zh) * | 2016-12-14 | 2019-04-19 | 中国电子科技集团公司第二十研究所 | 一种卫星导航干扰抑制算法 |
CN108829903A (zh) * | 2017-11-09 | 2018-11-16 | 北京广利核***工程有限公司 | 判定fpga冗余设计的代码与综合后电路一致性的方法和*** |
CN108829903B (zh) * | 2017-11-09 | 2021-11-05 | 北京广利核***工程有限公司 | 判定fpga冗余设计的代码与综合后电路一致性的方法和*** |
CN111591248A (zh) * | 2020-04-16 | 2020-08-28 | 杭州勒格智能设备有限公司 | 车门控制通信***及方法 |
CN113075701A (zh) * | 2021-03-25 | 2021-07-06 | 陕西凌云电器集团有限公司 | 一种具有差分定位和抗干扰功能的卫星信号接收方法及模块 |
CN113075701B (zh) * | 2021-03-25 | 2023-04-18 | 陕西凌云电器集团有限公司 | 具有差分定位和抗干扰功能的卫星信号接收方法及模块 |
CN116028770A (zh) * | 2023-01-18 | 2023-04-28 | 珠海微度芯创科技有限责任公司 | 适用于实、复协方差矩阵的特征值分解硬件实现方法 |
CN116028770B (zh) * | 2023-01-18 | 2024-01-12 | 珠海微度芯创科技有限责任公司 | 适用于实、复协方差矩阵的特征值分解硬件实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105137454B (zh) | 2017-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105137454A (zh) | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 | |
CN102798840B (zh) | 外辐射源雷达宽带信道化接收***及fpga实现方法 | |
Djurovic et al. | A hybrid CPF-HAF estimation of polynomial-phase signals: Detailed statistical analysis | |
CN107290730B (zh) | 互耦条件下双基地mimo雷达角度估算方法 | |
CN102590829B (zh) | 用于卫星导航***的完全自适应陷波器及其陷波方法 | |
CN105204006A (zh) | 基于子空间的干扰噪声协方差矩阵重构的波束形成方法 | |
CN105589080B (zh) | 一种基于记忆矩阵准求逆的卫星导航抗干扰的装置及其抗干扰方法 | |
CN104898113A (zh) | 一种多输入多输出雷达波形设计方法 | |
CN105259550A (zh) | 基于压缩感知的多输入多输出雷达二维角度估计方法 | |
CN107219496B (zh) | 一种改进的相关干涉仪鉴相方法 | |
CN108896954A (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN102055512B (zh) | 天线组阵中基于循环自相关的信号相位差估计装置与方法 | |
CN105842713A (zh) | 一种基于sft的ins辅助bds信号快速捕获方法 | |
CN103763230A (zh) | 一种改进的自适应盲源分离方法 | |
CN105572473B (zh) | 高分辨率线性时频分析方法 | |
CN103338024B (zh) | 天线组阵中时延的互补卡尔曼滤波装置与方法 | |
CN109001687A (zh) | 基于广义旁瓣相消结构的机载雷达空时自适应滤波方法 | |
CN103323667A (zh) | 贝塞尔函数与虚拟阵列相结合的sfm信号的参数估计方法 | |
CN105721078B (zh) | 一种通用群时延纹波自动化测试方法 | |
CN109581277A (zh) | 一种基于压缩感知理论的四维天线阵doa估计方法 | |
CN113376570A (zh) | 一种基于信道化处理的相关干涉仪实时测向方法 | |
CN105510706A (zh) | 一种高精度欠采样测频方法 | |
CN105282067A (zh) | 一种复数域盲源分离方法 | |
CN110413939B (zh) | 一种基于原子范数的到达角估计方法 | |
CN104914451B (zh) | 一种块Toeplitz矩阵低复杂度求逆的空时抗干扰方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170808 |