CN105071914B - 基于数字调频广播的时钟同步方法和调频广播接收机 - Google Patents
基于数字调频广播的时钟同步方法和调频广播接收机 Download PDFInfo
- Publication number
- CN105071914B CN105071914B CN201510446399.0A CN201510446399A CN105071914B CN 105071914 B CN105071914 B CN 105071914B CN 201510446399 A CN201510446399 A CN 201510446399A CN 105071914 B CN105071914 B CN 105071914B
- Authority
- CN
- China
- Prior art keywords
- frequency
- broadcast
- time
- clock
- broadcast receiver
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000001360 synchronised effect Effects 0.000 claims abstract description 44
- 230000002123 temporal effect Effects 0.000 claims description 39
- 238000011156 evaluation Methods 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 17
- 230000000875 corresponding effect Effects 0.000 description 30
- 238000001228 spectrum Methods 0.000 description 20
- 239000011159 matrix material Substances 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 11
- 238000004891 communication Methods 0.000 description 8
- 238000005070 sampling Methods 0.000 description 8
- 238000000926 separation method Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 230000002269 spontaneous effect Effects 0.000 description 4
- 230000017105 transposition Effects 0.000 description 4
- 230000003071 parasitic effect Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003042 antagnostic effect Effects 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L7/00—Arrangements for synchronising receiver with transmitter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于数字调频广播的时钟同步方法，该方法包括：调频广播接收机接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；调频广播接收机解调出无线数据广播信号中时间同步信号；调频广播接收机获取其接收到时间同步信号的本地时间信息；调频广播接收机根据自身的本地位置信息、时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率。本发明还公开一种调频广播接收机。本发明通过FM调频广播实现伪卫星***或类似伪卫星***中时钟的同步，即通过调频广播发射机将标准时间通过FM频段的数据广播发送至对应的调频广播接收机，解决了伪卫星***或类似伪卫星***的时钟同步过于依赖卫星导航的技术问题。
Description
技术领域
本发明涉及无线导航技术领域，尤其涉及一种基于数字调频广播的时钟同步方法和调频广播接收机。
背景技术
随着人类社会对卫星导航的日益依赖，现有全球卫星导航***(北斗、全球定位***(Global Positioning System，GPS)等，统称为全球卫星导航***(Global NavigationSatellite System，GNSS))信号弱的缺点日益突出：在低层楼宇、丛林和峡谷环境下难以定位；城市高层楼宇内基本无法定位；以隐私保护为名义干扰无线电导航的事件越来越多，进一步降低了导航***的日常可用性；在战争和对抗环境下，无线电导航对抗措施越来越成熟。故此，GNSS的可用性日益遭到质疑。
在这个背景下，世界上主要国家都已着手开展寻找GNSS增强和备份方案的工作，基本出发点都在于寻找和利用GNSS以外的更多导航源，达到增强或者备份现有GNSS的目的。其中，对卫星导航进行备份的应用最为广泛的是地面伪卫星技术。伪卫星(Pseudo-Satellite或Pseudolite，缩写为PL)，是布设于地面上发射某种定位信号的发射器，通常都是发射类似于GPS的信号。理论上，在任何威胁信号无效的地方，伪卫星都可以取代卫星进行导航定位，例如可以室内、地下停车场和隧道等场景进行导航。
但是，伪卫星***或类似伪卫星***面临的基本问题是需要不停地进行多个伪卫星节点之间的时钟同步。与GPS不同，伪卫星通常装备的是温度补偿晶振时钟，精度不高，会产生时钟漂移，在采样时间里与参考站的标准时间信号不能精确同步。但在定位模型中，***中的所有的伪卫星必须要保持同步。现有技术中，一般通过导航卫星来现实伪卫星的时钟同步，但是，如果伪卫星的时钟同步依赖于卫星导航，当卫星导航失效时，伪卫星的时钟精度无法到达导航要求而导致伪卫星失效，导致伪卫星***的时钟过同步于依赖卫星导航。
发明内容
本发明的主要目的在于提供一种基于数字调频广播的时钟同步方法和调频广播接收机，旨在解决伪卫星***或类似伪卫星***的时钟同步过于依赖卫星导航的技术问题。
为实现上述目的，本发明提供的一种基于数字调频广播的时钟同步方法，所述基于数字调频广播的时钟同步方法包括：
调频广播接收机接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；
所述调频广播接收机解调出所述无线数据广播信号中时间同步信号；
所述调频广播接收机获取其接收到所述时间同步信号的本地时间信息；
所述调频广播接收机根据自身的本地位置信息、所述时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率。
优选地，所述时间同步信号包括测距码、所述调频广播发射机广播所述测距码的远端时间信息和所述调频广播发射机的远端位置信息，
所述调频广播接收机根据所述时间同步信号、本地时间信息和本地位置信息，同步该调频广播接收机本地时钟的时刻和频率的步骤包括：
调频广播接收机根据所述本地时间信息，估算该调频广播接收机接收到所述测距码的估算时间信息以及所述调频广播接收机与调频广播发射机之间的相对频差；
调频广播接收机根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差；
调频广播接收机根据所述时间差和相对频差，同步该调频广播接收机本地时钟的时刻和频率。
优选地，所述所述调频广播接收机获取其接收到所述时间同步信号的本地时间信息的步骤包括；
所述调频广播接收机生成与所述调频广播发射机所发送的测距码的脉冲信号波形相同的本地测距码；
所述调频广播接收机获取所述本地测距码的脉冲波形与所述调频广播发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值；
所述调频广播接收机根据所述本地时间信息，估算该调频广播接收机接收到所述测距码的估算时间信息以及所述调频广播接收机与调频广播发射机之间的相对频差的步骤包括：
所述调频广播接收机根据所述调频广播发射机所发送的测距码的各个脉冲与第其一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
优选地，所述调频广播接收机根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差的步骤包括：
所述调频广播接收机根据所述远端位置信息和本地位置信息，计算得出信号自所述调频广播发射机到调频广播接收机的传播时间；
所述调频广播接收机根据所述远端时间信息、估计时间信息和传输时间，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差。
优选地，所述调频广播接收机根据所述时间差和相对频差，调整自身时钟的时刻和频率的步骤包括：
优选地，所述调频广播接收机根据所述时间差和相对频差，调整自身时钟的时刻和频率的步骤包括：
所述调频广播接收机根据所述相对频差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的频率；
所述调频广播接收机根据所述时间差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的时刻。
此外，为实现上述目的，本发明还提供一种调频广播接收机，其特征在于，所述调频广播接收机包括：
信号接收模块，用于接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；
解调模块，用于解调出所述无线数据广播信号中时间同步信号；
本地时间获取模块，用于获取所述调频广播接收机接收到所述时间同步信号的本地时间信息；
同步模块，用于根据自身的本地位置信息、所述时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率。
优选地，所述时间同步信号包括测距码、所述调频广播发射机广播所述测距码的远端时间信息和所述调频广播发射机的远端位置信息，
所述同步模块包括：
估算单元，用于根据所述本地时间信息，估算该调频广播接收机接收到所述测距码的估算时间信息以及所述调频广播接收机与调频广播发射机之间的相对频差；
时间差计算单元，用于根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差；
同步单元，用于根据所述时间差和相对频差，同步该调频广播接收机本地时钟的时刻和频率。
优选地，所述本地时间获取模块包括；
脉冲生成单元，用于生成与所述调频广播发射机所发送的测距码的脉冲信号波形相同的本地测距码；
时钟值获取单元，用于获取所述本地测距码的脉冲波形与所述调频广播发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值；
所述估算单元，还用于根据所述调频广播发射机所发送的测距码的各个脉冲与其第一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
优选地，所述时间差计算单元还用于：
根据所述远端位置信息和本地位置信息，计算得出信号自所述调频广播发射机到调频广播接收机的传播时间；
根据所述远端时间信息、估计时间信息和传输时间，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差。
优选地，所述同步单元还用于：
根据所述相对频差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的频率；
根据所述时间差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的时刻。
本发明通过调频广播接收机接收与标准时间同步的调频广播发射机广播的信号中的无线数据广播信号，然后提取无线数据广播信号中的时间同步信号以及获取该调频广播接收机的本地时间信息和本地位置信息，最后根据本地位置信息、时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率，从而将无线数据广播信号寄生在传统的FM调频信号的空置资源中，由于FM频段的数据广播具有传播距离远、绕射和传输能力强的特点，从而本发明通过FM调频广播实现伪卫星***或类似伪卫星***中时钟的同步，即通过调频广播发射机将标准时间通过FM频段的数据广播发送至对应的调频广播接收机以完成接收机时钟的同步，解决了伪卫星***或类似伪卫星***的时钟同步过于依赖卫星导航的技术问题。
附图说明
图1为本发明基于数字调频广播的时钟同步方法对应***的一个实施例的结构示意图；
图2为本发明基于数字调频广播的时钟同步方法第一实施例的流程示意图；
图3为调频信号与寄生数字信号的带内频谱关系图；
图4为调频信号与寄生数字信号的带外频谱关系图；
图5为本发明基于数字调频广播的时钟同步方法中调频数据广播信号D(t)的信号格式；
图6为根据本发明基于数字调频广播的时钟同步方法中测量脉冲到达时间方法的性能示意图；
图7为本发明调频广播接收机第一实施例的功能模块示意图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
为了更好的理解本发明基于数字调频广播的时钟同步方法，如图1所示，运用本发明基于数字调频广播的时钟同步方法的模型包括至少一个调频广播发射机A0和多个调频广播接收机(例如三个，A1、A2、A3)，调频广播接收机的时钟同步流程相似，为方便叙述，以下各实施例以调频广播接收机A1为例。
本发明提供一种基于数字调频广播的时钟同步方法，在本发明基于数字调频广播的时钟同步方法的第一实施例中，参照图2，基于数字调频广播的时钟同步方法包括：
步骤S10，调频广播接收机接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；
调频广播发射机A0将自己的时间同步到标准时间。一般来讲，最常用的标准时间是世界协调时，同步的方法是利用一个授时型全球导航卫星接收机接收全球导航卫星***的信号，利用授时型全球导航卫星将调频广播发射机A0的时间同步到标准时间。当全球导航卫星***不可用时，调频广播发射机A0可以依赖于自身时钟进行工作，此时，调频广播发射机A0的时间与标准时间会存在误差，但只要误差处于可允许的范围，仍然可视为“处于同步状态”。不同应用对可允许的范围的要求不同。伪卫星***来讲，通常***调频广播发射机A0的时间与世界协调时误差不超过1秒。
调频广播发射机A0(以下简称发射机A0)除广播传统调频信号FM(t)外，还同时广播一种寄生于调频广播的无线数据广播信号D(t)，该信号带有测距码、发射机A0广播测距码的远端时间信息和发射机A0的远端位置信息。因此发射机A0发射的信号可以表示为RF(t)＝FM(t)+D(t)，其中RF(t)为发射机A0发射的全部信号。例如现有技术中，将模拟声音广播信号和数字声音广播信号合成为一路模拟数字混合信号，共用一个模拟声音广播频道，混合调制信号放大后，经天馈***发射，得到预定频谱模式的射频信号。在同一频段内模拟声音广播信号和数字声音广播信号叠加、同时传送。
D(t)信号的频谱设计以不显著影响一般收音机的对现有调频广播的接收音质为出发点。图3给出了全频带频谱形式，也称作带内形式，可以看到D(t)的频谱B2与FM(t)的频谱B1存在交叠。图4给出了一种非全频带频谱形式，也称作带外形式，可以看到D(t)的频谱B2与FM(t)的频谱B1不存在交叠。这两种形式各有优缺点，在实际***中均可见。美国HDRadio采用了图4所示的带外形式，而中国的CDRadio则采用了如图3的带内形式。
一般地，D(t)信号实际能够实现几十Kbps的数据通信功能。根据本发明的D(t)的信号格式如图5所示。根据本发明，我们在D(t)信号中周期的放置测距码，相邻测距码之间为数据块儿。数据块的调制方式通常为正交频分复用(OFDM)信号。具体地，D(t)可以在时间t＝[T0+nTp,T0+nTp+Tpn]时间内发送测距码c(t-nT)，其中T0为第一个测距码的起始时刻，Tp为相邻测距码出现的周期，Tpn为测距码持续时长。测距码可以是伪随机码、Gold码等，且在本发明中，测距码不是连续的，而是在相邻的测距码之间存储数据块。
步骤S20，所述调频广播接收机解调出所述无线数据广播信号中时间同步信号；
调频广播接收机A1将D(t)从RF(t)中解调出来。如果D(t)的频谱与调频信号FM(t)不重叠，则采用滤波方法就可以分开。否则D(t)的频谱与调频信号FM(t)存在重叠，则需要使用带内分离方法，如专利CN201510011206一种带内同频数模音频广播信号的分离方法就公开了一种具体的分离方法。然后调频广播接收机A1对分离后的D(t)去除载波，并进行数字化采样得到D(n)。在数字通信中，去除载波是常规技术。
步骤S30，所述调频广播接收机获取其接收到所述时间同步信号的本地时间信息；
调频广播接收机A1在接收到发射机A0发送的时间同步信号时，记录当时调频广播接收机A1的本地时间信息。接收机A1在本端生成与发射机A0所发送的测距码的脉冲信号波形相同的本地测距码，并且接收机A1获取本地测距码的脉冲波形与发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值(即本地时间信息，此本地时钟值为多个)。
步骤S40，所述调频广播接收机根据自身的本地位置信息、所述时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率。
时间同步信号包括测距码、所述调频广播发射机广播所述测距码的远端时间信息和所述调频广播发射机的远端位置信息，
步骤S40包括：
步骤S41，调频广播接收机根据所述本地时间信息，估算该调频广播接收机接收到所述测距码的估算时间信息以及所述调频广播接收机与调频广播发射机之间的相对频差；
所述调频广播接收机根据所述调频广播发射机所发送的测距码的各个脉冲与第其一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
调频广播接收机A1生成与发射机A0信号c(t)波形相同的本地测距码CL(n)。将CL(n)与D(n)进行滑动相关运算。根据测距码的相关性，当D(n)和CL(n)相位对齐时，此时相关结果最大，即得到一个相关峰。假设在任一给定时间段Tsel内，所述滑动相关运算获得了M个相关峰，第m个相关峰是这一时间段内第m测距码所对应的相关峰，假设第m个相关峰对应的本地时钟值为λm。取其中M乘以rate(一个给定的百分比)个最大的相关峰进行下述估值运算：
待估计的值为：
(1)第一个相关峰对应的脉冲的起始到达本地时刻的最优估值θ；
(2)发射机A0和接收机A1(此处以接收机A1为例)时钟的频差系数k。
已知：
任意第m个脉冲相对于第1个脉冲的标称延迟为tm(t1＝0)，有
T＝[t1,t2,…tm.....,tM]T
其中[]T表示矩阵转置。
现有观测值：
M个脉冲到达时间(本地时钟值)Φ＝[λ1,λ2,λ3,..λm....,λM]T，假设各脉冲到达时间测量噪声为ζ＝[υ1,υ2,υ3,..υm...,υM]T。
观测方程：
第m个相关峰对应的本地时钟值为λm，即第m观测方程可表示为：
λm＝(k+1)tm+θ+υm
或
λm-tm＝k tm+θ+υm
其中k是待求的发射和接收时钟的频差系数；θ是待求的第一个相关峰对应的脉冲的起始到达本地时刻的最优估值，υm是第m个脉冲到达时间测量噪声
矩阵表示：上述观测方程可以表示矩阵形式。
令X＝[k,θ]T
则上述问题可表述为：
Φ-T＝Γ*X+ζ
其中
即
用最小二乘求解，有：
X＝(ΓTΓ)-1ΓT(Φ-T)
X＝[k,θ]T中的最优估值θ和频差系数k即为所求值。其中，θ为所求测距码到达时刻，频差系数k为与自身时钟调频广播发射机A0的相对频差。
需要解释的是，获取λm的技术在通信中是成熟技术，一般取相关峰最大的那个采样点以及超前滞后的相邻采样点(一共三个采样点)进行计算。
步骤S42，调频广播接收机根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差；
接收机A1根据其本地位置信息和发射机A0的远端位置信息，得出接收机A1与发射机A0之间的距离S，且调频广播信号在空气中传输速度接近光速，从而求得信号自发射机A0到接收机A1的传播时间；从而接收机根据远端时间信息、估计时间信息(即上述测距码到达时刻θ)和传输时间，计算得出接收机A1与发射机A0时钟之间的时间差。
步骤S43，调频广播接收机根据所述时间差和相对频差，同步该调频广播接收机本地时钟的时刻和频率。
接收机A1通常在数字电路上用数字秒计数器维护本地时间，该数字电路有一个工作时钟，该工作时钟的频率存在波动，但是频率标称值是已知的(频率标称值设为Fdef)，接收机A1默认启动时，接收机A1认为当数字秒计数器每增加Fdef时认为时间过了一秒，即秒溢出一次。当检测到接收机A1与发射机A0存在时间差时，将秒计时器的值(即接收机A1的时刻)相应的减少时间差乘以频率标称值；当检测到相对频差(即频差系数k)时，将该数字电路的工作时钟实际频率调整为Fdef/(1+K)。
在本实施例中，调频广播接收机接收与标准时间同步的调频广播发射机广播的信号中的无线数据广播信号，然后提取无线数据广播信号中的时间同步信号以及获取该调频广播接收机的本地时间信息和本地位置信息，最后根据本地位置信息、时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率，从而将无线数据广播信号寄生在传统的FM调频信号的空置资源中，由于FM频段的数据广播具有传播距离远、绕射和传输能力强的特点，从而本发明通过FM调频广播实现伪卫星***或类似伪卫星***中时钟的同步，即通过调频广播发射机将标准时间通过FM频段的数据广播发送至对应的调频广播接收机以完成接收机时钟的同步，解决了伪卫星***或类似伪卫星***的时钟同步过于依赖卫星导航的技术问题。
进一步地，在本发明基于数字调频广播的时钟同步方法第一实施例的基础上，提出基于数字调频广播的时钟同步方法的第二实施例，在第二实施例中，步骤S30包括：
步骤S31，所述调频广播接收机生成与所述调频广播发射机所发送的测距码的脉冲信号波形相同的本地测距码；
步骤S32，所述调频广播接收机获取所述本地测距码的脉冲波形与所述调频广播发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值；
调频广播接收机A1生成与发射机A0信号c(t)波形相同的本地测距码CL(n)。将CL(n)与D(n)进行滑动相关运算。根据测距码的相关性，当D(n)和CL(n)相位对齐时，此时相关结果最大，即得到一个相关峰，记录每个相关峰产生时的本地时钟值(即当前的时刻)。
步骤S41包括：
步骤S411，所述调频广播接收机根据所述调频广播发射机所发送的测距码的各个脉冲与第其一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
调频广播接收机A1生成与发射机A0信号c(t)波形相同的本地测距码CL(n)。将CL(n)与D(n)进行滑动相关运算。根据测距码的相关性，当D(n)和CL(n)相位对齐时，此时相关结果最大，即得到一个相关峰。假设在任一给定时间段Tsel内，所述滑动相关运算获得了M个相关峰，第m个相关峰是这一时间段内第m测距码所对应的相关峰，假设第m个相关峰对应的本地时钟值为λm。取其中M乘以rate(一个给定的百分比)个最大的相关峰进行下述估值运算：
待估计的值为：
(1)第一个相关峰对应的脉冲的起始到达本地时刻的最优估值θ；
(2)发射机A0和接收机A1(此处以接收机A1为例)时钟的频差系数k。
已知：
任意第m个脉冲相对于第1个脉冲的标称延迟为tm(t1＝0)，有
T＝[t1,t2,…tm.....,tM]T
其中[]T表示矩阵转置。
现有观测值：
M个脉冲到达时间(本地时钟值)Φ＝[λ1,λ2,λ3,..λm....,λM]T，假设各脉冲到达时间测量噪声为ζ＝[υ1,υ2,υ3,..υm...,υM]T。
观测方程：
第m个相关峰对应的本地时钟值为λm，即第m观测方程可表示为：
λm＝(k+1)tm+θ+υm
或
λm-tm＝k tm+θ+υm
其中k是待求的发射和接收时钟的频差系数；θ是待求的第一个相关峰对应的脉冲的起始到达本地时刻的最优估值，υm是第m个脉冲到达时间测量噪声
矩阵表示：上述观测方程可以表示矩阵形式。
令X＝[k,θ]T
则上述问题可表述为：
Φ-T＝Γ*X+ζ
其中
即
用最小二乘求解，有：
X＝(ΓTΓ)-1ΓT(Φ-T)
X＝[k,θ]T中的最优估值θ和频差系数k即为所求值。其中，θ为所求测距码到达时刻，频差系数k为与自身时钟调频广播发射机A0的相对频差。需要解释的是，获取λm的技术在通信中是成熟技术，一般取相关峰最大的那个采样点以及超前滞后的相邻采样点(一共三个采样点)进行计算。
已有的导航技术中(例如GPS导航、北斗导航等)，带宽一般达到2Mhz以上，且测距码是连续的，通过码相位跟踪就可以达到亚米级的码相位分辨率，而在本发明中测距码不是连续的，相邻的测距码之间存储有数据块，D(t)信号的带宽只有几百KHz(带宽窄)，通常采用1Msps的采样时钟去采样，此时时钟的分辨率只有1us，这样的同步误差显然太大了。
在第二实施例中，通过综合使用多个相关峰估算出第一个相关峰对应的脉冲的起始到达本地时刻的最优估值θ，以及发射机A0和接收机A1(此处以接收机A1为例)时钟的频差系数k，从而对接收机A1的本地时钟时刻和频率进行调整，以与发射机A0同步，实现了在频率比较低的采样时钟和窄带信号的情况下取得了高精度的时间同步，参照图6，图6给出了使用48个脉冲，脉冲使用率rate分别为0.6，0.8，1时最优估值θ均方误差，经测量同步的精度通常可以达到10纳秒量级，其中，图6中的横坐标Eb/N0为接收到的无线电信号的每比特数据的能量与噪声功率密度的比，是本领域中表示接收到的无线电信号质量的通用方法；图6的纵坐标RMSE为θ的估计值与真实值之间的方均根误差，单位为纳秒，是本领域中表示估值质量的通用方法。
在本发明基于数字调频广播的时钟同步方法第二实施例的基础上，提出基于数字调频广播的时钟同步方法的第三实施例，在第三实施例中，步骤S42包括：
步骤S421，所述调频广播接收机根据所述远端位置信息和本地位置信息，计算得出信号在所述调频广播发射机和调频广播接收机的传播时间；
设发射机A0的远端位置信息为Ptx＝[xtx,ytx,ztx]，接收机A1的本地位置信息为Prx＝[xrx,yrx,zrx]，计算信号自发射机A0到接收机A1的传播时间Tprop＝|Ptx-Prx|/C，其中C为光速。
步骤S422，所述调频广播接收机根据所述远端时间信息、估计时间信息和传输时间，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差。
远端时间信息包括在预设时间内(如Tsel内)第一个脉冲对应的发射机A0时间为θtx，接收机接收到第一个脉冲对应的本地时间为θrx＝θ，发射机A0与接收机A1时钟的时间差TE＝θtx+Tprop-θrx。
在本实施例中，通过根据接收机的本地位置信息和发射机的远端位置信息计算得出信号在发射机和接收机之间的传播时间，再根据发射机的远端时间信息、接收机的估计时间信息和传输时间，计算得出接收机时钟与发射机时钟之间的时间差，以一种简单的计算方式获得接收机时钟与发射机时钟之间的时间差。
在本发明基于数字调频广播的时钟同步方法第三实施例的基础上，提出基于数字调频广播的时钟同步方法的第四实施例，在第四实施例中，
步骤43包括：
步骤S431，所述调频广播接收机根据所述相对频差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的频率；
步骤S432，所述调频广播接收机根据所述时间差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的时刻。
在本实施例中，调频广播接收机通常在数字电路上用数字秒计数器来维护本地时间，该数字电路有一个工作时钟，其频率存在波动，但是频率标称值(单位：MHz)是已知的，设定该频率标称值为Fdef。接收机A1默认启动时，若接收机A1的数字秒计时器每增加Fdef则时间过了一秒，即秒溢出一次。当检测到发射机A0与接收机A1存在时间差TE时，则将秒计数器的值相应的减少{TE*Fdef}(需要四舍五入，也可能是负值)。当检测到频差系数k时，说明该数字电路的一个工作时钟实际频率为Fdef/(1+k)，因此调频广播接收机A1应调整秒溢出机制，变为当数字秒计数器达到频率标称值增加了Fdef/(1+K)(需要四舍五入)时认为时间过了一秒。
本发明提供一种调频广播接收机，在本发明调频广播接收机的第一实施例中，参照图7，调频广播接收机包括：
信号接收模块10，用于接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；
调频广播发射机A0将自己的时间同步到标准时间。一般来讲，最常用的标准时间是世界协调时，同步的方法是利用一个授时型全球导航卫星接收机接收全球导航卫星***的信号，利用授时型全球导航卫星将调频广播发射机A0的时间同步到标准时间。当全球导航卫星***不可用时，调频广播发射机A0可以依赖于自身时钟进行工作，此时，调频广播发射机A0的时间与标准时间会存在误差，但只要误差处于可允许的范围，仍然可视为“处于同步状态”。不同应用对可允许的范围的要求不同。伪卫星***来讲，通常***调频广播发射机A0的时间与世界协调时误差不超过1秒。
调频广播发射机A0(以下简称发射机A0)除广播传统调频信号FM(t)外，还同时广播一种寄生于调频广播的无线数据广播信号D(t)，该信号带有测距码、发射机A0广播测距码的远端时间信息和发射机A0的远端位置信息。因此发射机A0发射的信号可以表示为RF(t)＝FM(t)+D(t)，其中RF(t)为发射机A0发射的全部信号。例如现有技术中，将模拟声音广播信号和数字声音广播信号合成为一路模拟数字混合信号，共用一个模拟声音广播频道，混合调制信号放大后，经天馈***发射，得到预定频谱模式的射频信号。在同一频段内模拟声音广播信号和数字声音广播信号叠加、同时传送。
D(t)信号的频谱设计以不显著影响一般收音机的对现有调频广播的接收音质为出发点。图3给出了全频带频谱形式，也称作带内形式，可以看到D(t)的频谱B2与FM(t)的频谱B1存在交叠。图4给出了一种非全频带频谱形式，也称作带外形式，可以看到D(t)的频谱B2与FM(t)的频谱B1不存在交叠。这两种形式各有优缺点，在实际***中均可见。美国HDRadio采用了图4所示的带外形式，而中国的CDRadio则采用了如图3的带内形式。
一般地，D(t)信号实际能够实现几十Kbps的数据通信功能。根据本发明的D(t)的信号格式如图5所示。根据本发明，我们在D(t)信号中周期的放置测距码，相邻测距码之间为数据块儿。数据块的调制方式通常为正交频分复用(OFDM)信号。具体地，D(t)可以在时间t＝[T0+nTp,T0+nTp+Tpn]时间内发送测距码c(t-nT)，其中T0为第一个测距码的起始时刻，Tp为相邻测距码出现的周期，Tpn为测距码持续时长。测距码可以是伪随机码、Gold码等，且在本发明中，测距码不是连续的，而是在相邻的测距码之间存储数据块。
解调模块20，用于解调出所述无线数据广播信号中时间同步信号；
调频广播接收机A1的解调模块20将D(t)从RF(t)中解调出来。如果D(t)的频谱与调频信号FM(t)不重叠，则采用滤波方法就可以分开。否则D(t)的频谱与调频信号FM(t)存在重叠，则需要使用带内分离方法，如专利CN201510011206一种带内同频数模音频广播信号的分离方法就公开了一种具体的分离方法。然后调频广播接收机A1的解调模块20对分离后的D(t)去除载波，并进行数字化采样得到D(n)。在数字通信中，去除载波是常规技术。
本地时间获取模块30，用于获取所述调频广播接收机接收到所述时间同步信号的本地时间信息；
调频广播接收机A1的本地时间获取模块30在接收到发射机A0发送的时间同步信号时，记录当时调频广播接收机A1的本地时间信息。接收机A1在本端生成与发射机A0所发送的测距码的脉冲信号波形相同的本地测距码，并且接收机A1获取本地测距码的脉冲波形与发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值(即本地时间信息，此本地时钟值为多个)。
同步模块40，用于根据自身的本地位置信息、所述时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率。
时间同步信号包括测距码、所述调频广播发射机广播所述测距码的远端时间信息和所述调频广播发射机的远端位置信息，
所述同步模块40包括：
估算单元41，用于根据所述本地时间信息，估算该调频广播接收机接收到所述测距码的估算时间信息以及所述调频广播接收机与调频广播发射机之间的相对频差；
所述调频广播接收机根据所述调频广播发射机所发送的测距码的各个脉冲与第其一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
调频广播接收机A1生成与发射机A0信号c(t)波形相同的本地测距码CL(n)。将CL(n)与D(n)进行滑动相关运算。根据测距码的相关性，当D(n)和CL(n)相位对齐时，此时相关结果最大，即得到一个相关峰。假设在任一给定时间段Tsel内，所述滑动相关运算获得了M个相关峰，第m个相关峰是这一时间段内第m测距码所对应的相关峰，假设第m个相关峰对应的本地时钟值为λm。取其中M乘以rate(一个给定的百分比)个最大的相关峰进行下述估值运算：
待估计的值为：
(1)第一个相关峰对应的脉冲的起始到达本地时刻的最优估值θ；
(2)发射机A0和接收机A1(此处以接收机A1为例)时钟的频差系数k。
已知：
任意第m个脉冲相对于第1个脉冲的标称延迟为tm(t1＝0)，有
T＝[t1,t2,…tm.....,tM]T
其中[]T表示矩阵转置。
现有观测值：
M个脉冲到达时间(本地时钟值)Φ＝[λ1,λ2,λ3,..λm....,λM]T，假设各脉冲到达时间测量噪声为ζ＝[υ1,υ2,υ3,..υm...,υM]T。
观测方程：
第m个相关峰对应的本地时钟值为λm，即第m观测方程可表示为：
λm＝(k+1)tm+θ+υm
或
λm-tm＝k tm+θ+υm
其中k是待求的发射和接收时钟的频差系数；θ是待求的第一个相关峰对应的脉冲的起始到达本地时刻的最优估值，υm是第m个脉冲到达时间测量噪声
矩阵表示：上述观测方程可以表示矩阵形式。
令X＝[k,θ]T
则上述问题可表述为：
Φ-T＝Γ*X+ζ
其中
即
用最小二乘求解，有：
X＝(ΓTΓ)-1ΓT(Φ-T)
X＝[k,θ]T中的最优估值θ和频差系数k即为所求值。其中，θ为所求测距码到达时刻，频差系数k为与自身时钟调频广播发射机A0的相对频差。
需要解释的是，获取λm的技术在通信中是成熟技术，一般取相关峰最大的那个采样点以及超前滞后的相邻采样点(一共三个采样点)进行计算。
时间差计算单元42，用于根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差；
接收机A1的时间差计算单元42根据其本地位置信息和发射机A0的远端位置信息，得出接收机A1与发射机A0之间的距离S，且调频广播信号在空气中传输速度接近光速，从而求得信号自发射机A0到接收机A1的传播时间；从而接收机根据远端时间信息、估计时间信息(即上述测距码到达时刻θ)和传输时间，计算得出接收机A1与发射机A0时钟之间的时间差。
同步单元43，用于根据所述时间差和相对频差，同步该调频广播接收机本地时钟的时刻和频率。
接收机A1的同步单元43通常在数字电路上用数字秒计数器维护本地时间，该数字电路有一个工作时钟，该工作时钟的频率存在波动，但是频率标称值是已知的(频率标称值设为Fdef)，接收机A1默认启动时，接收机A1认为当数字秒计数器每增加Fdef时认为时间过了一秒，即秒溢出一次。当检测到接收机A1与发射机A0存在时间差时，将秒计时器的值(即接收机A1的时刻)相应的减少时间差乘以频率标称值；当检测到相对频差(即频差系数k)时，将该数字电路的工作时钟实际频率调整为Fdef/(1+K)。
在本实施例中，调频广播接收机的信号接收模块10接收与标准时间同步的调频广播发射机广播的信号中的无线数据广播信号，然后解调模块20提取无线数据广播信号中的时间同步信号以及本地时间获取模块30获取该调频广播接收机的本地时间信息和本地位置信息，最后同步模块40根据本地位置信息、时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率，从而将无线数据广播信号寄生在传统的FM调频信号的空置资源中，由于FM频段的数据广播具有传播距离远、绕射和传输能力强的特点，从而本发明通过FM调频广播实现伪卫星***或类似伪卫星***中时钟的同步，即通过调频广播发射机将标准时间通过FM频段的数据广播发送至对应的调频广播接收机以完成接收机时钟的同步，解决了伪卫星***或类似伪卫星***的时钟同步过于依赖卫星导航的技术问题。
进一步地，在本发明调频广播接收机第一实施例的基础上，提出调频广播接收机的第二实施例，在第二实施例中，所述本地时间获取模块30包括；
脉冲生成单元31，用于生成与所述调频广播发射机所发送的测距码的脉冲信号波形相同的本地测距码；
时钟值获取单元32，用于获取所述本地测距码的脉冲波形与所述调频广播发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值；
调频广播接收机A1的脉冲生成单元31生成与发射机A0信号c(t)波形相同的本地测距码CL(n)。时钟值获取单元32将CL(n)与D(n)进行滑动相关运算。根据测距码的相关性，当D(n)和CL(n)相位对齐时，此时相关结果最大，即得到一个相关峰，记录每个相关峰产生时的本地时钟值(即当前的时刻)。
所述估算单元41，还用于根据所述调频广播发射机所发送的测距码的各个脉冲与其第一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
调频广播接收机A1生成与发射机A0信号c(t)波形相同的本地测距码CL(n)。将CL(n)与D(n)进行滑动相关运算。估算单元41根据测距码的相关性，当D(n)和CL(n)相位对齐时，此时相关结果最大，即得到一个相关峰。假设在任一给定时间段Tsel内，所述滑动相关运算获得了M个相关峰，第m个相关峰是这一时间段内第m测距码所对应的相关峰，假设第m个相关峰对应的本地时钟值为λm。取其中M乘以rate(一个给定的百分比)个最大的相关峰进行下述估值运算：
待估计的值为：
(1)第一个相关峰对应的脉冲的起始到达本地时刻的最优估值θ；
(2)发射机A0和接收机A1(此处以接收机A1为例)时钟的频差系数k。
已知：
任意第m个脉冲相对于第1个脉冲的标称延迟为tm(t1＝0)，有
T＝[t1,t2,…tm.....,tM]T
其中[]T表示矩阵转置。
现有观测值：
M个脉冲到达时间(本地时钟值)Φ＝[λ1,λ2,λ3,..λm....,λM]T，假设各脉冲到达时间测量噪声为ζ＝[υ1,υ2,υ3,..υm...,υM]T。
观测方程：
第m个相关峰对应的本地时钟值为λm，即第m观测方程可表示为：
λm＝(k+1)tm+θ+υm
或
λm-tm＝k tm+θ+υm
其中k是待求的发射和接收时钟的频差系数；θ是待求的第一个相关峰对应的脉冲的起始到达本地时刻的最优估值，υm是第m个脉冲到达时间测量噪声
矩阵表示：上述观测方程可以表示矩阵形式。
令X＝[k,θ]T
则上述问题可表述为：
Φ-T＝Γ*X+ζ
其中
即
用最小二乘求解，有：
X＝(ΓTΓ)-1ΓT(Φ-T)
X＝[k,θ]T中的最优估值θ和频差系数k即为所求值。其中，θ为所求测距码到达时刻，频差系数k为与自身时钟调频广播发射机A0的相对频差。需要解释的是，获取λm的技术在通信中是成熟技术，一般取相关峰最大的那个采样点以及超前滞后的相邻采样点(一共三个采样点)进行计算。
已有的导航技术中(例如GPS导航、北斗导航等)，带宽一般达到2Mhz以上，且测距码是连续的，通过码相位跟踪就可以达到亚米级的码相位分辨率，而在本发明中测距码不是连续的，相邻的测距码之间存储有数据块，D(t)信号的带宽只有几百KHz(带宽窄)，通常采用1Msps的采样时钟去采样，此时时钟的分辨率只有1us，这样的同步误差显然太大了。
在第二实施例中，通过估算单元41综合使用多个相关峰估算出第一个相关峰对应的脉冲的起始到达本地时刻的最优估值θ，以及发射机A0和接收机A1(此处以接收机A1为例)时钟的频差系数k，从而对接收机A1的本地时钟时刻和频率进行调整，以与发射机A0同步，实现了在频率比较低的采样时钟和窄带信号的情况下取得了高精度的时间同步，参照图6，图6给出了使用48个脉冲，脉冲使用率rate分别为0.6，0.8，1时最优估值θ均方误差，经测量同步的精度通常可以达到10纳秒量级，其中，图6中的横坐标Eb/N0为接收到的无线电信号的每比特数据的能量与噪声功率密度的比，是本领域中表示接收到的无线电信号质量的通用方法；图6的纵坐标RMSE为θ的估计值与真实值之间的方均根误差，单位为纳秒，是本领域中表示估值质量的通用方法。
进一步地，在本发明调频广播接收机第二施例的基础上，提出调频广播接收机的第三实施例，在第三实施例中，所述时间差计算单元42还用于：
根据所述远端位置信息和本地位置信息，计算得出信号自所述调频广播发射机到调频广播接收机的传播时间；
设发射机A0的远端位置信息为Ptx＝[xtx,ytx,ztx]，接收机A1的本地位置信息为Prx＝[xrx,yrx,zrx]，时间差计算单元42计算信号自发射机A0到接收机A1的传播时间Tprop＝|Ptx-Prx|/C，其中C为光速。
根据所述远端时间信息、估计时间信息和传输时间，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差。
远端时间信息包括在预设时间内(如Tsel内)第一个脉冲对应的发射机A0时间为θtx，接收机接收到第一个脉冲对应的本地时间为θrx＝θ，发射机A0与接收机A1时钟的时间差TE＝θtx+Tprop-θrx。
在本实施例中，通过时间差计算单元42根据接收机的本地位置信息和发射机的远端位置信息计算得出信号在发射机和接收机之间的传播时间，再根据发射机的远端时间信息、接收机的估计时间信息和传输时间，计算得出接收机时钟与发射机时钟之间的时间差，以一种简单的计算方式获得接收机时钟与发射机时钟之间的时间差。
进一步地，在本发明调频广播接收机第三施例的基础上，提出调频广播接收机的第四实施例，在第四实施例中，所述同步单元43还用于：
根据所述相对频差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的频率；
根据所述时间差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的时刻。
在本实施例中，调频广播接收机通常在数字电路上用数字秒计数器来维护本地时间，该数字电路有一个工作时钟，其频率存在波动，但是频率标称值(单位：MHz)是已知的，设定该频率标称值为Fdef。接收机A1默认启动时，若接收机A1的数字秒计时器每增加Fdef则时间过了一秒，即秒溢出一次。当检测到发射机A0与接收机A1存在时间差TE时，则将秒计数器的值相应的减少{TE*Fdef}(需要四舍五入，也可能是负值)。当检测到频差系数k时，说明该数字电路的一个工作时钟实际频率为Fdef/(1+k)，因此调频广播接收机A1应调整秒溢出机制，变为当数字秒计数器达到频率标称值增加了Fdef/(1+K)(需要四舍五入)时认为时间过了一秒。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例的方法。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (6)
1.一种基于数字调频广播的时钟同步方法，其特征在于，所述基于数字调频广播的时钟同步方法包括：
调频广播接收机接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；
所述调频广播接收机解调出所述无线数据广播信号中时间同步信号；所述时间同步信号包括测距码、所述调频广播发射机广播所述测距码的远端时间信息和所述调频广播发射机的远端位置信息；
所述调频广播接收机生成与所述调频广播发射机所发送的测距码的脉冲信号波形相同的本地测距码；
所述调频广播接收机获取所述本地测距码的脉冲波形与所述调频广播发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值；
所述调频广播接收机根据所述调频广播发射机所发送的测距码的各个脉冲与第其一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差；
调频广播接收机根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差；
调频广播接收机根据所述时间差和相对频差，同步该调频广播接收机本地时钟的时刻和频率。
2.如权利要求1所述的基于数字调频广播的时钟同步方法，其特征在于，所述调频广播接收机根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差的步骤包括：
所述调频广播接收机根据所述远端位置信息和本地位置信息，计算得出信号自所述调频广播发射机到调频广播接收机的传播时间；
所述调频广播接收机根据所述远端时间信息、估计时间信息和传输时间，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差。
3.如权利要求1所述的基于数字调频广播的时钟同步方法，其特征在于，所述调频广播接收机根据所述时间差和相对频差，调整自身时钟的时刻和频率的步骤包括：
所述调频广播接收机根据所述相对频差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的频率；
所述调频广播接收机根据所述时间差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的时刻。
4.一种调频广播接收机，其特征在于，所述调频广播接收机包括：
信号接收模块，用于接收与标准时间同步的调频广播发射机所广播的无线数据广播信号；
解调模块，用于解调出所述无线数据广播信号中时间同步信号；所述时间同步信号包括测距码、所述调频广播发射机广播所述测距码的远端时间信息和所述调频广播发射机的远端位置信息；
本地时间获取模块，用于获取所述调频广播接收机接收到所述时间同步信号的本地时间信息；
同步模块，用于根据自身的本地位置信息、所述时间同步信号和本地时间信息，同步该调频广播接收机本地时钟的时刻和频率；
所述同步模块包括：
估算单元，用于根据所述本地时间信息，估算该调频广播接收机接收到所述测距码的估算时间信息以及所述调频广播接收机与调频广播发射机之间的相对频差；
时间差计算单元，用于根据所述远端位置信息、本地位置信息、远端时间信息和估算时间信息，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差；
同步单元，用于根据所述时间差和相对频差，同步该调频广播接收机本地时钟的时刻和频率；
所述本地时间获取模块包括；
脉冲生成单元，用于生成与所述调频广播发射机所发送的测距码的脉冲信号波形相同的本地测距码；
时钟值获取单元，用于获取所述本地测距码的脉冲波形与所述调频广播发射机所发送的测距码的脉冲波形相位对齐时的本地时钟值；
所述估算单元，还用于根据所述调频广播发射机所发送的测距码的各个脉冲与其第一个脉冲的标称延迟和所述本地时钟值，估算该调频广播接收机接收到所述测距码的第一个脉冲的估算时钟值以及所述调频广播接收机与调频广播发射机之间的相对频差。
5.如权利要求4所述的调频广播接收机，其特征在于，所述时间差计算单元还用于：
根据所述远端位置信息和本地位置信息，计算得出信号自所述调频广播发射机到调频广播接收机的传播时间；
根据所述远端时间信息、估计时间信息和传输时间，计算得出所述调频广播接收机时钟与调频广播发射机时钟之间的时间差。
6.如权利要求4所述的调频广播接收机，其特征在于，所述同步单元还用于：
根据所述相对频差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的频率；
根据所述时间差和所述调频广播接收机的频率标称值，调整所述调频广播接收机时钟的时刻。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510446399.0A CN105071914B (zh) | 2015-07-27 | 2015-07-27 | 基于数字调频广播的时钟同步方法和调频广播接收机 |
| PCT/CN2016/085079 WO2017016321A1 (zh) | 2015-07-27 | 2016-06-07 | 基于数字调频广播的时钟同步方法和调频广播接收机 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510446399.0A CN105071914B (zh) | 2015-07-27 | 2015-07-27 | 基于数字调频广播的时钟同步方法和调频广播接收机 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105071914A CN105071914A (zh) | 2015-11-18 |
| CN105071914B true CN105071914B (zh) | 2018-02-02 |
Family
ID=54501202
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510446399.0A Active CN105071914B (zh) | 2015-07-27 | 2015-07-27 | 基于数字调频广播的时钟同步方法和调频广播接收机 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN105071914B (zh) |
| WO (1) | WO2017016321A1 (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105071914B (zh) * | 2015-07-27 | 2018-02-02 | 深圳思凯微电子有限公司 | 基于数字调频广播的时钟同步方法和调频广播接收机 |
| CN106773985A (zh) * | 2016-12-29 | 2017-05-31 | 西北核技术研究所 | 一种用于远距离多点控制的高精度顺序控制单元和方法 |
| CN107589429B (zh) * | 2017-08-14 | 2020-05-01 | 深圳思凯微电子有限公司 | 基于调频数据广播的定位方法、装置、***及存储介质 |
| IL277231B2 (en) * | 2018-03-14 | 2024-03-01 | Locata Corp | Method and device for synchronizing a tracking network |
| CN109565772A (zh) * | 2018-10-31 | 2019-04-02 | 深圳市汇顶科技股份有限公司 | 时间同步方法、设备及存储介质 |
| KR102036649B1 (ko) * | 2019-06-17 | 2019-10-25 | 주식회사 에스씨솔루션 | 누수관리 시스템 및 이를 이용한 누수위치 예측 방법 |
| CN113037413B (zh) * | 2021-03-11 | 2023-04-21 | 成都德芯数字科技股份有限公司 | 一种调频同步方法、装置、***、设备、存储介质 |
| CN113132046B (zh) * | 2021-03-25 | 2022-10-11 | 中国电子科技集团公司第五十四研究所 | 一种基于锁模光频梳的共视法时间同步装置及方法 |
| CN114326362B (zh) * | 2021-12-30 | 2023-08-08 | 国网思极位置服务有限公司 | 一种寄生于调频电台的无线授时***和方法 |
| CN116847452B (zh) * | 2023-08-03 | 2024-03-15 | 大有期货有限公司 | 卫星信号校时同步*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102957642A (zh) * | 2011-08-24 | 2013-03-06 | 上海凯芯微电子有限公司 | 一种无线数据接收***及其接收方法 |
| CN103630915A (zh) * | 2012-08-24 | 2014-03-12 | 陈曦 | 一种利用数字调频广播进行导航定位的方法 |
| CN103647628A (zh) * | 2013-11-28 | 2014-03-19 | 陈辉 | 一种时间同步方法、装置及*** |
| CN104062895A (zh) * | 2014-06-26 | 2014-09-24 | 桂林电子科技大学 | 伪卫星时间同步及其定位方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5389039B2 (ja) * | 2007-11-02 | 2014-01-15 | ノヴァテル インコーポレイテッド | 正確な時間および周波数をネットワークにわたって分配するためのシステムおよび方法 |
| US8810452B2 (en) * | 2010-05-24 | 2014-08-19 | Trueposition, Inc. | Network location and synchronization of peer sensor stations in a wireless geolocation network |
| US20130093619A1 (en) * | 2011-10-07 | 2013-04-18 | Electronics And Telecommunications Research Institute | Apparatus and method for indoor positioning |
| CN102355633A (zh) * | 2011-10-11 | 2012-02-15 | 中国科学院软件研究所 | 一种基于fm广播数据***的时钟校准方法 |
| CN105071914B (zh) * | 2015-07-27 | 2018-02-02 | 深圳思凯微电子有限公司 | 基于数字调频广播的时钟同步方法和调频广播接收机 |
-
2015
- 2015-07-27 CN CN201510446399.0A patent/CN105071914B/zh active Active
-
2016
- 2016-06-07 WO PCT/CN2016/085079 patent/WO2017016321A1/zh active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102957642A (zh) * | 2011-08-24 | 2013-03-06 | 上海凯芯微电子有限公司 | 一种无线数据接收***及其接收方法 |
| CN103630915A (zh) * | 2012-08-24 | 2014-03-12 | 陈曦 | 一种利用数字调频广播进行导航定位的方法 |
| CN103647628A (zh) * | 2013-11-28 | 2014-03-19 | 陈辉 | 一种时间同步方法、装置及*** |
| CN104062895A (zh) * | 2014-06-26 | 2014-09-24 | 桂林电子科技大学 | 伪卫星时间同步及其定位方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105071914A (zh) | 2015-11-18 |
| WO2017016321A1 (zh) | 2017-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105071914B (zh) | 基于数字调频广播的时钟同步方法和调频广播接收机 | |
| US6963306B2 (en) | Position location and data transmission using pseudo digital television transmitters | |
| Rabinowitz et al. | A new positioning system using television synchronization signals | |
| EP2384449B1 (en) | Local clock frequency calibration using low earth orbit (leo) satellites | |
| CN104076371B (zh) | 利用卫星***的因特网热点定位 | |
| EP0902300B1 (en) | Telecommunications-assisted satellite positioning system | |
| US7042396B2 (en) | Position location using digital audio broadcast signals | |
| EP0902554B1 (en) | An auxiliary system for assisting a wireless terminal in determining its position from signals transmitted from a navigation satellite | |
| JP3240393B2 (ja) | 世界測位システム用トラッキングシステム | |
| CN100381835C (zh) | 基于数字广播电视信号的无线电组合定位方法 | |
| US8233091B1 (en) | Positioning and time transfer using television synchronization signals | |
| US20030052822A1 (en) | Position location using broadcast digital television signals | |
| US20070063893A1 (en) | Spot Locator | |
| US7432852B2 (en) | Method of synchronizing base stations of a terrestrial cellular communication network | |
| JP2011503582A (ja) | ネットワークにわたって位置を定めるためのシステム | |
| CN106908817A (zh) | 辅助导航定位方法及*** | |
| US20040201779A1 (en) | Symbol clock recovery for the ATSC digital television signal | |
| CN103630915B (zh) | 一种利用数字调频广播进行导航定位的方法 | |
| JP2009025292A (ja) | テレビジョン信号の監視ユニット | |
| KR101475036B1 (ko) | ＧＰＳ 및 Ｇａｌｉｌｅｏ 위성 신호의 주피크 특성을 이용한 위성 신호 멀티패스 추적 시스템 | |
| CN205384374U (zh) | 一种都市环境中的组合定位装置 | |
| JP5181344B2 (ja) | ナビゲーション信号送信装置、受信機、ナビゲーション信号の生成方法およびナビゲーション信号生成プログラム | |
| KR100766484B1 (ko) | 위성 이용 세계 국가 표준시 다중 방송방법 | |
| US7161533B2 (en) | Method and a system for positioning, and an electronic device | |
| CN107678047B (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 | ||
| CB02 | Change of applicant information |
Address after: 518000 Shenzhen integrated circuit design and application Industrial Park, No. 1089 Tea Light Road, Shenzhen, Guangdong 508-1, Nanshan District, China Applicant after: SHENZHEN SKYCASTER MICRONICS Inc. Address before: 518000 design & Application Industrial Park, Shenzhen, Nanshan District, Shenzhen, 1089, Guang Guang road, Guangdong, 508 Applicant before: SHENZHEN SKYCASTER MICRONICS Inc. |
|
| COR | Change of bibliographic data | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhang Guanghua Inventor after: Li Li Inventor after: Peng Tieyan Inventor before: Chen Xi Inventor before: Zhang Guanghua Inventor before: Li Li Inventor before: Peng Tieyan |