CN112073076B - 一种基于短时傅里叶变换的空频自适应处理抗干扰方法 - Google Patents
一种基于短时傅里叶变换的空频自适应处理抗干扰方法 Download PDFInfo
- Publication number
- CN112073076B CN112073076B CN202010762611.5A CN202010762611A CN112073076B CN 112073076 B CN112073076 B CN 112073076B CN 202010762611 A CN202010762611 A CN 202010762611A CN 112073076 B CN112073076 B CN 112073076B
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- interference
- frequency
- noise
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 title claims abstract description 13
- 230000003044 adaptive effect Effects 0.000 title claims description 17
- 239000013598 vector Substances 0.000 claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 239000004576 sand Substances 0.000 claims description 22
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000002452 interceptive effect Effects 0.000 claims description 8
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000002904 solvent Substances 0.000 claims description 2
- 230000036039 immunity Effects 0.000 claims 1
- 230000001629 suppression Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/12—Neutralising, balancing, or compensation arrangements
- H04B1/123—Neutralising, balancing, or compensation arrangements using adaptive balancing or compensation means
- H04B1/126—Neutralising, balancing, or compensation arrangements using adaptive balancing or compensation means having multiple inputs, e.g. auxiliary antenna for receiving interfering signal
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Signal Processing (AREA)
- Data Mining & Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Discrete Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Noise Elimination (AREA)
Abstract
本发明公开了一种基于短时傅里叶变换的空频自适应处理抗干扰方法,包括获取阵列接收信号;将传感器阵列接收到的时域信号转换到时频域;将时频域划分为期望信号区域和干扰噪声区域;去除区域中的噪声点后在干扰区域中分离各个干扰信号;计算期望信号和各干扰信号的导向矢量,重构各频点的干扰加噪声协方差矩阵;通过MVDR算法计算各频点的自适应权矢量;本发明的方法运算简单,同时具有良好的鲁棒性,能够获得很好的干扰抑制效果。
Description
技术领域
本发明涉及自适应抗干扰技术领域,主要涉及一种基于短时傅里叶变换的空频自适应处理抗干扰方法。
背景技术
抗干扰算法有多种形式,早些时候广泛使用的是利用阵列天线在空域上进行的空域自适应处理(SAP),SAP算法计算复杂度较低但所能抑制的干扰数量受到阵元数量限制。为了能使阵列抑制更多的干扰,在SAP算法基础上又发展了空时自适应处理(STAP),STAP在SAP的基础上,在每个通道后面添加若干个时域抽头,使得阵列拥有更多的自由度来抑制干扰信号,STAP在不增加阵元数量同时提高了抗干扰能力,但STAP的计算复杂度较高。为了降低STAP的计算复杂度,又提出了一种空时自适应处理的次优算法即SFAP,其通过对接收信号做快速傅里叶变换(FFT)将宽带信号划分为若干频带,再在各频带上进行SAP。对于阵元数为M、延迟节点数为N的STAP算法,求解自适应权值要进行MN×MN的矩阵求逆运算,而SFAP是在各频带分别求取权值只需要进行M×M的矩阵求逆运算,从而降低了计算复杂度。
傅里叶变换是时域到频域互相转化的工具,但它没有时间定位或时间局域化的能力,为了解决这一问题,人们又在其基上发展了短时傅里叶变换(STFT)。STFT是一种典型的线性时频分析方法,它将时间信号变换到时频域,揭示了各频率分量随时间的变化趋势,把时域和频域有机地结合起来对信号进行分析。通过STFT同样可以对信号子带划分,在各子带上进行空域滤波完成SFAP。
发明内容
发明目的:本发明提供了一种基于短时傅里叶变换的空频自适应处理抗干扰方法,通过STFT将接收到的宽带信号变换至时频域,在时频域上分离各输入信号后估计各信号的延时τ,利用τ计算信号的导向矢量用以IPNC矩阵重构及波束形成。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种基于短时傅里叶变换的空频自适应处理抗干扰方法,包括如下步骤:
步骤S1、获取阵列接收信号S;所述阵列接收信号S为大小Ne×N的样本矩阵;Ne是阵列天线的阵元个数,N为采样点数;
步骤S2、采用短时傅里叶变换将传感器阵列接收到的时域信号转换到时频域,即对阵列接收信号S的各行数据分别进行短时傅里叶变换;
步骤S3、将时频域划分为期望信号区域Ωs和干扰噪声区域Ωi+n;去除区域中的噪声点;在去除噪声后的干扰区域Ωi中分离各个干扰信号;
步骤S4、分别计算期望信号和各干扰信号的导向矢量a(f);重构各频点的干扰加噪声协方差矩阵;
步骤S5、通过MVDR算法计算各频点的自适应权矢量ω(f)。
进一步地,所述步骤S2中采用短时傅里叶变换后的转换信号表示如下:
其中,ai(f)为第i个信号的导向矢量,K表示有K+1个接收信号,Si(t,f)为第i个信号的时频域表示;时频域中单信源点处的信号可表示为:
X(t,f)=ai(t,f)Si(t,f)
信号导向矢量ai(f)如下式所示:
其中τi为第i个信号的延时矢量,且τi,1=0;所述时频域中单信源点处的信号表示变形如下式:
其中X(t,f,1)为第一个通道数据的快速傅里叶变换结果。
进一步地,所述步骤S3中将时频域划分为期望信号区域Ωs和干扰噪声区域Ωi+n、去除噪声点并分离各个干扰信号的具体方法如下:
步骤S3.1、求取各时频点的信号导向矢量a(t,f);
步骤S3.2、根据期望信号的标称导向矢量a0(f)及给定的误差范围,与a(t,f)做比较后将时频域分成期望信号区域Ωs和干扰噪声区域Ωi+n;具体如下:
其中Δ为用于衡量各时频点导向矢量ai(t,f)和标称导向矢量a0(f)的相似度的参数,a(t,f)为信号导向矢量;
定义误差导向矢量ae(f)为:
其中τe为误差延时矢量;根据标称导向矢量a0(f)和误差导向矢量ae(f)设置阈值如下:
对时频域进行如下划分:
步骤S3.3、设定噪声门限,根据门限大小去除Ωi+n中的噪声时频点得到期望信号区域Ωs和干扰噪声区域Ωi;计算Ωs和Ωi中各频点处的信号延时,根据延时对Ωi中的干扰信号进行分离;
设置噪声门限Δ2,在每个时间切片上根据噪声门限Δ2剔除噪声点;对某一时间切片(ts,f)∈Ωi:
满足上式的时频点(ts,fs)即为噪声,从Ωi中剔除;在Ωs和Ωi中求出各时频点处信号的延时如下:
将Ωi中各时频点处的τi(t,f)根据大小分成J类,由于信号重叠部分较少,所以包含时频点点数较多的K类区域为各干扰信号所在区域剩余部分为信号重叠区域Ωo。分别将Ωs和中的τi(t,f)求和取平均,即可得到期望信号延时估计值τs和各干扰信号延时估计值τk。
进一步地,所述步骤S4中分别计算期望信号和各干扰信号的导向矢量a(f)并重构各频点的干扰加噪声协方差矩阵具体方法如下:
根据步骤S3获取的期望信号延时估计值τs和干扰信号的延时估计值τk计算期望信号和干扰信号的导向矢量估计值如下:
利用时频域上的频率切片计算各频点的协方差矩阵如下:
其中N为频率切片长度;利用Rs(f)计算Capon空间谱:
则干扰信号的协方差矩阵如下所示:
因此噪声协方差矩阵Rn(f):
Rn(f)=σn 2I
其中I为单位矩阵;σn 2为噪声功率估计值可由Rs(f)的最小特征值得到,对Rs(f)做特征值分解有:
其中λi和ui分别表示第i个特征值及该特征值对应的特征向量,则有σn 2=λmin;
进一步地,所述步骤S5中自适应权矢量ω(f)如下所示:
其中Ri+n -1(f)为步骤S4中求得的干扰加噪声协方差矩阵的逆矩阵;as(f)为期望信号的导向矢量;[·]H为共轭转置算子。
有益效果:本发明具备以下优点:
相比于传统SFAP算法,使用STFT代替FFT,提高了数据利用率,节省了数据存储资源;能够利用时频谱简洁地重构干扰加噪声协方差矩阵矩阵,使得MVDR算法能够获得更好的抗干扰效果;最后计算过程中不需要过多的先验知识,更加便捷。
附图说明
图1是本发明的一种基于STFT的SFAP抗干扰方法流程图;
图2是天线阵元的分布图;
图3是接收信号的时频谱;
图4a是分离出的期望信号时频谱;
图4b是分离出的第一个二次调频干扰信号时频谱;
图4c是分离出的第二个二次调频干扰信号时频谱;
图5a是由权矢量画出的二维方向图;
图5b是二维方向图在中心频率处的剖面图;
图6是入射期望信号和自适应滤波输出结果的对比图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示的一种基于短时傅里叶变换的空频自适应处理抗干扰方法,包括如下步骤:
步骤S1、获取阵列接收信号S;所述阵列接收信号S为大小Ne×N的样本矩阵;Ne是阵列天线的阵元个数,N为采样点数。
步骤S2、采用短时傅里叶变换将传感器阵列接收到的时域信号转换到时频域,即对阵列接收信号S的各行数据分别进行短时傅里叶变换。其中变换后的转换信号表示如下:
其中,ai(f)为第i个信号的导向矢量,K表示有K+1个接收信号,Si(t,f)为第i个信号的时频域表示;时频域中单信源点处的信号可表示为:
X(t,f)=ai(t,f)Si(t,f)
信号导向矢量ai(f)如下式所示:
其中τi为第i个信号的延时矢量,且τi,1=0;所述时频域中单信源点处的信号表示变形如下式:
其中X(t,f,1)为第一个通道数据的快速傅里叶变换结果。
步骤S3、将时频域划分为期望信号区域Ωs和干扰噪声区域Ωi+n;去除区域中的噪声点;在去除噪声后的干扰区域Ωi中分离各个干扰信号。具体地,
步骤S3.1、求取各时频点的信号导向矢量a(t,f);
步骤S3.2、根据期望信号的标称导向矢量a0(f)及给定的误差范围,与a(t,f)做比较后将时频域分成期望信号区域Ωs和干扰噪声区域Ωi+n;具体如下:
其中Δ为用于衡量各时频点导向矢量ai(t,f)和标称导向矢量a0(f)的相似度的参数,a(t,f)为信号导向矢量;
定义误差导向矢量ae(f)为:
其中τe为误差延时矢量;根据标称导向矢量a0(f)和误差导向矢量ae(f)设置阈值如下:
对时频域进行如下划分:
步骤S3.3、设定噪声门限,根据门限大小去除Ωi+n中的噪声时频点得到期望信号区域Ωs和干扰噪声区域Ωi;计算Ωs和Ωi中各频点处的信号延时,根据延时对Ωi中的干扰信号进行分离;
设置噪声门限Δ2,在每个时间切片上根据噪声门限Δ2剔除噪声点;对某一时间切片(ts,f)∈Ωi:
满足上式的时频点(ts,fs)即为噪声,从Ωi中剔除;在Ωs和Ωi中求出各时频点处信号的延时如下:
将Ωi中各时频点处的τi(t,f)根据大小分成J类,由于信号重叠部分较少,所以包含时频点点数较多的K类区域为各干扰信号所在区域剩余部分为信号重叠区域Ωo。分别将Ωs和中的τi(t,f)求和取平均,即可得到期望信号延时估计值τs和各干扰信号延时估计值τk;根据信号延时的特性,根据τk的大小将Ωi中的干扰信号进行分离。
步骤S4、分别计算期望信号和各干扰信号的导向矢量a(f);重构各频点的干扰加噪声协方差(IPNC)矩阵;
根据步骤S3获取的期望信号延时估计值τs和干扰信号的延时估计值τk计算期望信号和干扰信号的导向矢量估计值如下:
利用时频域上的频率切片计算各频点的协方差矩阵如下:
其中N为频率切片长度;利用Rs(f)计算Capon空间谱:
则干扰信号的协方差矩阵如下所示:
因此噪声协方差矩阵Rn(f):
Rn(f)=σn 2I
其中I为单位矩阵;σn 2为噪声功率估计值可由Rs(f)的最小特征值得到,对Rs(f)做特征值分解有:
其中λi和ui分别表示第i个特征值及该特征值对应的特征向量,则有σn 2=λmin;
步骤S5、通过MVDR算法计算各频点的自适应权矢量ω(f)如下所示:
其中Ri+n -1(f)为步骤S4中求得的干扰加噪声协方差矩阵的逆矩阵;as(f)为期望信号的导向矢量;[·]H为共轭转置算子。
下面提供一份仿真实验进一步说明本发明提供的抗干扰方法的有效性,选取仿真平台为Matlab平台,仿真实验中天线阵元的分布图如图2所示,仿真参数如下表1所示:
表1仿真实验参数表
图3是接收信号的时频谱,图4a是分离出的期望信号,图4b、图4c是分离出的干扰信号。图5a是根据自适应权值画出的二维方向图。可以看到分别在两个干扰来向的对应频率范围内产生了零陷。图5b是二维方向图在中心频率处的剖面图,可以看到在干扰来向上产生了较深的零陷的同时保证了期望信号未在其真实来向上被当作干扰信号抑制掉。
图6是入射期望信号和自适应滤波输出结果的对比图。可以看出本发明提出的方法具有良好的鲁棒性,能够有效抑制宽带干扰信号取得良好的抗干扰效果。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种基于短时傅里叶变换的空频自适应处理抗干扰方法,其特征在于,包括如下步骤:
步骤S1、获取阵列接收信号S;所述阵列接收信号S为大小M*P的样本矩阵;M是阵列天线的阵元个数,P为采样点数;
步骤S2、采用短时傅里叶变换将传感器阵列接收到的时域信号转换到时频域,即对阵列接收信号S的各行数据分别进行短时傅里叶变换;
步骤S3、将时频域划分为期望信号区域Ωs和干扰噪声区域Ωi+n;去除区域中的噪声点;在去除噪声后的干扰区域Ωi中分离各个干扰信号;具体地,
步骤S3.1、求取各时频点的信号导向矢量a(t,f);
步骤S3.2、根据期望信号的标称导向矢量a0(f)及给定的误差范围,与a(t,f)做比较后将时频域分成期望信号区域Ωs和干扰噪声区域Ωi+n;具体如下:
其中△为用于衡量各时频点导向矢量ai(t,f)和标称导向矢量a0(f)的相似度的参数,a(t,f)为信号导向矢量;
定义误差导向矢量ae(f)为:
其中τe为误差延时矢量;根据标称导向矢量a0(f)和误差导向矢量ae(f)设置阈值如下:
对时频域进行如下划分:
步骤S3.3、设定噪声门限,根据门限大小去除Ωi+n中的噪声时频点得到期望信号区域Ωs和干扰噪声区域Ωi;计算Ωs和Ωi中各频点处的信号延时,根据延时对Ωi中的干扰信号进行分离;
设置噪声门限△2,在每个时间切片上根据噪声门限△2剔除噪声点;对某一时间切片(ts,f)∈Ωi:
其中X表示采用短时傅里叶变换后的转换信号;满足上式的时频点(ts,fs)即为噪声点,从Ωi中剔除;在Ωs和Ωi中求出各时频点处信号的延时如下:
将Ωi中各时频点处的τi(t,f)根据大小分成J类,由于信号重叠部分较少,所以包含时频点点数较多的K类区域为各干扰信号所在区域剩余部分为信号重叠区域Ωo;分别将Ωs和中的τi(t,f)求和取平均,即可得到期望信号延时估计值τs和各干扰信号延时估计值τk;
步骤S4、分别计算期望信号和各干扰信号的导向矢量a(f);重构各频点的干扰加噪声协方差矩阵;具体地,根据步骤S3获取的期望信号延时估计值τs和干扰信号的延时估计值τk计算期望信号和干扰信号的导向矢量估计值如下:
利用时频域上的频率切片计算各频点的协方差矩阵如下:
其中N为频率切片长度;利用Rs(f)计算Capon空间谱:
则干扰信号的协方差矩阵如下所示:
因此噪声协方差矩阵Rn(f):
Rn(f)=σn 2I
其中I为单位矩阵;σn 2为噪声功率估计值可由Rs(f)的最小特征值得到,对Rs(f)做特征值分解有:
其中λi和ui分别表示第i个特征值及该特征值对应的特征向量,则有σn 2=λmin;
其中K代表干扰信号的个数;
步骤S5、通过MVDR算法计算各频点的自适应权矢量ω(f)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010762611.5A CN112073076B (zh) | 2020-07-31 | 2020-07-31 | 一种基于短时傅里叶变换的空频自适应处理抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010762611.5A CN112073076B (zh) | 2020-07-31 | 2020-07-31 | 一种基于短时傅里叶变换的空频自适应处理抗干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112073076A CN112073076A (zh) | 2020-12-11 |
CN112073076B true CN112073076B (zh) | 2021-10-22 |
Family
ID=73657657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010762611.5A Active CN112073076B (zh) | 2020-07-31 | 2020-07-31 | 一种基于短时傅里叶变换的空频自适应处理抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112073076B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944757A (zh) * | 2012-11-23 | 2013-02-27 | 中国人民解放军国防科学技术大学 | 一种抗环境干扰的电磁辐射发射现场测试方法 |
CN107167778A (zh) * | 2017-07-03 | 2017-09-15 | 电子科技大学 | 一种基于协方差矩阵重构和导向矢量估计的稳健波束形成方法 |
CN109743078A (zh) * | 2019-01-11 | 2019-05-10 | 河南宏泰控飞信息技术有限公司 | 基于阵列天线的无人机跳频遥控信号侦收方法与设备 |
CN110113085A (zh) * | 2019-05-23 | 2019-08-09 | 中国人民解放军战略支援部队信息工程大学 | 一种基于协方差矩阵重构的波束形成方法及*** |
CN111245490A (zh) * | 2019-11-08 | 2020-06-05 | 潘峰 | 宽带信号提取方法、装置及电子设备 |
-
2020
- 2020-07-31 CN CN202010762611.5A patent/CN112073076B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944757A (zh) * | 2012-11-23 | 2013-02-27 | 中国人民解放军国防科学技术大学 | 一种抗环境干扰的电磁辐射发射现场测试方法 |
CN107167778A (zh) * | 2017-07-03 | 2017-09-15 | 电子科技大学 | 一种基于协方差矩阵重构和导向矢量估计的稳健波束形成方法 |
CN109743078A (zh) * | 2019-01-11 | 2019-05-10 | 河南宏泰控飞信息技术有限公司 | 基于阵列天线的无人机跳频遥控信号侦收方法与设备 |
CN110113085A (zh) * | 2019-05-23 | 2019-08-09 | 中国人民解放军战略支援部队信息工程大学 | 一种基于协方差矩阵重构的波束形成方法及*** |
CN111245490A (zh) * | 2019-11-08 | 2020-06-05 | 潘峰 | 宽带信号提取方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112073076A (zh) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mestre et al. | Finite sample size effect on minimum variance beamformers: Optimum diagonal loading factor for large arrays | |
CN105302936B (zh) | 基于相关计算和协方差矩阵重构的自适应波束形成方法 | |
Boashash et al. | Robust multisensor time–frequency signal processing: A tutorial review with illustrations of performance enhancement in selected application areas | |
CN109490850B (zh) | 主瓣干扰下宽带阵列自适应波束形成方法 | |
CN106506430B (zh) | 一种基于压缩感知技术的补偿峰均比非线性失真的新算法 | |
EP1932102A2 (en) | A method and apparatus for blind source separation | |
US9123348B2 (en) | Sound processing device | |
US20110149719A1 (en) | Method for separating blind signal and apparatus for performing the same | |
CN111273237B (zh) | 基于空域矩阵滤波和干扰对消的强干扰抑制方法 | |
Reddy et al. | Unambiguous speech DOA estimation under spatial aliasing conditions | |
CN110320535A (zh) | 一种基于小波包变换与空域处理的卫星导航接收机多干扰抑制方法 | |
CN113537102B (zh) | 一种微震信号的特征提取方法 | |
CN115236594B (zh) | 一种适用于极地脉冲噪声环境的冰下声源定位方法 | |
CN115001517A (zh) | 一种无线通信接收噪声抑制方法 | |
CN108845294B (zh) | 基于子带线性多约束最小方差准则的宽带发射自适应波束形成方法 | |
CN112073076B (zh) | 一种基于短时傅里叶变换的空频自适应处理抗干扰方法 | |
CN105610748B (zh) | 一种频率分段的通道均衡方法 | |
CN113805139A (zh) | 基于聚焦变换的宽带信号稀疏表示波达方向估计方法 | |
CN107255809B (zh) | 一种基于宽带聚焦矩阵的阻塞阵波束形成方法 | |
Gershman et al. | Coherent wideband DOA estimation of multiple FM signals using spatial time-frequency distributions | |
CN113625222B (zh) | 一种基于时延差范数加权的宽带信号高分辨方位估计方法 | |
CN115808659A (zh) | 基于低复杂度不确定集积分的稳健波束形成方法及*** | |
CN113422629B (zh) | 一种协方差矩阵重构自适应波束形成方法及*** | |
CN114647931A (zh) | 一种基于期望信号消除和空间谱估计的稳健波束形成方法 | |
CN114047481A (zh) | 一种基于子空间正交性的稳健自适应波束形成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |