Comments
Description
Transcript
スードライトの基礎知識 pdfファイル版
スードライトの基礎知識 INDOOR POSITIONINGの可能性 平成15年3月12日 測位衛星技術株式会社 スードライト(Pseudolites:PL)とは • • • • 地上に設置する“疑似的(Pseudo)な衛星(Satellite)” 1970年代GPS衛星打上前の地上テストにおいて使用 GPS衛星信号と互換性のある信号を発生 基本的にGPSと同じ周波数を使用(L1:1575.42MHz) • 衛星補強システムの1つとして期待 ⇒RTCM SC104でも信号仕様を定義 • ユーザ受信機に改良が必要 背景の確認 ●ポジショニングの多様化 ・カーナビゲーションの完全普及 ・携帯電話の測位機能(GPSOne,PHSによるポジショニングなど) ・ロケーションサービスの開始(ココセコム,イマドコサービス他) ●衛星測位環境の充実 ・GPS:第3周波数L5の追加、静止衛星からのGPS信号発信(オーバーレイ) ・2004年から欧州によるガリレオシステムSTART ・2005年ロシアGLONASSシステム完成 ●MTSATからWAAS信号発信 ●国土地理院の電子基準局⇒VRSサービス開始準備 ●ITRF測地系への統一(GPS利用の座標系) 【ポジショニングニーズ】 室内から屋外でのシームレスなポジショニング技術の確立 室内から屋外でのシームレスなポジショニング技術の確立 1.携帯電話などを対象としたシームレスナビゲーション⇒メートルレベル 2.INDOORにおける高精度ポジショニング ⇒センチメートルレベル Pseudolite技術への期待 実験システム系統図 GPS衛星 GPS信号L1 GPSアンテナ GPSアンテナ IN400 スードライト スードライト信号L1 GPSアンテナ 特定小電力無線機 RTCM ディファレンシャル補正値 スードライト対応 GPS受信機 IN2400 DGPS基地局 IN400 スードライト スードライトの利点 1. スードライトにより、GPS衛星の数が“追加” =>衛星数の増加による信頼性の向上 2. スードライト配置によるDOPの向上 =>測位精度の改善 3. 移動体における急激なジオメトリーの変化 =>RTKにおける高速ambiguity決定が可能 ・α=90°で正確な ambiguityの決定が可能 ・αが大きくなるにつれて 遠近比が問題となる Pseudolite 利用イメージ 室内外シームレスな測位 都市部でも測位可能! ITS 自動走行 歩行者ITS バリアフリー ロボット制御 実用化へのプロセス インテグリティーノーティクス社代理権確保 認証されたスードライト開発 ●インフラ設備 ●インフラ設備 実用評価試験 独自開発 普及機生産 ●端末機 実用評価実験 評価受信機 小型チップ 2003∼2004 携帯電話他情報端末への挿入 次世代携帯電話 ボード受信機 Pseudolite & RTK 先行事例∼瞬間FIX∼ Outrider TUAV スタンフォード大学 室内におけるスードライト研究例 スードライト発信機 プレートの裏側 ・商用6チャンネルGPS受信機を改良 ・4つのアンテナ ・19.2kbpsRFモデム Robot & Target Vehicles 【スードライト】 ・GPS衛星として6個のスードライトを配置 ・1つのスードライトをマスタとして、タイミングデータ を受信機の同期のために使用。 その他のアプリケーション • 露天掘り等における擬似衛星の利用 • その他GPS衛星の確保が難しい場所にお ける擬似衛星による補完 基礎検証実験概要 スードライトの配置(位置関係や距離)、設置場所、マルチパス、 遠近問題などによってどのように利用可能性や測位精度が影響されるかについて基礎的な実験を行う。 その結果をもってPL設置に関する各種パラメータのガイダンスを策定する。 利用可能性・精度の検証 基礎データの取得 室内での利用を 想定した際の課題 特に室内での利用に際しての課題(時刻同期、補正データ 配信、アンビギュイティの解決など)に関して、これらをクリアす るための基礎実験。 電波伝播特性に関する 基礎データ収集 利用可能範囲や測位精度の推定シミュレーションにおいて必要な GPS、およびスードライトの電波伝播特性(減衰、反射、吸収、透過 など)に関する基礎データを収集するための基礎実験。 遠近問題(Near/Far Problem) 1 10 • 地上におけるGPS信号レベルはほぼ一定(-130dBm) • 信号レベルは距離の2乗に逆比例 近い(Near)ゾーン⇒ PL信号がGPS信号を妨害 遠い( Far )ゾーン⇒ PL信号は受信不能 • 連続的な信号(非パルス方式)⇒遠近比は約10:1 遠近問題への対応法 (下記の複数を同時に採用) • TDMA方式. 方式 ランダムまたは特定のサイクルレートでパルス信号を 生成し出力。 • FDMA方式.GPSのL1帯から周波数をオフセットさせる。 L:1575.42MHz±1.023MHzを用いる。比較的最近の 技術で受信機側でより多くの改修が必要。 • CDMA方式. 異なるゴールドコードを利用. GPSのC/Aコード1023種 類の中でGPS、WAAS等で使用されているものを除く 712種類の中から利用。 パルス方式の基本概念(パルス duty cycleの考え方) PL2 PL2 PL1 PL2 PL2 PL1 PL2 PL1 PL2 24chips 1ms 1023chips ● PL1のパルスDuty⇒7.5% = 十分受信できるDutyサイクル 2.5%×3 ● 1コード長(1023chips)当たりに3パルスが望ましい。 ● 例 PL×4個 ⇒ 30% ● C/Aコード(1023chips)のコードパターンは、変わらない! パルス方式の基本概念 ・なぜ、GPS信号が受信できるのか? ・なぜ、PL(パルス)でコード判読できるのか? 1パルス=24chips 101001101 2.5%×3=7%でのパルス送信! この部分は、 GPS信号はできない! 1024chips 通常10スロット(RTCMで規定) GPS受信機では、10スロットのパターンを受け取る中で コードパタンを一致させる。 GPS信号妨害への対処法 • パルス方式 • コードは時間により拡散されている(TDMA) PL有効ゾーン 重なるPLゾーン数 <10~15でなければなら ない オーバーラップするPLゾーン数には限界あり カバレッジを犠牲にせず、 オーバーラップするPL数を制限するためには… • Nearゾーンを物理的に隔離 →高層ビルやタワーの屋上に設置 • セル分割方式 (携帯電話と同様) それぞれのセルに ~10台のPLを配置 スードライト基礎干渉実験 【目的】 スードライトの基礎検証実験を行なうにあたって、スードライトが最大出力で送出した 場合においても、実験エリア外の一般GPSユーザーに影響を及ぼさない安全 距離の算出と検証を行う。 敷地外 安全距離 実験エリア GPS信号強度 A.通常電界⇒−110dBm B.弱電界 ⇒−130dBm アッテネーター α スードライト出力 MAX;+3dBm スードライト 信号合成器 β 安全距離 ユーザーGPS a.測量用受信機 b. 普及型受信機 【実験手順】 境界線 1.GPS信号強度をA.Bそれぞれの入力レベルに対して ⇒アッテネータβを減衰量を段階的に小さくし、ユーザーGPSが干渉した 時のアッテネータ−量(減衰量)を測定する。 2.連続送出とパルス送出の両方式について測定を行なう。 3.その時の減衰量からスードライト最大出力時の離れるべき距離を算出する。 同軸伝送路実験概要 屋上GPSアンテナ スードライト出力 MAX;+3dBm 固定減衰器 信号合成器 −4dB −7dB GPS信号強度 A.通常電界⇒−110dBm B.弱電界 ⇒−130dBm 4分配器 2分配器 可変アッテネータ スードライト スペクトラム アナライザー 測量GPS受信機 一般GPS受信機 −7dB 0dB ∼ −100dB (安全距離の算出) 自由空間伝播損失=20log(4πD/λ) 測定した各自由空間伝播損失を 上記値を計算式に当てはめ、安全距離Dを算出する。 測量用GPS受信機:アシュッテック Z−Survyor2周波受信機 一般GPS受信機 :ガーミン社 GPSⅢ一周波受信機 (入力GPS信号強度について) 通常GPS受信機は、−110dBm∼−120dBmで受信できるように設計されて いる。従って、今回の実験では 通常受信を−110dBm、弱い受信感度時を−130dBmとして実験を行なった。 パルスDutyの状況 Duty100%連続波(1024chips) Duty10%(102chips) Duty20%(205chips) Duty5%(25chips) 実験機器 1 IN400スードライト IN1200 スードライト対応1周波受信機 2周波測量機 マゼラン社 Z−surveyor ガーミン GPSⅢ 1周波受信機 実際に設置を仮定した時のPL送信出力 −52.3dBm PL1 −62.7dBm PL2 12m 40m 116m 35m −130dBm 20° 2周波干渉強度 仮定出力PL1 干渉距離 仮定出力PL2 干渉距離 連続波 -97.3dBm -52.3dBm 2.7m -62.71dBm 0.82m 20%パルス -95.7dBm -52.3dBm 2.25m -62.71dBm 0.67m 10%パルス -83.0dBm -52.3dBm 0.5m -62.71dBm 0.16m 2.5%パルス 干渉なし -52.3dBm 対象外 -62.71dBm 対象外 GPS/PL テスト結果 • パルス方式による遠近比の向上 : – PLから18km離れた航空機においてPL信号を問題なく受信 – PLから18cmの距離でハンドヘルド受信機はGPS信号を受信 (IntegriNautics社による実験) ⇒遠近比:100,000:1 実験機器構成 GPS/PL用基準局 GPS用基準局 GPS L1/L2アンテナ GPS L1アンテナ PL基準局 IN2400BAS GPS Diff 基準局 Z-12 RTCM Type1 SS無線 (GPS/PL補正用) 特小無線 (GPS DGPS補正用) 2分配器 評価用受信機群(移動局) GPS RTK基準局 Z-12 DBEN GPS L1/L2アンテナ 特小無線 (GPS RTK補正用) 4分配器 GPS/PL 1周波 受信機 スードライト対応 ( IN2400Rov) 表示・記録用PC SS無線 (GPS/PL補正用) 特小無線 (GPS RTK補正用) 特小無線 (GPS DGPS補正用) GPS 2周波 搬送波受信機 測量・高精度用 (Ashtech Z-12Extreme) 表示・記録用PC GPS 1周波 搬送波受信機 高精度ナビ・コントロール (Ashtech G-12) 表示・記録用PC GPS 1周波 コード受信機 カーナビゲーション用 (Ashtech G-8) 表示・記録用PC Z12 (RTK-GPS) PL出力小、PL全て使用 水平方向測位結果 計測値(fix解) 計測値(float解) 真値 スタート地点 25 20 PRN 36 15 PRN 33 BAS PRN 地上PseudoLite 屋上PseudoLite 40 10 PRN 35 5 y[m] (北向きが正) PRN 39 0 -5 -10 -15 -20 -25 -30 -35 -40 -45 -30 -25 -20 -15 -10 -5 x[m] (東向きが正) 0 5 10 IN2400 (Pseudolite) PL出力小、PL屋上2個と40番を使用 水平方向測位結果 IN2400 (Pseudolite) PL出力小、PL地上3個を使用 25 計測値 真値 スタート地点 鉛直方向測位結果 0.8 計測値 真値 0.6 20 高さ[m] 0.4 15 0.2 0.0 10 -0.2 計測値の真値に対する残さ(水平方向は真値からの距離) 0.40 5 水平方向 鉛直方向 0.30 0.20 0 残さ[m] 0.10 -0.20 rms(水平):0.0417m rms(鉛直):0.1154m -0.30 -0.40 DOP 8 HDOP VDOP 7 6 y[m] (北向きが正) 0.00 -0.10 -5 -10 -15 DOP 5 4 -20 3 2 1 -25 0 全衛星数 GPS衛星数 PL数 衛星数 12 10 -30 衛星数 8 6 -35 4 2 -40 0 -20 0 20 40 60 80 100 120 140 160 180 水平方向移動距離[m] -45 -30 -25 -20 -15 -10 -5 x[m] (東向きが正) 0 5 10 実験会場(Θ館) 実験構成(PL送信アンテナ) PRN34 PRN34 PRN33 PRN33 PRN35 PRN36 PRN36 PRN35 実験構成(リピータ) 屋外GPS 受信アンテナ リピータ 送信アンテナ リピータ 屋外GPS 実験構成(BASE・ROV) ROV位置(P41) BASE P41 P5 BASE設置位置 ROV位置(P5) 屋内実験 • GPS信号の再放射をする前提上、受信機時刻の誤差を 補正する技術を利用する。 • 開発したスードライトシステムは、2つのスードライトを利 用してコードによる2次元的な屋外測位を可能にしている。 • ユーザーの高さと屋外アンテナの座標は補正情報(メッ セージ?)の一部として伝達される。 • 開発されたシステムはGSSシングルチャンネルシュミレー ターを利用しての検証もなされた。 • ソフトは同梱のコンピューターのリアルタイムモードとポス トプロセシングモードで実行された。 2つのスードライトのみで受信機の位置を割り出す等式 • ( XRA - XS )2+ (YRA -YS )2 +(ZRA- ZS )2 + TL + ( X - XRT )2+ (Y -YRT )2 +(Z- ZRT )2 + TR + AtmErr +TS = P02 • ( X – X1 )2+ (Y –Y1 )2 +(Z- Z1 )2 + TR + TPL1 = P1 2 • ( X – X2 )2+ (Y –Y2 )2 +(Z- Z2 )2 + TR + TPL2 = P2 2 • Z=H TPL2 ,TPL2 , AtmErr, TS - removed by DD X,Y,Z,TR - ? 屋内測位における問題点 • マルチパス (より好条件の設置場所の為、実験結果はルーバーで20-25mの誤差、 参考用受信機でもメーターレベルの誤差が生じた。) 2重解 (GPS衛星では起こり得ない、2つの解を等式は提示する。) • • 初期座標 – 非線形の擬似距離方程式を真値近くで線形化する必要がある。 – PLまでの距離は衛星までのそれより遥かに近くにあり、従って、初期座標誤差に よる方向コサイン誤差の方が重要である。 [δ X]= ([A]T[A]) -1 [A]T[δr] GPS Satellite PL <- 20 km -> True position True position <- 20 m -> Linearization point Linearization point - 既存の解決方法: 基準点から測位を開始する – 実生活の用途に不適当 -基準点が比較的重要でない測位計算方法が別途必要である。分析的解決案 に基づく特別な演算方法を利用した。 Example of the Analytical Solution, where TR = 0, Z = H • • • • • • x1 = [ -1/2*(-1/(-8*p2_y*p1_y-8*p1_x*p2_x+4*p2_y^2+4*p1_y^2+4*p2_x^2+4*p1_x^2)*(-4*p2_y^2*p1_y-4*p2_y*d2+4*p2_y^3-8*p1_x*p2_x*p2_y+4*p2_y*p1_x^2-4*p1_y*d1+4*p2_x^2*p2_y+4*p2_x^2*p1_y8*p1_x*p2_x*p1_y-4*p2_y*p1_y^2+4*p1_x^2*p1_y+4*p2_y*d1+4*d2*p1_y+4*p1_y^3+4*(12*p1_x^2*p2_x^2*d2+8*p2_y*d2*p1_x*p2_x*p1_y-8*p2_y^3*p1_y*p1_x*p2_x+12*p2_y^2*p1_y^2*p1_x*p2_x4*p2_y^2*d2*p1_x*p2_x-4*p2_y*d2*p2_x^2*p1_y-4*p2_y*d2*p1_x^2*p1_y-16*p1_x*p2_x^3*p2_y*p1_y+24*p1_x^2*p2_x^2*p2_y*p1_y-16*p1_x^3*p2_x*p2_y*p1_y-4*p1_x*p2_x*p2_y^2*d18*p1_x*p2_x*p2_y*p1_y^3-4*p2_y*p1_x^2*p1_y*d1-4*p1_y*d1*p2_x^2*p2_y-4*p1_y^2*d1*p1_x*p2_x-4*p1_x*p2_x*p1_y^2*d2-4*p1_x*p2_x*d2*d1+4*p2_y^3*p1_y*p1_x^2+4*p2_y^3*p1_y*p2_x^26*p2_y^2*p1_y^2*p2_x^2-6*p2_y^2*p1_y^2*p1_x^2+2*p2_y^2*d2*p1_x^2+2*p2_y^2*d2*p2_x^2+2*p2_y^4*p1_x*p2_x12*p1_x^2*p2_x^2*p2_y^2+8*p1_x^3*p2_x*p2_y^2+8*p1_x*p2_x^3*p2_y^2+4*p2_y*p1_x^4*p1_y+2*p2_y^2*p1_x^2*d1+4*p2_y*p1_x^2*p1_y^3+2*p1_y^2*d1*p2_x^2+2*p1_y^2*d1*p1_x^2+4*p2_x^4*p2_y*p1_y +2*p2_x^2*p2_y^2*d1+4*p2_x^2*p2_y*p1_y^3+8*p2_x^3*p1_y^2*p1_x12*p2_x^2*p1_y^2*p1_x^2+2*p2_x^2*p1_y^2*d2+8*p1_x^3*p2_x*p1_y^2+2*p1_x*p2_x*p1_y^4+2*p1_x^2*p1_y^2*d2+12*p1_x^2*p2_x^2*d1-8*p1_x*p2_x^3*d2-8*p1_x*p2_x^3*d1-p2_y^4*p1_x^2p2_y^4*p2_x^2-2*p2_y^2*p1_x^4-2*p2_x^4*p2_y^2-2*p2_x^4*p1_y^2-p2_x^2*p1_y^4-2*p1_x^4*p1_y^2-p1_x^2*p1_y^4+20*p1_x^3*p2_x^3-15*p1_x^2*p2_x^415*p1_x^4*p2_x^2+6*p1_x*p2_x^5+6*p1_x^5*p2_x+2*p2_x^4*d2+2*p2_x^4*d1-p2_x^2*d2^2-p2_x^2*d1^2-p1_x^2*d2^2+2*p1_x^4*d2+2*p1_x^4*d1-p1_x^2*d1^2-p2_x^6p1_x^6+8*p1_x*p2_x*p2_y*p1_y*d1+2*p1_x*p2_x*d2^2-8*p1_x^3*p2_x*d2-8*p1_x^3*p2_x*d1+2*p1_x*p2_x*d1^2+2*p2_x^2*d2*d1+2*p1_x^2*d2*d1)^(1/2))*p2_y+1/(-8*p2_y*p1_y8*p1_x*p2_x+4*p2_y^2+4*p1_y^2+4*p2_x^2+4*p1_x^2)*(-4*p2_y^2*p1_y-4*p2_y*d2+4*p2_y^3-8*p1_x*p2_x*p2_y+4*p2_y*p1_x^2-4*p1_y*d1+4*p2_x^2*p2_y+4*p2_x^2*p1_y-8*p1_x*p2_x*p1_y4*p2_y*p1_y^2+4*p1_x^2*p1_y+4*p2_y*d1+4*d2*p1_y+4*p1_y^3+4*(12*p1_x^2*p2_x^2*d2+8*p2_y*d2*p1_x*p2_x*p1_y-8*p2_y^3*p1_y*p1_x*p2_x+12*p2_y^2*p1_y^2*p1_x*p2_x-4*p2_y^2*d2*p1_x*p2_x4*p2_y*d2*p2_x^2*p1_y-4*p2_y*d2*p1_x^2*p1_y-16*p1_x*p2_x^3*p2_y*p1_y+24*p1_x^2*p2_x^2*p2_y*p1_y-16*p1_x^3*p2_x*p2_y*p1_y-4*p1_x*p2_x*p2_y^2*d1-8*p1_x*p2_x*p2_y*p1_y^34*p2_y*p1_x^2*p1_y*d1-4*p1_y*d1*p2_x^2*p2_y-4*p1_y^2*d1*p1_x*p2_x-4*p1_x*p2_x*p1_y^2*d2-4*p1_x*p2_x*d2*d1+4*p2_y^3*p1_y*p1_x^2+4*p2_y^3*p1_y*p2_x^2-6*p2_y^2*p1_y^2*p2_x^26*p2_y^2*p1_y^2*p1_x^2+2*p2_y^2*d2*p1_x^2+2*p2_y^2*d2*p2_x^2+2*p2_y^4*p1_x*p2_x12*p1_x^2*p2_x^2*p2_y^2+8*p1_x^3*p2_x*p2_y^2+8*p1_x*p2_x^3*p2_y^2+4*p2_y*p1_x^4*p1_y+2*p2_y^2*p1_x^2*d1+4*p2_y*p1_x^2*p1_y^3+2*p1_y^2*d1*p2_x^2+2*p1_y^2*d1*p1_x^2+4*p2_x^4*p2_y*p1_y +2*p2_x^2*p2_y^2*d1+4*p2_x^2*p2_y*p1_y^3+8*p2_x^3*p1_y^2*p1_x12*p2_x^2*p1_y^2*p1_x^2+2*p2_x^2*p1_y^2*d2+8*p1_x^3*p2_x*p1_y^2+2*p1_x*p2_x*p1_y^4+2*p1_x^2*p1_y^2*d2+12*p1_x^2*p2_x^2*d1-8*p1_x*p2_x^3*d2-8*p1_x*p2_x^3*d1-p2_y^4*p1_x^2p2_y^4*p2_x^2-2*p2_y^2*p1_x^4-2*p2_x^4*p2_y^2-2*p2_x^4*p1_y^2-p2_x^2*p1_y^4-2*p1_x^4*p1_y^2-p1_x^2*p1_y^4+20*p1_x^3*p2_x^3-15*p1_x^2*p2_x^415*p1_x^4*p2_x^2+6*p1_x*p2_x^5+6*p1_x^5*p2_x+2*p2_x^4*d2+2*p2_x^4*d1-p2_x^2*d2^2-p2_x^2*d1^2-p1_x^2*d2^2+2*p1_x^4*d2+2*p1_x^4*d1-p1_x^2*d1^2-p2_x^6p1_x^6+8*p1_x*p2_x*p2_y*p1_y*d1+2*p1_x*p2_x*d2^2-8*p1_x^3*p2_x*d2-8*p1_x^3*p2_x*d1+2*p1_x*p2_x*d1^2+2*p2_x^2*d2*d1+2*p1_x^2*d2*d1)^(1/2))*p1_y+p2_x^2+p2_y^2-d2-p1_x^2-p1_y^2+d1)/(p2_x+p1_x)]; x2 = [ -1/2*(-1/(-8*p2_y*p1_y-8*p1_x*p2_x+4*p2_y^2+4*p1_y^2+4*p2_x^2+4*p1_x^2)*(-4*p2_y^2*p1_y-4*p2_y*d2+4*p2_y^3-8*p1_x*p2_x*p2_y+4*p2_y*p1_x^2-4*p1_y*d1+4*p2_x^2*p2_y+4*p2_x^2*p1_y8*p1_x*p2_x*p1_y-4*p2_y*p1_y^2+4*p1_x^2*p1_y+4*p2_y*d1+4*d2*p1_y+4*p1_y^3-4*(12*p1_x^2*p2_x^2*d2+8*p2_y*d2*p1_x*p2_x*p1_y-8*p2_y^3*p1_y*p1_x*p2_x+12*p2_y^2*p1_y^2*p1_x*p2_x4*p2_y^2*d2*p1_x*p2_x-4*p2_y*d2*p2_x^2*p1_y-4*p2_y*d2*p1_x^2*p1_y-16*p1_x*p2_x^3*p2_y*p1_y+24*p1_x^2*p2_x^2*p2_y*p1_y-16*p1_x^3*p2_x*p2_y*p1_y-4*p1_x*p2_x*p2_y^2*d18*p1_x*p2_x*p2_y*p1_y^3-4*p2_y*p1_x^2*p1_y*d1-4*p1_y*d1*p2_x^2*p2_y-4*p1_y^2*d1*p1_x*p2_x-4*p1_x*p2_x*p1_y^2*d2-4*p1_x*p2_x*d2*d1+4*p2_y^3*p1_y*p1_x^2+4*p2_y^3*p1_y*p2_x^26*p2_y^2*p1_y^2*p2_x^2-6*p2_y^2*p1_y^2*p1_x^2+2*p2_y^2*d2*p1_x^2+2*p2_y^2*d2*p2_x^2+2*p2_y^4*p1_x*p2_x12*p1_x^2*p2_x^2*p2_y^2+8*p1_x^3*p2_x*p2_y^2+8*p1_x*p2_x^3*p2_y^2+4*p2_y*p1_x^4*p1_y+2*p2_y^2*p1_x^2*d1+4*p2_y*p1_x^2*p1_y^3+2*p1_y^2*d1*p2_x^2+2*p1_y^2*d1*p1_x^2+4*p2_x^4*p2_y*p1_y +2*p2_x^2*p2_y^2*d1+4*p2_x^2*p2_y*p1_y^3+8*p2_x^3*p1_y^2*p1_x12*p2_x^2*p1_y^2*p1_x^2+2*p2_x^2*p1_y^2*d2+8*p1_x^3*p2_x*p1_y^2+2*p1_x*p2_x*p1_y^4+2*p1_x^2*p1_y^2*d2+12*p1_x^2*p2_x^2*d1-8*p1_x*p2_x^3*d2-8*p1_x*p2_x^3*d1-p2_y^4*p1_x^2p2_y^4*p2_x^2-2*p2_y^2*p1_x^4-2*p2_x^4*p2_y^2-2*p2_x^4*p1_y^2-p2_x^2*p1_y^4-2*p1_x^4*p1_y^2-p1_x^2*p1_y^4+20*p1_x^3*p2_x^3-15*p1_x^2*p2_x^415*p1_x^4*p2_x^2+6*p1_x*p2_x^5+6*p1_x^5*p2_x+2*p2_x^4*d2+2*p2_x^4*d1-p2_x^2*d2^2-p2_x^2*d1^2-p1_x^2*d2^2+2*p1_x^4*d2+2*p1_x^4*d1-p1_x^2*d1^2-p2_x^6p1_x^6+8*p1_x*p2_x*p2_y*p1_y*d1+2*p1_x*p2_x*d2^2-8*p1_x^3*p2_x*d2-8*p1_x^3*p2_x*d1+2*p1_x*p2_x*d1^2+2*p2_x^2*d2*d1+2*p1_x^2*d2*d1)^(1/2))*p2_y+1/(-8*p2_y*p1_y8*p1_x*p2_x+4*p2_y^2+4*p1_y^2+4*p2_x^2+4*p1_x^2)*(-4*p2_y^2*p1_y-4*p2_y*d2+4*p2_y^3-8*p1_x*p2_x*p2_y+4*p2_y*p1_x^2-4*p1_y*d1+4*p2_x^2*p2_y+4*p2_x^2*p1_y-8*p1_x*p2_x*p1_y4*p2_y*p1_y^2+4*p1_x^2*p1_y+4*p2_y*d1+4*d2*p1_y+4*p1_y^3-4*(12*p1_x^2*p2_x^2*d2+8*p2_y*d2*p1_x*p2_x*p1_y-8*p2_y^3*p1_y*p1_x*p2_x+12*p2_y^2*p1_y^2*p1_x*p2_x-4*p2_y^2*d2*p1_x*p2_x4*p2_y*d2*p2_x^2*p1_y-4*p2_y*d2*p1_x^2*p1_y-16*p1_x*p2_x^3*p2_y*p1_y+24*p1_x^2*p2_x^2*p2_y*p1_y-16*p1_x^3*p2_x*p2_y*p1_y-4*p1_x*p2_x*p2_y^2*d1-8*p1_x*p2_x*p2_y*p1_y^34*p2_y*p1_x^2*p1_y*d1-4*p1_y*d1*p2_x^2*p2_y-4*p1_y^2*d1*p1_x*p2_x-4*p1_x*p2_x*p1_y^2*d2-4*p1_x*p2_x*d2*d1+4*p2_y^3*p1_y*p1_x^2+4*p2_y^3*p1_y*p2_x^2-6*p2_y^2*p1_y^2*p2_x^26*p2_y^2*p1_y^2*p1_x^2+2*p2_y^2*d2*p1_x^2+2*p2_y^2*d2*p2_x^2+2*p2_y^4*p1_x*p2_x12*p1_x^2*p2_x^2*p2_y^2+8*p1_x^3*p2_x*p2_y^2+8*p1_x*p2_x^3*p2_y^2+4*p2_y*p1_x^4*p1_y+2*p2_y^2*p1_x^2*d1+4*p2_y*p1_x^2*p1_y^3+2*p1_y^2*d1*p2_x^2+2*p1_y^2*d1*p1_x^2+4*p2_x^4*p2_y*p1_y +2*p2_x^2*p2_y^2*d1+4*p2_x^2*p2_y*p1_y^3+8*p2_x^3*p1_y^2*p1_x12*p2_x^2*p1_y^2*p1_x^2+2*p2_x^2*p1_y^2*d2+8*p1_x^3*p2_x*p1_y^2+2*p1_x*p2_x*p1_y^4+2*p1_x^2*p1_y^2*d2+12*p1_x^2*p2_x^2*d1-8*p1_x*p2_x^3*d2-8*p1_x*p2_x^3*d1-p2_y^4*p1_x^2p2_y^4*p2_x^2-2*p2_y^2*p1_x^4-2*p2_x^4*p2_y^2-2*p2_x^4*p1_y^2-p2_x^2*p1_y^4-2*p1_x^4*p1_y^2-p1_x^2*p1_y^4+20*p1_x^3*p2_x^3-15*p1_x^2*p2_x^415*p1_x^4*p2_x^2+6*p1_x*p2_x^5+6*p1_x^5*p2_x+2*p2_x^4*d2+2*p2_x^4*d1-p2_x^2*d2^2-p2_x^2*d1^2-p1_x^2*d2^2+2*p1_x^4*d2+2*p1_x^4*d1-p1_x^2*d1^2-p2_x^6p1_x^6+8*p1_x*p2_x*p2_y*p1_y*d1+2*p1_x*p2_x*d2^2-8*p1_x^3*p2_x*d2-8*p1_x^3*p2_x*d1+2*p1_x*p2_x*d1^2+2*p2_x^2*d2*d1+2*p1_x^2*d2*d1)^(1/2))*p1_y+p2_x^2+p2_y^2-d2-p1_x^2-p1_y^2+d1)/(p2_x+p1_x)]; y1 = [ 1/2/(-8*p2_y*p1_y-8*p1_x*p2_x+4*p2_y^2+4*p1_y^2+4*p2_x^2+4*p1_x^2)*(-4*p2_y^2*p1_y-4*p2_y*d2+4*p2_y^3-8*p1_x*p2_x*p2_y+4*p2_y*p1_x^2-4*p1_y*d1+4*p2_x^2*p2_y+4*p2_x^2*p1_y8*p1_x*p2_x*p1_y-4*p2_y*p1_y^2+4*p1_x^2*p1_y+4*p2_y*d1+4*d2*p1_y+4*p1_y^3+4*(12*p1_x^2*p2_x^2*d2+8*p2_y*d2*p1_x*p2_x*p1_y-8*p2_y^3*p1_y*p1_x*p2_x+12*p2_y^2*p1_y^2*p1_x*p2_x4*p2_y^2*d2*p1_x*p2_x-4*p2_y*d2*p2_x^2*p1_y-4*p2_y*d2*p1_x^2*p1_y-16*p1_x*p2_x^3*p2_y*p1_y+24*p1_x^2*p2_x^2*p2_y*p1_y-16*p1_x^3*p2_x*p2_y*p1_y-4*p1_x*p2_x*p2_y^2*d18*p1_x*p2_x*p2_y*p1_y^3-4*p2_y*p1_x^2*p1_y*d1-4*p1_y*d1*p2_x^2*p2_y-4*p1_y^2*d1*p1_x*p2_x-4*p1_x*p2_x*p1_y^2*d2-4*p1_x*p2_x*d2*d1+4*p2_y^3*p1_y*p1_x^2+4*p2_y^3*p1_y*p2_x^26*p2_y^2*p1_y^2*p2_x^2-6*p2_y^2*p1_y^2*p1_x^2+2*p2_y^2*d2*p1_x^2+2*p2_y^2*d2*p2_x^2+2*p2_y^4*p1_x*p2_x12*p1_x^2*p2_x^2*p2_y^2+8*p1_x^3*p2_x*p2_y^2+8*p1_x*p2_x^3*p2_y^2+4*p2_y*p1_x^4*p1_y+2*p2_y^2*p1_x^2*d1+4*p2_y*p1_x^2*p1_y^3+2*p1_y^2*d1*p2_x^2+2*p1_y^2*d1*p1_x^2+4*p2_x^4*p2_y*p1_y +2*p2_x^2*p2_y^2*d1+4*p2_x^2*p2_y*p1_y^3+8*p2_x^3*p1_y^2*p1_x12*p2_x^2*p1_y^2*p1_x^2+2*p2_x^2*p1_y^2*d2+8*p1_x^3*p2_x*p1_y^2+2*p1_x*p2_x*p1_y^4+2*p1_x^2*p1_y^2*d2+12*p1_x^2*p2_x^2*d1-8*p1_x*p2_x^3*d2-8*p1_x*p2_x^3*d1-p2_y^4*p1_x^2p2_y^4*p2_x^2-2*p2_y^2*p1_x^4-2*p2_x^4*p2_y^2-2*p2_x^4*p1_y^2-p2_x^2*p1_y^4-2*p1_x^4*p1_y^2-p1_x^2*p1_y^4+20*p1_x^3*p2_x^3-15*p1_x^2*p2_x^415*p1_x^4*p2_x^2+6*p1_x*p2_x^5+6*p1_x^5*p2_x+2*p2_x^4*d2+2*p2_x^4*d1-p2_x^2*d2^2-p2_x^2*d1^2-p1_x^2*d2^2+2*p1_x^4*d2+2*p1_x^4*d1-p1_x^2*d1^2-p2_x^6p1_x^6+8*p1_x*p2_x*p2_y*p1_y*d1+2*p1_x*p2_x*d2^2-8*p1_x^3*p2_x*d2-8*p1_x^3*p2_x*d1+2*p1_x*p2_x*d1^2+2*p2_x^2*d2*d1+2*p1_x^2*d2*d1)^(1/2))]; y2 = [ 1/2/(-8*p2_y*p1_y-8*p1_x*p2_x+4*p2_y^2+4*p1_y^2+4*p2_x^2+4*p1_x^2)*(-4*p2_y^2*p1_y-4*p2_y*d2+4*p2_y^3-8*p1_x*p2_x*p2_y+4*p2_y*p1_x^2-4*p1_y*d1+4*p2_x^2*p2_y+4*p2_x^2*p1_y8*p1_x*p2_x*p1_y-4*p2_y*p1_y^2+4*p1_x^2*p1_y+4*p2_y*d1+4*d2*p1_y+4*p1_y^3-4*(12*p1_x^2*p2_x^2*d2+8*p2_y*d2*p1_x*p2_x*p1_y-8*p2_y^3*p1_y*p1_x*p2_x+12*p2_y^2*p1_y^2*p1_x*p2_x4*p2_y^2*d2*p1_x*p2_x-4*p2_y*d2*p2_x^2*p1_y-4*p2_y*d2*p1_x^2*p1_y-16*p1_x*p2_x^3*p2_y*p1_y+24*p1_x^2*p2_x^2*p2_y*p1_y-16*p1_x^3*p2_x*p2_y*p1_y-4*p1_x*p2_x*p2_y^2*d18*p1_x*p2_x*p2_y*p1_y^3-4*p2_y*p1_x^2*p1_y*d1-4*p1_y*d1*p2_x^2*p2_y-4*p1_y^2*d1*p1_x*p2_x-4*p1_x*p2_x*p1_y^2*d2-4*p1_x*p2_x*d2*d1+4*p2_y^3*p1_y*p1_x^2+4*p2_y^3*p1_y*p2_x^26*p2_y^2*p1_y^2*p2_x^2-6*p2_y^2*p1_y^2*p1_x^2+2*p2_y^2*d2*p1_x^2+2*p2_y^2*d2*p2_x^2+2*p2_y^4*p1_x*p2_x12*p1_x^2*p2_x^2*p2_y^2+8*p1_x^3*p2_x*p2_y^2+8*p1_x*p2_x^3*p2_y^2+4*p2_y*p1_x^4*p1_y+2*p2_y^2*p1_x^2*d1+4*p2_y*p1_x^2*p1_y^3+2*p1_y^2*d1*p2_x^2+2*p1_y^2*d1*p1_x^2+4*p2_x^4*p2_y*p1_y +2*p2_x^2*p2_y^2*d1+4*p2_x^2*p2_y*p1_y^3+8*p2_x^3*p1_y^2*p1_x12*p2_x^2*p1_y^2*p1_x^2+2*p2_x^2*p1_y^2*d2+8*p1_x^3*p2_x*p1_y^2+2*p1_x*p2_x*p1_y^4+2*p1_x^2*p1_y^2*d2+12*p1_x^2*p2_x^2*d1-8*p1_x*p2_x^3*d2-8*p1_x*p2_x^3*d1-p2_y^4*p1_x^2p2_y^4*p2_x^2-2*p2_y^2*p1_x^4-2*p2_x^4*p2_y^2-2*p2_x^4*p1_y^2-p2_x^2*p1_y^4-2*p1_x^4*p1_y^2-p1_x^2*p1_y^4+20*p1_x^3*p2_x^3-15*p1_x^2*p2_x^415*p1_x^4*p2_x^2+6*p1_x*p2_x^5+6*p1_x^5*p2_x+2*p2_x^4*d2+2*p2_x^4*d1-p2_x^2*d2^2-p2_x^2*d1^2-p1_x^2*d2^2+2*p1_x^4*d2+2*p1_x^4*d1-p1_x^2*d1^2-p2_x^6p1_x^6+8*p1_x*p2_x*p2_y*p1_y*d1+2*p1_x*p2_x*d2^2-8*p1_x^3*p2_x*d2-8*p1_x^3*p2_x*d1+2*p1_x*p2_x*d1^2+2*p2_x^2*d2*d1+2*p1_x^2*d2*d1)^(1/2))]; 屋内PL実験 PL Ref.St. Rover • 訂正データーはPL時間訂正を含み, PL座標は室内座標である。訂正メッ セージの一部であるPL座標とPRNは手動でパネルを通してソフトに導入 する。 • 正確性:コード測位で2∼3mの精度(絶対精度) スムージング測位で10cm程度の精度(相対測位精度) • 演算法は初期座標の情報を必要としない。 Indoor PL test Roof RMS= Submeter PL / Re-transmitter RS RMS=1m Rover RMS=10m • • 屋外アンテナからの衛星信号及び屋内に伝達された同じ衛星信号からのみ計算され る測位位置の比較によって、トランスミッターはマルチパスのレベルのモニタリングを 可能にする。 実験において、下記の原因による大きなマルチパスエラーが生じた。 – 坂道 – 金属製の椅子 アンテナ比較 • Antenna C vs. B • Antenna A vs. B • Antenna C sees only 3 PL • Multipath level Antenna C 定点観測 •5 meter 誤差 •Multipath 誤差 •測位が安定しない •(初期位置は正確) Antenna A 動的観測 • 同じコースを周回 • No フィルター • Pseudorange のみ • No true initialization • 1-2 meter のシフト • Cos方向の誤算 • Carrier phase を用いれば 改善する。(予想) • Carrier phase データは、 乱れていた。 Antenna A 動的観測 (II) • Position chart – 1st yellow – 2nd red – 3rd green • Number – green PL+GPS – GPS NTTコミュニケーションズ・コンソーシアム 社会実験システム構成 システム構成 同期用GPSアンテナ スードライト発信機 スードライトアンテナ スードライト受信機 PL/GPSアンテナ 骨伝導スピーカ 基準局 SS無線機 SS無線機 電子コンパス テンキーボード 電子コンパス モニター背負子 テンキーボード SS無線機 ノートパソコン PL/GPSアンテナ PCカード USBハブ 骨伝導スピーカ バッテリ PL/GPS受信機 シリアルUSB変換 実験ルート 東京海上ビル 新丸ビル 実験ルート 皇居 郵船ビル 70m 丸ビル 東京駅 180m スードライトアンテナ GPSアンテナ 基準局 N 実験機器システムの構成の検討 実験機器の小型化 ¾ノートPCを軽量でかつ高性能のものを用意 ¾アプリケーションの高速化にも取り組み ¾スードライト受信器も小型・軽量化 つくば技術検証実験 地図データの収集・整理 空間データ作成 歩行者同様に、地物も高い位置精度を求めら れる。 道路台帳附図(縮尺1:1500)アナログ図面から ベクトルデータ化(ライン、ポリゴン)、データ構 造化 N 東京海上ビ ル 至皇居 和田倉門交差 点 横断歩道 地片 国 道1 号 新 丸 ビル 歩 道 地片 都道行幸通 り 至東京 駅 郵 船 ビル 地片ネットワ ーク 丸 ビ ル 実験機器の設置および調整 丸ビル 東京海上ビル 郵船ビル 丸ビル 新丸ビル *PL01 1:PRN No 33(1745) 2:DATA Source Internal 3:DATA RATE 50Hz 4:パルスパターン RTCM 5:Clock 内部 6:同期 外部 7:Gprime Mode Off 8:パルス遅延 0chips 9:パルス幅 24chips 10:周波数オフセット 0Hz *PL02 1:PRN No 34(1713) 2:DATA Source Internal 3:DATA RATE 50Hz 4:パルスパターン RTCM 5:Clock 内部 6:同期 外部 7:Gprime Mode Off 8:パルス遅延 24chips 9:パルス幅 24chips 10:周波数オフセット 0Hz 実験機器の設置および調整 新丸ビル 東京海上ビル 郵船ビル 丸ビル 新丸ビル *PL05 1:PRN No 39(0160) 2:DATA Source Internal 3:DATA RATE 50Hz 4:パルスパターン RTCM 5:Clock 内部 6:同期 外部 7:Gprime Mode Off 8:パルス遅延 96chips 9:パルス幅 24chips 10:周波数オフセット 0Hz *PL06 1:PRN No 40(0701) 2:DATA Source Internal 3:DATA RATE 50Hz 4:パルスパターン RTCM 5:Clock 内部 6:同期 外部 7:Gprime Mode Off 8:パルス遅延 120chips 9:パルス幅 24chips 10:周波数オフセット 0Hz 実験機器の設置および調整 IN2400 PL Data Base情報 PRN No. E/D E(m) N(m) U(m) N/A N/A Uoffset 丸の内駐車場出入口(基準局) PL= 33 0 -41.2480 東京海上ビル 郵船ビル 丸ビル 新丸ビル PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= PL= 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 0 -86.3926 0 -142.5537 0 -183.8711 255 0.0000 255 0.0000 0 -12.8862 0 -79.1212 255 110.1406 255 0.0000 255 0.0000 255 -79.0830 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 255 0.0000 -65.5176 -53.3096 -43.2988 -31.5635 0.0000 0.0000 6.2507 24.1347 -13.0195 0.0000 0.0000 24.2793 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 26.4736 26.4482 58.7432 58.7432 0.0000 0.0000 26.3440 26.3280 28.4180 0.0000 0.0000 26.3281 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.4004 -0.4004 0.0000 0.0000 -1.0449 -0.4004 -0.4004 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 モニタテストの模様 今後の測位技術(Pseudolite)研究の方向性 一般ユーザーに対して INDOOR-GPSの守備範囲を超えた部分 =地下街、トンネルその他(DEEP INDOOR) でのポジショニング 数10m 精度 高精度ユーザーに対して OUTDOOR & INDOOR 限られたエリア内の高精度測位を確実に確保 数cm 精度 共通部分(送信機) 送信フォーマット提唱 受信機開発 3年間で具体化までのプロセスが必要。 特化して推進する組織が必要。 周波数 法整備 先端衛星測位研究センター開設 屋上アンテナ GLONASS MAI 工作台 会議テーブル モスクワ航空大学 ルビジウム発振器 東京湾海事VRSネットワークセンター 受付 武漢 UTI 室内スードライト実験室 カルガリー大学 GPS スードライト送信機 基地局受信機 東京湾GISデータセンター ノッティンガム大学 GALILEO 最先端衛星測位 研究センター パーツその他 保管BOX 26000 NU スペクトラムアナライザー GPS基礎評価研究部 GPSシュミレーター プリンタ 海上移動体動態監視センター 武漢大学 国際的な産学協同プロジェクト リアルタイム情報収集サーバー 東京湾海事 海事情報解析サーバー データサーバー カンファレンスルーム WC リフレッシュコーナー 18000 準天頂衛星システムを利用した衛星測位シ ステムの姿 (1) GPS衛星と準天頂衛星を併用 (2008年以降) 現状(GPS衛星のみ) 準天頂衛星 GPS衛星:①∼④ ② ① GPS衛星:①∼④ ③ ④ ② 測位信号 ① × 通信・放送サービス 測位信号 誤差補正情報 測位信号 × 信号遮断 ③ 都市部、山間部 でも有利な仰角 ④ × 測位信号 信号遮断 信号遮断 山陰、ビル陰等で 4機補足ができないと 測位不能となる FM多重放送 DGPS 電子 基準点 測位精度: GPS単独測位 約10m ∼ 誤差補正情報 DGPS補正 約3 ∼ 5m程度 サービスエリア: FMラジオ受信地域 測位誤差 データ 測位精度: サブメータ、デシメータ、センチメータ サービス等級に依存 通信・放送 送信局 サービスエリア: 準天頂衛星: 全国一円及び沿岸部 高仰角からの高精度測位 通信・放送の複合サービスを提供 GPS衛星1機の代替機能を提供