CN106656106B - 一种频域自适应匹配滤波器方法 - Google Patents
一种频域自适应匹配滤波器方法 Download PDFInfo
- Publication number
- CN106656106B CN106656106B CN201611215146.3A CN201611215146A CN106656106B CN 106656106 B CN106656106 B CN 106656106B CN 201611215146 A CN201611215146 A CN 201611215146A CN 106656106 B CN106656106 B CN 106656106B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- matched filter
- filter
- obtains
- signal
- 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
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0043—Adaptive algorithms
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供的是一种频域自适应匹配滤波器方法。(1)对受噪声污染的接收信号x(t)进行处理,得到受噪声污染的接收信号x(t)和经典匹配滤波器脉冲响应函数h(t)的傅里叶变换结果X(ω)和H(ω);(2)根据步骤(1)得到的傅里叶变换结果X(ω)和H(ω),在频域进行经典匹配滤波器处理,得到匹配滤波器频域输出Y(ω);(3)根据步骤(2)中的输出Y(ω),在频域分别对其实部Re[Y(ω)]和虚部Im[Y(ω)]进行自适应线谱增强处理,得到频域自适应匹配滤波器的频域输出YALE(ω);(4)求出步骤(3)中频域输出YALE(ω)的傅里叶逆变换,得到频域自适应匹配滤波器的时域输出结果yALE(τ)。本发明的方法可以在比经典匹配滤波器的检测阈更低的信噪比下更可靠地工作,使得从强干扰背景中检测弱目标成为可能,从而提高探测距离。
Description
技术领域
本发明涉及的是一种匹配滤波器技术,是一种适用于通信、雷达、声呐等电子信息系统的信号检测的方法。
背景技术
在通信、雷达、声呐等电子信息系统中,尽管有各种各样的信号检测方法,匹配滤波器仍然是最常采用的经典检测器,它是在带限白噪声干扰条件下检测确知信号(所谓理想条件)输出信噪比最大的最佳线性滤波器。
从匹配滤波器诞生至今,学界从未停止对它的研究,但由于匹配滤波器早已被证明是理想条件下的最佳滤波器,人们对于它的改进多是针对特定干扰背景下的研究,然而从匹配滤波器本身的特性来改进其性能的研究却鲜有人问津。是否可以在理想条件下寻求到比最佳滤波器更佳的处理器呢?在线性运算及定常滤波范围内,匹配滤波器无疑是最佳的,不可能找到更佳的处理器。然而跳出定常滤波器的范畴,却有可能在自适应领域寻求到比匹配滤波器检测能力更强的检测方法。
在理想条件下,比最佳处理器更佳,这在检测领域无疑是重大的挑战,在主动声呐、雷达及通信中都将有重要的应用价值,为工程应用提供性能更高、更可靠的检测技术。
发明内容
本发明的目的是提供一种性能更强,适用于通信、雷达、声呐等电子信息系统的频域自适应匹配滤波器方法。
本发明的目的是这样实现的:
(1)对受噪声污染的接收信号x(t)进行处理,得到受噪声污染的接收信号x(t)和经典匹配滤波器脉冲响应函数h(t)的傅里叶变换结果X(ω)和H(ω);
(2)根据步骤(1)得到的傅里叶变换结果X(ω)和H(ω),在频域进行经典匹配滤波器处理,得到匹配滤波器频域输出Y(ω);
(3)根据步骤(2)中的输出Y(ω),在频域分别对其实部Re[Y(ω)]和虚部Im[Y(ω)]进行自适应线谱增强处理,得到频域自适应匹配滤波器的频域输出YALE(ω);
(4)求出步骤(3)中频域输出YALE(ω)的傅里叶逆变换,得到频域自适应匹配滤波器的时域输出结果yALE(τ)。
本发明还可以包括:
1、所述受噪声污染的接收信号为带限白色噪声干扰条件下的除到达时间及幅度未知外确知的宽带线性调频信号,x(t)表示为:
x(t)=s(t)+n(t)
其中s(t)为确知宽带线性调频信号,n(t)为零均值平稳加性带限白噪声,t表示时间,ω表示角频率。
2、匹配滤波器频域输出Y(ω)表示为:
Y(ω)=X(ω)·H(ω)
其中H(ω)是匹配滤波器的传输函数,为目标信号s(t)频谱的复共轭,若匹配滤波器输出信号瞬时功率信噪比峰值位置为t0时刻,则
其中S(ω)为目标信号s(t)的傅里叶变换。
与现有技术相比,本发明的有益效果是:
本发明结合频域自适应线谱增强技术,公开了一种频域自适应匹配滤波器方法,使匹配滤波器的检测性能得到了明显改善,与经典匹配滤波器相比,本发明方法处理增益显著提高,检测阈降低约5dB,接收机工作特性曲线明显优于匹配滤波器,从而证明了频域自适应匹配滤波器各个方面的性能都显著优于经典匹配滤波器。因此应用本发明的方法可以在比经典匹配滤波器的检测阈更低的信噪比下更可靠地工作,使得从强干扰背景中检测弱目标成为可能,从而提高探测距离。
附图说明
图1为本发明方法的原理框图;
图2为自适应线谱增强器的原理框图;
图3a为信噪比-10dB时采用经典匹配滤波器的归一化处理结果图;图3b为信噪比-10dB时采用本发明方法的归一化处理结果图;
图4为采用经典匹配滤波器和本发明方法的信噪比增益曲线对比图;
图5为虚警概率千分之一时采用经典匹配滤波器和本发明方法的检测概率曲线对比图;
图6为信噪比-16dB时采用经典匹配滤波器和本发明方法的ROC曲线对比图。
具体实施方式
下面结合附图对本发明进一步说明。
结合图1具体说明本发明方法的实施流程:
(1)分别求出受噪声污染的接收信号x(t)和经典匹配滤波器脉冲响应函数h(t)的傅里叶变换X(ω)和H(ω),接收信号x(t)为带限白色噪声干扰条件下的除到达时间及幅度未知外确知的宽带线性调频信号,x(t)可表示为:
x(t)=s(t)+n(t)
其中s(t)为确知宽带线性调频信号,n(t)为零均值平稳加性带限白噪声,t表示时间,ω表示角频率;
(2)根据步骤(1)中的傅里叶变换结果X(ω)和H(ω),在频域进行经典匹配滤波器处理,得到匹配滤波器频域输出Y(ω),Y(ω)的计算方法为:
Y(ω)=X(ω)·H(ω)
其中H(ω)是匹配滤波器的传输函数,为目标信号s(t)频谱的复共轭,若匹配滤波器输出信号瞬时功率信噪比峰值位置为t0时刻,则
其中S(ω)为目标信号s(t)的傅里叶变换;
(3)根据步骤(2)中的输出Y(ω),在频域分别对其实部Re[Y(ω)]和虚部Im[Y(ω)]进行自适应线谱增强处理,得到频域自适应匹配滤波器的频域输出YALE(ω);
(4)求出步骤(3)中频域输出YALE(ω)的傅里叶逆变换,得到频域自适应匹配滤波器的时域输出结果yALE(τ)。
步骤(3)中,频域自适应线谱增强处理原理参见图2,Δ=mτ0为一延迟线,τ0为采样周期,共延迟m个采样周期,Δ应使得延迟后信号中的噪声分量和x(ω)中的互不相关,同时,Δ应使得延迟后有用信号分量和x(ω)中的仍然相关。W(ω)为频域自适应线谱增强器的权系数矢量,采用LMS算法迭代求出,ε(ω)为估计误差,用于调节W(ω),使其达到最佳。
通过与经典匹配滤波器处理效果对比,说明本发明方法所取得的有益效果。
信噪比-10dB时采用经典匹配滤波器和本发明方法的归一化处理结果图如图3a-图3b所示,从图3a-图3b的对比中可以看出采用本发明方法处理后的输出信噪比明显高于匹配滤波器,在输入信噪比为-10dB时本发明方法与匹配滤波器相比取得了更高的处理增益。
图4为采用经典匹配滤波器和本发明方法的信噪比增益曲线对比图,每个输入信噪比下样本数为1000个,取输出信噪比均值与输入信噪比的差值为系统增益。图示输入信噪比范围内,本发明方法的系统增益曲线显著高于经典匹配滤波器。
图5为虚警概率千分之一时采用经典匹配滤波器和本发明方法的检测概率曲线对比图,从图中结果可知,在要求检测概率一定的情况下,本发明方法的检测阈比匹配滤波器降低约5dB,因此可以在信噪比更低的条件下工作。在一定背景干扰下,本发明方法成功检测信号的能力比经典匹配滤波器更强。
图6为信噪比-16dB时采用经典匹配滤波器和本发明方法的ROC曲线对比图,从图中可以看出,本发明方法的ROC曲线明显优于经典匹配滤波器。
本发明结合频域自适应线谱增强技术,公开了一种频域自适应匹配滤波器方法,使经典匹配滤波器的各项检测性能得到了明显改善,应用本发明的方法可以在比经典匹配滤波器的检测阈更低的信噪比下更可靠地工作,并且本发明方法算法简单,容易实现,有望广泛应用于通信、雷达、声呐等电子信息系统的信号检测领域。
Claims (3)
1.一种频域自适应匹配滤波器方法,其特征是:
(1)对受噪声污染的接收信号x(t)进行处理,得到受噪声污染的接收信号x(t)和经典匹配滤波器脉冲响应函数h(t)的傅里叶变换结果X(ω)和H(ω),t表示时间,ω表示角频率;
(2)根据步骤(1)得到的傅里叶变换结果X(ω)和H(ω),在频域进行经典匹配滤波器处理,得到匹配滤波器频域输出Y(ω);
(3)根据步骤(2)中的输出Y(ω),在频域分别对其实部Re[Y(ω)]和虚部Im[Y(ω)]进行自适应线谱增强处理,得到频域自适应匹配滤波器的频域输出YALE(ω);
(4)求出步骤(3)中频域输出YALE(ω)的傅里叶逆变换,得到频域自适应匹配滤波器的时域输出结果yALE(t)。
2.根据权利要求1所述的频域自适应匹配滤波器方法,其特征是所述受噪声污染的接收信号为带限白色噪声干扰条件下的除到达时间及幅度未知外确知的宽带线性调频信号,x(t)表示为:
x(t)=s(t)+n(t)
其中s(t)为确知宽带线性调频信号,n(t)为零均值平稳加性带限白噪声。
3.根据权利要求1或2所述的频域自适应匹配滤波器方法,其特征是匹配滤波器频域输出Y(ω)表示为:
Y(ω)=X(ω)·H(ω)
其中H(ω)是匹配滤波器的传输函数,为目标信号s(t)频谱的复共轭,若匹配滤波器输出信号瞬时功率信噪比峰值位置为t0时刻,则
其中S(ω)为目标信号s(t)的傅里叶变换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611215146.3A CN106656106B (zh) | 2016-12-26 | 2016-12-26 | 一种频域自适应匹配滤波器方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611215146.3A CN106656106B (zh) | 2016-12-26 | 2016-12-26 | 一种频域自适应匹配滤波器方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106656106A CN106656106A (zh) | 2017-05-10 |
CN106656106B true CN106656106B (zh) | 2019-04-23 |
Family
ID=58827788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611215146.3A Active CN106656106B (zh) | 2016-12-26 | 2016-12-26 | 一种频域自适应匹配滤波器方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106656106B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108535863B (zh) | 2018-03-30 | 2020-10-27 | 武汉华星光电技术有限公司 | 一种盖板的处理方法、控制器、盖板处理装置及存储介质 |
CN112462352B (zh) * | 2020-10-30 | 2022-10-18 | 哈尔滨工程大学 | 一种适用于低信噪比条件下的线谱增强方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1383612A (zh) * | 2000-05-16 | 2002-12-04 | 株式会社鹰山 | 匹配滤波器电路 |
CN101588329A (zh) * | 2008-05-19 | 2009-11-25 | 北京大学深圳研究生院 | 扩频通信多相位匹配滤波器捕获电路及其设计方法 |
CN103744066A (zh) * | 2014-01-09 | 2014-04-23 | 西安电子科技大学 | 数字正交鉴相与匹配滤波的优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2512395A (en) * | 2013-03-28 | 2014-10-01 | Sony Corp | Transmitter and method of transmitting |
-
2016
- 2016-12-26 CN CN201611215146.3A patent/CN106656106B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1383612A (zh) * | 2000-05-16 | 2002-12-04 | 株式会社鹰山 | 匹配滤波器电路 |
CN101588329A (zh) * | 2008-05-19 | 2009-11-25 | 北京大学深圳研究生院 | 扩频通信多相位匹配滤波器捕获电路及其设计方法 |
CN103744066A (zh) * | 2014-01-09 | 2014-04-23 | 西安电子科技大学 | 数字正交鉴相与匹配滤波的优化方法 |
Non-Patent Citations (3)
Title |
---|
Digital signal processing;Z Hussain 等;《Springer》;20110101;第44-48页 |
The Compressive Matched Filter;Mark A.Davenport等;《Rice University Department of Electrical and Copmputer Engineering Technical Report TREE 0610》;20061101;第1-16页 |
一种新的匹配滤波器工程实现方法;郝宗杰等;《电声技术》;20090215;第33卷(第5期);第45-47页 |
Also Published As
Publication number | Publication date |
---|---|
CN106656106A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772457B (zh) | 一种基于空时频架构的卫星导航抗干扰方法 | |
US9651670B2 (en) | Method and apparatus for adjusting threshold weight of frequency domain to improve anti-jamming performance and eliminate jamming in GPS system | |
CA2954635C (en) | Improved signal detection and characterization | |
CN106656106B (zh) | 一种频域自适应匹配滤波器方法 | |
CN109725297B (zh) | 一种基于回波极化度的有源转发式干扰识别方法 | |
US7295145B2 (en) | Selective-sampling receiver | |
Meller et al. | Processing of noise radar waveforms using block least mean squares algorithm | |
CN107479050B (zh) | 基于对称谱特性和次对称特性的目标检测方法和装置 | |
KR101908455B1 (ko) | 저전력 주파수 변조 연속파 레이더 시스템 및 그 제어 방법 | |
CN110568415B (zh) | 混合高斯模型下基于Arctan函数的信号检测方法 | |
CN112748429B (zh) | 一种快速噪声对消滤波方法 | |
CN106533394B (zh) | 一种基于自适应滤波器幅频响应的高精度频率估计方法 | |
CN117805756A (zh) | 一种杂波中低虚警率噪声干扰识别方法及系统 | |
Shin et al. | Weak radar signal detection based on variable band selection | |
Zhou et al. | The Influence of Automatic Gain Control on Narrowband Frequency Domain GPS Anti-Jamming Receiver | |
KR101303192B1 (ko) | 수동 소나시스템 및 이의 데몬 처리 향상기법 | |
CN107976659A (zh) | 基于S2method时频分析的目标极点特征提取与检测方法 | |
CN103901421A (zh) | 基于对角减载的水声阵列smi-mvdr空间谱估计方法 | |
CN104793197A (zh) | 基于ifft频谱相除法和梯度自适应格型滤波的直达波抑制方法 | |
Xia et al. | Weak signal detection based on FRFT-domain adaptive filtering algorithm | |
Bvrley et al. | A receiver for multicarrier logarithmic warped frequency domain code waveforms | |
US11698436B2 (en) | Radar devices and methods using multicarrier logarithmic warped frequency domain code waveforms | |
Erdogan et al. | Digital chirp rate adaptation for optimum FMCW interception using Wigner-Hough transform | |
CN111948629A (zh) | 一种高稳健性大多普勒单频脉冲信号检测方法 | |
CN116973835B (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 |