CN102707157B - 一种基于功率谱的单频脉冲信号参数估值方法 - Google Patents
一种基于功率谱的单频脉冲信号参数估值方法 Download PDFInfo
- Publication number
- CN102707157B CN102707157B CN201210155036.8A CN201210155036A CN102707157B CN 102707157 B CN102707157 B CN 102707157B CN 201210155036 A CN201210155036 A CN 201210155036A CN 102707157 B CN102707157 B CN 102707157B
- Authority
- CN
- China
- Prior art keywords
- signal
- pure
- spectrum
- power
- tau
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于功率谱的单频脉冲信号参数估值方法,该方法包括:提取待估值的单频脉冲信号功率谱,将所述功率谱的频率轴映射到时间轴,获得功率时间信号;将所述功率时间信号进行傅里叶变换,分别得到所述功率时间信号的幅度谱和相位谱;根据所述单频脉冲信号的脉冲宽度与所述幅度谱之间的关系,计算所述脉冲宽度;根据所述单频脉冲信号的载频与所述相位谱之间的关系,计算所述载频,实现单频脉冲信号的估值。本发明解决了现有技术中不能准确地估计出低信噪比环境下的单频脉冲信号参数问题,适用于单频脉冲信号功率谱由于噪声影响造成的谱峰扭曲、斜变以及主瓣过零点被淹没的情况。
Description
技术领域
本发明涉及信号处理技术,特别涉及一种基于功率谱的单频脉冲信号参数估值方法。
背景技术
信号参数的估值历来是通信领域研究和应用的热点,对于军用电子对抗以及民用频谱监管都具有极其重要的意义。估值是指对所测出的信号分析,获取信号的有用参数,以便确定信号用途和恢复信号。对于单频脉冲信号参数估值算法,目前相对简单的参数估值方案是调制参数估值的趋势和主流,但是算法的抗噪性能差。
单频脉冲信号是雷达及声纳中常用的信号形式,是指在发送的每个脉冲内载频是固定不变,具有方便、快捷、易于实现的优点。可以表示为:
其中:A是信号强度,τ是脉冲宽度,fc是载波频率,是信号初始相位。
信号功率谱表示为:
P(f)=A2τ2sinc2[τ(f-fc)]
常用的单频脉冲信号参数有载频、脉冲宽度。信号载频能表征信号的辐射频段,可以反映雷达的功能和用途,对提高声纳跟踪测量***的定位精度与改善跟踪效果具有重要意义。它一般不随时间地点变化而变化,是一个较为固定的分选参数。脉冲宽度是脉冲结束时间和到达时间之间的时间宽度,它是描述脉冲特征的一个重要参数。
目前常用的单频脉冲信号估值方法包括瞬时自相关法和过零点检测法,但都存在抗噪性能差的缺点,不适合在复杂环境下进行估计。
发明内容
本发明目的在于:克服已有单频脉冲信号估值技术存在的不能准确地估计出低信噪比环境下信号参数的问题,适用于单频脉冲信号功率谱由于噪声影响造成的谱峰扭曲、斜变以及主瓣过零点被淹没的情况。
本发明的技术方案是:本发明提供了一种基于功率谱的单频脉冲信号参数估值方法,其特征在于,该方法包括:
步骤1、提取待估值的单频脉冲信号功率谱P(f),将所述功率谱的频率轴映射到时间轴,获到功率谱的强度与时间的关系,称为功率时间信号P(t),将所述功率时间信号进行傅里叶变换,分别得到所述功率时间信号的幅度谱|F′(x)|和相位谱
其中:将所述功率谱P(f)=A2τ2sinc2[τ(f-fc)]的频率轴按照f=t的原则映射到时间轴,得到所述功率时间信号P(t),A为所述单频脉冲信号强度,表达式如下所示:
P(t)=A2τ2sinc2[τ(t-fc)]
对所述功率时间信号P(t)进行傅里叶变换得到F(x),x为变换后的横轴分量, F(x)表达式如下表示:
由于傅里叶变换的频移特性F(x)关于纵轴对称,取x≥0的部分F′(x)进行分析,表达式如下所示:
得到所述功率时间信号的幅度谱|F′(x)|,表达式如下所示:
得到所述功率时间信号的相位谱表达式如下所示:
步骤2、根据所述单频脉冲信号的脉冲宽度τ与所述功率时间信号的幅度谱|F′(x)|之间的关系,计算所述单频脉冲信号的脉冲宽度τ;
其中:根据所述功率时间信号幅度谱的表达式,计算所述功率时间信号幅度谱的首个过零点x0,所述首个过零点x0即为所述单频脉冲信号的脉冲宽度τ=x0;
步骤3、根据所述单频脉冲信号的载频fc与所述功率时间信号的相位谱之间的关系,计算所述单频脉冲信号的载频fc;
其中:根据所述功率时间信号相位谱的表达式,计算所述功率时间信号的相位谱中第一个极大值的横坐标x1和最后一个极大值的横坐标xN,N为极大值个数,所述两个极大值的横坐标x1和xN的差为Δx=xN-x1,所述单频脉冲信号的载频fc的表达式如下所示:
本发明的有益效果是:本发明通过提取待估值的单频脉冲信号功率谱,将所述功率谱的频率轴映射到时间轴,获得功率时间信号;将所述功率时间信号进行傅里叶变换,分别得到所述功率时间信号的幅度谱和相位谱;根据所述单频脉冲信号的脉冲宽度与所述幅度谱之间的关系,计算所述脉冲宽度;根据所述单频脉冲信号的载频与所述相位谱之间的关系,计算所述载频,实现单频脉冲信号的估值。本发明解决了现有技术中不能准确地估计出低信噪比环境下的单频脉冲信号参数问题,适用于单频脉冲信号功率谱由于噪声影响造成的谱峰扭曲、斜变以及主瓣过零点被淹没的情况。
附图说明
图1为本发明提供的基于功率谱的单频脉冲信号参数估值的流程图;
图2为本发明实施例功率时间信号的示意图;
图3为本发明实施例功率时间信号的幅度谱;
图4为本发明实施例功率时间信号的相位谱;
图5为本发明实施例信噪比-10dB到10dB的单频脉冲信号脉冲宽度估值误差曲线;
图6为本发明实施例信噪比-10dB到10dB的单频脉冲信号载频估值误差曲线。
具体实施方式
以下将参照图1-6对本发明的具体实施方式进行说明。
如图1所示,本发明实施例进行基于功率谱的单频脉冲信号参数估值方法包括下列步骤:
步骤1、提取待估值的单频脉冲信号功率谱P(f),将所述功率谱的频率轴映射到时间轴,得到功率谱的强度与时间的关系,称为功率时间信号P(t);将所述功率时间信号进行傅里叶变换,分别得到所述功率时间信号的幅度谱|F′(x)|和相位谱
其中,步骤1中所述功率时间信号的幅度谱和相位谱的确定方法具体包括:
将所述功率谱P(f)=A2τ2sinc2[τ(f-fc)]的频率轴按照f=t的原则映射到时间轴,得到所述功率时间信号P(t),A为所述单频脉冲信号强度,表达式如下所示:
P(t)=A2τ2sinc2[τ(t-fc)]
对所述功率时间信号P(t)进行傅里叶变换得到F(x),x为变换后的横轴分量, F(x)表达式如下表示:
由于傅里叶变换的频移特性F(x)关于纵轴对称,取x≥0的部分F′(x)进行分析,,表达式如下所示:
得到所述功率时间信号的幅度谱|F′(x)|,表达式如下所示:
得到所述功率时间信号的相位谱表达式如下所示:
步骤2、根据所述单频脉冲信号的脉冲宽度τ与所述功率时间信号的幅度谱|F′(x)|之间的关系,计算所述单频脉冲信号的脉冲宽度τ;
其中,步骤2中计算所述单频脉冲信号脉冲宽度的方法具体包括:
根据所述功率时间信号幅度谱的表达式,计算所述功率时间信号幅度谱的首个过零点x0,所述首个过零点x0即为所述单频脉冲信号的脉冲宽度τ=x0。
需要说明的是:
确定单频脉冲信号的脉冲宽度与所述功率时间信号的幅度谱之间关系的方法具体包括:
功率时间信号的幅度谱表达式如下所示:
根据矩形函数 可知,当x≥τ时|F′(x)|=0,由此得到幅度谱的首个过零点x0为单频脉冲信号的脉冲宽度τ=x0。
步骤3、根据所述单频脉冲信号的载频fc与所述功率时间信号的相位谱之间的关系,计算所述单频脉冲信号的载频fc。
其中,步骤3中计算所述单频脉冲信号载频的方法具体包括:
根据所述功率时间信号相位谱的表达式,计算所述功率时间信号的相位谱中第一个极大值的横坐标x1和最后一个极大值的横坐标xN,N为极大值个数,所述两个极大值的横坐标x1和xN的差为Δx=xN-x1,所述单频脉冲信号的载频fc的表达式如下所示:
需要说明的是:
确定单频脉冲信号的载频与所述功率时间信号的相位谱之间关系的方法具体包括:
功率时间信号的相位谱表达式如下所示:
根据相位谱的表达式可以得出相位谱的频率即为信号的载频,由此得到单频脉冲信号的载频fc的表达式如下所示:
本发明以单频脉冲信号为实施例对本发明方法进行详细说明。
本发明实施例可以应用到军用电子对抗以及民用频谱监管领域,在低信噪比环境下准确地估计出信号参数,可以反映雷达的功能和用途,对提高声纳跟踪测量***的定位精度与改善跟踪效果具有重要意义。
本发明以下内容用缩写CW表示单频脉冲信号,用缩写SNR表示信噪比。其中,载频fc=500MHz,脉冲宽度τ=1μs的CW信号得到的功率时间信号如图2所示,对功率时间信号进行傅里叶变换,分别得到功率时间信号的幅度谱如图3所示和相位谱如图4所示。
本发明在MATLAB6.5环境下估计CW信号参数。
CW信号载频fc=500MHz,脉冲宽度τ=1μs,SNR变化范围-10dB—10dB。本发明对CW信号的参数估值结果如表1、2所示;脉冲宽度估计误差曲线如图5所示,载频估计误差曲线如图6所示。由以上结果可知,当SNR在-10dB-10dB,脉冲宽度估值误差小于0.019,载频估值误差均为0,由此分析可知,在低信噪比情况下本发明方法仍能对CW信号进行参数估值,且精度较高。
表1信噪比变化下参数估值结果
表2信噪比变化下参数估值结果
Claims (1)
1.一种基于功率谱的单频脉冲信号参数估值方法,其特征在于,该方法包括:
步骤1、提取待估值的单频脉冲信号功率谱P(f),将所述功率谱的频率轴映射到时间轴,获得功率谱的强度与时间的关系,称为功率时间信号P(t),将所述功率时间信号进行傅里叶变换,分别得到所述功率时间信号的幅度谱|F′(x)|和相位谱
其中:将所述功率谱P(f)=A2τ2sinc2[τ(f-fc)]的频率轴按照f=t的原则映射到时间轴,得到所述功率时间信号P(t),A为所述单频脉冲信号强度,τ表示脉冲宽度,fc表示载波频率,表达式如下所示:
P(t)=A2τ2sinc2[τ(t-fc)]
对所述功率时间信号P(t)进行傅里叶变换得到F(x),x为变换后的横轴分量, F(x)表达式如下表示:
由于傅里叶变换的频移特性F(x)关于纵轴对称,取x≥0的部分F′(x)进行分析,表达式如下所示:
得到所述功率时间信号的幅度谱|F′(x)|,表达式如下所示:
得到所述功率时间信号的相位谱表达式如下所示:
步骤2、根据所述单频脉冲信号的脉冲宽度τ与所述功率时间信号的幅度谱|F′(x)|之间的关系,计算所述单频脉冲信号的脉冲宽度τ;
其中:根据所述功率时间信号幅度谱的表达式,计算所述功率时间信号幅度谱的首个过零点x0,所述首个过零点x0即为所述单频脉冲信号的脉冲宽度τ=x0;
步骤3、根据所述单频脉冲信号的载频fc与所述功率时间信号的相位谱之间的关系,计算所述单频脉冲信号的载频fc;
其中:根据所述功率时间信号相位谱的表达式,计算所述功率时间信号的相位谱中第一个极大值的横坐标x1和最后一个极大值的横坐标xN,N为极大值个数,所述两个极大值的横坐标x1和xN的差为Δx=xN-x1,所述单频脉冲信号的载频fc的表达式如下所示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210155036.8A CN102707157B (zh) | 2012-05-18 | 2012-05-18 | 一种基于功率谱的单频脉冲信号参数估值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210155036.8A CN102707157B (zh) | 2012-05-18 | 2012-05-18 | 一种基于功率谱的单频脉冲信号参数估值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102707157A CN102707157A (zh) | 2012-10-03 |
CN102707157B true CN102707157B (zh) | 2015-10-28 |
Family
ID=46900090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210155036.8A Expired - Fee Related CN102707157B (zh) | 2012-05-18 | 2012-05-18 | 一种基于功率谱的单频脉冲信号参数估值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102707157B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103063909B (zh) * | 2012-12-18 | 2015-10-28 | 天津理工大学 | 一种基于功率谱的线性调频信号参数估值方法 |
CN103441975B (zh) * | 2013-08-30 | 2016-12-28 | 天津理工大学 | 一种基于功率谱的二相编码信号参数估值方法 |
CN108469602B (zh) * | 2018-03-29 | 2021-09-07 | 东南大学 | 一种基于谱特征提取的脉冲信号类型自动判别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1395400A (zh) * | 2002-07-05 | 2003-02-05 | 清华大学 | 无需先验知识的常用数字信号调制方式的自动识别方法 |
CN101743731A (zh) * | 2007-11-05 | 2010-06-16 | 英国广播公司 | 用于信号发现的方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010037189A1 (en) * | 2000-01-20 | 2001-11-01 | Dan Onu | Method of estimating phase noise spectral density and jitter in a periodic signal |
US7778785B2 (en) * | 2008-02-14 | 2010-08-17 | Advantest Corporation | Signal-to-noise ratio measurement for discrete waveform |
-
2012
- 2012-05-18 CN CN201210155036.8A patent/CN102707157B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1395400A (zh) * | 2002-07-05 | 2003-02-05 | 清华大学 | 无需先验知识的常用数字信号调制方式的自动识别方法 |
CN101743731A (zh) * | 2007-11-05 | 2010-06-16 | 英国广播公司 | 用于信号发现的方法和装置 |
Non-Patent Citations (3)
Title |
---|
A new approach of pitch detection based on morphology filter and wavelet transform;zhang J,Wang H P;《Proc. of the Second International Workshop on Education Technology and Computer Science》;20101231;第751-753页 * |
基于能量分离算法的多分量线性调频信号频率估计;刘柳等;《天津理工大学学报》;20110430;第27卷(第2期);第10-12页 * |
非参数模型谱估计算法性能比较;张佳明等;《天津理工大学学报》;20090831;第25卷(第4期);第19-22页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102707157A (zh) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901429B (zh) | 基于稀疏孔径的机动目标逆合成孔径雷达成像方法 | |
CN107942322B (zh) | 多径利用雷达时间反演目标距离估计方法 | |
CN102323575B (zh) | Pd雷达在微弱信号检测过程中的距离走动校正方法 | |
CN103941089B (zh) | 基于dft的正弦信号频率估计方法 | |
CN103063909B (zh) | 一种基于功率谱的线性调频信号参数估值方法 | |
CN104678372B (zh) | 正交频分复用雷达超分辨距离与角度值联合估计方法 | |
CN107612587A (zh) | 一种用于跳频非合作通信中跳频信号的参数估计方法 | |
Zheng et al. | ISAR imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis | |
CN103364783B (zh) | 基于单通道sar的运动目标径向速度无模糊估计方法 | |
CN101702701B (zh) | 极低信噪比下频率偏移的估计与补偿方法 | |
CN103278807B (zh) | 双通道欠采样线扫频脉冲信号的时延估计方法 | |
CN103856431B (zh) | 基于快速傅里叶变换的mpsk调制下开环载波同步算法 | |
CN109412644B (zh) | 一种直扩msk信号的多普勒频率估计方法 | |
US20150285905A1 (en) | Sub-nyquist radar processing using doppler focusing | |
CN102841337A (zh) | 一种天波超视距雷达回波信号去非线性相位污染方法 | |
CN102707157B (zh) | 一种基于功率谱的单频脉冲信号参数估值方法 | |
CN104219761A (zh) | 一种基于最大斜率的超宽带无线定位方法 | |
CN102279396B (zh) | 一种基于分数阶Fourier变换的宽带线性调频脉冲测距方法 | |
CN106330342A (zh) | 一种低计算复杂度的水声通信多普勒因子估计方法 | |
CN106154257A (zh) | 基于FFT与apFFT的精密测量雷达二次测频方法 | |
CN106569182B (zh) | 基于最小熵的相位编码信号载频估计方法 | |
CN116643250B (zh) | 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法 | |
Wu et al. | Frequency estimation algorithm for ranging of millimeter wave LFMCW radar | |
CN104901718A (zh) | 基于直接序列扩频信号载波频率测量的多普勒估计方法 | |
CN104320360B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151028 Termination date: 20200518 |