WO2018124250A1 - 測位装置、測位方法及びプログラム - Google Patents

測位装置、測位方法及びプログラム Download PDF

Info

Publication number
WO2018124250A1
WO2018124250A1 PCT/JP2017/047138 JP2017047138W WO2018124250A1 WO 2018124250 A1 WO2018124250 A1 WO 2018124250A1 JP 2017047138 W JP2017047138 W JP 2017047138W WO 2018124250 A1 WO2018124250 A1 WO 2018124250A1
Authority
WO
WIPO (PCT)
Prior art keywords
positioning
reception
receiving
antenna
signal
Prior art date
Application number
PCT/JP2017/047138
Other languages
English (en)
French (fr)
Inventor
俊雄 阿部
洋二 毛利
Original Assignee
スペースリンク株式会社
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by スペースリンク株式会社 filed Critical スペースリンク株式会社
Priority to JP2018559619A priority Critical patent/JPWO2018124250A1/ja
Publication of WO2018124250A1 publication Critical patent/WO2018124250A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
    • G01C21/28—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • 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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/22—Multipath-related issues
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • 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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • 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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • 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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/36—Constructional details or hardware or software details of the signal processing chain relating to the receiver frond end

Definitions

  • the present invention relates to a positioning device, a positioning method, and a program.
  • positioning by a GPS satellite in an object such as a flying object is generally performed by installing a GPS antenna above the object to be measured.
  • the antenna In the case of one antenna, the antenna is inclined due to the movement, inclination, etc. of the flying object. Visibility may be reduced. Then, there is a possibility that the positioning signal of the GPS satellite whose visibility is lowered cannot be received, the number of receiving satellites is reduced, and the positioning cannot be performed.
  • a positioning system has been developed in which a plurality of antennas and a receiver are arranged above an object to be positioned, and a positioning calculation of the receiver is performed using a radio signal received by each antenna.
  • the technique regarding the flying object provided with the positioning function by GPS is described in patent document 1, for example.
  • An object of the present invention is to perform positioning using a positioning satellite more appropriately.
  • positioning using a positioning satellite can be performed more appropriately.
  • FIG. 1 is a block diagram illustrating a configuration of a multi-antenna positioning apparatus 100 according to the embodiment.
  • a GPS signal is received as a positioning signal and the influence of multipath and interference waves is reduced using two antennas.
  • the multi-antenna positioning device 100 is configured as a flying object (here, a drone) will be described.
  • the multi-antenna positioning apparatus 100 includes an antenna (first antenna 101-1 and second antenna 101-2), front end (first RF front end 102-1,..., Second RF front end 102-2), clock generation.
  • Device 121 a reception processing unit (channel 1 reception processing unit ch 1 to channel 16 reception processing unit ch 16), a pseudo-range availability determination unit 118, and a positioning calculation unit 119.
  • the reception processing unit includes a circuit having a similar configuration of a number of channels (here, channel 1 to channel 16) corresponding to the number of positioning satellites that can be captured simultaneously. In FIG. 1, only the channel 1 reception processing unit ch1 is shown in detail.
  • the channel 1 reception processing unit ch1 includes a first mixer 104-1, an accumulator 105-1, a correlation value calculation unit 107-1, a carrier tracking unit 109-1, a code tracking unit 110-1, and a second mixer 113.
  • a pseudo-range calculation unit 116-1 which are realized by a logic circuit such as ASIC (Application Specific Integrated Circuit) that performs a predetermined calculation, or software executed by a CPU (Central Processing Unit).
  • Software to be executed by the CPU includes a method in which the software is executed by being loaded in a ROM (Read Only Memory), a method in which the software is stored in an external storage device such as a hard disk device, and is loaded onto a RAM (Random Access Memory) at the time of execution.
  • the channel 2 reception processing unit ch2 and subsequent channels have the same configuration as the channel 1 reception processing unit ch1.
  • Each channel can be realized on the same computer by parallel processing by a logic circuit or time division processing. Note that at least 4 channels are required for the calculation of positioning, but in practice, more channels, for example, 12 channels or more are usually mounted.
  • the first antenna 101-1 and the second antenna 101-2 are installed on a flying body (drone) that requires positioning, and receive a positioning signal electromagnetic wave to generate a high-frequency signal.
  • the first antenna 101-1 and the second antenna 101-2 are installed at positions symmetrical with respect to the central axis.
  • FIG. 2 is a configuration diagram showing an arrangement example of antennas installed on the flying object in the first embodiment.
  • the first antenna 101-1 and the second antenna 101-2 are installed on the flying object at an angle of 180 ° from the central axis.
  • the positioning signal is obtained by performing spectrum spreading on a carrier wave using a spreading code of a pseudo random noise code.
  • a C / A (Coarse / Acquisition) code is used as the pseudo random noise code.
  • the positioning signal includes navigation data necessary for positioning calculation.
  • the first RF front end 102-1 and the second RF front end 102-2 are connected to the first antenna 101-1 and the second antenna 101-2, respectively.
  • a common clock generated by the clock generator 121 is input to the first RF front end 102-1 and the second RF front end 102-2 and received by the first antenna 101-1 and the second antenna 101-2, respectively.
  • the high-frequency signal which is the positioning signal, is converted into a digital signal by the first RF front end 102-1 and the second RF front end 102-2. By converting to a digital signal, the subsequent processing can be executed by digital signal processing by a logic circuit or a CPU of a computer.
  • the first RF front end 102-1 and the second RF front end 102-2 amplify the high-frequency signal, down-convert the frequency of the high-frequency signal, and the down-converted high-frequency signal is input from the clock generator 121.
  • the channel 1 reception processing unit CH1 is a circuit that calculates a pseudo distance from a digitized positioning signal.
  • the first mixer 104-1,..., And the first mixer 104-16 of each channel respectively convert the digitized positioning signals 103-1,. , ..., despreading by mixing with 114-16 and restoring positioning signals, but with logic circuits and logic that performs predetermined operations such as addition and multiplication realized on a computer by software is there.
  • the ch1 replica signal 114-1 may be a combination of the carrier replica signal 111-1 to be matched with the carrier wave of the positioning signal and the code replica signal 112-1 to be matched with the code included in the positioning signal.
  • the ch1 accumulator 105-1 is configured to output a correlation value by performing accumulation processing for a predetermined time on the digitized positioning signal 103-1 subjected to despreading.
  • Correlation value calculation result 106-1 is an amount including the phase difference from replica signal 114-1 in the information, and correlation value calculation unit 107-1 receives correlation value calculation result 106-1 and performs a predetermined calculation process. I do.
  • the carrier tracking unit 109-1 generates a carrier replica signal 111-1 to be matched with the carrier wave of the positioning signal based on the correlation value calculation result 108-1 inputted thereto, and the correlation value calculation result 108-1 In this configuration, the phase of the carrier replica signal 111-1 is changed so that the size becomes maximum.
  • the code tracking unit 110-1 Based on the correlation value calculation result 108-1 input thereto, the code tracking unit 110-1 generates a code replica signal 112-1 to be matched with the code included in the positioning signal, and the correlation value calculation result 108- In this configuration, the phase of the code replica signal 112-1 is changed so that the magnitude of 1 is maximized, and the code delay time 115-1 of the code included in the received positioning signal is output.
  • the pseudo distance calculation unit 116-1 is configured to calculate the pseudo distance 117-1 based on the code delay time 115-1.
  • the pseudo distance 117-1 is output from the reception processing unit of each channel.
  • the positioning calculation unit 119-1 is configured to perform positioning calculation based on the pseudo distances 117-1 of a plurality of channels at the receiver time. The other channels operate in the same manner as ch1.
  • Each of the plurality of positioning satellites generates and transmits a positioning signal obtained by modulating a carrier wave with a code of a corresponding channel.
  • the first antenna 101-1 and the second antenna 101-2 installed on the top of the flying object receive positioning signals. The strength and phase of the positioning signal received by each antenna varies.
  • the positioning signal received by the first antenna 101-1 is sent to the first RF front end 102-1, where it is amplified and down-converted to a frequency suitable for digital processing, and is generated by the clock generator 121.
  • a / D conversion is performed on the basis of a clock common to all the antenna systems to obtain a digitized positioning signal 103-1.
  • the positioning signal received by the second antenna 101-2 is also sent to the second RF front end 102-2 in the same manner, and becomes a digitized positioning signal 103-2.
  • the digitized positioning signals 103-1 and 103-2 are A / D converted based on a common clock to become a digitized positioning signal 103-2.
  • the digitized positioning signal 103-1 is sent to the channel 1 reception processing unit ch1,...,
  • the channel 8 reception processing unit ch8, and the digitized positioning signal 103-2 is sent to the channel 9 reception processing unit ch9,. 16 is sent to the reception processing unit ch16.
  • the same processing is performed based on the code included in the positioning signal of the positioning satellite received by each channel. Below, the process performed by channel 1 reception process part CH1 is demonstrated.
  • the first mixer 104-1 mixes the input digitized positioning signal 103-1 with the common replica signal 114-1 by multiplication by digital processing or the like.
  • the mixed signal is subjected to a cumulative process for a predetermined time by the correlation value calculation unit 107-1.
  • Correlation value calculation result 106-1 is output from correlation value calculation section 107-1.
  • An equation model of these correlation values is shown in Equation 1.
  • P is the signal strength
  • ⁇ is the code phase difference between the received positioning signal and the replica signal 114-1
  • ⁇ f and ⁇ are the carrier frequency differences between the received positioning signal and the replica signal 114-1, respectively.
  • the phase difference is shown in Equation 1.
  • R represents the autocorrelation of codes included in the receiver signal and the replica signal, and is a maximum value when both ⁇ and ⁇ f are zero.
  • a carrier replica signal 111-1 and a code replica signal 112-1 are generated so that R always has a maximum value.
  • the carrier tracking unit 109-1 generates the carrier replica signal 111-1 by changing the phase of the signal of the carrier waveform so that ⁇ f becomes zero.
  • the code tracking unit 110-1 generates the code replica signal 112-1 by changing the phase of the signal of the waveform of the code corresponding to the tracking positioning satellite so that ⁇ becomes zero.
  • the second mixer 113-1 mixes the carrier replica signal 111-1 and the code replica signal 112-1, and generates a replica signal 114-1.
  • the replica signal 114-1 is controlled so that R is maximized.
  • the correlation value calculation result 106-1 is calculated by a predetermined calculation in the correlation value calculation unit 107-1, and is output as the correlation value calculation result 108-1.
  • the correlation value calculation result 108-1 is input to the carrier tracking unit 109-1 and the code tracking unit 110-1.
  • the code tracking unit 110-1 calculates the code delay time 115-1 of the received code with respect to an arbitrary receiver time t, and calculates the code delay time 115-1 To 116-1.
  • the pseudo distance calculation unit 116-1 calculates the pseudo distance 117-1 at the receiver time t from the code delay time 115-1 and outputs it.
  • the process until the pseudo-position 117-1 is output with the digitized positioning signal 103-1 as an input is called a channel.
  • the reception processing unit of each channel tracks a signal from a different GPS satellite and generates a pseudo distance 117-1 corresponding to the satellite.
  • the plurality of pseudo distances 117-1 generated in synchronization with the receiver time t is input to the pseudo distance availability determination unit 118-1.
  • FIG. 3 is a flowchart of the pseudo distance determination process for two antennas in the first embodiment.
  • the difference is smaller than a predetermined value X that is considered to be unaffected by multipath and interference waves, it is input to the positioning calculation unit 119 (step S2).
  • the positioning calculation unit 119 performs a positioning calculation from these pseudo distances and outputs a positioning result 120.
  • step S3 the pseudorange availability determination unit 118-1 Are not used for positioning (step S3).
  • step S3 the pseudo distance determination process ends.
  • a specific calculation method of the positioning calculation a known technique can be used. For example, it is described in the following document. (“Fundamental GPS Basic Concept / Positioning Principle / Signal and Receiver", Original by Pratap Misra and Per Enge, Translated by the Japan Institute of Navigation GPS, pp. 159-163, Shoyo Bunko)
  • the positioning calculation unit obtains the pseudo distance 117-1 from the antenna 101-1, and the pseudo distance 117-2 from the antenna 101-2. If there is no influence of multipath and interference wave, the absolute value of the difference in the pseudo distance obtained by subtracting the pseudo distance 117-2 from the pseudo distance 117-1 as shown in the flowchart of FIG. . When the difference between the absolute values of the pseudoranges is smaller than the predetermined value X, the pseudorange of the satellite is adopted for the positioning calculation.
  • FIG. 4 is a configuration diagram showing an arrangement example of antennas installed on the flying object in the second embodiment when there are three antennas.
  • the first antenna 101-1, the second antenna 101-2, and the third antenna 101-3 are arranged at an angle of 120 ° from each other in the flying object.
  • FIG. 1 is used in the second embodiment.
  • FIG. 5 is a configuration diagram showing an arrangement example of the antennas installed on the flying object in the second embodiment when the number of antennas is four. Furthermore, FIG.
  • N is a natural number of 3 or more.
  • N is a natural number of 3 or more.
  • N is a natural number of 3 or more.
  • the first antenna 101-1, the second antenna 101-2,..., The Mth antenna 101-M (M is a natural number of 3 or more) are respectively the first front end 102-1 and the second front end 102. -2, ..., the M-th front end 102-M is connected in a pair.
  • the first mixer 104-1, the second mixer 104-2,..., The Nth mixer 104-N are connected to the subsequent stage, and the first correlation value calculation unit 107- is further connected to the subsequent stage. 1, the second correlation value calculation unit 107-2,..., And the Nth correlation value calculation unit 107-N are connected.
  • the configuration other than ch1 is the same as that of ch1.
  • the positioning calculation unit obtains the pseudo distance 117-1 from the antenna 101-1, and obtains the pseudo distance 117-2 from the antenna 101-2.
  • the antenna 101-2 to obtain a pseudorange 117-2
  • the antenna 101-3 to obtain a pseudorange 117-3.
  • the absolute value of the difference between the pseudo distance 117-1 and the pseudo distance 117-2 is compared with a predetermined value X (step S11), and the absolute value of the difference between the pseudo distance 117-1 and the pseudo distance 117-2 is the predetermined value X. If smaller, the pseudo distance 117-1 or the pseudo distance 117-2 is used for positioning (step S12). If the absolute value of the difference between the pseudo distance 117-1 and the pseudo distance 117-2 is equal to or greater than a predetermined value X, one or both of the pseudo distance 117-1 and the pseudo distance 117-2 may be affected by multipath or interference. Then, the absolute value of the difference between the pseudo distance 117-2 and the pseudo distance 117-3 is compared with a predetermined value X (step S13).
  • the pseudo distance 117-2 or the pseudo distance 117-3 is used for positioning (step S14). If the absolute value of the difference between the pseudo distance 117-2 and the pseudo distance 117-3 is equal to or greater than a predetermined value X, one or both of the pseudo distance 117-2 and the pseudo distance 117-3 may be affected by multipath or interference. Then, the absolute value of the difference between the pseudo distance 117-3 and the pseudo distance 117-1 is compared with a predetermined value X (step S15).
  • step S16 When the absolute value of the difference between the pseudo distance 117-3 and the pseudo distance 117-1 is smaller than the predetermined value X, the pseudo distance 117-3 or the pseudo distance 117-1 is used for positioning (step S16), but the pseudo distance 117 is used. If the absolute value of the difference between ⁇ 3 and the pseudorange 117-1 is equal to or greater than the predetermined value X, the pseudorange of this satellite is not used for the positioning calculation. After step S12, step S14, and step S16, the pseudo distance determination process ends.
  • FIG. 8 is a flowchart of pseudo distance determination processing when N antennas are mounted in the second embodiment.
  • the absolute value of the difference between the pseudo distance 117-1 and the pseudo distance 117-2 is compared with a predetermined value X (step S21), and the absolute value of the difference between the pseudo distance 117-1 and the pseudo distance 117-2 is the predetermined value X.
  • the pseudo distance 117-1 or the pseudo distance 117-2 is used for positioning (step S22). If the absolute value of the difference between the pseudo distance 117-1 and the pseudo distance 117-2 is equal to or greater than a predetermined value X, one or both of the pseudo distance 117-1 and the pseudo distance 117-2 may be affected by multipath or interference. Then, the absolute value of the difference between the pseudo distance 117-2 and the pseudo distance 117-3 is compared with a predetermined value X (step S23). If the absolute value of the difference between the pseudo distance 117-2 and the pseudo distance 117-3 is smaller than the predetermined value X, the pseudo distance 117-2 or the pseudo distance 117-3 is used for positioning (step S24).
  • the absolute value of the difference between the pseudo distance 117-2 and the pseudo distance 117-3 is equal to or greater than a predetermined value X, one or both of the pseudo distance 117-2 and the pseudo distance 117-3 may be affected by multipath or interference. Then, the absolute value of the difference between the pseudo distance 117-3 and the pseudo distance 117-4 is compared with a predetermined value X, and the same processing is executed (steps S25 and S26). Thereafter, the same processing is repeated up to the pseudo distance N (steps S27 and S28). However, when the absolute value of the difference between the pseudo distances received from adjacent antennas continues to the pseudo distance N, the absolute value of the difference between 117-1 and the pseudo distance 117-3 is a predetermined value.
  • the pseudo distance 117-3 or the pseudo distance 117-1 is used for positioning.
  • the absolute value of the difference between the pseudorange 117-3 and the pseudorange 117-1 is equal to or greater than the predetermined value X, the pseudorange of this satellite is not used for positioning calculation.
  • the same satellite is received by a plurality of antennas, and the pseudo distance calculated based on the received signal obtained from the receiving antenna affected by the multipath and interference wave is calculated as many as the number of antennas. If pseudoranges calculated from multiple antennas within the same satellite do not fall within a certain error, it is assumed that the pseudorange of the antenna is affected by multipath or interference, and that pseudorange is used for positioning calculation. By not putting in, the effect that a highly accurate and reliable positioning result can be obtained is acquired.
  • the present invention can be used for a flying object other than a drone, an agricultural construction machine, or a positioning sensor for automatic driving of a car.
  • a gyro sensor may be provided. Then, when the attitude of the multi-antenna positioning apparatus 100 is changed, a receiving antenna located on the upper side (upper side) is selected from the detection result of the gyro sensor, and calculation is performed based on the received signal obtained from the selected receiving antenna. Positioning can be performed using the pseudorange. In this case, the multi-antenna positioning apparatus 100 can perform positioning using the received signal of the receiving antenna estimated to have better quality. Therefore, positioning using the positioning satellite can be performed more appropriately.
  • a reception antenna with a better reception state may be selected based on the reception state of each reception antenna. That is, from the reception state of each reception antenna, it is determined which reception signal of which reception antenna is of higher quality, and using the pseudo distance calculated based on the reception signal obtained from the reception antenna of higher quality Positioning can be performed.
  • Various parameters in satellite positioning can be used as an indicator of the reception state used when judging the quality of the received signal. For example, the carrier phase (degree of coincidence with the carrier replica signal, etc.), signal strength A change (fluctuation state based on a history of signal intensity in a predetermined period) or the like can be used. With such a configuration, positioning using a positioning satellite can be performed using a reception signal of a receiving antenna with higher reliability without adding a special configuration.
  • a receiving antenna used for positioning is selected from three or more receiving antennas provided in the multi-antenna positioning apparatus 100, two receiving antennas can be sequentially compared as in the above-described embodiment. Priority can be given to comparison between receiving antennas that are spatially separated. For example, when the number of receiving antennas provided in the multi-antenna positioning apparatus 100 is four, two receiving antennas at diagonal positions are preferentially compared rather than two adjacent receiving antennas. Thereby, when the attitude
  • a receiving antenna to be used for positioning from three or more receiving antennas provided in the multi-antenna positioning apparatus 100, a comparison of three or more receiving antennas is performed at the same time instead of comparing two receiving antennas each. It may be done. Then, based on the dispersion state of the pseudo distance calculated based on the reception signals obtained from these three or more reception antennas, a reception antenna estimated to be suitable for positioning is selected, and the reception antennas are selected.
  • the pseudo-range at may be used for positioning.
  • the pseudo distances of the three receiving antennas are within the preset threshold range Y, and the pseudo distances of the remaining one receiving antenna are out of the threshold range Y. If the value is an arbitrary value, any one of the reception antennas within the threshold range Y may be selected, and the pseudo distance at the reception antenna may be used for positioning. Thereby, it is possible to perform positioning by comparing a plurality of receiving antennas provided in the multi-antenna positioning apparatus 100 in a lump and selecting a receiving antenna more quickly.
  • the present invention can be implemented by appropriately combining the above-described embodiment and each modification.
  • the processing in the above-described embodiment can be executed by either hardware or software. That is, it is only necessary that the multiple antenna positioning device 100 has a function capable of executing the above-described processing, and what functional configuration and hardware configuration are used to realize this function is not limited to the above-described example.
  • a program constituting the software is installed on a computer from a network or a storage medium.
  • the storage medium for storing the program includes a removable medium distributed separately from the apparatus main body, or a storage medium incorporated in the apparatus main body in advance.
  • the removable medium is composed of, for example, a magnetic disk, an optical disk, a magneto-optical disk, or the like.
  • the optical disk is composed of, for example, a CD-ROM (Compact Disk-Read Only Memory), a DVD (Digital Versatile Disk), a Blu-ray Disc (registered trademark), and the like.
  • the magneto-optical disk is constituted by an MD (Mini-Disk) or the like.
  • the storage medium incorporated in advance in the apparatus main body is constituted by, for example, a ROM or a hard disk in which a program is stored.

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Automation & Control Theory (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

測位衛星を用いた測位をより適切に行うこと。　測位装置は、複数のアンテナと、１または複数の受信回路と、前記アンテナ及び前記受信回路を含む受信経路を介して複数の衛星からの測位信号を受信する受信部と、同一の前記衛星からの測位信号について、少なくとも一部が異なる前記受信経路を介して受信された複数の受信信号に基づいて、測位に用いる前記受信信号を判断する判断部と、前記判断部によって測位に用いると判断された前記受信信号に基づいて、現在位置を算出する位置算出部と、を備える。

Description

測位装置、測位方法及びプログラム
　本発明は、測位装置、測位方法及びプログラムに関する。
　従来、飛翔体等の物体におけるＧＰＳ衛星による測位は、一般に測位を行う物体の上部にＧＰＳアンテナを設置するが、一つのアンテナの場合、飛翔体の移動、傾き等により、アンテナが傾き、ＧＰＳ衛星の可視性が低下する場合がある。すると、可視性が低下したＧＰＳ衛星の測位信号を受信できず、受信衛星が減り、測位を行えない可能性がある。
　この問題を解決するため、測位を行う物体の上部に複数のアンテナと受信機を配置し、各アンテナで受信された無線信号により受信機の測位計算を行う測位システムが開発されている。
　なお、ＧＰＳによる測位機能を備えた飛翔体に関する技術は、例えば、特許文献１に記載されている。
特開２０１７－１２４６８９号公報
　しかしながら、アンテナを複数にした場合、ビルによる電波の反射波によるマルチパスや妨害波が入力されることが考えられ、その結果、測位結果の誤差が大きくなり間違った測位結果を招くことが考えられる。
　本発明の課題は、測位衛星を用いた測位をより適切に行うことである。
　上記課題を解決するため、本発明の一態様に係る測位装置は、
　複数のアンテナと、
　１または複数の受信回路と、
　前記アンテナ及び前記受信回路を含む受信経路を介して複数の衛星からの測位信号を受信する受信部と、
　同一の前記衛星からの測位信号について、少なくとも一部が異なる前記受信経路を介して受信された複数の受信信号に基づいて、測位に用いる前記受信信号を判断する判断部と、
　前記判断部によって測位に用いると判断された前記受信信号に基づいて、現在位置を算出する位置算出部と、
　を備えることを特徴とする。
　本発明によれば、測位衛星を用いた測位をより適切に行うことができる。
実施形態に係る複数アンテナ測位装置１００の構成を示すブロック図である。 第１の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。 第１の実施形態における２アンテナ時の擬似距離判断処理のフローチャート図である。 アンテナを３つとした場合の第２の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。 アンテナを４つとした場合の第２の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。 アンテナをＮ個とした場合の第２の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。 第２の実施形態における３アンテナ時の擬似距離判断処理のフローチャート図である。 第２の実施形態におけるＮ個のアンテナ搭載時の擬似距離判断処理のフローチャート図である。
［第１の実施形態］
　以下、本発明の実施形態について、図面を参照して説明する。
　図１は、実施形態に係る複数アンテナ測位装置１００の構成を示すブロック図である。
　本実施形態においては、測位信号としてＧＰＳ信号を受信し、２個のアンテナを使用してマルチパス、妨害波の影響を低減する例を示す。
　なお、以下の実施形態においては、複数アンテナ測位装置１００を飛翔体（ここではドローンとする）として構成した例について説明する。
　初めに、複数アンテナ測位装置１００の構成について説明する。
　複数アンテナ測位装置１００は、アンテナ（第１アンテナ１０１－１、第２アンテナ１０１－２）、フロントエンド（第１ＲＦフロントエンド１０２－１、・・・、第２ＲＦフロントエンド１０２－２）、クロック発生器１２１、受信処理部（チャンネル１受信処理部ｃｈ１～チャンネル１６受信処理部ｃｈ１６）、擬似距離可否判断部１１８、及び測位演算部１１９から構成される。
　受信処理部は、同時に捕捉可能な測位衛星の数に対応した数のチャンネル（ここではチャンネル１からチャンネル１６）の同様の構成を有する回路を含む。図１においては、これらのうち、チャンネル１受信処理部ｃｈ１のみ詳細な構成を示す。
　チャンネル１受信処理部ｃｈ１は、第１混合器１０４－１、累算器１０５－１、相関値演算部１０７－１、搬送波追尾部１０９－１、コード追尾部１１０－１、第２混合器１１３－１、擬似距離演算部１１６－１から構成され、これらは所定の演算を行うＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）等のロジック回路や、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）で実行されるソフトウェアによって実現される。ＣＰＵにより実行させるソフトウェアは、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）に搭載して実行させる方法、ハードディスク装置等の外部記憶装置に記憶し、実行時にＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）上にローディングする方法等もある。チャンネル２受信処理部ｃｈ２以降も、チャンネル１受信処理部ｃｈ１と同様の構成を有している。なお、それぞれのチャンネルは、ロジック回路による並列処理や、時分割処理により同一のコンピュータ上で実現することができる。なお、測位の計算には、最低でも４チャンネルが必要であるが、実際には、それより多くのチャンネル、例えば１２チャンネル以上が実装されることが通常である。
　第１アンテナ１０１－１、第２アンテナ１０１－２は、測位が必要とされる飛翔体（ドローン）の機体に設置され、測位信号の電磁波を受信して高周波信号を発生する。飛翔体には、第１アンテナ１０１－１、第２アンテナ１０１－２を、その中心軸に対して対称な位置に設置する。
　図２は、第１の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。
　飛翔体に第１アンテナ１０１－１、第２アンテナ１０１－２を中心軸上から１８０°の角度を空けて設置する。測位信号は、搬送波を擬似ランダム雑音コードの拡散符号によってスペクトラム拡散を行ったものである。擬似ランダム雑音コードとしては、Ｃ／Ａ（Ｃｏａｒｓｅ／Ａｃｑｕｉｓｉｔｉｏｎ）コードを使用する。また、測位信号には、測位演算に必要な航法データが含まれている。
　第１アンテナ１０１－１、第２アンテナ１０１－２には、それぞれ、第１ＲＦフロントエンド１０２－１、第２ＲＦフロントエンド１０２－２が接続される。第１ＲＦフロントエンド１０２－１、第２ＲＦフロントエンド１０２－２には、クロック発生器１２１で発生された共通のクロックが入力され、それぞれ、第１アンテナ１０１－１、第２アンテナ１０１－２で受信された測位信号である高周波信号を次の第１ＲＦフロントエンド１０２－１、第２ＲＦフロントエンド１０２－２でデジタル信号に変換する。デジタル信号に変換することにより、以降の処理を、ロジック回路やコンピュータのＣＰＵによるデジタル信号処理で実行可能となる。
　第１ＲＦフロントエンド１０２－１、第２ＲＦフロントエンド１０２－２は、高周波信号を増幅し、前記の高周波信号の周波数をダウンコンバートし、ダウンコンバートされた高周波信号を、クロック発生器１２１から入力された共通のクロックを基準としてＡ／Ｄ（Ａｎａｌｏｇ／Ｄｉｇｉｔａｌ）変換し、その後、サンプリングして量子化することによって、それぞれ、デジタル化された測位信号１０３－１、１０３－２、・・・、１０３－１６を生成する。それらのデジタル化測位信号１０３－１及び１０３－２は、共通のクロックによってサンプリングされ、観測時刻は同期している。
　チャンネル１受信処理部ＣＨ１は、デジタル化された測位信号から擬似距離を算出する回路である。各ｃｈの第１混合器１０４－１、・・・、第１混合器１０４－１６は、それぞれ、デジタル化測位信号１０３－１、・・・、１０３－１６を各ｃｈのレプリカ信号１１４－１、・・・、１１４－１６と混合することによって逆拡散を行い、測位信号の復元を行うが、ロジック回路や、ソフトウェアによってコンピュータ上に実現する、加算、乗算等の所定の演算を行うロジックである。
　ｃｈ１のレプリカ信号１１４－１は、測位信号の搬送波と一致させるべき搬送波レプリカ信号１１１－１と、測位信号に含まれるコードと一致させるべきコードレプリカ信号１１２－１を合成したものであり得る。ｃｈ１の累算器１０５－１は、逆拡散を行ったデジタル化測位信号１０３－１に対して一定時間の累積処理を行うことによって、相関値を出力する構成である。相関値演算結果１０６－１は、レプリカ信号１１４－１との位相差を情報に含む量であり、相関値演算部１０７－１は、相関値演算結果１０６－１を入力し、所定の演算処理を行う。
　搬送波追尾部１０９－１は、そこに入力される相関値演算結果１０８－１に基づいて、測位信号の搬送波と一致させるべき搬送波レプリカ信号１１１－１を生成し、相関値演算結果１０８－１の大きさが極大になるように搬送波レプリカ信号１１１－１の位相を変化させる構成である。
　コード追尾部１１０－１は、そこに入力される相関値演算結果１０８－１に基づいて、測位信号に含まれるコードと一致させるべきコードレプリカ信号１１２－１を生成し、相関値演算結果１０８－１の大きさが極大になるようにコードレプリカ信号１１２－１の位相を変化させると共に、受信した測位信号に含まれるコードのコード遅延時間１１５－１を出力する構成である。
　擬似距離演算部１１６－１は、コード遅延時間１１５－１に基づいて、擬似距離１１７－１を計算する構成である。擬似距離１１７－１は、それぞれのチャンネルの受信処理部から出力される。
　測位演算部１１９－１は、受信機時刻における複数のチャンネルの擬似距離１１７－１に基づき、測位演算を行う構成である。他のｃｈもｃｈ１と同様の動作となる。
　次に、本実施形態に係る複数アンテナ測位装置１００の動作について説明する。
　複数の測位衛星のそれぞれは、それに対応するチャンネルのコードで搬送波を変調した測位信号を生成して送信している。飛翔体の上部に設置された第１アンテナ１０１－１及び第２アンテナ１０１－２は、測位信号を受信する。それぞれのアンテナが受信する測位信号の強度と位相は変動する。第１アンテナ１０１－１で受信された測位信号は、第１ＲＦフロントエンド１０２－１に送られ、そこで、増幅された後にデジタル処理に適した周波数までダウンコンバートされ、クロック発生器１２１で発生させられたすべてのアンテナ系に共通のクロックに基づいてＡ／Ｄ変換されてデジタル化測位信号１０３－１となる。第２アンテナ１０１－２で受信された測位信号も、同様に第２ＲＦフロントエンド１０２－２に送られ、デジタル化測位信号１０３－２となる。デジタル化測位信号１０３－１及び１０３－２は、共通のクロックに基づいてＡ／Ｄ変換されてデジタル化測位信号１０３－２となる。デジタル化測位信号１０３－１は、チャンネル１受信処理部ｃｈ１、・・・、チャンネル８受信処理部ｃｈ８に送られ、デジタル化測位信号１０３－２はチャンネル９受信処理部ｃｈ９、・・・、チャンネル１６受信処理部ｃｈ１６に送られる。それぞれのチャンネルの受信処理部では、それぞれが受信する測位衛星の測位信号に含まれるコードに基づいて、同様の処理が実施される。以下では、チャンネル１受信処理部ＣＨ１で行われる処理について説明する。
　第１混合器１０４－１は、入力されたデジタル化測位信号１０３－１に共通のレプリカ信号１１４－１をデジタル処理による乗算等によって混合する。混合された信号は、相関値演算部１０７－１によって一定時間の累積処理が行われる。相関値演算部１０７－１から、相関値演算結果１０６－１が出力される。それらの相関値の数式モデルを式１に示す。式１において、Pは信号強度、Δτは受信した測位信号とレプリカ信号１１４－１の間のコード位相差、Δｆとφはそれぞれ受信した測位信号とレプリカ信号１１４－１の間の搬送波の周波数差と位相差である。
Figure JPOXMLDOC01-appb-M000001
　なお、式１において、Ｒは受信機信号とレプリカ信号に含まれるコードの自己相関を表し、ΔτとΔｆとが共にゼロであるとき最大値となる。ＧＰＳの測位信号では、Ｒが常に最大値となるような搬送波レプリカ信号１１１－１とコードレプリカ信号１１２－１とが生成される。このために、搬送波追尾部１０９－１は、搬送波の波形の信号の位相をΔｆがゼロになるように変化させることによって搬送波レプリカ信号１１１－１を生成する。そして、コード追尾部１１０－１は、追尾している測位衛星に対応するコードの波形の信号の位相をΔτがゼロになるように変化させることによってコードレプリカ信号１１２－１を生成する。次に、第２混合器１１３－１が、搬送波レプリカ信号１１１－１とコードレプリカ信号１１２－１を混合してレプリカ信号１１４－１を生成する。このレプリカ信号１１４－１は、Ｒが最大になるように制御されているものである。
　式（１）において、Ｒはその最大値で正規化されているとする。このとき、ＧＰＳの測位信号が追尾されている状態、つまりΔτとΔｆが共にゼロであるとき、相関値演算結果１０６－１は、式（２）で表される。式（２）において、ＰＡとφＡはそれぞれ相関値演算結果１０６－１の信号強度と位相差である。
Figure JPOXMLDOC01-appb-M000002
　相関値演算結果１０６－１は、相関値演算部１０７－１において所定の演算により計算され、相関値演算結果１０８－１として出力される。
　相関値演算結果１０８－１は、搬送波追尾部１０９－１およびコード追尾部１１０－１に入力される。コード追尾部１１０－１は、受信された測位信号に含まれるコードを追尾しているとき、任意の受信機時刻ｔに対する受信コードのコード遅延時間１１５－１を算出し、それを擬似距離演算部１１６－１に出力する。擬似距離演算部１１６－１は、コード遅延時間１１５－１から受信機時刻ｔにおける擬似距離１１７－１を算出して出力する。
　デジタル化測位信号１０３－１を入力として擬似距離１１７－１を出力するまでの処理は、チャンネルと呼ばれる。各チャンネルの受信処理部では、それぞれ異なるＧＰＳ衛星からの信号を追尾し、その衛星に対応した擬似距離１１７－１の生成を行う。受信機時刻tに同期して生成された複数の擬似距離１１７－１は、擬似距離可否判断部１１８－１に入力される。
　図３は、第１の実施形態における２アンテナ時の擬似距離判断処理のフローチャート図である。
　擬似距離可否判断部１１８－１において、同一衛星で複数アンテナにより得られた擬似距離でそれぞれの擬似距離１１７－１（＝Ａ）及び擬似距離１１７－２（＝Ｂ）の差の絶対値を所定の値Ｘと比較し（ステップＳ１）、その差がマルチパス、妨害波の影響を受けていないと考えられる所定の値Ｘより小さい場合に測位演算部１１９に入力される（ステップＳ２）。測位演算部１１９では、これら擬似距離より測位演算が行われ、測位結果１２０が出力される。なお、同一衛星で複数アンテナにより得られた擬似距離の差がマルチパス、妨害波の影響を受けていないと考えられる所定値以上である場合には、擬似距離可否判断部１１８－１は、これらの擬似距離を測位に使用しない（ステップＳ３）。ステップＳ２及びステップＳ３の後、擬似距離判断処理は終了となる。
　測位演算の具体的演算方法に関しては、公知の技術を用いることができる。例えば次の文献に記載がある。（「精説GPS　基本概念・測位原理・信号と受信機」，Pratap Misra and Per Enge原著，日本航海学会GPS研究会訳，pp. 159-163，正陽文庫）
　測位演算部でアンテナ１０１－１から擬似距離１１７－１が得られ、アンテナ１０１－２から擬似距離１１７－２が得られる。マルチパス、妨害波の影響を受けていなければ、図３のフローチャートのように擬似距離１１７－１から擬似距離１１７－２を引いた擬似距離の差の絶対値はアンテナの位置から近い値となる。その擬似距離の絶対値の差が所定値Ｘより小さい場合、この衛星の擬似距離を測位演算に採用する。
［第２の実施形態］
　第１の実施形態では、２つのアンテナで受信される測位信号を用いているが、３つ以上のアンテナが利用できる場合、これらの３つ以上のアンテナで受信される測位信号を用いることができる。
　図４は、アンテナを３つとした場合の第２の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。第１アンテナ１０１－１、第２アンテナ１０１－２、第３アンテナ１０１－３は、飛翔体においてお互いに１２０°の角度を空けて配置される。第１の実施形態と同様に第２の実施形態でも図１を使用する。
　また、図５は、アンテナを４つとした場合の第２の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。
　さらに、図６は、アンテナをＮ個（Ｎは３以上の自然数）とした場合の第２の実施形態における飛翔体に設置したアンテナの配置例を示す構成図である。
　第２の実施形態において、アンテナが一般的にＮ個設置されている場合、複数アンテナ測位装置１００のアンテナから相関値演算部までの構成は、Ｎ組設置される。
　第１アンテナ１０１－１、第２アンテナ１０１－２、・・・、第Ｍアンテナ１０１－Ｍ（Ｍは３以上の自然数）には、それぞれ、第１フロントエンド１０２－１、第２フロントエンド１０２－２、・・・、第Ｍフロントエンド１０２－Ｍが一対に接続される。それらの後段に、第１混合器１０４－１、第２混合器１０４－２、・・・、第Ｎ混合器１０４－Ｎが接続され、それらのさらに後段に、第１相関値演算部１０７－１、第２相関値演算部１０７－２、・・・、第Ｎ相関値演算部１０７－Ｎが接続される。ｃｈ１以外についても、ｃｈ１の同様の構成である。
　図１では測位演算部でアンテナ１０１－１から擬似距離１１７－１が得られ、アンテナ１０１－２から擬似距離１１７－２が得られるが、同様に、図５では、３アンテナの場合、アンテナ１０１－１から擬似距離１１７－１が得られ、アンテナ１０１－２から擬似距離１１７－２が得られ、アンテナ１０１－３から擬似距離１１７－３が得られる。
　図７は、第２の実施形態における３アンテナ時の擬似距離判断処理のフローチャート図である。
　３アンテナの場合の各アンテナについて、図７のように同一衛星で各アンテナから受信した擬似距離１１７－１（＝Ａ）、擬似距離１１７－２（＝Ｂ）、擬似距離１１７－３（＝Ｃ）で測位に使用できる擬似距離を判断する。擬似距離１１７－１と擬似距離１１７－２の差の絶対値を所定の値Ｘと比較し（ステップＳ１１）、擬似距離１１７－１と擬似距離１１７－２の差の絶対値が所定の値Ｘより小さい場合、擬似距離１１７－１または擬似距離１１７－２を測位に利用する（ステップＳ１２）。もし擬似距離１１７－１と擬似距離１１７－２の差の絶対値が所定の値Ｘ以上の場合、擬似距離１１７－１、擬似距離１１７－２の片方または両方がマルチパスまたは妨害波の影響を受けていると判断し、次に擬似距離１１７－２と擬似距離１１７－３の差の絶対値を所定の値Ｘと比較する（ステップＳ１３）。擬似距離１１７－２と擬似距離１１７－３の差の絶対値が所定の値Ｘより小さい場合、擬似距離１１７－２または擬似距離１１７－３を測位に利用する（ステップＳ１４）。もし擬似距離１１７－２と擬似距離１１７－３の差の絶対値が所定の値Ｘ以上の場合、擬似距離１１７－２、擬似距離１１７－３の片方または両方がマルチパスまたは妨害波の影響を受けていると判断し、次に擬似距離１１７－３と擬似距離１１７－１の差の絶対値を所定の値Ｘと比較する（ステップＳ１５）。擬似距離１１７－３と擬似距離１１７－１の差の絶対値が所定の値Ｘより小さい場合、擬似距離１１７－３または擬似距離１１７－１を測位に利用する（ステップＳ１６）が、擬似距離１１７－３と擬似距離１１７－１の差の絶対値が所定の値Ｘ以上の場合、この衛星の擬似距離を測位計算に利用しない。ステップＳ１２、ステップＳ１４及びステップＳ１６の後、擬似距離判断処理は終了となる。
（アンテナＮ個の場合の擬似距離判断処理）
　図８は、第２の実施形態におけるＮ個のアンテナ搭載時の擬似距離判断処理のフローチャート図である。
　図８のフローチャートのように、同一衛星で各アンテナから受信した擬似距離１１７－１（＝Ａ１）、擬似距離１１７－２（＝Ａ２）、・・・、擬似距離１１７－Ｎ（＝ＡＮ）で測位に使用できる擬似距離を判断する。擬似距離１１７－１と擬似距離１１７－２の差の絶対値を所定の値Ｘと比較し（ステップＳ２１）、擬似距離１１７－１と擬似距離１１７－２の差の絶対値が所定の値Ｘより小さい場合、擬似距離１１７－１または擬似距離１１７－２を測位に利用する（ステップＳ２２）。もし擬似距離１１７－１と擬似距離１１７－２の差の絶対値が所定の値Ｘ以上の場合、擬似距離１１７－１、擬似距離１１７－２の片方または両方がマルチパスまたは妨害波の影響を受けていると判断し、次に擬似距離１１７－２と擬似距離１１７－３の差の絶対値を所定の値Ｘと比較する（ステップＳ２３）。擬似距離１１７－２と擬似距離１１７－３の差の絶対値が所定の値Ｘより小さい場合、擬似距離１１７－２または擬似距離１１７－３を測位に利用する（ステップＳ２４）。もし擬似距離１１７－２と擬似距離１１７－３の差の絶対値が所定の値Ｘ以上の場合、擬似距離１１７－２、擬似距離１１７－３の片方または両方がマルチパスまたは妨害波の影響を受けていると判断し、次に擬似距離１１７－３と擬似距離１１７－４の差の絶対値を所定の値Ｘと比較し、同様の処理を実行する（ステップＳ２５，Ｓ２６）。以降、疑似距離Ｎまで同様の処理を繰り返す（ステップＳ２７，Ｓ２８）。ただし、隣接する番号のアンテナから受信した擬似距離の差の絶対値が所定値Ｘより大きい場合が擬似距離Ｎまで続く場合、１１７－１と擬似距離１１７－３の差の絶対値が所定の値Ｘより小さい場合、擬似距離１１７－３または擬似距離１１７－１を測位に利用する。なお、擬似距離１１７－３と擬似距離１１７－１の差の絶対値が所定の値Ｘ以上の場合、この衛星の擬似距離を測位計算に利用しない。ステップＳ２２、ステップＳ２４、ステップＳ２６及びステップＳ２８の後、擬似距離判断処理は終了となる。
　本発明によれば、同一衛星を複数アンテナにより受信し、マルチパス、妨害波の影響を受けた受信アンテナから得られる受信信号を基に算出された擬似距離をアンテナの数ほど計算する。もし同一衛星で複数のアンテナから算出した擬似距離がある誤差以内に入っていなければ、そのアンテナの擬似距離はマルチパス、または妨害波の影響を受けていると考え、その擬似距離を測位演算に入れないことにより、高精度で信頼性の高い測位結果を得られることができるという効果が得られる。
　応用例として、ドローン以外の飛翔体、農業建機、または自動車の自動運転の測位センサー等に利用可能である。
［変形例１］
　複数アンテナ測位装置１００をドローン等の飛翔体、その他、姿勢変化を行う移動体として構成する場合、ジャイロセンサを備えることとしてもよい。
　そして、複数アンテナ測位装置１００が姿勢変化した場合に、ジャイロセンサの検出結果から、より上側（上空側）に位置する受信アンテナを選択し、選択された受信アンテナから得られる受信信号を基に算出された擬似距離を用いて測位を行うことができる。
　この場合、複数アンテナ測位装置１００において、品質がより良好であると推定される受信アンテナの受信信号を用いて、測位を行うことができる。
　したがって、測位衛星を用いた測位をより適切に行うことが可能となる。
［変形例２］
　複数アンテナ測位装置１００において、測位に用いる受信アンテナを選択する場合、受信アンテナ毎の受信状態の良否に基づいて、より受信状態が良好な受信アンテナを選択することとしてもよい。
　即ち、各受信アンテナの受信状態から、いずれの受信アンテナの受信信号がより高品質であるかを判断し、より高品質な受信アンテナから得られる受信信号を基に算出された擬似距離を用いて測位を行うことができる。受信信号の品質を判断する場合に用いられる受信状態の指標としては、衛星測位における各種パラメータを用いることが可能であるが、例えば、搬送波位相（搬送波レプリカ信号との一致度合等）、信号強度の変化（所定期間の信号強度の履歴に基づく変動状態）等を用いることができる。
　このような構成により、特別な構成を追加することなく、より信頼性の高い受信アンテナの受信信号を用いて、測位衛星を用いた測位を行うことができる。
［変形例３］
　複数アンテナ測位装置１００に備えられた３個以上の受信アンテナから測位に用いる受信アンテナを選択する場合、上述の実施形態のように、２個の受信アンテナを順次比較することができ、この場合に、空間的に離れて配置された受信アンテナ同士の比較を優先することができる。
　例えば、複数アンテナ測位装置１００に備えられている受信アンテナが４個の場合、隣接する２個の受信アンテナではなく、対角の位置にある２個の受信アンテナを優先して比較する。
　これにより、複数アンテナ測位装置１００の姿勢が水平から傾斜した場合等に、衛星信号の受信状態が異なると推定される２個の受信アンテナの受信信号を優先して比較することができる。
　そのため、衛星信号の受信状態が異なる状況が発生した場合に、より早い段階で、測位に用いることが適切でない受信アンテナが存在することを検出できる。
　また、複数アンテナ測位装置１００に備えられた３個以上の受信アンテナから測位に用いる受信アンテナを選択する場合、２個ずつの受信アンテナの比較に代えて、３個以上の受信アンテナの比較を同時に行うこととしてもよい。そして、これら３個以上の受信アンテナから得られる受信信号を基に算出された擬似距離の分散状態に基づいて、測位に用いることが適切であると推定される受信アンテナを選択し、その受信アンテナにおける擬似距離を測位に用いることとしてもよい。
　例えば、４個の受信アンテナの擬似距離のうち、３個の受信アンテナにおける擬似距離が予め設定された閾値範囲Ｙに収まっており、残りの１個の受信アンテナにおける擬似距離が閾値範囲Ｙから外れた値である場合、閾値範囲Ｙに収まっているいずれかの受信アンテナを選択し、その受信アンテナにおける擬似距離を測位に用いることとしてもよい。
　これにより、複数アンテナ測位装置１００に備えられた複数の受信アンテナを一括して比較し、より速やかに受信アンテナを選択して、測位を行うことができる。
　なお、本発明は、上述の実施形態に限定されるものではなく、本発明の目的を達成できる範囲での変形、改良等は本発明に含まれるものである。
　例えば、上記実施形態及び各変形例を適宜組み合わせて、本発明を実施することが可能である。
　上述の実施形態における処理は、ハードウェア及びソフトウェアのいずれにより実行させることも可能である。
　即ち、上述の処理を実行できる機能が複数アンテナ測位装置１００に備えられていればよく、この機能を実現するためにどのような機能構成及びハードウェア構成とするかは上述の例に限定されない。
　上述の処理をソフトウェアにより実行させる場合には、そのソフトウェアを構成するプログラムが、コンピュータにネットワークや記憶媒体からインストールされる。
　プログラムを記憶する記憶媒体は、装置本体とは別に配布されるリムーバブルメディア、あるいは、装置本体に予め組み込まれた記憶媒体等で構成される。リムーバブルメディアは、例えば、磁気ディスク、光ディスク、または光磁気ディスク等により構成される。光ディスクは、例えば、ＣＤ－ＲＯＭ（Ｃｏｍｐａｃｔ　Ｄｉｓｋ－Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ），ＤＶＤ（Ｄｉｇｉｔａｌ　Ｖｅｒｓａｔｉｌｅ　Ｄｉｓｋ），Ｂｌｕ－ｒａｙ　Ｄｉｓｃ（登録商標）等により構成される。光磁気ディスクは、ＭＤ（Ｍｉｎｉ－Ｄｉｓｋ）等により構成される。また、装置本体に予め組み込まれた記憶媒体は、例えば、プログラムが記憶されているＲＯＭやハードディスク等で構成される。
　以上、本発明のいくつかの実施形態について説明したが、これらの実施形態は、例示に過ぎず、本発明の技術的範囲を限定するものではない。本発明はその他の様々な実施形態を取ることが可能であり、さらに、本発明の要旨を逸脱しない範囲で、省略や置換等種々の変更を行うことができる。これら実施形態やその変形は、本明細書等に記載された発明の範囲や要旨に含まれると共に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
ｃｈ１　チャンネル１受信処理部
ｃｈ２　チャンネル２受信処理部
ｃｈ３　チャンネル３受信処理部
ｃｈ４　チャンネル４受信処理部
ｃｈ５　チャンネル５受信処理部
ｃｈ６　チャンネル６受信処理部
ｃｈ７　チャンネル７受信処理部
ｃｈ８　チャンネル８受信処理部
ｃｈ９　チャンネル１受信処理部
ｃｈ１０　チャンネル２受信処理部
ｃｈ１１　チャンネル３受信処理部
ｃｈ１２　チャンネル４受信処理部
ｃｈ１３　チャンネル５受信処理部
ｃｈ１４　チャンネル６受信処理部
ｃｈ１５　チャンネル７受信処理部
ｃｈ１６　チャンネル８受信処理部
１００　複数アンテナ測位装置
１０１－１～１０１－Ｍ　第１アンテナ～第Ｍアンテナ
１０２－１～１０２－Ｍ　第１ＲＦフロントエンド～第ＭＲＦフロントエンド
１０３－１～１０３－Ｎ　デジタル化測位信号
１０４－１～１０４－Ｎ　ｃｈ１からｃｈＮの第１混合器
１０５－１～１０５－Ｎ　ｃｈ１からｃｈＮの相関値演算部
１０６－１～１０６－Ｎ　ｃｈ１からｃｈＮの相関値
１０７－１～１０７－Ｎ　ｃｈ１からｃｈＮの相関値演算部
１０８－１～１０８－Ｎ　ｃｈ１からｃｈＮの相関値演算結果
１０９－１～１０９－Ｎ　ｃｈ１からｃｈＮの搬送波追尾部
１１０－１～１１０－Ｎ　ｃｈ１からｃｈＮのコード追尾部
１１１－１～１１１－Ｎ　ｃｈ１からｃｈＮの搬送波レプリカ信号
１１２－１～１１２－Ｎ　ｃｈ１からｃｈＮのコードレプリカ信号
１１３－１～１１３－Ｎ　ｃｈ１からｃｈＮの第２混合器
１１４－１～１１４－Ｎ　ｃｈ１からｃｈＮのレプリカ信号
１１５－１～１１５－Ｎ　ｃｈ１からｃｈＮのコード遅延時間
１１６－１～１１６－Ｎ　ｃｈ１からｃｈＮの擬似距離演算部
１１７－１～１１７－Ｎ　ｃｈ１からｃｈＮの擬似距離
１１８　擬似距離可否判断部

Claims (11)

  1. 　複数のアンテナと、
    　１または複数の受信回路と、
    　前記アンテナ及び前記受信回路を含む受信経路を介して複数の衛星からの測位信号を受信する受信部と、
    　同一の前記衛星からの測位信号について、少なくとも一部が異なる前記受信経路を介して受信された複数の受信信号に基づいて、測位に用いる前記受信信号を判断する判断部と、
    　前記判断部によって測位に用いると判断された前記受信信号に基づいて、現在位置を算出する位置算出部と、
    　を備えることを特徴とする測位装置。
  2. 　前記判断部は、前記受信信号の受信状態に基づいて、測位に用いる前記受信信号を判断することを特徴とする請求項１に記載の測位装置。
  3. 　前記判断部は、前記複数の受信信号から算出される擬似距離の差に基づいて、測位に用いる前記受信信号を判断することを特徴とする請求項１または２に記載の測位装置。
  4. 　前記判断部は、前記複数の受信信号から算出される擬似距離が予め設定された閾値よりも大きい場合、当該複数の受信信号を測位に用いないように判断することを特徴とする請求項３に記載の測位装置。
  5. 　前記複数のアンテナは、互いに離間して配置され、
    　前記判断部は、前記複数のアンテナのうち、より大きく離間した２つの前記アンテナを介して受信した前記受信信号から算出される前記擬似距離の差に基づいて、測位に用いる前記受信信号を判断することを特徴とする請求項３または４に記載の測位装置。
  6. 　前記判断部は、前記複数のアンテナを介して受信した前記受信信号のうち、受信状態が予め設定された条件に適合する複数の前記受信信号のいずれかを測位に用いるように判断することを特徴とする請求項１から５のいずれか１項に記載の測位装置。
  7. 　前記判断部は、前記受信信号における搬送波位相に基づいて、測位に用いる前記受信信号を判断することを特徴とする請求項１から６のいずれか１項に記載の測位装置。
  8. 　前記判断部は、前記受信信号における信号強度の変化に基づいて、測位に用いる前記受信信号を判断することを特徴とする請求項１から７のいずれか１項に記載の測位装置。
  9. 　前記測位装置の姿勢を検出する姿勢検出部を備え、
    　前記判断部は、前記姿勢検出部によって検出された前記測位装置の姿勢に基づいて、より上側に位置する前記アンテナを含む前記受信経路による前記受信信号を測位に用いると判断することを特徴とする請求項１から８のいずれか１項に記載の測位装置。
  10. 　複数のアンテナと、
    　１または複数の受信回路と、を備える測位装置で実行される測位方法であって、
    　前記アンテナ及び前記受信回路を含む受信経路を介して複数の衛星からの測位信号を受信する受信ステップと、
    　同一の前記衛星からの測位信号について、少なくとも一部が異なる前記受信経路を介して受信された複数の受信信号に基づいて、測位に用いる前記受信信号を判断する判断ステップと、
    　前記判断ステップにおいて測位に用いると判断された前記受信信号に基づいて、現在位置を算出する位置算出ステップと、
    　を含むことを特徴とする測位方法。
  11. 　複数のアンテナと、
    　１または複数の受信回路と、を備える測位装置を制御するコンピュータに、
    　前記アンテナ及び前記受信回路を含む受信経路を介して複数の衛星からの測位信号を受信する受信機能と、
    　同一の前記衛星からの測位信号について、少なくとも一部が異なる前記受信経路を介して受信された複数の受信信号に基づいて、測位に用いる前記受信信号を判断する判断機能と、
    　前記判断機能によって測位に用いると判断された前記受信信号に基づいて、現在位置を算出する位置算出機能と、
    　を実現させることを特徴とするプログラム。
PCT/JP2017/047138 2016-12-27 2017-12-27 測位装置、測位方法及びプログラム WO2018124250A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018559619A JPWO2018124250A1 (ja) 2016-12-27 2017-12-27 測位装置、測位方法及びプログラム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016254158 2016-12-27
JP2016-254158 2016-12-27

Publications (1)

Publication Number Publication Date
WO2018124250A1 true WO2018124250A1 (ja) 2018-07-05

Family

ID=62709540

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/047138 WO2018124250A1 (ja) 2016-12-27 2017-12-27 測位装置、測位方法及びプログラム

Country Status (2)

Country Link
JP (1) JPWO2018124250A1 (ja)
WO (1) WO2018124250A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6445207B1 (ja) * 2018-07-06 2018-12-26 Ｔｅｒｒａ Ｄｒｏｎｅ株式会社 情報収集装置及びこれを搭載した無人航空機
WO2021229886A1 (ja) * 2020-05-11 2021-11-18 ソニーセミコンダクタソリューションズ株式会社 受信装置、および、受信装置の制御方法
US11378697B2 (en) * 2019-10-30 2022-07-05 Qualcomm Incorporated Multipath mitigation for multiband GNSS receiver

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08240654A (ja) * 1995-03-03 1996-09-17 Toyota Motor Corp グローバルポジショニングシステムを用いた位置標定装置
US20060125689A1 (en) * 2004-12-10 2006-06-15 Narayan Anand P Interference cancellation in a receive diversity system
JP2009192900A (ja) * 2008-02-15 2009-08-27 Nikon Corp カメラ
JP2010164339A (ja) * 2009-01-13 2010-07-29 Toyota Motor Corp Ｇｎｓｓ受信装置及び測位方法
JP2011027518A (ja) * 2009-07-24 2011-02-10 Casio Computer Co Ltd 測位装置及びプログラム
JP2011160072A (ja) * 2010-01-29 2011-08-18 Softbank Mobile Corp 信号処理装置及び信号処理プログラム
JP2014153085A (ja) * 2013-02-05 2014-08-25 Railway Technical Research Institute 車両位置計測方法、および車両位置計測システム

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001311767A (ja) * 2000-04-28 2001-11-09 Clarion Co Ltd Ｇｐｓ受信機及びｇｐｓ受信機におけるマルチパス判定方法
JP2005207815A (ja) * 2004-01-21 2005-08-04 Mitsubishi Electric Corp マルチパス信号検出装置
JP2011007683A (ja) * 2009-06-26 2011-01-13 Seiko Epson Corp マルチパス信号信頼度判定方法、コード位相誤差算出方法、コード位相算出方法、位置算出方法及びマルチパス信号信頼度判定装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08240654A (ja) * 1995-03-03 1996-09-17 Toyota Motor Corp グローバルポジショニングシステムを用いた位置標定装置
US20060125689A1 (en) * 2004-12-10 2006-06-15 Narayan Anand P Interference cancellation in a receive diversity system
JP2009192900A (ja) * 2008-02-15 2009-08-27 Nikon Corp カメラ
JP2010164339A (ja) * 2009-01-13 2010-07-29 Toyota Motor Corp Ｇｎｓｓ受信装置及び測位方法
JP2011027518A (ja) * 2009-07-24 2011-02-10 Casio Computer Co Ltd 測位装置及びプログラム
JP2011160072A (ja) * 2010-01-29 2011-08-18 Softbank Mobile Corp 信号処理装置及び信号処理プログラム
JP2014153085A (ja) * 2013-02-05 2014-08-25 Railway Technical Research Institute 車両位置計測方法、および車両位置計測システム

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6445207B1 (ja) * 2018-07-06 2018-12-26 Ｔｅｒｒａ Ｄｒｏｎｅ株式会社 情報収集装置及びこれを搭載した無人航空機
WO2020008706A1 (ja) * 2018-07-06 2020-01-09 Ｔｅｒｒａ Ｄｒｏｎｅ株式会社 情報収集装置及びこれを搭載した無人航空機
JP2020010185A (ja) * 2018-07-06 2020-01-16 Ｔｅｒｒａ Ｄｒｏｎｅ株式会社 情報収集装置及びこれを搭載した無人航空機
US11378697B2 (en) * 2019-10-30 2022-07-05 Qualcomm Incorporated Multipath mitigation for multiband GNSS receiver
US11686853B2 (en) 2019-10-30 2023-06-27 Qualcomm Incorporated Multipath mitigation for multiband GNSS receiver
WO2021229886A1 (ja) * 2020-05-11 2021-11-18 ソニーセミコンダクタソリューションズ株式会社 受信装置、および、受信装置の制御方法

Also Published As

Publication number Publication date
JPWO2018124250A1 (ja) 2019-08-08

Similar Documents

Publication Publication Date Title
US11474258B2 (en) System for determining a physical metric such as position
KR102572948B1 (ko) 로컬 오실레이터를 사용하여 생성된 로컬 신호의 주파수 또는 위상을 정정하는 방법 및 시스템
US10321430B2 (en) Method, apparatus, computer program, chip set, or data structure for correlating a digital signal and a correlation code
US11137500B2 (en) Method, apparatus, computer program, chip set, or data structure for correlating a digital signal and a correlation code
US9612337B2 (en) Antenna beam forming for tracking a transmitter signal
US9581700B2 (en) Method and apparatus tracking global navigation satellite system (GNSS)
WO2018124250A1 (ja) 測位装置、測位方法及びプログラム
CN109313271B (zh) 一种用于使数字信号和相关码相关的方法、装置、计算机程序、芯片组或数据结构
KR20150130989A (ko) 위치 추정에서 멀티패스 효과를 방지하기 위해 송신 다양화를 제공하는 시스템 및 방법
US9720097B2 (en) Multi-antenna radio-navigation signals reception device
JP6165479B2 (ja) ダイバーシチ受信装置、ダイバーシチ受信方法、及びダイバーシチ受信プログラム
KR101271257B1 (ko) 위성항법 신호 추적 장치, 위치 정보 획득 장치, 그리고 위성항법 신호 추적 방법
JP2013246024A (ja) Ｇｐｓモジュール及びｇｐｓによる位置測定方法並びにコンピュータプログラム
JPWO2017022390A1 (ja) 測位装置、測位方法、および、測位プログラム
US20190011568A1 (en) Device for combining multiple satellite navigation signals, and signal processing device comprising same
Lennen et al. Deep Integration of GNSS Receiver Functions into High Volume Platforms

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17886720

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2018559619

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17886720

Country of ref document: EP

Kind code of ref document: A1