CN115356906B - 一种线性光学采样的双阈值拟合方法及时间偏差估计方法 - Google Patents
一种线性光学采样的双阈值拟合方法及时间偏差估计方法 Download PDFInfo
- Publication number
- CN115356906B CN115356906B CN202210826088.7A CN202210826088A CN115356906B CN 115356906 B CN115356906 B CN 115356906B CN 202210826088 A CN202210826088 A CN 202210826088A CN 115356906 B CN115356906 B CN 115356906B
- Authority
- CN
- China
- Prior art keywords
- los
- threshold
- pulse
- signal
- check
- 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
- 230000003287 optical effect Effects 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005070 sampling Methods 0.000 title claims abstract description 23
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 210000001520 comb Anatomy 0.000 claims description 8
- 230000003321 amplification Effects 0.000 claims description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 5
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 7
- 230000000630 rising effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04D—APPARATUS OR TOOLS SPECIALLY DESIGNED FOR MAKING OR MAINTAINING CLOCKS OR WATCHES
- G04D7/00—Measuring, counting, calibrating, testing or regulating apparatus
- G04D7/002—Electrical measuring and testing apparatus
- G04D7/003—Electrical measuring and testing apparatus for electric or electronic clocks
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Communication System (AREA)
Abstract
本发明公开了一种线性光学采样的双阈值拟合方法及时间偏差估计方法。本发明拟合方法为：1)对待拟合LOS信号进行检测；2)当检测到LOS信号电压值高于脉冲检测阈值Vthreshold时，记录下当前时刻tth1；3)当检测到LOS信号电压值高于脉冲达峰判定阈值Vcheck时，记录下当前时刻tcheck；当再次检测到LOS信号电压值低于Vthreshold时，记录当前时刻tth2；4)如果满足tcheck＜tth2，则(tth1,tth2)区间对应一个完整的LOS脉冲信号；如果不满足tcheck＜tth2，则继续检测，直到LOS信号电压值再次低于Vthreshold时，记录当前时刻tth3，此时(tth1,tth3)区间对应一个完整的LOS脉冲信号；5)根据完整的LOS脉冲信号对应的时间区间拟合得到LOS信号峰值对应的中心时刻t。
Description
技术领域
本发明涉及信息科学通信中高精度时间测量领域，尤其涉及一种线性光学采样的双阈值拟合方法及时间偏差估计方法。
背景技术
时间是极为重要的基本物理量，它不仅在日常生产生活中有着重要的应用，更在军事、前沿研究等领域扮演着重要角色。得益于高精度原子钟、光钟的飞速发展，时间信号已经成为了目前测量精度最高的物理量之一，被广泛用作参考物理量。
其中，时间信号的高精度测量主要指的是高精度时间偏差测量。目前比较成熟的方案有双混频时间偏差测量方案，其基本的思想是通过引入第三个公共频率源，并与两台待测频率源存在一个微小的频率差，来实现两台待测频率源微小时间差的放大。放大后的时间差容易被一般的时间探测器探测得到，之后再除以放大倍数，即可得到高精度的两台待测频率源的时间偏差。
近年来，研究人员又提出了线性光学采样的方案(也称作双光梳技术)来实现高精度的时间偏差测量。其基本思想也是通过引入第三台光梳，其频率与待测光梳的频率fr相差一个小量Δfr，来实现两台待测光梳脉冲信号时间间隔的放大。对于得到的第三台光梳分别与两台待测光梳进行线性光学采样所得到的线性光学采样信号(LOS信号)，需要对其进行拟合，以找出中心时刻，从而精确地计算出放大后的时间间隔，然后除以放大倍数，即fr/Δfr，从而得到精确的两台待测光梳的时间偏差。
另一方面，由于实际的应用场景是测量出时间偏差后，实时地进行两台待测光梳的同步，因此对于采集得到的LOS信号这一脉冲序列，需要不断地进行高精度的拟合操作，由于其中的脉冲脉宽远小于脉冲周期，因此对于整个周期进行拟合会极大地浪费资源、降低速度。
发明内容
针对现有技术中存在的技术问题，本发明的目的在于提出一种线性光学采样的双阈值拟合方法及时间偏差估计方法。本发明提出了双阈值探测的方法对LOS脉冲信号进行拟合，来实现对于拟合区间的合理精准选取。
本发明的技术方案为：
一种线性光学采样的双阈值拟合方法，其步骤包括：
1)对于待拟合的LOS信号进行检测；其中，参考光梳的频率与待测光梳的频率fr相差设定量Δfr；将参考光梳与待测光梳进行线性光学采样得到所述LOS信号；
2)当检测到LOS信号电压值高于脉冲检测阈值Vthreshold时，记录下当前时刻tth1作为LOS信号的开始；
3)当检测到LOS信号电压值高于脉冲达峰判定阈值Vcheck时，记录下当前时刻tcheck；当再次检测到LOS信号电压值低于脉冲检测阈值Vthreshold时，记录下当前时刻tth2；
4)如果满足tcheck＜tth2，则认为(tth1,tth2)区间对应一个完整的LOS脉冲信号；如果不满足tcheck＜tth2，则继续检测，直到LOS信号电压值再次低于脉冲检测阈值Vthreshold时，记录下当前时刻tth3，此时(tth1,tth3)区间对应一个完整的LOS脉冲信号；
5)根据完整的LOS脉冲信号对应的时间区间拟合得到LOS信号峰值对应的中心时刻t。一种线性光学采样的双阈值时间偏差估计方法，其步骤包括:
1)对于待拟合的两列LOS信号分别进行检测；其中，将参考光梳与第一待测光梳进行线性光学采样所得LOS信号记为第一列LOS信号，将参考光梳与第二待测光梳进行线性光学采样所得LOS信号记为第二列LOS信号；参考光梳的频率与待测光梳的频率fr相差设定量Δfr；
2)当检测到第一列LOS信号电压值高于脉冲检测阈值Vthreshold时，记录下当前时刻tth1作为第一列LOS信号的开始；
3)当检测到第一列LOS信号电压值高于脉冲达峰判定阈值Vcheck时，记录下当前时刻tcheck；当再次检测到第一列LOS信号电压值低于脉冲检测阈值Vthreshold时，记录下当前时刻tth2；
4)如果满足tcheck＜tth2，则认为(tth1,tth2)区间对应一个完整的LOS脉冲信号；如果不满足tcheck＜tth2，则继续检测，直到第一列LOS信号电压值再次低于脉冲检测阈值Vthreshold时，记录下当前时刻tth3，此时(tth1,tth3)区间对应一个完整的LOS脉冲信号；
5)根据完整的LOS脉冲信号对应的时间区间拟合得到第一列LOS信号峰值对应的中心时刻t1；
6)按照步骤2～5)的方法拟合得到第二列LOS信号峰值对应的中心时刻t2；
7)根据中心时刻t1、t2的差值得出放大后的时间间隔，将所述时间间隔除以对应的放大倍数得到两台待测光梳的时间偏差。
进一步的，获取连续n个LOS脉冲信号对应的时间区间(tth1-nts,tth2+nts)；步骤5)中，根据(tth1-nts,tth2+nts)拟合得到LOS信号峰值对应的中心时刻t；其中，ts表示单次采样的采样时间，n为拟合点数。
进一步的，LOS信号中的噪底电压为5mV，脉冲检测阈值Vthreshold为20mV，脉冲达峰判定阈值Vcheck为200mV。
进一步的，利用高斯拟合得到LOS信号峰值对应的中心时刻。
与现有技术相比，本发明的积极效果为：
通过引入双阈值判定，第一去除了未出现LOS脉冲时的大量无用数据，节省了计算资源，提升了拟合速度；第二脉冲检测阈值和脉冲达峰阈值的双阈值判定，准确定位了LOS脉冲开始和结束时刻，使得拟合区间更为精确，利用高斯拟合能够最终获得更为精确的拟合峰中心时刻，提高时间偏差测量精度。
附图说明
图1为本发明原理图。
图2为本发明时间偏差估计方法流程图。
具体实施方式
为使本发明的上述特征和优点更明显易懂，下文特举实施例，并配合附图作详细说明。
采集到的信号为周期性的脉冲信号，如图1所示。其中脉冲的重复频率为1kHz(周期为1ms)，脉冲宽度为1μs，占整个周期的1/1000。利用高精度采集卡进行电压数据的连续采集，本发明方法流程如图2所示，考虑整个过程：
1.在未采集到LOS脉冲信号时，测到的电压值很小，一般为***的噪底(假设5mV)；
2.突然检测到电压值高于脉冲检测阈值Vthreshold时，例如：20mV，此时记录下该时刻tth1，对应某一LOS脉冲信号的开始；
3.当检测到电压值高于脉冲达峰判定阈值Vcheck时，例如：200mV，此时记录下时刻tcheck；
4.当再次检测到电压值低于脉冲检测阈值Vthreshold时，例如：20mV，此时记录下该时刻tth2，对应某一LOS脉冲信号的结束；
5.如果满足tcheck＜tth2，则认为(tth1,tth2)区间对应一个完整的LOS脉冲信号；
6.如果5中的条件不满足，说明tth2仍然处于脉冲的上升沿，实际是脉冲上升沿的电压波动，继续检测，直到电压值再次低于脉冲检测阈值Vthreshold时，记录下该时刻tth3，对应某一LOS脉冲信号的结束，此时(tth1,tth3)区间对应一个完整的LOS脉冲信号；若tcheck小于tth3，则进行下一步操作；
7.根据具体实验，在拟合过程中具体要多选择一些数据来进行拟合。实际拟合区间选择为：(tth1-nts,tth2+nts)。如果出现6中的情况，这一区间为(tth1-nts,tth3+nts)。其中ts表示单次采样的采样时间,n为拟合点数。
8.在拟合区间内进行高斯拟合，获得拟合后LOS脉冲信号峰值对应的中心时刻。如前所述，第三台光梳分别与两台待测光梳进行线性光学采样后得到了两列LOS信号，分别通过上述拟合找出相应的两个中心时刻，两者相减得出放大后的时间间隔，然后除以放大倍数得到精确的两台待测光梳的时间偏差。
以上实施例仅用以说明本发明的技术方案而非对其进行限制，本领域的普通技术人员可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明的精神和范围，本发明的保护范围应以权利要求所述为准。
Claims (8)
1.一种线性光学采样的双阈值拟合方法，其步骤包括：
1)对于待拟合的LOS信号进行检测；其中，参考光梳的频率与待测光梳的频率fr相差设定量Δfr；将参考光梳与待测光梳进行线性光学采样得到所述LOS信号；
2)当检测到LOS信号电压值高于脉冲检测阈值Vthreshold时，记录下当前时刻tth1作为LOS信号的开始；
3)当检测到LOS信号电压值高于脉冲达峰判定阈值Vcheck时，记录下当前时刻tcheck；当再次检测到LOS信号电压值低于脉冲检测阈值Vthreshold时，记录下当前时刻tth2；
4)如果满足tcheck＜tth2，则认为(tth1,tth2)区间对应一个完整的LOS脉冲信号；如果不满足tcheck＜tth2，则继续检测，直到LOS信号电压值再次低于脉冲检测阈值Vthreshold时，记录下当前时刻tth3，此时(tth1,tth3)区间对应一个完整的LOS脉冲信号；
5)根据完整的LOS脉冲信号对应的时间区间拟合得到LOS信号峰值对应的中心时刻t。
2.根据权利要求1所述的方法，其特征在于，获取连续n个LOS脉冲信号对应的时间区间(tth1-nts,tth2+nts)；步骤5)中，根据(tth1-nts,tth2+nts)拟合得到LOS信号峰值对应的中心时刻t；其中，ts表示单次采样的采样时间，n为拟合点数。
3.根据权利要求1所述的方法，其特征在于，LOS信号中的噪底电压为5mV，脉冲检测阈值Vthreshold为20mV，脉冲达峰判定阈值Vcheck为200mV。
4.根据权利要求1或2或3所述的方法，其特征在于，利用高斯拟合得到LOS信号峰值对应的中心时刻t。
5.一种线性光学采样的双阈值时间偏差估计方法，其步骤包括:
1)对于待拟合的两列LOS信号分别进行检测；其中，将参考光梳与第一待测光梳进行线性光学采样所得LOS信号记为第一列LOS信号，将参考光梳与第二待测光梳进行线性光学采样所得LOS信号记为第二列LOS信号；参考光梳的频率与待测光梳的频率fr相差设定量Δfr；
2)当检测到第一列LOS信号电压值高于脉冲检测阈值Vthreshold时，记录下当前时刻tth1作为第一列LOS信号的开始；
3)当检测到第一列LOS信号电压值高于脉冲达峰判定阈值Vcheck时，记录下当前时刻tcheck；当再次检测到第一列LOS信号电压值低于脉冲检测阈值Vthreshold时，记录下当前时刻tth2；
4)如果满足tcheck＜tth2，则认为(tth1,tth2)区间对应一个完整的LOS脉冲信号；如果不满足tcheck＜tth2，则继续检测，直到第一列LOS信号电压值再次低于脉冲检测阈值Vthreshold时，记录下当前时刻tth3，此时(tth1,tth3)区间对应一个完整的LOS脉冲信号；
5)根据完整的LOS脉冲信号对应的时间区间拟合得到第一列LOS信号峰值对应的中心时刻t1；
6)按照步骤2～5)的方法拟合得到第二列LOS信号峰值对应的中心时刻t2；
7)根据中心时刻t1、t2的差值得出放大后的时间间隔，将所述时间间隔除以对应的放大倍数得到两台待测光梳的时间偏差。
6.根据权利要求5所述的方法，其特征在于，获取连续n个LOS脉冲信号对应的时间区间(tth1-nts,tth2+nts)；步骤5)中，根据(tth1-nts,tth2+nts)拟合得到LOS信号峰值对应的中心时刻t；其中，ts表示单次采样的采样时间，n为拟合点数。
7.根据权利要求5所述的方法，其特征在于，LOS信号中的噪底电压为5mV，脉冲检测阈值Vthreshold为20mV，脉冲达峰判定阈值Vcheck为200mV。
8.根据权利要求5或6或7所述的方法，其特征在于，利用高斯拟合得到LOS信号峰值对应的中心时刻。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210826088.7A CN115356906B (zh) | 2022-07-13 | 2022-07-13 | 一种线性光学采样的双阈值拟合方法及时间偏差估计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210826088.7A CN115356906B (zh) | 2022-07-13 | 2022-07-13 | 一种线性光学采样的双阈值拟合方法及时间偏差估计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115356906A CN115356906A (zh) | 2022-11-18 |
| CN115356906B true CN115356906B (zh) | 2024-04-30 |
Family
ID=84031927
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210826088.7A Active CN115356906B (zh) | 2022-07-13 | 2022-07-13 | 一种线性光学采样的双阈值拟合方法及时间偏差估计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115356906B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5159491A (zh) * | 1974-11-20 | 1976-05-24 | Hewlett Packard Yokogawa | |
| CN1270391A (zh) * | 1999-04-09 | 2000-10-18 | 三洋电机株式会社 | 信号记录设备 |
| JP2001267894A (ja) * | 2000-03-14 | 2001-09-28 | Nippon Telegr & Teleph Corp <Ntt> | 電圧比較回路 |
| US6377082B1 (en) * | 2000-08-17 | 2002-04-23 | Agere Systems Guardian Corp. | Loss-of-signal detector for clock/data recovery circuits |
| US6489803B1 (en) * | 2001-07-03 | 2002-12-03 | Silicon Laboratories, Inc. | Use of dual hysteresis modes in determining a loss of signal output indication |
| CN106253977A (zh) * | 2016-08-22 | 2016-12-21 | 青岛海信宽带多媒体技术有限公司 | Los告警判决门限的调整方法及光模块 |
| WO2018086949A1 (de) * | 2016-11-11 | 2018-05-17 | Leoni Kabel Gmbh | Verfahren und messanordnung zur überwachung einer leitung |
| CN108873008A (zh) * | 2018-06-12 | 2018-11-23 | 天津大学 | 一种基于双光梳干涉的水下距离高精度测量方法 |
| CN113534106A (zh) * | 2021-06-03 | 2021-10-22 | 北京大学 | 一种微腔光梳激光器、测距装置及测距方法 |
| CN114513251A (zh) * | 2022-01-04 | 2022-05-17 | 北京大学 | 一种基于双混频的连续相位测量方法 |
-
2022
- 2022-07-13 CN CN202210826088.7A patent/CN115356906B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5159491A (zh) * | 1974-11-20 | 1976-05-24 | Hewlett Packard Yokogawa | |
| CN1270391A (zh) * | 1999-04-09 | 2000-10-18 | 三洋电机株式会社 | 信号记录设备 |
| JP2001267894A (ja) * | 2000-03-14 | 2001-09-28 | Nippon Telegr & Teleph Corp <Ntt> | 電圧比較回路 |
| US6377082B1 (en) * | 2000-08-17 | 2002-04-23 | Agere Systems Guardian Corp. | Loss-of-signal detector for clock/data recovery circuits |
| US6489803B1 (en) * | 2001-07-03 | 2002-12-03 | Silicon Laboratories, Inc. | Use of dual hysteresis modes in determining a loss of signal output indication |
| CN106253977A (zh) * | 2016-08-22 | 2016-12-21 | 青岛海信宽带多媒体技术有限公司 | Los告警判决门限的调整方法及光模块 |
| WO2018086949A1 (de) * | 2016-11-11 | 2018-05-17 | Leoni Kabel Gmbh | Verfahren und messanordnung zur überwachung einer leitung |
| CN108873008A (zh) * | 2018-06-12 | 2018-11-23 | 天津大学 | 一种基于双光梳干涉的水下距离高精度测量方法 |
| CN113534106A (zh) * | 2021-06-03 | 2021-10-22 | 北京大学 | 一种微腔光梳激光器、测距装置及测距方法 |
| CN114513251A (zh) * | 2022-01-04 | 2022-05-17 | 北京大学 | 一种基于双混频的连续相位测量方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于CNN-BiLSTM的NLOS/LOS信号识别模型;李建霖等;《第十三届卫星导航年会论文集》;20220525;全文 * |
| 量子干涉及其在光存储中的应用;郭弘;《科技导报》;20061231;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115356906A (zh) | 2022-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101086510A (zh) | 改进型数字触发器 | |
| CN104545870A (zh) | 用于心率检测的床垫和心率检测方法 | |
| CN101504431A (zh) | 随机脉冲时间序列的纳秒级在线检测*** | |
| CN106443203A (zh) | 一种脉冲信号检测***及方法 | |
| CN106353788B (zh) | 模拟数字混合式脉冲幅度分析器及其分析方法 | |
| CN108896879B (zh) | 基于局部放电信号特征的诊断图谱相位开窗参数调整方法 | |
| CN104502684A (zh) | 一种全数字化峰值到达时刻鉴别方法 | |
| CN107544072B (zh) | 一种预置波形匹配的高精度测距***及方法 | |
| CN106855633B (zh) | 一种同步提取惯性卫星组合量测数据的方法 | |
| CN110596439B (zh) | 一种数字触发检测方法 | |
| CN110554364B (zh) | 一种雷达信号脉冲到达时间测量精度的测试***及方法 | |
| CN103368676A (zh) | 一种基于周期性脉冲信号的异地数据同步采集方法及*** | |
| CN101136624B (zh) | 脉冲信号识别装置及方法 | |
| GB2590556A (en) | Systems and methids of QT interval analysis | |
| CN106556780A (zh) | 一种局部放电类型确定方法及*** | |
| CN115356906B (zh) | 一种线性光学采样的双阈值拟合方法及时间偏差估计方法 | |
| CN107976480A (zh) | 一种飞行时间质谱仪的高精度时间测量方法 | |
| CN103698779A (zh) | 卫星授时机及导航接收机测试站 | |
| CN105403765A (zh) | 一种用于脉冲功率测量的稳幅环路装置 | |
| CN107525686A (zh) | 一种航天器多遥测参数分析窗口划分方法 | |
| CN114236594A (zh) | 一种核脉冲信号数字三角-梯形双通道成形方法 | |
| CN212254401U (zh) | 一种等效采样的光纤分布式温度测量装置 | |
| CN104977978B (zh) | 一种用于***时钟不同步信号的采集及处理方法 | |
| CN109633243B (zh) | 一种基于多相位采样的束流信号峰值幅度精确提取方法 | |
| CN101661010A (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 |