CN109490853A - 一种线性调频脉冲信号中心频率处谱线值确定方法 - Google Patents

一种线性调频脉冲信号中心频率处谱线值确定方法 Download PDF

Info

Publication number
CN109490853A
CN109490853A CN201710814543.0A CN201710814543A CN109490853A CN 109490853 A CN109490853 A CN 109490853A CN 201710814543 A CN201710814543 A CN 201710814543A CN 109490853 A CN109490853 A CN 109490853A
Authority
CN
China
Prior art keywords
spectral line
value
centre frequency
frequency
data flow
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.)
Granted
Application number
CN201710814543.0A
Other languages
English (en)
Other versions
CN109490853B (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.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
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 Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201710814543.0A priority Critical patent/CN109490853B/zh
Publication of CN109490853A publication Critical patent/CN109490853A/zh
Application granted granted Critical
Publication of CN109490853B publication Critical patent/CN109490853B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种线性调频脉冲信号中心频率处谱线值确定方法，本发明方法首先构建谱线值确定***，所构建的***包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块。数据流生成模块对线性调频脉冲信号进行采样，生成数据流；频域点索引值生成模块获取中心频率处的频域点索引值；中心频率处谱线值获取模块使用数据流和频域点索引值计算中心频率处谱线值。本发明方法解决了传统的信号特定位置谱线值确定方法计算量大的问题。

Description

一种线性调频脉冲信号中心频率处谱线值确定方法
技术领域
本发明涉及一种信号谱线值确定方法，特别是一种线性调频脉冲信号中心频率处谱线值确定方法。
背景技术
在分析主被动雷达复合导引头产生的线性调频脉冲信号时，需要对线性调频脉冲信号中心频率处的谱线值进行确定。传统的信号特定位置谱线值确定方法为：对线性调频脉冲信号进行采样，形成采样序列；对采样序列进行快速傅立叶变换处理，将采样序列从时域变换到频域，得到频域数据；计算出特定位置处所对应的频域点索引值；在频域数据中找到频域点索引值所对应的谱线值。此方法易实现，但在进行快速傅立叶变换处理时，需要计算所有数据点的谱线值才能得到特定位置的谱线值，也就是说，计算全部数据点的谱线值和只计算特定位置这一个频点的谱线值，所需的运算量是一样的，很明显，使用FFT的传统信号特定位置谱线值确定方法浪费的运算量太多。
发明内容
本发明目的在于提供一种线性调频脉冲信号中心频率处谱线值确定方法，解决传统方法计算量大的问题。
一种线性调频脉冲信号中心频率处谱线值确定方法，其具体步骤为：
第一步搭建谱线值确定***
谱线值确定***，包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块。
数据流生成模块的功能为：对线性调频脉冲信号进行采样，生成数据流；
频域点索引值生成模块的功能为：获取中心频率处的频域点索引值；
中心频率处谱线值获取模块的功能为：使用数据流和频域点索引值计算中心频率处谱线值。
第二步数据流生成模块对线性调频脉冲信号进行采样，生成数据流
数据流生成模块在脉冲宽度的时间范围内对线性调频脉冲信号x(t)进行采样，生成数据流x(n)，采样频率为fs，t为时间参数，n为时域点索引，n＝1,2,…,N，N为数据流长度，N＝fs·τ，τ为线性调频脉冲信号的脉冲宽度。
第三步频域点索引值生成模块获取中心频率处的频域点索引值
频域点索引值生成模块获取中心频率处所对应的频域点索引值k0：k0＝f0·N/fs，其中f0为线性调频脉冲信号的中心频率。
第四步中心频率处谱线值获取模块通过数据流和频域点索引值获取中心频率处谱线值
中心频率处谱线值获取模块使用数据流x(n)和频域点索引值k0进行递推，获取递推运算结果y(n)：
递推的初始条件为y(-1)＝0。其中y(n-1)为递推运算结果y(n)上一个数据点的值；WN为旋转因子，WN＝e-j(2π/N)，j表示e-j(2π/N)为复数表示形式；为复数固定值；表示将与y(n-1)进行复数相乘。
中心频率处谱线值获取模块通过递推获取y(N)的值，y(N)的值就是中心频率处的谱线值。
更优的，其中第二步中采样频率fs≥4(f0+B/2)，其中f0为线性调频脉冲信号的中心频率，B为线性调频脉冲信号的调制带宽，f0+B/2为线性调频脉冲信号的最大频率。
一种线性调频脉冲信号中心频率处谱线值确定***，包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块。数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块顺次连接。数据流生成模块对线性调频脉冲信号进行采样，生成数据流；频域点索引值生成模块获取中心频率处的频域点索引值；中心频率处谱线值获取模块使用数据流和频域点索引值获取中心频率处谱线值。
本发明方法从窄带滤波器组的数学模型角度考虑，使用线性调频脉冲信号数据流和中心频率处的频域索引值，并通过递推运算最终确定线性调频脉冲信号中心频率处谱线值。本方法的运算效率高于使用FFT的传统信号特定位置谱线值确定方法。目前，本方法已成功应用在末制导的SAR成像导引头试验样机上，试验结果表明：确定的线性调频脉冲信号中心频率处谱线值准确，满足***使用要求。
具体实施方式
实施例1
一种线性调频脉冲信号中心频率处谱线值确定方法，其具体步骤为：
第一步搭建谱线值确定***
谱线值确定***，包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块。
数据流生成模块的功能为：对线性调频脉冲信号进行采样，生成数据流；
频域点索引值生成模块的功能为：获取中心频率处的频域点索引值；
中心频率处谱线值获取模块的功能为：使用数据流和频域点索引值计算中心频率处谱线值。
第二步数据流生成模块对线性调频脉冲信号进行采样，生成数据流
数据流生成模块在脉冲宽度的时间范围内对线性调频脉冲信号x(t)进行采样，生成数据流x(n)，采样频率为fs，fs≥4(f0+B/2)，其中f0为线性调频脉冲信号的中心频率，B为线性调频脉冲信号的调制带宽，f0+B/2为线性调频脉冲信号的最大频率。t为时间参数，n为时域点索引，n＝1,2,…,N，N为数据流长度，N＝fs·τ，τ为线性调频脉冲信号的脉冲宽度。
第三步频域点索引值生成模块获取中心频率处的频域点索引值
频域点索引值生成模块获取中心频率处所对应的频域点索引值k0：k0＝f0·N/fs，其中f0为线性调频脉冲信号的中心频率。
第四步中心频率处谱线值获取模块通过数据流和频域点索引值获取中心频率处谱线值
中心频率处谱线值获取模块使用数据流x(n)和频域点索引值k0进行递推，获取递推运算结果y(n)：
递推的初始条件为y(-1)＝0。其中y(n-1)为递推运算结果y(n)上一个数据点的值；WN为旋转因子，WN＝e-j(2π/N)，j表示e-j(2π/N)为复数表示形式；为复数固定值；表示将与y(n-1)进行复数相乘。
中心频率处谱线值获取模块通过递推获取y(N)的值，y(N)的值就是中心频率处的谱线值。
实施例2
一种线性调频脉冲信号中心频率处谱线值确定***，包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块。数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块顺次连接。数据流生成模块对线性调频脉冲信号进行采样，生成数据流；频域点索引值生成模块获取中心频率处的频域点索引值；中心频率处谱线值获取模块使用数据流和频域点索引值获取中心频率处谱线值。

Claims (3)

1.一种线性调频脉冲信号中心频率处谱线值确定方法，其特征在于具体步骤为：
第一步 搭建谱线值确定***
谱线值确定***，包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块；
数据流生成模块的功能为：对线性调频脉冲信号进行采样，生成数据流；
频域点索引值生成模块的功能为：获取中心频率处的频域点索引值；
中心频率处谱线值获取模块的功能为：使用数据流和频域点索引值计算中心频率处谱线值；
第二步 数据流生成模块对线性调频脉冲信号进行采样，生成数据流
数据流生成模块在脉冲宽度的时间范围内对线性调频脉冲信号x(t)进行采样，生成数据流x(n)，采样频率为fs，t为时间参数，n为时域点索引，n＝1,2,…,N，N为数据流长度，N＝fs·τ，τ为线性调频脉冲信号的脉冲宽度；
第三步 频域点索引值生成模块获取中心频率处的频域点索引值
频域点索引值生成模块获取中心频率处所对应的频域点索引值k0：k0＝f0·N/fs，其中f0为线性调频脉冲信号的中心频率；
第四步 中心频率处谱线值获取模块通过数据流和频域点索引值获取中心频率处谱线值
中心频率处谱线值获取模块使用数据流x(n)和频域点索引值k0进行递推，获取递推运算结果y(n)：
递推的初始条件为y(-1)＝0。其中y(n-1)为递推运算结果y(n)上一个数据点的值；WN为旋转因子，WN＝e-j(2π/N)，j表示e-j(2π/N)为复数表示形式；为复数固定值；表示将与y(n-1)进行复数相乘；
中心频率处谱线值获取模块通过递推获取y(N)的值，y(N)的值就是中心频率处的谱线值。
2.如权利要求1所述的线性调频脉冲信号中心频率处谱线值确定方法，其特征在于：所述第二步中采样频率fs≥4(f0+B/2)，其中f0为线性调频脉冲信号的中心频率，B为线性调频脉冲信号的调制带宽，f0+B/2为线性调频脉冲信号的最大频率。
3.一种线性调频脉冲信号中心频率处谱线值确定***，其特征在于包括：数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块；
数据流生成模块、频域点索引值生成模块和中心频率处谱线值获取模块顺次连接；
数据流生成模块对线性调频脉冲信号进行采样，生成数据流；频域点索引值生成模块获取中心频率处的频域点索引值；中心频率处谱线值获取模块使用数据流和频域点索引值获取中心频率处谱线值。
CN201710814543.0A 2017-09-10 2017-09-10 一种线性调频脉冲信号中心频率处谱线值确定方法 Active CN109490853B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710814543.0A CN109490853B (zh) 2017-09-10 2017-09-10 一种线性调频脉冲信号中心频率处谱线值确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710814543.0A CN109490853B (zh) 2017-09-10 2017-09-10 一种线性调频脉冲信号中心频率处谱线值确定方法

Publications (2)

Publication Number Publication Date
CN109490853A true CN109490853A (zh) 2019-03-19
CN109490853B CN109490853B (zh) 2022-11-22

Family

ID=65688785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710814543.0A Active CN109490853B (zh) 2017-09-10 2017-09-10 一种线性调频脉冲信号中心频率处谱线值确定方法

Country Status (1)

Country Link
CN (1) CN109490853B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115060355A (zh) * 2022-04-12 2022-09-16 东南大学 一种基于线性调频脉冲的谐振子品质因数测量方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100225524A1 (en) * 2009-03-03 2010-09-09 Wieslaw Jerzy Szajnowski Spectral analysis and fmcw automotive radar utilizing the same
CN102353838A (zh) * 2011-06-30 2012-02-15 南京航空航天大学 一种用fpga芯片实现的快速高精度频率测量方法
CN102624661A (zh) * 2012-03-13 2012-08-01 电子科技大学 基于谱峰突出指数的谱峰个数识别方法
CN103412287A (zh) * 2013-09-01 2013-11-27 西安电子科技大学 基于lvd的线性调频信号参数估计方法
CN106291516A (zh) * 2016-07-27 2017-01-04 河海大学 一种声呐应答式干扰的剔除方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100225524A1 (en) * 2009-03-03 2010-09-09 Wieslaw Jerzy Szajnowski Spectral analysis and fmcw automotive radar utilizing the same
CN102353838A (zh) * 2011-06-30 2012-02-15 南京航空航天大学 一种用fpga芯片实现的快速高精度频率测量方法
CN102624661A (zh) * 2012-03-13 2012-08-01 电子科技大学 基于谱峰突出指数的谱峰个数识别方法
CN103412287A (zh) * 2013-09-01 2013-11-27 西安电子科技大学 基于lvd的线性调频信号参数估计方法
CN106291516A (zh) * 2016-07-27 2017-01-04 河海大学 一种声呐应答式干扰的剔除方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张书仙等: "M-Rife频率估计法在测控信号载波捕获中的应用", 《遥测遥控》 *
张庆国: "基于FFT并实时修正的Rife频率估计算法", 《信号处理》 *
钱云襄等: "线性调频连续波信号参数估计算法", 《现代雷达》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115060355A (zh) * 2022-04-12 2022-09-16 东南大学 一种基于线性调频脉冲的谐振子品质因数测量方法
CN115060355B (zh) * 2022-04-12 2024-03-26 东南大学 一种基于线性调频脉冲的谐振子品质因数测量方法

Also Published As

Publication number Publication date
CN109490853B (zh) 2022-11-22

Similar Documents

Publication Publication Date Title
US20200011911A1 (en) High-precision frequency measuring system and method
CN105229431A (zh) 具有改进的距离确定的料位测量
CN106546817B (zh) 一种具有反馈功能的频率估计和能量估计方法
CN109342813B (zh) 一种基于dft和二分法的正弦信号频率估计方法
CN109490852B (zh) 一种线性调频脉冲信号调频斜率极性确定方法
CN103823177B (zh) 基于窗函数设计的滤波器性能检测方法和***
CN103744066A (zh) 数字正交鉴相与匹配滤波的优化方法
CN109490853A (zh) 一种线性调频脉冲信号中心频率处谱线值确定方法
CN107688167B (zh) 一种多时宽线性调频脉冲压缩信号幅度包络曲线生成方法
CN105745550B (zh) 模拟至信息转换器
CN102901951A (zh) 基于gpu的雷达信号脉内特征实时分析实现方案
CN110808929A (zh) 相减策略的实复转换式信噪比估计算法
US7444365B2 (en) Non-linear digital rank filtering of input signal values
CN114114231A (zh) 一种基于fft-czt的雷达测距方法
CN115630271A (zh) 信号频率估计方法、装置、设备及存储介质
CN108152795A (zh) 一种宽带线性调频脉冲信号初始频率估计方法
CN108169715B (zh) 一种同相和正交通道相位不平衡度确定方法及***
Ni et al. Polyphase-modulated radar signal recognition based on time-frequency amplitude and phase features
CN106597100B (zh) 一种广域电网动态频率的插值fft估计方法
US11531639B2 (en) Pipelined cognitive signal processor
CN112162153A (zh) 基于相位直线拟合的正弦波信号频率估计方法
CN108535542B (zh) 一种寻峰鉴相方法
CN102420792B (zh) 基于变长-迭代运算的单音信号频率估计方法
US10036769B2 (en) Apparatus for efficient measurement of tone frequency, amplitude, and phase
CN109672489A (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