測位装置、測位装置の制御方法、測位装置の制御プログラム、測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体
【課題】衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる測位装置等を提供すること。
【解決手段】測位装置20が発生するレプリカ測位基礎符号と測位基礎符号の相関値の最大値に対応する受信周波数であるピーク周波数F0を特定するピーク周波数特定手段と、ピーク周波数F0よりも低い周波数である低周波数F1と、ピーク周波数F0よりも高い周波数である高周波数F2を算出する参照周波数算出手段と、低周波数F1に対応する相関値P1と、高周波数F2に対応する相関値P2を算出する参照相関値算出手段と、ピーク周波数F0に対応する相関値P1及びピーク周波数F0と、低周波数F1に対応する相関値P2及び低周波数F1と、高周波数F2に対応する相関値P3及び高周波数F2とに基づいて、補正後ピーク周波数Frを算出する補正後ピーク周波数算出手段等を有する。
【解決手段】測位装置20が発生するレプリカ測位基礎符号と測位基礎符号の相関値の最大値に対応する受信周波数であるピーク周波数F0を特定するピーク周波数特定手段と、ピーク周波数F0よりも低い周波数である低周波数F1と、ピーク周波数F0よりも高い周波数である高周波数F2を算出する参照周波数算出手段と、低周波数F1に対応する相関値P1と、高周波数F2に対応する相関値P2を算出する参照相関値算出手段と、ピーク周波数F0に対応する相関値P1及びピーク周波数F0と、低周波数F1に対応する相関値P2及び低周波数F1と、高周波数F2に対応する相関値P3及び高周波数F2とに基づいて、補正後ピーク周波数Frを算出する補正後ピーク周波数算出手段等を有する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、測位衛星からの電波を利用する測位装置、測位装置の制御方法、測位装置の制御プログラム、測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体に関するものである。
【背景技術】
【0002】
従来、衛星航法システムである例えば、GPS(Global Positioning System)を利用してGPS受信機の現在位置を測位する測位システムが実用化されている。
このGPS受信機は、GPS衛星の軌道等を示す航法メッセージ(概略衛星軌道情報:アルマナック、精密衛星軌道情報:エフェメリス等を含む)に基づいて、GPS衛星からの電波(以後、衛星電波と呼ぶ)に乗せられている擬似雑音符号(以後、PN(Psuedo random noise code)符号と呼ぶ)の一つであるC/A(Clear and AcquisionまたはCoarse and Access)コードを受信する。C/Aコードは、測位の基礎となる符号である。
GPS受信機は、そのC/AコードがどのGPS衛星から発信されたものであるかを特定したうえで、例えば、そのC/Aコードの発信時刻と受信時刻に基づいて、GPS衛星とGPS受信機の距離(擬似距離)を算出する。そして、GPS受信機は、3個以上のGPS衛星についての擬似距離と、各GPS衛星の衛星軌道上の位置に基づいて、GPS受信機の位置を測位するようになっている(特開平10−339772号公報等参照)。
GPS受信機は、受信したC/AコードとGPS受信機が有しているレプリカのC/Aコードとの間で符号同期を行い、最大の相関値を示す位相(以下、コードフェーズと呼ぶ)を算出する。GPS受信機は、このコードフェーズを利用して、上述の擬似距離を算出することができる。
上述のC/Aコードは衛星電波に乗せられているから、上述の符号同期を正確に行うためには、符号同期と、受信した衛星電波のキャリア周波数(IF(中間)キャリア周波数)とGPS受信機内部の周波数の同期(以下「周波数同期」と呼ぶ)を行う必要がある。
衛星電波の信号強度が強く、例えば、1ミリ秒(ms)という短時間ごとに相関結果(コヒーレント結果)を出力することができる場合には、そのコヒーレント結果に基づいて周波数を修正するPLL(Phase Locked Loop)を構成することによって、周波数同期を行うことができる(例えば、特開2003−98244段落0020参照)。
しかし、衛星電波の強度が弱い場合には、PLLによる周波数同期ができず、いずれ符号同期もできなくなる。
これに対して、本来のIFキャリア周波数を予想して予想IFキャリア周波数を設定し、その予想IFキャリア周波数よりも所定の値だけ高い周波数及び低い周波数における信号レベルの差が低減するようにして、予想IFキャリア周波数を真のIFキャリア周波数に近づける技術が提案されている(例えば、特許文献1)。
【特許文献1】特開2003−255036
【発明の開示】
【発明が解決しようとする課題】
【0003】
上述の技術においては、予想IFキャリア周波数を適切に定める必要がある。ところが、衛星電波の信号強度が極めて弱い場合には、予想IFキャリア周波数を適切に定めることができない場合があるという問題がある。
【0004】
そこで、本発明は、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる測位装置、測位装置の制御方法、測位装置の制御プログラム、測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的とする。
【課題を解決するための手段】
【0005】
前記目的は、第1の発明によれば、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置であって、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号の相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定手段と、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出手段と、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出手段と、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出手段と、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信手段と、を有することを特徴とする測位装置により達成される。
【0006】
第1の発明の構成によれば、前記測位装置は、前記ピーク周波数特定手段を有するから、前記ピーク周波数を特定することができる。
また、前記測位装置は、前記補正後ピーク周波数算出手段を有するから、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数によって規定される第1点と、前記低周波数に対応する前記相関値及び前記低周波数によって規定される第2点と、前記高周波数に対応する前記相関値及び前記高周波数によって規定される第3点の3点に基づいて、補正後ピーク周波数を算出することができる。
前記レプリカ測位基礎符号の位相を固定した場合、前記相関値と受信周波数(IFキャリア周波数)との関係を示すグラフは、理論上、相関値の最大値に対応する点を頂点とする二等辺三角形を描く。そして、上述の第1点はその頂点付近に位置し、第2点及び第3点はそれぞれ異なる斜辺に位置する。そして、第2点及び第3点のいずれか一方は、第1点と同一の斜辺に位置するから、斜辺の傾きを特定することができる。二等辺三角形においては、一方の斜辺の傾きを特定することができれば、他方の斜辺の傾きも特定することができる。そして、二つの斜辺が交わった点が頂点である。そしてこの頂点に対応する周波数が、上述の前記補正後ピーク周波数である。
上述のように、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を適切に定めることができない場合であっても、前記ピーク周波数は必ず1つ存在する。そして、前記ピーク周波数を特定すれば、周波数サーチのサーチ幅による制限を受けることなく、前記補正後ピーク周波数算出手段によって前記補正後ピーク周波数を算出することができる。
そして、前記測位装置は、前記電波受信手段を有するから、前記補正後ピーク周波数を使用して、前記電波を受信することができる。このため、精度よく前記相関値を算出することができ、現在位置を精度良く算出することができる。
これにより、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる。
【0007】
第2の発明は、第1の発明の構成において、前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行う受信周波数制御手段を有することを特徴とする測位装置である。
である。
【0008】
第2の発明の構成によれば、前記測位装置は、前記受信周波数制御手段を有するから、前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行うことができる。
これにより、前記電波の信号強度が所定の強度範囲の場合に、前記受信周波数を継続的に前記電波のIFキャリア周波数に近づけることができる。
【0009】
第3の発明は、第2の発明の構成において、前記補正後ピーク周波数算出手段と前記受信周波数制御手段とが平行して作動することを特徴とする請求項2に記載の測位装置である。
【0010】
第3の発明の構成によれば、前記電波の信号強度が所定の強度よりも大きい場合には、前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行うことができる。また、前記電波の信号強度が所定の強度よりも小さい場合に、前記補正後ピーク周波数を使用して、前記電波を受信することができる。
このため、前記受信強度が所定の強度よりも大きい状態からから小さい状態に移行した場合に、連続的に精度良く測位を行うことができる。
【0011】
第4の発明は、第1の発明乃至第3の発明のいずれかの構成において、前記発信源は、測位衛星であることを特徴とする測位装置である。
【0012】
前記目的は、第5の発明によれば、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、を有することを特徴とする測位装置の制御方法である。
【0013】
第5の発明の構成によれば、第1の発明の構成と同様に、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる。
【0014】
前記目的は、第6の発明によれば、コンピュータに、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、を実行させることを特徴とする測位装置の制御プログラムによって達成される。
【0015】
前記目的は、第7の発明によれば、コンピュータに、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、を実行させることを特徴とする測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体によって達成される。
【発明を実施するための最良の形態】
【0016】
以下、この発明の好適な実施の形態を添付図面等を参照しながら、詳細に説明する。
尚、以下に述べる実施の形態は、本発明の好適な具体例であるから、技術的に好ましい種々の限定が付されているが、本発明の範囲は、以下の説明において特に本発明を限定する旨の記載がない限り、これらの態様に限られるものではない。
【0017】
図1は、本発明の実施の形態の端末20等を示す概略図である。
図1に示すように、端末20は、複数の測位衛星である例えば、GPS衛星12a,12b,12c及び12dから、電波S1,S2,S3及びS4を受信することができる。電波S1等は電波の一例である。GPS衛星12a等は、発信源の一例でもある。
電波S1等には各種のコード(符号)が乗せられている。そのうちの一つがC/Aコードである。このC/Aコードは、1,023チップ(chip)から構成されている。そして、このC/Aコードは、1.023Mbpsのビット率、1,023bit(=1msec)のビット長の信号である。このC/Aコードは、測位基礎符号の一例である。そして、端末20は、現在位置を測位する測位装置の一例である。
端末20は、自動車15に搭載されており、自動車15の移動に連れて移動しつつ、現在位置を測位するようになっている。
【0018】
端末20は、例えば、3個以上の異なるGPS衛星12a等からのC/Aコードを受信して、現在位置を測位することができるようになっている。
端末20は、まず、受信したC/AコードがどのGPS衛星に対応するものかを特定する。次に、端末20は、相関処理によって、受信したC/Aコードの位相(以下「コードフェーズ」と呼ぶ)を算出する。続いて、端末20は、そのコードフェーズを使用して、各GPS衛星12a等と端末20との距離(以後、擬似距離と呼ぶ)を算出する。続いて、現在時刻における各GPS衛星12a等の衛星軌道上の位置と、上述の擬似距離に基づいて、現在位置の測位演算を行うことができるように構成されている。
C/Aコードは、電波S1等に乗せられているから、端末20が、電波S1等を受信するときの受信周波数が不正確であると、相関処理によって算出するコードフェーズの精度も劣化する。GPS衛星12a等はその衛星軌道上を移動しているから、この受信周波数は変化し続けているが、電波S1等の信号強度が強い場合には、電波S1等を利用したPLLによって周波数同期を確保し続けることができる。
ところが、電波S1等の信号強度が極めて弱い場合には、PLLは有効に機能しない。また、電波S1等の信号強度が極めて弱い場合には、電波S1等のIFキャリア周波数を的確に予想することも困難である。
この点、端末20は、以下に説明するように、電波S1等の信号強度が極めて弱い場合
に、IFキャリア周波数を予想することなく、精度よく現在位置を測位することができるようになっている。
【0019】
端末20は例えば、携帯電話機、PHS(Personal Handy−phone System)、PDA(Personal Digital Assistance)等であるが、これらに限らない。
また、GPS衛星12a等は、4個に限らず、3個でもよいし、5個以上でもよい。
【0020】
(端末20の主なハードウエア構成について)
図2は、端末20の主なハードウエア構成を示す概略図である。
図2に示すように、端末20は、コンピュータを有し、コンピュータは、バス22を有する。バス22には、CPU(Central Processing Unit)24、記憶装置26等が接続されている。記憶装置26は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)等である。
また、バス22には、外部記憶装置28が接続されている。外部記憶装置28は、例えば、HDD(Hard Disc Drive)である。
また、バス22には、電源装置30、入力装置32、GPS装置34、表示装置47、時計48が接続されている。
【0021】
(GPS装置34の構成について)
図3は、GPS装置34の構成を示す概略図である。
図3に示すように、GPS装置34は、RF部35とベースバンド部36で構成される。
RF部35は、アンテナ35aで電波S1等を受信する。そして、増幅器であるLNA35bが、電波S1に乗せられているC/Aコード等の信号を増幅する。そして、ミキサ35cが、信号の周波数を中間(IF)キャリア周波数にダウンコンバートする。そして、直交(IQ)検波器35dが信号をIQ分離する。続いて、ADコンバータ35e1及び35e2が、IQ分離された信号をそれぞれデジタル信号に変換するように構成されている。
ベースバンド部36は、RF部35からデジタル信号に変換されたIFキャリア周波数の信号を受信する。
ベースバンド部36の相関部37は、入力したデジタル信号を例えば、10ミリ秒(ms)の時間において同期積算し、その積算結果とレプリカC/Aコードとの相関をとる処理であるコヒーレントを行う。相関部37は、NCO38、コードジェネレータ39及び相関器40から構成される。コードジェネレータ39は、NCO38が発生するクロックのタイミングでレプリカC/Aコードを発生する。相関器40はC/AコードとレプリカC/Aコードとの相関をとり、位相の特定及び相関値の算出を行う。相関部37には、周波数及びレプリカC/Aコードの位相を設定することができる。
信号積算器41は、相関部37から出力された相関値を積算する処理であるインコヒーレントを行う。
コードフェーズ検出器42は、相関部37の出力値と信号積算器41からの出力値から、コードフェーズを検出する。
【0022】
上述のように、相関処理は、コヒーレントと、インコヒーレントから構成される。
コヒーレントは、相関部37が、受信したC/AコードとレプリカC/Aコードの相関をとる処理である。
例えば、コヒーレント時間が20msecであれば、20msecの時間において同期積算したC/AコードとレプリカC/Aコードとの相関値等を算出する。コヒーレント処理の結果、相関をとった位相と、相関値が出力される。
インコヒーレントは、コヒーレント結果の相関値を積算することによって、インコヒーレント値を算出する処理である。
相関処理の結果、コヒーレント処理で出力された位相と、インコヒーレント値が出力される。相関値Pはインコヒーレント値である。
電波S1等の信号強度が十分に強い場合には、位相検出器43が相関器40から位相情報を取得して、NCO38に供給し、PLLを構成することができる。この結果、IFキャリア周波数と同期した周波数において、レプリカC/Aコードを発生させることができる。具体的には、相関値Pが最大になるように、受信周波数の制御が行われる。
位相検出器43、相関器40及びNCO38は、受信周波数制御手段の一例である。
【0023】
(端末20の主なソフトウエア構成について)
図4は、端末20の主なソフトウエア構成を示す概略図である。
図4に示すように、端末20は、各部を制御する制御部100、図2のGPS装置34に対応するGPS部102、時計48に対応する計時部104、各種プログラムを格納する第1記憶部110、各種情報を格納する第2記憶部150を有する。
【0024】
図4に示すように、端末20は、第2記憶部150に、衛星軌道情報152を格納している。衛星軌道情報152は、アルマナック152a及びエフェメリス152bを含む。アルマナック152aは、すべてのGPS衛星12a等の概略の軌道を示す情報である。エフェメリス152bは、各GPS衛星12a等の精密な軌道を示す情報である。
端末20は、アルマナック152a及びエフェメリス152bを、測位のために使用する。
【0025】
図4に示すように、端末20は、第2記憶部150に、初期位置情報154を格納している。初期位置情報154は、端末20の現在の初期位置P0を示す情報である。初期位置Q0は、例えば、前回測位時の測位位置である。
【0026】
図4に示すように、端末20は、第1記憶部110に、観測可能衛星算出プログラム112を格納している。観測可能衛星算出プログラム112は、制御部100が、アルマナック152aを参照して、計時部104によって計測した現在時刻において初期位置Q0から観測可能なGPS衛星12a等を示す観測可能衛星情報156を生成するためのプログラムである。
制御部100は、生成した観測可能衛星情報156を第2記憶部150に格納する。
【0027】
図4に示すように、端末20は、第1記憶部110に、第1推定周波数算出プログラム114を格納している。第1推定周波数算出プログラム114は、制御部100が、電波S1等ごとのIFキャリア周波数の予測値である第1推定周波数Aを算出するためのプログラムである。第1推定周波数Aは、端末20が、例えば、現在時刻においてGPS衛星12aからの電波S1を受信するときの電波S1のIFキャリア周波数の予測値である。
【0028】
図5は、第1推定周波数算出プログラム114の説明図である。
図5に示すように、第1推定周波数Aは、発信周波数H1からドップラー偏移H2を加えた周波数である。発信周波数H1は、電波S1等がGPS衛星12a等から発信されるときの周波数である例えば、1.5GHzと、ミキサ35cによるダウンコバート率によって決定される既知の値である。ドップラー偏移H2は、GPS衛星12a等と端末20の相対移動によって生じる周波数偏移であり、常に変動している。ドップラー偏移H2は、端末20の初期位置P0とエフェメリス152bによって、算出することができる。
制御部100は、第1推定周波数Aを示す第1推定周波数情報158を第2記憶部150に格納する。
【0029】
ただし、端末20の位置が正確な現在位置ではなくて初期位置Q0であること、GPS衛星12a等と端末20は常に相対移動を行っていることによって、算出したドップラー偏移H1は真のドップラー偏移とは乖離する可能性がある。
このため、第1推定周波数Aは、真のIFキャリア周波数とは乖離するのが通常である。
【0030】
図4に示すように、端末20は、第1記憶部110に、第1相関プログラム116を格納している。第1相関プログラム116は、制御部100が、GPS衛星12a等から受信したC/AコードとレプリカC/Aコードとの相関値を算出し、さらに、C/Aコードの位相(コードフェーズ)である第1位相CP1を算出するためのプログラムである。
なお、第1位相CP1は、C/Aコードの位相であり、レプリカC/Aコードの位相でもある。
【0031】
図6は、第1相関プログラム116の説明図である。
図6(a)に示すように、制御部100はベースバンド部36によって、C/Aコードの1チップを例えば、等間隔で分割して、相関処理を行う。C/Aコードの1チップは、例えば、32等分される。すなわち、32分の1チップの位相幅(第1位相幅W1)間隔で相関処理を行う。制御部100が相関処理を行うときの第1位相幅W1間隔の位相を第1サンプリング位相SC1と呼ぶ。
第1位相幅W1は、アンテナ35aに入力する信号の信号強度が−155dBm以上である場合に、相関最大値Pmaxを検出することができる位相幅として規定されている。32分の1チップの位相幅であれば、信号強度が−155dBm以上であれば弱電界であっても、相関最大値Pmaxを検出することができることがシミュレーションによって明らかになっている。
【0032】
図6(b)に示すように、制御部100は、推定周波数A±100kHzの範囲を、100Hz単位で変動させながら相関処理を行う。各周波数ごとに、コードフェーズCPを第1位相幅W1で変動させて、相関最大値Pmaxを算出することができる周波数及びコードフェーズを特定する。
なお、測位開始時においては、レプリカC/Aコードを0から1023チップまで変動させる。
そして、一度相関最大値Pmaxに対応するコードフェーズ及び周波数を特定すると、その後は、相関最大値Pmaxに対応するコードフェーズ及び周波数を中心にして、測位開始時よりも狭い範囲において、信号S1等のサーチを行う。例えば、制御部100は、既に算出している第1測位位相CP1を中心に、±256チップの位相範囲をサーチする。また、周波数については、相関最大値Pmaxに対応する周波数を中心に、±1.0kHzの範囲を100Hz単位でサーチする。この条件を、第1トラッキング条件と呼ぶ。
【0033】
図6(c)に示すように、ベースバンド部36からは、2チップ分の位相C1乃至C64に対応する相関値Pが出力される。各位相C1乃至C64が、第1サンプリング位相SC1である。
Pnoiseに対するPmaxの比をSNRと呼ぶ。Pnoiseは環境雑音の信号レベルである。PmaxはGPS衛星12a等からの信号レベルである。
信号S1等の信号強度が弱い状態においては、図6(c)SNR1は、比較的小さい。
【0034】
制御部100は、相関値Pmaxに対応する第1位相CP1を特定する。
制御部100は、第1位相CP1を示す第1位相情報160を第2記憶部位150に格納する。
SNR1が小さいほど、第1位相CP1の精度は低い。
なお、上述の第1相関プログラム116に基づく端末20の動作を、第1相関処理と呼ぶ。
【0035】
図4に示すように、端末20は、第1記憶部110に、第1測位プログラム118を格納している。第1測位プログラム118は、制御部100が、3個以上のGPS衛星12a等に対応する第1位相CP1に基づいて、現在位置を測位して測位位置Q1を算出するためのプログラムである。
【0036】
図7は、測位方法を示す概念図である。
図7に示すように、例えば、GPS衛星12aと端末20との間には、複数のC/Aコードが連続的に並んでいると観念することができる。そして、GPS衛星12aと端末20との間の距離は、C/Aコードの長さの整数倍とは限らないから、コード端数C/Aaが存在する。つまり、GPS衛星12aと端末20との間には、C/Aコードの整数倍の部分(C/Aコードがn個(nは整数)並んでいる部分)と、端数部分(コード端数C/Aa)が存在する。C/Aコードの整数倍の部分とコード端数C/Aaの合計の長さが擬似距離である。端末20は、この擬似距離を使用して測位を行う。
GPS衛星12aの軌道上の位置はエフェメリス152bを使用して算出可能である。そして、GPS衛星12aの軌道上の位置と初期位置Q0との距離を算出すれば、C/Aコードの整数倍の部分を特定することができる。
そして、図7に示すように、レプリカC/Aコードの位相を例えば、矢印X1方向に移動させながら、相関処理を行を行う。
相関値が最大になった位相がコード端数C/Aaである。そして、このコード端数C/Aaが、第1位相CP1である。
【0037】
制御部100は、3個以上のGPS衛星12a等に対応する第1位相CP1に基づいて、各GPS衛星12a等と端末20との擬似距離を算出する。そして、各GPS衛星12a等の軌道上の位置はエフェリス152bによって算出する。そして、3個以上のGPS衛星12a等の軌道上の位置と、擬似距離に基づいて、現在位置を測位して測位位置Q1を算出する。
制御部100は、測位位置Q1を示す第1測位位置情報162を第2記憶部150に格納する。
【0038】
図4に示すように、端末20は、第1記憶部110に、測位位置出力プログラム120を格納している。測位位置出力プログラム120は、制御部100が、測位位置Q1又は後述の測位位置Q2を表示装置47に表示するためのプログラムである。
【0039】
図4に示すように、端末20は、第1記憶部110に、第2相関プログラム122を格納している。第2相関プログラム122は、制御部100が相関処理を行い、相関値Pとコードフェ−ズCPを算出するためのプログラムである。
制御部100は、相関値Pとコードフェ−ズCP示す第2相関情報164を第2記憶部150に格納する。
【0040】
図4に示すように、端末20は、第1記憶部110に、ピーク周波数特定プログラム124を格納している。ピーク周波数特定プログラム124と制御部100は、ピーク周波数特定手段の一例である。
【0041】
図8は、ピーク周波数特定プログラム124の説明図である。
図8に示すように、制御部100は、相関最大値Pmaxに対応する周波数をピーク周波数F0として特定する。ピーク周波数F0はピーク周波数の一例である。
このピーク周波数F0は、端末20が例えば、100Hzの幅でサーチを行った結果であるから、受信した電波S1等の真のIFキャリア周波数と、最大約50Hzの乖離を生じている。
制御部100は、ピーク周波数F0を示すピーク周波数情報166を第2記憶部150に格納する。
【0042】
図4に示すように、端末20は、第1記憶部110に、参照周波数算出プログラム126を格納している。参照周波数算出プログラム126と制御部100は、参照周波数算出手段の一例である。
制御部100は参照周波数算出プログラム126に基づいて、ピーク周波数F0よりも100Hz低い周波数F1と、ピーク周波数F0よりも100Hz高い周波数F2を算出する。制御部100は、周波数F1及び周波数F2を示す参照周波数情報168を第2記憶部150に格納する。周波数F1は低周波数の一例である。周波数F2は高周波数の一例である。
ピーク周波数F0と周波数F1との周波数差分と、ピーク周波数F0と周波数F2との周波数差分が等しくなるように、周波数F1と周波数F1が規定される。本実施の形態では、周波数差分は100Hzに設定されている。
なお、本実施の形態とは異なり、周波数差分は100Hzに限らない。
【0043】
図4に示すように、端末20は、第1記憶部110に、参照相関値算出プログラム128を格納している。参照相関値算出プログラム128と制御部100は、参照相関値算出手段の一例である。
制御部100は参照相関値算出プログラム128に基づいて、周波数F1に対応する相関値P1と、周波数F2に対応する相関値P2を算出する。具体的には、制御部100は、第2相関情報164を参照して、相関値P1及び相関値P2を算出する。
制御部100は、相関値P1及び相関値P2を示す参照相関値情報170を第2記憶部150に格納する。
【0044】
図4に示すように、端末20は、第1記憶部110に、第2推定周波数算出プログラム130を格納している。第2推定周波数算出プログラム130は、制御部100が、ピーク周波数F0と相関ピーク値Pmax(P0)、周波数F1と相関値P1、及び、周波数F2と相関値P2とに基づいて、第2推定周波数Frを算出するためのプログラムである。第2推定周波数Frは、補正後ピーク周波数の一例である。第2推定周波数算出プログラム130と制御部100は、補正後ピーク周波数算出手段の一例である。
【0045】
図9及び図10は、第2推定周波数算出プログラム130の説明図である。
相関値Pと周波数Fを示すグラフは、図9及び図10に示すように、二等辺三角形を描く。
図9(a)及び図10(a)に示すように、ピーク周波数F0と相関ピーク値P0によって点G0が規定される。周波数F1と相関値P1によって点G1が規定される。そして、周波数F2と相関値P2とによって点G2が規定される。
【0046】
図9(a)及び図9(b)に示すように、相関値P1が相関値P2よりも小さい場合には、点G0と点G1は傾きa(aは正の数)の同一直線上にある。点G0と点G1を結んだ直線が直線L1である。
そして、点G2は、傾きが−aの直線上にある。傾きが−aで、かつ、点G2を通る直線が直線L2である。
そして、直線L1と直線L2の交点が二等辺三角形の頂点Hである。頂点Hに対応する周波数が、第2推定周波数Frである。図9(b)の連立方程式1を解くと、未知数Fr、Pr及び傾きaを算出することができる。
【0047】
図10(a)及び図10(b)に示すように、相関値P1が相関値P2よりも大きい場合には、点G0と点G2は傾き−a(aは正の数)の同一直線上にある。点G0と点G2を結んだ直線が直線L2である。
そして、点G1は、傾きがaの直線上にある。傾きがaで、かつ、点G1を通る直線が直線L1である。
そして、直線L1と直線L2の交点が二等辺三角形の頂点Hである。頂点Hに対応する周波数が、第2推定周波数Frである。図10(b)の連立方程式2を解くと、未知数Fr、Pr及び傾きaを算出することができる。
【0048】
なお、相関値P1が相関値P2と等しい場合には、ピーク周波数F0が第2推定周波数Frである。
制御部100は、第2推定周波数Frを示す第2推定周波数情報172を第2記憶部150に格納する。
この第2推定周波数Frは、周波数Fのサーチステップである100Hzの制限を受けないから、精度の高い情報である。すなわち、ピーク周波数F0よりも、真のIFキャリア周波数に近い。
【0049】
図4に示すように、端末20は、第1記憶部110に、第2位相決定プログラム132を格納している。第2位相決定プログラム132は、制御部100が、第2推定周波数Frを使用して、電波S1等を受信して、相関処理を行い、測位を行うための第2位相CP2を算出するためのプログラムである。第2位相決定プログラム132と制御部100は、電波受信手段の一例である。
【0050】
図11は、第2位相決定プログラム132の説明図である。
図11の相関グラフにおけるSNR2は、図6(c)のグラフにおけるSNR1よりも大きい。これは、第2推定周波数Frが、真のIFキャリア周波数に非常に近いからである。
このため、相関最大値Pmaxに対応する位相である第2位相CP2は精度の高い位相情報である。
制御部100は、第2位相CP2を示す第2位相情報174を第2記憶部150に格納する。
上述の第2相関プログラム122、ピーク周波数特定プログラム124、参照周波数算出プログラム126、参照相関値算出プログラム128及び第2推定周波数算出プログラム130、第2位相決定プログラム132に基づく端末20の動作を、第2相関処理と呼ぶ。
【0051】
図4に示すように、端末20は、第1記憶部110に、第2測位プログラム134を格納している。第2測位プログラム134は、制御部100が、3個以上のGPS衛星12aについての第2位相CP2を使用して、測位を行い、測位位置Q2を算出するためのプログラムである。
制御部100は、測位位置Q2を示す第2測位位置情報176を第2記憶部150に格納する。
【0052】
図4に示すように、端末20は、第1記憶部110に、信号強度評価プログラム136を格納している。
信号強度評価プログラム136は、アンテナ35aに入力する信号の信号強度SPを評価するためのプログラムである。アンテナ35aに入力する信号の信号強度SPは、相関値から推定することができる。
制御部100は、例えば、信号強度SPが−138dBm以上である場合には、第1相関処理を行って、測位位置Q1を算出する。
そして、制御部100は、信号強度SPが−142dBm以下である場合には、第2相関処理を行って、測位位置Q2を算出する。
そして、制御部100は、信号強度SPが−142dBmより大きく−138dBm未満である場合には、第1相関処理と第2相関処理を平行して実施する。そして、制御部100は、第1位相CP1を使用して、測位位置Q1を算出する。
【0053】
端末20は、以上のように構成されている。
上述のように、端末20は、ピーク周波数F0(図4参照)を特定することができる。
また、端末20は、第2推定周波数Fr(図4参照)を算出することができる。
レプリカC/Aコードの位相を固定した場合、相関値と受信周波数(IFキャリア周波数)との関係を示すグラフは、図9に示すように、相関値の最大値に対応する点を頂点とする二等辺三角形を描く。そして、ピーク周波数F0に対応する点G0はその頂点H付近に位置し、ピーク周波数F0の前後の周波数F1及びF2にそれぞれ対応する点G1及びG2はそれぞれ異なる斜辺に位置する。そして、点G1及び点G2のいずれか一方は、点G1と同一の斜辺に位置するから、斜辺の傾きaを特定することができる。二等辺三角形においては、一方の斜辺の傾きを特定することができれば、他方の斜辺の傾きも特定することができる。そして、二つの斜辺が交わった点が頂点Hである。そしてこの頂点Hに対応する周波数が、上述の第2推定周波数Frである。
上述のように、電波S1等の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることができない場合であっても、ピーク周波数F0は必ず1つ存在する。そして、ピーク周波数F0を特定すれば、第2推定周波数Frを算出することができる。
そして、端末20は、第2推定周波数Frを使用して、電波S1等を受信することができる。このため、精度よく相関値Pを算出することができ、現在位置を精度よく算出することができる。
これにより、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる。
【0054】
また、端末20は、PLLによって、レプリカC/Aコードと受信したC/Aコードとのコヒーレント値が最大になるように受信周波数の制御を行うことができる。
これにより、電波S1等の信号強度が所定の強度範囲の場合に、PLLを有効に機能させ、受信周波数を継続的に電波S1等のIFキャリア周波数に近づけることができる。
【0055】
また、端末20は、電波S1等の信号強度が所定範囲の場合には、上述の第1相関処理と第2相関処理を平行して行うことができる。このため、信号強度SPが所定の強度よりも大きい状態からから小さい状態に移行した場合に、連続的に精度良く測位を行うことができる。
【0056】
以上が本実施の形態に係る端末20の構成であるが、以下、その動作例を主に図12及び図13を使用して説明する。
図12及び図13は端末20の動作例を示す概略フローチャートである。
【0057】
まず、端末20は、各GPS衛星12a等について、エフェメリス152bと初期位置Q0から、推定周波数Aを算出する(図12のステップST1)。
続いて、端末20は、第1相関処理を行う(ステップST2)。
【0058】
続いて、端末20は、信号強度SPを判断する(ステップST3)。
端末20は、ステップST3において、信号強度SPが−138dBm以上であると判断すると、第1相関処理を継続し(ステップST4A)、第1位相CP1を使用して現在位置を測位して、測位位置Q1を算出する(ステップST5A)。
続いて、端末20は、測位位置Q1を出力する(ステップST6A)。
続いて、端末20は、測位が測位規定回数である例えば、10回に達したか否かを判断する(ステップST7)。
端末20は、測位が測位規定回数に達したと判断した場合には測位を終了する。
端末20は、測位が測位規定回数に達していない判断した場合にはステップST3以降を実施する。
【0059】
端末20は、ステップST3において、信号強度SPが−142dBm以下であると判断すると、第1相関処理を停止し、第2相関処理を行う(ステップST4B)。
第2相関処理においては、端末20は、まず、ピーク周波数F0(図4参照)を特定する(図13のステップST101)。このステップST101は、ピーク周波数特定ステップの一例である。
続いて、端末20は、周波数F1及びF2(図4参照)を算出する(ステップST102)。このステップST102は、参照周波数算出ステップの一例である。
【0060】
続いて、端末20は、相関値P2及びP3(図4参照)を算出する(ステップST103)。このステップST103は、参照相関値算出ステップの一例である。
続いて、端末20は、第2推定周波数Fr(図4参照)を算出する(ステップST104)。このステップST104は、補正後ピーク周波数算出ステップの一例である。
【0061】
続いて、端末20は、第2位相CP2を算出し、その第2位相CP2を使用して現在位置を測位して、測位位置Q2を算出する(図12のステップST5B)。
続いて、端末20は、測位位置Q2を出力し(ステップST6B)、ステップST7を実施する。
【0062】
端末20は、ステップST3において、信号強度SPが−142dBmより大きく−138dBm未満であると判断すると、第1相関処理と第2相関処理を平行実施する(図12のステップST4C)。
続いて、端末20は、第1位相CP1を使用して現在位置を測位して、測位位置Q1を算出する(ステップST5C)。
続いて、端末20は、測位位置Q1を出力し(ステップST6C)、ステップST7を実施する。
端末20は、測位が測位規定回数に達していない判断した場合にはステップST3以降を実施する。再度のステップST3において、端末20は、信号強度SPが−138dBm以下であると判断すると、ステップST4Bに進む。ここで、第1相関処理と第2相関処理が平行して継続しているから、第1相関処理を停止し、直ちに第2相関処理を実施することができる。これは、第1相関処理においてPLLが機能しなくなった後に、第2相関処理を開始するのではなくて、信号強度SPが−142dBm以下に低下するおそれがある中程度の状態(信号強度SPが−142dBmより大きく−138dBm未満の状態)で、第2相関処理を継続しておくことを意味する。このため、第2相関処理において、新たに広範囲の周波数及び位相をサーチする必要がないから、迅速にステップST5B以下を実施することができる。
また、中程度の状態(信号強度SPが−142dBmより大きく−138dBm未満の状態)は、信号強度SPが−138dBm以上になる可能性がある状態でもある。予め第1相関処理を継続しておくことによって、信号強度が−138dBm以上になった場合に、直ちに第1相関処理だけを実施する状態に移行することができる。
【0063】
本発明は、上述の各実施の形態に限定されない。さらに、上述の各実施の形態は、相互に組み合わせて構成するようにしてもよい。
【図面の簡単な説明】
【0064】
【図1】本発明の実施の形態の端末等を示す概略図である。
【図2】端末の主なハードウエア構成を示す概略図である。
【図3】GPS装置の構成を示す概略図である。
【図4】端末の主なソフトウエア構成を示す概略図である。
【図5】第1推定周波数算出プログラムの説明図である。
【図6】第1相関プログムの説明図である。
【図7】測位方法を示す概念図である。
【図8】ピーク周波数特定プログラムの説明図である。
【図9】第2推定周波数算出プログラムの説明図である。
【図10】第2推定周波数算出プログラムの説明図である。
【図11】第2位相決定プログラムの説明図である。
【図12】端末の動作例を示す概略フローチャートである。
【図13】端末の動作例を示す概略フローチャートである。
【符号の説明】
【0065】
12a,12b,12c,12d・・・GPS衛星、20・・・端末、34・・・GPS装置、112・・・観測可能衛星算出プログラム、114・・・第1推定周波数算出プログラム、116・・・第1相関プログラム、118・・・第1測位プログラム、120・・・測位位置出力プログラム、122・・・第2相関プログラム、124・・・ピーク周波数特定プログラム、126・・・参照周波数算出プログラム、128・・・参照相関値算出プログラム、130・・・第2推定周波数算出プログラム、132・・・第2位相特定プログラム、134・・・第2測位プログラム、136・・・信号強度評価プログラム
【技術分野】
【0001】
本発明は、測位衛星からの電波を利用する測位装置、測位装置の制御方法、測位装置の制御プログラム、測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体に関するものである。
【背景技術】
【0002】
従来、衛星航法システムである例えば、GPS(Global Positioning System)を利用してGPS受信機の現在位置を測位する測位システムが実用化されている。
このGPS受信機は、GPS衛星の軌道等を示す航法メッセージ(概略衛星軌道情報:アルマナック、精密衛星軌道情報:エフェメリス等を含む)に基づいて、GPS衛星からの電波(以後、衛星電波と呼ぶ)に乗せられている擬似雑音符号(以後、PN(Psuedo random noise code)符号と呼ぶ)の一つであるC/A(Clear and AcquisionまたはCoarse and Access)コードを受信する。C/Aコードは、測位の基礎となる符号である。
GPS受信機は、そのC/AコードがどのGPS衛星から発信されたものであるかを特定したうえで、例えば、そのC/Aコードの発信時刻と受信時刻に基づいて、GPS衛星とGPS受信機の距離(擬似距離)を算出する。そして、GPS受信機は、3個以上のGPS衛星についての擬似距離と、各GPS衛星の衛星軌道上の位置に基づいて、GPS受信機の位置を測位するようになっている(特開平10−339772号公報等参照)。
GPS受信機は、受信したC/AコードとGPS受信機が有しているレプリカのC/Aコードとの間で符号同期を行い、最大の相関値を示す位相(以下、コードフェーズと呼ぶ)を算出する。GPS受信機は、このコードフェーズを利用して、上述の擬似距離を算出することができる。
上述のC/Aコードは衛星電波に乗せられているから、上述の符号同期を正確に行うためには、符号同期と、受信した衛星電波のキャリア周波数(IF(中間)キャリア周波数)とGPS受信機内部の周波数の同期(以下「周波数同期」と呼ぶ)を行う必要がある。
衛星電波の信号強度が強く、例えば、1ミリ秒(ms)という短時間ごとに相関結果(コヒーレント結果)を出力することができる場合には、そのコヒーレント結果に基づいて周波数を修正するPLL(Phase Locked Loop)を構成することによって、周波数同期を行うことができる(例えば、特開2003−98244段落0020参照)。
しかし、衛星電波の強度が弱い場合には、PLLによる周波数同期ができず、いずれ符号同期もできなくなる。
これに対して、本来のIFキャリア周波数を予想して予想IFキャリア周波数を設定し、その予想IFキャリア周波数よりも所定の値だけ高い周波数及び低い周波数における信号レベルの差が低減するようにして、予想IFキャリア周波数を真のIFキャリア周波数に近づける技術が提案されている(例えば、特許文献1)。
【特許文献1】特開2003−255036
【発明の開示】
【発明が解決しようとする課題】
【0003】
上述の技術においては、予想IFキャリア周波数を適切に定める必要がある。ところが、衛星電波の信号強度が極めて弱い場合には、予想IFキャリア周波数を適切に定めることができない場合があるという問題がある。
【0004】
そこで、本発明は、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる測位装置、測位装置の制御方法、測位装置の制御プログラム、測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的とする。
【課題を解決するための手段】
【0005】
前記目的は、第1の発明によれば、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置であって、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号の相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定手段と、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出手段と、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出手段と、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出手段と、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信手段と、を有することを特徴とする測位装置により達成される。
【0006】
第1の発明の構成によれば、前記測位装置は、前記ピーク周波数特定手段を有するから、前記ピーク周波数を特定することができる。
また、前記測位装置は、前記補正後ピーク周波数算出手段を有するから、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数によって規定される第1点と、前記低周波数に対応する前記相関値及び前記低周波数によって規定される第2点と、前記高周波数に対応する前記相関値及び前記高周波数によって規定される第3点の3点に基づいて、補正後ピーク周波数を算出することができる。
前記レプリカ測位基礎符号の位相を固定した場合、前記相関値と受信周波数(IFキャリア周波数)との関係を示すグラフは、理論上、相関値の最大値に対応する点を頂点とする二等辺三角形を描く。そして、上述の第1点はその頂点付近に位置し、第2点及び第3点はそれぞれ異なる斜辺に位置する。そして、第2点及び第3点のいずれか一方は、第1点と同一の斜辺に位置するから、斜辺の傾きを特定することができる。二等辺三角形においては、一方の斜辺の傾きを特定することができれば、他方の斜辺の傾きも特定することができる。そして、二つの斜辺が交わった点が頂点である。そしてこの頂点に対応する周波数が、上述の前記補正後ピーク周波数である。
上述のように、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を適切に定めることができない場合であっても、前記ピーク周波数は必ず1つ存在する。そして、前記ピーク周波数を特定すれば、周波数サーチのサーチ幅による制限を受けることなく、前記補正後ピーク周波数算出手段によって前記補正後ピーク周波数を算出することができる。
そして、前記測位装置は、前記電波受信手段を有するから、前記補正後ピーク周波数を使用して、前記電波を受信することができる。このため、精度よく前記相関値を算出することができ、現在位置を精度良く算出することができる。
これにより、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる。
【0007】
第2の発明は、第1の発明の構成において、前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行う受信周波数制御手段を有することを特徴とする測位装置である。
である。
【0008】
第2の発明の構成によれば、前記測位装置は、前記受信周波数制御手段を有するから、前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行うことができる。
これにより、前記電波の信号強度が所定の強度範囲の場合に、前記受信周波数を継続的に前記電波のIFキャリア周波数に近づけることができる。
【0009】
第3の発明は、第2の発明の構成において、前記補正後ピーク周波数算出手段と前記受信周波数制御手段とが平行して作動することを特徴とする請求項2に記載の測位装置である。
【0010】
第3の発明の構成によれば、前記電波の信号強度が所定の強度よりも大きい場合には、前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行うことができる。また、前記電波の信号強度が所定の強度よりも小さい場合に、前記補正後ピーク周波数を使用して、前記電波を受信することができる。
このため、前記受信強度が所定の強度よりも大きい状態からから小さい状態に移行した場合に、連続的に精度良く測位を行うことができる。
【0011】
第4の発明は、第1の発明乃至第3の発明のいずれかの構成において、前記発信源は、測位衛星であることを特徴とする測位装置である。
【0012】
前記目的は、第5の発明によれば、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、を有することを特徴とする測位装置の制御方法である。
【0013】
第5の発明の構成によれば、第1の発明の構成と同様に、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる。
【0014】
前記目的は、第6の発明によれば、コンピュータに、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、を実行させることを特徴とする測位装置の制御プログラムによって達成される。
【0015】
前記目的は、第7の発明によれば、コンピュータに、複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、を実行させることを特徴とする測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体によって達成される。
【発明を実施するための最良の形態】
【0016】
以下、この発明の好適な実施の形態を添付図面等を参照しながら、詳細に説明する。
尚、以下に述べる実施の形態は、本発明の好適な具体例であるから、技術的に好ましい種々の限定が付されているが、本発明の範囲は、以下の説明において特に本発明を限定する旨の記載がない限り、これらの態様に限られるものではない。
【0017】
図1は、本発明の実施の形態の端末20等を示す概略図である。
図1に示すように、端末20は、複数の測位衛星である例えば、GPS衛星12a,12b,12c及び12dから、電波S1,S2,S3及びS4を受信することができる。電波S1等は電波の一例である。GPS衛星12a等は、発信源の一例でもある。
電波S1等には各種のコード(符号)が乗せられている。そのうちの一つがC/Aコードである。このC/Aコードは、1,023チップ(chip)から構成されている。そして、このC/Aコードは、1.023Mbpsのビット率、1,023bit(=1msec)のビット長の信号である。このC/Aコードは、測位基礎符号の一例である。そして、端末20は、現在位置を測位する測位装置の一例である。
端末20は、自動車15に搭載されており、自動車15の移動に連れて移動しつつ、現在位置を測位するようになっている。
【0018】
端末20は、例えば、3個以上の異なるGPS衛星12a等からのC/Aコードを受信して、現在位置を測位することができるようになっている。
端末20は、まず、受信したC/AコードがどのGPS衛星に対応するものかを特定する。次に、端末20は、相関処理によって、受信したC/Aコードの位相(以下「コードフェーズ」と呼ぶ)を算出する。続いて、端末20は、そのコードフェーズを使用して、各GPS衛星12a等と端末20との距離(以後、擬似距離と呼ぶ)を算出する。続いて、現在時刻における各GPS衛星12a等の衛星軌道上の位置と、上述の擬似距離に基づいて、現在位置の測位演算を行うことができるように構成されている。
C/Aコードは、電波S1等に乗せられているから、端末20が、電波S1等を受信するときの受信周波数が不正確であると、相関処理によって算出するコードフェーズの精度も劣化する。GPS衛星12a等はその衛星軌道上を移動しているから、この受信周波数は変化し続けているが、電波S1等の信号強度が強い場合には、電波S1等を利用したPLLによって周波数同期を確保し続けることができる。
ところが、電波S1等の信号強度が極めて弱い場合には、PLLは有効に機能しない。また、電波S1等の信号強度が極めて弱い場合には、電波S1等のIFキャリア周波数を的確に予想することも困難である。
この点、端末20は、以下に説明するように、電波S1等の信号強度が極めて弱い場合
に、IFキャリア周波数を予想することなく、精度よく現在位置を測位することができるようになっている。
【0019】
端末20は例えば、携帯電話機、PHS(Personal Handy−phone System)、PDA(Personal Digital Assistance)等であるが、これらに限らない。
また、GPS衛星12a等は、4個に限らず、3個でもよいし、5個以上でもよい。
【0020】
(端末20の主なハードウエア構成について)
図2は、端末20の主なハードウエア構成を示す概略図である。
図2に示すように、端末20は、コンピュータを有し、コンピュータは、バス22を有する。バス22には、CPU(Central Processing Unit)24、記憶装置26等が接続されている。記憶装置26は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)等である。
また、バス22には、外部記憶装置28が接続されている。外部記憶装置28は、例えば、HDD(Hard Disc Drive)である。
また、バス22には、電源装置30、入力装置32、GPS装置34、表示装置47、時計48が接続されている。
【0021】
(GPS装置34の構成について)
図3は、GPS装置34の構成を示す概略図である。
図3に示すように、GPS装置34は、RF部35とベースバンド部36で構成される。
RF部35は、アンテナ35aで電波S1等を受信する。そして、増幅器であるLNA35bが、電波S1に乗せられているC/Aコード等の信号を増幅する。そして、ミキサ35cが、信号の周波数を中間(IF)キャリア周波数にダウンコンバートする。そして、直交(IQ)検波器35dが信号をIQ分離する。続いて、ADコンバータ35e1及び35e2が、IQ分離された信号をそれぞれデジタル信号に変換するように構成されている。
ベースバンド部36は、RF部35からデジタル信号に変換されたIFキャリア周波数の信号を受信する。
ベースバンド部36の相関部37は、入力したデジタル信号を例えば、10ミリ秒(ms)の時間において同期積算し、その積算結果とレプリカC/Aコードとの相関をとる処理であるコヒーレントを行う。相関部37は、NCO38、コードジェネレータ39及び相関器40から構成される。コードジェネレータ39は、NCO38が発生するクロックのタイミングでレプリカC/Aコードを発生する。相関器40はC/AコードとレプリカC/Aコードとの相関をとり、位相の特定及び相関値の算出を行う。相関部37には、周波数及びレプリカC/Aコードの位相を設定することができる。
信号積算器41は、相関部37から出力された相関値を積算する処理であるインコヒーレントを行う。
コードフェーズ検出器42は、相関部37の出力値と信号積算器41からの出力値から、コードフェーズを検出する。
【0022】
上述のように、相関処理は、コヒーレントと、インコヒーレントから構成される。
コヒーレントは、相関部37が、受信したC/AコードとレプリカC/Aコードの相関をとる処理である。
例えば、コヒーレント時間が20msecであれば、20msecの時間において同期積算したC/AコードとレプリカC/Aコードとの相関値等を算出する。コヒーレント処理の結果、相関をとった位相と、相関値が出力される。
インコヒーレントは、コヒーレント結果の相関値を積算することによって、インコヒーレント値を算出する処理である。
相関処理の結果、コヒーレント処理で出力された位相と、インコヒーレント値が出力される。相関値Pはインコヒーレント値である。
電波S1等の信号強度が十分に強い場合には、位相検出器43が相関器40から位相情報を取得して、NCO38に供給し、PLLを構成することができる。この結果、IFキャリア周波数と同期した周波数において、レプリカC/Aコードを発生させることができる。具体的には、相関値Pが最大になるように、受信周波数の制御が行われる。
位相検出器43、相関器40及びNCO38は、受信周波数制御手段の一例である。
【0023】
(端末20の主なソフトウエア構成について)
図4は、端末20の主なソフトウエア構成を示す概略図である。
図4に示すように、端末20は、各部を制御する制御部100、図2のGPS装置34に対応するGPS部102、時計48に対応する計時部104、各種プログラムを格納する第1記憶部110、各種情報を格納する第2記憶部150を有する。
【0024】
図4に示すように、端末20は、第2記憶部150に、衛星軌道情報152を格納している。衛星軌道情報152は、アルマナック152a及びエフェメリス152bを含む。アルマナック152aは、すべてのGPS衛星12a等の概略の軌道を示す情報である。エフェメリス152bは、各GPS衛星12a等の精密な軌道を示す情報である。
端末20は、アルマナック152a及びエフェメリス152bを、測位のために使用する。
【0025】
図4に示すように、端末20は、第2記憶部150に、初期位置情報154を格納している。初期位置情報154は、端末20の現在の初期位置P0を示す情報である。初期位置Q0は、例えば、前回測位時の測位位置である。
【0026】
図4に示すように、端末20は、第1記憶部110に、観測可能衛星算出プログラム112を格納している。観測可能衛星算出プログラム112は、制御部100が、アルマナック152aを参照して、計時部104によって計測した現在時刻において初期位置Q0から観測可能なGPS衛星12a等を示す観測可能衛星情報156を生成するためのプログラムである。
制御部100は、生成した観測可能衛星情報156を第2記憶部150に格納する。
【0027】
図4に示すように、端末20は、第1記憶部110に、第1推定周波数算出プログラム114を格納している。第1推定周波数算出プログラム114は、制御部100が、電波S1等ごとのIFキャリア周波数の予測値である第1推定周波数Aを算出するためのプログラムである。第1推定周波数Aは、端末20が、例えば、現在時刻においてGPS衛星12aからの電波S1を受信するときの電波S1のIFキャリア周波数の予測値である。
【0028】
図5は、第1推定周波数算出プログラム114の説明図である。
図5に示すように、第1推定周波数Aは、発信周波数H1からドップラー偏移H2を加えた周波数である。発信周波数H1は、電波S1等がGPS衛星12a等から発信されるときの周波数である例えば、1.5GHzと、ミキサ35cによるダウンコバート率によって決定される既知の値である。ドップラー偏移H2は、GPS衛星12a等と端末20の相対移動によって生じる周波数偏移であり、常に変動している。ドップラー偏移H2は、端末20の初期位置P0とエフェメリス152bによって、算出することができる。
制御部100は、第1推定周波数Aを示す第1推定周波数情報158を第2記憶部150に格納する。
【0029】
ただし、端末20の位置が正確な現在位置ではなくて初期位置Q0であること、GPS衛星12a等と端末20は常に相対移動を行っていることによって、算出したドップラー偏移H1は真のドップラー偏移とは乖離する可能性がある。
このため、第1推定周波数Aは、真のIFキャリア周波数とは乖離するのが通常である。
【0030】
図4に示すように、端末20は、第1記憶部110に、第1相関プログラム116を格納している。第1相関プログラム116は、制御部100が、GPS衛星12a等から受信したC/AコードとレプリカC/Aコードとの相関値を算出し、さらに、C/Aコードの位相(コードフェーズ)である第1位相CP1を算出するためのプログラムである。
なお、第1位相CP1は、C/Aコードの位相であり、レプリカC/Aコードの位相でもある。
【0031】
図6は、第1相関プログラム116の説明図である。
図6(a)に示すように、制御部100はベースバンド部36によって、C/Aコードの1チップを例えば、等間隔で分割して、相関処理を行う。C/Aコードの1チップは、例えば、32等分される。すなわち、32分の1チップの位相幅(第1位相幅W1)間隔で相関処理を行う。制御部100が相関処理を行うときの第1位相幅W1間隔の位相を第1サンプリング位相SC1と呼ぶ。
第1位相幅W1は、アンテナ35aに入力する信号の信号強度が−155dBm以上である場合に、相関最大値Pmaxを検出することができる位相幅として規定されている。32分の1チップの位相幅であれば、信号強度が−155dBm以上であれば弱電界であっても、相関最大値Pmaxを検出することができることがシミュレーションによって明らかになっている。
【0032】
図6(b)に示すように、制御部100は、推定周波数A±100kHzの範囲を、100Hz単位で変動させながら相関処理を行う。各周波数ごとに、コードフェーズCPを第1位相幅W1で変動させて、相関最大値Pmaxを算出することができる周波数及びコードフェーズを特定する。
なお、測位開始時においては、レプリカC/Aコードを0から1023チップまで変動させる。
そして、一度相関最大値Pmaxに対応するコードフェーズ及び周波数を特定すると、その後は、相関最大値Pmaxに対応するコードフェーズ及び周波数を中心にして、測位開始時よりも狭い範囲において、信号S1等のサーチを行う。例えば、制御部100は、既に算出している第1測位位相CP1を中心に、±256チップの位相範囲をサーチする。また、周波数については、相関最大値Pmaxに対応する周波数を中心に、±1.0kHzの範囲を100Hz単位でサーチする。この条件を、第1トラッキング条件と呼ぶ。
【0033】
図6(c)に示すように、ベースバンド部36からは、2チップ分の位相C1乃至C64に対応する相関値Pが出力される。各位相C1乃至C64が、第1サンプリング位相SC1である。
Pnoiseに対するPmaxの比をSNRと呼ぶ。Pnoiseは環境雑音の信号レベルである。PmaxはGPS衛星12a等からの信号レベルである。
信号S1等の信号強度が弱い状態においては、図6(c)SNR1は、比較的小さい。
【0034】
制御部100は、相関値Pmaxに対応する第1位相CP1を特定する。
制御部100は、第1位相CP1を示す第1位相情報160を第2記憶部位150に格納する。
SNR1が小さいほど、第1位相CP1の精度は低い。
なお、上述の第1相関プログラム116に基づく端末20の動作を、第1相関処理と呼ぶ。
【0035】
図4に示すように、端末20は、第1記憶部110に、第1測位プログラム118を格納している。第1測位プログラム118は、制御部100が、3個以上のGPS衛星12a等に対応する第1位相CP1に基づいて、現在位置を測位して測位位置Q1を算出するためのプログラムである。
【0036】
図7は、測位方法を示す概念図である。
図7に示すように、例えば、GPS衛星12aと端末20との間には、複数のC/Aコードが連続的に並んでいると観念することができる。そして、GPS衛星12aと端末20との間の距離は、C/Aコードの長さの整数倍とは限らないから、コード端数C/Aaが存在する。つまり、GPS衛星12aと端末20との間には、C/Aコードの整数倍の部分(C/Aコードがn個(nは整数)並んでいる部分)と、端数部分(コード端数C/Aa)が存在する。C/Aコードの整数倍の部分とコード端数C/Aaの合計の長さが擬似距離である。端末20は、この擬似距離を使用して測位を行う。
GPS衛星12aの軌道上の位置はエフェメリス152bを使用して算出可能である。そして、GPS衛星12aの軌道上の位置と初期位置Q0との距離を算出すれば、C/Aコードの整数倍の部分を特定することができる。
そして、図7に示すように、レプリカC/Aコードの位相を例えば、矢印X1方向に移動させながら、相関処理を行を行う。
相関値が最大になった位相がコード端数C/Aaである。そして、このコード端数C/Aaが、第1位相CP1である。
【0037】
制御部100は、3個以上のGPS衛星12a等に対応する第1位相CP1に基づいて、各GPS衛星12a等と端末20との擬似距離を算出する。そして、各GPS衛星12a等の軌道上の位置はエフェリス152bによって算出する。そして、3個以上のGPS衛星12a等の軌道上の位置と、擬似距離に基づいて、現在位置を測位して測位位置Q1を算出する。
制御部100は、測位位置Q1を示す第1測位位置情報162を第2記憶部150に格納する。
【0038】
図4に示すように、端末20は、第1記憶部110に、測位位置出力プログラム120を格納している。測位位置出力プログラム120は、制御部100が、測位位置Q1又は後述の測位位置Q2を表示装置47に表示するためのプログラムである。
【0039】
図4に示すように、端末20は、第1記憶部110に、第2相関プログラム122を格納している。第2相関プログラム122は、制御部100が相関処理を行い、相関値Pとコードフェ−ズCPを算出するためのプログラムである。
制御部100は、相関値Pとコードフェ−ズCP示す第2相関情報164を第2記憶部150に格納する。
【0040】
図4に示すように、端末20は、第1記憶部110に、ピーク周波数特定プログラム124を格納している。ピーク周波数特定プログラム124と制御部100は、ピーク周波数特定手段の一例である。
【0041】
図8は、ピーク周波数特定プログラム124の説明図である。
図8に示すように、制御部100は、相関最大値Pmaxに対応する周波数をピーク周波数F0として特定する。ピーク周波数F0はピーク周波数の一例である。
このピーク周波数F0は、端末20が例えば、100Hzの幅でサーチを行った結果であるから、受信した電波S1等の真のIFキャリア周波数と、最大約50Hzの乖離を生じている。
制御部100は、ピーク周波数F0を示すピーク周波数情報166を第2記憶部150に格納する。
【0042】
図4に示すように、端末20は、第1記憶部110に、参照周波数算出プログラム126を格納している。参照周波数算出プログラム126と制御部100は、参照周波数算出手段の一例である。
制御部100は参照周波数算出プログラム126に基づいて、ピーク周波数F0よりも100Hz低い周波数F1と、ピーク周波数F0よりも100Hz高い周波数F2を算出する。制御部100は、周波数F1及び周波数F2を示す参照周波数情報168を第2記憶部150に格納する。周波数F1は低周波数の一例である。周波数F2は高周波数の一例である。
ピーク周波数F0と周波数F1との周波数差分と、ピーク周波数F0と周波数F2との周波数差分が等しくなるように、周波数F1と周波数F1が規定される。本実施の形態では、周波数差分は100Hzに設定されている。
なお、本実施の形態とは異なり、周波数差分は100Hzに限らない。
【0043】
図4に示すように、端末20は、第1記憶部110に、参照相関値算出プログラム128を格納している。参照相関値算出プログラム128と制御部100は、参照相関値算出手段の一例である。
制御部100は参照相関値算出プログラム128に基づいて、周波数F1に対応する相関値P1と、周波数F2に対応する相関値P2を算出する。具体的には、制御部100は、第2相関情報164を参照して、相関値P1及び相関値P2を算出する。
制御部100は、相関値P1及び相関値P2を示す参照相関値情報170を第2記憶部150に格納する。
【0044】
図4に示すように、端末20は、第1記憶部110に、第2推定周波数算出プログラム130を格納している。第2推定周波数算出プログラム130は、制御部100が、ピーク周波数F0と相関ピーク値Pmax(P0)、周波数F1と相関値P1、及び、周波数F2と相関値P2とに基づいて、第2推定周波数Frを算出するためのプログラムである。第2推定周波数Frは、補正後ピーク周波数の一例である。第2推定周波数算出プログラム130と制御部100は、補正後ピーク周波数算出手段の一例である。
【0045】
図9及び図10は、第2推定周波数算出プログラム130の説明図である。
相関値Pと周波数Fを示すグラフは、図9及び図10に示すように、二等辺三角形を描く。
図9(a)及び図10(a)に示すように、ピーク周波数F0と相関ピーク値P0によって点G0が規定される。周波数F1と相関値P1によって点G1が規定される。そして、周波数F2と相関値P2とによって点G2が規定される。
【0046】
図9(a)及び図9(b)に示すように、相関値P1が相関値P2よりも小さい場合には、点G0と点G1は傾きa(aは正の数)の同一直線上にある。点G0と点G1を結んだ直線が直線L1である。
そして、点G2は、傾きが−aの直線上にある。傾きが−aで、かつ、点G2を通る直線が直線L2である。
そして、直線L1と直線L2の交点が二等辺三角形の頂点Hである。頂点Hに対応する周波数が、第2推定周波数Frである。図9(b)の連立方程式1を解くと、未知数Fr、Pr及び傾きaを算出することができる。
【0047】
図10(a)及び図10(b)に示すように、相関値P1が相関値P2よりも大きい場合には、点G0と点G2は傾き−a(aは正の数)の同一直線上にある。点G0と点G2を結んだ直線が直線L2である。
そして、点G1は、傾きがaの直線上にある。傾きがaで、かつ、点G1を通る直線が直線L1である。
そして、直線L1と直線L2の交点が二等辺三角形の頂点Hである。頂点Hに対応する周波数が、第2推定周波数Frである。図10(b)の連立方程式2を解くと、未知数Fr、Pr及び傾きaを算出することができる。
【0048】
なお、相関値P1が相関値P2と等しい場合には、ピーク周波数F0が第2推定周波数Frである。
制御部100は、第2推定周波数Frを示す第2推定周波数情報172を第2記憶部150に格納する。
この第2推定周波数Frは、周波数Fのサーチステップである100Hzの制限を受けないから、精度の高い情報である。すなわち、ピーク周波数F0よりも、真のIFキャリア周波数に近い。
【0049】
図4に示すように、端末20は、第1記憶部110に、第2位相決定プログラム132を格納している。第2位相決定プログラム132は、制御部100が、第2推定周波数Frを使用して、電波S1等を受信して、相関処理を行い、測位を行うための第2位相CP2を算出するためのプログラムである。第2位相決定プログラム132と制御部100は、電波受信手段の一例である。
【0050】
図11は、第2位相決定プログラム132の説明図である。
図11の相関グラフにおけるSNR2は、図6(c)のグラフにおけるSNR1よりも大きい。これは、第2推定周波数Frが、真のIFキャリア周波数に非常に近いからである。
このため、相関最大値Pmaxに対応する位相である第2位相CP2は精度の高い位相情報である。
制御部100は、第2位相CP2を示す第2位相情報174を第2記憶部150に格納する。
上述の第2相関プログラム122、ピーク周波数特定プログラム124、参照周波数算出プログラム126、参照相関値算出プログラム128及び第2推定周波数算出プログラム130、第2位相決定プログラム132に基づく端末20の動作を、第2相関処理と呼ぶ。
【0051】
図4に示すように、端末20は、第1記憶部110に、第2測位プログラム134を格納している。第2測位プログラム134は、制御部100が、3個以上のGPS衛星12aについての第2位相CP2を使用して、測位を行い、測位位置Q2を算出するためのプログラムである。
制御部100は、測位位置Q2を示す第2測位位置情報176を第2記憶部150に格納する。
【0052】
図4に示すように、端末20は、第1記憶部110に、信号強度評価プログラム136を格納している。
信号強度評価プログラム136は、アンテナ35aに入力する信号の信号強度SPを評価するためのプログラムである。アンテナ35aに入力する信号の信号強度SPは、相関値から推定することができる。
制御部100は、例えば、信号強度SPが−138dBm以上である場合には、第1相関処理を行って、測位位置Q1を算出する。
そして、制御部100は、信号強度SPが−142dBm以下である場合には、第2相関処理を行って、測位位置Q2を算出する。
そして、制御部100は、信号強度SPが−142dBmより大きく−138dBm未満である場合には、第1相関処理と第2相関処理を平行して実施する。そして、制御部100は、第1位相CP1を使用して、測位位置Q1を算出する。
【0053】
端末20は、以上のように構成されている。
上述のように、端末20は、ピーク周波数F0(図4参照)を特定することができる。
また、端末20は、第2推定周波数Fr(図4参照)を算出することができる。
レプリカC/Aコードの位相を固定した場合、相関値と受信周波数(IFキャリア周波数)との関係を示すグラフは、図9に示すように、相関値の最大値に対応する点を頂点とする二等辺三角形を描く。そして、ピーク周波数F0に対応する点G0はその頂点H付近に位置し、ピーク周波数F0の前後の周波数F1及びF2にそれぞれ対応する点G1及びG2はそれぞれ異なる斜辺に位置する。そして、点G1及び点G2のいずれか一方は、点G1と同一の斜辺に位置するから、斜辺の傾きaを特定することができる。二等辺三角形においては、一方の斜辺の傾きを特定することができれば、他方の斜辺の傾きも特定することができる。そして、二つの斜辺が交わった点が頂点Hである。そしてこの頂点Hに対応する周波数が、上述の第2推定周波数Frである。
上述のように、電波S1等の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることができない場合であっても、ピーク周波数F0は必ず1つ存在する。そして、ピーク周波数F0を特定すれば、第2推定周波数Frを算出することができる。
そして、端末20は、第2推定周波数Frを使用して、電波S1等を受信することができる。このため、精度よく相関値Pを算出することができ、現在位置を精度よく算出することができる。
これにより、衛星電波の信号強度が極めて弱い場合において、予想IFキャリア周波数を定めることを要せずに、精度良く測位することができる。
【0054】
また、端末20は、PLLによって、レプリカC/Aコードと受信したC/Aコードとのコヒーレント値が最大になるように受信周波数の制御を行うことができる。
これにより、電波S1等の信号強度が所定の強度範囲の場合に、PLLを有効に機能させ、受信周波数を継続的に電波S1等のIFキャリア周波数に近づけることができる。
【0055】
また、端末20は、電波S1等の信号強度が所定範囲の場合には、上述の第1相関処理と第2相関処理を平行して行うことができる。このため、信号強度SPが所定の強度よりも大きい状態からから小さい状態に移行した場合に、連続的に精度良く測位を行うことができる。
【0056】
以上が本実施の形態に係る端末20の構成であるが、以下、その動作例を主に図12及び図13を使用して説明する。
図12及び図13は端末20の動作例を示す概略フローチャートである。
【0057】
まず、端末20は、各GPS衛星12a等について、エフェメリス152bと初期位置Q0から、推定周波数Aを算出する(図12のステップST1)。
続いて、端末20は、第1相関処理を行う(ステップST2)。
【0058】
続いて、端末20は、信号強度SPを判断する(ステップST3)。
端末20は、ステップST3において、信号強度SPが−138dBm以上であると判断すると、第1相関処理を継続し(ステップST4A)、第1位相CP1を使用して現在位置を測位して、測位位置Q1を算出する(ステップST5A)。
続いて、端末20は、測位位置Q1を出力する(ステップST6A)。
続いて、端末20は、測位が測位規定回数である例えば、10回に達したか否かを判断する(ステップST7)。
端末20は、測位が測位規定回数に達したと判断した場合には測位を終了する。
端末20は、測位が測位規定回数に達していない判断した場合にはステップST3以降を実施する。
【0059】
端末20は、ステップST3において、信号強度SPが−142dBm以下であると判断すると、第1相関処理を停止し、第2相関処理を行う(ステップST4B)。
第2相関処理においては、端末20は、まず、ピーク周波数F0(図4参照)を特定する(図13のステップST101)。このステップST101は、ピーク周波数特定ステップの一例である。
続いて、端末20は、周波数F1及びF2(図4参照)を算出する(ステップST102)。このステップST102は、参照周波数算出ステップの一例である。
【0060】
続いて、端末20は、相関値P2及びP3(図4参照)を算出する(ステップST103)。このステップST103は、参照相関値算出ステップの一例である。
続いて、端末20は、第2推定周波数Fr(図4参照)を算出する(ステップST104)。このステップST104は、補正後ピーク周波数算出ステップの一例である。
【0061】
続いて、端末20は、第2位相CP2を算出し、その第2位相CP2を使用して現在位置を測位して、測位位置Q2を算出する(図12のステップST5B)。
続いて、端末20は、測位位置Q2を出力し(ステップST6B)、ステップST7を実施する。
【0062】
端末20は、ステップST3において、信号強度SPが−142dBmより大きく−138dBm未満であると判断すると、第1相関処理と第2相関処理を平行実施する(図12のステップST4C)。
続いて、端末20は、第1位相CP1を使用して現在位置を測位して、測位位置Q1を算出する(ステップST5C)。
続いて、端末20は、測位位置Q1を出力し(ステップST6C)、ステップST7を実施する。
端末20は、測位が測位規定回数に達していない判断した場合にはステップST3以降を実施する。再度のステップST3において、端末20は、信号強度SPが−138dBm以下であると判断すると、ステップST4Bに進む。ここで、第1相関処理と第2相関処理が平行して継続しているから、第1相関処理を停止し、直ちに第2相関処理を実施することができる。これは、第1相関処理においてPLLが機能しなくなった後に、第2相関処理を開始するのではなくて、信号強度SPが−142dBm以下に低下するおそれがある中程度の状態(信号強度SPが−142dBmより大きく−138dBm未満の状態)で、第2相関処理を継続しておくことを意味する。このため、第2相関処理において、新たに広範囲の周波数及び位相をサーチする必要がないから、迅速にステップST5B以下を実施することができる。
また、中程度の状態(信号強度SPが−142dBmより大きく−138dBm未満の状態)は、信号強度SPが−138dBm以上になる可能性がある状態でもある。予め第1相関処理を継続しておくことによって、信号強度が−138dBm以上になった場合に、直ちに第1相関処理だけを実施する状態に移行することができる。
【0063】
本発明は、上述の各実施の形態に限定されない。さらに、上述の各実施の形態は、相互に組み合わせて構成するようにしてもよい。
【図面の簡単な説明】
【0064】
【図1】本発明の実施の形態の端末等を示す概略図である。
【図2】端末の主なハードウエア構成を示す概略図である。
【図3】GPS装置の構成を示す概略図である。
【図4】端末の主なソフトウエア構成を示す概略図である。
【図5】第1推定周波数算出プログラムの説明図である。
【図6】第1相関プログムの説明図である。
【図7】測位方法を示す概念図である。
【図8】ピーク周波数特定プログラムの説明図である。
【図9】第2推定周波数算出プログラムの説明図である。
【図10】第2推定周波数算出プログラムの説明図である。
【図11】第2位相決定プログラムの説明図である。
【図12】端末の動作例を示す概略フローチャートである。
【図13】端末の動作例を示す概略フローチャートである。
【符号の説明】
【0065】
12a,12b,12c,12d・・・GPS衛星、20・・・端末、34・・・GPS装置、112・・・観測可能衛星算出プログラム、114・・・第1推定周波数算出プログラム、116・・・第1相関プログラム、118・・・第1測位プログラム、120・・・測位位置出力プログラム、122・・・第2相関プログラム、124・・・ピーク周波数特定プログラム、126・・・参照周波数算出プログラム、128・・・参照相関値算出プログラム、130・・・第2推定周波数算出プログラム、132・・・第2位相特定プログラム、134・・・第2測位プログラム、136・・・信号強度評価プログラム
【特許請求の範囲】
【請求項1】
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置であって、
前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号の相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定手段と、
前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出手段と、
前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出手段と、
前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出手段と、
前記補正後ピーク周波数を使用して、前記電波を受信する電波受信手段と、
を有することを特徴とする測位装置。
【請求項2】
前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行う受信周波数制御手段を有することを特徴とする請求項1に記載の測位装置。
【請求項3】
前記補正後ピーク周波数算出手段と前記受信周波数制御手段とが平行して作動することを特徴とする請求項2に記載の測位装置。
【請求項4】
前記発信源は、測位衛星であることを特徴とする請求項1乃至請求項3のいずれかに記載の測位装置。
【請求項5】
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、
前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、
前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、
前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、
前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、
を有することを特徴とする測位装置の制御方法。
【請求項6】
コンピュータに、
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、
前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、
前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、
前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、
前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、
を実行させることを特徴とする測位装置の制御プログラム。
【請求項7】
コンピュータに、
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、
前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、
前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、
前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、
前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、
を実行させることを特徴とする測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項1】
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置であって、
前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号の相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定手段と、
前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出手段と、
前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出手段と、
前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出手段と、
前記補正後ピーク周波数を使用して、前記電波を受信する電波受信手段と、
を有することを特徴とする測位装置。
【請求項2】
前記レプリカ測位基礎符号と前記測位基礎符号とのコヒーレント値が最大になるように受信周波数の制御を行う受信周波数制御手段を有することを特徴とする請求項1に記載の測位装置。
【請求項3】
前記補正後ピーク周波数算出手段と前記受信周波数制御手段とが平行して作動することを特徴とする請求項2に記載の測位装置。
【請求項4】
前記発信源は、測位衛星であることを特徴とする請求項1乃至請求項3のいずれかに記載の測位装置。
【請求項5】
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、
前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、
前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、
前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、
前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、
を有することを特徴とする測位装置の制御方法。
【請求項6】
コンピュータに、
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、
前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、
前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、
前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、
前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、
を実行させることを特徴とする測位装置の制御プログラム。
【請求項7】
コンピュータに、
複数の発信源からの電波に乗せられている測位基礎符号を使用して、現在位置を測位する測位装置が、前記測位装置が発生するレプリカ測位基礎符号と前記測位基礎符号との相関値の最大値に対応する受信周波数であるピーク周波数を特定するピーク周波数特定ステップと、
前記測位装置が、前記ピーク周波数よりも低い周波数である低周波数と、前記ピーク周波数よりも高い周波数である高周波数を算出する参照周波数算出ステップと、
前記測位装置が、前記低周波数に対応する前記相関値と、前記高周波数に対応する前記相関値を算出する参照相関値算出ステップと、
前記測位装置が、前記ピーク周波数に対応する前記相関値及び前記ピーク周波数と、前記低周波数に対応する前記相関値及び前記低周波数と、前記高周波数に対応する前記相関値及び前記高周波数とに基づいて、補正後ピーク周波数を算出する補正後ピーク周波数算出ステップと、
前記測位装置が、前記補正後ピーク周波数を使用して、前記電波を受信する電波受信ステップと、
を実行させることを特徴とする測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【公開番号】特開2007−256111(P2007−256111A)
【公開日】平成19年10月4日(2007.10.4)
【国際特許分類】
【出願番号】特願2006−81533(P2006−81533)
【出願日】平成18年3月23日(2006.3.23)
【出願人】(000002369)セイコーエプソン株式会社 (51,324)
【Fターム(参考)】
【公開日】平成19年10月4日(2007.10.4)
【国際特許分類】
【出願日】平成18年3月23日(2006.3.23)
【出願人】(000002369)セイコーエプソン株式会社 (51,324)
【Fターム(参考)】
[ Back to top ]