(57)【要約】  本発明は、ＩＭＵに基づくｕｗｂ測位異常値処理方法を開示し、前記ＩＭＵに基づくｕｗｂ測位異常値処理方法は、ドローンのＩＭＵデータを取得するステップＳ１と、ドローンのオイラー角を計算するステップＳ２と、Ｘ、Ｙ、Ｚ軸の加速度及び△ｔ時間内のＸ、Ｙ、Ｚ軸上の変位を計算し、ドローンの変位値Ｑを推定するステップＳ３と、ｕｗｂ測定値について異常検出を行うステップＳ４と、異常なｕｗｂ測定値に対してデータ融合を行い、ｕｗｂ測定値を補正するステップＳ５と、を含む。本発明は、従来の生データの前処理過程における異常値の判別が不正確であること、及び異常値補正の精度などの問題を解決し、測定ハードウェアデータの正確性を向上させる。さらに、本発明はハードウェア機器を追加する必要がなく、複雑な環境でデータを測定するときのｕｗｂセンサの安定性を高め、また、後続のｕｗｂ測位解に有効で正確な測定出力データを提供し、ｕｗｂ測位の安定性、効率及びリアルタイム性を高める。【選択図】図１
2022500617.tif 000002
【特許請求の範囲】
【請求項１】
ＩＭＵに基づくｕｗｂ測位異常値処理方法であって、
ドローンのＩＭＵデータを取得するステップＳ１と、
加速度計とジャイロスコープを組み合わせて、ドローンのオイラー角を計算するステップＳ２と、
ステップＳ２で得られたドローンのオイラー角から、Ｘ、Ｙ、Ｚ軸の加速度及び△ｔ時間内のＸ、Ｙ、Ｚ軸上の変位を計算し、ドローンの変位値Ｑを推定するステップＳ３と、
ステップＳ３で得られた変位値Ｑ及びｕｗｂセンサの距離データから、ｕｗｂ測定値について異常検出を行い、ｕｗｂ測定値が異常である場合、ステップＳ５に入るステップＳ４と、
異常なｕｗｂ測定値に対してデータ融合を行い、ｕｗｂ測定値を補正するステップＳ５と、を含む、ことを特徴とするＩＭＵに基づくｕｗｂ測位異常値処理方法。
【請求項２】
前記ステップＳ２では、具体的には、
オイラー角及びジャイロスコープのドリフトを推定し、つまり、
【数４８】
2022500617.tif 000050
であり、ｂｇ，ｋはジャイロスコープのドリフト量、θ、φはそれぞれピッチ角、ロール角である。）のようにし、
【数４９】
2022500617.tif 000051
のため（ωｂはドローンの角速度、Ｒは測地座標系に対するボディ座標系の回転行列である。）、第３列を抽出し、つまり、
【数５０】
2022500617.tif 000052
であり、
【数５１】
2022500617.tif 000053
ここで、ドリフト誤差ｂｇを減算し、次にジャイロスコープ測定ホワイトノイズｗｇを減算し、ドリフトノイズモデリングにはホワイトガウスノイズとすることができ、
【数５２】
2022500617.tif 000054
はジャイロスコープのノイズであり、離散化して、
【数５３】
2022500617.tif 000055
（式中、Ｔｓはサンプリング時間である。）を得て、
単純化して状態変数を分離し、
【数５４】
2022500617.tif 000056
を得て、
単純化には、Ｘｋ＝Φｋ−１Ｘｋ−１＋Γｋ−１ｗｋ−１とし、
ここで、ｋ−１時刻でのシステム遷移行列は、
【数５５】
2022500617.tif 000057
であり、ｋ−１時刻での状態ベクトルは、
【数５６】
2022500617.tif 000058
であり、観測モデルは、
【数５７】
2022500617.tif 000059
であり、ここで、
【数５８】
2022500617.tif 000060
であり、Ｋ時刻での観測行列は、
【数５９】
2022500617.tif 000061
であり、前記のように、単純化したプロセスモデルと観測モデルは、
【数６０】
2022500617.tif 000062
になり、プロセスモデルと観測モデルを得た後、初期化
【数６１】
2022500617.tif 000063
（ここで、Ｅ（Ｘ０）は零時刻での数学的期待値、Ｅ［（Ｘ０−Ｅ（Ｘ０））（Ｘ０−Ｅ（Ｘ０））Ｔ］は零時刻の誤差共分散の計算を示す。）を行い、
状態推定予測Ｘｋ｜ｋ−１＝Φｋ−１Ｘｋ−１を行い、
誤差共分散予測：
Ｐｋ｜ｋ−１＝Φｋ−１Ｐｋ−１｜ｋ−１Φｋ−１Ｔ＋Γｋ−１Ｑｋ−１Γｋ−１Ｔ（Ｑｋ−１はシステムノイズ共分散行列を示す。）を行い、
カルマンゲイン行列は、
【数６２】
2022500617.tif 000064
（ここで、Ｈｋ＝［Ｏ３×３−ｇＩ］。ＲＴは測定ノイズ共分散行列である。）であり、
状態推定更新：
Ｘｋ｜ｋ＝Ｘｋ｜ｋ−１＋Ｋｋ（Ｚ−ＨｋＸｋ｜ｋ−１）
（ここで、Ｚは観測データである。）を行い、
誤差共分散更新：
Ｐｋ｜ｋ＝（Ｉ−ＫｋＨｋ）Ｐｋ｜ｋ−１
（Ｉは単位行列である。）を行い、状態ベクトルＸｋから正確なオイラー角を得る、ことを特徴とする請求項１に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法。
【請求項３】
前記ステップＳ３では、Ｘ、Ｙ、Ｚ軸の加速度の計算は、具体的には、
【数６３】
2022500617.tif 000065
（ここで、ｘｋ（ｉ，ｊ）はｘｋの第ｉ行、第ｊ列のデータを表す。）から、ピッチ角θとロール角φを求め、ドローンの飛行速度が非常に小さく、すなわちその瞬時加速度を無視できるため、測定した加速度がｇ重力加速度に近く、測定した加速度の値ａｂから、
【数６４】
2022500617.tif 000066
ことを特徴とする請求項２に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法。
【請求項４】
前記ステップＳ３では、△ｔ時間内のＸ、Ｙ、Ｚ軸上の変位の計算は、式
【数６５】
2022500617.tif 000067
であり、ここで、ｖｔは現在速度であり、ｖｔ＝（Ｐｔ−Ｐｔ−１）／Δｔは三次元ベクトルであり、Ｐｔはｔ時刻のｕｗｂ測位システムの位置ソリューションによる位置であり、前時刻の平均速度を現在時間の速度とし、また、
【数６６】
2022500617.tif 000068
はΔｔ時間内に零に近いので、
Δｓ＝ｖｔΔｔ＋ａｔΔｔであり、
式中、ａはＩＭＵソリューションにより得られる正確な加速度のナビゲーション座標系での加速度であり、ｔはシステムの動作時間であり、上記式よりｔ＋Δｔ時刻の推定変位ｓｔ＋Δｓを得て、
ここで、
【数６７】
2022500617.tif 000069
であり、Δｓ１、Δｓ２、Δｓ３はそれぞれサンプリング時間Δｔ内のナビゲーション座標系でのＸ、Ｙ、Ｚ軸の変位である、ことを特徴とする請求項３に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法。
【請求項５】
前記ステップＳ３では、ドローンの変位値Ｑの計算式は、
【数６８】
2022500617.tif 000070
であり、
ｘ、ｙ、ｚは前回測定した測位値であり、（ａ＊１、ａ＊２、ａ＊３）は＊番目の基地局の位置値である、ことを特徴とする請求項４に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法。
【請求項６】
前記ステップＳ４では、異常値を検出して判断する際に、式は
【数６９】
2022500617.tif 000071
はそれぞれｕｗｂの現在の測定距離値と前の時刻の測定値であり、ＱはΔｔ時間内の変位値であり、ｍｅａｎはｎ個の測定データの標準偏差値であり、ｋは比重パラメータであり、具体的な環境パラメータに応じて調整することができ、式
【数７０】
2022500617.tif 000072
を測定異常値とする、ことを特徴とする請求項５に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法。
【請求項７】
前記ステップＳ５ｕｗｂでは、測定値を補正する式は、
【数７１】
2022500617.tif 000073
であり、
ここで、ｓｔｄはｎ個の測定データの平均値である、ことを特徴とする請求項６に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
発明は、ｕｗｂ測位の技術分野に関し、特にＩＭＵに基づくｕｗｂ測位異常値処理方法に関する。
【背景技術】
【０００２】
Ｕｗｂ測位技術は、１ＧＨｚ以上の周波数帯域幅を用いた無線キャリア通信技術であり、正弦波キャリアではなく、ナノ秒スケールの非正弦波ナローパルスを用いてデータを伝送するため、占めるスペクトル範囲が広く、無線通信を使用するにもかかわらず、そのデータ伝送速度が数百メガビット／ｓ以上に達する。ｕｗｂは高速データ伝送、低消費電力、高安全性、高測位精度の特徴があるため、さまざまな分野の室内測位に広く使用されており、本特許のｕｗｂの応用背景は、室内でドローンに高精度の測位データを提供することである。室内環境の複雑さ及び各ｕｗｂハードウェアの違いのため、ｕｗｂ距離測定入力値の生データにはジャンプ、ロスが生じ、実際のデータが損失されてしまい、この場合、ｕｗｂ測位に連続的かつ正確な測定データ値を提供するために、このような異常データ値を処理・充填する必要がある。
【０００３】
現在、異常値の常の処理方式としては、異常値を含むデータを削除し、異常値を欠損値として、欠損値処理方法で処理し、平均値で補正するか、処理しないなどがあるが、以上のような異常値処理方法では、ｕｗｂ測位データの処理過程に一定の欠陥があり、たとえば、平均値で補正する方法では、ドローンが大きな加速度で飛行する場合、平均値で異常値を代えると、明らかなデータの歪みが発生し、一方、ほかの処理方法のいずれも、後続のアルゴリズムによる推定に誤差をもたらす。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明の目的は、従来技術の欠陥を解決し、ＩＭＵに基づくｕｗｂ測位異常値処理方法を提供することであり、４つのｕｗｂ基地局でｕｗｂ測位を行う例では、従来の生データ前処理過程における異常値の判別が不正確であること、及び異常値補正の精度などの問題を解決し、測定ハードウェアデータの正確性を向上させる。本発明は、ハードウェア機器を追加する必要がなく、複雑な環境でデータを測定するときのｕｗｂセンサの安定性を高め、また、後続のｕｗｂ測位解に有効で正確な測定出力データを提供し、ｕｗｂ測位の安定性、効率及びリアルタイム性を高める。
【課題を解決するための手段】
【０００５】
上記目的を達成させるために、本発明による技術案は以下のとおりである。
ＩＭＵに基づくｕｗｂ測位異常値処理方法であって、
ドローンのＩＭＵデータを取得するステップＳ１と、
加速度計とジャイロスコープを組み合わせて、ドローンのオイラー角を計算するステップＳ２と、
ステップＳ２で得られたドローンのオイラー角から、Ｘ、Ｙ、Ｚ軸の加速度及び△ｔ時間内のＸ、Ｙ、Ｚ軸上の変位を計算し、ドローンの変位値Ｑを推定するステップＳ３と、
ステップＳ３で得られた変位値Ｑ及びｕｗｂセンサの距離データから、ｕｗｂ測定値について異常検出を行い、ｕｗｂ測定値が異常である場合、ステップＳ５に入るステップＳ４と、
異常なｕｗｂ測定値に対してデータ融合を行い、ｕｗｂ測定値を補正するステップＳ５と、を含む。
【０００６】
さらに、前記ステップＳ２では、具体的には、
オイラー角及びジャイロスコープのドリフトを推定し、つまり、（式
【数１】
2022500617.tif 000003
であり、ｂｇ，ｋはジャイロスコープのドリフト量、θ、φはそれぞれピッチ角、ロール角である。）のようにし、
【数２】
2022500617.tif 000004
であり、
【数３】
2022500617.tif 000005
離散化して、
【数４】
2022500617.tif 000006
を得て、単純化には、Ｘｋ＝Φｋ−１Ｘｋ−１＋Γｋ−１ｗｋ−１とし、
ここで、ｋ−１時刻でのシステム遷移行列は、
【数５】
2022500617.tif 000007
であり、
ｋ−１時刻でのシステムノイズ行列は、
【数６】
2022500617.tif 000008
であり、
ｋ−１時刻での状態ベクトルは、
【数７】
2022500617.tif 000009
であり、
ｋ−１時刻でのシステムノイズは、
【数８】
2022500617.tif 000010
であり、
観測モデルは、
【数９】
2022500617.tif 000011
であり、
ここで、
【数１０】
2022500617.tif 000012
であり、
Ｋ時刻での観測行列は、
【数１１】
2022500617.tif 000013
であり、
前記のように、単純化したプロセスモデルと観測モデルは、
【数１２】
2022500617.tif 000014
になり、
プロセスモデルと観測モデルを得た後、初期化
【数１３】
2022500617.tif 000015
（ここで、Ｅ（Ｘ０）は零時刻での数学的期待値、Ｅ［（Ｘ０−Ｅ（Ｘ０））（Ｘ０−Ｅ（Ｘ０））Ｔ］は零時刻の誤差共分散の計算を示す。）を行い、
状態推定予測Ｘｋ｜ｋ−１＝Φｋ−１Ｘｋ−１を行い、
誤差共分散予測：
Ｐｋ｜ｋ−１＝Φｋ−１Ｐｋ−１｜ｋ−１Φｋ−１Ｔ＋Γｋ−１Ｑｋ−１Γｋ−１Ｔ（Ｑｋ−１はシステムノイズ共分散行列を示す。）を行い、
カルマンゲイン行列は、
【数１４】
2022500617.tif 000016
（ここで、Ｈｋ＝［Ｏ３×３−ｇＩ］。ＲＴは測定ノイズ共分散行列である。）であり、
状態推定更新：
Ｘｋ｜ｋ＝Ｘｋ｜ｋ−１＋Ｋｋ（Ｚ−ＨｋＸｋ｜ｋ−１）
（ここで、Ｚは観測データである。）を行い、
誤差共分散更新：
Ｐｋ｜ｋ＝（Ｉ−ＫｋＨｋ）Ｐｋ｜ｋ−１
（Ｉは単位行列である。）を行い、状態ベクトルＸｋから正確なオイラー角を得る。
さらに、前記ステップＳ３では、Ｘ、Ｙ、Ｚ軸の加速度の計算は、具体的には、
【数１５】
2022500617.tif 000017
（ここで、ｘｋ（ｉ，ｊ）はｘｋの第ｉ行、第ｊ列のデータを表す。）から、ピッチ角θとロール角φを求め、ドローンの飛行速度が非常に小さく、すなわちその瞬時加速度を無視できるため、測定した加速度がｇ重力加速度に近く、測定した加速度の値ａｂから、
【数１６】
2022500617.tif 000018
【０００７】
さらに、前記ステップＳ３では、△ｔ時間内のＸ、Ｙ、Ｚ軸上の変位の計算は、式
【数１７】
2022500617.tif 000019
であり、
ここで、ｖｔは現在速度であり、ｖｔ＝（Ｐｔ−Ｐｔ−１）／Δｔは三次元ベクトルであり、Ｐｔはｔ時刻のｕｗｂ測位システムの位置ソリューションによる位置であり、前時刻の平均速度を現在時間の速度とし、また、
【数１８】
2022500617.tif 000020
がΔｔ時間内に零に近いので、
Δｓ＝ｖｔΔｔ＋ａｔΔｔであり、であり、
式中、ａはＩＭＵソリューションにより得られる正確な加速度のナビゲーション座標系での加速度であり、ｔはシステムの動作時間であり、上記式よりｔ＋Δｔ時刻の推定変位ｓｔ＋Δｓを得て、
ここで、
【数１９】
2022500617.tif 000021
であり、Δｓ１、Δｓ２、Δｓ３はそれぞれサンプリング時間Δｔ内のナビゲーション座標系でのＸ、Ｙ、Ｚ軸の変位である。
【０００８】
さらに、前記ステップＳ３では、ドローンの変位値Ｑの計算式は、
【数２０】
2022500617.tif 000022
であり、
ｘ、ｙ、ｚは前回測定した測位値であり、（ａ＊１、ａ＊２、ａ＊３）は＊番目の基地局の位置値である。
【０００９】
さらに、前記ステップＳ４では、異常値を検出して判断する際に、式は
【数２１】
2022500617.tif 000023
はそれぞれｕｗｂの現在の測定距離値と前の時刻の測定値であり、ＱはΔｔ時間内の変位値であり、ｍｅａｎはｎ個の測定データの標準偏差値であり、ｋは比重パラメータであり、具体的な環境パラメータに応じて調整することができ、式
【数２２】
2022500617.tif 000024
を測定異常値とする。
【００１０】
さらに、前記ステップＳ５ｕｗｂでは、測定値を補正する式は、
【数２３】
2022500617.tif 000025
であり、
ここで、ｓｔｄはｎ個の測定データの平均値である。
【発明の効果】
【００１１】
従来技術に比べて、本技術案の原理及び利点は以下のとおりである。
１．本技術案では、ドローンに備えるＩＭＵセンサと組み合わせてドローンのΔｔ時間内の変位を予測することにより、ｕｗｂで測定された距離値の変動範囲を効果的に予測し、データが異常であるか否かを正確に判別できる。
２．本技術案では、異常値データが判断され、且つドローンのΔｔ時間内の変位が推定されることにより、異常値をより正確に補正し、生データ値を効果的に回復できる。
【００１２】
本発明の実施例又は従来技術の技術案をより明瞭に説明するために、以下、実施例又は従来技術の説明に必要なサービスを簡単に説明するが、勿論、以下の説明における図面は、本発明のいくつかの実施例に過ぎず、当業者であれば、創造的な努力を必要とせずに、これら図面に基づいてほかの図面を得ることができる。
【図面の簡単な説明】
【００１３】
【図１】本発明のＩＭＵに基づくｕｗｂ測位異常値処理方法の原理のフローチャートである。
【発明を実施するための形態】
【００１４】
以下、特定の実施例にて本発明をさらに説明する。
【００１５】
図１に示すように、本実施例に記載のＩＭＵに基づくｕｗｂ測位異常値処理方法は、ステップＳ１〜ステップＳ５を含む。
【００１６】
Ｓ１、ドローンのＩＭＵデータを取得する。
【００１７】
Ｓ２、加速度計とジャイロスコープを組み合わせて、ドローンのオイラー角を計算し、具体的には、
オイラー角及びジャイロスコープのドリフトを推定し、つまり、
【数２４】
2022500617.tif 000026
であり、ｂｇ，ｋはジャイロスコープのドリフト量、θ、φはそれぞれピッチ角、ロール角である。）のようにし、
【数２５】
2022500617.tif 000027
であり、
【数２６】
2022500617.tif 000028
はジャイロスコープのノイズであり、
離散化して、
【数２７】
2022500617.tif 000029
（式中、Ｔｓはサンプリング時間である。）を得て、
単純化して状態変数を分離し、
【数２８】
2022500617.tif 000030
を得て、
単純化には、Ｘｋ＝Φｋ−１Ｘｋ−１＋Γｋ−１ｗｋ−１とし、
ここで、ｋ−１時刻でのシステム遷移行列は、
【数２９】
2022500617.tif 000031
であり、
ｋ−１時刻でのシステムノイズ行列は、
【数３０】
2022500617.tif 000032
であり、
ｋ−１時刻での状態ベクトルは、
【数３１】
2022500617.tif 000033
であり、
ｋ−１時刻でのシステムノイズは、
【数３２】
2022500617.tif 000034
であり、
観測モデルは、
【数３３】
2022500617.tif 000035
であり、
ここで、
【数３４】
2022500617.tif 000036
であり、
Ｋ時刻での観測行列は、
【数３５】
2022500617.tif 000037
であり、
前記のように、単純化したプロセスモデルと観測モデルは、
【数３６】
2022500617.tif 000038
になり、
プロセスモデルと観測モデルを得た後、初期化
【数３７】
2022500617.tif 000039
（ここで、Ｅ（Ｘ０）は零時刻での数学的期待値、Ｅ［（Ｘ０−Ｅ（Ｘ０））（Ｘ０−Ｅ（Ｘ０））Ｔ］は零時刻の誤差共分散の計算を示す。）を行い、
状態推定予測Ｘｋ｜ｋ−１＝Φｋ−１Ｘｋ−１を行い、
誤差共分散予測：
Ｐｋ｜ｋ−１＝Φｋ−１Ｐｋ−１｜ｋ−１Φｋ−１Ｔ＋Γｋ−１Ｑｋ−１Γｋ−１Ｔ（Ｑｋ−１はシステムノイズ共分散行列を示す。）を行い、
カルマンゲイン行列は、
【数３８】
2022500617.tif 000040
（ここで、Ｈｋ＝［Ｏ３×３−ｇＩ］。ＲＴは測定ノイズ共分散行列である。）であり、
状態推定更新：
Ｘｋ｜ｋ＝Ｘｋ｜ｋ−１＋Ｋｋ（Ｚ−ＨｋＸｋ｜ｋ−１）
（ここで、Ｚは観測データである。）を行い、
誤差共分散更新：
Ｐｋ｜ｋ＝（Ｉ−ＫｋＨｋ）Ｐｋ｜ｋ−１
（Ｉは単位行列である。）を行い、状態ベクトルＸｋから正確なオイラー角を得る。
【００１８】
Ｓ３、ステップＳ２で得られたドローンのオイラー角から、Ｘ、Ｙ、Ｚ軸の加速度を計算し、具体的には、
【数３９】
2022500617.tif 000041
（ここで、ｘｋ（ｉ，ｊ）はｘｋの第ｉ行、第ｊ列のデータを表す。）から、ピッチ角θとロール角φを求め、ドローンの飛行速度が非常に小さく、すなわちその瞬時加速度を無視できるため、測定した加速度がｇ重力加速度に近く、測定した加速度の値ａｂから、
【数４０】
2022500617.tif 000042
△ｔ時間内のＸ、Ｙ、Ｚ軸上の変位を計算し、具体的には、式
【数４１】
2022500617.tif 000043
であり、
ここで、ｖｔは現在速度であり、ｖｔ＝（Ｐｔ−Ｐｔ−１）／Δｔは三次元ベクトルであり、Ｐｔはｔ時刻のｕｗｂ測位システムの位置ソリューションによる位置であり、前時刻の平均速度を現在時間の速度とし、また、
【数４２】
2022500617.tif 000044
はΔｔ時間内に零に近いので、
Δｓ＝ｖｔΔｔ＋ａｔΔｔであり、
式中、ａはＩＭＵソリューションにより得られる正確な加速度のナビゲーション座標系での加速度であり、ｔはシステムの動作時間であり、上記式よりｔ＋Δｔ時刻の推定変位ｓｔ＋Δｓを得て、
ここで、
【数４３】
2022500617.tif 000045
であり、Δｓ１、Δｓ２、Δｓ３はそれぞれサンプリング時間Δｔ内のナビゲーション座標系でのＸ、Ｙ、Ｚ軸の変位である。
ドローンの変位値Ｑを以下の式により推定し、
【数４４】
2022500617.tif 000046
であり、
ｘ、ｙ、ｚは前回測定した測位値であり、（ａ＊１、ａ＊２、ａ＊３）は＊番目の基地局の位置値である。
【００１９】
Ｓ４、ステップＳ３で得られた変位値Ｑ及びｕｗｂセンサの距離データから、ｕｗｂ測定値について異常検出を行い、ｕｗｂ測定値が異常である場合、ステップＳ５に入り、
異常値を検出して判断する際に、式は
【数４５】
2022500617.tif 000047
はそれぞれｕｗｂの現在の測定距離値と前の時刻の測定値であり、ＱはΔｔ時間内の変位値であり、ｍｅａｎはｎ個の測定データの標準偏差値であり、ｋは比重パラメータであり、具体的な環境パラメータに応じて調整することができ、式
【数４６】
2022500617.tif 000048
を測定異常値とする。
【００２０】
Ｓ５、異常なｕｗｂ測定値に対してデータ融合を行い、ｕｗｂ測定値を補正し、測定値
【数４７】
2022500617.tif 000049
であり、
ここで、ｓｔｄはｎ個の測定データの平均値である。
【００２１】
４つのｕｗｂ基地局がｕｗｂ測位を行う本実施例では、従来の生データの前処理過程における異常値の判別が不正確であること、及び異常値補正の精度などの問題を解決し、測定ハードウェアデータの正確性を向上させる。さらに、本実施例ハードウェア機器を追加する必要がなく、複雑な環境でデータを測定するときのｕｗｂセンサの安定性を高め、また、後続のｕｗｂ測位解に有効で正確な測定出力データを提供し、ｕｗｂ測位の安定性、効率及びリアルタイム性を高める。
【００２２】
以上の実施例は、本発明の好適な実施例に過ぎず、これらにより本発明の実施範囲を制限するものではなく、このため、本発明の形態、原理に基づいて行われる変化はすべて本発明の特許範囲に含まれるものとする。
【図１】
2022500617.tif 000074
【国際調査報告】
2022500617.tif 000075
2022500617.tif 000076
2022500617.tif 000077
2022500617.tif 000078