CN105004921B - 基于李沙育图形的异频高精度频率测量方法 - Google Patents
基于李沙育图形的异频高精度频率测量方法 Download PDFInfo
- Publication number
- CN105004921B CN105004921B CN201510485837.4A CN201510485837A CN105004921B CN 105004921 B CN105004921 B CN 105004921B CN 201510485837 A CN201510485837 A CN 201510485837A CN 105004921 B CN105004921 B CN 105004921B
- Authority
- CN
- China
- Prior art keywords
- frequency
- lissajouf
- signal
- measured signal
- module
- 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
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000003550 marker Substances 0.000 claims description 19
- 230000000737 periodic effect Effects 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000003786 synthesis reaction Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004891 communication Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 8
- 238000000691 measurement method Methods 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 2
- 235000002918 Fraxinus excelsior Nutrition 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 239000002956 ash Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Landscapes
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
Abstract
本发明公开了一种基于李沙育图形的异频高精度频率测量方法，依次包括以下步骤：A：将参考信号和被测信号送入FPGA模块，然后将被测信号的粗测频率通过FPGA模块生成参考频率信号；B：参考频率信号和被测信号同时输入数字示波器产生动态李沙育图形；C：通过高分辨率CCD图像采集设备采集动态李沙育图形，并对动态李沙育图形的翻转周期进行计算，利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号；D：上位机利用调整过的参考频率信号以及动态李沙育图形的翻转周期，计算出被测模块输出的被测信号的频率。本发明能够准确测量被测信号的频率值，系统稳定可靠，测量精度高。
Description
技术领域
本发明涉及一种频率测量方法，尤其涉及一种基于李沙育图形的异频高精度频率测量方法。
背景技术
目前,传统的频率测量方法以及基于异频相位处理的频率测量方法，其精度受制于检测器件的分辨率。对于一般的TTL门电路器件来说，当输入信号为窄脉冲且其持续接近Tpd(门电路延时时间一般为ns量级)时，输出的变化跟不上输入的变化，无法保证正确的逻辑输出，导致产生相位重合点模糊区，因此分辨率难以达到ns量级。李沙育图形的形状特征由两频率信号之间的相位关系所决定，理论上不存在测量误差，其精度主要受制于CCD图像采集卡的采样速度，因此能够达到现有频率测量电路所无法实现的高精度。
发明内容
本发明的目的是提供一种基于李沙育图形的异频高精度频率测量方法，能够通过李沙育图形的翻转周期计算出参考信号和被测信号之间的微小频差，从而准确测量被测信号的频率值，本发明无需混频、倍频等复杂频率变换电路，电路结构简单，操作方便，成本低廉，系统稳定可靠，测量精度高。
本发明采用下述技术方案：
基于李沙育图形的异频高精度频率测量方法，依次包括以下步骤：
A：将参考模块输出的参考信号送入FPGA模块，同时将被测模块输出的被测信号送入FPGA模块；然后结合参考信号，利用FPGA模块对被测信号进行粗测得到被测信号的粗测频率，然后将被测信号的粗测频率通过FPGA模块生成参考频率信号；
B：将FPGA模块生成的参考频率信号和被测模块输出的被测信号同时输入数字示波器，利用数字示波器产生动态李沙育图形；
C：通过高分辨率CCD图像采集设备将示波器上显示的动态李沙育图形实时采集到上位机中，利用上位机对动态李沙育图形的翻转周期进行计算，并将计算结果反馈至FPGA模块，利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号；
D：上位机利用调整过的参考频率信号以及动态李沙育图形的翻转周期，计算出被测模块输出的被测信号的频率，被测信号的频率计算公式如下：
fx＝f′x+Δf
其中，Δf是被测信号与参考频率信号的频率偏差；A是参考频率信号与被测信号标称值的最大频率公约数，A≡1；τ是动态李沙育图形的翻转周期值，fx是被测信号的频率，f′x是被测信号频率的粗测频率。
所述的步骤A中，利用FPGA内部的计数器，在同一闸门时间内对参考信号和被测信号进行计数，以获得参考信号和被测信号的计数值，然后结合参考信号的周期，经FPGA内部配置的MCU对计数值进行处理，从而获得被测信号的粗测频率。
所述的步骤A中，被测信号的粗测频率被送入FPGA内部配置的控制模块，以产生用于对FPGA内部配置的频率合成模块进行控制的控制频率信号，在控制频率信号的控制下，FPGA内部配置的频率合成模块产生与被测信号的粗测频率相同数值的参考频率信号。
所述的步骤C中，高分辨率CCD图像采集设备是指具有高清成像功能的工业摄像机。
所述的步骤C中，翻转周期结果通过串行通信反馈到FPGA模块，利用FPGA模块对参考频率信号进行调整。如果翻转周期的值在0.5s至30s范围内，能够通过高分辨率CCD图像采集并利用上位机计算出来，则直接得到被测频率的粗测值偏差并显示被测频率值；如果翻转周期小于0.5s且不能通过高分辨率CCD图像采集并利用上位机计算出来，则将翻转周期结果通过串行通信反馈到FPGA内部配置的频率合成模块和控制模块，以控制参考信号频率产生步进式变化，并实时观察动态李沙育图形，直至动态李沙育图形稳定并使翻转周期在0.5s至30s范围内，从而便于通过CCD图像采集及计算。上述方法即可实现利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号。
所述的步骤C中，动态李沙育图形的翻转周期计算方法如下：
C1：对于李沙育图形中的第i帧Fi(m,n)，m和n分别表示第i帧图像的行和列；忽略第i帧Fi(m,n)中的色彩和饱和度信息，只保留亮度信息，并将亮度信息转变为具有256个灰度级的灰度图像I(m,n)；
C2：对灰度图像I(m,n)进行分割得到二值图像J(m,n)：
其中，Threshold为阈值，阈值根据图像采集质量进行设定；
C3：将二值图像J(m,n)中位于李沙育图形上的像素和背景上的像素分别用1和0表示，则二值图像J(m,n)的归一化面积为
其中，NAi表示第i帧图像的归一化面积，M×N表示李沙育图形视频的分辨率；
C4：对李沙育图形视频中所有的帧均进行步骤C1至步骤C3的处理，得到归一化面积实验曲线NA(i)，1≤i≤NOF，其中NOF表示李沙育图形视频中的帧数；
C5：将归一化面积实验曲线NA(i)所产生的重合程度最大的帧设为翻转周期的标记帧，第i个标记帧与第i+2个标记帧之间的时间间隔即为李沙育图形的翻转周期值，i为任意自然数；
设某段李沙育图形视频中总共有N个标记帧，第l个标记帧为李沙育图形视频中的第kl帧，则可以得到N-2个翻转周期值，那么第p个翻转周期值为：
其中，τp为第p个翻转周期值，kp+2为第p+2个标记帧，kp为第p个标记帧，FR为李沙育图形视频的帧率。
所述的步骤C中，上位机采用LabVIEW虚拟仪器。
本发明基于李沙育图形的翻转周期计算出参考信号和被测信号之间的微小频差，最终获得被测信号的频率值。本发明无需混频、倍频等复杂的频率变换电路，不存在±1个字的频率计数误差，以软件算法处理和硬件实现为主，电路结构简单，操作方便，成本低廉，系统稳定可靠，测量精度高。
附图说明
图1为本发明的流程示意图。
具体实施方式
以下结合附图和实施例对本发明作以详细的描述：
如图1所示，本发明所述的基于李沙育图形的异频高精度频率测量方法，依次包括以下步骤：
A：将参考模块输出的参考信号送入FPGA模块，同时将被测模块输出的被测信号送入FPGA模块；然后结合参考信号，利用FPGA模块对被测信号进行粗测得到被测信号的粗测频率，然后将被测信号的粗测频率通过FPGA模块生成参考频率信号；
粗测频率通过直接测频法获得。利用FPGA内部的计数器，在同一闸门时间内对参考信号和被测信号进行计数，以获得参考信号和被测信号的计数值，然后结合参考信号的周期，经FPGA内部配置的MCU对计数值进行处理，从而获得被测信号的粗测频率。由于直接测频法在测频过程中存在±1字的计数误差，故所得频率称为粗测频率。
参考频率信号的生成通过FPGA内部配置的频率合成模块和控制模块完成。被测信号的粗测频率被送入FPGA内部配置的控制模块，以产生用于对FPGA内部配置的频率合成模块进行控制的控制频率信号，在控制频率信号的控制下，FPGA内部配置的频率合成模块产生与被测信号的粗测频率相同数值的参考频率信号，参考频率信号与被测频率信号客观上具有微小的频率偏差，使参考频率信号与被测频率信号的直接相位比对结果更具有规律性。
B：将FPGA模块生成的参考频率信号和被测模块输出的被测信号同时输入数字示波器，利用数字示波器产生动态李沙育图形；
C：通过高分辨率CCD图像采集设备将示波器上显示的动态李沙育图形实时采集到上位机中，利用上位机对动态李沙育图形的翻转周期进行计算，并将计算结果反馈至FPGA模块，利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号；
所述的高分辨率CCD图像采集设备是指具有高清成像功能的工业摄像机且分辨率大于1000万像素。
翻转周期结果通过串行通信反馈到FPGA模块，利用FPGA模块对参考频率信号进行调整。如果翻转周期的值在0.5s至30s范围内，能够通过高分辨率CCD图像采集并利用上位机计算出来，则直接得到被测频率的粗测值偏差并显示被测频率值；如果翻转周期小于0.5s且不能通过高分辨率CCD图像采集并利用上位机计算出来，则将翻转周期结果通过串行通信反馈到FPGA内部配置的频率合成模块和控制模块，以控制参考信号频率产生步进式变化，并实时观察动态李沙育图形，直至动态李沙育图形稳定并使翻转周期在0.5s至30s范围内，从而便于通过CCD图像采集及计算。上述方法即可实现利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号。
步骤C中，动态李沙育图形的翻转周期计算方法如下：
C1：对于李沙育图形中的第i帧Fi(m,n)，m和n分别表示第i帧图像的行和列；忽略第i帧Fi(m,n)中的色彩和饱和度信息，只保留亮度信息，并将亮度信息转变为具有256个灰度级的灰度图像I(m,n)；
C2：对灰度图像I(m,n)进行分割得到二值图像J(m,n)：
其中，Threshold为阈值，阈值根据图像采集质量进行设定；
C3：将二值图像J(m,n)中位于李沙育图形上的像素和背景上的像素分别用1和0表示，则二值图像J(m,n)的归一化面积为
其中，NAi表示第i帧图像的归一化面积，M×N表示李沙育图形视频的分辨率；
C4：对李沙育图形视频中所有的帧均进行步骤C1至步骤C3的处理，得到归一化面积实验曲线NA(i)，1≤i≤NOF，其中NOF表示李沙育图形视频中的帧数；
C5：将归一化面积实验曲线NA(i)所产生的重合程度最大的帧设为翻转周期的标记帧，第i个标记帧与第i+2个标记帧之间的时间间隔即为李沙育图形的翻转周期值，i为任意自然数；
设某段李沙育图形视频中总共有N个标记帧，第l个标记帧为李沙育图形视频中的第kl帧，则可以得到N-2个翻转周期值，那么第p个翻转周期值为：
其中，τp为第p个翻转周期值，kp+2为第p+2个标记帧，kp为第p个标记帧，FR为李沙育图形视频的帧率。
D：上位机利用调整过的参考频率信号以及动态李沙育图形的翻转周期，计算出被测模块输出的被测信号的频率。
被测信号的频率计算如下：
fx＝f′x+Δf
其中，Δf是被测信号与参考频率信号的频率偏差；A是参考频率信号与被测信号标称值的最大频率公约数，本发明中A≡1；τ是动态李沙育图形的翻转周期值，fx是被测信号的频率；f′x是被测信号频率的粗测频率。
本发明中，上位机也可采用LabVIEW虚拟仪器。
下表为使用本发明所述的基于李沙育图形的异频高精度频率测量方法的实测结果：
综上，传统的李沙育图形只是单纯用于观察两频率信号之间的倍数关系，而本发明基于李沙育图形的翻转周期计算出参考信号和被测信号之间的微小频差，最终获得被测信号的频率值。本发明无需混频、倍频等复杂的频率变换电路，不存在±1个字的频率计数误差，以软件算法处理和硬件实现为主，电路结构简单，操作方便，成本低廉，系统稳定可靠，测量精度高。
Claims (6)
1.基于李沙育图形的异频高精度频率测量方法，其特征在于：依次包括以下步骤：
A：将参考模块输出的参考信号送入FPGA模块，同时将被测模块输出的被测信号送入FPGA模块；然后结合参考信号，利用FPGA模块对被测信号进行粗测得到被测信号的粗测频率，然后将被测信号的粗测频率通过FPGA模块生成参考频率信号；
B：将FPGA模块生成的参考频率信号和被测模块输出的被测信号同时输入数字示波器，利用数字示波器产生动态李沙育图形；
C：通过高分辨率CCD图像采集设备将示波器上显示的动态李沙育图形实时采集到上位机中，利用上位机对动态李沙育图形的翻转周期进行计算，并将计算结果反馈至FPGA模块，利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号；
其中，动态李沙育图形的翻转周期计算方法如下：
C1：对于李沙育图形中的第i帧Fi(m,n)，m和n分别表示第i帧图像的行和列；忽略第i帧Fi(m,n)中的色彩和饱和度信息，只保留亮度信息，并将亮度信息转变为具有256个灰度级的灰度图像I(m,n)；
C2：对灰度图像I(m,n)进行分割得到二值图像J(m,n)：
其中，Threshold为阈值，阈值根据图像采集质量进行设定；
C3：将二值图像J(m,n)中位于李沙育图形上的像素和背景上的像素分别用1和0表示，则二值图像J(m,n)的归一化面积为
其中，NAi表示第i帧图像的归一化面积，M×N表示李沙育图形视频的分辨率；
C4：对李沙育图形视频中所有的帧均进行步骤C1至步骤C3的处理，得到归一化面积实验曲线NA(i)，1≤i≤NOF，其中NOF表示李沙育图形视频中的帧数；
C5：将归一化面积实验曲线NA(i)所产生的重合程度最大的帧设为翻转周期的标记帧，第i个标记帧与第i+2个标记帧之间的时间间隔即为李沙育图形的翻转周期值，i为任意自然数；设某段李沙育图形视频中总共有N个标记帧，第l个标记帧为李沙育图形视频中的第kl帧，则可以得到N-2个翻转周期值，那么第p个翻转周期值为：
其中，τp为第p个翻转周期值，kp+2为第p+2个标记帧，kp为第p个标记帧，FR为李沙育图形视频的帧率；
D：上位机利用调整过的参考频率信号以及动态李沙育图形的翻转周期，计算出被测模块输出的被测信号的频率，被测信号的频率计算公式如下：
其中，Δf是被测信号与参考频率信号的频率偏差；A是参考频率信号与被测信号标称值的最大频率公约数，A≡1；τ是动态李沙育图形的翻转周期值，fx是被测信号的频率，fx′是被测信号频率的粗测频率。
2.根据权利要求1所述的基于李沙育图形的异频高精度频率测量方法，其特征在于：所述的步骤A中，利用FPGA内部的计数器，在同一闸门时间内对参考信号和被测信号进行计数，以获得参考信号和被测信号的计数值，然后结合参考信号的周期，经FPGA内部配置的MCU对计数值进行处理，从而获得被测信号的粗测频率。
3.根据权利要求1所述的基于李沙育图形的异频高精度频率测量方法，其特征在于：所述的步骤A中，被测信号的粗测频率被送入FPGA内部配置的控制模块，以产生用于对FPGA内部配置的频率合成模块进行控制的控制频率信号，在控制频率信号的控制下，FPGA内部配置的频率合成模块产生与被测信号的粗测频率相同数值的参考频率信号。
4.根据权利要求1所述的基于李沙育图形的异频高精度频率测量方法，其特征在于：所述的步骤C中，高分辨率CCD图像采集设备是指具有高清成像功能的工业摄像机。
5.根据权利要求1所述的基于李沙育图形的异频高精度频率测量方法，其特征在于：所述的步骤C中，翻转周期结果通过串行通信反馈到FPGA模块，利用FPGA模块对参考频率信号进行调整； 如果翻转周期的值在0.5s至30s范围内，能够通过高分辨率CCD图像采集并利用上位机计算出来，则直接得到被测频率的粗测值偏差并显示被测频率值；如果翻转周期小于0.5s且不能通过高分辨率CCD图像采集并利用上位机计算出来，则将翻转周期结果通过串行通信反馈到FPGA内部配置的频率合成模块和控制模块，以控制参考信号频率产生步进式变化，并实时观察动态李沙育图形，直至动态李沙育图形稳定并使翻转周期在0.5s至30s范围内，从而便于通过CCD图像采集及计算；上述方法即可实现利用FPGA模块对参考频率信号进行调整，使参考频率信号能够实时跟踪被测模块输出的被测信号。
6.根据权利要求1所述的基于李沙育图形的异频高精度频率测量方法，其特征在于：所述的步骤C中，上位机采用LabVIEW虚拟仪器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510485837.4A CN105004921B (zh) | 2015-08-10 | 2015-08-10 | 基于李沙育图形的异频高精度频率测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510485837.4A CN105004921B (zh) | 2015-08-10 | 2015-08-10 | 基于李沙育图形的异频高精度频率测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105004921A CN105004921A (zh) | 2015-10-28 |
| CN105004921B true CN105004921B (zh) | 2017-08-15 |
Family
ID=54377664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510485837.4A Expired - Fee Related CN105004921B (zh) | 2015-08-10 | 2015-08-10 | 基于李沙育图形的异频高精度频率测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105004921B (zh) |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB532392A (en) * | 1939-07-21 | 1941-01-23 | Standard Telephones Cables Ltd | Improvements in or relating to electric frequency testing |
| JPS6053854A (ja) * | 1983-09-02 | 1985-03-27 | Anritsu Corp | 周波数測定装置 |
-
2015
- 2015-08-10 CN CN201510485837.4A patent/CN105004921B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105004921A (zh) | 2015-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104048744B (zh) | 一种非接触式的基于影像的实时在线振动测量方法 | |
| CN107274453A (zh) | 一种结合标定与校正的摄像机三维测量装置、系统及方法 | |
| CN103837741A (zh) | 基于fpga的等精度频率测试系统及其设计方法 | |
| CN104219464B (zh) | 一种采样位置自适应调整的ccd视频信号处理系统 | |
| CN104991118B (zh) | 一种高分辨异频信号频率测量系统和测量方法 | |
| CN102591996B (zh) | 用于多部量测仪器的数据同步架构及方法 | |
| CN102928677A (zh) | 一种纳米级脉冲信号采集方法 | |
| CN106662472A (zh) | 偏差估计装置和方法以及故障诊断装置和方法 | |
| CN108318053A (zh) | 一种空间光学遥感相机成像时刻标定精度测量方法及系统 | |
| CN105336635B (zh) | Cd-sem装置的校正方法、应用cd-sem装置的方法及cd-sem装置 | |
| CN116124393A (zh) | 一种离轴测量时桥梁多点动挠度测量方法和装置 | |
| CN201935964U (zh) | 石英晶体振荡器仪表智能化自动检测装置 | |
| CN109191386A (zh) | 一种基于BPNN的快速Gamma校正方法及装置 | |
| CN104155664A (zh) | 一种星载接收机校时功能测试系统及方法 | |
| CN205593617U (zh) | 具有显示流量数据识别的lng加气机计量检定标准装置 | |
| CN108594147A (zh) | 一种模拟信号和数字信号同步采集及同步时间差校准方法 | |
| CN104678338A (zh) | 数字示波器及自检方法、自动测试系统 | |
| CN105004921B (zh) | 基于李沙育图形的异频高精度频率测量方法 | |
| CN103487692A (zh) | 一种触摸屏的质量检测方法 | |
| US10548559B2 (en) | Method for calibrating working plane of medical detection apparatus | |
| CN204795439U (zh) | 用于大畸变图像采集系统的预畸变测试卡成型系统 | |
| CN109709071A (zh) | 基于合成rgb相移彩色条纹背景纹影的流动显示装置和方法 | |
| CN116243583A (zh) | 用于tdl-tdc的神经网络测量校准系统和方法 | |
| CN103513103A (zh) | 物理量测量设备和物理量测量方法 | |
| TWI494549B (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 | ||
| GR01 | Patent grant | ||
| 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: 20170815 |