一种线性调频信号参数分析的综合方法

合集下载

线性调频信号的时频分析研究

线性调频信号的时频分析研究随着通信技术的发展，线性调频信号（Linear Frequency Modulation，LFM）在通信系统中得到了广泛的应用。

线性调频信号是一种在一段时间内频率线性变化的信号，其具有宽带、抗多径衰落、抗高噪声等特点，因此适用于高分辨率雷达、超声定位、地震勘探等领域。

为了更好地理解和设计线性调频信号的应用系统，对其进行时频分析研究是非常重要的。

时频分析是一种将信号在时间和频率域上进行联合分析的方法，可以提供关于信号特性的更详细的信息。

对于线性调频信号而言，时频分析可以帮助我们获得信号的调频特性和调制参数。

下面将介绍几种常见的时频分析方法，以及它们在线性调频信号研究中的应用。

STFT是一种将信号在时间和频率上进行分析的方法，它通过将信号分成多个小时间窗口，并对每个窗口进行傅里叶变换，得到该窗口内信号的频谱信息。

STFT可以提供线性调频信号的瞬时频率信息，帮助我们理解信号的调频特性。

2. Wigner-Ville分布（Wigner-Ville Distribution，WVD）WVD是一种采用时频联合分析的方法，它通过计算信号的瞬时相位和瞬时幅度，得到信号在时频上的分布。

WVD可以提供线性调频信号的瞬时频率和瞬时频谱信息，有助于我们研究信号的调频参数和调频性质。

3. 希尔伯特-黄变换（Hilbert-Huang Transform，HHT）此外，还有一些其他的时频分析方法，如连续小波变换（Continuous Wavelet Transform，CWT）、自适应滤波器（Adaptive Filter），它们在线性调频信号研究中也有一定的应用。

通过将这些方法相互结合，可以更好地理解线性调频信号的时频特性和调制参数。

在线性调频信号的时频分析研究中，我们可以分析信号的频谱特性、瞬时频率变化、调制参数等。

通过这些分析，我们可以了解信号是否具有带宽限制特性、频率变化规律，以及在特定调制参数下，信号的传输性能如何。

一种新的线性调频脉冲信号参数估计算法

ｃｈｉｒｐｒａｄａｒｓｉｇｎａｌｓ．ＡｃｈｉｒｐｓｉｇｎａｌｐａｒａｍｅｔｅｒｅｓｔｉｍａｔｉｏｎｍｅｔｈｏｄｂａｓｅｄｏｎｐｅｒｉｏｄＷＶＨＴ（ＰＷＶＨＴ）ｉｓ
ＷＶＨＴｉｍａｇｅｗｈｅｎｔｈｅｔｒａｄｉｔｉｏｎａｌＷｉｇｎｅｒ — Ｖｉｌｌｅ－Ｈｏｕｇｈｔｒａｎｓｆｏｒｍ（ＷＶＨＴ）ｉｓｕｓｅｄｔｏｄｅｔｅｃｔａｎｄｅｓｔｉｍａｔｅ
Ａｂｓｔｒａｃｔ：ＴｈｅｏｕｔｐｕｔＳＮＲａｔｔｈｅｔｒｕｅｐａｒａｍｅｔｅｒｖａｌｕｅｓｄｏｅｓｎｏｔｉｎｃｒｅａｓｅａｎｄｔｈｅｍｕｌｔｉｐｌｅｐｅａｋｓｏｃｃｕｒｉｎ
样本信号信噪比的关系，采用ＭＡＴＬＡＢ对ＰＷＶＨＴ进行线性调频脉冲信号参数估计的误差进行了
仿真。仿真结果表明：在对多线性调频脉冲分量信号进行参数估计时，ＰＷＶＨＴ的性能更优。
关键词：信息处理技术；线性调频脉冲信号；参数估计；Ｗｉｇｎｅｒ．Ｖｉｌｌｅ．Ｈｏｕｇｈ变换；周期Ｗｉｇｎｅｒ —
ｓｕｐｅｒｉｏｒｔｏＷＶＨＴｆｏｒｅｓｔｉｍａｔｉｎｇａｃｈｉｒｐｓｉｇｎａｌｔｈａｔｉｓｃｏｍｐｏｓｅｄｏｆｍｕｌｔｉｐｌｅｐｕｌｓｅｓ．

一种新的线性调频脉冲信号参数估计算法

一种新的线性调频脉冲信号参数估计算法陈磊;陈殿仁;刘颖【摘要】在采用传统Wigner-Ville-Hough变换(WVHT)进行线性调频脉冲信号检测和参数估计的过程中,当观测数据具有多个线性调频脉冲周期时,在估计参数真值处的输出信噪比不能因为观察时间的增加而变好,反而会在时频变换图中出现多个峰值,影响信号参数的估计性能.为此,研究了一种基于周期WVHT(PWVHT)的线性调频脉冲信号参数估计方法,给出通用信号的PWVHT公式并推导了线性调频脉冲信号的PWVHT公式,分析PWVHT输出信号的信噪比与观测时间和样本信号信噪比的关系,采用MATLAB对PWVHT进行线性调频脉冲信号参数估计的误差进行了仿真.仿真结果表明:在对多线性调频脉冲分量信号进行参数估计时,PWVHT的性能更优.【期刊名称】《兵工学报》【年(卷),期】2014(035)002【总页数】7页(P207-213)【关键词】信息处理技术;线性调频脉冲信号;参数估计;Wigner-Ville-Hough变换;周期Wigner-Ville-Hough变换【作者】陈磊;陈殿仁;刘颖【作者单位】长春理工大学电信学院,吉林长春130022;长春理工大学电信学院,吉林长春130022;长春理工大学电信学院,吉林长春130022【正文语种】中文【中图分类】TN974近年来,在雷达信号处理领域,线性调频(LFM)脉冲信号参数检测和估计受到广泛关注[1]。

多种信号处理技术被应用于LFM脉冲信号的的参数检测和估计算法中,包含功率谱估计、Wigner-Ville分布、Choi-Williams分布等方法的非线性时频变换技术,用来产生信号的时频图像,进而分析信号瞬时频率随时间的变换情况[2]。

邓兵等[3]分析了基于分数阶傅里叶变换的LFM脉冲时延估计特性。

刘锋等[4]提出了一种基于联合Wigner-Ville分布-随机Hough变换改进算法的LFM信号参数快速估计方法;陈鹏等[5]研究了一种基于离散分数阶傅立叶变换的水下动目标线性调频回波检测算法;罗洁思等[6]研究了一种基于多尺度LFM基信号稀疏分解的多分量LFM信号检测方法。

基于综合算法的多分量线性调频信号参数估计方法及系统[发明专利]

专利名称：基于综合算法的多分量线性调频信号参数估计方法及系统
专利类型：发明专利
发明人：黄理杰,钱江
申请号：CN202010653459.7
申请日：20200708
公开号：CN111766444A
公开日：
20201013
专利内容由知识产权出版社提供
摘要：本发明公开一种基于STFT、神经网络和Radon变换的多分量线性调频信号参数估计方法，首先求多分量LFM信号的短时傅里叶变换(STFT)；用神经网络判断LFM信号分量数量；然后通过统一频率峰值高度消除信号强度影响；最后通过逐次消去法估计每个LFM分量的起始频率和调频斜率(用到Radon变换)；估计LFM信号的幅度。

本发明通过神经网络判断LFM信号分量数量以及通过统一频率峰值高度消除信号强度影响，可以直接从STFT图判断LFM信号分量数量，而不是从Radon平面判断LFM信号分量数量，判断更为准确；不会受到交叉项影响，也不会受到不同信号分量幅度不同或同一信号分量幅度随时间变化的影响。

申请人：电子科技大学
地址：611731 四川省成都市高新区(西区)西源大道2006号
国籍：CN
代理机构：重庆航图知识产权代理事务所(普通合伙)
代理人：胡小龙
更多信息请下载全文后查看。

线性调频连续波信号检测与参数估计算法

线性调频连续波信号检测与参数估计算法朱文涛;苏涛;杨涛;郑纪彬;张龙【摘要】针对长时间积累的方法较难在线性调频连续波信号的检测和参数估计中应用的问题，该文基于归一化变窗长相干平均法，提出一种联合帧间相关法与循环平稳法的线性调频连续波信号检测与参数估计算法。

该方法首先利用归一化变窗长相干平均法实现噪声方差的一致性，然后在此基础上利用帧间相关法实现周期的精确估计，最后利用循环平稳法估计信号的相位参数。

该算法以较低的计算复杂度实现信号的长时间积累，在低信噪比条件下具有较好的估计性能。

仿真结果验证了该方法的有效性。

%To deal with the issue that the long time accumulation method is hard to be applied to the detection and parameter estimation of Linear Frequency Modulated Continuous Wave (LFMCW) signal, a novel algorithm combining the joint inter-frame correlation method and cyclostationary method is proposed to detect LFMCW signal and estimate its parameters. Firstly, consistency of noise variance is completed by the normalized window-variant coherent average method. Then, on this basis, the period is exactly estimated with inter-frame correlation method. Finally, cyclostationary method is utilized to estimate the phase parameters of LFMCW signal. The proposed algorithm realizes long time accumulation with low computational complexity. For low Signal-to-Noise Ratio (SNR), this method still achieves good estimation performance. Simulation results verify the effectiveness of the proposed method.【期刊名称】《电子与信息学报》【年(卷),期】2014(000)003【总页数】7页(P552-558)【关键词】信号检测;参数估计;线性调频连续波信号;归一化变窗长相干平均;帧间相关;循环平稳【作者】朱文涛;苏涛;杨涛;郑纪彬;张龙【作者单位】西安电子科技大学雷达信号处理国家重点实验室西安 710071;西安电子科技大学雷达信号处理国家重点实验室西安 710071;西安电子科技大学雷达信号处理国家重点实验室西安 710071;西安电子科技大学雷达信号处理国家重点实验室西安 710071;西安电子科技大学雷达信号处理国家重点实验室西安710071【正文语种】中文【中图分类】TN957.51调频连续波(FMCW)雷达由于结构简单、体积小、距离分辨率高、无距离盲区、成本低、低功耗和低截获等优点，在军用导航、战场侦察与地面成像等领域得到越来越广泛的应用。

一种多径线性调频信号的参数估计方法

数重构出源信号，并将其和观测信号作循环互相关变换，再根据信号循环频率域上的多径信号
特征，估计出最强路径信号分量的时延和衰减因子参数，后重构最强路径分量，然最后将其从
观测信号中滤去，这样逐次分离出各径信号，而有效地抑制了强路径信号分量对弱路径信号从
ｎｎｓｒｃｎｔｕｔｄｕｉｇｔｅｅｔｔａａｅｅｄｉｃｔｒｍｅｏｓｒｉｇｓｇａｓｅｔｉｅｏｓｒｃｅｓｎｓｉｅｐｍｔｒａｓｕｏｔｂｅｎｉｌ，ｗｈｃａｈａｍｄｒｓｎｆｈｖｎｉｈｃｎｓｐｒｔａｈｐｔｉａｄｒｄｃｅｅｆｃｆｓｏｇｐｔｏｏｅｔｕｎｗａｏｏｎｎ．Ｆ — ｅａａｅｅｃａｈｓｇｌａｅｕｅｔｆｔｔｎａｃｍｐｎｎｐｅｋｃｍｐｅｔｉｎｎｈｅｏｒｈｏｎｌ，ｓｍｅｃｍｐｔｒｓｍｕａｉｎｒｉｅｌａｙｏｏｕｅｉｌｔｓａｅｇｖｎ．ｏ
（．ｅｉｉｔｙＲｐｅｅｔｉｅａｍｎ，ｅｉ０４；．ｉｉｄ１ｔａ，ｅｉ０９）１ＢｉｇＭｌｒｅｒｎｔｅＤｐｒｅｔＢｉｇ１０１２ｊｎｉａｓａｖｔｊ０ｎｒｎＲ．ｓＹｒＢｉｇ１０４ａｇｌｄｊ０ｎ
ＡｂｔａｔＡｅｍｅｈｄｆｒｅｔｔｇｐｒｍｅｅｓｆｒＬＭｙｃｃｉＣ＇Ｓｃｒｅａｉｎｔｎｆｒｃｎ・ｓｒｃ：ｎｗｔｏｓｍａｉａａｔｒＦｂｙｌＩＳ・ｏｌｔｒｓｍｏ－ｏｉｎｏｃＯ－ｏａｏｂｎｎｉｉｎｌｅａａｉｎｔｈｏｏｙｉｒｓｎｅｉｉｇｗｔｓｇａｐｒｔｃｎｌｇｐｅｅｔｄ．Ｆｒｔｈｅｃａａｔｒａａｔｒｏｍｉｅｉａｈｓｏｅｓｉ，ｔｈｒｃｅｒｍｅｅｆｓｐｓｅｔｄｓｇｌｔｎａｅｅｔｔｄｂｙｌｏｒｌｔｎｔｎｆｒｄｔｅｅｔｄｓａｓｒｃｎｔｃｅｙｔｅｅｔｔｄｒｓｍａｅｙｃｃｉｃｒｅａｉｒｓｏａｍｔｉｌｉｅｏｓｕｔｄｂｓｉｅｉｃｏａｍｎｈｉｅｎｇｒｈａｍｐｒｍｅｅ．Ａｄｔｅ．ｔｅｐｒｍｅｅｓｏｅｓｏｇｓｉｎｏｏｅｔｒｓｉｔｄｂｙｌｒｓ — ａａｔｒｓｎｎｈａａｔｒｆｔｎｅｔｇａｃｍｐｎｎｅｅｔｈｈｔｒｓｌａａｍｅｙｃｃｉｃｏｓｃ

线性调频信号的时频分析研究(优选.)

最新文件---------------- 仅供参考--------------------已改成-----------word文本 --------------------- 方便更改赠人玫瑰，手留余香。

线性调频信号的时频分析研究摘要线性调频信号是其中一类有代表性的非平稳信号，该信号广泛存在于雷达、声纳、语音、地球物理和生物医学信号处理中。

对于这类频率随时间变化的信号，传统的时间域和频率域的分析方法都不能够全面的反映信号的特征，而时频分析是分析和处理非平稳信号的有力工具。

利用时频分布，可以对各种信号进行分析、处理，提取信号在特定时间特定频率所具有的特征信息。

文中介绍了线性调频信号的定义及特性，描述了短时傅里叶变换，Wigner—Ville 分布，Wigner—Hough分布三种时频分析方法。

通过对时频分析方法的原理介绍,运用MATLAB 中的工具箱，对一个线性调频信号进行时频分析的MATLAB仿真。

通过对几种时频分析方法比对分析和基于MATLAB信号降噪的仿真实验，验证几种分析方法的优越性。

关键词：线性调频信号时频分析短时Chirp-Fourier变换 Wigner—Ville分布Wigner—Hough分布Linear FM signal time-frequency analysisAbstractIn modern signal processing, linear frequency modulation signal is one representative of non-stationary signals, the signal is widespread in radar, sonar, speech, and geophysics, and biomedical signal processing. Such frequency time-varying signal, the traditional time domain and frequency domain analysis methods are not able to fully reflect the characteristics of the signal, but when the frequency analysis is a powerful tool for analysis and processing of non-stationary signals. Using time-frequency distribution to analyze a variety of signal processing, extract the signal characteristics with a specific frequency at aspecific time.This paper introduces the definition and characteristics of the linear FM signal, describes the short-term Chirp-Fourier Transform, Gabor distribution ，Wigner-Ville distribution of two kinds of time-frequency analysis. By the principle of time-frequency analysis method, the use of the toolbox in MATLAB, MATLAB simulation of time-frequency analysis of a linear FM signal. By frequency analysis of several methods of analysis and MATLAB-based signal to noise simulation and validation of several advantages of the method.Key words: LFM signal Time-frequency analysis Wigner-Ville distribution Discrete Chirp-Fourier transform目录1 绪论 (3)1.1 课题背景及研究意义 (3)1.2 国内外发展状况 (6)1.3本论文的主要内容 (7)2 线性调频信号 (8)2.1 线性调频信号的定义 (9)2.2线性调频信号的特点 (9)2.3 线性调频信号的仿真 (11)3 线性调频信号的时频分析方法研究 (15)3.1时频分析的定义 (15)3.2时频分析基本思想 (16)3.3 时频分析方法的介绍和仿真 (16)3.3.1 短时傅里叶变换 (16)3.3.2 Winger—Ville分布变换结果 (24)3.3.3 W-H变换结果 (31)4 结论 (34)附录 (36)参考文献 (42)致谢 (44)1 绪论本章介绍了本文的研究背景和意义，概述了线性调频信号和时频分析理论及应用的研究进展和现状，给出了全文的内容安排。

一种基于HHT的线性调频信号分析与参数估计

第39卷第3期应用科技V ol.39, No.3 2012年6月 Applied Science and Technology Jun. 2012 doi:10.3969/j.issn.1009-671X.201112015基于HHT的线性调频信号分析与参数估计甄浩川1，申丽然2，张晓林11. 哈尔滨工程大学信息与通信工程学院, 黑龙江哈尔滨 1500012. 大连海事大学信息科学技术学院, 辽宁大连 116026摘要：基于Hilbert-Huang变换（HHT），提出一种有效的线性调频（LFM）信号的分析与参数估计的方法. 首先对LFM信号进行HHT得到其Hilbert谱，并根据能量的准则提取其中的主成分，完成LFM信号分析. 对利用能量型主成分提取法得到的LFM信号主成分进行最小二乘直线的拟合，计算直线的斜率与截距，得到LFM信号的参数的估计值. 实验证明，能量型主成分提取方法，在较高的信噪比范围内具有一定的LFM信号估计效果.关键词：线性调频信号; 经验模态分解; Hilbert-Huang变换; 主成分提取; 本征模态函数; 参数估计中图分类号：X75 文献标志码：A 文章编号：1009-671X（2012）03-0045-06HHT-Based linear frequency modulated signal analysisand parameter estimationZHEN Haochuan1, SHEN Liran2, ZHANG Xiaolin11. College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China2. College of Information Science and Technology, Dalian Maritime University, Dalian 116026, ChinaAbstract:Based on Hilbert-Huang transform (HHT), an effective method is proposed to analyze the linear frequency modulated (LFM) signal and estimate the parameters. The signal was processed through HHT first to get its Hilbert spectrum. According to a criterion of energy, the principal components were extracted for analysis of LFM signal. The principal components of LFM signal obtained by the energy-oriented principal component extraction method were processed by the least-squares line fitting method, then the gradient and intercept of the line were computed to get the estimated value of the parameters of the LFM signal. Experiments prove that the energy-oriented principal component extraction has a good effect on LFM signal estimation within a higher SNR range.Keywords:LFM signal; empirical mode decomposition（EMD）; HHT; principal component extraction; intrinsic mode function（IMF）; parameter estimation线性调频信号（LFM）作为一种常见的非平稳信号，具有较大的时宽和带宽积、频率随着时间线性变化等特点，广泛应用于雷达、声纳和地质探测等领域. 因此，LFM的参数估计是一项同时具有理论价值与实用价值的研究课题，向来备受关注. 有很多方法可以对信号的瞬时频率进行估计，但是均有一定的局限性. 例如相位差分法只适合于几乎没有噪声的情况；用线性回归的方法估计得到的瞬时频率[1]在信噪比（SNR）较低的时候效果也会恶化甚至无效.Hilbert-Huang变换[2]（HHT）是由N.E.Huang于1998年提出的一种适合于非平稳非线性信号处理的收稿日期：2011-12-08.作者简介：甄浩川（1986-），男，硕士研究生，主要研究方向：语音信号处理. E-mail: zhenhaochuan@.新理论. 其核心在于经验模态分解（EMD）部分，无需任何先验知识，完全自适应的根据信号本身的局部特征，得到本征模态函数（IMF），再对IMF进行Hilbert变换，即可得到被分解信号的Hilbert谱. IMF 的获得保证了瞬时频率的物理意义及完整性，也保证了Hilbert谱的实用性. HHT自出现以来，就受到了学术界的广泛关注，从其最早使用于地震数据和海洋潮汐数据的分析，扩展到现在的语音增强[3]、语音端点检测[4]、合成孔径雷达[5]和故障检测[6]等领域.文中对LFM信号的HHT分析，从提取主成分的角度对文献[7]所提方法称之为直接型主成分提取法. 通过进行深入研究和分析，进而提出能量型主成分提取法，提取LFM信号Hilbert谱主成分，实现对LFM·46·应用科技第39卷信号的分析. 能量型主成分提取法可以作为一种有效的LFM 信号参数估计方法使用，且其在各种信噪比环境下均具有较好的LFM 信号参数估计效果.1 Hilbert-Huang 变换HHT 分为2大步骤：EMD 和Hilbert 变换. 1.1 经验模态分解（EMD ）EMD 的处理过程是简单的，它最主要的目的是将信号分解成一些平稳的、线性的数据序列. 每一个序列被称为一个本征模态函数（IMF ）. 作为HHT 重要且创新的部分，EMD 是服务于Hilbert 变换的，即IMF 是最适合进行Hilbert 变换的. 因此关于IMF 的定义，需要满足以下2个条件：1）在全部数据集中，极值点的数目与过零点的数目必须相等或者至多相差一个；2）在任意点，上下包络的均值为零，这里的上下包络分别是由局部最大值和局部最小值定义的.给定一个信号()s t ，通过IMF 的定义，将EMD 算法表述如下：1）提取信号()s t 的全部的局部极值点. 2）分别通过对全部的局部极大值点和局部极小值点使用三次样条插值函数，得到上包络线()max e t 和下包络线()min e t . 计算上下包络线的均值()()()()max min 2.m t e t e t =+ （1）3）计算1h :()1()h s t m t =−. 如果1h 满足IMF 的2个条件，就得到了1个IMF 、1h .4）如果1h 不是1个IMF ，需要将1h 视为待处理的信号，并且重复步骤1）~3）. 这个过程就叫做筛选. 筛选不会停止，除非第k 次的时候，1k h 满足了IMF 的2个条件，就得到了第1个IMF ，记做1c : 11k c h =. 实际上，严格满足IMF 定义的2个条件是非常不容易的，因为Huang 等人提出利用一种标准差SD 作为另一个判据，满足SD 的1k h 也可以称之为1个IMF ：()21(1)1201(1)()().()Tk k t k h t h t SD h t −=−⎡⎤−=⎢⎥⎢⎥⎣⎦∑ （2） 5）将1c 从数据中剔除：()11r s t c =−.6）定义1r 为新的数据，使其通过步骤1）~5）获得2c ;7）最终，全部的IMF 被筛选出来，筛选后的余项为n r ，原始信号可以表达为各IMF 与余项的和1(),ni n i X t c r ==+∑ （3）至此，信号被分解完毕.可以看出，信号在EMD 分解过程之前和之中几乎无需人为干涉，完全由数据驱动，得到的分解结果是自适应的，是反应信号内在物理意义的，尤其是IMF 保证了Hilbert 变换中瞬时频率的物理意义.1.2 Hilbert 变换对所有IMF 分别进行Hilbert 变换，求得各IMF 的瞬时频率和瞬时幅度，再联合时间轴，得到一种时频分析方式.对于给定信号()x t ，由各IMF 构成：()().Ni i x t c t ==∑1（4）对()x t 做Hilbert 变换，并将式（4）代入得 ()()()()()()()()()()()d d ππ...d πd d ...d πππd .πNi i N N N Ni i i i c x y t t t c c c t t t c c c t t t c d t t τττττττττττττττττττττττττ======−−⎛⎞ +++=⎜⎟−−−⎝⎠+++=−−− =−∑∫∫∫∫∫∫∑∑∫11212111111111（5）相当于对各IMF 分量()i c t 做Hilbert 变换；所以，各IMF 分量的解析信号为()()()j ,i i i e t c t d t=+ （6）各IMF 分量的瞬时幅度为()i a t = （7）各IMF 分量的瞬时相位为()()()arctan ,i i i d t t c t θ⎛⎞=⎜⎟⎝⎠（8）则各IMF 分量的瞬时频率为()d ;d i i t tθω=（9）因此原信号可以表示为()()j ()Re e.i N t i ix t a t θ=∑ （10）式中Re 表示取实部. 又由于瞬时频率的定义，此时可以将信号()x t 写做()()()()Re exp j d ,Ni i ix t a t t t ω=∑∫ （11）并将其看作既是时间的函数，又是瞬时频率的函数.第3期甄浩川，等：基于HHT 的线性调频信号分析与参数估计 ·47·这样构成的时频平面称为Hilbert 谱. 严格意义上讲，只有经过EMD 分解后，利用IMF 做Hilbert 变换得到的谱才叫Hilbert-Huang 谱.2 LFM 信号的Hilbert 谱主成分提取与参数估计设LFM 信号的模型为()()()20exp j2π12.x t A f t at =+⎡⎤⎣⎦ (12)式中：A 是幅度，0f 和a 分别是其调频斜率和初始频率，文中重点关注0f 和a 的估计，默认1A =. 仿真信号的采样频率为1kHz s f = ，调频斜率050Hz f = ，初始频率250Hz/s a = ，采样点数为1 000个，时间长度为1s .直接型主成分提取法是由文献[7]提出，即对LFM 信号进行EMD 分解得到众IMF ，将第一阶IMF 在Hilbert 谱上进行最小二乘直线拟合，得到LFM 信号的估计. 但文献[7]只给出了在无噪声情况下的实验结果. 鉴于此，对直接型主成分提取法做更进一步的讨论，从原理上进行分析和实验，并提出具有估计能力的能量型主成分提取法.2.1 直接型主成分提取法图1是LFM 信号的时域波形，图2是LFM 信号经过EMD 分解后得到的各IMF.图1 LFM 信号的时域波形学者Rilling 等[8]通过对EMD 算法的研究和其与高斯白噪声的相关实验，得到EMD 分解与滤波器组等价的结论. 这种滤波器组将输入信号由高频到低频，逐层分解并输出，滤波器组的参数由信号自身的时域局部表现决定，因此是自适应的信号分解方法. 可以看出，第一阶IMF 在时域的波形直观看上去与LFM 信号的时域波形有些相似之处，频率逐渐增快. 图3是LFM 信号的Hilbert 谱，图4是对LFM 信号进行HHT 后的Hilbert 谱.（a ）第1至第5阶IMF（b ）第6阶至第7阶IMF图2 LFM信号经EMD 分解后得到的IMF图3 LFM信号的Hilbert 谱图4 LFM 信号的Hilbert-Huang 谱·48·应用科技第39卷可以从LFM 信号的数学模型推导出，其瞬时频率与时间是呈线性关系的，但经过HHT 后却出现了一些LFM 原本并不含有的能量非常微弱的瞬时频率成分，这些成分都是偏离那条理论的直线的而且自身无线性可言. 关于瞬时频率，Cohen [9]在其书中提到瞬时频率的定义方法，可能会导致所求瞬时频率可能不是理论频谱中的频率甚至是负值. 经过EMD 分解得到的IMF 保证了瞬时频率非负，但是无法避免另一个问题的出现. 仔细观察图3可以发现，LFM 信号在Hilbert 谱上是有很小一部分散乱的点，但都聚集在直线附近，发散的不严重. 而图4所示经过EMD 分解后的Hilbert 谱中，散乱的点分布在更广的频谱范围. 这样的结果是因为EMD 分解过程使得Cohen 所提的现象得到了扩散. 虽然直觉上认为这是不好的现象，但是文中却恰恰是研究并利用这一特性. 将2幅Hilbert 谱对比，发现只有其中由第一阶IMF 构成的分量从频率覆盖范围看与LFM 信号本身的Hilbert 谱最为接近，且其所含能量最高，瞬时频率点也比较聚集在直线附近；因此，完全有理由相信第一阶IMF 就是使用EMD 算法得到的LFM 信号在Hilbert 谱的主成分，文中称之为直接型主成分提取法，也就是文献[7]所提方法. 这种主成分的获得是通过对时域波形和Hilbert 谱的直接观察得到的.图5 无噪声时直接型提取到的LFM 信号主成分图6 无噪声时利用直接型提取到的主成分拟合的直线图5是从纯净的LFM 中提取出来的第一阶IMF 的Hilbert 谱，将这些瞬时频率点拟合后的直线如图6所示，效果非常不错，=2496422a .ˆ、0=502147f ˆ.，与真实值的误差都非常小. 但是进一步实验发现，直接型主成分提取法几乎不具备抗噪性能.图7是SNR=10dB 时，对LFM 信号进行HHT 后的Hilbert 谱. 图8是使用直接型主成分提取法得到的LFM 的主成分的Hilbert 谱.图7 SNR=10dB 时LFM 信号的Hilbert-Huang 谱图8 SNR=10 dB 时直接型提取到的LFM 信号主成分图9 SNR=10 dB 时利用直接型提取到的主成分拟合的直线图9是拟合得到的直线，此时，估计效果恶劣，=64.6261a -ˆ，0=2848324f ˆ.，与真实值相差甚远. 将图5和8对比可以知道，在加入高斯白噪声后，即便是SNR 不低的情况下，直接型主成分提取法得到的第3期甄浩川，等：基于HHT 的线性调频信号分析与参数估计 ·49·信号主成分已经被高斯白噪声严重的污染，瞬时频率成分扩散到了更宽的频谱上，这直接导致直线拟合不准确，估计性能急剧恶化. 2.2 能量型主成分提取考虑EMD 分解的特性，即每个分量都具备与原信号相同的时间轴长度；因此，各IMF 均保留了原信号在时间轴上的信息并分别在不同频带给予表达. 观察发现，并不是每个IMF 在每一时刻的瞬时频率都落在了纯净LFM 信号在Hilbert 谱中的那条直线上，而且Huang [2]也指出各IMF 不会在同一时间位置包含相同的瞬时频率，那么可以认为在每一时刻只有一个IMF 对LFM 信号的时频性质起主要作用. 由此，提出一种以幅度最大即能量最高为判据的主成分提取法：对于同一时刻τ，选出各IMF 在时刻τ对应的瞬时幅度最大值k A ，()arg max ,j jt k A t τ== （13）并将其对应的瞬时频率k f 记录下来. 将所有这样的瞬时幅度的集合及其对应的瞬时频率的集合，作为LFM 信号的主成分，这种方法称为能量型主成分提取法.图10 SNR=10 dB 时能量型提取到的LFM 信号主成分图11 SNR=10 dB 时利用能量型提取到的主成分拟合的直线图10是SNR=10dB 时，使用能量型主成分提取得到LFM 的主成分的Hilbert 谱. 图11是拟合得到的直线，此时效果非常不错，0=2526399=480699a f 、ˆ..ˆ，与真实值的误差非常小，证明能量型主成分提取法具备抗噪性能，性能的分析将在下一节详述.3 实验与分析仍使用之前的仿真信号模型及参数，并加入不同SNR 的加性高斯白噪声，对文中所提方法的估计效果进行实验并分析. 为了评价2种方法对调频斜率a 和起始频率0f 的估计效果，定义估计偏差δ，在某SNR 条件下为()11ˆ=100%.N i i N δααα=⎛⎞−×⎜⎟⎝⎠∑ （14）式中：α是参数的设定值，ˆα是参数的估计值，N 为在某SNR 条件下对某一参数的重复估计次数，在这里设定100N =. 另外，文中所有的直线拟合均采用的是最小二乘直线拟合，利用该准则拟合得到的直线在坐标系上与数据点之间的距离平方和最小. 理想的高斯白噪声是覆盖全部频谱的，其功率谱密度是常数，而且与信号是不相关的；因此，最小二乘准则使得主分量中离散的数据点聚集起来，在很大程度上降低了高斯白噪声的影响.图12 斜率的估计偏差图13 起始频率的估计偏差图12、13是2种方法对2个参数在不同SNR 下的估计效果，步长为1 dB. 通过求取直线的斜率与截·50·应用科技第39卷距来完成LFM信号对于起始频率0f和调频斜率a的估计，部分SNR下的估计结果列于表1、2.表1 调频斜率的估计（标准250 Hz/s，100次的平均值）方法 5 dB 10 dB 15 dB 20 dB直接型-76.337 0 -56.618 214.228 0 111.373 7能量型211.425 0 250.862 3253.069 0 252.169 3表2 初始频率的估计（标准50 Hz/s，100次的平均值）方法 5 dB 10 dB 15 dB 20 dB直接型323.372 1 281.573 9217.412 8 145.552 6能量型87.268 5 49.463 047.420 3 48.310 9可以看出，在5~20 dB，虽然直接型提取法的起始频率估计误差曲线随着SNR的升高，趋势是下降的，但是无论对f还是a的估计误差都相对过大. 相比之下，能量型提取法则效果好得多，2条误差曲线均是随SNR的升高而下降，且均在SNR为10 dB位置附近估计性能达到最佳. 当SNR>10dB时，虽然误差有所上升，但依然是很小的.同时可以看出，随着SNR的降低，直接型主成分提取法效果不断恶化. 实验发现，在SNR较低的时候，如0 dB，由于所获得的主成分的点过于分散，甚至于人眼都已经很难判断LFM信号是否存在；此外，该方法提取出的主成分随着SNR的降低，其主成分的意义也在下降，表现是LFM在某些时刻的理论频率值出现在了2阶甚至3阶的IMF中，换句话说，集中在1阶IMF中的主成分在随着SNR的下降而不断向高阶IMF扩散，这也是提出能量型主成分提取法的主要依据. 能量型主成分提取法的效果也在随着SNR的降低而降低，但是其所获得的主成分的点相对来说是集中的，且拟合得到的直线更接近纯净的LFM信号的理论时频特性.4结束语介绍了HHT变换，并研究了LFM信号的HHT 分析，通过对文献[7]方法的研究和分析，得出通过EMD得到LFM信号主成分的这样一种本质思想，进而提出了一种LFM信号的主成分提取方法并进行信号的分析，且具备一定的参数估计性能. 在相同SNR 条件下：SNR非常高时，能量型提取法和直接型提取法的估计效果相差不多；随着SNR的降低，能量型提取法的估计效果依然较好，而直接型提取法却逐渐失效. 此外，加性高斯白噪声的随机性，对能量型主成分提取法的影响使得本方法缺乏一定的稳定性.综上所述，直接型主成分提取法基本只适用于对纯净LFM的分析；能量型主成分提取法却具备较好的估计性能，适用于较宽的信噪比范围，具有一定的实用价值.参考文献:[1] TRETTER S. Estimating the frequency of a noisy sinusoid bylinear regression[J]. IEEE Transactions on Information Theory, 1985, 31(6): 832-835.[2] HUANG N E, SHEN Z, LONG S R, et al. The empirical modedecomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of Royal Society of London, 1998(A454): 903-995.[3] KHALDI K, TURKI-HADJ A M, BOUDRAA A O. A newEMD denoising approach dedicated to voiced speech signals[C]// International Conference on Signals, Circuits andSystems. Monastir, Tunisia, 2008: 1-5.[4] LU Z M, LIU B S, SHEN L R. Speech endpoint detection instrong noisy environment based on the Hilbert-Huang transform[C]// International Conference on Mechatronics andAutomation. Changchun, China, 2009: 4322-4326.[4] CHEN S H, ZHANG R H, SU H B, et al. SAR andmultispectral image fusion using generalized IHS transform based on a trous wavelet and EMD decompositions[J].Sensors Journal, 2010, 10(3): 737-745.[6] ESPINOSA A G, ROSERO J A, CUSIDO J, et al. Faultdetection by means of Hilbert-Huang transform of the statorcurrent in a PMSM with demagnetization[J]. IEEE Transactions on Energy Conversion, 2010, 25(2): 312-318. [7] 崔华. 一种新的线性调频信号的瞬时频率估计方法[J]. 计算机应用研究, 2008, 25(8): 2532-2533, 2536.[8] FLANDRIN P, RILLING G. GONCALVES P. Empirical modedecomposition as a filter bank[J]. IEEE Signal Processing Letters, 2004, 11(2): 112-114.[9] COHEN L. Time-frequency analysis[M]. Englewood:Prentice-Hall, 1995: 39-45.。

一种估计线性调频信号参数的方法

一种估计线性调频信号参数的方法
彭榕辉;周荫清
【期刊名称】《遥测遥控》
【年(卷),期】2002(023)002
【摘要】合成孔径雷达实际发射的脉冲信号并非标准的线性调频信号,而是带有高次相位误差的调频信号,在高次相位误差较大的情况下,将对成像质量产生影响,因此,必须对发射信号的线性调频项和高次相位误差进行测试.文中提出利用一维相位解缠和多项式拟合测试发射脉冲信号参数的方法,仿真结果表明,具有较高的估计精度和较好的实用性.
【总页数】4页(P12-15)
【作者】彭榕辉;周荫清
【作者单位】北京航空航天大学电子工程系,北京,100083;北京航空航天大学电子工程系,北京,100083
【正文语种】中文
【中图分类】TN91
【相关文献】
1.一种多径线性调频信号的参数估计方法 [J], 史建锋;张锦政;马强
2.一种低复杂度线性调频信号参数估计算法 [J], 熊竹林;刘策伦;安建平
3.一种新型线性调频信号参数快速估计算法 [J], 姜晓斐;李弋鹏;贾燕婕
4.一种新的基于分数阶Fourier变换的线性调频信号的参数估计方法 [J], 陈蓉;汪
一鸣
5.一种基于新误差标准的ANF线性调频信号参数估计方法 [J], 赵红梅;崔艳因版权原因，仅展示原文概要，查看原文内容请购买。

一种线性调频信号参数提取方法

一种线性调频信号参数提取方法李彦;姜军;朱丰;柏又青【摘要】Aimed at the parameter extraction of Linear Frequency Modulation (LFM)signal,an novel pa-rameter extraction method is proposed in this paper based on Erosion transform and Hough transform. First,LFM signal is described by using B-distribution analysis method.Onthe basis of this,the Erosion transform is utilized to highlight the characters of radar LFM signal,and then Hough transform is em-ployed to extract LFM signal parameters.The simulation indicates that the method can not only extract the parameters of LFM signal in precision,but also has a good anti-noise performance.%针对线性调频(LFM)信号参数提取问题，提出了一种基于Erosion变换和 Hough 变换的LFM信号参数提取方法。

利用B-分布分析方法描述LFM信号，在此基础上，进一步利用Ero-sion变换突出雷达LFM信号的特点，利用Hough变换提取LFM信号参数。

仿真分析表明：该方法能够较精确的提取出LFM信号参数，误差不超过1%，并具有良好的抗噪性能。

【期刊名称】《空军工程大学学报（自然科学版）》【年(卷),期】2014(000)006【总页数】4页(P49-52)【关键词】LFM信号;参数提取;B-分布;Erosion变换;Hough变换【作者】李彦;姜军;朱丰;柏又青【作者单位】空军工程大学理学院，陕西西安，710051;空军工程大学理学院，陕西西安，710051;空军工程大学信息与导航学院，陕西西安，710077;空军工程大学理学院，陕西西安，710051【正文语种】中文【中图分类】TN953线性调频(Linear Frequency Modulation, LFM)信号是一种常见的雷达信号形式，在雷达成像和目标识别中具有广泛的应用[1]。