CN109270831A - Bpm短波多频点定时系统 - Google Patents
Bpm短波多频点定时系统 Download PDFInfo
- Publication number
- CN109270831A CN109270831A CN201811387186.5A CN201811387186A CN109270831A CN 109270831 A CN109270831 A CN 109270831A CN 201811387186 A CN201811387186 A CN 201811387186A CN 109270831 A CN109270831 A CN 109270831A
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- module
- utc
- signals
- 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
Links
- 238000001914 filtration Methods 0.000 claims abstract description 19
- 238000005070 sampling Methods 0.000 claims description 27
- 238000001514 detection method Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 15
- 230000003044 adaptive effect Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 9
- 230000002596 correlated effect Effects 0.000 claims description 7
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000000875 corresponding effect Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 230000011664 signaling Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000008571 general function Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 239000005433 ionosphere Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241001672694 Citrus reticulata Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/14—Setting the time according to the time information carried or implied by the radio signal the radio signal being a telecommunication standard signal, e.g. GSM
- G04R20/16—Tuning or receiving; Circuits therefor
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/14—Setting the time according to the time information carried or implied by the radio signal the radio signal being a telecommunication standard signal, e.g. GSM
- G04R20/18—Decoding time data; Circuits therefor
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R40/00—Correcting the clock frequency
- G04R40/06—Correcting the clock frequency by computing the time value implied by the radio signal
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Synchronisation In Digital Transmission Systems (AREA)
- Electric Clocks (AREA)
Abstract
本发明提供了一种BPM短波多频点定时系统，通道选择单元接收BPM授时台发送的四个频点授时信号，将各频点授时信号均变为固定载频点2.5MHz数字信号；解调搜索单元对接收到的数字信号进行解调和低通滤波，将得到的正交时号信号与本地1KHz正交时号信号进行四路折叠匹配相关，利用相关值和门限值对两个相关峰进行检测，得到两个相关峰最大值的位置送入定时输出模块；定时输出单元识别出UTC秒时号和分时号，将UTC秒时号进行滤波处理防止误判，对本地秒信号、分信号进行同步。本发明采用全数字接收方式，在任何地点均可接收BPM短波授时信号，具有多频点接收、自动化程度高、定时精度高、全数字化等特点。
Description
技术领域
本发明涉及一种定时系统，属于无线电授时和定时技术领域。
背景技术
BPM短波授时系统作为大科学装置之一，由中国科学院国家授时中心承担发播任务，通过电离层一次或多次反射进行时间信号传递，可覆盖我国陆地和近海海域，授时精度为1ms，具有覆盖范围广、发射设备简单、接收成本低等优点。其主要功能是进行定时和校频，也应用于通信、无线电波传播、气象和电离层等相关科学研究。
BPM短波授时系统每天以四个频率(2.5MHz±5KHz、5MHz±5KHz、10MHz±5KHz、15MHz±10KHz)交替发播标准时间和标准频率信号，通过调幅方式发播协调世界时UTC时号、世界时UT1时号、无调制载波及BPM呼号。四个频率的发播时间如表1所示，各频点发播程序如表2所示。
表1BPM授时系统发播时间
| 发射频率/MHz | UTC时间 | 北京时间 |
| 2.5 | 07:30—01：00 | 15:30—09：00 |
| 5.0 | 00:00—24：00 | 00:00—24：00 |
| 10.0 | 00:00—24：00 | 00:00—24：00 |
| 15.0 | 01:00—09：00 | 09:00—17：00 |
表2BPM授时系统各频点发播程序
| 1)59<sup>m</sup>00<sup>s</sup>～00<sup>m</sup>00<sup>s</sup> 6)29<sup>m</sup>00<sup>s</sup>～30<sup>m</sup>00<sup>s</sup> | BPM呼号(1分钟) |
| 2)00<sup>m</sup>00<sup>s</sup>～10<sup>m</sup>00<sup>s</sup> 7)30<sup>m</sup>00<sup>s</sup>～40<sup>m</sup>00<sup>s</sup> | UTC时号(10分钟) |
| 3)10<sup>m</sup>00<sup>s</sup>～15<sup>m</sup>00<sup>s</sup> 8)40<sup>m</sup>00<sup>s</sup>～45<sup>m</sup>00<sup>s</sup> | 无调制载波(5分钟) |
| 4)15<sup>m</sup>00<sup>s</sup>～25<sup>m</sup>00<sup>s</sup> 9)45<sup>m</sup>00<sup>s</sup>～55<sup>m</sup>00<sup>s</sup> | UTC时号(10分钟) |
| 5)25<sup>m</sup>00<sup>s</sup>～29<sup>m</sup>00<sup>s</sup> 7)55<sup>m</sup>00<sup>s</sup>～59<sup>m</sup>00<sup>s</sup> | UT1时号(4分钟) |
UTC和UT1时号均采用1KHz的标准音频调制，起始时刻为正弦波的零相位，UTC时号包含10ms的秒时号和300ms的分时号，UT1包含100ms的秒时号和300ms的分时号，具体格式如下所示：
式中，m(t)为发播时号信号，f0频率1KHz；
BPM呼号在每小时的59m00s～00m00s和29m00s～30m00s发播，前40s为莫尔斯电码：—····――·——，后20s为女声普通话语音广播：BPM标准时间标准频率发波台。同时为避免与周边国家短波授时台(如日本的JTY和印度的ATA等授时台)的相互干扰，BPM的UTC时号超前UTC(NTSC)20ms发播。
目前，基于全球卫星导航系统的定时终端使用比较广泛，但卫星导航系统存在信号功率低，容易受到遮挡和干扰等问题，因此短波授时仍作为重要授时手段使用。短波主要靠天波进行传播，传播信道属于变参信号，信号很容易受到电离层变化和其他无线电系统干扰的影响，且时号不是连续发播，接收信号时有时无。现在我国大部分用户使用的短波接收机依然是模拟技术，多数采用传统的超外差调谐的方式和分立元件构成，例如1987年范荣美、吕彩田、李为丰在论文中提出的《BPM定时仪研究》及国家授时中心研制的PO23BPM短波定时仪，这些设备存在结构复杂、精度低、数字化程度低、需要手动选择接收频点等问题。1988年范荣美在专利“中点检测技术及短波时号定时仪”中提出利用中点检测技术进行时号检测，其要求接收时号是对称的，且检测效率不高。随着数字化技术的发展和应用，出现了一些数字化短波接收手段，例如西安电子科技大学马卓、张伟、刘维、杜栓义在专利提出的“利用短波授时信号校准本地时钟的方法”利用傅里叶变化技术实现校准本地时钟，但需要根据发播程序确定接收信号时间。总结现有BPM定时终端，不同程度上存在自动化程度低、时号检测效率低、数字化程度不高、需要手动选择接收频点和调节信号增益等问题。因此研制多频点接收，全数字化接收，自动化程度高的BPM定时接收机，对我国BPM授时系统的建设和发展以及研制新型BPM定时接收机有着非常重要的意义。
发明内容
为了克服现有技术的不足，本发明提供一种BPM多频点定时系统，采用全数字接收方式，在任何地点均可接收BPM短波授时信号，输出标准的UTC秒信号和分信号，具有多频点接收、自动化程度高、定时精度高、全数字化等特点，可满足我国BPM短波用户对BPM定时系统低成本和全自动的需求。
本发明解决其技术问题所采用的技术方案是：一种BPM短波多频点定时系统，包括通道选择单元、解调搜索单元和定时输出单元。
所述的通道选择单元接收BPM授时台发送的四个频点授时信号，实时计算出其中信号功率最大的一个频点授时信号，通过AD采样将各频点授时信号均变为固定载频点2.5MHz数字信号，然后经过带通滤波器送入解调搜索单元；所述的解调搜索单元对接收到的数字信号进行解调和低通滤波，得到两路1KHz的正交时号信号，将得到的正交时号信号与本地1KHz正交时号信号进行四路折叠匹配相关，相关结果求取平方和，经过相关延迟相减得到两个相关值峰，然后通过短时自适应门限方法计算出门限值，最后利用相关值和门限值对两个相关峰进行检测，得到两个相关峰最大值的位置送入定时输出模块；所述的定时输出单元首先通过信号相关峰的出现的位置和理论位置关系对信号进行识别，识别出UTC秒时号和分时号，接着将UTC秒时号进行滤波处理防止误判，滤波后结果进行时延计算得到UTC秒时号准确位置，通过UTC秒时号初始位置对本地秒信号进行同步，同时通过秒信号计数和识别的分时号对本地分信号进行同步。
所述的通道选择单元包括时钟基准模块、采样时钟产生模块、AD采样模块、FIR带通滤波器模块、信号强度检测模块和通道选择模块；首先通道选择模块自动按顺序从四个频点中设置接收频点，给出通道号；时钟基准模块接收来自外部的时钟源，通过锁相环产生系统所需的工作时钟；采样时钟产生模块根据通道选择模块设定的通道号通过工作时钟分频产生四个频点的采样时钟，采样后各频点信号均为2.5MHz；AD采样模块通过采样时钟将所选频点的模拟信号转换成数字信号，完成对信号的采样、量化和编码；FIR带通滤波器模块对AD采样信号进行滤波，滤除带外噪声和干扰，并将滤波结果送入信号强度检测模块；信号强度检测模块通过快速傅里叶变化对接收信号进行计算，计算出对应通道号2.5MHz频点信号功率；最后，通道选择模块根据每一个频点信号功率判断四个频点中信号功率最大的一个作为系统接收信号；通道选择单元按照设定周期重复计算功率最大的信号。
所述解调搜索单元包括解调模块、低通滤波器模块、相关延迟相减模块、自适应门限模块和自动搜索模块；解调模块使用本地ROM表中存储的一周2.5MHz正交载波与FIR滤波器输出信号进行混频，得到两路正交混频信号；低通滤波模块分别对两路正交混频信号进行滤波，滤除高频分量和噪声，得到两路1KHz正交时号信号，并通过累加平均的方式进行降采样；相关延迟相减模块采用四路折叠匹配相关器将两路1KHz正交信号与本地ROM表中存储的10ms正交1KHz信号分别进行匹配相关得到四个相关值，匹配相关时长为10ms，将四个相关值平方和相加后得到相关值R(t)，相关值R(t)延迟10ms得到延迟相关值R(t+10ms)，通过相减R(t)-R(t+10ms)得到相关延迟相减结果R′(t)；自适应门限模块按照设定周期对相关延迟相减结果R′(t)累加求平均计算出信号检测的门限值，平均结果乘以8得到动态自适应门限V；自动搜索模块在时钟的驱动下计数作为搜索基准计数，通过比较R′(t)和门限V，连续判断一个搜索基准计数内相关延迟相减结果R′(t)中最大的两个极值点，记录两个极值点在一个搜索基准计数内计数的位置N1和N2，并在一个搜索基准计数结束时将位置值N1和N2送入定时输出单元，然后检测下一个基准计数内最大极值点位置。
所述定时输出单元包括时号识别模块、峰值滤波模块、时延计算模块、UTC秒信号产生模块和UTC分信号产生模块；时号识别模块接收解调搜索单元送来的位置N1和N2，通过判断N1和N2的位置差ΔN＝N2-N1对时号进行识别，若ΔN为f*10ms，则判定为UTC秒时号，N1为UTC秒时号相关位置；若ΔN为f*990ms，也判定为UTC秒时号，索引值N2为UTC秒时号相关位置，满足以上两种条件时产生UTC秒时号标识和索引值N(N为N1或N2)；若ΔN为f*300ms或f*700ms，则判定为分时号，产生分信号标识；若ΔN为f*100ms或f*900ms，则判定为UT1秒时号；时号识别模块中连续5次识别出UTC秒时号，且前后一次N的值误差在±10内，则认为UTC秒时号识别正确；将识别的UTC秒时号位置N送入峰值滤波模块中；峰值滤波模块将时号识别模块中UTC秒时号相关位置N送入存储长度为20的移位寄存器，然后进行从大到小排序，取中间两个值的平均值作为秒时号相关峰位置N′，从而完成UTC秒时号位置的滤波和误判值剔除，得到秒时号相关峰位置；当1分钟内未成功识别出UTC秒时号，清除移位寄存器，重新开始峰值滤波模块；时延计算模块计算接收系统定时信号相对于发射台的总定时时延Δt＝tr-t20ms+td，式中tr为接收终端处理时延，定时系统通过接收BPM授时模拟器信号后输出的1PPS信号与模拟源输出的1PPS信号的时差测量均值获取；t20ms为发播UTC秒时号超前UTC(NTSC)20ms；td为传播时延；υd为电磁波传输速度；将秒信号时延换算为秒信号调整量NΔt＝Δt×f，UTC秒时号出现的位置为N′-NΔt；UTC秒信号产生模块是由100KHz时钟产生一个1s的计数，秒信号起始位置为自动搜索模块中UTC秒时号在1s基准计数中的位置N′-NΔt，实现了UTC秒信号的同步；UTC分信号产生模块在UTC秒信号的触发下进行秒计数，计数范围为0～59，计数为0时输出UTC分信号，当时号识别模块中产生分信号标识，则UTC秒信号和分信号重合，将秒计数置为0，实现分信号的同步。
本发明的有益效果是：
(1)本发明可接收BPM短波台发出的四个频点，通过自适应选频接收信号强度最好的一个频点信号，无需手动选择频点；
(2)本发明采用相关延迟相减方法获取信号相关峰，信号识别准确率高，无需通过发播时间判定此刻是否接收信号；
(3)本发明采用短时自适应门限技术，无需调节接收信号增益；
(4)本发明对时号峰值检测结果进行滤波获得秒信号位置，具有较高的定时精度；
(5)本发明采用全数字解调和处理方式，在一般FPGA芯片上就可实现，集成化程度高可节省硬件电路资源，可靠性较高，能够，满足BPM短波用户的精度要求。
附图说明
图1是本发明总体功能框架图
图2是BPM短波授时系统时号波形示意图
图3是本发明信号处理流程框图
图4是一路折叠匹配相关器结构图
图5是时号经过四路折叠匹配相关平方相加后仿真图
图6是时号相关延迟相减结果仿真图
图7是BPM定时系统总定时延迟示意图
图8是在5MHz频点上实测的时差图
具体实施方式
下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。
本发明提供一种BPM短波多频点定时系统，用于BPM短波授时信号的数字化接收和处理，为用户提供标准的UTC秒信号和分信号。
BPM多频点定时系统，包括通道选择单元，解调搜索单元和定时输出单元。通道选择单元是接收BPM授时台发送的四个频点授时信号，实时计算出其中信号最好的一个频点授时信号，通过AD采样各频点均变为固定载频点2.5MHz数字信号，然后经过带通滤波，将其送入解调搜索单元；解调搜索单元对接收到的数字信号进行解调和低通滤波，得到两路1KHz的正交时号信号，接着与本地1KHz正交时号信号进行四路折叠匹配相关，相关结果求取平方和，在经过相关延迟相减得到两个相关值峰，然后通过短时自适应门限方法计算出门限值，最后利用相关值和门限值对两个相关峰进行检测，得到两个相关峰最大值的位置送入定时输出模块；定时输出单元首先通过信号相关峰的出现的位置和理论位置关系对信号进行识别，识别出UTC秒时号和分时号，接着将UTC秒时号进行滤波处理防止误判，滤波后结果进行时延计算得到UTC秒时号准确位置，通过UTC秒时号初始位置对本地秒信号进行同步，同时通过秒信号计数和识别的分时号对本地分信号进行同步。
所述的通道选择单元包括时钟基准模块、采样时钟产生模块、AD采样模块、FIR带通滤波器模块、信号强度检测模块和通道选择模块。首先通道选择模块自动按顺序从四个频点(2.5MHz、5MHz、10Mhz、15MHz)中设置接收频点，给出通道号；时钟基准模块接收来自外部的时钟源，通过锁相环产生系统所需的工作时钟；采样时钟产生模块根据通道选择模块设定的通道号通过工作时钟分频产生四个频点(2.5MHz、5MHz、10Mhz、15MHz)的采样时钟(12.5MHz、7.5MHz、12.5MHz、12.5MHz)，采样后各频点信号均为2.5MHz；AD采样模块通过采样时钟将所选频点的模拟信号转换成数字信号，完成对信号的采样、量化和编码；FIR带通滤波器模块对AD采样信号进行滤波，滤除带外噪声和干扰，FIR滤波器带宽一般取8KHz，并将滤波结果送入信号强度检测模块。信号强度检测模块通过快速傅里叶变化(2048点)对接收信号进行计算，计算出对应通道号2.5MHz频点信号功率；最后，通道选择模块根据每一个频点信号功率判断四个频点中信号功率最大的一个作为系统接收信号。通道选择单元根据信号发播格式对上述步骤每5分钟进行一次。
所述解调搜索单元包括解调模块、低通滤波器模块、相关延迟相减模块、自适应门限模块和自动搜索模块。解调搜索单元是利用通道选择单元中FIR滤波器输出的信号，频点为2.5MHz，工作时钟为12.5MHz。首先解调模块使用本地ROM表中存储的一周2.5MHz正交载波与FIR滤波器输出信号进行混频，得到两路正交混频信号；低通滤波模块分别对两路正交混频信号进行滤波，滤除高频分量和噪声，得到两路1KHz正交时号信号，带宽为±3kHz，并通过累加平均的方式进行降采样率至100KHz(累加次数为125次)；相关延迟相减模块采用四路折叠匹配相关器将两路1KHz正交信号与本地ROM表中存储的10ms正交1KHz信号分别进行匹配相关得到四个相关值，匹配相关时长为10ms，将四个相关值平方和相加后得到相关值R(t)，相关值R(t)延迟10ms得到延迟相关值R(t+10ms)，通过相减R(t)-R(t+10ms)得到相关延迟相减结果R′(t)。自适应门限模块通过延迟相减结果R′(t)计算出信号检测的门限值，计算方法是每2秒对相关延迟相减结果R′(t)累加求平均，平均结果乘以8得到动态自适应门限V，此过程连续重复执行；自动搜索模块在100KHz时钟的驱动下，循环产生时长为1s的计数，作为搜索基准计数，通过比较R′(t)和门限V，连续判断1s基准计数内相关延迟相减结果R′(t)中最大的两个极值点，记录两个极值点在1s内计数的位置N1和N2，并在1s计数结束时将位置值N1和N2送入定时输出单元，然后检测下一个1s基准计数内最大极值点位置。
所述定时输出单元包括时号识别模块、峰值滤波模块、时延计算模块、UTC秒信号产生模块和UTC分信号产生模块。时号识别模块接收解调搜索单元送来的位置N1和N2，通过判断N1和N2的位置差ΔN＝N2-N1对时号进行识别，判断过程如下描述：
(1)UTC秒时号判断：若ΔN为f*10ms，则判定为UTC秒时号，N1为UTC秒时号相关位置；若ΔN为f*990ms，也判定为UTC秒时号，索引值N2为UTC秒时号相关位置，满足以上两种条件时产生UTC秒时号标识和索引值N(N为N1或N2)；
(2)UTC/UT1分时号判断：若ΔN为f*300ms或f*700ms，则判定为分时号，产生分信号标识。
(3)UT1秒时号判断：若ΔN为f*100ms或f*900ms，则判定为UT1秒时号；
时号识别模块中连续5次识别出UTC秒时号，且前后一次N的值误差在±10内，则认为UTC秒时号识别正确。将识别的UTC秒时号位置N送入峰值滤波模块中。峰值滤波模块将时号识别模块中UTC秒时号相关位置N送入存储长度为20的移位寄存器，然后进行从大到小排序，取中间两个值的平均值作为秒时号相关峰位置N′，从而完成UTC秒时号位置的滤波和误判值剔除，得到较为准确的秒时号相关峰位置。因为短波发播不连续且信号传播不稳定，可能长时间搜索不到信号，因此当1分钟内未成功识别出UTC秒时号，应清除移位寄存器，重新开始峰值滤波模块；时延计算模块计算接收系统定时信号相对于发射台的总定时时延，其计算公式为：
Δt＝tr-t20ms+td
式中tr为接收终端处理时延，定时系统可通过接收BPM授时模拟器信号后输出的1PPS信号与模拟源输出的1PPS信号的时差测量均值获取；t20ms为发播UTC秒时号超前UTC(NTSC)20ms，为固定值；td为传播时延，可采用经验估算公式td＝D/υd进行估算，其中D为收发两地间的大地基线长度，υd为电磁波传输速度，取285000km/s。时号相关处理时钟为100kHz，秒时号分辨率为1/f＝10us，因此需要将秒信号时延换算为秒信号调整量NΔt＝Δt×f，因此UTC秒时号出现的位置为N′-NΔt。
UTC秒信号产生模块是由100KHz时钟产生一个1s的计数，秒信号起始位置为自动搜索模块中UTC秒时号在1s基准计数中的位置N′-NΔt，这样就实现了UTC秒信号的同步；UTC分信号产生模块在UTC秒信号的触发下进行秒计数，计数范围为0～59，计数为0时输出UTC分信号，当时号识别模块中产生分信号标识，则UTC秒信号和分信号重合，将秒计数置为0，实现分信号的同步。
参照图1总体功能框架图，本发明公开了BPM短波多频点定时系统，主要包括通道选择单元，解调搜索单元和定时输出单元。系统计算出四个频点信号的强度，选择出最强的短波信号进行解调，采用相关延迟相减和自适应门限技术对语音信号进行搜索捕获，通过判断信号相关峰值的出现的位置进行时号识别，对识别位置结果进行滤波，并进行时延修正，从而实现UTC秒信号和分信号的同步。依据图2BPM短波时号示意图和图3信号处理流程框，实现方案包括如下：
步骤1：信号采集，时钟基准产生整个系统工作时钟12.5MHz，根据短波信号四个频点(2.5MHz、5MHz、10MHz、15MHz)使用进行数字采集，分别产生其采样时钟为12.5MHz、7.5MHz、12.5MHz和12.5MHz，采样后信号频率统一为2.5Mhz。
步骤2：通道选择，通道选择器按照接收顺序设置对应频点的通道号，并产生该频点信号的采样频率，通过AD芯片对该频点信号进行采集，将采样后的数字信号通过FIR带通滤波器，滤除带外噪声和干扰，滤波器采用128阶的hamming窗，带宽为±3kHz。然后将FIR输出信号进行快速傅里叶变化(FFT)，FFT时钟速率为12.5MHz，点数为2048点；通过FFT计算出信号功率谱，检测每个通道2.5MHz频点上信号强度，通过比较选择信号功率最大的一个作为接收频点。根据短波信号发播程序，该步骤每5分钟重复一次。
步骤3：信号解调，本地ROM表存储采样时钟为12.5MHz，频率为2.5MHz，长度为1个周期的标准数字正交单载波信号，在DDS驱动下循环产生本地ROM读取地址，生成正交支路的载波。其中DDS工作时钟为12.5MHz，累加器字长为2^32，频率控制字为858993459。正交支路载波分别与步骤2得到的2.5MHz信号相乘，通过低通滤波器，得到两路正交语音调制信号。其中低通滤波器采用128阶的Taylor窗，工作时钟为100kHz，带宽2kHz。最后通过累加平均的方式将两路正交语音信号降采样至100KHz(累加次数为125次)。
步骤4：时号相关，利用步骤3获取的两路正交语音信号，与本地存储的10ms正交语音信号进行四路折叠匹配滤波，设计折叠匹配滤波器如图4所示，数据通道折叠次数为10次，每折中ROM长度为100，本地语音信号数据ROM长度为10ms，在工作时钟的驱动下循环相关，通过四路折叠匹配相关结果取平方和得到相关值，UTC/UT1时号经过匹配相关后波形如图5所示。然后利用相关延迟相减方法获取延迟相减相关结果。
步骤5：信号捕获，利用步骤4中的延迟相减相关结果进行2秒累加求平均，平均结果乘以8得到动态自适应门限；同时通过100KHz时钟产生1s的基准计数，范围为0～99999，检测出1s内两个相关结果的极大值，记录极大值位置N1和N2，此步骤在1s基准计数内进行，计数结束后重新开始检测。
步骤6：时号识别，利用步骤5极大值位置N1和N2进行时号识别，计算ΔN＝N2-N1，结合图6中UTC/UT1分秒时号相关峰出现的位置关系，按照技术方案中时号识别模块中判断条件进行识别时号，给出UTC秒时号和分时号标识，记录UTC秒时号出现位置N(N1或N2)。
步骤7：峰值滤波，峰值滤波通过对步骤6中判定出UTC秒时号标识，对秒时号相关峰值位置送至移位寄存器，存储长度20次，然后按从大到小排序，取中间两个值的平均值作为接收UTC秒时号位置。当1分钟内未识别出UTC秒时号时，应清除移位寄存器，重新开始峰值滤波模块。
步骤8：秒信号调整量计算，图7为定时系统总定时延迟示意图，按照技术方案中时延计算模块所述方法计算出UTC秒时号在1s基准计数中的位置N′-NΔt。
步骤9：UTC秒信号产生，UTC秒信号是用100KHz时钟进行循环计数产生的。首先根据步骤8计算的UTC秒时号位置，在1s基准计数的N′-NΔt，调整本地秒计数为0，对本地秒信号进行校准；当未检测到UTC秒时号时通过计数器维持秒信号产生。图8是在临潼接收5MHz频点信号，秒信号与UTC(NTSC)秒信号比较实测的时差结果。
步骤10：UTC分信号产生，UTC分信号使用UTC秒信号累加计数(0-59)产生，当步骤6中识别出分时号，清除本地秒触发计数为0，从而对本地UTC分信号进行校准，当没有分信号标识时通过秒信号累加60次进行维持产生。
Claims (4)
1.一种BPM短波多频点定时系统，包括通道选择单元、解调搜索单元和定时输出单元，其特征在于：所述的通道选择单元接收BPM授时台发送的四个频点授时信号，实时计算出其中信号功率最大的一个频点授时信号，通过AD采样将各频点授时信号均变为固定载频点2.5MHz数字信号，然后经过带通滤波器送入解调搜索单元；所述的解调搜索单元对接收到的数字信号进行解调和低通滤波，得到两路1KHz的正交时号信号，将得到的正交时号信号与本地1KHz正交时号信号进行四路折叠匹配相关，相关结果求取平方和，经过相关延迟相减得到两个相关值峰，然后通过短时自适应门限方法计算出门限值，最后利用相关值和门限值对两个相关峰进行检测，得到两个相关峰最大值的位置送入定时输出模块；所述的定时输出单元首先通过信号相关峰的出现的位置和理论位置关系对信号进行识别，识别出UTC秒时号和分时号，接着将UTC秒时号进行滤波处理防止误判，滤波后结果进行时延计算得到UTC秒时号准确位置，通过UTC秒时号初始位置对本地秒信号进行同步，同时通过秒信号计数和识别的分时号对本地分信号进行同步。
2.根据权利要求1所述的BPM短波多频点定时系统，其特征在于：所述的通道选择单元包括时钟基准模块、采样时钟产生模块、AD采样模块、FIR带通滤波器模块、信号强度检测模块和通道选择模块；首先通道选择模块自动按顺序从四个频点中设置接收频点，给出通道号；时钟基准模块接收来自外部的时钟源，通过锁相环产生系统所需的工作时钟；采样时钟产生模块根据通道选择模块设定的通道号通过工作时钟分频产生四个频点的采样时钟，采样后各频点信号均为2.5MHz；AD采样模块通过采样时钟将所选频点的模拟信号转换成数字信号，完成对信号的采样、量化和编码；FIR带通滤波器模块对AD采样信号进行滤波，滤除带外噪声和干扰，并将滤波结果送入信号强度检测模块；信号强度检测模块通过快速傅里叶变化对接收信号进行计算，计算出对应通道号2.5MHz频点信号功率；最后，通道选择模块根据每一个频点信号功率判断四个频点中信号功率最大的一个作为系统接收信号；通道选择单元按照设定周期重复计算功率最大的信号。
3.根据权利要求1所述的BPM短波多频点定时系统，其特征在于：所述解调搜索单元包括解调模块、低通滤波器模块、相关延迟相减模块、自适应门限模块和自动搜索模块；解调模块使用本地ROM表中存储的一周2.5MHz正交载波与FIR滤波器输出信号进行混频，得到两路正交混频信号；低通滤波模块分别对两路正交混频信号进行滤波，滤除高频分量和噪声，得到两路1KHz正交时号信号，并通过累加平均的方式进行降采样；相关延迟相减模块采用四路折叠匹配相关器将两路1KHz正交信号与本地ROM表中存储的10ms正交1KHz信号分别进行匹配相关得到四个相关值，匹配相关时长为10ms，将四个相关值平方和相加后得到相关值R(t)，相关值R(t)延迟10ms得到延迟相关值R(t+10ms)，通过相减R(t)-R(t+10ms)得到相关延迟相减结果R′(t)；自适应门限模块按照设定周期对相关延迟相减结果R′(t)累加求平均计算出信号检测的门限值，平均结果乘以8得到动态自适应门限V；自动搜索模块在时钟的驱动下计数作为搜索基准计数，通过比较R′(t)和门限V，连续判断一个搜索基准计数内相关延迟相减结果R′(t)中最大的两个极值点，记录两个极值点在一个搜索基准计数内计数的位置N1和N2，并在一个搜索基准计数结束时将位置值N1和N2送入定时输出单元，然后检测下一个基准计数内最大极值点位置。
4.根据权利要求1所述的BPM短波多频点定时系统，其特征在于：所述定时输出单元包括时号识别模块、峰值滤波模块、时延计算模块、UTC秒信号产生模块和UTC分信号产生模块；时号识别模块接收解调搜索单元送来的位置N1和N2，通过判断N1和N2的位置差ΔN＝N2-N1对时号进行识别，若ΔN为f*10ms，则判定为UTC秒时号，N1为UTC秒时号相关位置；若ΔN为f*990ms，也判定为UTC秒时号，索引值N2为UTC秒时号相关位置，满足以上两种条件时产生UTC秒时号标识和索引值N(N为N1或N2)；若ΔN为f*300ms或f*700ms，则判定为分时号，产生分信号标识；若ΔN为f*100ms或f*900ms，则判定为UT1秒时号；时号识别模块中连续5次识别出UTC秒时号，且前后一次N的值误差在±10内，则认为UTC秒时号识别正确；将识别的UTC秒时号位置N送入峰值滤波模块中；峰值滤波模块将时号识别模块中UTC秒时号相关位置N送入存储长度为20的移位寄存器，然后进行从大到小排序，取中间两个值的平均值作为秒时号相关峰位置N′，从而完成UTC秒时号位置的滤波和误判值剔除，得到秒时号相关峰位置；当1分钟内未成功识别出UTC秒时号，清除移位寄存器，重新开始峰值滤波模块；时延计算模块计算接收系统定时信号相对于发射台的总定时时延Δt＝tr-t20ms+td，式中tr为接收终端处理时延，定时系统通过接收BPM授时模拟器信号后输出的1PPS信号与模拟源输出的1PPS信号的时差测量均值获取；t20ms为发播UTC秒时号超前UTC(NTSC)20ms；td为传播时延；υd为电磁波传输速度；将秒信号时延换算为秒信号调整量NΔt＝Δt×f，UTC秒时号出现的位置为N′-NΔt；UTC秒信号产生模块是由100KHz时钟产生一个1s的计数，秒信号起始位置为自动搜索模块中UTC秒时号在1s基准计数中的位置N′-NΔt，实现了UTC秒信号的同步；UTC分信号产生模块在UTC秒信号的触发下进行秒计数，计数范围为0～59，计数为0时输出UTC分信号，当时号识别模块中产生分信号标识，则UTC秒信号和分信号重合，将秒计数置为0，实现分信号的同步。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811387186.5A CN109270831B (zh) | 2018-11-21 | 2018-11-21 | Bpm短波多频点定时系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811387186.5A CN109270831B (zh) | 2018-11-21 | 2018-11-21 | Bpm短波多频点定时系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109270831A true CN109270831A (zh) | 2019-01-25 |
| CN109270831B CN109270831B (zh) | 2020-09-22 |
Family
ID=65190922
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811387186.5A Active CN109270831B (zh) | 2018-11-21 | 2018-11-21 | Bpm短波多频点定时系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109270831B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111142366A (zh) * | 2019-12-19 | 2020-05-12 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种新的短波授时方法 |
| CN111638640A (zh) * | 2020-05-27 | 2020-09-08 | 中国科学院国家授时中心 | 一种bpm短波授时信号模拟器 |
| CN112099339A (zh) * | 2020-07-27 | 2020-12-18 | 季华实验室 | Utc数据接收器及计时器 |
| CN112269312A (zh) * | 2020-10-21 | 2021-01-26 | 四川汉星航通科技有限公司 | 一种无线授时信号自动对准装置 |
| CN112612319A (zh) * | 2021-02-01 | 2021-04-06 | 西安空天电子技术有限公司 | 一种bpm短波时号产生器 |
| CN116755319A (zh) * | 2023-08-18 | 2023-09-15 | 青岛鼎信通讯股份有限公司 | 电能表时钟系统及控制方法、装置、计算机设备、介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008149028A2 (fr) * | 2007-05-22 | 2008-12-11 | France Telecom | Synchronisation en phase de noeuds dans un reseau de telecommunications |
| CN105791204A (zh) * | 2016-04-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 一种交流irig-b码的全数字包络检波解调方法与装置 |
| CN107171761A (zh) * | 2017-04-19 | 2017-09-15 | 西安电子科技大学 | 利用短波授时信号校准本地时钟的方法 |
| CN206948336U (zh) * | 2017-06-09 | 2018-01-30 | 西安宏泰时频技术有限公司 | 一种小型bpm短波授时控制系统及控制电路板 |
| CN108333916A (zh) * | 2018-03-15 | 2018-07-27 | 中国科学院国家授时中心 | 一种基于Chirp信号的BPM短波授时信号与定时方法 |
-
2018
- 2018-11-21 CN CN201811387186.5A patent/CN109270831B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008149028A2 (fr) * | 2007-05-22 | 2008-12-11 | France Telecom | Synchronisation en phase de noeuds dans un reseau de telecommunications |
| CN105791204A (zh) * | 2016-04-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 一种交流irig-b码的全数字包络检波解调方法与装置 |
| CN107171761A (zh) * | 2017-04-19 | 2017-09-15 | 西安电子科技大学 | 利用短波授时信号校准本地时钟的方法 |
| CN206948336U (zh) * | 2017-06-09 | 2018-01-30 | 西安宏泰时频技术有限公司 | 一种小型bpm短波授时控制系统及控制电路板 |
| CN108333916A (zh) * | 2018-03-15 | 2018-07-27 | 中国科学院国家授时中心 | 一种基于Chirp信号的BPM短波授时信号与定时方法 |
Non-Patent Citations (1)
| Title |
|---|
| 郭美军等: "BPM授时信号传输损耗的仿真建模", 《宇航计测技术》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111142366A (zh) * | 2019-12-19 | 2020-05-12 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种新的短波授时方法 |
| CN111638640A (zh) * | 2020-05-27 | 2020-09-08 | 中国科学院国家授时中心 | 一种bpm短波授时信号模拟器 |
| CN112099339A (zh) * | 2020-07-27 | 2020-12-18 | 季华实验室 | Utc数据接收器及计时器 |
| CN112099339B (zh) * | 2020-07-27 | 2021-08-24 | 季华实验室 | Utc数据接收器及计时器 |
| CN112269312A (zh) * | 2020-10-21 | 2021-01-26 | 四川汉星航通科技有限公司 | 一种无线授时信号自动对准装置 |
| CN112612319A (zh) * | 2021-02-01 | 2021-04-06 | 西安空天电子技术有限公司 | 一种bpm短波时号产生器 |
| CN112612319B (zh) * | 2021-02-01 | 2024-05-03 | 西安空天电子技术有限公司 | 一种bpm短波时号产生器 |
| CN116755319A (zh) * | 2023-08-18 | 2023-09-15 | 青岛鼎信通讯股份有限公司 | 电能表时钟系统及控制方法、装置、计算机设备、介质 |
| CN116755319B (zh) * | 2023-08-18 | 2023-11-21 | 青岛鼎信通讯股份有限公司 | 电能表时钟系统及控制方法、装置、计算机设备、介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109270831B (zh) | 2020-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109270831A (zh) | Bpm短波多频点定时系统 | |
| CN101666869B (zh) | 微弱卫星导航信号二次捕获方法与装置 | |
| CN107171761B (zh) | 利用短波授时信号校准本地时钟的方法 | |
| CN103033824A (zh) | 高性能导航卫星空间信号质量评估方法 | |
| MXPA05000474A (es) | Metodo y aparato de mitigacion de correlacion cruzada para utilizarse en un receptor del sistema de posicionamiento global. | |
| CN105743612B (zh) | 实时盲解调大频移短时突发信号的方法 | |
| CN103199887A (zh) | 一种适用于直接扩频信号的通用捕获方法 | |
| CN107819534B (zh) | 导航信息广播方法、发射机、接收方法、接收机和系统 | |
| CN105704081B (zh) | 多个偏移载波调幅信号的优选和解调方法 | |
| WO2015013602A2 (en) | Accurate timing distribution by high-frequency radio | |
| CN110988924A (zh) | 一种基于阵列天线抗干扰bds/gps接收机 | |
| CN106656400A (zh) | Pcma系统时延抖动情况下的精确频差快速捕获方法 | |
| CN1009790B (zh) | 短波无线电网络的通信方法和设备 | |
| CN107621643B (zh) | 一种适用于导航信号质量评估的相关域参数精确解算方法 | |
| Chen et al. | Evaluation of binary offset carrier signal capture algorithm for development of the digital health literacy instrument | |
| CN107037457A (zh) | 一种基于Inmarsat系统的星基增强接收机 | |
| CN106802424A (zh) | 一种多频卫星导航接收机快速引导跟踪方法及装置 | |
| CN109633715B (zh) | 一种基于并行码相位捕获gps和北斗三代信号的方法 | |
| CN105204037B (zh) | 一种长码扩频测距信号相关损失测试方法 | |
| CN110149197A (zh) | 一种用于时钟同步系统的高精度同步方法及同步系统 | |
| Berger et al. | Seismic detectors: the state of the art | |
| Kouba et al. | Ionospheric drift measurements: Skymap points selection | |
| Beckmann | Interference detection and suppression in Loran-C receivers | |
| Hu et al. | Time Information Transmission Based on FM Broadcast Signal | |
| Malyshev et al. | Algorithm for separating GNSS signals into components |
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 |