CN111487462A - 一种超快速测频方法 - Google Patents

一种超快速测频方法 Download PDF

Info

Publication number
CN111487462A
CN111487462A CN202010314738.0A CN202010314738A CN111487462A CN 111487462 A CN111487462 A CN 111487462A CN 202010314738 A CN202010314738 A CN 202010314738A CN 111487462 A CN111487462 A CN 111487462A
Authority
CN
China
Prior art keywords
frequency measurement
frequency
sampling points
sequence
signal
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
CN202010314738.0A
Other languages
English (en)
Other versions
CN111487462B (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.)
8511 Research Institute of CASIC
Original Assignee
8511 Research Institute of CASIC
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 8511 Research Institute of CASIC filed Critical 8511 Research Institute of CASIC
Priority to CN202010314738.0A priority Critical patent/CN111487462B/zh
Publication of CN111487462A publication Critical patent/CN111487462A/zh
Application granted granted Critical
Publication of CN111487462B publication Critical patent/CN111487462B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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)

Abstract

本发明提供了一种超快速测频方法，包括以下步骤：构建N个采样点的测频查找表；采用高速AD进行信号采集和单比特量化；对系统收集的采样点进行滑动，对每次滑动后收集的N个采样点作为地址，在查找表中获取采样点的信号频率；对查表获取的信号频率进行平均，作为测量结果。

Description

一种超快速测频方法
技术领域
本发明涉及一种电子对抗技术，特别是一种超快速测频方法。
背景技术
如何快速对目标频率进行快速测量是电子对抗技术领域内的关键技术，随着电子对抗技术的不断发展，在电子战接收机领域内的新技术、新方法层出不穷，研究人员提出了具有高处理带宽、高处理速度的单比特测频方法，单比特测频方法基于单比特接收机技术，该技术通过采用位数较少的ADC、简化DFT运算的核函数ejωn(使其值为±1、±j)的方式，从而简化了乘法运算，提高了运算速度。但上述方法还是无法避免大量的加减法运算，影响了处理的实时性。
发明内容
本发明的目的在于提供一种超快速测频方法。
实现本发明目的的技术方案为：一种超快速测频方法，包括以下步骤：
构建N个采样点的测频查找表；采用高速AD进行信号采集和单比特量化；对系统收集的采样点进行滑动，对每次滑动后收集的N个采样点作为地址，在查找表中获取采样点的信号频率；对查表获取的信号频率进行平均，作为测量结果fout。
进一步地，测频查找表的构建具体包括以下步骤：
步骤11，对N点二进制序列进行去直流，即将0值转换为-1，然后对序列进行补零，形成长度为2a(a为正整数)的序列s(n)；
步骤12，对序列s(n)进行FFT运算，求出信号频谱S(k)，在S(k)中横坐标为[0,2a-1]的范围内寻找幅度最大点的位置kmax，并根据采样频率fs与FFT点数的相应关系，计算出测频结果(第一奈奎斯特区)fout；
步骤13，遍历所有N点二进制序列，以N点二进制序列为地址，测频结果为内容，生成测频查找表。
进一步地，
Figure BDA0002459106860000011
进一步地，
Figure BDA0002459106860000012
其中M为滑动的次数，kimax为第i次滑动后寻找幅度最大点的位置。
进一步地，经处理后的待测单比特信号为
Figure BDA0002459106860000021
其中，s(t)为输入的射频信号。
本发明与现有技术相比，其显著优点为：(1)通过查表的方式可以在极短的时间内给出采集信号的频率测量结果；(2)配合微波分频器可快速测量高频信号的频率；(3)利用快速测频结果引导接收机其他组件或算法模块工作。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为一种超快速测频方法流程图。
图2为高速AD采样序列与查找表地址对应关系示例示意图。
图3为采样滑动处理示意图。
图4为实例测频仿真结果示意图。
具体实施方式
结合图1和图2，本发明的一种超快速测频方法，包括以下步骤：
步骤1，构建N个采样点的测频查找表；
步骤2，采用高速AD进行信号采集和单比特量化，假设输入的射频信号可表示为s(t)，经处理后的待测单比特信号可表示为
Figure BDA0002459106860000022
步骤2，利用系统收集的N个采样点作为地址，查找步骤1构建的查找表获取信号频率；
步骤3，利用系统收集的N个采样点作为地址，查找步骤1构建的查找表获取信号频率；
步骤4，对系统收集采样点进行滑动，产生新的N个采样点作为地址，查找步骤1构建的查找表获取信号频率；
步骤5，多次重复步骤4，对步骤3和步骤4的查表结果进行平均，作为测量结果，即
Figure BDA0002459106860000031
本发明的N个采样点测频查找表的构建方法，包括以下步骤：
步骤11，对N点二进制序列进行去直流，即将0值转换为-1，然后对序列进行补零，形成长度为2a(a为正整数)的序列s(n)；
步骤12，对序列s(n)进行FFT运算，求出信号频谱S(k)，在S(k)中横坐标为[0,2a-1]的范围内寻找幅度最大点的位置kmax，并根据采样频率fs与FFT点数的相应关系，计算出测频结果(第一奈奎斯特区)fout。
Figure BDA0002459106860000032
步骤13，遍历所有N点二进制序列，以N点二进制序列为地址，测频结果为内容，生成测频查找表。
实施例一
结合图1至图4，各相关参数的取值为：N＝15，a＝9，nstep＝2，M＝64，fs＝2500MHz，输入信号频率范围为168.75～1125MHz，仿真信噪比为12dB。根据本专利公开的一种超快速测频方法，其步骤具体如下：
步骤1，构建15个采样点的测频查找表；
步骤11，对15点二进制序列进行去直流，即将0值转换为-1，然后对序列进行补零，形成长度为做2a＝512点的序列s(n)；
步骤12，对序列s(n)进行FFT运算，求出信号频谱S(k)，在S(k)中横坐标为[0,255]的范围内寻找幅度最大点的位置kmax，并根据采样频率fs与FFT点数的相应关系，计算出测频结果(第一奈奎斯特区)fout。
Figure BDA0002459106860000033
步骤13，遍历所有15点二进制序列，以15点二进制序列为地址，测频结果为内容，生成测频查找表。
步骤2，采用高速AD进行信号采集和单比特量化，产生0/1取值的二进制采样序列；
步骤3，利用系统收集的15个采样点作为地址，查找步骤1构建的查找表获取信号频率；
步骤4，将步骤3的采样窗口滑动2个采样点，产生新的15个采样点作为地址，查找步骤1构建的查找表获取信号频率；
步骤5，重复63次步骤4，对步骤3和步骤4的查表结果进行平均，作为测量结果。
实例测频时间分析：实际采用数据点数为143个，数据点采样耗时57.2ns，以156.25MHz工作时钟为例需要10个时钟周期，Rom表读出延迟1个时钟周期，64次数据相加采用8路串行操作5个时钟周期，也即17个工作时钟后可以输出测频索引，总耗时17*6.4＝108.8ns。

Claims (5)

1.一种超快速测频方法，其特征在于，包括以下步骤：
构建N个采样点的测频查找表；
采用高速AD进行信号采集和单比特量化；
对系统收集的采样点进行滑动，对每次滑动后收集的N个采样点作为地址，在查找表中获取采样点的信号频率；
对查表获取的信号频率进行平均，作为测量结果
Figure FDA0002459106850000011
2.根据权利要求1所述的方法，其特征在于，测频查找表的构建具体包括以下步骤：
步骤11，对N点二进制序列进行去直流，即将0值转换为-1，然后对序列进行补零，形成长度为2a(a为正整数)的序列s(n)；
步骤12，对序列s(n)进行FFT运算，求出信号频谱S(k)，在S(k)中横坐标为[0,2a-1]的范围内寻找幅度最大点的位置kmax，并根据采样频率fs与FFT点数的相应关系，计算出测频结果(第一奈奎斯特区)fout；
步骤13，遍历所有N点二进制序列，以N点二进制序列为地址，测频结果为内容，生成测频查找表。
3.根据权利要求2所述的方法，其特征在于，
Figure FDA0002459106850000012
4.根据权利要求3所述的方法，其特征在于，
Figure FDA0002459106850000013
其中M为滑动的次数，kimax为第i次滑动后寻找幅度最大点的位置。
5.根据权利要求1所述的方法，其特征在于，经处理后的待测单比特信号为
Figure FDA0002459106850000014
其中，s(t)为输入的射频信号。
CN202010314738.0A 2020-04-21 2020-04-21 一种超快速测频方法 Active CN111487462B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010314738.0A CN111487462B (zh) 2020-04-21 2020-04-21 一种超快速测频方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010314738.0A CN111487462B (zh) 2020-04-21 2020-04-21 一种超快速测频方法

Publications (2)

Publication Number Publication Date
CN111487462A true CN111487462A (zh) 2020-08-04
CN111487462B CN111487462B (zh) 2022-05-13

Family

ID=71812923

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010314738.0A Active CN111487462B (zh) 2020-04-21 2020-04-21 一种超快速测频方法

Country Status (1)

Country Link
CN (1) CN111487462B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114062782A (zh) * 2021-10-21 2022-02-18 中国电子科技集团公司第二十九研究所 适用于宽带射频信号谱估计的2比特采样量化系统及方法
CN114089290A (zh) * 2021-11-12 2022-02-25 扬州宇安电子科技有限公司 一种快速高效测频引导噪声干扰方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4928105A (en) * 1988-05-25 1990-05-22 Telemus Electronic Systems, Inc. Intrapulse radar receiver
JP2005181177A (ja) * 2003-12-22 2005-07-07 Kokusai Denki Engineering:Kk 周波数検出方法
CN101057150A (zh) * 2004-09-13 2007-10-17 菲尔特罗尼克公开有限公司 用于测量接收的信号的频率的方法和设备
US20100069030A1 (en) * 2006-02-10 2010-03-18 Thierry Briand Frequency measuring broadband digital receiver
CN110068729A (zh) * 2019-04-22 2019-07-30 南京磐能电力科技股份有限公司 一种信号相量计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4928105A (en) * 1988-05-25 1990-05-22 Telemus Electronic Systems, Inc. Intrapulse radar receiver
JP2005181177A (ja) * 2003-12-22 2005-07-07 Kokusai Denki Engineering:Kk 周波数検出方法
CN101057150A (zh) * 2004-09-13 2007-10-17 菲尔特罗尼克公开有限公司 用于测量接收的信号的频率的方法和设备
US20100069030A1 (en) * 2006-02-10 2010-03-18 Thierry Briand Frequency measuring broadband digital receiver
CN110068729A (zh) * 2019-04-22 2019-07-30 南京磐能电力科技股份有限公司 一种信号相量计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王坤达: "基于实时数字信号处理的宽带单比特瞬时测频接收技术", 《舰船电子对抗》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114062782A (zh) * 2021-10-21 2022-02-18 中国电子科技集团公司第二十九研究所 适用于宽带射频信号谱估计的2比特采样量化系统及方法
CN114062782B (zh) * 2021-10-21 2023-04-28 中国电子科技集团公司第二十九研究所 适用于宽带射频信号谱估计的2比特采样量化系统及方法
CN114089290A (zh) * 2021-11-12 2022-02-25 扬州宇安电子科技有限公司 一种快速高效测频引导噪声干扰方法

Also Published As

Publication number Publication date
CN111487462B (zh) 2022-05-13

Similar Documents

Publication Publication Date Title
WO2018188228A1 (zh) 高精度频率测量系统及方法
CN111487462B (zh) 一种超快速测频方法
CN105572499B (zh) 一种基于矢量网络分析仪的眼图生成方法
CN112737711B (zh) 一种基于自适应噪声基底估计的宽带载波检测方法
CN105223482B (zh) 局部放电特高频信号波形的小波分解二值去噪方法
CN102109555B (zh) 具有数字中频信号处理系统的频谱分析仪及其实现方法
CN103956756B (zh) 一种电力系统低频振荡模态辨识方法
CN103457603A (zh) 一种基于平均频谱测试adc动态参数的方法
CN103983849B (zh) 一种实时高精度的电力谐波分析方法
CN108957175B (zh) 基于改进的hht算法的电能质量扰动识别方法
CN113269169A (zh) 一种轴承故障检测方法和装置
CN114374446B (zh) 一种基于脉冲信号的幅相特性测量方法
CN102200550B (zh) 一种用于高精度测量相位差的延迟正交数字中频鉴相方法
CN111865865A (zh) 适用于高灵敏度星载ads-b接收机的频偏及相偏估计方法
Roy et al. Fast OMP algorithm and its FPGA implementation for compressed sensing‐based sparse signal acquisition systems
CN108760890B (zh) 一种使用信号能谱系数监测脆性材料损伤演化过程的方法
CN112764000A (zh) 基于分位点的对数正态纹理海杂波幅度模型参数估计方法
CN105652254A (zh) 室外场rcs测量方法及系统
CN103095249A (zh) 一种中值滤波电路及方法
CN115438790B (zh) 量子态信息处理系统、量子测控系统、量子计算机
CN113406453B (zh) 一种基于mcu的prpd/prps图谱数据处理方法及检测装置
CN109085545A (zh) 基于fista的压缩感知信号实时重构方法
CN101964635B (zh) 数字信号的自动增益控制方法
CN107576842B (zh) 一种宽频带同步采样方法
Qi et al. Research on EMD improvement methods

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