JP2005108246A - センサフュージングを利用した無人移動体の位置推定方法、その装置およびプログラムを記録したコンピュータ読み取り可能な記録媒体。 - Google Patents
センサフュージングを利用した無人移動体の位置推定方法、その装置およびプログラムを記録したコンピュータ読み取り可能な記録媒体。 Download PDFInfo
- Publication number
- JP2005108246A JP2005108246A JP2004288232A JP2004288232A JP2005108246A JP 2005108246 A JP2005108246 A JP 2005108246A JP 2004288232 A JP2004288232 A JP 2004288232A JP 2004288232 A JP2004288232 A JP 2004288232A JP 2005108246 A JP2005108246 A JP 2005108246A
- Authority
- JP
- Japan
- Prior art keywords
- sensor
- signal
- function
- estimation
- channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000009826 distribution Methods 0.000 claims abstract description 9
- 230000008859 change Effects 0.000 claims abstract description 5
- 230000004927 fusion Effects 0.000 claims description 33
- 230000003595 spectral effect Effects 0.000 claims description 27
- 238000012545 processing Methods 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 12
- 238000005314 correlation function Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 9
- 238000012821 model calculation Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims 1
- 238000000354 decomposition reaction Methods 0.000 description 24
- 230000000875 corresponding effect Effects 0.000 description 14
- 230000008569 process Effects 0.000 description 14
- 238000013459 approach Methods 0.000 description 12
- 230000014509 gene expression Effects 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 10
- 230000008901 benefit Effects 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000013215 result calculation Methods 0.000 description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000007257 malfunction Effects 0.000 description 2
- 238000007500 overflow downdraw method Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000001373 regressive effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000005309 stochastic process Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000002730 additional effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/02—Control of position or course in two dimensions
-
- 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/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
Abstract
【課題】 センサフュージングを利用した無人移動体の位置推定方法およびその装置を提供すること。
【解決手段】 無人移動体の位置推定を検知できる少なくとも2つ以上のセンサにおいて所定のパラメータ値を測定する第1段階と、前記パラメータ値を選択的に結合する第2段階と、前記パラメータの予想される範囲内でその変化を検出する第3段階と推定とエラー分布を利用して、センサデータと所定のデータ偏差とにより表現される無人移動体の位置を推定する第4段階とを含むセンサフュージングを利用した無人移動体の位置推定方法を提供する。
【選択図】 図3
【解決手段】 無人移動体の位置推定を検知できる少なくとも2つ以上のセンサにおいて所定のパラメータ値を測定する第1段階と、前記パラメータ値を選択的に結合する第2段階と、前記パラメータの予想される範囲内でその変化を検出する第3段階と推定とエラー分布を利用して、センサデータと所定のデータ偏差とにより表現される無人移動体の位置を推定する第4段階とを含むセンサフュージングを利用した無人移動体の位置推定方法を提供する。
【選択図】 図3
Description
本発明は無人移動体の位置推定方法およびその装置に係り、より詳細にはセンサフュージングを利用した位置推定方法およびその装置に関する。
現状のセンサデータフュージョン技術では、特定の問題に対して的確な解決手法を提供できていない。データフュージョンによる解決手法を定義する場合、研究者は、例えばカルマンフィルタスキームのような公知のカーネルを利用できるものの、専用の手法を構築する必要がある。データを組み合わせるためにこのようなアプローチを用いるシステムは、計算が極端に複雑になることがあり、そのコストも高いものとなり、実現は極めて難しい。
一連の評価者は、現実的な計画としてこれらのいくつかを組み合わせること、および、評価者がリアルタイムで処理しなければならない制約があることを提案してきた。センサデータフュージョンでの主要なアプローチとしては、一般的なカルマンフィルタリング(KF:ordinary Kalman Filtering)、拡張カルマンフィルタリング(EKF:Extended Kalman Filtering)、共分散共通部(CI:Covariance Intersection)、隠れマルコフモデル(HMM:Hidden Markov Models)、部分観測マルコフ決定過程(POMDP:Partially Observable Markov Decision Process)およびベイズネットワークによる手法(Bayesian Networks solution)などがある。
これらの手法の使用には、それぞれ制約および限界がある。主な制約は、分散に依存するモデルを使用する必要があるということであり、EKFの場合には、相互相関の積を計算する必要があり、POMDPの場合には、いくつかのプロセスにおける、以前の状態と現在の状態との間の低い関連を分析する必要がある。それにより、センシング構造を構築する幾つかのよく知られた手法がある。よく知られたセンシング構造は、分散化されたフュージョン構造(a decentralized fusion structure)、分配されたフュージョン構造(a distributed fusion structure)、連合化されたフュージョン構造(a federated fusion structure)、階層化されたフュージョン構造(a hierarchical fusion structure)である。そして、これらのフュージョン構造は、それぞれ幾つかの長所および短所を有している。
分散化され、分配されたフュージョン構造は、拡大縮小可能(scalable)であり、生存可能(survivable)であり、そしてモジュール化可能であるという長所がある。しかしながら、エラー推定がフュージョンチャンネル(a fusion channel)に依存しているという短所を有している。
連合化され、階層化されたフュージョン構造は、フュージョンカスケード(fusion cascade)ごとに回帰的なエラー推定が可能であり、そして、モジュール化できるという長所がある。しかしながら、拡大縮小性が低く(non-scalable)、生存率(survivable)が低いという短所がある。
連合化され、階層化されたフュージョン構造は、フュージョンカスケード(fusion cascade)ごとに回帰的なエラー推定が可能であり、そして、モジュール化できるという長所がある。しかしながら、拡大縮小性が低く(non-scalable)、生存率(survivable)が低いという短所がある。
移動ロボット分野におけるセンサデータフュージョンでは、幾つかの主なアプローチが用いられる。これまでは、EKFが主な状態推定技術であった。EKFは、状態遷移における1次テーラー近似と、観測された状態軌道に関する観測方程式とに基づいている。したがって、EKFの適用は、要求される導関数(derivatives)が存在するという仮定に基づいて、適当な労力により求められる。テーラー線形近似は、多くの場合において、正確な表現(representation)を十分に提供することができない。そして、あまりに粗い近似により、重大な偏り(biases)または収束問題(convergent problems)が多く現れる。
反復、高次フィルタリングおよび統計的な線形近似のような、EKFよりさらに複雑な推定手法を利用可能である。より進んだ手法は、一般的に推定の正確性を向上させるが、具現するためにはあまりにも複雑であり、かつ計算負担が増加してしまう。
したがって、本発明が解決しようとする技術的課題は、拡張可能であり、生存可能であり、そして、モジュール化可能なセンサフュージングを利用した無人移動体の位置推定方法およびその装置を提供することである。
前記した技術的課題を解決するためになされた本発明に係るセンサフュージングを利用した無人移動体の位置推定方法は、無人移動体の位置推定を検知できる少なくとも2つ以上のセンサにおいて所定のパラメータ値を測定する第1段階と、前記パラメータ値を選択的に結合する第2段階と、前記パラメータの予想される範囲内でその変化を検出する第3段階と推定とエラー分布を利用して、センサデータと所定のデータ偏差とにより表現される無人移動体の位置を推定する第4段階とを含むことを特徴としている。
また、前記第1段階は、ソース信号を受信する段階と、前記ソース信号を、高速フーリエ変換(FFT:Fast Fourier Transformation)を使用して周波数領域信号に変換し、スペクトル密度関数を計算する段階と、
多項式を、スペクトルおよび信号に依存した表現で表して、対応する相関関数を計算し、係数を求める段階とを含むことを特徴としている。
多項式を、スペクトルおよび信号に依存した表現で表して、対応する相関関数を計算し、係数を求める段階とを含むことを特徴としている。
前記した技術的課題を解決するためになされた本発明に係るセンサフュージングを利用した無人移動体の位置推定装置は、センサハードウェア階層と対応するソフトウェア階層とから構成され、前記ソフトウェア階層は、電源部のあるセンサを具備し、電源部のあるセンサを提供し、制御信号シーケンス提供を担当し、前記センサからの加工されていないデータを抽出して前処理階層に受け渡すことを担当するセンサチャンネル部と、フュージョンアルゴリズムを通じてプロセシングのための交差および自動相関チャンネルを含む交差プロダクトを計算し、チャンネルパラメータのためのエラーフィードバックを支援し、信号処理表現に対するエラー推定を取得する交差チャンネルモデル計算/フィードバック支援部と、直交加重関数の線形結合を生成し、信号キー特性に対応する推定信号の表現に対する加重関数のセットを生成し、エラー推定方程式を用いて、エラー補償に対する適切な規則を生成する推定分解部と、前記推定分解部において生成された加重関数のセットを、分解された加重係数のセットおよびそれに対応する測定された信号値に対する分散された任意値の推定セットに重畳させ、フュージングされた信号推定の最終結果を取得して、エラー最小化関数を推定する推定重畳部と、最終結果の計算に関する必要な情報を抽出し、前記無人移動体の位置および現状に応じた位置測定に関するキーとなる特徴を抽出し、最終結果を環境状態と関連付け、前記無人移動体の位置に関する、スケーリングおよび調整されていない情報を取得する、最終結果計算部とを含むことを特徴としている。
また、前記センサチャンネル部は、高速フーリエ変換を用いて信号データを処理することによって、スペクトル領域で信号を分析することを特徴としている。
また、前記センサチャンネル部は、スペクトル関数の状態を追跡し、センサチャンネルの状態を予測および分析し、自己回帰法および最小二乗誤差法により、スペクトル関数を多項式に当てはめ、前記センサチャンネルの抽象モデルを用いて前記センサチャンネルのキーパラメータを取得し、環境条件によって所定時間、前記センサチャンネルのチューニングを行うことを特徴としている。
さらに、前記交差チャンネルモデル計算/フィードバック支援部は、相関関数の計算のために、積分畳み込みによる加工されていない信号の変換、スペクトル関数およびパワースペクトル関数を利用し、前記スペクトル関数および前記パワースペクトル関数により、信号チャンネルの交差雑音加重値を決定し、信号スペクトル関数を分析し、初期段階で環境に関する情報を抽出し、交差関連結果、エラー最小化フィードバックのサポートおよびセンサチャンネルのキー周波数を取得することを特徴としている。
そして、本発明は、前記したセンサフュージングを利用した無人移動体の位置推定方法を、コンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体を提供する。
本発明によれば、拡大縮小可能であり、いかなる環境条件であっても容易に拡張または圧縮することができる。そして、生存可能であって、もしセンサソースの1つがなくなるか、誤作動すると、システム全体に被害がおよぶことなく、指数的に関連したエラー推定が減少するだけである。また、モジュール化可能であって、どの種類のセンサが、どの種類のセンシングをつかさどっているかを容易に特定できる。また、フュージョンチャンネルごとにエラー推定および出力補正が可能である。したがって、すべてのセンサソースは、次のレベルのデータフュージョンに対する非回帰的なエラー推定および警告の能力を、それぞれ有している。
以下、本発明の実施の形態を添付した図面を参照して詳細に説明する。
1.導入
本発明は、オブジェクト様式(object-like)の階層化構造アプローチによる、新しいセンサデータフュージョンの技術を提供する。この技術は、カルーネン・レーベ分解方法の多次元拡張により得られた非線形変形の近似化に基づいている。このアプローチの背景となる原理は、通常のフィルタリング技術とは異なっている。このカルーネン・レーベ分解方法のために、補間公式の導関数は不要である。あらかじめ定義された方程式さえも不要である。これはセンサ信号から計算されたスペクトル関数に基づいた自動復帰多項式フィッティング原理によるものである。もちろん、多項式の次数に対する上限を制限する必要はある。本発明の実施形態がテーラー近似化に基づいたフィルタのように複雑であっても、計算に対する負担は極めて小さい。また、推定エラーが分散されるという仮定のもとで、本発明のフュージョンの技術は、エラーの補償により、エラー計算をさらに正確なものとすることができる。本発明によるフュージョンの技術のエントリーポイントに対する深いフィードバックに基づいた最小化によって、他のフィルタリング技術(テーラー近似化を含み)より高精度なエラー推定が可能である。
本発明は、オブジェクト様式(object-like)の階層化構造アプローチによる、新しいセンサデータフュージョンの技術を提供する。この技術は、カルーネン・レーベ分解方法の多次元拡張により得られた非線形変形の近似化に基づいている。このアプローチの背景となる原理は、通常のフィルタリング技術とは異なっている。このカルーネン・レーベ分解方法のために、補間公式の導関数は不要である。あらかじめ定義された方程式さえも不要である。これはセンサ信号から計算されたスペクトル関数に基づいた自動復帰多項式フィッティング原理によるものである。もちろん、多項式の次数に対する上限を制限する必要はある。本発明の実施形態がテーラー近似化に基づいたフィルタのように複雑であっても、計算に対する負担は極めて小さい。また、推定エラーが分散されるという仮定のもとで、本発明のフュージョンの技術は、エラーの補償により、エラー計算をさらに正確なものとすることができる。本発明によるフュージョンの技術のエントリーポイントに対する深いフィードバックに基づいた最小化によって、他のフィルタリング技術(テーラー近似化を含み)より高精度なエラー推定が可能である。
2.一般的なアプローチ
信号処理のための分解方法およびその長所を説明する。一般的な信号処理方法のうち1つは、係数を有する明確な周期関数のセットで信号を表現することである。この方法の利点は、信号を量および質のパラメータで容易に説明できるということである。このようなアプローチにより、周波数領域において信号が検討されることが公知となっている(スペクトル表現)。
信号処理のための分解方法およびその長所を説明する。一般的な信号処理方法のうち1つは、係数を有する明確な周期関数のセットで信号を表現することである。この方法の利点は、信号を量および質のパラメータで容易に説明できるということである。このようなアプローチにより、周波数領域において信号が検討されることが公知となっている(スペクトル表現)。
本発明における信号の周波数領域における表現は、キー周波数とセンサチャンネルの一般的なピクチャとを示す。センサフュージョン技術で最も広く使われる構造を分析してみると、信号をあらかじめ分析する方法や構造は使われないことが分かる。このような技術が産業応用分野で広範囲に適用されているにもかかわらず、移動ロボットへの応用では、ほぼ適用されることはない。ソース品質の分析によって、現在の方法(信号表現方法)の信頼性はよく知られている。ソース品質の分析を通じて、チャンネル状態の予測をモニタリングおよび診断することが可能である。
SLAM(Simultaneous Localization and Mapping)またはセルフナビゲーション技術について、ロボットシステムの認識装置において発生する主要な問題のひとつはセンサ信号の処理、とりわけ、センサデータフュージョンである。しかし、もし信号がなくなるか、またはノイズによって妨害されると、センサデータフュージョンパートに入力される値は妨害され、その結果、妨害された測定結果が、センサデータフュージョンパートから出力され、データ処理の最終段階において誤った位置と方向情報が出力される。このため、ソース信号のモニタリングと診断とを容易に具現する軽くて強固な技術を用いることが必要である。したがって、結合または複合化されたセンサデータフュージョン方法が提案されている。
本発明では、ソース信号の前処理およびデータフュージョンに対応する幾つかの適当な階層がある。明確な理解のために、本発明による方法を詳しく説明する。
次のようなソース信号の前処理に対する一般的な構造が提案される。
(1)ソース信号を受信
(2)高速フーリエ変換(FFT:Fast Fourier Transformation)を使用して前記ソース信号を周波数領域信号に変換
(3)スペクトル密度関数を計算
(4)スペクトルおよび信号に依存する表現として、多項式に当てはめる(処理中の信号の分析−チャンネルの安定性と品質)
(5)対応する相関(共変動)関数および対応する係数を算出
(6)本発明のカーネルに相当する分解方法を実行。
(7)予測とエラー推定モデルとを計算
(1)ソース信号を受信
(2)高速フーリエ変換(FFT:Fast Fourier Transformation)を使用して前記ソース信号を周波数領域信号に変換
(3)スペクトル密度関数を計算
(4)スペクトルおよび信号に依存する表現として、多項式に当てはめる(処理中の信号の分析−チャンネルの安定性と品質)
(5)対応する相関(共変動)関数および対応する係数を算出
(6)本発明のカーネルに相当する分解方法を実行。
(7)予測とエラー推定モデルとを計算
次に、前記した方法の各段階を説明する。第1に、スペクトル関数を示す代表的な多項式を求める。第2に、この多項式を利用してT−R領域における多項式のルート分布を用いて信号の品質の分析が可能である。プロセス(プロセスの代表信号)の状態および条件を表す変換関数を取得することが主な要件である場合に、このようなアプローチは極めて有用である。そして、キー周波数(プロセスの主要で、かつ特徴的な周波数)を取得することができ、ハードウェアのどの部分が信号処理に影響を与えるかを分析することが可能である。
第3に、相関(共変動)関数は、2種の方式で計算される。1)ソース、加工されていない信号を用いる−この方法は、ソースの元の形態を与える。2)スペクトル関数を用いる−この方法は、周波数領域から相関の形態を与える。
次に、本発明の数学的な背景で、幾つかのキーとなる事項を検討する。
次に、本発明の数学的な背景で、幾つかのキーとなる事項を検討する。
3.定義および説明
本発明の要点を定義する。簡略に説明するために、1次元である場合を考える。この方法は容易にn次元に拡張することができる。次元の効果によって独立チャンネルの数が想定される。
本発明の要点を定義する。簡略に説明するために、1次元である場合を考える。この方法は容易にn次元に拡張することができる。次元の効果によって独立チャンネルの数が想定される。
本発明の基本概念は、明確な制約条件下で確率的(任意の)なプロセスのような観測可能なプロセスを表現するアルゴリズムにある。本発明の主な原理は、非周期の確率的なプロセスを、相関関係のない係数を有する一連の直交関数に分解することである。同時に、直交関数に分解する過程の間、エラー最小化方法を使用する。このエラー最小化は、クロスチャンネルおよびインチャンネルのノイズおよびエラーの抑制に非常に効果的である。その結果、前記した方法で得られる結果は、必要な情報の抽出のために容易に使われる。データ分析のための付加的な性質は、前記スペクトル関数を概観するのに用いられる。次に、前記した方法を、順を追って説明する。
3.1 推定の定義
ソース信号の定義は、時間に相関する関数と考えられる。本発明において、環境状態を記述する結果関数を、明確かつ確実に特定することが必要である。そのため、信号システム(SS:Signal System)の統計に基づく推定値
は、演算子F(X(t))パラメータの決定により求められ、この推定値は、物理的に測定された状態座標SS・X(t)Rqを用いてインジケータ
の
により求められる。
ソース信号の定義は、時間に相関する関数と考えられる。本発明において、環境状態を記述する結果関数を、明確かつ確実に特定することが必要である。そのため、信号システム(SS:Signal System)の統計に基づく推定値
は、演算子F(X(t))パラメータの決定により求められ、この推定値は、物理的に測定された状態座標SS・X(t)Rqを用いてインジケータ
の
により求められる。
ここで、1次元の場合(p=q=1)を考え、線形推定の生成および適用の側面から説明する。
ここで、a=a(t)は、ソース信号を分析して決定関数を得ることを意味している。X=X(t)は[0,τ]確率過程の連続的な平均値の2乗であり、ソース信号の偏差のように表現される。bは、任意の値であり、[0,τ]は、SS関数の周期であり、T∈[0,τ]は、測定のために決定された時間を意味している。
任意値のセットY、X(t)、t∈[0,τ]のあらゆる有限次元の分布は、一様分布であり、数式1の線形推定でのaとbパラメータはエラー伝播の最小値
から得られ、数式2の最小値を意味する。
数式2よると、加重関数aは普通t∈[0、T]の関数のクラスに属する。このような関数のクラスは、推定的な決定、または、または任意値セットY,X(t),t∈[0,τ]の過去の分析に基づいて選択される。
L2[0,T]を固定すれば、数式2は数式3のようになる。
数式3からJは、b=b0および
で最小値を有する。ここで、b0=b0(a)は数式4により決定される。
E(estimation)を用いて任意値YとX(t)とをセンタリングした後、
が、得られる。数式6を考える。
数式4を数式1および数式2に代入し、数式5および数式6を考えると、数式7、数式8そして数式9が得られる。
数式7および数式8から、Yの推定およびエラー推定εが次式から得られる。
数式10は、b=b0の場合の、数式1のバイアスのない推定の特性を表している。
aに依存する数式9の関数から、J=J(a)と考えられる。この点から、よく定義されたエラー最小化関数と決定された関数との間の関係が得られる。前記のように、エラー推定システムと一定した定義との間の依存性は、無視されて避けられる。
任意値のセットY、X(t)、t∈[0,τ]のあらゆる有限次元の分布は、一様分布であり、数式1の線形推定でのaとbパラメータはエラー伝播の最小値
から得られ、数式2の最小値を意味する。
L2[0,T]を固定すれば、数式2は数式3のようになる。
で最小値を有する。ここで、b0=b0(a)は数式4により決定される。
aに依存する数式9の関数から、J=J(a)と考えられる。この点から、よく定義されたエラー最小化関数と決定された関数との間の関係が得られる。前記のように、エラー推定システムと一定した定義との間の依存性は、無視されて避けられる。
3.2 分解
相関関数およびクロス相関関数への旧来のアプローチにより、数式11および数式12が得られる。
数式9は、数式13のように書換えられる。
パラメータaを決定するために、加重関数のクラスを考えるにあたって、加重関数のクラスおよびJ=J(a)の効果的な最小化アルゴリズムが表現されなくてはならない。
基本的に、仮想技術領域での他の課題のようにこの問題を解決するために、[0,T]でプライベート相関関数R(t,s)の直交システム[φi:1≦i≦∞]が必要となる。
そして、カルーネン・ローベ直交分解は数式15で表される。
ここでξは実数または複素数であり、数式16のように定義される。
非周期の任意のプロセスは、非相関の任意の係数を有するフーリエ級数で表現できないが、非相関係数を有する直交関数[φi:1≦i≦∞]の級数に拡張される。
数式15は、[0,T]において、2乗平均値に集束する。直交システム[φi:1≦i≦∞]は、L2[0,T]におよぶ。結果的に、それぞれの加重関数a∈L2[0,T]は、ある程度の正確性を有し(L2[0,T]空間メトリックスで)、φi関数の有限セットの線形的な組合せで近似化することが可能である。
システム[φi:1≦i≦∞]はL2[0,T]で直交であるので、数式17のように示すことができる。
ここで、δijはクロネッカーのデルタである。
数式12および数式15ないし数式17から数式18が得られる。
プライベート相関関数R(t,s)の非負の決定からプライベート値(λi≧0)を考える。数式17と前記数式から数式18は数式19のように表現される。
数式17および数式18は、数式15における直交分解の性質を示している。任意のプロセスx=x(t)が集中するために、数式16と数式19から数式20を得る。
前記の式のように、同一のプロセスX=X(t)を考えると、数式16、数式19および数式20から数式15の係数ξiは均一に分布した任意の値から独立し、ξi∈N(0,λi)である。数式15の分解において、正のλiに対応する要素だけが重要であることが明らかである。
相関関数およびクロス相関関数への旧来のアプローチにより、数式11および数式12が得られる。
基本的に、仮想技術領域での他の課題のようにこの問題を解決するために、[0,T]でプライベート相関関数R(t,s)の直交システム[φi:1≦i≦∞]が必要となる。
数式15は、[0,T]において、2乗平均値に集束する。直交システム[φi:1≦i≦∞]は、L2[0,T]におよぶ。結果的に、それぞれの加重関数a∈L2[0,T]は、ある程度の正確性を有し(L2[0,T]空間メトリックスで)、φi関数の有限セットの線形的な組合せで近似化することが可能である。
システム[φi:1≦i≦∞]はL2[0,T]で直交であるので、数式17のように示すことができる。
数式12および数式15ないし数式17から数式18が得られる。
3.3 結合および重畳
この時点で分解の最終結果は、出力関数の推定およびエラー推定−システム品質の決定に対する最小化関数−により示される。mの値を固定させ、数式21のように加重関数を考える。
数式6と数式9に戻って数式22を受け入れ、
数式15、数式17および数式19を基礎として数式23と数式24を得られる。
数式21および数式22から、もしも、iでλi=0であれば、ρi=0であり、したがって、数式24は加重関数aのαiに対して独立している。これがλi(1≦i≦m)を正数である理由である。加重関数の係数
を考えると、数式25が得られる。
数式25は、数式21のクラスでJ=J(a)に対する最小値を提供し、
は数式26のように示される。
数式26を数式24に代入すれば、数式27が得られる。
これから数式21に注目すれば、数式28の結論が得られる。
ここでcor(ξi,y)は任意値ξiとyとの相関係数である。
数式23と数式26によれば、数式24のα0での応答である推定値
は数式29のように与えられる。
伝播エラー
の分散
は数式27か数式28から得られる。
実際には、任意値Yの特性である
と、相関関数r(t)およびR(t,s)とは常にあらかじめ定義されるものではない。
ここで、nは、観察の番号であり、YνおよびXν(t)は、観察v(1≦v≦n)に応答するYとX(t)に対する任意の値である。
この時点で分解の最終結果は、出力関数の推定およびエラー推定−システム品質の決定に対する最小化関数−により示される。mの値を固定させ、数式21のように加重関数を考える。
を考えると、数式25が得られる。
は数式26のように示される。
数式23と数式26によれば、数式24のα0での応答である推定値
は数式29のように与えられる。
の分散
は数式27か数式28から得られる。
実際には、任意値Yの特性である
と、相関関数r(t)およびR(t,s)とは常にあらかじめ定義されるものではない。
3.4 分析
前記した数式は、SSに基づいた移動装置の現在位置と方向の伝播および推定システムに対する数学的なツールである。
これらの結果は、多次元のケースに容易に拡張でき、推定値の多様なパラメータを考慮して、プロセスの幾つかの特性の交差関係および交差ファンクションの解析を可能にする。
前記した数式は、SSに基づいた移動装置の現在位置と方向の伝播および推定システムに対する数学的なツールである。
これらの結果は、多次元のケースに容易に拡張でき、推定値の多様なパラメータを考慮して、プロセスの幾つかの特性の交差関係および交差ファンクションの解析を可能にする。
4.完全な分解アルゴリズム
本発明によるマルチセンサデータフュージョンは、次のように具現される。第1に、プロセスを初期化する。第2に、ソース信号は、デトレンド(detrend)されて集約される。第3に、アナログおよびデジタルごとのアプローチによって、カルーネン・ローベ分解が行われる。第4に、
と
とを計算する。第5には、J(a0)(28)を用いて、エラー推定値を計算する。そして、所定の時間内で、推定値を更新し、エラー推定に対する最小化関数を計算する。そして、前記の第2の段階を繰り返す。
これらの手順は、フュージョン信号処理と予測に対するクローズ−ループシーケンスを提供する。
本発明によるマルチセンサデータフュージョンは、次のように具現される。第1に、プロセスを初期化する。第2に、ソース信号は、デトレンド(detrend)されて集約される。第3に、アナログおよびデジタルごとのアプローチによって、カルーネン・ローベ分解が行われる。第4に、
と
とを計算する。第5には、J(a0)(28)を用いて、エラー推定値を計算する。そして、所定の時間内で、推定値を更新し、エラー推定に対する最小化関数を計算する。そして、前記の第2の段階を繰り返す。
これらの手順は、フュージョン信号処理と予測に対するクローズ−ループシーケンスを提供する。
5.オブジェクト様式のセミレベル情報のフュージョン
前記した内容により、(図1および図2に示すように、)拡大縮小が可能で、生存可能で、かつモジュール化可能な特性を有し、フュージョンチャンネルごとに、出力の修正およびエラー推定が可能なフュージョンシステムを構築できる。拡大縮小が可能な特性によって、フュージョンシステムは、環境条件に応じて、容易に拡張または圧縮される。また、生存可能な特性により、センサソースの1つがなくなるか、または誤作動した場合であっても、システム全体に障害が及ぶことなく、関連したエラー推定が減少するだけである。また、モジュール化が可能な特性により、如何なる種類のセンサが如何なる種類のセンシングに対して関連があるかを容易に識別できる。フュージョンチャンネル当たりエラー推定および出力補正が可能である。これによってあらゆるセンサソースは次のレベルのデータフュージョンに対する自身の非回帰的なエラー推定および警告能力を有する。
前記した内容により、(図1および図2に示すように、)拡大縮小が可能で、生存可能で、かつモジュール化可能な特性を有し、フュージョンチャンネルごとに、出力の修正およびエラー推定が可能なフュージョンシステムを構築できる。拡大縮小が可能な特性によって、フュージョンシステムは、環境条件に応じて、容易に拡張または圧縮される。また、生存可能な特性により、センサソースの1つがなくなるか、または誤作動した場合であっても、システム全体に障害が及ぶことなく、関連したエラー推定が減少するだけである。また、モジュール化が可能な特性により、如何なる種類のセンサが如何なる種類のセンシングに対して関連があるかを容易に識別できる。フュージョンチャンネル当たりエラー推定および出力補正が可能である。これによってあらゆるセンサソースは次のレベルのデータフュージョンに対する自身の非回帰的なエラー推定および警告能力を有する。
センサフュージングを利用した無人移動体の位置推定方法は、次のような段階よりなる。(1)無人移動体の位置の推定を感知するために、複数のモデルパラメータを利用して、センサデータを動的に観測する。(2)パラメータ値を選択的に結合する。(3)所定の予測される範囲内で、前記パラメータ値の変化を検出する。(4)推定およびエラー分布を利用して、センサデータと所定のデータ偏差とにより表現される無人移動体の位置を推定する。
前記した(1)段階で、動的に観測されたパラメータ値は、センサからのリアルタイムの情報である。前記センサによって、同じ動的モデル(モデル成分の分解で表現されるセンサ伝達関数)を構成できる。このセンサ伝達関数を使用して、どのパラメータがより重要であるかを判断できる。このために、あらゆるパラメータ値を計算することが必要である。つまり、各パラメータ加重係数が計算される。加重係数を計算するために自動回帰分析手法を使用できる。数式14ないし数式19はチャンネル(センサ)分解モデルを示す。
前記した(2)段階で、どのパラメータが入力されたかを決定する。センサフュージョン方法を利用した適切な計算のためにセンサパラメータモデル間の関連性を決定する必要がある。このために、交差相関を計算する必要がある。これはモデルの自由度と全体モデルの他のパートとの間の相関度を決定する標準的な手法である。この分析により、どのチャンネル(センサモデル)が、物体の位置およびエラーの推定に効果的であるかを決定できる。
また、エラー分布とそれに対応するチャンネルエラー補償を分析するために、この結果を使用できる。数式9、13、24および28はエラー推定を示し、数式14ないし19は分解モデルを示し、数式21ないし24と27はエラー推定と分解モデル間の連結を示している。
前記した(3)段階では、チャンネル(センサ)のモデルおよびそれに対応するパラメータを、適切なモデル性能のために追跡する必要がある。このために、J(a)をリアルタイムで追跡することが必要である。これは従来技術によれば、膨大なデータ量およびそのフローのために難しい作業である。しかし、本発明では分解モデルを使用するため、リアルタイムで計算可能である。つまり、数式28を利用してリアルタイムでパラメータ値の変化量を追跡できる。
前記した(4)段階で、アルゴリズムに対する最終方程式、すなわち数式28ないし数式33は、チャンネル(センサ)モデル、分解モデルおよび交差相関分析の実行後に生成される。これらの数式は、推定とエラー分散を得るための主な結果である。
図3は、本発明の構成を示すブロック図であって、センサチャンネル部300、交差チャンネルモデル計算/フィードバック支援部320、推定分解部340、推定重畳部360および最終結果計算部380を含んで構成されている。
センサチャンネル部300は、センサハードウェア階層とそれに対応するソフトウェア階層で構成される。ソフトウェア階層は、センサに電源および制御信号シーケンスを供給し、センサからの加工されていないデータを抽出して前処理階層に受け渡すことを担当する。この時点で信号と関連したモデルが構成される。このようなモデルの構成方法を次に説明する。
(1)高速フーリエ変換(FFT)を用いて信号データをプロセシングすることによってスペクトル領域で信号を分析できる。第1に、スペクトル関数形態の状態を追跡し、センサチャンネルの状態を予測または分析することが可能である。第2に、最小平均2乗誤差方法のアプローチによりスペクトル関数で自動回帰方法を用いて多項式をフィットすることが可能である。この方法の長所は、信号関連モデルを用いてインプロセス信号のモニタリングおよび分析を容易にできることである。これにより、診断のような信号チャンネルプロセシングを行える。
(2)信号チャンネルに、チャンネルパラメータブロックのモデルを導入する。チャンネルパラメータに対する柔軟性あるフィードバックを提供するブロックの主なアイディアは、チャンネルパラメータチューニングを支援する。なぜなら、あらゆるディバイスの演算サイクルの間に、インプロセスまたは適切に機能させるオフラインチューニングを行うことが必要であるためである。そのため、もしも、チャンネルの抽象的なモデルを取得することが可能であれば、そのチャンネルのキーパラメータを取得することが可能である。そして、後で環境条件によって所定の時間、そのチャンネルのチューニングを行うことが可能である。
交差チャンネルモデル計算/フィードバック支援部320は、フュージョンアルゴリズムのプロセシングのための交差積を計算する。このために、交差および自動相関チャンネルのような交差関連プロダクトが必要である。チャンネルパラメータのチューニングためのエラーフィードバックの支援として、信号処理方法による理由で、全体信号処理の状況の表現のためのエラー推定を取得するために、幾つかのポイントを規定することが必要である。第1に、相関関数の計算のための2種類の方法がある。1つは積分畳み込みを経て一般的な加工のされていない信号の変換を用いる方法である。そして、もう1つはスペクトル関数およびパワースペクトル関数を用いる方法である。第2に、前記方法は信号チャンネルで簡単な相関関数の計算だけでなく、交差雑音加重値を決定することが可能である。信号のスペクトル関数を分析することによって、初期段階で、環境の重要な特徴に関する情報を抽出できる。これにより、前記交差チャンネルモデル計算/フィードバック支援部320において、交差関連した結果、エラー最小化フィードバック支援、および、センサチャンネルのキー周波数を取得できる。
推定分解部340は、直交加重関数の線形組合せを作りうる。信号キー特性と相応する数学的背景に助けられて推定信号の表現に対した加重関数のセットを作る。エラー推定方程式も考慮せねばならない。そういう方程式でセンサチャンネル部300でエラー補償に対する適切な規則を得る。次のパートは最適信号処理に対する推定計算と最小化方程式とを考慮する。
推定重畳部360は、分解推定された加重関数のセットを、分解加重係数のセットおよび測定された信号値に対する分散任意値の推定の対応するセットに重畳する。この時点で、フュージングされた信号推定の最終結果が得られる。また、エラー最小化関数を推定することも必要である。
前記した最終結果計算部380は、移動ディバイス位置情報の抽出およびエラー関連データ分析し、最終結果の計算と関連した必要な情報を抽出する。移動ディバイスの位置および現状に応じてキー特徴が抽出される。前記方法の次の段階は、環境状態の最終結果を容易に相関させることである。その結果、図4に示したように、移動ディバイスの位置について、スケーリングされず、かつ調整されていない情報を取得する。
一方、この他にセンサからの信号処理に対して次のようなシーケンスを考える。図5に示すように、信号を受信して、システムの加重関数を用いて信号を処理する。(時間移動Ts≦40msであるリアルタイムバッファ)FFTを用いて信号を処理し、信号のスペクトル関数を取得する。スペクトル関数の範囲内で、信号の定性的性質を分析することが可能である。これは加重周波数、スペクトル範囲、形態とタイプおよび全体システムのどのパートが、スペクトルでの周波数部分を定義しているのかを含んでいる。スペクトル関数に対応して、自動回帰モデルを取得して、ルート分布T−R領域によりスペクトルの性質を分析することが可能である。このような分布のタイプおよび種類は、モデル分解階層を表している。このような分析によって、システム全体の幾つかのパラメータ間の連結性および関連性が取得される(例えば、差動ロボットモデルでの速度と方向パラメータ間の関係)。このような自動回帰モデルを取得した後、リアルタイムのモニタリングと各センシングチャンネルが診断することによって簡潔で、かつ柔軟な数学的ソフトウェアセットを構成することができる。
本発明は無人移動物体、例えばロボットの位置推定およびエラー補正に有効に適用することができる。
300 センサチャンネル部
320 交差チャンネルモデル計算/フィードバック支援部
340 推定分解部
360 推定重畳部
380 最終結果計算部
320 交差チャンネルモデル計算/フィードバック支援部
340 推定分解部
360 推定重畳部
380 最終結果計算部
Claims (7)
- (a)少なくとも2つ以上のセンサから無人移動体の位置推定を検知できる所定のパラメータ値を測定する段階と、
(b)前記パラメータ値を選択的に結合する段階と、
(c)所定の予想される範囲内で、前記パラメータ値の変化を検出する段階と、
(d)推定とエラー分布を利用して、センサデータと所定のデータ偏差とにより表現される無人移動体の位置を推定する段階と、
を含むことを特徴とするセンサフュージングを利用した無人移動体の位置推定方法。 - 前記段階(a)は、
ソース信号を受信する段階と、
前記ソース信号を、高速フーリエ変換(FFT:Fast Fourier Transformation)を使用して周波数領域信号に変換し、スペクトル密度関数を計算する段階と、
多項式を、スペクトルおよび信号に依存した表現で表して、対応する相関関数を計算し、係数を求める段階と、
を含むことを特徴とする請求項1に記載のセンサフュージングを利用した無人移動体の位置推定方法。 - 複数のセンサを具備し、前記センサからの加工されていないデータを抽出して前処理階層に受け渡すセンサチャンネル部と、
所定のフュージョンアルゴリズムを通じて、プロセシングのための交差および自動相関チャンネルを含む交差プロダクトを計算し、チャンネルパラメータのためのエラーフィードバックを支援し、信号処理表現に対するエラー推定を取得する交差チャンネルモデル計算/フィードバック支援部と、
直交加重関数の線形結合を生成し、信号キー特性に対応する推定信号表現に対する加重関数のセットを生成し、エラー推定方程式を用いて、エラー補償に対する規則を生成する推定分解部と、
前記推定分解部において生成された加重関数を、分解された加重係数のセットおよびそれに対応する測定された信号値に対する分散された任意値の推定セットに重畳させ、フュージングされた信号推定の最終結果を取得し、エラー最小化関数を推定する推定重畳部と、
最終結果の計算に関する必要な情報を抽出し、前記無人移動体の位置および現状に応じた位置測定に関するキーとなる特徴を抽出し、最終結果を環境状態と関連付け、前記無人移動体の位置に関する、スケーリングおよび調整されていない情報を取得する、最終結果計算部と、
を含むことを特徴とするセンサフュージングを利用した無人移動体の位置推定装置。 - 前記センサチャンネル部は、
高速フーリエ変換を用いて信号データを処理することによって、スペクトル領域で信号を分析すること、
を特徴とする請求項3に記載のセンサフュージングを利用した無人移動体の位置推定装置。 - 前記センサチャンネル部は、
スペクトル関数の状態を追跡し、センサチャンネルの状態を予測および分析し、自己回帰法および最小二乗誤差法により、スペクトル関数に多項式を当てはめ、前記センサチャンネルの抽象モデルを用いて前記センサチャンネルのキーパラメータを取得し、環境条件によって所定時間に前記センサチャンネルのチューニングを行うこと、
を特徴とする請求項4に記載のセンサフュージングを利用した無人移動体の位置推定装置。 - 前記交差チャンネルモデル計算/フィードバック支援部は、
相関関数計算のために、積分畳込みによる加工されていない信号の変換、スペクトル関数およびパワースペクトル関数を利用し、前記スペクトル関数および前記パワースペクトル関数により、信号チャンネルの交差雑音加重値を決定し、信号スペクトル関数を分析し、初期段階で環境に関する情報を抽出し、交差関連結果、エラー最小化フィードバックのサポートおよびセンサチャンネルのキー周波数を取得すること、
を特徴とする請求項3に記載のセンサフュージングを利用した無人移動体の位置推定装置。 - 請求項1または請求項2に記載のセンサフュージングを利用した無人移動体の位置推定方法を、コンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020030068073A KR100707168B1 (ko) | 2003-09-30 | 2003-09-30 | 센서 퓨징을 이용한 무인 이동체의 위치 추정 방법 및 장치 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005108246A true JP2005108246A (ja) | 2005-04-21 |
Family
ID=34545534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004288232A Pending JP2005108246A (ja) | 2003-09-30 | 2004-09-30 | センサフュージングを利用した無人移動体の位置推定方法、その装置およびプログラムを記録したコンピュータ読み取り可能な記録媒体。 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20070233336A1 (ja) |
JP (1) | JP2005108246A (ja) |
KR (1) | KR100707168B1 (ja) |
CN (1) | CN1645283A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012248032A (ja) * | 2011-05-27 | 2012-12-13 | Fujitsu Ltd | 地図処理方法及びプログラム、並びにロボットシステム |
JP2012531890A (ja) * | 2009-06-30 | 2012-12-13 | アストン ユニバーシティ | 生物細胞の特性又は挙動を特徴づけるためにバイオセンシングするマイクロプレートの振動 |
JP2016110576A (ja) * | 2014-12-10 | 2016-06-20 | 株式会社豊田中央研究所 | 自己位置推定装置及び自己位置推定装置を備えた移動体 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2004294651A1 (en) * | 2003-10-21 | 2005-06-16 | Proxy Aviation Systems, Inc. | Methods and apparatus for unmanned vehicle control |
US7395156B2 (en) * | 2005-06-23 | 2008-07-01 | Raytheon Company | System and method for geo-registration with global positioning and inertial navigation |
US7840287B2 (en) * | 2006-04-13 | 2010-11-23 | Fisher-Rosemount Systems, Inc. | Robust process model identification in model based control techniques |
DE102007029299B4 (de) * | 2007-06-22 | 2011-12-22 | Fraba Ag | Optischer Sensor für Positionieraufgaben |
US20090240390A1 (en) * | 2008-03-21 | 2009-09-24 | Nenad Nenadic | System and method for component monitoring |
CN103162678B (zh) * | 2013-03-11 | 2015-05-13 | 哈尔滨工业大学 | 批量mems陀螺信息融合方法 |
US9037403B2 (en) * | 2013-03-26 | 2015-05-19 | Toyota Motor Engineering & Manufacturing North America, Inc. | Intensity map-based localization with adaptive thresholding |
US9755869B2 (en) | 2013-11-18 | 2017-09-05 | Bae Systems Information And Electronic Systems Integration Inc. | Process for tunnelized cyclostationary to achieve low-energy spectrum sensing |
WO2016123032A1 (en) * | 2015-01-26 | 2016-08-04 | Batten George W Jr | Floor patterns for navigation corrections |
US9802620B2 (en) * | 2015-12-18 | 2017-10-31 | GM Global Technology Operations LLC | Position error estimate and implementation for autonomous driving |
US10666473B2 (en) | 2016-09-27 | 2020-05-26 | Bae Systems Information And Electronic Systems Integration Inc. | Techniques for implementing a portable spectrum analyzer |
KR102463176B1 (ko) | 2017-10-16 | 2022-11-04 | 삼성전자주식회사 | 위치 추정 장치 및 방법 |
CN111324172B (zh) * | 2018-12-13 | 2021-11-23 | 北京小米松果电子有限公司 | 遥杆的校准方法、装置、电子设备和存储介质 |
CN112598019B (zh) * | 2020-10-28 | 2024-04-02 | 南京航空航天大学 | 异类振动信号时空多互相关分析和自适应加权融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001083232A (ja) * | 1999-08-16 | 2001-03-30 | Raytheon Co | 受動ターゲット位置決定装置 |
-
2003
- 2003-09-30 KR KR1020030068073A patent/KR100707168B1/ko not_active IP Right Cessation
-
2004
- 2004-09-24 CN CNA2004100995016A patent/CN1645283A/zh active Pending
- 2004-09-29 US US10/951,908 patent/US20070233336A1/en not_active Abandoned
- 2004-09-30 JP JP2004288232A patent/JP2005108246A/ja active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012531890A (ja) * | 2009-06-30 | 2012-12-13 | アストン ユニバーシティ | 生物細胞の特性又は挙動を特徴づけるためにバイオセンシングするマイクロプレートの振動 |
JP2012248032A (ja) * | 2011-05-27 | 2012-12-13 | Fujitsu Ltd | 地図処理方法及びプログラム、並びにロボットシステム |
JP2016110576A (ja) * | 2014-12-10 | 2016-06-20 | 株式会社豊田中央研究所 | 自己位置推定装置及び自己位置推定装置を備えた移動体 |
Also Published As
Publication number | Publication date |
---|---|
KR20050031809A (ko) | 2005-04-06 |
US20070233336A1 (en) | 2007-10-04 |
CN1645283A (zh) | 2005-07-27 |
KR100707168B1 (ko) | 2007-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2005108246A (ja) | センサフュージングを利用した無人移動体の位置推定方法、その装置およびプログラムを記録したコンピュータ読み取り可能な記録媒体。 | |
Santos et al. | An evaluation of 2D SLAM techniques available in robot operating system | |
JP6068468B2 (ja) | 予測及び予知のための逐次的カーネル回帰モデリング方法 | |
Puig et al. | Passive robust fault detection using interval observers: Application to the DAMADICS benchmark problem | |
Xiong et al. | Detection of satellite attitude sensor faults using the UKF | |
JP2020017293A (ja) | 検出データからマシンの挙動のための動的モデルを形成する方法 | |
Baghernezhad et al. | Computationally intelligent strategies for robust fault detection, isolation, and identification of mobile robots | |
Gning et al. | Mixture of uniform probability density functions for non linear state estimation using interval analysis | |
Morales-Menendez et al. | Real-time monitoring of complex industrial processes with particle filters | |
JP2014524094A (ja) | パターン・シーケンスを持つカーネル回帰モデリングを用いる監視システム | |
Huang et al. | Adaptive recursive decentralized cooperative localization for multirobot systems with time-varying measurement accuracy | |
US7787696B2 (en) | Systems and methods for adaptive sampling and estimating a systematic relationship between a plurality of points | |
Zhu et al. | Robust particle filter for state estimation using measurements with different types of gross errors | |
Zhang et al. | Neural network based uncertainty prediction for autonomous vehicle application | |
Heravi et al. | Long term trajectory prediction of moving objects using gaussian process | |
Hébert et al. | Probabilistic map learning: Necessity and difficulties | |
Seignez et al. | Experimental vehicle localization by bounded-error state estimation using interval analysis | |
CN111623796B (zh) | 一种基于信息融合的钢轨里程估计方法 | |
Martinez et al. | Feature uncertainty estimation in sensor fusion applied to autonomous vehicle location | |
Angelov et al. | Evolving inferential sensors in the chemical process industry | |
Lenac et al. | Revival of filtering based SLAM? Exactly sparse delayed state filter on Lie groups | |
Jatana | Predictive maintenance of industrial robots | |
Huang et al. | Centralized fusion of unscented kalman filter based on huber robust method for nonlinear moving target tracking | |
CN111679651A (zh) | 用于故障引起的变结构变参数系统的辨识方法及系统 | |
Lin et al. | Visual tracking with sensing dynamics compensation using the Expectation-Maximization algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051219 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080212 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080701 |