# US 11,442,173 B2 · GNSS-based timing synchronization > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 11,442,173 B2 - Document ID: US11442173B2 - Status: Granted - Kind: B2 - Application number: 16/261,718 - Filed: 2019-01-30 - Published: 2020-07-30 - Granted: 2022-09-13 - Earliest claimed: 2019-01-30 - Assignees: Mitsubishi Electric Research Laboratories, Inc. - Applicants: Mitsubishi Electric Research Laboratories, Inc. - Inventors: Kyeong-Jin Kim; Sriramya Bhamidipati; Hongbo Sun - Web page: https://patentgrep.com/patent/US11442173B2 ## CPC - G01S 19/215: {issues related to spoofing} - G01S 19/36: relating to the receiver frond end - G01S 19/39: the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO - G01S: RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES ## Abstract A time synchronization system an input interface configured to receive pseudorange measurements from a set of antennas of at least one Global Navigation Satellite System (GNSS) receiving system that operates the antennas at known positions according to a clock with an unknown time bias to receive signals from satellites in sight of the antennas and an output interface configured to output the determined time bias. The time synchronization system also includes a processor configured to compare the pseudorange measurements of different antennas to produce an antenna-specific clock offset for each of the antennas, correct the pseudorange measurements of the antennas according to the antenna-specific clock offsets of corresponding antennas to produce corrected pseudorange measurements, and determine the time bias of the clock using the corrected pseudorange measurements. ## Claims ### Claim 1 A time synchronization system, comprising: an input interface configured to receive pseudorange measurements from a set of antennas of at least one Global Navigation Satellite System (GNSS) receiving system that operates the set of antennas at known positions according to a clock with an unknown time bias to receive signals from satellites in sight of each antenna of the set of antennas; a processor configured to: compare the pseudorange measurements of the set of antennas to produce an antenna-specific clock offset for each antenna of the set of antennas, wherein to produce the antenna-specific clock offset for each antenna of the set of antennas, the processor is configured to: determine, for each pair of antennas from the set of antennas, a difference pseudorange residual between pseudorange measurements of a pair of antennas of the set of antennas, wherein a pseudorange residual of a first antenna of the set of antennas removes known components from the pseudorange measurement of the first antenna, and wherein the difference pseudorange residual between the pair of antennas of the set of antennas removes common components in the pseudorange measurements of the pair of antennas of the set of antennas, such that the difference pseudorange residual includes a difference between the antenna-specific clock offsets of the pair of antennas of the set of antennas subject to noise; and compare the difference pseudorange residuals between different pairs of antennas of the set of antennas to produce the antenna-specific clock offset for each antenna of the set of antennas; correct the pseudorange measurements of the set of antennas according to the antenna-specific clock offset for each antenna of the set of antennas to produce corrected pseudorange measurements; and determine the time bias of the clock using the corrected pseudorange measurements; and an output interface configured to output the determined time bias. ### Claim 2 (depends on 1) The time synchronization system of claim 1, wherein different antennas of the set of antennas are in sight of different combinations of the satellites, and wherein at least one first antenna of the set of antennas has a first antenna-specific clock offset different from a second antenna-specific clock offset of a second antenna of the set of antennas. ### Claim 3 (depends on 1) The time synchronization system of claim 1, wherein the processor is configured to compare the difference pseudorange residuals while ignoring the noise or by sampling the noise on a known distribution of the noise. ### Claim 4 (depends on 1) The time synchronization system of claim 1, wherein the processor is configured to: for each pair of antennas of the set of antennas, probabilistically determine the difference pseudorange residuals using a Gaussian model for a distribution of difference pseudorange residual according to pseudorange measurements from all satellites in sight of respective pair of antennas of the set of antennas; and sample the distribution of the difference pseudorange residual of each pair of antennas to produce the antenna-specific clock offset for each antenna of the set of antennas by comparing the sampled distribution of the difference pseudorange residual of each pair of antennas. ### Claim 5 (depends on 1) The time synchronization system of claim 1, wherein to produce the antenna-specific clock offset for the first antenna of the set of antennas, the processor is configured to determine a probability of a first antenna-specific clock offset of the first antenna with respect to probabilities of antenna-specific clock offsets of a plurality of antennas of the set of antennas other than the first antenna as a marginal distribution of the first antenna-specific clock offset of the first antenna conditioned on a joint posteriori distribution of the antenna-specific clock offsets of all antennas of the set of antennas. ### Claim 6 (depends on 5) The time synchronization system of claim 5, wherein the antenna-specific clock offsets of all antennas of the set of antennas are determined based on corresponding marginal distributions of the antenna-specific clock offsets. ### Claim 7 (depends on 5) The time synchronization system of claim 5, wherein a second antenna-specific clock offset for a second antenna k of the set of antennas is determined using a belief propagation as a probabilistic belief bt(αk) at time t on a marginal distribution of an antenna-specific clock offset αk using a product of a prior distribution of the probabilistic belief bt(αk) and incoming belief messages from all antennas of the set of antennas except the second antenna k. ### Claim 8 (depends on 7) The time synchronization system of claim 7, wherein the belief propagation updates the probabilistic belief bt(αk) of the probability of the antenna-specific clock offset αk based on a computation of message from a factor node, fk,n, to variable nodes and a computation of message, from a factor node, fk, to a variable of the factor node, fk. ### Claim 9 (depends on 1) The time synchronization system of claim 1, wherein the processor executes a Kalman filter configured to determine the time bias of the clock using the corrected pseudorange measurements of the set of antennas. ### Claim 10 (depends on 1) The time synchronization system of claim 1, wherein the set of antennas includes multiple directional antennas. ### Claim 11 (depends on 1) The time synchronization system of claim 1, wherein the at least one GNSS receiving system includes a set of GNSS receiving systems, wherein each GNSS receiving system of the set of GNSS receiving systems includes one or multiple antennas of the set of antennas. ### Claim 12 (depends on 11) The time synchronization system of claim 11, wherein the processor randomly selects a subset of GNSS receiving systems from the set of GNSS receiving systems to perform time synchronization. ### Claim 13 (depends on 1) The time synchronization system of claim 1, further comprising: a phasor measurement unit configured to operate using the clock with the determined time bias. ### Claim 14 A time synchronization method, wherein the method uses a processor coupled with stored instructions implementing the method, wherein the instructions, when executed by the processor carry out steps of the method, comprising: receiving pseudorange measurements from a set of antennas of at least one Global Navigation Satellite System (GNSS) receiving system that operates the set of antennas at known positions according to a clock with an unknown time bias to receive signals from satellites in sight of each antenna of the set of antennas; comparing the pseudorange measurements of the set of antennas to produce an antenna-specific clock offset for each antenna of the set of antennas, wherein a first antenna-specific clock offset for an antenna k of the set of antennas is determined using a belief propagation as a probabilistic belief bt(αk) at time t on a marginal distribution of an antenna-specific clock offset αk using a product of a prior distribution of the probabilistic belief bt(αk) and incoming belief messages from all antennas of the set of antennas except the antenna k; correcting the pseudorange measurement of the set of antennas according to the antenna-specific clock offset for each antenna of the set of antennas to produce corrected pseudorange measurements; determining the time bias of the clock using the corrected pseudorange measurements; and correcting the time of the clock using the determined time bias. ### Claim 15 (depends on 14) The time synchronization method of claim 14, wherein different antennas of the set of antennas are in sight of different combinations of the satellites, and wherein, at least one first antenna of the set of antennas has a second antenna-specific clock offset different from a third antenna-specific clock offset of a second antenna of the set of antennas. ### Claim 16 (depends on 14) The time synchronization method of claim 14, further comprising: determining a probability of a second antenna-specific clock offset of a first antenna with respect to probabilities of antenna-specific clock offsets of a plurality of antennas of the set of antennas other than the first antenna as a marginal distribution of the second antenna-specific clock offset of the first antenna conditioned on a joint posteriori distribution of the antenna-specific clock offsets of all antennas of the set of antennas. ### Claim 17 A non-transitory computer readable storage medium embodied thereon a program executable by a processor for performing a method, the method comprising: receiving pseudorange measurements from a set of antennas of at least one Global Navigation Satellite System (GNSS) receiving system that operates the set of antennas at known positions according to a clock with an unknown time bias to receive signals from satellites in sight of each antenna of the set of antennas; comparing the pseudorange measurements of the set of antennas to produce an antenna-specific clock offset for each antenna of the set of antennas by: determining, for each pair of antennas from the set of antennas, a difference pseudorange residual between pseudorange measurements of a pair of antennas of the set of antennas, wherein a pseudorange residual of a first antenna of the set of antennas removes known components from the pseudorange measurement of the first antenna, and wherein the difference pseudorange residual between the pair of antennas of the set of antennas removes common components in the pseudorange measurements of the pair of antennas of the set of antennas, such that the difference pseudorange residual includes a difference between the antenna-specific clock offsets of the pair of antennas of the set of antennas subject to noise; and comparing the difference pseudorange residuals between different pairs of antennas of the set of antennas to produce the antenna-specific clock offset for each antenna of the set of antennas; correcting the pseudorange measurement of the set of antennas according to the antenna-specific clock offset for each antenna of the set of antennas to produce corrected pseudorange measurements; determining the time bias of the clock using the corrected pseudorange measurements; and outputting the determined time bias. ### Claim 18 (depends on 17) The medium of claim 17, wherein the method comprising: determining a probability of a first antenna-specific clock offset of the first antenna with respect to probabilities of antenna-specific clock offsets of a plurality of antennas of the set of antennas other than the first antenna as a marginal distribution of the first antenna-specific clock offset of the first antenna conditioned on a joint posteriori distribution of the antenna-specific clock offsets of all antennas of the set of antennas. ### Claim 19 (depends on 17) The medium of claim 17, wherein the antenna-specific clock offsets of all antennas of the set of antennas are determined based on corresponding marginal distributions of the antenna-specific clock offsets. ### Claim 20 (depends on 19) The medium of claim 19, wherein a second antenna-specific clock offset for a second antenna k of the set of antennas is determined using a belief propagation as a probabilistic belief bt(αk) at time t on a marginal distribution of an antenna-specific clock offset αk using a product of a prior distribution of the probabilistic belief bt(αk) and incoming belief messages from all antennas of the set of antennas except the second antenna k.