CN108362939B - 一种线性调频信号的频域参数测量方法 - Google Patents
一种线性调频信号的频域参数测量方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 86
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 8
- 238000013501 data transformation Methods 0.000 claims abstract description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements 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):
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):
可知,X′(k)为X(k)的线性加权;
S302.根据推导结果,在频域对X(k)进行内插,以得到所需的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号。
其中,所述步骤S4包括以下子步骤:对于各个数据段,利用傅里叶变换的时域延时性质,将对应的频谱信号X′(k)上增加时域延时信息,得到:
DFT{}表示离散傅里叶变换,t表示时域窗口选取的第t个数据段,q表示时域窗口选取数据段的时间间隔,(t-1)q即为各数据段的延时信息;式中,
本发明的有益效果是:本发明利用短时频谱综合恢复出原始信号频谱,拟合出了原始信号更多的频谱信息,提高频率分辨率,进而提高带宽、起始频率的测量精度。
附图说明
图1为窄数据窗口内信号频谱分析的示意图;
图2为原始信号频谱示意图;
图3为本发明的方法流程图;
图4为本发明实施例中的频谱恢复示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图3所示,一种线性调频信号的频域参数测量方法,包括以下步骤:
S1.时域数据获取:对信号进行时域划窗,将信号划分为多个时域的数据段x(n),n=0,1,…,N-1,每个数据段包含N个离散数据;
S2.数据变换:对每一个数据段进行离散傅里叶变换,得到各个数据段对应的子带频谱X(k):
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):
在k≠m时,X′(k)依然为X(k)的线性加权,只是权值较为复杂而已。
S302.根据推导结果,在频域对X(k)进行内插,以得到所需的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号。
所述步骤S4包括以下子步骤:
对于各个数据段,利用傅里叶变换的时域延时性质,将对应的频谱信号X′(k)上增加时域延时信息,得到:
DFT{}表示离散傅里叶变换,t表示时域窗口选取的第t个数据段,q表示时域窗口选取数据段的时间间隔,(t-1)q即为各数据段的延时信息;式中,
在本申请的实施例中,利用本申请的方法为一个带宽为40MHz的LFM信号综合恢复出的完整信号频谱,得到的结果如图4所示,对比图2可知,恢复出的频谱与原始频谱近似一致,由此可以精确测量出带宽、起始频率参数。
综上,本发明利用短时频谱综合恢复出原始信号频谱,拟合出了原始信号更多的频谱信息,提高频率分辨率,进而提高带宽、起始频率的测量精度。
Claims (2)
1.一种线性调频信号的频域参数测量方法,其特征在于:包括以下步骤:
S1.时域数据获取:对信号进行时域划窗,将信号划分为多个时域的数据段x(n),n=0,1,…,N-1,每个数据段包含N个离散数据;
S2.数据变换:对每一个数据段进行离散傅里叶变换,得到各个数据段对应的子带频谱X(k):
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):
可知,X′(k)为X(k)的线性加权;
S302.根据推导结果,在频域对X(k)进行内插,以得到所需的频谱信号X′(k),所述频谱信号X′(k)为包含r*N个数据的频域信号;
S4.频谱时延:利用傅里叶变换的时域延时性质,在每一个数据段对应的频谱信号X′(k)上增加时域延时信息;
S5.频谱综合:将带有时域延时信息的所有频谱信号X′(k)进行频谱综合,恢复原始信号的完整频谱;
S6.频域参数测量:利用恢复得到的完整频谱,对原始信号进行频域参数测量。
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)
| 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)
| 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 | 湖南大学 | 一种基于时域准同步的动态信号相量测量方法 |
-
2018
- 2018-01-31 CN CN201810098772.1A patent/CN108362939B/zh active Active
Patent Citations (10)
| 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)
| 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 |















