[go: up one dir, main page]

CN108362939B - 一种线性调频信号的频域参数测量方法 - Google Patents

一种线性调频信号的频域参数测量方法 Download PDF

Info

Publication number
CN108362939B
CN108362939B CN201810098772.1A CN201810098772A CN108362939B CN 108362939 B CN108362939 B CN 108362939B CN 201810098772 A CN201810098772 A CN 201810098772A CN 108362939 B CN108362939 B CN 108362939B
Authority
CN
China
Prior art keywords
frequency
signal
spectrum
frequency spectrum
time domain
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
Application number
CN201810098772.1A
Other languages
English (en)
Other versions
CN108362939A (zh
Inventor
杨振宇
马爽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Tiger Microwave Technology Co Ltd
Original Assignee
Chengdu Tiger Microwave Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Tiger Microwave Technology Co Ltd filed Critical Chengdu Tiger Microwave Technology Co Ltd
Priority to CN201810098772.1A priority Critical patent/CN108362939B/zh
Publication of CN108362939A publication Critical patent/CN108362939A/zh
Application granted granted Critical
Publication of CN108362939B publication Critical patent/CN108362939B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种线性调频信号的频域参数测量方法,包括以下步骤:S1.时域数据获取:对信号进行时域划窗,将信号划分为多个数据段;S2.数据变换:对各数据段进行离散傅里叶变换,得到各数据段对应的子带频谱;S3.频谱内插:对各个子带频谱进行内插处理,得到各数据段对应的点数扩展后的频谱信号;S4.频谱时延:在点数扩展后的频谱信号上增加时域延时信息;S5.频谱综合:将带有时域延时信息的所有频谱信号进行频谱综合,恢复原始信号的完整频谱;S6.频域参数测量:利用恢复得到的完整频谱,进行频域参数测量。本发明利用短时频谱综合恢复出原始信号频谱,拟合出了原始信号更多的频谱信息,提高频率分辨率,进而提高带宽、起始频率的测量精度。

Description

一种线性调频信号的频域参数测量方法
技术领域
本发明涉及频域参数测量,具体是涉及一种线性调频信号的频域参数测量方法。
背景技术
在对空域信号侦测的过程中,为及时得到截获信号脉冲描述字中的频域参数,如带宽、起始频率或中心频率,短时傅里叶变换(STFT)是一种常用的处理方法;通过对信号在时域滑窗,将信号划分为时间较短的数据段,再对各短时信号进行DFT(离散傅里叶变换)得到子带频谱,分析其对应的频域信息,但是对于调制信号例如线性调频信号(LFM),各短时数据段的频谱并不能准确反映原始完整信号的频域信息,图1为对窄数据窗口内信号进行频谱分析的示意图(短时信号频谱),图2为原始信号频谱,对比两图可以看出,由各段频谱仅反映原始信号的部分频谱信息,且对各子带频谱参数的测量受频率分辨率等的影响,频率测量误差较大,由各子带测量参数进一步获取原始信号的带宽、中心频率等参数必定会带来较大误差。
发明内容
本发明的目的在于克服现有技术的不足,提供一种线性调频信号的频域参数测量方法,利用短时频谱综合恢复出原始信号频谱,拟合出了原始信号更多的频谱信息,提高频率分辨率,进而提高带宽、起始频率的测量精度。
本发明的目的是通过以下技术方案来实现的:一种线性调频信号的频域参数测量方法,包括以下步骤:
S1.时域数据获取:对信号进行时域划窗,将信号划分为多个时域的数据段x(n),n=0,1,…,N-1,每个数据段包含N个离散数据;
S2.数据变换:对每一个数据段进行离散傅里叶变换,得到各个数据段对应的子带频谱X(k):
Figure BDA0001565774480000011
S3.频谱内插:对各个数据段对应的子带频谱X(k)进行内插处理,得到各数据段对应的点数扩展后的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号,r为正整数;
S4.频谱时延:利用傅里叶变换的时域延时性质,在每一个数据段对应的频谱信号X′(k)上增加时域延时信息;
S5.频谱综合:将带有时域延时信息的所有频谱信号X′(k)进行频谱综合,恢复原始信号的完整频谱;
S6.频域参数测量:利用恢复得到的完整频谱,对原始信号进行频域参数测量。
其中,所述步骤S3包括以下子步骤:
S301.设x(n)中补(r-1)*N个零,r为正整数,得到r*N点的新序列x′(n),则通过离散傅里叶变换,推导出x′(n)对应的频谱信号X′(k):
Figure BDA0001565774480000021
可知,X′(k)为X(k)的线性加权;
S302.根据推导结果,在频域对X(k)进行内插,以得到所需的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号。
其中,所述步骤S4包括以下子步骤:对于各个数据段,利用傅里叶变换的时域延时性质,将对应的频谱信号X′(k)上增加时域延时信息,得到:
Figure BDA0001565774480000022
DFT{}表示离散傅里叶变换,t表示时域窗口选取的第t个数据段,q表示时域窗口选取数据段的时间间隔,(t-1)q即为各数据段的延时信息;式中,
Figure BDA0001565774480000023
Figure BDA0001565774480000024
即为增加延时信息后的频谱信号。
本发明的有益效果是:本发明利用短时频谱综合恢复出原始信号频谱,拟合出了原始信号更多的频谱信息,提高频率分辨率,进而提高带宽、起始频率的测量精度。
附图说明
图1为窄数据窗口内信号频谱分析的示意图;
图2为原始信号频谱示意图;
图3为本发明的方法流程图;
图4为本发明实施例中的频谱恢复示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图3所示,一种线性调频信号的频域参数测量方法,包括以下步骤:
S1.时域数据获取:对信号进行时域划窗,将信号划分为多个时域的数据段x(n),n=0,1,…,N-1,每个数据段包含N个离散数据;
S2.数据变换:对每一个数据段进行离散傅里叶变换,得到各个数据段对应的子带频谱X(k):
Figure BDA0001565774480000031
S3.频谱内插:对各个数据段对应的子带频谱X(k)进行内插处理,得到各数据段对应的点数扩展后的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号,r为正整数;
S4.频谱时延:利用傅里叶变换的时域延时性质,在每一个数据段对应的频谱信号X′(k)上增加时域延时信息;
S5.频谱综合:将带有时域延时信息的所有频谱信号X′(k)进行频谱综合,拟合恢复原始信号的完整频谱;
S6.频域参数测量:利用恢复得到的完整频谱,对原始信号进行频域参数测量。
其中,所述步骤S3包括以下子步骤:
S301.设x(n)中补(r-1)*N个零,r为正整数,得到r*N点的新序列x′(n),则通过离散傅里叶变换,推导出x′(n)对应的频谱信号X′(k):
Figure BDA0001565774480000032
可知,在k=m时,X′(k)为X(k)的线性加权,可进一步简化为
Figure BDA0001565774480000033
在k≠m时,X′(k)依然为X(k)的线性加权,只是权值较为复杂而已。
S302.根据推导结果,在频域对X(k)进行内插,以得到所需的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号。
所述步骤S4包括以下子步骤:
对于各个数据段,利用傅里叶变换的时域延时性质,将对应的频谱信号X′(k)上增加时域延时信息,得到:
Figure BDA0001565774480000041
DFT{}表示离散傅里叶变换,t表示时域窗口选取的第t个数据段,q表示时域窗口选取数据段的时间间隔,(t-1)q即为各数据段的延时信息;式中,
Figure BDA0001565774480000042
Figure BDA0001565774480000043
即为增加延时信息后的频谱信号。
在本申请的实施例中,利用本申请的方法为一个带宽为40MHz的LFM信号综合恢复出的完整信号频谱,得到的结果如图4所示,对比图2可知,恢复出的频谱与原始频谱近似一致,由此可以精确测量出带宽、起始频率参数。
综上,本发明利用短时频谱综合恢复出原始信号频谱,拟合出了原始信号更多的频谱信息,提高频率分辨率,进而提高带宽、起始频率的测量精度。

Claims (2)

1.一种线性调频信号的频域参数测量方法,其特征在于:包括以下步骤:
S1.时域数据获取:对信号进行时域划窗,将信号划分为多个时域的数据段x(n),n=0,1,…,N-1,每个数据段包含N个离散数据;
S2.数据变换:对每一个数据段进行离散傅里叶变换,得到各个数据段对应的子带频谱X(k):
Figure FDA0002459983290000011
S3.频谱内插:对各个数据段对应的子带频谱X(k)进行内插处理,得到各数据段对应的点数扩展后的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号,r为正整数;
所述步骤S3包括以下子步骤:
S301.设x(n)中补(r-1)*N个零,r为正整数,得到r*N点的新序列x′(n),则通过离散傅里叶变换,推导出x′(n)对应的频谱信号X′(k):
Figure FDA0002459983290000012
可知,X′(k)为X(k)的线性加权;
S302.根据推导结果,在频域对X(k)进行内插,以得到所需的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号;
S4.频谱时延:利用傅里叶变换的时域延时性质,在每一个数据段对应的频谱信号X′(k)上增加时域延时信息;
S5.频谱综合:将带有时域延时信息的所有频谱信号X′(k)进行频谱综合,恢复原始信号的完整频谱;
S6.频域参数测量:利用恢复得到的完整频谱,对原始信号进行频域参数测量。
2.根据权利要求1所述的一种线性调频信号的频域参数测量方法,其特征在于:所述步骤S4包括以下子步骤:
对于各个数据段,利用傅里叶变换的时域延时性质,将对应的频谱信号X′(k)上增加时域延时信息,得到:
Figure FDA0002459983290000021
DFT{}表示离散傅里叶变换,t表示时域窗口选取的第t个数据段,q表示时域窗口选取数据段的时间间隔,(t-1)q即为各数据段的延时信息;式中,
Figure FDA0002459983290000022
Figure FDA0002459983290000023
即为增加延时信息后的频谱信号。
CN201810098772.1A 2018-01-31 2018-01-31 一种线性调频信号的频域参数测量方法 Active CN108362939B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810098772.1A CN108362939B (zh) 2018-01-31 2018-01-31 一种线性调频信号的频域参数测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810098772.1A CN108362939B (zh) 2018-01-31 2018-01-31 一种线性调频信号的频域参数测量方法

Publications (2)

Publication Number Publication Date
CN108362939A CN108362939A (zh) 2018-08-03
CN108362939B true CN108362939B (zh) 2020-06-23

Family

ID=63007713

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810098772.1A Active CN108362939B (zh) 2018-01-31 2018-01-31 一种线性调频信号的频域参数测量方法

Country Status (1)

Country Link
CN (1) CN108362939B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110007140A (zh) * 2019-04-16 2019-07-12 北京遥感设备研究所 一种雷达设备lfm信号调频斜率测试方法及系统
CN117118536B (zh) * 2023-10-25 2023-12-19 南京派格测控科技有限公司 调频稳定性的确定方法、装置、设备及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003012416A1 (en) * 2001-08-01 2003-02-13 Liposcience, Inc. Method of determining presence and concentration of lipoprotein x in blood plasma and serum
CN101566684A (zh) * 2009-06-08 2009-10-28 武汉大学 雷达数字信号处理方法及装置
CN101776747A (zh) * 2010-01-15 2010-07-14 山东大学 基于自适应重采样的雷达信号压缩方法
CN101833035A (zh) * 2010-04-19 2010-09-15 天津大学 线性调频信号参数估计方法及其实施装置
CN102055411A (zh) * 2010-12-22 2011-05-11 成都凯腾四方数字广播电视设备有限公司 基于多通道反馈的功率放大器线性化校正电路及方法
CN102308546A (zh) * 2011-07-26 2012-01-04 华为技术有限公司 一种多载波光信号的接收方法和装置
CN102680948A (zh) * 2012-05-15 2012-09-19 东南大学 一种线性调频信号调频率和起始频率估计方法
CN103226190A (zh) * 2013-03-19 2013-07-31 南京大学 一种基于真延时的宽带相控阵波束形成方法
CN103532656A (zh) * 2013-08-08 2014-01-22 北京理工大学 基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法
CN103869162A (zh) * 2014-03-05 2014-06-18 湖南大学 一种基于时域准同步的动态信号相量测量方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003012416A1 (en) * 2001-08-01 2003-02-13 Liposcience, Inc. Method of determining presence and concentration of lipoprotein x in blood plasma and serum
CN101566684A (zh) * 2009-06-08 2009-10-28 武汉大学 雷达数字信号处理方法及装置
CN101776747A (zh) * 2010-01-15 2010-07-14 山东大学 基于自适应重采样的雷达信号压缩方法
CN101833035A (zh) * 2010-04-19 2010-09-15 天津大学 线性调频信号参数估计方法及其实施装置
CN102055411A (zh) * 2010-12-22 2011-05-11 成都凯腾四方数字广播电视设备有限公司 基于多通道反馈的功率放大器线性化校正电路及方法
CN102308546A (zh) * 2011-07-26 2012-01-04 华为技术有限公司 一种多载波光信号的接收方法和装置
CN102680948A (zh) * 2012-05-15 2012-09-19 东南大学 一种线性调频信号调频率和起始频率估计方法
CN103226190A (zh) * 2013-03-19 2013-07-31 南京大学 一种基于真延时的宽带相控阵波束形成方法
CN103532656A (zh) * 2013-08-08 2014-01-22 北京理工大学 基于分数阶傅里叶域信道化的宽带线性调频信号多假目标干扰方法
CN103869162A (zh) * 2014-03-05 2014-06-18 湖南大学 一种基于时域准同步的动态信号相量测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LFM signal detection using the origin moment of fractional spectrum;Yin Zhiping;《2008 9th International Conference on Signal Processing》;20081029;全文 *
Parameters estimation of LFM signal based on fractional order cross spectrum;Kaihui Ding;《Proceedings of 2011 International Conference on Electronic & Mechanical Engineering and Information Technology》;20110814;全文 *
连续波背景下多信号处理方法;黄慧慧;《数据采集与处理》;20060930;第21卷(第3期);全文 *

Also Published As

Publication number Publication date
CN108362939A (zh) 2018-08-03

Similar Documents

Publication Publication Date Title
Oberlin et al. The Fourier-based synchrosqueezing transform
CN108362939B (zh) 一种线性调频信号的频域参数测量方法
CN119622236B (zh) 一种基于峰度最大化与连续性准则的自适应时频分析方法
CN106772270A (zh) 一种雷达回波信号的采样方法及重构方法
CN111077371B (zh) 一种提高相位测量精度的方法和装置
CN103487669A (zh) 基于任意频率信号间相位特征处理的相位噪声测量方法
CN107966687B (zh) 基于部分自相关谱的mimo雷达信号调制类型识别方法
CN104422956A (zh) 一种基于稀疏脉冲反演的高精度地震谱分解方法
CN117130048B (zh) 一种时频谱的时间同步压缩提取变换方法
Sekhar et al. Adaptive spectrogram vs. adaptive pseudo-Wigner–Ville distribution for instantaneous frequency estimation
CN113126055A (zh) Mimo雷达信号处理方法、装置、终端设备及存储介质
CN117373484B (zh) 一种基于特征变换的开关柜声纹故障检测方法
Pal et al. Soft-thresholding for spectrum sensing with coprime samplers
CN117706213A (zh) 一种电磁继电器抗电磁脉冲干扰测试方法、装置和设备
WO2023226572A1 (zh) 特征表示的提取方法、装置、设备、介质及程序产品
Qi et al. On timing skews of multicoset samplers in compressive spectrum sensing for millimeter-wave
CN119646393B (zh) 非平稳信号处理方法、装置、电子设备及存储介质
CN110687362A (zh) 一种实时相噪测试方法
RU2813684C1 (ru) Способ и устройство измерения спектра и кепстральных параметров информационных акустических сигналов телерадиовещания
Wang et al. Fast algorithm on parameter estimation of wideband LFM signal based on down-chirp and CS
Saulig et al. Nonstationary signals information content estimation based on the local Rényi entropy in the time-frequency domain
CN117033908A (zh) 一种基于多尺度小波的实时时频分析方法
Liang et al. A Structure of MWC Based on Overlapping Window Method for Wideband LFM
CN110888916B (zh) 一种多通道间的数据处理方法、装置、设备和存储介质
Yuehua et al. The method of range profile for step frequency MMW radar based on wavelet transform power spectrum estimator

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