知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-07-26
(45)【発行日】2023-08-03
(54)【発明の名称】無線イベント検出及び監視のための方法、装置及びシステム
(51)【国際特許分類】
H04W 4/16 20090101AFI20230727BHJP
H04B 17/30 20150101ALI20230727BHJP
H04L 27/18 20060101ALI20230727BHJP
H04W 4/33 20180101ALI20230727BHJP
H04W 24/10 20090101ALI20230727BHJP
【ＦＩ】
H04W4/16
H04B17/30
H04L27/18 A
H04W4/33
H04W24/10
【請求項の数】
18
【外国語出願】
(21)【出願番号】P 2019212505
(22)【出願日】2019-11-25
(62)【分割の表示】P 2018529656の分割
【原出願日】2016-12-09
(65)【公開番号】P2020065259
(43)【公開日】2020-04-23
【審査請求日】2019-12-24
【審判番号】
【審判請求日】2022-08-23
(31)【優先権主張番号】62/265,155
(32)【優先日】2015-12-09
(33)【優先権主張国・地域又は機関】US
(31)【優先権主張番号】62/307,081
(32)【優先日】2016-03-11
(33)【優先権主張国・地域又は機関】US
(31)【優先権主張番号】62/316,850
(32)【優先日】2016-04-01
(33)【優先権主張国・地域又は機関】US
(31)【優先権主張番号】62/383,235
(32)【優先日】2016-09-02
(33)【優先権主張国・地域又は機関】US
(31)【優先権主張番号】62/411,504
(32)【優先日】2016-10-21
(33)【優先権主張国・地域又は機関】US
(73)【特許権者】
【識別番号】517016392
【氏名又は名称】オリジン ワイヤレス， インコーポレイテッド
【氏名又は名称原語表記】ＯＲＩＧＩＮ ＷＩＲＥＬＥＳＳ， ＩＮＣ．
【住所又は居所原語表記】Ｓｕｉｔｅ １０７０，７５００ Ｇｒｅｅｎｗａｙ Ｃｅｎｔｅｒ Ｄｒｉｖｅ，Ｇｒｅｅｎｂｅｌｔ，ＭＡＲＹＬＡＮＤ ２０７７０ Ｕ．Ｓ．Ａ．
(74)【代理人】
【識別番号】100076428
【弁理士】
【氏名又は名称】大塚 康徳
(74)【代理人】
【識別番号】100115071
【弁理士】
【氏名又は名称】大塚 康弘
(74)【代理人】
【識別番号】100112508
【弁理士】
【氏名又は名称】高柳 司郎
(74)【代理人】
【識別番号】100116894
【弁理士】
【氏名又は名称】木村 秀二
(74)【代理人】
【識別番号】100130409
【弁理士】
【氏名又は名称】下山 治
(74)【代理人】
【識別番号】100134175
【弁理士】
【氏名又は名称】永川 行光
(74)【代理人】
【識別番号】100166660
【弁理士】
【氏名又は名称】吉田 晴人
(72)【発明者】
【氏名】シュー， キンイ
(72)【発明者】
【氏名】チャン， フェン
(72)【発明者】
【氏名】チェン， チェン
(72)【発明者】
【氏名】ワン， ベイベイ
(72)【発明者】
【氏名】サファー， ゾルタン
(72)【発明者】
【氏名】ハン， イ
(72)【発明者】
【氏名】ライ， ハン－コック デュック
(72)【発明者】
【氏名】リウ， ケー． ジェイ． レイ
【合議体】
【審判長】中木 努
【審判官】圓道 浩史
【審判官】齋藤 哲
(56)【参考文献】
【文献】特開２００８－０２６３１０（ＪＰ，Ａ）
【文献】特開２０１４－２２８２９１（ＪＰ，Ａ）
【文献】特開２００５－０９１０２６（ＪＰ，Ａ）
【文献】国際公開第２０１６／０１１４３３（ＷＯ，Ａ２）
【文献】Ｋｕｎ ＱＩＡＮ ｅｔ ａｌ．，“ＰＡＤＳ： Ｐａｓｓｉｖｅ Ｄｅｔｅｃｔｉｏｎ ｏｆ Ｍｏｖｉｎｇ Ｔａｒｇｅｔｓ ｗｉｔｈ Ｄｙｎａｍｉｃ Ｓｐｅｅｄ ｕｓｉｎｇ ＰＨＹ Ｌａｙｅｒ Ｉｎｆｏｒｍａｔｉｏｎ”，２０１４ ２０ｔｈ ＩＥＥＥ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｐａｒａｌｌｅｌ ａｎｄ Ｄｉｓｔｒｉｂｕｔｅｄ Ｓｙｓｔｅｍｓ （ＩＣＰＡＤＳ），２０１４年１２月，ＤＯＩ： １０．１１０９／ＰＡＤＳＷ．２０１４．７０９７７８４
(58)【調査した分野】(Int.Cl.，ＤＢ名)
H04B 7/24- 7/26
H04W 4/00-99/00
G01S 7/00- 7/42, 13/00-13/95
3GPP TSG RAN WG1-4
3GPP TSG SA WG1-4
3GPP TSG CT WG1,4
(57)【特許請求の範囲】
【請求項1】
会場内で移動する物体を監視する監視システムのチャネル情報を用意するための方法であって、前記方法は、プロセッサ、メモリ及び前記監視システムの命令のセットを用いて、前記会場内の無線マルチパスチャネルを通じて前記監視システムの無線送信機から前記監視システムの無線受信機へ送信された少なくとも１つのプロービング信号に基づいて、チャネル情報（ＣＩ）の第１及び第２セットを取得することと、ここで、各ＣＩは、時間領域又は周波数領域におけるチャネル状態情報を含み、
ＣＩの前記第１及び第２セットの位相歪みを補正、補償、低減又は除去することでＣＩの前記第１及び第２セットを前処理することによって、ＣＩの前記第１及び第２セットをクリーニングすることと、
クリーニングされたＣＩの前記第１及び第２セットを解析することと、
クリーニングされたＣＩの前記第１及び第２セットの前記解析に基づいて、第１期間において前記物体の動作に対して取得されてクリーニングされた、ＣＩの前記第１セットと、ＣＩの前記第１セットよりも後の時間に取得されてクリーニングされた、ＣＩの前記第２セットとの時間反転共振強度を計算することによって、前記会場内における前記物体の前記動作を監視することと、
を含む方法。
【請求項2】
請求項１に記載の方法であって、サンプリング周波数オフセット（ＳＦＯ）及びシンボルタイミングオフセット（ＳＴＯ）により生じる、ＣＩの前記第１及び第２セットにおける線形位相オフセットを推定すること、
搬送周波数オフセット（ＣＦＯ）及び共通位相オフセットにより生じる、ＣＩの前記第１及び第２セットにおける初期位相歪みを推定すること、及び
ＣＩの前記第１及び第２セットにおける残留位相オフセットを除去すること、
のうちの少なくとも１つによって、ＣＩの前記第１及び第２セットをクリーニングすることを更に含む、方法。
【請求項3】
請求項２に記載の方法であって、前記線形位相オフセットを推定することは、ＣＩと前記ＣＩのシフトされたバージョンとの要素ごとの積に基づいて、線形位相係数を取得することを含む、方法。
【請求項4】
請求項１に記載の方法であって、前記会場内における前記物体の前記動作を監視することは、前記物体の前記動作に関連付けられたイベントの監視、前記イベントの検出、前記物体の前記動作に関連付けられた対象イベントの監視、屋内イベントの監視、正常イベントの監視、異常イベントの監視、ユーザ定義イベントの監視、前記物体の監視、前記動作に関連付けられた状態の監視、前記動作に関連付けられた状況の監視、前記会場内の人間の監視、前記人間の前記動作の監視、前記動作に関連付けられた状態遷移の検出、前記動作の検出、前記人間の動作の検出、あるロケーションにおける前記人間の動作の検出、あるロケーション付近の移動の検出、「人間が動いている」ことの検出、「人間が動いていない」ことの検出、存在の検出、不存在の検出、閉塞の検出、変化の検出、前記人間の動きの監視、人間の移動の監視、前記物体の呼吸動作の監視、継続的な移動の監視、人間の呼吸と人間の歩行と人間の運転とのうちの少なくとも１つの監視、前記人間の認識、前記物体の認識、前記イベントの認識、及び、人間と物体と扉と窓と部屋と領域と間取りとのうちの少なくとも１つに関連付けられたイベントの監視、のうちの少なくとも１つを含む、方法。
【請求項5】
請求項１に記載の方法であって、クリーニングされたＣＩの前記第１及び第２セットに基づく前記監視に基づいて、警報、信号、警報信号、報告、インジケーション、及び警告のうちの少なくとも１つをユーザへ送ることを更に含む、方法。
【請求項6】
請求項１に記載の方法であって、前記会場は、領域、監視領域、家屋、環境、オフィス、店舗、工場、ホテルの部屋、博物館、教室、倉庫、自動車、トラック、バス、船舶、列車、航空機、移動住宅、洞窟及びトンネル、のうちの少なくとも１つを含む、方法。
【請求項7】
請求項１に記載の方法であって、前記少なくとも１つのプロービング信号は、ハンドシェイク、交換、ＷｉＦｉパケットの交換、常時送出されるプローブ信号、複数のサブキャリアを有する直交周波数分割多重（ＯＦＤＭ）信号、及び、プリアンブルとショートプリアンブルとロングプリアンブルと符号とペイロードとトレーニング系列と有効トレーニング系列と疑似ランダム系列とのうちの少なくとも１つを有するフレーム、のうちの少なくとも１つを含む、方法。
【請求項8】
会場内で移動する物体を監視する監視システムのチャネル情報を用意するための装置であって、無線受信機と、
前記無線受信機と通信可能に接続されたプロセッサと、
前記プロセッサと通信可能に接続されたメモリと、
前記メモリに格納された命令のセットと、を備え、前記命令のセットは、前記プロセッサによって実行されると前記システムに、
前記会場内の無線マルチパスチャネルを通じて前記監視システムの無線送信機から前記無線受信機へ送信された少なくとも１つのプロービング信号に基づいて、チャネル情報（ＣＩ）の第１及び第２セットを取得することと、ここで、各ＣＩは、時間領域又は周波数領域におけるチャネル状態情報（ＣＳＩ）を含み、
ＣＩの前記第１及び第２セットの位相歪みを補正、補償、低減又は除去することでＣＩの前記第１及び第２セットを前処理することによって、ＣＩの前記第１及び第２セットをクリーニングすることと、
クリーニングされたＣＩの前記第１及び第２セットを解析することと、
クリーニングされたＣＩの前記第１及び第２セットの前記解析に基づいて、第１期間において前記物体の動作に対して取得されてクリーニングされたＣＩの前記第１セットと、前記第１期間より後の第２期間において取得されてクリーニングされたＣＩの前記第２セットとの時間反転共振強度を計算することによって、前記会場内における前記物体の前記動作を監視することと、
を行わせる、装置。
【請求項9】
請求項８に記載の装置であって、前記命令のセットは前記システムに更に、サンプリング周波数オフセット（ＳＦＯ）及びシンボルタイミングオフセット（ＳＴＯ）により生じる、ＣＩの前記第１及び第２セットにおける線形位相オフセットを推定すること、
搬送周波数オフセット（ＣＦＯ）及び共通位相オフセットにより生じる、ＣＩの前記第１及び第２セットにおける初期位相歪みを推定すること、及び
ＣＩの前記第１及び第２セットにおける残留位相オフセットを除去すること、
のうちの少なくとも１つによって、ＣＩの前記第１及び第２セットをクリーニングすることを行わせる、装置。
【請求項10】
請求項９に記載の装置であって、前記線形位相オフセットを推定することは、ＣＩと前記ＣＩのシフトされたバージョンとの要素ごとの積に基づいて、線形位相係数を取得することを含む、装置。
【請求項11】
請求項８に記載の装置であって、前記会場内における前記物体の前記動作を監視することは、前記物体の前記動作に関連付けられたイベントの監視、前記物体の前記動作に関連付けられた対象イベントの監視、屋内イベントの監視、正常イベントの監視、異常イベントの監視、ユーザ定義イベントの監視、前記物体の監視、前記動作に関連付けられた状態の監視、前記動作に関連付けられた状況の監視、前記会場内の人間の監視、前記人間の前記動作の監視、前記動作に関連付けられた状態遷移の検出、前記動作の検出、前記人間の動作の検出、あるロケーションにおける前記人間の動作の検出、あるロケーション付近の移動の検出、「人間が動いている」ことの検出、「人間が動いていない」ことの検出、存在の検出、不存在の検出、閉塞の検出、変化の検出、前記人間の動きの監視、人間の移動の監視、前記物体の呼吸動作の監視、継続的な移動の監視、人間の呼吸と人間の歩行と人間の運転とのうちの少なくとも１つの監視、前記人間の認識、前記物体の認識、前記イベントの認識、及び、人間と物体と扉と窓と部屋と領域と間取りとのうちの少なくとも１つに関連付けられたイベントの監視、のうちの少なくとも１つを含む、装置。
【請求項12】
請求項８に記載の装置であって、前記命令のセットは前記システムに更に、クリーニングされたＣＩの前記第１及び第２セットに基づく前記監視に基づいて、警報、信号、警報信号、報告、インジケーション、及び警告のうちの少なくとも１つをユーザへ送ることを行わせる、装置。
【請求項13】
請求項８に記載の装置であって、前記会場は、領域、監視領域、家屋、環境、オフィス、店舗、工場、ホテルの部屋、博物館、教室、倉庫、自動車、トラック、バス、船舶、列車、航空機、移動住宅、洞窟及びトンネル、のうちの少なくとも１つを含む、装置。
【請求項14】
請求項８に記載の装置であって、前記少なくとも１つのプロービング信号は、ハンドシェイク、交換、ＷｉＦｉパケットの交換、常時送出されるプローブ信号、複数のサブキャリアを有する直交周波数分割多重（ＯＦＤＭ）信号、及び、プリアンブルとショートプリアンブルとロングプリアンブルと符号とペイロードとトレーニング系列と有効トレーニング系列と疑似ランダム系列とのうちの少なくとも１つを有するフレーム、のうちのうちの少なくとも１つを含む、装置。
【請求項15】
会場内で移動する物体を監視する監視システムのチャネル情報を用意するための方法であって、プロセッサ、メモリ及び前記監視システムの命令のセットを用いて、前記会場内の無線マルチパスチャネルを通じて前記監視システムの無線送信機から前記監視システムの無線受信機へ送信された少なくとも１つのプロービング信号に基づいて、チャネル情報（ＣＩ）の第1及び第２セットを取得することと、
ここで、各ＣＩは、時間領域又は周波数領域におけるチャネル状態情報（ＣＳＩ）を含み、
ＣＩの前記第１及び第２セットの処理に基づいてＣＩの前記第１及び第２セットのダイナミクスを解析することと、
ここで、ＣＩの前記第１及び第２セットは、リモート処理される、クラウドリソースにより処理される、ローカルで処理される、及び、ＣＩの前記第１及び第２セットを生成するチップ又はシステムと一体化した又は近接したリソースにより処理される、のうちの少なくとも１つであり、
前記ＣＩの前記第１及び第２セットについての前記ダイナミクスの前記解析に基づいて、第１期間において前記物体の動作に対して取得されて、位相歪みを補正、補償、低減又は除去することでクリーニングされたＣＩの前記第１セットと、前記第１期間より後の第２期間において取得されて、位相歪みを補正、補償、低減又は除去することでクリーニングされたＣＩの前記第２セットとの間の時間反転共振強度を計算することによって、前記会場内における前記物体の前記動作を監視することと、
を含む方法。
【請求項16】
請求項１５に記載の方法であって、サンプリング周波数オフセット（ＳＦＯ）及びシンボルタイミングオフセット（ＳＴＯ）により生じる、ＣＩの前記第１及び第２セットにおける線形位相オフセットを推定すること、
搬送周波数オフセット（ＣＦＯ）及び共通位相オフセットにより生じる、ＣＩの前記第１及び第２セットにおける初期位相歪みを推定すること、及び
ＣＩの前記第１及び第２セットにおける残留位相オフセットを除去すること、
のうちの少なくとも１つによって、ＣＩの前記セットをクリーニングすることを更に含む、方法。
【請求項17】
請求項１６に記載の方法であって、前記線形位相オフセットを推定することは、ＣＩと前記ＣＩのシフトされたバージョンとの要素ごとの積に基づいて、線形位相係数を取得することを含む、方法。
【請求項18】
請求項１５に記載の方法であって、クリーニングされたＣＩの前記第１及び第２セットに基づく前記監視に基づいて、警報、信号、警報信号、報告、インジケーション、及び警告のうちの少なくとも１つをユーザへ送ることを更に含み、
前記会場を監視することは、イベントの監視、対象イベントの監視、前記会場におけるイベントの監視、屋内イベントの監視、正常イベントの監視、異常イベントの監視、ユーザ定義イベントの監視、物体の監視、前記物体の状態の監視、前記物体の状況の監視、前記会場内の人間の監視、状態遷移の検出、動作の検出、人間の動作の検出、あるロケーションにおける人間の動作の検出、あるロケーション付近の移動の検出、「人間が動いている」ことの検出、「人間が動いていない」ことの検出、存在の検出、不存在の検出、閉塞の検出、変化の検出、人間の動きの監視、人間の移動の監視、呼吸の監視、継続的な移動の監視、人間の呼吸と人間の歩行と人間の運転とのうちの少なくとも１つの監視、前記人間の認識、及び、人間と物体と扉と窓と部屋と領域と間取りとのうちの少なくとも１つに関連付けられたイベントの監視、のうちの少なくとも１つを含み、
前記少なくとも１つのプロービング信号は、ハンドシェイク、交換、ＷｉＦｉパケットの交換、常時送出されるプローブ信号、複数のサブキャリアを有する直交周波数分割多重（ＯＦＤＭ）信号、及び、プリアンブルとショートプリアンブルとロングプリアンブルと符号とペイロードとトレーニング系列と有効トレーニング系列と疑似ランダム系列とのうちの少なくとも１つを有するフレーム、のうちのうちの少なくとも１つを含む、方法。
【発明の詳細な説明】
【技術分野】
【0001】
関連出願の相互参照
本願は、２０１６年１２月９日に提出された米国特許出願第６２／２６５，１５５号、２０１６年３月１１日に提出された米国特許出願第６２／３０７，０８１号（００３３Ｐ０１）、２０１６年４月１日に提出された米国特許出願第６２／３１６，８５０号（００３３Ｐ０２）、２０１６年９月２日に提出された米国特許出願第６２／３８３，２３５号（００３６Ｐ０１）、及び２０１６年１０月２１日に提出された米国特許出願第６２／４１１，５０４号の優先権を主張する。上記の出願の全体の内容が本明細書において援用される。
本発明は、無線セキュリティ／監視の用途に関する。更に詳細には、本発明は、無線信号を使用する時間反転技術に基づいてイベント検出、動作検出を行う方法、装置及びシステムに関する。
本願は、２０１６年１２月９日に提出された米国特許出願第６２／２６５，１５５号、２０１６年３月１１日に提出された米国特許出願第６２／３０７，０８１号（００３３Ｐ０１）、２０１６年４月１日に提出された米国特許出願第６２／３１６，８５０号（００３３Ｐ０２）、２０１６年９月２日に提出された米国特許出願第６２／３８３，２３５号（００３６Ｐ０１）、及び２０１６年１０月２１日に提出された米国特許出願第６２／４１１，５０４号の優先権を主張する。上記の出願の全体の内容が本明細書において援用される。
本発明は、無線セキュリティ／監視の用途に関する。更に詳細には、本発明は、無線信号を使用する時間反転技術に基づいてイベント検出、動作検出を行う方法、装置及びシステムに関する。
【背景技術】
【0002】
現在の屋内監視システムは主に、屋内の物体の状態を監視するための接触センサと、人間の動作の存在を検出するためのパッシブ赤外線センサ（ＰＩＲ）に依存する。いくつかの例において、会場をセキュリティ保護するためには、複数のセンサを配置する必要があり、動作を検出するには見通し内（ＬＯＳ）環境が必要とされる。既存のシステムは、ユーザフィードバックの補助なしで学習や自己の更新を行わない。現在、大部分のホーム／オフィスセキュリティシステムは、制御パネル、接触センサ及びＰＩＲで構成される。接触センサは、接触が確立又は切断されたかを検出するために磁気又は電流を使用する装置であり、扉又は窓の開閉動作を監視するために扉及び窓に一般に設置されている。ＰＩＲセンサは、受信した赤外線波の変化を追跡することにより、環境内で熱を発生させる動く物体を検出する。いくつかの例において、接触センサは扉に取り付けられる必要があり、ＰＩＲセンサは動く物体を検出するために見通しを必要とするという点において、各センサの監視有効範囲は制限される。更に、会場を保護する際に検出するべきイベントが多い場合、コストは大幅に増加する。
【0003】
近年、無線信号を使用して屋内のイベント及び動作を検出する検知が注目を集めている。受信無線周波数（ＲＦ）信号は伝搬環境により変化する場合があるため、デバイスフリー屋内検知システムは、受信ＲＦ信号の変化を通じて環境内の動作を監視できる。屋内イベント検出のために信号送信中の変動を識別するのに使用できるＲＦ信号の特徴の例は、受信信号強度（ＲＳＳ）及びチャネル状態情報（ＣＳＩ）を含む。環境変化の影響を受けやすいため、ＲＳＳインジケータ（ＲＳＳＩ）は、屋内の動作を示すため及びそれを更に認識するために適用できる。チャネル状態情報は現在、いくつかの市販の装置において使用可能である。
【0004】
デバイスフリー屋内監視システムにおける技術の別のカテゴリは、対象の反射を使用して対象を追跡するレーダ技術を採用する。レーダ技術は、超広帯域（ＵＷＢ）検知を使用することにより、異なるパスを通る無線信号の飛行時間（ＴｏＦ）におけるサブナノ秒の遅延を識別できる。
【0005】
しかし、屋内モニタシステムの上述した技術は制限を有する。第１に、スカラとしての受信信号強度インジケータは自由度１しか有さず且つマルチパス効果の影響を強く受けるため、異なる屋内イベント及び物体を区別するための受信信号強度インジケータの分解能は低い。受信信号強度インジケータに基づくシステムにおいて、屋内検出の性能は、複数の検知装置又はアンテナを配置するコストを費やして保証される。更に、大部分のチャネル状態情報に基づく屋内検知システムは、チャネル状態情報の振幅のみに依存し、位相情報はどれほど有用であっても廃棄される。その一方で、異なるＴｏＦ情報を取得するために、レーダに基づく技術は１ＧＨｚ以上の帯域幅を使用して環境を検知するが、この帯域幅は市販のＷｉＦｉ装置で実現できず、センサから取得された結果は屋内イベントの種類を検出するために更なる労力を必要とする場合が多い。
【発明の概要】
【0006】
一般的態様において、複数のイベントを検出するための装置が提供される。前記装置は、プロセッサと、命令のセットを格納する記憶装置と、を備え、前記命令は、前記プロセッサによって実行された場合に前記プロセッサに、チャネル状態情報を分類し、かつ、前記チャネル状態情報の前記分類に基づいてイベントを検出するよう構成された少なくとも１つの分類器をトレーニングすることを行わせる。前記少なくとも１つの分類器の前記トレーニングは、前記装置によって検出される複数の既知のイベントのそれぞれについて、当該既知のイベントが会場内で発生する期間の間に、当該会場内の無線送信機と無線受信機との間の無線マルチパスチャネルのトレーニングチャネル状態情報（ＣＳＩ）の少なくとも１つのセットであって、前記無線送信機から前記無線マルチパスチャネルを介して前記無線受信機へ送信される少なくとも１つの第１プロービング信号から導出される、トレーニングチャネル状態情報の前記少なくとも１つのセットを取得すること、及び、前記複数の既知のイベントと、前記複数のイベントのそれぞれと関連付けられたトレーニングチャネル状態情報の前記少なくとも１つのセットとに基づいて、前記少なくとも１つの分類器をトレーニングすることを含む。前記プロセッサは更に、現在のイベントが発生している期間の間に、前記無線送信機と前記無線受信機との間の無線マルチパスチャネルの測定されたチャネル状態情報のセットであって、前記無線送信機から無線マルチパスチャネルを介して送信されて前記無線受信機で受信される少なくとも１つの第２プロービング信号から導出される、測定されたチャネル状態情報の前記セットを受信することと、前記複数の既知のイベントのいずれが前記現在のイベントにマッチングするかを判定するために、測定されたチャネル状態情報の前記セット又は測定されたチャネル状態情報の前記セットから導出されたデータの少なくとも１つに対して、前記分類器を適用することと、を行う。
【0007】
他の一般的態様において、推定されたチャネル状態情報（ＣＳＩ）に残留する位相オフセットを除去するための方法が提供される。前記方法は、デバイスから無線信号を受信することと、サンプリング周波数オフセット（ＳＦＯ）及びシンボルタイミングオフセット（ＳＴＯ）によって生じる、前記推定されたチャネル状態情報に含まれる線形位相オフセットを推定することと、搬送波周波数オフセット（ＣＦＯ）によって生じる、前記チャネル状態情報に含まれる共通位相オフセットを推定することと、いずれの参照チャネル状態情報も参照せずに、受信機におけるクリーンチャネル状態情報を取得することと、を含む。
【0008】
他の一般的態様において、ＷｉＦｉ信号を使用してイベントを監視するための方法が提供される。前記方法は、会場内で発生しうる複数の既知のイベントのそれぞれについて、当該既知のイベントが発生する少なくとも１つの期間の間に前記会場内で送信機から受信機へ送信される、ＩＥＥＥ８０２．１１標準規格に準拠したＷｉＦｉ信号の特性に関する第１情報を収集することと、前記会場を監視することと、を含み、前記会場を監視することは、第２期間の間に前記会場内で前記送信機から前記受信機へ送信されるＷｉＦｉ信号の特性に関する第２情報を判定することと、前記第２情報を前記第１情報と比較するか、又は前記第２情報から導出されるデータを前記第１情報から導出されるデータと比較することと、前記比較に基づいて、前記第２期間の間に、前記会場において前記複数の既知のイベントのいずれが発生しているかを判定することと、を含む。前記ＷｉＦｉ信号の特性に関する前記第１情報は、少なくとも、第１期間の間に前記送信機から前記受信機へ送信された前記ＷｉＦｉ信号から導出された、チャネル状態情報の第１セットを含み、前記ＷｉＦｉ信号の特性に関する前記第２情報は、少なくとも、前記第２期間の間に前記送信機から前記受信機へ送信された前記ＷｉＦｉ信号から導出された、チャネル状態情報の第２セットを含む。前記第２情報を前記第１情報と比較することは、チャネル状態情報の前記第２セットをチャネル状態情報の前記第１セットと比較することを含む。前記第２情報から導出された前記データを前記第１情報から導出された前記データを比較することは、チャネル状態情報の前記第２セットから導出されたデータをチャネル状態情報の前記第１セットから導出されたデータと比較することを含む。
【0009】
他の一般的態様において、ＷｉＦｉ信号を使用してイベントを監視するための方法が提供される。前記方法は、会場内で発生しうる複数の既知のイベントのそれぞれについて、当該既知のイベントが発生する少なくとも１つの期間の間に、ＩＥＥＥ８０２．１１標準規格に準拠した少なくとも１つのＷｉＦｉ信号を、前記会場内で送信機から受信機へ送信することと、前記受信機で受信されたＷｉＦｉ信号ごとに、当該ＷｉＦｉ信号からチャネル状態情報のセットを推定することと、ここで、チャネル状態情報の前記セットは、当該ＷｉＦｉ信号に関連付けられた前記既知のイベントに関連付けられており、主成分分析を使用して、前記複数の既知のイベントと、当該複数の既知のイベントに関連付けられたチャネル状態情報の前記セットとに基づいて、分類器をトレーニングすることと、を含む。前記方法は、監視段階の期間中に、第２期間の間に前記送信機から前記受信機へ少なくとも１つのＷｉＦｉ信号を送信することと、前記第２期間の間に前記受信機で受信された前記ＷｉＦｉ信号から、チャネル状態情報の第２セットを判定することと、前記第２期間の間に前記複数の既知のイベントのいずれが発生したかを判定するために、チャネル状態情報の前記第２セット又はチャネル状態情報の前記第２セットから導出されたデータに対して、前記分類器を適用することと、を含む。
【図面の簡単な説明】
【0010】
【0011】
【0012】
【0013】
【0014】
【0015】
【0016】
【0017】
【0018】
【0019】
【0020】
【0021】
【0022】
【0023】
【0024】
【0025】
【0026】
【0027】
【0028】
【0029】
【0030】
【0031】
【0032】
【0033】
【0034】
【0035】
【0036】
【0037】
【0038】
【0039】
【0040】
【0041】
【0042】
【0043】
【0044】
【0045】
【0046】
【0047】
【0048】
【0049】
【0050】
【0051】
【0052】
【0053】
【0054】
【0055】
【0056】
【0057】
【0058】
【0059】
【0060】
【0061】
【0062】
【0063】
【発明を実施するための形態】
【0064】
本明細書において、屋内会場であってもよく又は無線信号を反射して豊富なマルチパス環境を提供できる複数の表面を有する屋外会場であってもよい会場内のイベントを監視又は検出するシステム及び方法を説明する。会場は、例えば部屋、家屋、オフィス、店舗、工場、ホテルの部屋、博物館、教室、倉庫、自動車、トラック、バス、船舶、列車、航空機、移動住宅、洞窟又はトンネルであってもよい。例えば会場は、１つ又は複数の階を有する建物であってもよく、建物の一部分は地下であってもよい。建物の形状は、例えば円形、正方形、矩形、三角形又は不規則な形状であってもよい。これらは単なる例であり、本発明は他の種類の会場又は空間におけるイベントを検出するために使用可能である。
【0065】
送信機が会場においてマルチパスチャネルと呼ばれる複数のパスを介して受信機に無線信号を送出する場合、マルチパスチャネルは会場内の物体又は表面の配置等の会場の特徴の影響を受ける。従って、送信機及び受信機が会場内の固定位置に配置され且つ送信機が同一の波形を有する信号を受信機に繰り返し送出する場合、会場の特徴が時間と共に変化すると、受信機において受信される信号の波形も時間と共に変化する場合がある。受信機において受信された信号の波形の特徴を測定することにより、会場内で起きているイベントを推測することができる。
【0066】
いくつかの実現例において、システムはトレーニング段階で動作し、次に監視段階で動作する。トレーニング段階において、ユーザが監視することを意図する特定のイベントを表すために参照値を使用できるように特定のパラメータの参照値を確立するために、分類器がトレーニングされる。イベントは、例えば「扉が開く」、「扉が閉じる」、「窓が開く」、「窓が閉じる」、「室内で人間が動いていない」又は「室内で人間が動いている」であってもよい。これらは単なる例であり、システムは他の種類のイベントを検出することもできる。分類器は、イベントが会場内で発生している時に測定された無線信号から導出された情報を使用してトレーニングされる。監視段階において、会場内の無線信号が測定され、無線信号から導出された情報が分類器に提供される。測定された無線信号から導出された情報及び参照値に基づいて、分類器は会場内で発生したイベントを判定する。
【0067】
以下の説明において、環境固有シグニチャ及び時間反転無線システムの概要を提供する。次に、ベースライン法、統計的方法、機械学習法、時間反転共振強度の分散に基づく方法、アンテナ相関に基づく方法及び到着時間に基づく方法を使用する分類器等の分類器の種々の実現例を説明する。
【0068】
環境固有シグニチャ
図３５は、２つの送受信機２１０及び２１２を備える無線システム２０８の例示的な一実施形態を示す。本実施形態において、アンテナを備える送受信機Ａ２１０は無線信号２１４を発射し、当該信号は無線チャネル２１６を伝搬し、アンテナを備える送受信機Ｂ２１２にマルチパス無線信号２１８として到着する。例示的な実施形態において、少なくとも１つのアンテナは少なくとも１つの無線信号をチャネルに発射してもよく、少なくとも１つのアンテナは無線チャネルから信号を受信してもよい。実施形態において、送信アンテナ及び受信アンテナは互いに離間して配置されてもよく、いくつかの実施形態において、それらは同一場所に設置されてもよい。例えば装置、コンピュータ、モバイル装置及びアクセスポイント等は２つ以上のアンテナを備えてもよく、アンテナは送信アンテナ及び受信アンテナのいずれか又は双方として動作されてもよい。いくつかの実施形態において、少なくとも１つのアンテナは、チャネルへの無線信号の発射及びチャネルからのマルチパス信号の受信の双方に使用されてもよい単一アンテナであってもよい。実施形態において、アンテナは、異なる時間スロット、異なる周波数帯域、異なる方向及び／又は異なる偏波で信号を送受信してもよく、あるいは同一又は同様の時間、同一又は同様の周波数帯域、同一又は同様の方向及び／又は同一又は同様の偏波で信号を送受信してもよい。いくつかの実施形態において、アンテナ及び／又はアンテナを備える装置は、信号送信及び信号受信のタイミング、搬送周波数、方向及び／又は偏波を調節してもよい。
図３５は、２つの送受信機２１０及び２１２を備える無線システム２０８の例示的な一実施形態を示す。本実施形態において、アンテナを備える送受信機Ａ２１０は無線信号２１４を発射し、当該信号は無線チャネル２１６を伝搬し、アンテナを備える送受信機Ｂ２１２にマルチパス無線信号２１８として到着する。例示的な実施形態において、少なくとも１つのアンテナは少なくとも１つの無線信号をチャネルに発射してもよく、少なくとも１つのアンテナは無線チャネルから信号を受信してもよい。実施形態において、送信アンテナ及び受信アンテナは互いに離間して配置されてもよく、いくつかの実施形態において、それらは同一場所に設置されてもよい。例えば装置、コンピュータ、モバイル装置及びアクセスポイント等は２つ以上のアンテナを備えてもよく、アンテナは送信アンテナ及び受信アンテナのいずれか又は双方として動作されてもよい。いくつかの実施形態において、少なくとも１つのアンテナは、チャネルへの無線信号の発射及びチャネルからのマルチパス信号の受信の双方に使用されてもよい単一アンテナであってもよい。実施形態において、アンテナは、異なる時間スロット、異なる周波数帯域、異なる方向及び／又は異なる偏波で信号を送受信してもよく、あるいは同一又は同様の時間、同一又は同様の周波数帯域、同一又は同様の方向及び／又は同一又は同様の偏波で信号を送受信してもよい。いくつかの実施形態において、アンテナ及び／又はアンテナを備える装置は、信号送信及び信号受信のタイミング、搬送周波数、方向及び／又は偏波を調節してもよい。
【0069】
例示的な実施形態におけるアンテナは、電気出力又は電気信号を電波、マイクロ波、マイクロ波信号又は無線信号に変換し且つ電波、マイクロ波、マイクロ波信号又は無線信号を電気出力又は電気信号に変換するあらゆる種類の電気装置であってもよい。例えば少なくとも１つのアンテナは指向性アンテナ又は全方向性アンテナとして構成されてもよいが、それらに限定されない。少なくとも１つのアンテナは、何らかの種類のモノポールアンテナ、ダイポールアンテナ及びクワドラポールアンテナ等であってもよい。少なくとも１つのアンテナは、何らかの種類のループアンテナであってもよく且つ／又は１本のワイヤから形成されてもよい。少なくとも１つのアンテナは、パッチアンテナ、パラボラアンテナ、ホーンアンテナ、八木アンテナ、折り返しダイポールアンテナ、マルチバンドアンテナ、短波アンテナ、マイクロ波アンテナ、同軸アンテナ、メタマテリアルアンテナ、衛星アンテナ、誘電共振アンテナ、フラクタルアンテナ、ヘリカルアンテナ、等方性放射体、Ｊポールアンテナ、スロットアンテナ、マイクロストリップアンテナ、コンフォーマルアンテナ、皿型アンテナ、テレビアンテナ、ラジオアンテナ、ランダムワイヤアンテナ、セクタアンテナ、セルラアンテナ、スマートアンテナ及び傘形アンテナ等であってもよい。少なくとも１つのアンテナは、直線アレイアンテナ、フェーズドアレイアンテナ、反射アレイアンテナ及び指向性アレイアンテナ等のアンテナアレイの一部であってもよい。少なくとも１つのアンテナは、狭帯域アンテナ又は広帯域アンテナ、高利得アンテナ又は低利得アンテナ、調節可能アンテナ又はチューナブルアンテナ、あるいは固定アンテナであってもよい。あらゆる種類のアンテナは、本明細書中で説明するシステム、方法及び技術において使用されるように構成されてもよい。実施形態において、例示的なアンテナに関連する放射パターンは調整可能であってもよく、本明細書中で説明する例示的なシステム、方法及び技術の性能を向上するために調整されてもよい。
【0070】
実施形態において、電気信号は無線送信のために１つ以上のアンテナに適用されてもよく、処理するために１つ以上のアンテナから受信されてもよい。実施形態において、無線信号は電波又はマイクロ波であってもよい。実施形態において、無線信号はキロヘルツからテラヘルツの範囲のいずれかの搬送周波数を有してもよい。実施形態において、アンテナは、フィルタ、増幅器、スイッチ、モニタポート及びインピーダンス整合回路等のうちの少なくとも１つを備えてもよい。実施形態において、電気信号はアナログ回路網及び／又はデジタル回路網を使用して生成されてもよく、少なくとも１つのアンテナを駆動するために使用されてもよい。実施形態において、少なくとも１つのアンテナから受信された電気信号は、アナログ回路網及び／又はデジタル回路網を使用して処理されてもよい。本明細書中で開示される本発明の例示的な実施形態において、電気信号は、サンプリング、デジタル化、格納、比較、相関、時間反転、増幅、減衰、調節、補償、統合及び処理等が行われてもよい。
【0071】
本開示において、チャネルの特性をプロービングするために送信アンテナにより発射される信号は、場合によってはプローブ信号、チャンネルプローブ信号又はチャネルプローブ波形と呼ばれてもよい。図３６は、第１の装置２１０から広帯域無線チャネル２１６を介して第２の装置２１２へ送信される無線信号２１４を示す。チャネルプローブ信号２１４は、受信プローブ波形２１８とも呼ばれるものとして第２の装置２１２に到着してもよい。当該受信プローブ波形２１８は、少なくとも１つのアンテナと受信機電子機器のセットとを備える受信機により受信され且つ処理されてもよい。例示的な実施形態において、受信プローブ波形２１８を処理することは、装置２１０と装置２１２との間の広帯域チャネルに対する推定チャネル応答を生成してもよい。実施形態において、プローブ信号及び受信信号はデジタル信号に変換されるアナログ信号であってもよく（また、アナログ信号に変換されるデジタル信号であってもよく）、デジタル信号プロセッサ（ＤＳＰ），フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、ＡＲＭ（Advanced RISC Machine）プロセッサ、マイクロプロセッサ、コンピュータ及び特定用途向け集積回路（ＡＳＩＣ）等を使用して処理され且つ／又は生成されてもよい。
【0072】
時間領域において、通信リンクのチャネルインパルス応答は
としてモデル化可能である。式中、ｈｉ［ｋ］は長さＬを有するチャネルインパルス応答（ＣＩＲ）のｋ番目のタップであり、δ［］はディラックのデルタ関数である。尚、チャネル応答の時間領域表現ｈ及びチャネル応答の周波数領域表現Ｈはフーリエ変換により関連付けられる。
【0073】
例示的な実施形態において、チャネルインパルス応答が既知である場合、受信プローブ波形は、チャネルプローブ信号とチャネルインパルス応答とを畳み込むことにより予測されてもよい。チャネルインパルス応答又は推定チャネル応答は、実際のチャネルインパルス応答の近似値又は推定値であってもよい。例えば推定チャネル応答は、チャネルの「十分に正確な」推定値であると考えられるか又はチャネルの特定の特性を優先的にプロービングするように選択される特定のチャネル長に切り捨てられてもよい。更に、推定チャネル応答は、特定の用途に対して「十分に正確である」と判定された離散化信号の時間／振幅分解能を用いて受信プローブ波形の離散化近似値から導出されてもよい。推定チャネル応答は、実際のチャネル応答のフィルタリングされたバージョンであってもよく、チャネルの十分に正確な推定値であってもよい。何が「十分に正確である」かの判定は、用途、無線装置で使用されるハードウェア構成要素、装置の処理能力、装置の許容電力消費及びシステム性能の所望の精度等に依存してもよい。
【0074】
装置により送信されるプローブ信号が単一のパルス信号又はインパルス信号である場合、受信プローブ波形はチャネルインパルス応答の十分に正確な推定値であってもよく、推定チャネル応答を取得するために受信プローブ波形の受信、離散化及び格納以外の追加処理は殆ど不要であってもよい。装置により送信されるプローブ信号が単一のパルス信号又はインパルス信号以外の波形である場合、受信機は、推定チャネル応答を判定するために受信プローブ波形に追加処理を実行する必要がある場合がある。例示的な一実施形態において、受信器が受信プローブ波形を検出し且つ離散化してもよい。アナログデジタル（Ａ／Ｄ）変換器が離散化を実行するために使用されてもよい。実施形態において、デコンボリューション処理は、離散化された受信プローブ波形とチャネルプローブ信号の表現とを使用して推定チャネル応答を得てもよい。実施形態において、推定チャネル応答を得るために他の数学関数が使用されてもよい。チャネルインパルス応答（ＣＩＲ）は、本明細書においてチャネル応答（ＣＲ）、ＣＲ信号、ＣＩＲ信号、チャネルプローブ信号応答及び推定チャネル応答と呼ばれてもよい。チャネル応答は測定され且つ／又は計算されてもよく、且つ／又は測定及び計算の組み合わせにより生成されてもよい。本開示において、チャネル応答及び受信プローブ波形をロケーション固有シグニチャと呼ぶ場合もある。
【0075】
実施形態において、広帯域チャネルのチャネル応答の推定値の精度を向上又は低下するために、異なるチャネルプローブ信号が選択されてもよい。例示的な実施形態において、チャネルプローブ信号はパルス又はインパルスであってもよい。更に、チャネルプローブ信号は、規則的なパターン、任意のパターン又は不規則的なパターンを有する一連のパルスであってもよい。チャネルプローブ信号は波形であってもよい。波形は、実質的に矩形波形、二乗余弦波形、ガウス波形、ローレンツ波形、あるいは何らかの最適又は所望の方法でチャネルをプロービングするように設計された形状を有する波形であってもよい。例えばチャネルプローブ波形は、周波数がチャープされてもよく、あるいは何らかの最適又は所望の方法でチャネルをプロービングするように適合される周波数スペクトルを有してもよい。プローブ波形は、異なる中心周波数及び帯域幅を有する複数の波形であってもよい。プローブ波形は、振幅変調、位相変調、周波数変調、パルス位置変調又は偏波変調されてもよく、あるいは振幅、位相、周波数、パルス位置及び偏波の何らかの組み合わせで変調されてもよい。
【0076】
波形は、関連する通信チャネルを介して交換されることを意図してもよいデータストリームのビット持続時間とほぼ等しい時間幅を有してもよい。波形は、関連する通信チャネルを介して交換されることを意図するデータストリームのビット持続時間の略半分、略４分の１、略１０分の１、略１００分の１又はそれ以下の時間幅を有してもよい。プローブ信号／波形はデータパターンであってもよく、反復するデータパターンであってもよい。プローブ信号は、パケット及び／又はフレーミング情報、同期及び／又はクロック再生情報、ストリーム取り込み情報、装置ＩＤ／ネットワーク／リンク層動作情報を含んでもよい。プローブ信号は、動作環境、並びに／又はシステムの送信機及び／又は受信機内の電子構成要素に対して適合された周波数スペクトルを有してもよい。プローブ信号はチャネルインパルス応答の推定値であってもよく、あるいはチャネルインパルス応答の推定値の変更されたバージョンであってもよい。例えばプローブ信号は、推定チャネル応答の時間反転バージョンであってもよい。プローブ信号は、送信機及び／又は受信機内の特定の電子構成要素及び／又は特定の環境要素により生じる信号歪みを補償し且つ／又は強調するように設計されてもよい。
【0077】
チャネルプロービング信号の１つの例示的な種類は、周期パルス系列である。そのようなチャネルプロービング信号を用いる場合、受信プローブ波形は周期的チャネルパルス応答のノイズバージョンであってもよい。実施形態において、ノイズを抑制してチャネル応答を抽出するために時間平均化法を使用できる。
【0078】
いくつかの実施形態において、時間平均化法はチャネル応答の信頼できる測度を提供しない場合がある。チャネル応答推定を向上するためには、ノイズを抑制するために更に長いパルス系列を使用できる。システムの性能を更に向上するために、短い擬似ランダムパルス系列をチャネルプロービング信号として使用できる。そのような場合、受信プローブ波形は、擬似ランダム系列とチャネル応答との畳み込みであってもよい。
【0079】
実施形態において、プロービング信号として使用される擬似ランダム系列は、受信機が既知であってもよい。その場合、チャネル応答は、受信信号が擬似ランダム系列と畳み込まれる相関に基づく方法を使用して推定可能である。一般に、擬似ランダム系列の自己相関は、シンボル間干渉が生じて推定チャネル応答に誤差が存在する場合があるため、理想的なデルタ関数でない場合がある。実施形態において、シンボル間干渉によるそのような種類のチャネル推定誤差は、擬似ランダム系列でなく、自己相関関数に対する理想的なデルタ形状を有してもよい直交ゴレイ相補系列を使用することにより最小化又は回避されてもよい。
【0080】
実施形態において、無線装置は、ｆ１ＧＨｚの中心周波数を有する第１の無線信号を送信してもよい。実施形態において、第１の無線信号はチャネルプローブ信号、パルス信号、フレーム信号、擬似ランダムノイズ（ＰＮ）系列及びプリアンブル信号等であってもよい。実施形態において、無線信号の帯域幅は、約１０ＭＨｚ、約２０ＭＨｚ、約４０ＭＨｚ、約６０ＭＨｚ、約１２５ＭＨｚ、約２５０ＭＨｚ、約５００ＭＨｚ及び約１ＧＨｚ等であってもよい。実施形態において、無線装置はｆ２ＧＨｚの中心周波数を有する第２の無線信号を送出してもよい。実施形態において、第２の無線信号はチャネルプローブ信号、パルス信号、フレーム信号、ＰＮ系列及びプリアンブル信号等であってもよい。実施形態において、無線信号の帯域幅は、約１０ＭＨｚ、約２０ＭＨｚ、約４０ＭＨｚ、約６０ＭＨｚ、約１２５ＭＨｚ、約２５０ＭＨｚ、約５００ＭＨｚ及び約１ＧＨｚ等であってもよい。実施形態において、第１の無線信号の周波数スペクトル及び第２の無線信号の周波数スペクトルはオーバラップする周波数を含んでもよい。いくつかの実施形態において、２つの無線信号の間にオーバラップする周波数が存在しなくてもよい。いくつかの実施形態において、異なる無線信号の周波数スペクトルは、いわゆるガードバンド又はガードバンド周波数により分離されてもよい。第１の無線信号（例えば、周波数ｆ１における）を使用してプロービングされたチャネルに対するチャネル応答はＨｉｊ（ｆ１）と表されてもよい。第２の無線信号（例えば、プローブ周波数ｆ２における）を使用してプロービングされたチャネルに対するチャネル応答はＨｉｊ（ｆ２）と表されてもよい。実施形態において、３つ以上のプローブ周波数信号がチャネルをプロービングするために使用されてもよい。３つ以上のプローブ周波数信号は何らかのオーバラップする周波数を有してもよく、あるいはオーバラップする周波数を有さなくてもよい。
【0081】
実施形態において、無線装置は、異なる無線信号搬送周波数に同調して無線チャネルをプロービングするために、チャネル同調及び／又は周波数ホッピングを使用してもよい。いくつかの実施形態において、無線装置は、指定された周波数帯域内の異なるチャネルに同調して無線チャネルをプロービングしてもよい。例えば無線装置は、最初にＷｉＦｉ（ＩＥＥＥ８０２．１１）信号伝送帯域幅内の１つのチャネルに同調した後に無線帯域内の別のチャネルに同調してもよい。周波数同調は、１つのチャネルから次のチャネルに順次行われてもよいが、ＷｉＦｉ帯域内で１つのチャネルから別のチャネルにランダムにホップしてもよい。実施形態において、異なるチャネルは異なるチャネル帯域幅を有してもよい。実施形態において、プローブ信号を生成するため及び／又は受信信号内のチャネル情報を解析するために、あらゆる無線プロトコルが使用されてもよい。
【0082】
実施形態において、チャネルをプロービングするために複数のチャネルプローブ信号が使用されてもよい。いくつかの実現例において、同一のプローブ信号が複数回送出されてもよく、受信プローブ波形が平均化され且つ／又は比較されてもよい。例えばプローブ信号は２回、５回、１０回、３０回、５０回、１００回、５００回又は１０００回送出されてもよい。実施形態において、プローブ信号は１回送出されてもよく、２回～１０００回のいずれかの回数送出されてもよい。実施形態において、プローブ信号は１００１回以上送出されてもよい。例えば何らかの監視及びセキュリティ用途において、プローブ信号は継続的に送出されてもよい。空間を監視してプロービングするために、例えば毎秒１個、毎秒１０個、毎秒１００個等のプローブ信号が継続的に送出されてもよい。プローブ信号が継続的に送出される速度は、環境の変化が検出される必要のある速度により判定されてもよい。
【0083】
実施形態において、一部の受信プローブ波形のみが更なる処理に使用されてもよい。例えば一部の受信プローブ波形及び／又は推定チャネル応答が廃棄又は削除されてもよい。廃棄及び／又は削除される波形及び／又は応答は、他の受信波形及び／又は推定応答と十分に異なるため、外れ値（outlier）でありチャネルの十分に正確な表現ではないと考えられるものであってもよい。いくつかの実施形態において、異なるプローブ信号は異なる時間に送出されてもよく且つ／又は受信機からのフィードバックに応答して送出されてもよい。例えば送信機におけるプローブ信号は、受信プローブ波形、推定チャネル応答、並びに／又は受信プローブ波形及び／又は推定チャネル応答の類似性を向上するように調整されてもよい。実施形態において、送信機は少なくとも２つの異なるプローブ信号を送出してもよく、受信機は少なくとも２つの異なる受信プローブ波形のいずれか１つ、いくつか又は全てに基づいてチャネル応答を推定してもよい。実施形態において、プローブ信号は、以前に測定され且つ／又は計算されたチャネル応答のバージョン及び／又は測定され且つ／又は計算されたチャネル応答の時間反転バージョンであってもよい。
【0084】
本開示中で以下に更に詳細に説明するように、波形、シグニチャ及び／又は応答の類似性、マッチング又は相関は、仮想時間反転処理技術、時間反転共振強度、パターン認識及び／又はパターンマッチング、線形及び／又は非線形のサポートベクトルマシン及び／又はサポートベクトルネットワーク、機械学習、データマイニング、分類、統計的分類、タギング及びカーネルトリック（例えば、カーネル関数を適用するカーネル法）等を使用して判定されてもよい。
【0085】
実施形態において、受信プローブ波形を処理することは、受信信号のいずれか一部分を増幅又は減衰することを含んでもよい。実施形態において、チャネルは１回プロービングされてもよく、あるいはチャネルは２回以上プロービングされてもよい。実施形態において、複数の受信プローブ波形の測定、処理及び記録等が行われてもよい。実施形態において、いくつかのチャネル応答が他のチャネル応答と平均化されてもよい。実施形態において、いくつかのチャネル応答が廃棄されてもよく又は記録されなくてもよい。実施形態において、いくつかのチャネル応答が異なる環境条件で測定されて格納されてもよい。そのような格納された応答信号は、元の測定値に関連する環境条件を示すための参照信号として使用されてもよい。実施形態において、新規に測定されたチャネル応答は、以前に格納されたチャネル応答の中で新規に測定されたチャネル応答と最もマッチングするものを判定するために、多くの以前に格納されたチャネル応答と比較されてもよい。その場合、最も密接に相関するか又は最もマッチングする、以前に格納されたチャネル応答の環境パラメータが、新規に測定されたチャネル応答に関連付けられてもよい。例示的な実施形態において、環境条件は、温度、物体のロケーション又は配置、人間のロケーション及び／又は配置、物体の姿勢、人間の姿勢、アクセスポイントのロケーション及び／又は姿勢、端末装置、センサの位置及び／又は姿勢、信号反射体の位置及び／又は姿勢、信号散乱体の位置及び／又は姿勢、並びに信号減衰体の位置及び／又は姿勢等を含んでもよいが、それらに限定されなくてもよい。
【0086】
例示的な一実施形態において、推定チャネル応答は、特定の環境における２つの装置間のチャネル応答、あるいは会場又は特定の環境における装置と物体及び／又は構造との間のチャネル応答を表すため、環境固有波形及び／又はシグニチャであると考えられてもよい。図３６に示すように、装置２１０と装置２１２との間で送信された信号が伝搬する会場又は環境内の１つ以上の物体及び／又は構造及び／又は表面において１つ以上の動きが存在する場合、信号が伝搬する複数の伝搬パスのうちの少なくともいくつかが変化することによりチャネル応答が変化する可能性がある。推定チャネル波形の特性及びその変化量は、会場、環境及びシステム内のハードウェア構成要素に依存してもよい。
【0087】
時間反転無線システムの概要
時間反転無線システムの概要を以下に提供する。図１Ａを参照すると、時間反転システムは、無線信号を１回以上反射させる場合がある構造又は物体を有する環境において使用可能である。例えば会場１０２は、第１の部屋１０４及び第２の部屋１０６を有してもよい。第１の部屋１０４内の第１の装置１０８が第２の部屋１０６内の第２の装置１１０へ信号を送信する場合、信号は複数の方向に伝搬し、例えば１１２、１１４及び１１６である複数の伝搬パスを通って移動することにより第２の装置１１０に到達できる。複数の伝搬パスを通って移動する信号をマルチパス信号と呼ぶ。信号が伝搬パスを通って移動する際、信号は歪む場合がある。第２の装置１１０により受信されるマルチパス信号は、第１の装置１０８により送信された信号と非常に異なる場合がある。
時間反転無線システムの概要を以下に提供する。図１Ａを参照すると、時間反転システムは、無線信号を１回以上反射させる場合がある構造又は物体を有する環境において使用可能である。例えば会場１０２は、第１の部屋１０４及び第２の部屋１０６を有してもよい。第１の部屋１０４内の第１の装置１０８が第２の部屋１０６内の第２の装置１１０へ信号を送信する場合、信号は複数の方向に伝搬し、例えば１１２、１１４及び１１６である複数の伝搬パスを通って移動することにより第２の装置１１０に到達できる。複数の伝搬パスを通って移動する信号をマルチパス信号と呼ぶ。信号が伝搬パスを通って移動する際、信号は歪む場合がある。第２の装置１１０により受信されるマルチパス信号は、第１の装置１０８により送信された信号と非常に異なる場合がある。
【0088】
例えば図１Ｂを参照すると、第１の装置１０８がパルス信号を送出する場合、第２の装置１１０により受信される信号は波形１２０を有してもよい。波形１２０をチャネルインパルス応答（ＣＩＲ）と呼び、これは時間領域におけるチャネル状態情報（ＣＳＩ）を反映する。本明細書において、「チャネル状態情報」（ＣＳＩ）は通信リンクのチャネル特性を示す。これは、信号が送信機から受信機に伝搬する様子を説明する。チャネル状態情報は、受信チャネルサウンディング信号又は受信チャネルプロービング信号に基づくチャネル推定を使用して取得可能である。時間領域におけるチャネル状態情報は、チャネルインパルス応答（ＣＩＲ）と呼ばれることが多く、チャネル遅延スプレッドの間のチャネル減衰及び位相シフトを反映する複素値のベクトルにより表すことができる。周波数領域におけるチャネル状態情報は、チャネル周波数応答（ＣＦＲ）と呼ばれることが多く、例えばＷｉＦｉチャネル内のサブキャリアである各周波数ユニットにおけるチャネル減衰及び位相シフトを反映する複素値のベクトルにより表すことができる。
【0089】
図１Ｃを参照すると、時間反転波形（シグニチャ）１３０は、波形１２０を時間に関して反転することにより生成できる。第２の装置１１０が波形１３０を有する信号を送出する場合、信号は伝搬パス１１２、１１４及び１１６を含む種々の方向に伝搬し（インパルス信号の伝搬方向に対して逆方向に）、第１の装置１０８に到達する。第１の装置１０８において受信されるマルチパス信号は、第１の装置１０８から第２の装置１１０に以前に送出されたインパルス信号と同様のインパルス信号を形成してもよい。
【0090】
図１Ｂ及び図１Ｃに示す波形１２０及び１３０は単なる例である。時間反転システムにおける波形は、例えば環境及び送信される情報又はデータに依存して変化する場合がある。更に、第１の装置１０８から送出される初期信号は、あらゆる種類の信号であってもよく、インパルスは一例にすぎない。初期信号はあらゆる波形であってもよい。
【0091】
いくつかの実現例において、第２の装置１１０が第１の装置１０８へデータストリームを送信することを意図する場合、第２の装置１１０は装置１０８から受信した信号の正規化時間反転共役を基本送信波形（シグニチャ）として使用してもよい。第２の装置１１０は、基本波形（シグニチャ）においてデータストリームを符号化し、無線チャネルを介して信号を送信してもよい。装置１０８において受信される信号は、送信信号とチャネルインパルス応答との畳み込みに加法性白色ガウスノイズが加わったものとして説明されてもよい。送信信号はチャネルインパルス応答の時間反転バージョンに基づいて設計されているため、第１の装置１０８は、受信信号に対して単純な調整を実行してダウンサンプリングするだけで、第２の装置１１０により送信されたデータストリームを再生できる。
【0092】
いくつかの例において、送信装置又は基地局又はアクセスポイントは、２つ以上の受信装置に信号を同時に送出してもよい。送信信号は、複数の伝搬パスを通って各受信機に移動する。受信機が異なるロケーションに位置するため、マルチパス信号は異なる伝搬パスを通って受信機に到達する。送信機から送出される信号の波形を慎重に構成することにより、各受信機が自身に向けられたデータを十分に高品質で受信できるようにすることが可能である。
【0093】
チャネルプロービング
いくつかの実現例において、図３７を参照すると、チャネルプロービングは送受信機又はアクセスポイント（ＡＰ）５９０が端末装置５９２と通信する場合に実行される。チャネルプロービング段階において、端末装置５９２は、アップサンプリング部５９６を使用してインパルス信号５９４をαでアップサンプリングし、ルートレイズドコサイン（ＲＲＣ）フィルタＲＲＣB,fs５９８を使用してアップサンプル信号をフィルタリングする。この場合、Ｂは帯域幅であり、ｆｓ＝αＢはサンプルレートである。ビットレートを上げる目的の１つは、オーバサンプリングすることにより、データ送信段階において送信機と受信機との間にサンプル周波数オフセットが存在する場合に正確なサンプルを取得できることである。オーバサンプリングはデータ送信段階におけるステップの１つであってもよいため、チャネルを正確に推定するために、チャネルプロービング段階にオーバサンプリング及びダウンサンプリングを更に含むことができる。αの値はユーザが選択できる。例えばαは４に等しくなるように選択できるが、他の値であってもよい。ルートレイズドコサインフィルタは、帯域幅を制限するためにパルス整形を実行できる。異種時間反転システムにおいて、異なる種類の装置は異なる帯域幅を使用してもよいため、異なる種類の装置に対するルートレイズドコサインフィルタの帯域幅も異なる。
いくつかの実現例において、図３７を参照すると、チャネルプロービングは送受信機又はアクセスポイント（ＡＰ）５９０が端末装置５９２と通信する場合に実行される。チャネルプロービング段階において、端末装置５９２は、アップサンプリング部５９６を使用してインパルス信号５９４をαでアップサンプリングし、ルートレイズドコサイン（ＲＲＣ）フィルタＲＲＣB,fs５９８を使用してアップサンプル信号をフィルタリングする。この場合、Ｂは帯域幅であり、ｆｓ＝αＢはサンプルレートである。ビットレートを上げる目的の１つは、オーバサンプリングすることにより、データ送信段階において送信機と受信機との間にサンプル周波数オフセットが存在する場合に正確なサンプルを取得できることである。オーバサンプリングはデータ送信段階におけるステップの１つであってもよいため、チャネルを正確に推定するために、チャネルプロービング段階にオーバサンプリング及びダウンサンプリングを更に含むことができる。αの値はユーザが選択できる。例えばαは４に等しくなるように選択できるが、他の値であってもよい。ルートレイズドコサインフィルタは、帯域幅を制限するためにパルス整形を実行できる。異種時間反転システムにおいて、異なる種類の装置は異なる帯域幅を使用してもよいため、異なる種類の装置に対するルートレイズドコサインフィルタの帯域幅も異なる。
【0094】
フィルタリングされた信号は、デジタルアナログ変換器（ＤＡＣ）６００によりアナログベースバンド信号に変換される。ベースバンド信号は、電圧制御発振器（ＶＣＯ）６０２により高搬送周波数に変調され、バンドパスフィルタ６０４を通過し、電力増幅器６０６により増幅されて無線周波数（ＲＦ）信号が生成される。ＲＦ信号は、アンテナ６０８によりブロードキャストされる。ブロードキャスト信号は、プロファイルｈを有するマルチパスチャネル６１０を通過し、ここでノイズｎ'が信号に追加される。ブロードキャスト信号は、アクセスポイント５９０においてアンテナ６１２により取り込まれる。取り込まれた信号は、低ノイズ増幅器（ＬＮＡ）６１４により増幅され、バンドパスフィルタ６１６によりフィルタリングされ、電圧制御発振器６１８を使用してアナログベースバンド信号に変換される。ベースバンド信号は、アナログデジタル変換器（ＡＤＣ）６２０によりサンプルレートｆｓ＝αＢでサンプリングされ、デジタル信号が生成される。デジタル信号は、ルートレイズドコサインフィルタＲＲＣB,fs６２２によりフィルタリングされ、ダウンサンプリング部６２４により係数αでダウンサンプリングされる。インパルス信号が端末装置５９２により送信されたため、ダウンサンプリング信号はチャネル推定値
として扱われる。
【0095】
ルートレイズドコサインフィルタ５９８はデジタルフィルタであり、ルートレイズドコサインフィルタ６２２はアナログフィルタであり、双方は同一の帯域幅を有する。バンドパスフィルタ６０４は、端末装置５９２により使用される周波数帯域を定義する。従って、端末装置５９２が周波数ｆ１～周波数ｆ２の通過帯域を有するバンドパスフィルタ６０４を含む場合、端末装置５９２は周波数帯域ｆ１～ｆ２を使用するか又はそれに関連付けられると言う。図３７の例において、端末装置５９２のバンドパスフィルタ６０４の通過帯域はアクセスポイント５９０のバンドパスフィルタ６１６の通過帯域と同一である。
【0096】
アクセスポイント５９０は、本明細書中で説明する通信技術を使用して複数の他の装置と通信する基地局又は送受信機であってもよい。
【0097】
サンプリングレートｆｓ＝αＢの場合、離散チャネルインパルス応答を以下のように書くことができる。
式中、Ｔｓ＝１／（αＢ）である。完璧なチャネル推定値（チャネルプロービング段階においてノイズ及び干渉が無視される）を仮定すると、図３７の２つのルートレイズドコサインフィルタ５９８、６２２の間の同等のチャネルインパルス応答を以下のように書くことができる。
【0098】
図２６を参照すると、いくつかの実現例において、会場１４０におけるイベントを検出するために、送信機１４２はマルチパスチャネル１４６を介して受信機１４４に無線信号を送出する。送信機１４２はチャネルプロービング信号１４８を生成し、無線周波数変調１５０を実行して、関連するマルチパスチャネル１４６を通って受信機１４４に伝搬する無線周波数信号を生成する。受信機１４４は、送信された無線周波数信号を受信し、無線周波数復調１５２を実行して一意のチャネル応答（波形）１５４を再生し、これを記録する。送信機１４２に関連するチャネル応答１５４を受信すると、受信機１４４は、受信機１４４において受信したチャネル応答に基づいて、送信機１４２に対するシグニチャ波形１５６を計算する。文脈に依存して、用語「チャネルプロービング信号」は、無線周波数変調前のプローブ信号を示してもよく、あるいは無線周波数変調を適用後の信号を示してもよい。いくつかの実現例において、送信機１４２は図３７の装置５９２と同様の方法で実現可能であり、受信機１４４は図３７の装置５９０と同様の方法で実現可能である。
【0099】
受信機１４４は、チャネルプロービング信号１４８を使用してチャネル応答を推定できる。この目的を達成するために、受信機１４４は最初に、チャネル情報を知らない状態でマルチパスの影響を相殺してチャネルプロービング信号と同期する必要がある。一般的なチャネルプロービング信号は、プリアンブルと可能性のある確認ビットの２つの部分とを含む。受信機１４４は、マルチパスフェージングにより破損されないプリアンブルの基本特徴を抽出することにより、チャネルプロービング信号と同期する。その後、確認ビットを使用して、同期を更に検証する。同期が達成されると、受信機１４４はチャネルプロービング信号内の情報を使用して受信信号からチャネル応答を抽出する。
【0100】
【0101】
時間領域における時間反転パラメータ
時間領域におけるチャネル状態情報の表現がベクトルｈ＝［ｈ［１］，ｈ［２］，...，ｈ［Ｌ］］として表されると仮定する。式中、Ｌはチャネル状態情報のタップの数である。例えばｈ［１］、ｈ［２］、ｈ［Ｌ］等は、受信プロービング信号のデジタルサンプルであってもよい。サンプリングレートは、例えば１２５ＭＨｚであってもよい（他のサンプリングレートを使用することもできる）。時間領域における２つのチャネル状態情報サンプルｈ１及びｈ２の間の時空間共振の強度TR(ｈ1, ｈ2)を「時間反転共振強度（ＴＲＲＳ）」と呼び、これを以下のように定義する。
式中、「＊」は畳み込みを示し、ｌはタップインデックスを示し、ｇ２はｈ２のＴＲシグニチャである。すなわち、
ｇ２［ｋ］＝ｈ* ２［Ｌ－ｋ－１］， ｋ＝０，１，...，Ｌ－１ 式（２）
時間領域におけるチャネル状態情報の表現がベクトルｈ＝［ｈ［１］，ｈ［２］，...，ｈ［Ｌ］］として表されると仮定する。式中、Ｌはチャネル状態情報のタップの数である。例えばｈ［１］、ｈ［２］、ｈ［Ｌ］等は、受信プロービング信号のデジタルサンプルであってもよい。サンプリングレートは、例えば１２５ＭＨｚであってもよい（他のサンプリングレートを使用することもできる）。時間領域における２つのチャネル状態情報サンプルｈ１及びｈ２の間の時空間共振の強度TR(ｈ1, ｈ2)を「時間反転共振強度（ＴＲＲＳ）」と呼び、これを以下のように定義する。
ｇ２［ｋ］＝ｈ* ２［Ｌ－ｋ－１］， ｋ＝０，１，...，Ｌ－１ 式（２）
【0102】
２つのチャネル状態情報サンプルｈ１及びｈ２の間の時間反転共振強度が高いほど、ｈ１とｈ２との間の類似性が高いことを示し、２つのチャネル状態情報サンプルｈ１及びｈ２の間の時間反転共振強度が低いほど、ｈ１とｈ２との間の類似性が低いことを示す。環境が特定の状態にある場合の時間ｔ１において測定されたチャネル状態情報をｈ１が表し且つその後の時間ｔ２において測定されたチャネル状態情報をｈ２が表す場合、ｈ１とｈ２との間の時間反転共振強度は、時間ｔ２における環境と時間ｔ１における環境との類似度のインジケーションを提供する。
【0103】
チャネルインパルス応答（ＣＩＲ）は、時間領域において測定されるチャネル状態情報である。ｈ１とｈ*
２との内積は（ｈ１＊ｇ２）［Ｌ－１］と同等であるため、式（１）で定義された共振強度はｈ１とｈ２との間の相互相関係数の定義と同様である。しかし、式（１）の分子は畳み込み済み系列における最大絶対値である。本ステップの目的は、２つのチャネル状態情報推定値の間のあらゆる可能性のあるランダム位相歪みを補正することである。チャネル状態情報は通常、初期位相歪み及び線形位相歪みの影響を受ける。初期位相歪みは通常、搬送周波数オフセット（ＣＦＯ）及び共通位相オフセット（ＣＰＥ）により生じる。線形位相歪みは、シンボルタイミングオフセット（ＳＴＯ）及びサンプリング周波数オフセット（ＳＦＯ）により生じる場合がある。例えばチャネル状態情報の最初のいくつかのタップが損なわれるか又は異なる測定値に追加される場合があり、それにより２つのチャネル状態情報推定値の間に同期誤差が生じる。チャネル状態情報推定値における位相歪みを補正すると、時間反転共振強度は、精度を維持しつつ、チャネルインパルス応答間の一部又は実質的に全ての類似性を取り込むことができる。位相歪みを補正する方法を本明細書の後節で説明する。
【0104】
周波数領域における時間反転パラメータ
受信機側において、各チャネル状態サウンディングの後、それぞれに対するＬ×Ｍの未加工のチャネル状態情報行列を以下のように収集できる。
Ｈi ＝［ｈ(1), ｈ(2) ,... , ｈ(M)］, ∀i 式（３）
この場合、Ｍは送信機と受信機との間のリンクの数であり、Ｌはチャネル周波数応答（ＣＦＲ）における使用可能なサブキャリアの数である。本明細書において、各送信機は１つ以上のアンテナを有することができ、各受信機は１つ以上のアンテナを有することができる。送信機アンテナと受信機アンテナの各対は１つの通信リンクに対応する。従って、例えば送信機が３つのアンテナを有し且つ受信機が２つのアンテナを有する場合、このイベント検出器システムには３×２＝６個のリンクが存在すると言う。チャネル周波数応答は、周波数領域において測定されるチャネル状態情報である。
受信機側において、各チャネル状態サウンディングの後、それぞれに対するＬ×Ｍの未加工のチャネル状態情報行列を以下のように収集できる。
Ｈi ＝［ｈ(1), ｈ(2) ,... , ｈ(M)］, ∀i 式（３）
この場合、Ｍは送信機と受信機との間のリンクの数であり、Ｌはチャネル周波数応答（ＣＦＲ）における使用可能なサブキャリアの数である。本明細書において、各送信機は１つ以上のアンテナを有することができ、各受信機は１つ以上のアンテナを有することができる。送信機アンテナと受信機アンテナの各対は１つの通信リンクに対応する。従って、例えば送信機が３つのアンテナを有し且つ受信機が２つのアンテナを有する場合、このイベント検出器システムには３×２＝６個のリンクが存在すると言う。チャネル周波数応答は、周波数領域において測定されるチャネル状態情報である。
【0105】
２つのチャネル状態情報サンプルｈ１及びｈ２の間の周波数領域における時間反転時空間共振強度（ＴＲＲＳ） TR(ｈ1, ｈ2) を以下のように定義する。
式中、Ｌはチャネル状態情報ベクトルの長さであり、ｋはサブキャリアインデックスであり、ｇ２は以下のように取得されるｈ２の時間反転シグニチャである。
ｇ2［ｋ］＝ｈ* 2［ｋ］， ｋ=０, １, ...，Ｌ－１ 式（５）
ｇ2［ｋ］＝ｈ* 2［ｋ］， ｋ=０, １, ...，Ｌ－１ 式（５）
【0106】
式（４）の分子において、線形位相成分にわたる最大値探索動作が実行され、例えばＷｉＦｉである直交周波数分割多重（ＯＦＤＭ）に基づくシステムにおけるランダム線形位相オフセットが低減又は除去される。推定チャネル状態情報における線形位相オフセットは、発振器のサンプリングタイミングオフセット（ＳＴＯ）及びサンプリング周波数オフセット（ＳＦＯ）から生じる。
【0107】
２つの推定マルチパスプロファイル（時間領域におけるチャネルインパルス応答及び周波数領域におけるチャネル周波数応答）を比較する場合、それらは最初に時間反転空間にマッピングされ、そこで各々は１つの時間反転シグニチャとして表される。その場合、時間反転時空間共振強度は、マッピングされた時間反転空間における２つのマルチパスプロファイルの間の類似性を定量化するメトリックである。
【0108】
周波数領域におけるＴＲＲＳの統計的挙動
位相歪みを補正した後に推定チャネル周波数応答に位相オフセットが存在しない場合、多入力多出力（ＭＩＭＯ）システムにおける各リンクに対する式（４）の時間反転共振強度を以下のように計算できる。
式中、ｋはサブキャリアのインデックスであり、Ｌはサブキャリアの総数である。チャネルが定常性を有するという仮定に従うと、チャネル周波数応答（ＣＦＲ）ｈ０及びｈ１が同一のロケーション又は同一の状態から生じる場合、ｈ１を以下のようにモデル化できる。
ｈ１＝ｈ０＋ｎ 式（７）
式中、ｎはガウスノイズのベクトルであり、
Ｌはチャネル周波数応答の長さである。
位相歪みを補正した後に推定チャネル周波数応答に位相オフセットが存在しない場合、多入力多出力（ＭＩＭＯ）システムにおける各リンクに対する式（４）の時間反転共振強度を以下のように計算できる。
ｈ１＝ｈ０＋ｎ 式（７）
式中、ｎはガウスノイズのベクトルであり、
【0109】
【0110】
【0111】
【0112】
【0113】
【0114】
【0115】
時間反転イベント検出器：ベースライン法
段階１．オフライントレーニング
時間反転イベント検出器（ＴＲＥＤ）は、一意の屋内マルチパスプロファイル及び時間反転技術を利用して、例えば屋内イベントであるイベントを区別して検出する。オフライントレーニング段階において、あらゆる対象イベントのマルチパスプロファイルが収集され、時間反転空間における対応する時間反転シグニチャがデータベースに格納される。しかし、ノイズ及びチャネルフェージングにより、特定の状態からのチャネル状態情報は時間と共に若干変化する場合がある。この種の変動に対処するために、状態毎にいくつかの瞬間チャネル状態情報サンプルを収集してトレーニングセットを構築する。
段階１．オフライントレーニング
時間反転イベント検出器（ＴＲＥＤ）は、一意の屋内マルチパスプロファイル及び時間反転技術を利用して、例えば屋内イベントであるイベントを区別して検出する。オフライントレーニング段階において、あらゆる対象イベントのマルチパスプロファイルが収集され、時間反転空間における対応する時間反転シグニチャがデータベースに格納される。しかし、ノイズ及びチャネルフェージングにより、特定の状態からのチャネル状態情報は時間と共に若干変化する場合がある。この種の変動に対処するために、状態毎にいくつかの瞬間チャネル状態情報サンプルを収集してトレーニングセットを構築する。
【0116】
基本的に、イベント検出器は分類を使用する。分類は、カテゴリの所属が既知である観察値を含むデータのトレーニングセットに基づいて、新規の観察値がどのカテゴリ／クラスセットに属するかを識別することである。特に、具体的な一実現例において、分類を実現する処理又はアルゴリズムは分類器として知られる。その一例は、所定の電子メールを迷惑メールクラス又は非迷惑メールクラスに割り当てること、特定の症状の存在に基づいて所定の患者の観察された特徴が示す通りに患者に診断を割り当てること等である。多くの場合、個々の観察値は、特徴としても知られる定量化可能なプロパティのセットに解析される。いくつかの例において、特徴は分類可能で、順序を示すことができ、整数値又は実数値を有することができる。他の例において、分類器は、相似関数又は距離関数を使用して観察値と以前の観察値とを比較することにより機能する。更に別の例において、分類器は、入力データをカテゴリにマッピングする分類アルゴリズムにより実現される数学関数を更に参照する場合がある。分類器による決定が行われ且つそれが正確な決定であると確認されると、その後に収集される観察値は以前の観察値に追加され、それを用いて分類器／特徴を再トレーニングできる。
【0117】
開示されるイベント検出器において、トレーニングセットは、既知である所定の各イベントにおいて収集された受信プロービング信号からの推定ＣＳＩを含む。ベースライン法において、イベント検出器は、類似性スコアとして見なすことができるＴＲＲＳにより、監視段階で取得された推定ＣＳＩとトレーニング段階における観察値とを比較する。イベント検出に対する決定が行われ且つそれが正確な決定であると確認されると、監視段階において収集されたＣＳＩをデータベースに追加でき、それを用いて特徴を再トレーニングできる。
【0118】
特に、Ｓを状態セットとする各状態Ｓｉ∈Ｓに対して、対応する時間領域トレーニングチャネル状態情報、すなわちチャネルインパルス応答が推定され、以下のようにＨｉを形成する。
この場合、ベクトル変数の上付き文字＊は共役演算子を表す。Ｌはチャネル状態情報ベクトルの長さを示し、ｋはチャネル状態情報における時間領域タップのインデックスである。その場合、トレーニングデータベースｇはＧｉの集合である。
【0119】
段階２．オンライン監視
トレーニングデータベースｇを構築すると、時間反転イベント検出器はイベント検出を行える状態になる。監視段階の間、受信機は現在の推定チャネル状態情報とｇ内の時間反転シグニチャとをマッチングさせ続け、最も強い時間反転時空間共振を生成するものを探す。未知の監視中チャネル状態情報Ｈ~と状態Ｓｉ中のサンプルとの間の時間反転共振強度は、以下のように定義される。
トレーニングデータベースｇを構築すると、時間反転イベント検出器はイベント検出を行える状態になる。監視段階の間、受信機は現在の推定チャネル状態情報とｇ内の時間反転シグニチャとをマッチングさせ続け、最も強い時間反転時空間共振を生成するものを探す。未知の監視中チャネル状態情報Ｈ~と状態Ｓｉ中のサンプルとの間の時間反転共振強度は、以下のように定義される。
【0120】
【0121】
【0122】
ベースライン法の実験結果
ベースライン法を使用して実験が行われた。図２Ａの間取り図２２０に詳細な設定を示す。図中、異なるドットマークは送信機（ＴＸ）及び受信機（ＲＸ）の異なるロケーションを表す。実験において、Ｄ１～Ｄ８と示される複数の木製扉の開閉状態を検出する。丸いドット（例えば、２２２ａ、２２２ｂ、２２２ｃ、２２２ｄ）で示される送信機の各ロケーションは互いに１ｍ離間し、受信機の候補ロケーションは星（例えば、２２４ａ及び２２４ｂ）で示される。ＴＸ－ＲＸロケーションは、見通し内（ＬＯＳ）送信及び見通し外（ＮＬＯＳ）送信の双方を含む。
ベースライン法を使用して実験が行われた。図２Ａの間取り図２２０に詳細な設定を示す。図中、異なるドットマークは送信機（ＴＸ）及び受信機（ＲＸ）の異なるロケーションを表す。実験において、Ｄ１～Ｄ８と示される複数の木製扉の開閉状態を検出する。丸いドット（例えば、２２２ａ、２２２ｂ、２２２ｃ、２２２ｄ）で示される送信機の各ロケーションは互いに１ｍ離間し、受信機の候補ロケーションは星（例えば、２２４ａ及び２２４ｂ）で示される。ＴＸ－ＲＸロケーションは、見通し内（ＬＯＳ）送信及び見通し外（ＮＬＯＳ）送信の双方を含む。
【0123】
最初に、複数の屋内イベントを検出する際の時間反転イベント検出器の性能を調べる。更に、受信信号強度インジケータ（ＲＳＳＩ）に基づく屋内検出方法と時間反転イベント検出器との性能比較を調べる。受信信号強度インジケータは、受信機において受信された信号エネルギーを表すスカラである。無線信号が空中を移動する場合、異なる物体に遭遇し、送信信号の複数の変化したコピーが生成される。更に、送信中、異なる送信信号は屋内環境に応答して異なる減衰が生じる。従って、屋内環境の各スナップショット、すなわち各屋内イベントは異なる受信信号強度インジケータ値に対応してもよく、これは異なる屋内イベントを認識するために使用できる。トレーニング段階において、システムは異なる注目屋内イベントに対する受信信号強度インジケータを収集してデータベースに格納する。受信信号強度インジケータを使用して屋内イベントを分類する分類器は最近傍法に基づき、監視中の環境からの受信信号強度インジケータはトレーニング段階で収集された全ての屋内イベントの受信信号強度インジケータと比較される。屋内イベント検出に対する決定は、監視中の環境からの受信信号強度インジケータと最もマッチングする格納済み受信信号強度インジケータを見つけることにより行われる。
【0124】
受信信号強度インジケータに基づく方法を使用する実験において、受信機はロケーションＢ（２２４ａ）又はロケーションＣ（２２４ｂ）のいずれかに配置され、送信機のロケーションはドット（２２２ａ～２２２ｄ）により表される。それらは最も左側のロケーションから１ｍずつ離間し、それぞれを「軸１」～「軸４」と呼ぶ。全部で２つの受信機のロケーションと４つの送信機のロケーションがあり、すなわち８個のＴＸ－ＲＸ対がある。時間反転イベント検出器の目的は、Ｄ１～Ｄ８のうち、他の全ての扉が開いている中で閉まった木製扉を検出することである。
【0125】
時間反転イベント検出器及び受信信号強度インジケータに基づく方法の全体的な誤警報率及び検出率を表２及び表３に示す。
【0126】
見通し内の場合及び見通し外の場合の双方において、時間反転イベント検出器は受信信号強度インジケータに基づく方法より高い検出率及び低い誤警報率を達成し、性能が優れている。時間反転イベント検出器は、見通し外の場合に９６．９２％以上の検出率及び３．０８％以下の誤警報を得ることができ、見通し内の場合に９７．８９％以上の検出率及び２．１１％以下の誤警報を得ることができる。更に、受信機と送信機との間の距離が増加するほど、双方の方法の精度は向上する。
【0127】
図２Ｂを参照すると、時間反転イベント検出器に対する人間の移動の影響を調査するために、図中に網掛け部分２３０として示す領域を継続的に前後に移動する人間が０人、１人及び２人いる状態で実験を行った。また、送信機は丸印２３２で表すロケーションに配置され、受信機は六角形印２３４で表すロケーションに配置され、「Ｄ１」及び「Ｄ２」で示す隣接する２つの扉の状態を表４に示すように検出した。
【0128】
人間の移動により導入される干渉はマルチパス伝搬を変化させ、時間反転イベント検出器の監視段階における時間反転共振強度に変動を生じさせる。人間が会場内を動き回る場合、導入される干渉は変化し続け、各干渉の持続時間は短い。時間反転共振強度におけるバースト変動の影響に対処するために、「スライディングウィンドウ」と組み合わせた「多数決」法を採用して、検出結果を時間にわたり平滑化する。先行するＫ－１個の出力Ｓ*
k（ｋ＝ｔ－Ｋ＋１，...，ｔ－１）及び現在の結果S_t^*を有すると仮定した場合、全てのＳ*
k（ｋ＝ｔ－Ｋ＋１，...，ｔ－１）にわたり投票することにより、時間スタンプｔに対する決定が行われる。Ｋは平滑化用スライディングウィンドウのサイズを示す。Ｋ個の出力Ｓ*
kのうち最も多くの票を有するイベントが時間スタンプｔに対するイベントとして選択される。
【0129】
例えばＫ＝１１であり、時間反転イベント検出器の１１個の出力Ｓ*
k（ｋ＝ｔ－１０，．．．，ｔ）のセットが｛Ｅ１，Ｅ１，Ｅ２，Ｅ１，Ｅ１，Ｅ２，Ｅ２，Ｅ２，Ｅ２，Ｅ１，Ｅ２｝であるとする。本例において、５個の出力がイベントＥ１に対応し、６個の出力がイベントＥ２に対応し、イベントＥ２が時間スタンプｔに対するイベントとして選択される。時間スタンプｔ＋３において、時間反転イベント検出器の１１個の出力Ｓ*
k（ｋ＝ｔ－７，．．．，ｔ＋３）のセットが｛Ｅ１，Ｅ１，Ｅ２，Ｅ２，Ｅ２，Ｅ２，Ｅ１，Ｅ２，Ｅ１，Ｅ１，Ｅ１｝であるとする。本例において、６個の出力がイベントＥ１に対応し、５個の出力がイベントＥ２に対応するため、イベントＥ１が時間スタンプｔ＋３に対するイベントとして選択される。
【0130】
表５において、人間の移動（ＨＭ）がない場合と、１人の人間及び２人の人間により意図的に実行された継続的な人間の移動がある場合との平滑化アルゴリズム使用時又は未使用時の時間反転イベント検出器の平均精度を比較する。この場合、スライディングウィンドウの長さはＫ＝２０であり、これは０．１秒の持続時間を表す。第１に、時間反転イベント検出器の精度は、屋内イベント、送信機及び受信機のロケーションの付近で継続的な移動を実行する人間の数が増加するにつれて減少する。更に、採用された平滑化アルゴリズムは人間の動作に対する時間反転イベント検出器のロバスト性を向上し、平滑化を行わない場合と比較して精度を７％～９％向上する。その一方で、実験において、最も脆弱な状態は全ての扉が開いている状態「００」（表４を参照）であることが更に分かった。この状態で人間の移動が存在する場合、時間反転イベント検出器が「扉が閉じた」という誤報を生成する確率は他の状態より高い。これは、人間が扉のロケーションの付近を移動する場合、人間の身体が扉のロケーションにおける重度の閉鎖と見なされ、木製扉が閉じた状態と類似するため、マルチパスチャネル状態情報の変化が「扉が閉じた」と同様であるためである。
【0131】
時間反転イベント検出器：統計的方法
いくつかの実現例において、イベント検出器は、「開いている扉は無い」、「玄関扉が開いている」等である屋内イベントの発生を検出するように設計される。ここでは、統計値に基づくイベント検出器の詳細を開示する。統計値に基づくイベント監視システム（又は保護システム）は、オフライントレーニング段階及びオンライン学習段階である２つの動作段階を有する。トレーニングセットは、所定の既知の各イベントにおいて収集された受信プロービング信号からの推定ＣＳＩを含み、それらから各イベントに対する時間反転共振強度の関数における統計的分布を計算する。オンライン監視段階では、未知のイベントにおいて収集されたＣＳＩを使用する確率特徴が抽出され、未知のイベントは既知のイベントの１つに分類される。イベント検出の決定が行われ且つそれが正確な決定であると確認されると、監視段階において収集されたチャネル状態情報をトレーニングセットに追加でき、それを用いて特徴を再トレーニングできる。
いくつかの実現例において、イベント検出器は、「開いている扉は無い」、「玄関扉が開いている」等である屋内イベントの発生を検出するように設計される。ここでは、統計値に基づくイベント検出器の詳細を開示する。統計値に基づくイベント監視システム（又は保護システム）は、オフライントレーニング段階及びオンライン学習段階である２つの動作段階を有する。トレーニングセットは、所定の既知の各イベントにおいて収集された受信プロービング信号からの推定ＣＳＩを含み、それらから各イベントに対する時間反転共振強度の関数における統計的分布を計算する。オンライン監視段階では、未知のイベントにおいて収集されたＣＳＩを使用する確率特徴が抽出され、未知のイベントは既知のイベントの１つに分類される。イベント検出の決定が行われ且つそれが正確な決定であると確認されると、監視段階において収集されたチャネル状態情報をトレーニングセットに追加でき、それを用いて特徴を再トレーニングできる。
【0132】
段階１．オフライントレーニング
図１９を参照すると、オフライントレーニング段階２４０において、イベント監視システムは、監視される屋内イベント毎に同一のイベントプロファイルから収集されたチャネル周波数応答の間の時間反転共振強度の対数正規分布モデルを構築することを目的とする。
図１９を参照すると、オフライントレーニング段階２４０において、イベント監視システムは、監視される屋内イベント毎に同一のイベントプロファイルから収集されたチャネル周波数応答の間の時間反転共振強度の対数正規分布モデルを構築することを目的とする。
【0133】
特に、正常イベント及び異常イベントを含む屋内イベントのセットをＳとする各屋内イベントＳｉ∈Ｓに対して、対応する周波数領域チャネル周波数応答が受信機側で取得され、以下のように推定される。
Ｈｉ＝［ｈi (1), ｈi (2),..., ｈi (M)］, ｉ＝１,２,...,Ｎ 式（３２）
式中、ＮはＳのサイズであり、Ｍは送信機と受信機との間のリンクの数であり、Ｈｉの次元はＬ×Ｍである。Ｍは、送信機アンテナの数と受信機アンテナの数の積である。ここで、用語「正常イベント」は、ユーザ／環境が正常であると考え、警報を生成する必要がないと考えるイベントを示す。用語「異常イベント」は、ユーザ／環境が正常でない又は異常であると考え、警報を生成する必要があると考えるイベントを示す。
Ｈｉ＝［ｈi (1), ｈi (2),..., ｈi (M)］, ｉ＝１,２,...,Ｎ 式（３２）
式中、ＮはＳのサイズであり、Ｍは送信機と受信機との間のリンクの数であり、Ｈｉの次元はＬ×Ｍである。Ｍは、送信機アンテナの数と受信機アンテナの数の積である。ここで、用語「正常イベント」は、ユーザ／環境が正常であると考え、警報を生成する必要がないと考えるイベントを示す。用語「異常イベント」は、ユーザ／環境が正常でない又は異常であると考え、警報を生成する必要があると考えるイベントを示す。
【0134】
いくつかの実現例において、クラス内時間反転共振強度の統計値を学習するために、Ｈｉの少なくとも３００個の実現例が収集される。屋内イベントＳｉに対するチャネル状態情報の十分なサンプルを収集すると、以下のステップを使用して統計値が推定される。
【0135】
・位相オフセットの較正：チャネル状態情報の各実現例にランダム線形位相オフセットがあるため、何らかの位相較正法を適用して全ての実現例を整列させ、その後、式（６）を使用して時間反転共振強度を計算できる。
【0136】
・ＣＳＩ代表生成器：本ステップにおいて、リンク毎にトレーニングセット内の全ての屋内イベントに対して、同一の屋内イベントからの他の全てのチャネル周波数応答に最も類似する代表チャネル状態情報が見つけられる。最初に、屋内イベントＳｉに対して収集された全ての実現例の間のリンクｍに対する対の時間反転共振強度が計算される。時間反転共振強度は異なるチャネル周波数応答の間の類似性を定量化するのに使用できるため、自身とイベントＳｉにおけるリンクｍに対する他の全てのチャネル周波数応答との間の時間反転共振強度が事前に定義された類似性閾値を上回る数が最も多いチャネル周波数応答が、イベントＳｉにおけるリンクｍに対する代表チャネル状態情報として選択される。全てのリンク及び全ての屋内イベントに対して、これを繰り返す。その結果、全ての代表チャネル状態情報は以下のような集合を形成する。
式中、ｈｉ
(l)（ｎｌ）は、ｎｌ番目の実現例からのリンクｌに対するイベントＳｉのチャネル周波数応答を示す。ｎｌ（ｌ＝１，２，．．．，Ｍ）は実現例インデックスであり、これは全てのｌ及びＨｒｅｐ，ｉにおけるｉに対して異なってもよい。
【0137】
・対数正規パラメータの推定：代表チャネル状態情報が多入力多出力システムにおける全てのイベント及び全てのリンクに対して見つけられると、対数正規分布パラメータを推定できる。リンクｍ及びイベントＳｉに対して、代表チャネル状態情報Ｈ(m)
ｒｅｐ，ｉと他の全ての実現例ｈｉ
(m)（ｌ）との間の時間反転共振強度が式（６）を使用して計算され、以下のように示される。
【0138】
【0139】
トレーニングの最後に、全てのトレーニング済みイベントに対する代表チャネル状態情報Ｈｒｅｐの集合及び対数正規分布パラメータＱｒｅｐの集合を含むことにより、トレーニングデータベースが構築される。全てのトレーニング済みイベントは、警報が生成されない正常イベントＳnormal及び異常イベントが検出された場合に警報がユーザに報告される必要のある異常イベントＳabnormalである２つのグループに分類可能である。
【0140】
段階II．オンライン監視
監視段階において、図１９及び図２０に示すように、受信機は未知のイベントからチャネル状態情報をＨtest＝［ｈ(1) test, ｈ(2) test,..., ｈ(M) test］として継続的に収集する。統計値に基づくイベント検出器は、統計メトリックを計算して評価することにより、未知の監視されたチャネル状態情報に最も類似するトレーニングデータベース内のイベントを見つける。本明細書中、用語「監視されたチャネル状態情報」又は「検査されたチャネル状態情報」は、監視段階において収集されるチャネル状態情報を示す。詳細なステップを以下に列挙する。
監視段階において、図１９及び図２０に示すように、受信機は未知のイベントからチャネル状態情報をＨtest＝［ｈ(1) test, ｈ(2) test,..., ｈ(M) test］として継続的に収集する。統計値に基づくイベント検出器は、統計メトリックを計算して評価することにより、未知の監視されたチャネル状態情報に最も類似するトレーニングデータベース内のイベントを見つける。本明細書中、用語「監視されたチャネル状態情報」又は「検査されたチャネル状態情報」は、監視段階において収集されるチャネル状態情報を示す。詳細なステップを以下に列挙する。
【0141】
・ＴＲＲＳの計算：取得されたチャネル状態情報測定値Ｈtestが位相オフセットにより破損されるため、トレーニング段階で使用したのと同一の位相較正アルゴリズムをチャネル状態情報Ｈtestに最初に適用して、ランダム位相オフセットを除去する必要がある。その後、トレーニング済み屋内イベント毎に、代表チャネル状態情報と監視測定値との間の時間反転共振強度を式（６）により全てのリンクに対して計算する。これは以下のように示される。
【0142】
・統計メトリックの計算：時間反転共振強度が取得されると、Ｈtestと全てのトレーニング済みイベントＳｉとの間の統計メトリックが以下のように計算される。
式中、Ｆ(μ,σ)（ｘ）は、パラメータ（μ，σ）及び変数ｘを有する対数正規分布の累積分布関数（ＣＤＦ）である。関数Ｆ(μ,σ)（ｘ）が独立変数ｘ及びパラメータ（μ，σ）を有する累積分布関数であるため、これはｘと共に増加する。より大きいｘはより高いＦ(μ,σ)（ｘ）を示し、より高いＦ(μ,σ)（ｘ）はより大きいｘを示す。従って、より大きいＷi,testはより小さいＴＲ(m)
i,testセットを示し、これは、監視されたチャネル状態情報とトレーニング済みイベントＳｉのチャネル状態情報との間の類似性が低いため、未知のイベントの監視されたチャネル状態情報がトレーニング済みイベントＳｉに属する可能性が低いことを意味する。そのため、当該統計メトリックは、未知のイベントからの監視されたチャネル状態情報が何らかのトレーニング済みイベントＳｉに属する可能性の低さを反映する。
【0143】
・決定：提案される統計値に基づくイベント検出器は、２段階状態チェッカを用いて、監視中のチャネル状態情報測定値に対する決定を行う。詳細を以下に示す。
（ａ）段階１－正常状態チェッカ：環境が正常である（すなわち、正常イベントのうちの１つのみが発生する）か又は異常である（すなわち、現在の環境状態がトレーニング済み正常イベントのいずれにも属さない）かを検出するために、メトリック値Ｗnormal,testを使用して、事前に定義された閾値γnormalと比較する。
Ｗnormal,test＜γnormalの場合、監視中のチャネル状態情報測定値に対するシステムの決定はＤtest＝ｎｏｒｍａｌであり、決定手順は終了する。Ｗnormal,test＜γnormalでない場合、決定手順は次の段階である異常状態チェッカに進む。
（ｂ）段階２－異常状態チェッカ：本ステップでは、Ｈtestが収集される時間の間に何らのイベントが屋内環境において発生するため、どのトレーニング済み異常イベントが発生しているのかを判定することを目的とする。従って、メトリック値Ｗabnormal,testを使用して、事前に定義された閾値γabnormalと比較する。
その後、監視時間に対するシステムの決定が式（４３）のＷabnormal,testに基づいて以下のように行われる。
Ｄtest＝ｎは、トレーニング段階でトレーニングされていない何らかの異常イベントが発生していることを示す。
（ａ）段階１－正常状態チェッカ：環境が正常である（すなわち、正常イベントのうちの１つのみが発生する）か又は異常である（すなわち、現在の環境状態がトレーニング済み正常イベントのいずれにも属さない）かを検出するために、メトリック値Ｗnormal,testを使用して、事前に定義された閾値γnormalと比較する。
（ｂ）段階２－異常状態チェッカ：本ステップでは、Ｈtestが収集される時間の間に何らのイベントが屋内環境において発生するため、どのトレーニング済み異常イベントが発生しているのかを判定することを目的とする。従って、メトリック値Ｗabnormal,testを使用して、事前に定義された閾値γabnormalと比較する。
【0144】
統計的方法の実験結果
図３を参照すると、実験は一軒の家屋において行われ、間取り図２５０が図示される。間取り図２５０において、間取り図２５０中の矢印で示される扉及び窓の状態が統計値に基づくイベント検出器により監視される。試験家屋において、受信器のロケーション２５２が与えられ、２つの送信機のロケーション２５４ａ、２５４ｂが試験されて評価される。一方の送信機のロケーション２５４ａは玄関扉の領域に存在し、他方の送信機のロケーション２５４ｂは台所に存在する。
図３を参照すると、実験は一軒の家屋において行われ、間取り図２５０が図示される。間取り図２５０において、間取り図２５０中の矢印で示される扉及び窓の状態が統計値に基づくイベント検出器により監視される。試験家屋において、受信器のロケーション２５２が与えられ、２つの送信機のロケーション２５４ａ、２５４ｂが試験されて評価される。一方の送信機のロケーション２５４ａは玄関扉の領域に存在し、他方の送信機のロケーション２５４ｂは台所に存在する。
【0145】
受信機動作特徴（ＲＯＣ）性能が送信機の双方のロケーション２５４ａ、２５４ｂにおいて評価される。構成Ｉでは、送信機はロケーション２５４ａに配置される。構成IIでは、送信機はロケーション２５４ｂに配置される。図４Ａ及び図４Ｂは、検出率と誤警報率との間のトレードオフを表す受信機動作特徴の曲線プロットを示すグラフ２６０及びグラフ２７０である。検出率は、正確に分類されたクラスｉからの監視サンプルの百分率であり、誤警報率は、クラスｉに含まれないがクラスｉとして識別された監視サンプルの百分率である。
【0146】
送信機が玄関扉領域２５４ａに存在し且つ受信機がＡｌｉｃｅの部屋２５２の内側に配置される第１の構成において、全てのイベントがトレーニングデータベースに含まれる場合、受信機動作特徴性能は完全ではなく、図４Ｂのグラフ２７０に示すように、イベントｅ７の勉強部屋の扉２５６をｅ１と区別できない。これは、イベントｅ７が送信機及び受信機から遠すぎ、イベントｅ７により導入されたマルチパスの変化が小さすぎて、イベントｅ７と全ての扉が閉じているイベントｅ１とを分離できないためである。遠くのイベントをトレーニングデータベースから除外すると、図４Ａのグラフ２６０に示すように、受信機動作特徴性能は完全である。一般に、適切な監視性能を有するためには、イベントに対して監視されている物体（例えば、扉又は窓）は送信機又は受信機の少なくとも一方に対する見通し内パスを有するか又は送信機と受信機との間のリンクの近くに位置する必要がある。上記の基準を「適切な監視性能のための基準」と呼ぶ。構成Ｉにおいて、ｅ６及びｅ７以外の全てのイベントは適切な監視性能のための基準を満たす。
【0147】
構成Ｉにおける結果と同様に、構成IIの場合、ｅ２、ｅ８及びｅ９以外の全ての対象イベントは適切な監視性能のための基準を満たす。適切な監視性能のための基準を満たす対象イベントに対して、システムは図５Ａのグラフ２８０に示すように、誤警報率がほぼ０である完璧な検出精度を有する。図５Ｂを参照すると、グラフ２９０は、適切な監視性能のための基準を満たさないイベントに対して、システムが高い誤警報率を有することを示す。構成Ｉ及び構成IIの双方において、適切な監視性能のための基準を満たさないイベントは、「全ての扉が閉じている」と分類される。
【0148】
評価：動作
ここでは、構成Ｉのイベントリストの全ての対象イベントが意図的に実行され且つシステムが各イベントの間にチャネル状態情報を継続的に収集する監視環境がシミュレーションされる。３つの例を図６Ａ～図６Ｃに示す。図６Ａは、全ての扉及び窓が閉じている状況に対する決定を示すグラフ３００である。図６Ｂは、Ｂｏｂの部屋の扉が部屋の内側から開かれる状況に対する決定を示すグラフ３１０である。図６Ｃは、Ａｌｉｃｅの部屋の窓が外側から開かれる状況に対する決定を示すグラフ３２０である。グラフ３００、３１０、３２０において、ｘ軸は決定インデックス又は時間インデックスを表し、決定は一定期間にわたり周期的に行われる。ｙ軸は表６に対応する状態インデックスを表す。これらの例において、システムはイベントｅ１からイベントｅ４への移行（図６Ｂに示すように）及びイベントｅ１からイベントｅ９への移行（図６Ｃに示すように）の検出に成功している。
ここでは、構成Ｉのイベントリストの全ての対象イベントが意図的に実行され且つシステムが各イベントの間にチャネル状態情報を継続的に収集する監視環境がシミュレーションされる。３つの例を図６Ａ～図６Ｃに示す。図６Ａは、全ての扉及び窓が閉じている状況に対する決定を示すグラフ３００である。図６Ｂは、Ｂｏｂの部屋の扉が部屋の内側から開かれる状況に対する決定を示すグラフ３１０である。図６Ｃは、Ａｌｉｃｅの部屋の窓が外側から開かれる状況に対する決定を示すグラフ３２０である。グラフ３００、３１０、３２０において、ｘ軸は決定インデックス又は時間インデックスを表し、決定は一定期間にわたり周期的に行われる。ｙ軸は表６に対応する状態インデックスを表す。これらの例において、システムはイベントｅ１からイベントｅ４への移行（図６Ｂに示すように）及びイベントｅ１からイベントｅ９への移行（図６Ｃに示すように）の検出に成功している。
【0149】
家屋の外側のアクティビティに対する統計値に基づくイベント検出器のロバスト性をシミュレーションする。図７Ａを参照すると、グラフ３３０は玄関扉の外側を歩いている人間がいる時の決定結果を示す。出力された決定は、「全ての扉が閉じている」として全て正確に判定されている。図７Ｂを参照すると、グラフ３４０は裏口の扉の外側を歩いている人間がいる時の決定結果を示す。出力された決定は、「全ての扉が閉じている」として全て正確に判定されている。図７Ｃを参照すると、グラフ３５０は自動車が家屋の外側を走行している時の決定結果を示す。出力された決定は、「全ての扉が閉じている」として全て正確に判定されている。グラフ３３０、３４０及び３５０に示す決定結果は、システムが外側のダイナミクスに対して信頼性及び安定性を有することを示す。
【0150】
時間反転イベント検出器：機械学習法
屋内イベントの発生を検出するという、統計値に基づく方法と同一の目的を有する機械学習に基づくイベント検出器を以下に説明する。機械学習に基づく方法は、イベントを特徴空間にマッピングし、特徴空間における異なる屋内イベント間の距離を最大にしようとする。本方法は、計算の複雑性が高くなるが（統計値に基づく方法と比較して）、より高い決定性能を達成できる。機械学習に基づくイベント検出器は、図２１に示すように、オフライントレーニング段階４７０及びオンライン監視段階４７２で動作できる。
屋内イベントの発生を検出するという、統計値に基づく方法と同一の目的を有する機械学習に基づくイベント検出器を以下に説明する。機械学習に基づく方法は、イベントを特徴空間にマッピングし、特徴空間における異なる屋内イベント間の距離を最大にしようとする。本方法は、計算の複雑性が高くなるが（統計値に基づく方法と比較して）、より高い決定性能を達成できる。機械学習に基づくイベント検出器は、図２１に示すように、オフライントレーニング段階４７０及びオンライン監視段階４７２で動作できる。
【0151】
段階１．オフライントレーニング
統計値に基づく方法と同様に、トレーニング段階４７０において、監視される全ての対象イベントの屋内マルチパスプロファイルが収集され、サポートベクトルマシン（ＳＶＭ）及び特徴次元縮小法、すなわち主成分分析（ＰＣＡ）である分類器をトレーニングするために使用される。
統計値に基づく方法と同様に、トレーニング段階４７０において、監視される全ての対象イベントの屋内マルチパスプロファイルが収集され、サポートベクトルマシン（ＳＶＭ）及び特徴次元縮小法、すなわち主成分分析（ＰＣＡ）である分類器をトレーニングするために使用される。
【0152】
データの取得：最初に、Ｓが正常イベント及び異常イベントを含む屋内イベントのセットである各屋内イベントＳｉ∈Ｓに対して、対応する周波数領域チャネル周波数応答が取得され、受信機側において式（３２）と同様の形式で推定される。
【0153】
・位相オフセットの較正：チャネル状態情報の各実現例にランダム線形位相オフセットがあるため、位相較正法を適用して全ての実現例を個別に整列させた後、各チャネル周波数応答のエネルギーを正規化する。
【0154】
・初期特徴の生成：各イベントＳｉに対して、チャネル状態情報プロファイルＨｉが初期特徴として処理される。各リンクチャネル状態情報ｈi
(m)は複素値ベクトルであるため、機械学習用途で使用する前に処理する必要がある。複素値正規化チャネル状態情報の実数部及び虚数部を抽出して連結し、以下のような新規の特徴ベクトルを形成する。
【0155】
・特徴次元縮小：主成分分析において、直交変換を適用して、相関変数である可能性のあるものを含む元の特徴を各次元が無相関である新規の特徴空間に変換する。主成分分析は、全てのイベントに対して先行ステップで取得された新規に生成された特徴ベクトルを使用して学習される。いくつかの実現例において、大きい順に２００個の変換後成分が新規特徴として選択され、特徴の次元が縮小される。同一の直交変換が監視段階で使用される。尚、上述した２００という数は一例にすぎず、用途に応じて、大きい順に例えば１００個又は３００個である別の数の変換後成分を新規特徴として選択できる。
【0156】
・サポートベクトルマシンのトレーニング：より小さい次元の新規特徴が生成されると、それらは線形サポートベクトルマシンに提供され、異なるクラス（イベント）間のマージンが最大になるように重みベクトルがトレーニングされる。
【0157】
トレーニング段階の最後に、主成分分析に対する直交変換及び線形サポートベクトルマシン分類器が学習される。
【0158】
段階２．オンライン監視
未知のイベントから収集された各サンプルチャネル状態情報に対して、トレーニング段階で使用されたのと同一の前処理を最初に実行して、位相オフセット及び初期位相歪みを低減又は除去し、その後、チャネル状態情報の実数部及び虚数部を連結することにより特徴ベクトルが生成される。その後、学習された主成分分析の直交変換を特徴ベクトルに適用して、次元数を減少する。新規の特徴ベクトルは、トレーニング済み線形サポートベクトルマシンを使用して分類され且つラベル付けされる。イベント検出に関する決定が行われ且つそれが正確な決定であると確認されると、監視段階で収集されたチャネル状態情報をトレーニングセットに追加でき、それを用いてサポートベクトルマシンを再トレーニングできる。
未知のイベントから収集された各サンプルチャネル状態情報に対して、トレーニング段階で使用されたのと同一の前処理を最初に実行して、位相オフセット及び初期位相歪みを低減又は除去し、その後、チャネル状態情報の実数部及び虚数部を連結することにより特徴ベクトルが生成される。その後、学習された主成分分析の直交変換を特徴ベクトルに適用して、次元数を減少する。新規の特徴ベクトルは、トレーニング済み線形サポートベクトルマシンを使用して分類され且つラベル付けされる。イベント検出に関する決定が行われ且つそれが正確な決定であると確認されると、監視段階で収集されたチャネル状態情報をトレーニングセットに追加でき、それを用いてサポートベクトルマシンを再トレーニングできる。
【0159】
機械学習に基づく方法の実験結果
本節では、前節と同一のデータを使用する受信機動作特性性能及び運用テストを通じて、機械学習に基づくイベント検出器の性能を評価する。送信機が玄関扉領域に配置される構成Ｉを調べる。主成分分析の次元縮小後の分類特徴を形成するために、２０個の実数値主成分が選択される。
本節では、前節と同一のデータを使用する受信機動作特性性能及び運用テストを通じて、機械学習に基づくイベント検出器の性能を評価する。送信機が玄関扉領域に配置される構成Ｉを調べる。主成分分析の次元縮小後の分類特徴を形成するために、２０個の実数値主成分が選択される。
【0160】
適切な監視性能の基準を満たす対象イベントに対して、機械学習に基づくイベント検出器は、図８Ａのグラフ３６０に示すように、誤警報率がほぼ０である完璧な検出精度を有する。図８Ｂを参照すると、グラフ３７０は、適切な監視性能のための基準を満たさないイベントに対して、機械学習に基づくイベント検出器の誤警報率が高いことを示す。
【0161】
図９Ａは、全ての扉及び窓が閉じている状況に対して機械学習に基づくイベント検出器により行われた決定を示すグラフ３８０である。図９Ｂは、Ｂｏｂの部屋の扉が部屋の内側から開かれる状況に対する決定を示すグラフ３９０である。図９Ｃは、Ａｌｉｃｅの部屋の窓が外側から開かれる状況に対する決定を示すグラフ４００である。グラフ３８０、３９０及び４００において、ｘ軸は決定インデックス又は時間インデックスを表し、決定は一定期間にわたり周期的に行われる。ｙ軸は表６に対応する状態インデックスを表す。これらの例において、システムはイベントｅ１からイベントｅ４への移行を検出することに成功し（図９Ｂに示すように）、イベントｅ１からイベントｅ９への移行を検出することに成功した（図９Ｃに示すように）。
【0162】
上記の結果は、受信機動作性能評価及び運用テストの双方において、適切な監視性能のための基準を満たす対象イベントに対して誤警報率がほぼ０である完璧な検出率を達成する点で、機械学習に基づくイベント検出器が統計値に基づくイベント検出器と同様と性能を有すること、並びに運用テストにおいて対象イベントの発生が検出され且つ報告されることを示す。
【0163】
人間動作検出器：ＴＲＲＳ分散に基づく方法
図２２は、時間反転共振強度分散に基づく人間動作検出器４８０の処理ブロック図を示す。これは、時間反転共振強度の時系列の分散により定量化される大きいダイナミクスが検出されると即座に報告するように設計される。家屋内での人間の動作は、動作が送信機又は受信機の近くで起きる場合は特に、チャネルフェージング及びノイズにより発生する影響と比較して大きい変動を時間反転共振強度の時系列に導入する場合がある。時間反転共振強度分散に基づく人間動作検出器は、オフライントレーニング段階４８２及びオンライン監視段階４８４である２つの段階で動作できる。オフライントレーニング段階４８２において、マルチパス伝搬環境における異なる種類のダイナミクスにより導入される時間反転共振強度分散が調べられ、人間の動作がある場合と人間の動作がない場合のダイナミクスを判定する際に使用するために適切な閾値が取得される。オンライン監視段階４８４において、検出器４８０は測定されたチャネル状態情報を受信し、時間反転共振強度の時系列における分散を判定し、分散とトレーニング段階において確立された閾値とを比較し、家屋内に人間の動作が存在するか否かの判定を行う。
図２２は、時間反転共振強度分散に基づく人間動作検出器４８０の処理ブロック図を示す。これは、時間反転共振強度の時系列の分散により定量化される大きいダイナミクスが検出されると即座に報告するように設計される。家屋内での人間の動作は、動作が送信機又は受信機の近くで起きる場合は特に、チャネルフェージング及びノイズにより発生する影響と比較して大きい変動を時間反転共振強度の時系列に導入する場合がある。時間反転共振強度分散に基づく人間動作検出器は、オフライントレーニング段階４８２及びオンライン監視段階４８４である２つの段階で動作できる。オフライントレーニング段階４８２において、マルチパス伝搬環境における異なる種類のダイナミクスにより導入される時間反転共振強度分散が調べられ、人間の動作がある場合と人間の動作がない場合のダイナミクスを判定する際に使用するために適切な閾値が取得される。オンライン監視段階４８４において、検出器４８０は測定されたチャネル状態情報を受信し、時間反転共振強度の時系列における分散を判定し、分散とトレーニング段階において確立された閾値とを比較し、家屋内に人間の動作が存在するか否かの判定を行う。
【0164】
段階１．オフライントレーニング
トレーニング段階４８２において、システムは、監視領域内に人間の動作が存在しない場合の時間反転共振強度の時系列におけるダイナミクスと比較した人間の動作が存在する場合のダイナミクスの大きさを学習する必要がある。トレーニングにおける詳細なステップを以下に列挙する。
・データの取得：最初に、屋内環境の状況を２つのクラス、すなわち、環境が静的であるＳ０と、監視領域内で動いている人間が少なくとも１人いるＳ１に分ける。双方のクラスに対して、チャネル状態情報サンプルとも呼ばれるチャネル周波数応答が継続的に収集され、以下のように示される。
式中、ｔはチャネル周波数応答が取り込まれる期間の時間インデックスである。チャネル周波数応答が位相歪みにより破損する場合があるため、時間反転共振強度分散に基づく人間動作検出器が時間反転共振強度の時系列のダイナミクスを学習する前に、位相オフセットがそれぞれ個別に補償される。
・ダイナミクスの取得：静的状態Ｓ０及び動作状態Ｓ１の双方に対するチャネル周波数応答測定値の時系列が取得されると、時間ウィンドウ内で時間に沿って時間反転共振強度の分散値を追跡（トラッキング）することにより、時間反転共振強度の時系列のダイナミクスが定量化される。分散の時系列を取得するために、各状態Ｓｉ（ｉ＝０，１）に対して、長さＮ及びオーバラップＮ－１を有するスライディングウィンドウがＨｉ（ｔ）の時系列に適用され、ウィンドウ内で継続的に収集されたチャネル周波数応答をセグメント化する。例えばウィンドウ長Ｎとすると、チャネル周波数応答系列はＨｉ（ｔ０）～Ｈｉ（ｔ０＋Ｎ＊Ｔｓ）であり、Ｔｓはサンプリング時間であり、すなわち各チャネル状態情報の検知の時間間隔はＴｓである。その後、時間ｔ０から始まるウィンドウに対する時間反転共振強度系列が、ｔ０≦ｔ≦ｔ０＋Ｎ＊Ｔｓである
トレーニング段階４８２において、システムは、監視領域内に人間の動作が存在しない場合の時間反転共振強度の時系列におけるダイナミクスと比較した人間の動作が存在する場合のダイナミクスの大きさを学習する必要がある。トレーニングにおける詳細なステップを以下に列挙する。
・データの取得：最初に、屋内環境の状況を２つのクラス、すなわち、環境が静的であるＳ０と、監視領域内で動いている人間が少なくとも１人いるＳ１に分ける。双方のクラスに対して、チャネル状態情報サンプルとも呼ばれるチャネル周波数応答が継続的に収集され、以下のように示される。
・ダイナミクスの取得：静的状態Ｓ０及び動作状態Ｓ１の双方に対するチャネル周波数応答測定値の時系列が取得されると、時間ウィンドウ内で時間に沿って時間反転共振強度の分散値を追跡（トラッキング）することにより、時間反転共振強度の時系列のダイナミクスが定量化される。分散の時系列を取得するために、各状態Ｓｉ（ｉ＝０，１）に対して、長さＮ及びオーバラップＮ－１を有するスライディングウィンドウがＨｉ（ｔ）の時系列に適用され、ウィンドウ内で継続的に収集されたチャネル周波数応答をセグメント化する。例えばウィンドウ長Ｎとすると、チャネル周波数応答系列はＨｉ（ｔ０）～Ｈｉ（ｔ０＋Ｎ＊Ｔｓ）であり、Ｔｓはサンプリング時間であり、すなわち各チャネル状態情報の検知の時間間隔はＴｓである。その後、時間ｔ０から始まるウィンドウに対する時間反転共振強度系列が、ｔ０≦ｔ≦ｔ０＋Ｎ＊Ｔｓである
【0165】
【0166】
段階２．オンライン監視
監視段階において、システムが継続的に環境を感知してチャネル状態情報サンプルを収集するため、チャネル状態情報の時系列Ｈtest（ｔ），ｔ０≦ｔ≦ｔ０＋Ｎ＊Ｔｓ，∀ｔ０を取得するために、トレーニング段階と同一の長さＮ及びオーバラップ長さＮ－１を有するスライディングウィンドウが適用される。位相オフセット及び初期位相歪みを補償した後、時間反転共振強度系列が
として取得され、この時間ウィンドウ内の分散はσtest（ｔ０）である。家屋内に人間の動作があるか否かを判定するために、σtest（ｔ０）値が閾値γmotionと比較される。σtest（ｔ０）≧γmotionの場合のみ、システムは人間の動作の存在を報告する。
監視段階において、システムが継続的に環境を感知してチャネル状態情報サンプルを収集するため、チャネル状態情報の時系列Ｈtest（ｔ），ｔ０≦ｔ≦ｔ０＋Ｎ＊Ｔｓ，∀ｔ０を取得するために、トレーニング段階と同一の長さＮ及びオーバラップ長さＮ－１を有するスライディングウィンドウが適用される。位相オフセット及び初期位相歪みを補償した後、時間反転共振強度系列が
【0167】
時間反転共振強度分散に基づく人間動作検出器を使用することにより、誰かが家屋内を歩いている場合、システムはそれを検出できる。更に、動作が送信機又は受信機の近くで生じる場合、精度は高くなる。内側に金属を有する大きい物体が部屋の中に存在する場合があるため、無線信号（送信機から送出され、受信機により受信される）が動いている人間の身体に到達する前に大部分又はほぼ全ての電磁波を反射する金属の後方で人間の動作が発生する場合、動作検出は妥協される場合がある。
【0168】
図１０Ａを参照すると、グラフ４１０は、家屋の中を歩いている人間がいないことを人間動作検出器が正確に判定したことを示す。図１０Ｂを参照すると、グラフ４２０は、玄関扉領域を歩いている人間がいることを人間動作検出器が正確に判定したことを示す。
【0169】
図１０Ｃを参照すると、グラフ４３０は、決定４３２ａ、４３２ｂ及び４３２ｃで表される特定の期間に広範囲を歩いている人間がいること、並びに決定４３４ａ及び４３４ｂで表される特定の期間に広範囲を歩いている人間がいないことを人間動作検出器が正確に判定したことを示す。図１０Ｄを参照すると、グラフ４４０は、特定の期間に居間領域を歩いている人間がいること及び特定の期間に居間領域を歩いている人間がいないことを人間動作検出器が正確に判定したことを示す。図１０Ｅを参照すると、グラフ４５０は、特定の期間に台所を歩いている人間がいること及び特定の期間に台所を歩いている人間がいないことを人間動作検出器が正確に判定したことを示す。
【0170】
図１０Ｆを参照すると、グラフ４６０は、家屋の外側を歩いている人間がいる場合の決定結果を示す。人間動作検出器は、家屋の中を歩いている人間がいないことを正確に判定した。決定は、家屋の外側を歩いている人間の影響を受けなかった。
【0171】
平滑化のための時間ダイバーシチ
いくつかの例では、屋内環境において、無線送信におけるノイズ及び外のダイナミクスが存在し、受信機側における推定チャネル状態情報に影響を及ぼす場合がある。これらの種類の干渉は、保護システムの誤検出又は誤警報を引き起こす場合がある。いくつかの実現例において、これらの干渉は僅かであり且つ短時間である場合が多いため、時間と矛盾しない決定のみを出力するために時間ダイバーシチ平滑化法が提供される。時間ダイバーシチ平滑化法は、屋内イベントの一般的な持続時間が通常は数秒であるという仮定に依存する。
いくつかの例では、屋内環境において、無線送信におけるノイズ及び外のダイナミクスが存在し、受信機側における推定チャネル状態情報に影響を及ぼす場合がある。これらの種類の干渉は、保護システムの誤検出又は誤警報を引き起こす場合がある。いくつかの実現例において、これらの干渉は僅かであり且つ短時間である場合が多いため、時間と矛盾しない決定のみを出力するために時間ダイバーシチ平滑化法が提供される。時間ダイバーシチ平滑化法は、屋内イベントの一般的な持続時間が通常は数秒であるという仮定に依存する。
【0172】
いくつかの実現例において、時間ダイバーシチ平滑化アルゴリズムは、以下に説明するように２段階の多数決処理を使用する。
【0173】
・段階Ｉの多数決：当段階の多数決は、各単一チャネル状態情報推定値に基づいて生成された未加工の決定に直接適用される。長さＮ１及びオーバラップ長Ｏ１を有するスライディングウィンドウＳＷ１を用いて、インデックスｎの決定は
ＤMV1（ｎ）＝MajorityVote｛Ｄ（１＋（ｎ－１）＊Ｏ１），Ｄ（２＋（ｎ－１）＊Ｏ１），...，Ｄ（Ｎ１＋（ｎ－１）＊Ｏ１）｝により取得される。式中、Ｄ（ｎ）は、イベント検出器又は人間動作検出器のいずれかからのｎ番目のチャネル状態情報測定値に対する決定である。
ＤMV1（ｎ）＝MajorityVote｛Ｄ（１＋（ｎ－１）＊Ｏ１），Ｄ（２＋（ｎ－１）＊Ｏ１），...，Ｄ（Ｎ１＋（ｎ－１）＊Ｏ１）｝により取得される。式中、Ｄ（ｎ）は、イベント検出器又は人間動作検出器のいずれかからのｎ番目のチャネル状態情報測定値に対する決定である。
【0174】
・段階IIの多数決：段階Ｉからの決定ＤMV1（ｎ）が取得されると、第２段階の多数決において、長さＮ２及びオーバラップ長Ｏ２を有する第２のスライディングウィンドウＳＷ２が決定系列ＤMV1（ｎ），ｎ＝１，２，．．．に適用される。この場合、本ステップにおいて、時間ダイバーシチ平滑化を適用後の最終的な決定出力はＤfinal（ｎ）＝MajorityVote｛ＤMV1（１＋（ｎ－１）＊Ｏ２），ＤMV1（２＋（ｎ－１）＊Ｏ２），...，ＤMV1（Ｎ２＋（ｎ－１）＊Ｏ２）｝である。
【0175】
時間ダイバーシチ平滑化を適用することにより導入される時間遅延は、（Ｎ１＋（Ｎ２－１）＊Ｏ１）×Ｔｓ秒である。
【0176】
人間動作検出：アンテナ相関に基づく方法
いくつかの実現例において、アンテナ相関分散に基づく人間動作検出器は、アンテナ相関の分散により定量化される大きいダイナミクスが検出されると即座に報告するように設計される。人間の動作は、動作が送信機又は受信機の近くで起きる場合は特に、チャネルフェージング及びノイズにより発生する影響と比較して大きい変動をアンテナ相関に導入する場合がある。本方法において、人間の動作がある状況と人間の動作がない状況との間のアンテナ相関の分散におけるダイナミクスを区別するために、アンテナ相関の分散に基づく人間動作検出器が使用するのに適した閾値が取得される。
いくつかの実現例において、アンテナ相関分散に基づく人間動作検出器は、アンテナ相関の分散により定量化される大きいダイナミクスが検出されると即座に報告するように設計される。人間の動作は、動作が送信機又は受信機の近くで起きる場合は特に、チャネルフェージング及びノイズにより発生する影響と比較して大きい変動をアンテナ相関に導入する場合がある。本方法において、人間の動作がある状況と人間の動作がない状況との間のアンテナ相関の分散におけるダイナミクスを区別するために、アンテナ相関の分散に基づく人間動作検出器が使用するのに適した閾値が取得される。
【0177】
アンテナ相関を計算する前に、位相補償アルゴリズムを使用して、各リンクにおけるＣＳＩの位相歪みを補償する。図１１を参照すると、いくつかの実現例において、３アンテナＷｉＦｉ装置４９０が送信機として使用され、３アンテナＷｉＦｉ装置５００が受信機として使用される。尚、３アンテナＷｉＦｉ装置は例として使用されるにすぎず、送信機及び受信機は４個以上又は２個以下のアンテナを有することができる。例えば送信機及び受信機はそれぞれ、５個、１０個、２０個、３０個又は３１個以上のアンテナを有することができる。
【0178】
図１１において、Ｈｉ,ｊは、ｉ番目のＴＸアンテナとｊ番目のＲＸアンテナとの間のリンクにおけるチャネル状態情報推定値を表す。アンテナ相関γは、ＴＸアンテナ相関又はＲＸアンテナ相関であってもよく、以下のように計算可能である。
式中、ＭはＴＸアンテナの数を表し、ＮはＲＸアンテナの数を表す。
式（４７）及び式（４８）から、各時間インスタンスにおけるアンテナ相関γｔを図１２に示すように計算できる。
【0179】
背景ノイズ及びハードウェアの障害により、アンテナ相関γｔは多くのノイズを含む場合がある。従って、動作検出にγｔを使用する前に、γｔに対する後処理が必要な場合がある。最初に、以下の式（４９）に示すように、γｔの時間差を取得する。
Δγｔ＝γｔ－γｔ－１ 式（４９）
次に、ウィンドウサイズＷを有する移動平均ウィンドウがΔγｔに適用され、平均値及び分散は以下のように動作検出に使用できる。
Δγｔ＝γｔ－γｔ－１ 式（４９）
次に、ウィンドウサイズＷを有する移動平均ウィンドウがΔγｔに適用され、平均値及び分散は以下のように動作検出に使用できる。
【0180】
ここでは、実験は図１３に示すようなオフィスで行われ、４個の実験が実行された。送信機のロケーションを丸で示し、受信機のロケーションを三角形で示す。各実験において、矩形で示す人間が送信機と受信機との間の空間を繰り返し通過した。以下の実験に対して、ＲＸアンテナ相関を考察する。
【0181】
実験１：小さい部屋における見通し内の場合の動作：第１の実験では、送信機及び受信機の双方が見通し内の状況で小さい部屋に配置され、人間は送信機と受信機との間のリンクを繰り返し通過した。アンテナ相関を図１４のグラフ５１０に示す。図１４のグラフ５１０は、アンテナ相関が時間インデックスに沿って非常に多くのノイズを含むことを示す。従って、アンテナ相関の差分が式（４９）を使用して判定され、結果を図１５のグラフ５２０に示す。
【0182】
ウィンドウサイズＷ＝１秒を選択した場合、Δγｔの統計値を使用して、図１６及び図１７に示すように人間の動作を検出できる。図１６は、送信機と受信機との間の領域を歩いている人間がいる場合のΔγｔの平均を示すグラフ５３０である。図１７は、送信機と受信機との間の領域を歩いている人間がいる場合のΔγｔの分散を示すグラフ５４０である。図１６及び図１７のピークは、検出された人間の動作を表す。
【0183】
第２の実験では、送信機及び受信機をオフィスの入口に配置した。オフィスのカード鍵システムを使用して、入室履歴と歩いている人間がいる場合のＶａｒ（Δγｔ）とが、図１８のグラフ５５０に示すように、時間において整列される。
【0184】
尚、オフィスのカード鍵システムにより提供される入室履歴は、オフィスからの人間の退室を記録しない。従って、Ｖａｒ（Δγｔ）におけるいくつかの未知のピークは、オフィスから退室する人間により生じるものである場合がある。
【0185】
人間動作検出：ノイズレベルに基づく方法
いくつかの実現例において、線形位相オフセット及びランダム初期位相の問題を回避するために、チャネル周波数応答の振幅情報のみが使用される。説明の便宜上、時間ｉにおける（正規化された）チャネル周波数応答の振幅を以下のようなロングベクトルと示す。
いくつかの実現例において、線形位相オフセット及びランダム初期位相の問題を回避するために、チャネル周波数応答の振幅情報のみが使用される。説明の便宜上、時間ｉにおける（正規化された）チャネル周波数応答の振幅を以下のようなロングベクトルと示す。
【0186】
【0187】
【0188】
【0189】
誤警報要件αを与えた場合、閾値ηは次式により判定される。
式中、Ｑ（・）はＱ関数であり、すなわち標準正規累積分布関数の補数である。誤警報確率αが小さい場合、Ｑ-1（α）は大きく、それと同時にηも大きい。上記の説明は、観察ウィンドウＩの長さが上述したアルゴリズムの性能において重要な役割を果たすことを示す。ウィンドウのサイズＩが大きいほどＰ^
Ｎの分散は減少し、閾値ηをより小さい値に設定でき、それにより、式（５６）による観察ウィンドウの持続時間に相当する持続時間を有する動作に対する検出確率を向上できる。しかし、ウィンドウサイズが大きいほど、システムの応答時間は長くなり、持続時間が短い動作に対する検出確率は低下する。異なる目的又は用途に応じて、異なるウィンドウサイズが選択されてもよい。
【0190】
人間動作検出：到着時間（ＴｏＡ）に基づく方法
【0191】
【0192】
ダイバーシチを使用するＴＲＲＳの組み合わせ
上記で定義したように、同一イベント又は異なるイベントにおいて収集されたＣＳＩ間の類似性を特徴づけるために、時間反転共振強度を以下のように計算できる。
ｈ０及びｈ１は２つのフィンガープリントであり、Ｋは使用可能なサブキャリアの総数であり、Ｈ０［ｋ］及びＨ１［ｋ］はサブキャリアｋにおけるＣＳＩであり、ηは同期誤差を補償したｈ０とｈ１との間の変更された相互相関であり、Λ０、Λ１はそれぞれｈ０及びｈ１のチャネルエネルギーである。受信機及び送信機の無線周波数フロントエンド成分の不一致により受信機が送信機と完全に同期しない場合があるため、ｅ-jkφの追加の位相回転を使用して、ηの計算における同期誤差により発生した位相歪みを解消する。その場合、φは以下に説明するアルゴリズム１を使用して推定及び補償可能である。上記の式（６４）は、時間反転共振強度が０から１の範囲であることを示す。更に詳細には、時間反転共振強度が大きいほど、２つのＣＳＩ間の類似性が高いことを示し、従って、２つの関連するイベント間の類似性が高いことを示す。
上記で定義したように、同一イベント又は異なるイベントにおいて収集されたＣＳＩ間の類似性を特徴づけるために、時間反転共振強度を以下のように計算できる。
【0193】
いくつかの実現例において、イベント検出システム又は人間動作検出システムは、周波数ダイバーシチ及び／又は空間ダイバーシチを含むダイバーシチを利用することにより、有効帯域幅の増加を実現できる。周波数ダイバーシチは、異なるＷｉＦｉチャネルを使用することにより達成可能である。更に、周波数ダイバーシチは、直交周波数分割多重システムにおいてサブキャリアにより導入されるダイバーシチを含むことができる。大量のサブキャリア（例えばＩＥＥＥ８０２．１１ｎプロトコルを使用する４０ＭＨｚ帯域幅の場合、１１４個のサブキャリア）により、非常に安定した性能を得られる場合がある。
【0194】
無線通信システムに存在することのあるハードウェア不良を補償することにより、異なるダイバーシチを組み合わせた場合の時間反転共振強度を計算する方法を以下に説明する。一般に、時間反転共振強度は、２つのＣＳＩ間の差分を定量化するメトリックであると考えることができる。２つのＣＳＩは、２つの異なる物理的ロケーションから発生してもよいが、２つの仮想（論理的）ロケーション／状態から発生してもよい。前者の場合、ダイバーシチを有する時間反転共振強度は、ロケーション特定及び追跡に役立つことができる。後者の場合、ダイバーシチを有する時間反転共振強度は、呼吸監視、人間認識及び／又は保護／セキュリティシステム等の用途に使用可能である。
【0195】
現在のＷｉＦｉシステムには２つの異なるダイバーシチが存在し、すなわち周波数ダイバーシチ及び空間ダイバーシチが存在する。ＩＥＥＥ８０２．１１ｎによると、３５個のＷｉＦｉチャネルが、４０ＭＨｚの最大帯域幅を有する２．４ＧＨｚ及び５ＧＨｚの周波数帯域におけるＷｉＦｉ送信に割り当てられる。ＷｉＦｉチャネルが多数であることで周波数ダイバーシチが得られ、深いフェージング又は重大な干渉が発生した場合にＷｉＦｉ装置が周波数ホッピングを実行する機会が与えられる。それに対して、空間ダイバーシチは多入力多出力ＷｉＦｉ装置で利用可能であり、これはスペクトル効率を大幅に向上する高度な技術である。多入力多出力装置はＩＥＥＥ８０２．１１ｎ／ａｃの有用な構成要素になっており、多くの市販のＷｉＦｉ装置で広範に使用されている。ＷｉＦｉシステムの場合、双方の種類のダイバーシチを収集して、粒度がはるかに細かいフィンガープリントを提供することにより、僅か４０ＭＨｚの帯域幅で測定されたフィンガープリントと比較して曖昧さを低減することができる。
【0196】
図２３Ａ～図２３Ｃは、周波数ダイバーシチ及び空間ダイバーシチを個別に又は一緒に使用することにより広い有効帯域幅を生成するための一般的な原理を示す。図２３Ａは、ＷｉＦｉ装置５６０が１つのアンテナを有する例を示す。図２３Ｂは、ＷｉＦｉ装置５７０が４つのアンテナを有する例を示す。図２３Ｃは、ＷｉＦｉ装置５８０が２つのアンテナを有する例を示す。ＷｉＦｉ装置は複数のＷｉＦｉチャネルで動作できるため、周波数ホッピングを実行して異なるＷｉＦｉチャネルにおけるＣＳＩを取得することにより、周波数ダイバーシチを利用できる。図２３Ａに示すように、４つの異なるＷｉＦｉチャネルにおけるＣＳＩは連結され、広い有効帯域幅のフィンガープリントを形成する。周波数ダイバーシチは単アンテナＷｉＦｉ装置で利用可能であるが、周波数ホッピングを実行するには時間がかかる。時間効率のために、多アンテナＷｉＦｉ装置では空間ダイバーシチを利用できる。４つのアンテナを有する図２３ＢのＷｉＦｉ受信機５７０の場合、４つの受信アンテナにおけるＣＳＩを組み合わせて、広い有効帯域幅のフィンガープリントを形成できる。図２３Ｃの例では、周波数ダイバーシチ及び空間ダイバーシチの双方が使用され、２つのＷｉＦｉチャネルにおけるＣＳＩ及び２つの受信アンテナからのＣＳＩが組み合わされてフィンガープリントを形成する。
【0197】
ＷｉＦｉシステムの場合、空間ダイバーシチはアンテナリンクの数により判定されるが、周波数ダイバーシチは使用可能なＷｉＦｉチャネルの数に依存する。最大空間ダイバーシチをＳとし、最大周波数ダイバーシチをＦとし、各ＷｉＦｉチャネルの帯域幅をＷとすると、チャネル状態情報測定値は
と書くことができ、Ｓ×Ｆ×Ｗの有効帯域幅を有する粒度の細かいチャネル状態情報フィンガープリントを提供できる。
【0198】
【0199】
以下に示すアルゴリズム１は、
の計算を詳細に説明する。アルゴリズム１において、ステップ４～９は仮想リンク（ｓ，ｆ）におけるチャネルエネルギーを計算し、ステップ１０～１４は仮想リンク（ｓ，ｆ）における２つのフィンガープリントの変更された相互相関を計算する。各仮想リンクにおけるチャネルエネルギー及び変更された相互相関はそれぞれ、ステップ９及びステップ１５に示すように累算される。最後に、時間反転共振強度がステップ１８により取得される。ηs,fの計算は、サイズＮの離散フーリエ変換と同一の形式をとる
により近似され、従って、高速フーリエ変換により効率的に計算可能である。計算において大きいＮを使用することで、ηs,fのより正確な近似を得られる。
【0200】
アルゴリズム１において開示した時間反転共振強度の組み合わせに加えて、時間反転共振強度を組み合わせる他の方法が考えられる。更に簡潔に記述すると、チャネルｈ０,ｄとチャネルｈ１,ｄとの間のｄ番目のリンクにおける（正規化）共振強度は以下のように再定義される。
【0201】
結果として、各リンクｄに対して３つの特徴が得られる。
１．正規化共振強度TR（ｈ０,ｄ，ｈ１,ｄ）
２．リンクｄにおけるチャネルｈ０のチャネルエネルギー：Λ０,ｄ
３．リンクｄにおけるチャネルｈ１のチャネルエネルギー：Λ１,ｄ
１．正規化共振強度TR（ｈ０,ｄ，ｈ１,ｄ）
２．リンクｄにおけるチャネルｈ０のチャネルエネルギー：Λ０,ｄ
３．リンクｄにおけるチャネルｈ１のチャネルエネルギー：Λ１,ｄ
【0202】
次の計算ステップは、各リンクの（正規化）共振強度を組み合わせて一般的特徴を形成することである。時間反転共振強度を組み合わせる４つの方法を以下に説明する。
【0203】
最初に、チャネルｈ０とチャネルｈ１との間の共振強度の基本的組み合わせは次式である。
式中、Ｄはアンテナリンクの総数である。
【0204】
本方法において、各リンクは品質に関係なく平等に扱われる。この組み合わせアルゴリズムにより、位置特定及び認識における性能の評価基準を得られることが期待できる。
【0205】
第２に、チャネルｈ０とチャネルｈ１との間の時間反転共振強度のエネルギー加重組み合わせは次式である。
【0206】
前述の方法と異なり、エネルギー加重組み合わせ法は、チャネル利得の積Λ０,ｄΛ１,ｄに従って異なるリンクを組み合わせる。より高い品質（すなわち、より高いチャネル利得）を有するリンクは、他のリンクより加重される。一般に、エネルギー加重組み合わせ法は基本的組み合わせ法より高い性能を得られる。
【0207】
第３に、チャネルｈ０とチャネルｈ１との間の共振強度の全リンク連結組み合わせは次式である。
【0208】
本方法の背景にある考えは、各リンクチャネルを個別に整列させた後、使用可能なリンクにおける全てのＣＳＩを連結して、新規のチャネル状態情報を効果的に形成することである。
【0209】
第４に、チャネルｈ０とチャネルｈ１との間の共振強度の部分リンク連結組み合わせは次式である。
式中、ＮＲＸは受信アンテナの数であり、Ｄｉは受信アンテナｉを含むアンテナリンクのセットを表す。例えば３×３多入力多出力システムにおいて、アンテナリンク１、４、７はそれぞれ、送信アンテナ１と受信アンテナ１との間のリンク、送信アンテナ２と受信アンテナ１との間のリンク、送信アンテナ３と受信アンテナ１との間のリンクを表す。従ってＤ１＝｛１，４，７｝である。本方法の背景にある考えは、最初に同一の受信アンテナにおけるＣＳＩを連結して時間反転共振強度を計算して、全ての受信アンテナにおける計算された時間反転共振強度の平均を得ることである。
【0210】
チャネル状態情報の位相オフセットの補償
チャネル状態情報は、（ｉ）搬送周波数オフセット（ＣＦＯ）α及び共通位相オフセット（ＣＰＥ）βにより生じる初期位相歪みと、（ii）シンボルタイミングオフセット（ＳＴＯ）σ及びサンプリング周波数オフセット（ＳＦＯ）δにより生じる線形位相歪みによる影響を受ける場合がある。
チャネル状態情報は、（ｉ）搬送周波数オフセット（ＣＦＯ）α及び共通位相オフセット（ＣＰＥ）βにより生じる初期位相歪みと、（ii）シンボルタイミングオフセット（ＳＴＯ）σ及びサンプリング周波数オフセット（ＳＦＯ）δにより生じる線形位相歪みによる影響を受ける場合がある。
【0211】
ランダム性は、以下のようにしてチャネル状態情報フィンガープリントの位相品質を歪める。
式中、φｋ＝α＋δｋである。
【0212】
従って、
で定義される時間反転共振強度の計算は、位相歪みを補償するための「ｍａｘ」演算を含み、非常に複雑である。
【0213】
位相歪みを補償しつつ複雑さを軽減するために、タイプＩのチャネル状態情報フィンガープリント等のチャネル状態情報フィンガープリントを使用できる。
【0214】
例えば、以下のように書くことができる。
式中、Ｈ１,１［ｋ］は送信アンテナ１から受信アンテナ１において受信されたチャネル状態情報フィンガープリントを表し、Ｈ２,１［ｋ］は送信アンテナ２から受信アンテナ１において受信されたチャネル状態情報フィンガープリントを表す。
【0215】
各サブキャリアｋに対するアンテナ間フィンガープリントＧ１,２［ｋ］を以下のように形成できる。
上記は、シンボルタイミングオフセットσ及び共通位相オフセットβの影響が相殺されたことを示す。α及びδは残留した搬送周波数オフセット及びサンプリング周波数オフセットであるため、多くの場合に無視できる。
【0216】
３×３多入力多出力システムにおける各無線周波数チェーンに対して、３つの新規のチャネル状態情報フィンガープリントを形成できる。従って、全体で、Ｇ1,2,...,9を形成できる。更に一般に、
の新規のチャネル状態情報フィンガープリントを形成できる。
【0217】
次に、単純な内積を適用して共振強度を取得する。例えば較正後の２つのＣＳＩ間の共振強度（リンク毎）を以下のように計算する。
【0218】
別の実施形態において、以下のフィンガープリント（チャネル状態情報フィンガープリントII）を使用できる。ＴＸアンテナｉ及びＲＸアンテナｊに対して、以下のように書くことができる。
これは、線形位相シフトを完全に除去する。
【0219】
ＣＳＩ位相のサニタイズ：サブキャリアｋに対するグランドトゥルースチャネル状態情報を
式中、αは初期位相オフセットを表し、βは線形位相オフセットを表し、ｎｋはサブキャリアｋにおける推定雑音である。
【0220】
線形位相オフセットのクリーニング：最初に、線形位相オフセットαに注目する。いくつかの例において、αを除去する１つの方法は、全てのサブキャリアに沿って位相成分にわたり線形フィッティングを行い、推定線形位相オフセットとして勾配を抽出することである。これは解析では正しいが、ノイズの多いアンラッピングのため、線形フィッティング法は実際の実現例では正しくない場合がある。
【0221】
アンラッピングを用いる線形位相較正：未加工のチャネル状態情報位相は連結されて［－π，π］となり、これをラッピング位相φｗと呼ぶ。全てのサブキャリアに沿ってチャネル状態情報位相を線形フィッティングするためには、最初に各サブキャリアにおける位相を以下のようにアンラッピングする必要がある。
φ［ｋ］＝φ［ｋ－１］＋Ｗ｛φｗ［ｋ］－φｗ［ｋ－１］｝ 式（７９）
式中、φ［ｋ］はサブキャリアｋにおけるアンラッピング位相を示し、φｗ［ｋ］はラッピング位相である。演算Ｗはアンラッピング関数である。
φ［ｋ］＝φ［ｋ－１］＋Ｗ｛φｗ［ｋ］－φｗ［ｋ－１］｝ 式（７９）
式中、φ［ｋ］はサブキャリアｋにおけるアンラッピング位相を示し、φｗ［ｋ］はラッピング位相である。演算Ｗはアンラッピング関数である。
【0222】
ラッピング位相φ［ｋ］に更なる位相φｎ［ｋ］を追加するノイズのため、
であり、アンラッピング法はＷ｛φｗ［ｋ］＋φｎ［ｋ］－φｗ［ｋ－１］｝に劣化する場合がある。式（８０）におけるアンラッピング手順の誤トリガ又はミストリガにより、少なくとも|２π－φｎ［ｋ］|の誤りが振幅に導入される場合がある。
【0223】
従って、アンラッピング動作が破損したため、線形フィッティングにより取得される線形位相オフセットは正しくない場合がある。同一のロケーション（状態）からのＣＳＩに対する線形フィッティング方法の一例を図２４に示す。アンラッピング動作により誤りが導入される。
【0224】
実験値から、電力において
である場合、すなわち信号対ノイズ比（ＳＮＲ）が小さい場合、破損した位相アンラッピングを得る可能性が高いことが更にわかった。
【0225】
アンラッピングを用いない線形位相較正：位相をアンラッピングせず且つ特別なハードウェア実現例を用いずに線形位相オフセットを較正するために、以下の方法を使用する。
【0226】
【0227】
【0228】
式（８３）に示すように、真の線形位相オフセットβ以外に剰余項βｒｅｓが得られ、これは各ロケーションに対して以下のようにほぼ固定される。
式中、ｈ［ｋ］Ｈｋはグランドトゥルースチャネル状態情報である。各サブキャリアにおける位相が独立して無相関であるという仮定によると、βｒｅｓは非常に小さく、無視できる。
【0229】
線形位相較正方法の例を図２５に示す。
【0230】
初期位相オフセットのクリーニング：線形位相オフセットβ^がチャネル状態情報ベクトルに対して除去されると、次のステップは、初期位相オフセットαを除去することである。
【0231】
【0232】
線形位相オフセットのクリーニングと同様の方法を使用して、αパラメータが以下のように推定される。
【0233】
βｒｅｓと同様に、初期位相の剰余αｒｅｓは全てのサブキャリアにわたる共通位相として近似的に推定される。
この場合、ノイズの独立性が考慮される。各サブキャリアにおける位相が独立していると仮定すると、αｒｅｓは小さいはずである。
【0234】
上述したように、線形位相較正及び初期位相較正は並行して実現可能であり、双方は高速で効率的な配列演算である。上述した位相較正法は、従来技術のチャネル状態情報整列アルゴリズムで必要とされる参照チャネル状態情報を必要とせずに、各チャネル状態情報に対する位相を較正して補償する。実験結果から、ＣＳＩが同一のロケーション／状態から取得される場合、較正後のＣＳＩは全て一致し、適切に整列されている。
【0235】
上述したアルゴリズムに従って位相が較正され且つサニタイズされると、残りのチャネル状態情報は、位相歪みを有さないグランドトゥルースチャネル状態情報であると見ることができる。その結果、共振強度を計算するためのＦＦＴに基づく方法は不要になる。位相較正を行った後、取得されるＣＳＩは、ＡｏＡ又はＴｏＦを推定するために使用されてもよいグランドトゥルースチャネル状態情報を維持する。更に、較正されたＣＳＩの共振強度を計算するために、単純な内積法が使用される。この場合、共振強度は、送信中に伝搬する屋内環境に対する電磁波応答であるため、その物理的意味を維持する。
【0236】
図２６に関連して上述したように、会場１４０におけるイベントを検出するために、送信機１４２はマルチパスチャネル１４６を介して受信機１４４に無線信号を送出する。受信機１４４がチャネル特性を正確に推定するための送信機１４２と受信機１４４との間のハンドシェイク処理の例を以下に説明する。
【0237】
最小二乗（ＬＳ）に基づくチャネルプロービングハンドシェイク
チャネルプロービングのための最小二乗（ＬＳ）に基づくハンドシェイク方法を以下に説明する。本方法において、送信機１４２は受信機１４４にプローブ信号を送出する。受信機１４４は受信信号と同期し、受信信号に基づいてチャネル応答を推定する。
チャネルプロービングのための最小二乗（ＬＳ）に基づくハンドシェイク方法を以下に説明する。本方法において、送信機１４２は受信機１４４にプローブ信号を送出する。受信機１４４は受信信号と同期し、受信信号に基づいてチャネル応答を推定する。
【0238】
図２７を参照すると、いくつかの実現例において、送信機１４２により送信されるチャネルプロービング信号１６０のフレームは、４つの連続する部分、すなわちプリアンブル１６２、バーカ符号１６４、ペイロード１６６及びトレーニング系列１６８を含む。各フレームの長さは本方法では固定され、送信機１４２及び受信機１４４の双方により認識される。いくつかの例において、受信機１４４がより正確なチャネル情報及びタイミング情報を取得できるように、送信機１４２はチャネルプロービング信号の２つ以上のフレームを受信機１４４に送出してもよい。いくつかの実現例において、タイミングの精度はプリアンブル１６２及びバーカ符号１６４の長さを増加することにより向上でき、チャネル情報の精度はトレーニング系列１６８の長さを増加することにより向上できる。
【0239】
プリアンブル１６２は、例えばパルス位置変調（ＰＰＭ）を使用して変調される。パルス位置変調シンボル「０」及び「１」を図２８に示す。プリアンブル１６２は連続するパルス位置変調シンボル「０」を含む。パルス位置変調シンボル「０」において、ＬはＬ＝τｓ／δｔにより定義されるマルチパスチャネルの長さであり、τｓは遅延スプレッドであり、δｔはシステムのサンプリング周期である。実際、チャネル長が正確に認識されない場合でも、受信機１４４は環境の一般的な遅延スプレッド及びシステムのサンプリングレートに基づいてチャネル長を推定できる。本方法において、推定チャネル長が実際のチャネル長より長い場合、性能に影響はない。従って、受信機１４４はチャネル長を過大評価でき、誤りを生じることなくＬを実際のチャネル長より大きく設定できる。
【0240】
送信機１４２により送信されるバーカ符号１６４は、送信機１４２及び受信機１４４の双方により認識される符号である。バーカ符号１６４の長さはＬｂである。送信機１４２により送信されるペイロード１６６は、パルス位置変調を使用して変調される符号化情報であってもよい。情報は、トレーニング系列１６８のインデックスを含んでもよい。ペイロード１６６の長さはＬｐである。
【0241】
図２９を参照すると、トレーニング系列１６８は、ガードインターバルＡ１７０、有効トレーニング系列（ＥＴＳ）１７２及びガードインターバルＢ１７４を含む３つの連続する部分を有する。いくつかの例において、ガードインターバルＡ１７０は長さＬａを有する「０」ビットの系列であり、有効トレーニング系列１７２は長さＬｅを有する「０」ビット及び「１」ビットの系列であり、ガードインターバルＢ１７４は、ＬＢ≧Ｌである長さＬＢを有する「０」ビットの系列である。
【0242】
送信機１４２から送信されたチャネルプロービング信号は、マルチパスチャネル１４６を通って受信機１４４に到達する。受信信号は、ノイズにより汚染されたチャネルプロービング信号とチャネル応答との畳み込みである。受信機１４４は、受信信号を使用してチャネル応答を推定する。
【0243】
図３０を参照すると、チャネルプロービング信号を受信した場合、受信機１４４は、チャネルプロービング信号のプリアンブル１６２に対してスライディングウィンドウエネルギー検出を使用して同期インデックスを判定する。グラフ１８０は、受信機１４４において受信される信号の一例を示す。スライディングウィンドウエネルギー検出法は、ｌｐがＰＰＭシンボルの長さであるウィンドウサイズＬｗ＝ｌｐ＋Ｌを有する移動ウィンドウ１８２を使用する。ＥＩはｉ番目、（ｉ＋Ｌｗ）番目、（ｉ＋２Ｌｗ）番目、...、（ｉ＋（γ－１）Ｌｗ）番目のサンプルから開始するウィンドウにより累積された平均エネルギーであると定義する。この場合、γは平均ウィンドウの数であり、ｉ∈｛０，１，２，...，Ｌｗ－１｝である。同期インデックスｉｓはＥｉを最大化するインデックスｉであり、ｉ＝argmaxｉＥｉと書くことができる。この場合、インデックスｉは受信信号内のｉ番目のサンプルに対応する。同期インデックスｉｓが判定された後、他のインデックスは同期インデックスｉｓを参照として使用して判定される。
【0244】
図３１を参照すると、同期インデックスｉｓを取得した後、受信機１４４は同期を確認するためにバーカ符号１６４の検出を開始する。受信機１４４は、受信信号のインデックスｉｓ＋ｍＬｗからインデックスｉｓ＋ｍＬｗ＋Ｌｂ－１までの部分を復調し、復調された信号の部分と事前格納済みのバーカ符号系列とを比較する。この場合、Ｌｂはバーカ符号１６４の長さであり、ｍ∈｛０，１，．．．，Ｎ｝であり、Ｎはユーザにより設定される探索範囲である。いくつかの例において、Ｎはプリアンブル１６２内のパルス位置変調シンボルの数より大きく設定される。復調信号が事前格納済みバーカ符号とマッチングする場合、送信機は、受信機がチャネルプロービング段階にあり、ペイロード及びトレーニング系列が受信信号内のバーカ符号に後続することを保証される。更に、送信機は、バーカ符号１６４が開始する受信信号内の位置がマッチングインデックスｉｓ＋ｍＬｗであること及びバーカ符号１６４が終了する受信信号内の位置がインデックスｉｓ＋ｍＬｗ＋Ｌｂ－１であることを認識する。当該情報により、受信機１４４はチャネル応答を推定するための情報を取得できる。
【0245】
ペイロード１６６がバーカ符号１６４に後続するため、受信機１４４は、ペイロード１６６がインデックスｉｓ＋ｍＬｗ＋Ｌｂから開始すると仮定する。ペイロード１６６は固定長Ｌｄを有するため、受信機１４４は受信チャネルプロービング信号のインデックスｉｓ＋ｍＬｗ＋Ｌｂからインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄ－１までの部分を復調し、復調部分を復号化してトレーニング系列のインデックスｑを取得する。この場合、インデックスｑは受信信号内のｑ番目のサンプルを示さない。寧ろ、送信機１４２及び受信機１４４の双方により認識される２つ以上のトレーニング系列のセットが存在し、送信機１４２はセット内のいずれかのトレーニング系列を選択できる。インデックスｑは、トレーニング系列のセットの中で送信機１４２により選択されたトレーニング系列を示す。
【0246】
受信機１４４は、トレーニング系列が受信信号のインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄから開始すると仮定し、受信信号のインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄからインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄ＋Ｌｔ－１までの部分を使用してチャネル応答を推定する。この場合、Ｌｔ＝Ｌｅ＋Ｌである。ｕは受信信号のインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄからインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄ＋Ｌｔ－１までの部分を示すとする。受信機１４２は、ペイロード１６６から取得したインデックスｑに対応する変換行列Ｃｑを使用して、ｕを推定チャネル応答に変換する。この変換を以下のように書くことができる。
【0247】
変換行列Ｃｑは、受信機により使用されるトレーニング系列ｓｑにより以下のように構成されうる。
Ｃｑ＝ (Ｔ'Ｔ)-1Ｔ'
式中、Ｔはｓｑにより生成されるテプリッツ行列であり、送信機１４２により使用される有効トレーニング系列である。行列Ｔの１列目はｓｑｅ＝［ｓｑ；０］であり、０はＬ個の０で構成されるベクトルである。図３２を参照すると、テプリッツ行列Ｔ１９０の第１の列は、ｓｑｅから取得される第１の要素以外は全て０で構成される。行列Ｔ１９０は次元ｃｍ×ｃｎを有し、ｃｍ＝Ｌｅ＋Ｌ及びｃｎ＝Ｌａ＋τ＋Ｌである。パラメータτは、先行ステップにおける同期ミスマッチの最大値より大きい必要のある補償係数であり、すなわち、トレーニング系列が開始する受信信号におけるインデックスと受信機１４４により検出されたトレーニング系列１６８の始点であるインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄとの間の差分である。
Ｃｑ＝ (Ｔ'Ｔ)-1Ｔ'
式中、Ｔはｓｑにより生成されるテプリッツ行列であり、送信機１４２により使用される有効トレーニング系列である。行列Ｔの１列目はｓｑｅ＝［ｓｑ；０］であり、０はＬ個の０で構成されるベクトルである。図３２を参照すると、テプリッツ行列Ｔ１９０の第１の列は、ｓｑｅから取得される第１の要素以外は全て０で構成される。行列Ｔ１９０は次元ｃｍ×ｃｎを有し、ｃｍ＝Ｌｅ＋Ｌ及びｃｎ＝Ｌａ＋τ＋Ｌである。パラメータτは、先行ステップにおける同期ミスマッチの最大値より大きい必要のある補償係数であり、すなわち、トレーニング系列が開始する受信信号におけるインデックスと受信機１４４により検出されたトレーニング系列１６８の始点であるインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄとの間の差分である。
【0248】
図３３は、ＳＮＲ＝０ｄＢ及びチャネル長Ｌ＝２００である１０００個の実現例に対する同期ミスマッチヒストグラムの一例を示すグラフ２００である。ミスマッチ量はチャネル長と比較して小さく、何らかの特定の値により制限される。τが大きいほど多くのシステムリソースを使用するため、τの値は、同期誤差を制限できる可能な限り小さい値になるように選択されうる。
【0249】
推定チャネル応答が実際のチャネル応答のサイズであるＬより大きいサイズを有することが、
及び変換行列Ｃｑの次元からわかる。長いチャネル応答が推定される理由は、同期誤差により、推定チャネル応答は実際のチャネル応答がシフトされたバージョンであるためである。長いチャネル応答を推定することにより、有効チャネル応答が推定の中に含まれることが保証されるため、重要なタップが失われず、このことは、受信装置において送信波形をフォーカスするためにマルチパスチャネルを整合フィルタとして使用する時間反転システムにとって重要である。
【0250】
有効トレーニング系列が特に選択される場合、チャネル応答は変換行列Ｃｑを使用せずに容易に推定されうる。例えば有効トレーニング系列がＤ≧ＬであるＤにより分離されるインパルス列である場合、チャネル応答は受信信号内のサイズＤのウィンドウを平均化することにより単純に推定されうる。ｖｉが受信信号のインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄ＋ｉＤからインデックスｉｓ＋ｍＬｗ＋Ｌｂ＋Ｌｄ＋ｉＤ＋Ｄ－１までの部分を示すとすると、推定チャネル応答を以下のように書くことができる。
【0251】
相関に基づくチャネルプロービングハンドシェイク
チャネルプロービングに対する相関に基づくハンドシェイク方法を以下に説明する。本方法において、時間同期及びチャネル応答推定の双方のために単一系列が使用される。
チャネルプロービングに対する相関に基づくハンドシェイク方法を以下に説明する。本方法において、時間同期及びチャネル応答推定の双方のために単一系列が使用される。
【0252】
送信機１４２は、マルチパスチャネル１４６を介して受信機１４４へ擬似ランダム（ＰＮ）系列ｘを送信する。系列ｘは、例えば以下の式に示す確率質量関数（ｐｍｆ）を有する分散ランダム変数の列であってもよい。
【0253】
ここで使用される擬似ランダム符号の１つの重要な特徴は、長い擬似ランダム符号の自己相関関数がデルタ関数である傾向があることである。例えば図３４を参照して１０００の長さを有する擬似ランダム符号の自己相関関数を考えると、自己相関関数はｍ＝０において非常に高い値を有するが、他のロケーションにおいて大幅に抑制される。
【0254】
以下の説明において、複数の記号が使用される。Ｌはチャネル応答の長さであり、Ｎは擬似ランダム符号の長さであり、αはチャネルプロービングに対する閾値に関するパラメータである。
【0255】
以下の式に示すように、受信信号ｙは擬似ランダム符号ｘとチャネルｈとの畳み込み結果である。
【0256】
以下の式において定義される相互相関関数ｒｙｘ（ｍ）は、擬似ランダム符号ｘ及びチャネルｈの自己相関の畳み込みに等しい。
【0257】
相関による同期のためのアルゴリズムを以下に説明する。送信機は相互相関ｒyx（ｍ）を計算し、その最大振幅を見つける。
ｃｏｒｒmax＝ｍａｘ｜ｒyx（ｍ）｜
ｃｏｒｒmaxに基づいて、送信機は受信信号ｙの同期インデックスｉｓを以下のように見つける。
ｃｏｒｒmax＝ｍａｘ｜ｒyx（ｍ）｜
ｃｏｒｒmaxに基づいて、送信機は受信信号ｙの同期インデックスｉｓを以下のように見つける。
【0258】
パラメータαは実験により選択可能であり、０～１の範囲である。例えばシミュレーションで使用される超広帯域チャネルモデルによると、パラメータαは約０．２であるように選択される。インデックスｉｓを取得した後、以下のように相互相関関数ｒyx（ｍ）を切り捨てることにより推定チャネル応答
【0259】
例示的な一実施形態において、受信器１４４は無線チャネル１４６に対するチャネル状態情報（ＣＳＩ）を取り込むために使用される回路を含んでもよい。チャネル状態情報は、無線信号におけるサブキャリア毎の無線チャネル応答を推定するために使用されてもよい。チャネル状態情報は、無線チャネルに関する振幅情報及び位相情報の双方を含む複素ベクトルであってもよい。いくつかの実施形態において、チャネル振幅情報のみが判定され且つ／又はエクスポートされ且つ／又は処理されてもよい。いくつかの実施形態において、チャネル位相情報のみが判定され且つ／又はエクスポートされ且つ／又は処理されてもよい。好適な実施形態において、チャネル振幅情報及びチャネル位相情報の双方がイベント監視システムの用途の一部として判定され且つ／又はエクスポートされ且つ／又は処理されてもよい。実施形態において、チャネル状態情報信号は、マイクロプロセッサ及び／又はＤＳＰ等のプロセッサ、コンピュータ、ＦＰＧＡ及びＡＳＩＣ等において処理されてもよく、チャネル状態情報信号を生成するチップ又はシステムと一体化されるか又は非常に近接するリソース上でローカルに処理されてもよく、且つ／又はチャネル状態情報信号は信号処理を実行するために構成されてもよい距離計算装置においてリモートに処理されてもよく、且つ／又はチャネル状態情報信号は少なくともいくつかのクラウドリソースを使用して処理されてもよい。
【0260】
実施形態において、チャネル状態情報は、８０２．１１信号伝送プロトコルに対して存在するような正規化機構により判定されてもよい。いくつかの実施形態において、チャネル状態情報は、無線受信機においてパラメータを監視し且つ記録すること及び／又は受信機と送信機との間で渡されるプロービング信号又はトレーニング系列又はトレーニングプリアンブルを監視することにより判定されてもよい。実験室レベルでの実証において、チャネル状態情報を取り込むために使用された回路は、ＵＳＲＰ（National Instruments' Universal Software Radio Peripheral platform）及び／又はＩｎｔｅｌ５３００チップを含んだ。
【0261】
実施形態において、無線プロトコルパケット構成の特定の部分は、チャネル状態情報及び／又はチャネル周波数応答を推定するために使用されてもよい。例示的な８０２．１１ネットワークにおいて、ショートプリアンブルの後のロングプリアンブルがチャネル状態情報を判定するために使用されてもよい。例示的な８０２．１１通信システムは、直交周波数分割多重（ＯＦＤＭ）を使用してもよい。各直交周波数分割多重チャネルは、複数のサブキャリアを含んでもよい。直交周波数分割多重チャネルは、２１個以上のサブキャリア、３１個以上のサブキャリア、４１個以上のサブキャリア、４８個のサブキャリア、５６個のサブキャリア、６４個のサブキャリア又は１１４個のサブキャリア等を含んでもよい。実施形態において、直交周波数分割多重信号におけるサブキャリアのいくつかはヌルサブキャリアであってもよく、あるいはパイロットサブキャリアであってもよい。実施形態において、直交周波数分割多重シンボルにおけるサブキャリアのうちのいくつかのみがデータサブキャリアであってもよい。
【0262】
実施形態において、回路はサブキャリアのうちのいくつか又は全てに対するチャネル状態情報を提供してもよい。例えばいくつかの例示的な実施形態において、Ｉｎｔｅｌ５３００チップ及び／又はＵＳＲＰは直交周波数分割多重チャネルにおける３０個のサブキャリアに対するチャネル状態情報をエクスポートしてもよい。実施形態において、不完全なタイミング及び／又は周波数同期に起因するチャネル状態情報の変動を緩和するために、チャネル状態情報の更なる処理が実行されてもよい。実施形態において、１つのサブバンドからのチャネル状態情報がベクトルとして構成されてもよく、これを帯域内フィンガープリントと呼んでもよい。いくつかの実施形態において、１つのサブバンドからのチャネル状態情報は複数のチャネル状態情報判定値の平均であってもよく、いくつかの判定値がトリミング又は廃棄されてもよい。実施形態において、所定のチャネルに対するサブバンドチャネル状態情報は、１回、２回、１０回、３０回、５０回及び１００回等、測定され且つ／又は判定されてもよい。実施形態において、サブバンドチャネル状態情報のうちのいくつか又は全てが帯域内フィンガープリントを形成するために使用されてもよい。
【0263】
チャネル状態情報が特定の数のサブキャリアに対して報告される上述の例は一例にすぎず、限定することを全く意図しない。あらゆる数のサブキャリアからのチャネル状態情報が本明細書中で説明するイベント監視システムにおいて使用されてもよい。また、チャネル状態情報及び帯域内フィンガープリントを判定するために実行されてもよいチャネル測定の数又は交換されるＷｉＦｉパケットの数は例示であり、限定することを全く意図しない。例えばイベントを監視するイベント監視システムは、推定チャネル応答及びイベントを判定するために使用可能なプローブ信号を常時送出してもよい。
【0264】
実施形態において、複数の異なる中心周波数、周波数帯域幅及び／又はアンテナの帯域内フィンガープリントが組み合わされ且つ／又は連結されて複合フィンガープリント又は帯域間フィンガープリントにされてもよい。実施形態において、帯域間フィンガープリントが複合チャネル応答を形成してもよい。当該複合チャネル応答は、会場内で発生しているイベントを判定するために使用されてもよい。実施形態において、複合チャネル応答を使用するイベント監視システムの精度は、単一の帯域内フィンガープリントを使用して達成される精度より高い場合がある。
【0265】
例示的な一実施形態において、無線装置（例えば、送信機１４２及び受信機１４４）間の無線チャネルに対するチャネル状態情報を判定するために、８０２．１１フレーム内の２つのロングプリアンブルを使用できる。例えば一方の無線装置は無線ルータ／アクセスポイントであってもよく、無線ホットスポット（例えば、自宅又は職場の）をサポートしてもよい。無線ルータ／アクセスポイントは、本明細書中で説明するシステム及び方法と動作するように特に変更及び／又は構成されていない標準的な装置であってもよい。換言すると、本明細書中で説明するシステム及び方法は、イベント情報を判定するために、そのような既存の無線装置及びそれらの手法を利用できる。この例示的な実施形態において、第２の無線装置（例えば、コンピュータ、スマートフォン、タブレット等）がルータ／アクセスポイントの範囲内に存在する場合、第２の無線装置はルータ／アクセスポイントからＷｉＦｉ信号を受信してもよく、ＷｉＦｉ信号に含まれるロングプリアンブルを処理して第２の装置とルータ／アクセスポイントとの間の環境固有シグニチャを判定してもよい。いくつかの実現例において、第２の無線装置は、第２の無線装置がシステムと動作できるようにするソフトウェア及び／又はファームウェアを含んでもよい。いくつかの実現例において、第２の無線装置はホットスポットに接続されない。寧ろ、第２の無線装置は、ルータ／アクセスポイントから標準的なプリアンブル情報を受信できるパッシブ装置であってもよい。このように、第２の無線装置を「スニファ（ｓｎｉｆｆｅｒ）」と呼んでもよい。実施形態において、２つのロングプリアンブルは１つのフレームで、例えば直交周波数分割多重チャネル毎に送出される。実施形態において、直交周波数分割多重チャネルは、例えば２０ＭＨｚ又は４０ＭＨｚの帯域幅を有してもよい。
【0266】
上述したように、システムはチャネル状態情報を判定し且つチャネル状態情報と種々のイベントとを関連付けることができる。例えば１つの例において、第２の無線装置は、チャネル状態情報を格納し且つそれを特定のイベント（例えば、「扉が開いて窓が閉じた」）に関連付けることができる。特定のイベントがその後に発生する場合、第２の装置はインジケーション（例えば、警報）を提供してもよい。いくつかの実現例において、第２の無線装置が（例えば、ＷｉＦｉネットワークに加わることにより）インターネットにアクセスできる場合、第２の無線装置は取得した環境情報（例えば、イベントに関連付けられたチャネル状態情報）をアップロードしてもよい。環境情報は、そのような環境情報を格納するように構成される１つ以上のサーバ、ハードドライブ、記憶装置及び／又はデータベース等のリソース（例えば、クラウドベースのリソース）にアップロードされてもよい。
【0267】
いくつかの実現例において、種々のイベントに関連付けられたチャネル状態情報は、リソース（例えば、クラウド内のデータベース）にアップロード可能であり、情報はチャネル状態情報固有シグニチャと既知のイベントとを関連付けるために使用可能である。いくつかの例において、イベント分類器をトレーニングする装置は、イベントに関連付けられた物体（例えば、扉及び窓）を有する間取り図を示すためにユーザインタフェースを提示するディスプレイを有してもよい。ユーザが扉をタップした場合、扉が閉鎖位置と開放位置との間で切り替わってもよい。扉が開放位置にあると仮定する。ユーザは、ボタンをタップしてトレーニングセッションを開始してもよく、それによりシステムはチャネル状態情報を収集して「開扉」イベントに関連付けられたパラメータを認識するようにイベント分類器をトレーニングすることを開始する。「開扉」イベントに対するトレーニングが完了した後、ユーザは扉をタップして扉を閉鎖位置に切り替えて、ボタンをタップしてトレーニングセッションを開始してもよく、それによりシステムはチャネル状態情報を収集して「閉扉」イベントに関連付けられたパラメータを認識するようにイベント分類器をトレーニングすることを開始する。ユーザは、他のイベントに対する分類器のトレーニングを同様の方法で開始できる。
【0268】
環境監視システムにおける送信機及び受信機の各々は、１つ以上のプロセッサ及び１つ以上のコンピュータ可読媒体（例えば、ＲＡＭ、ＲＯＭ、ＳＤＲＡＭ、ハードディスク、光ディスク及びフラッシュメモリ）を含むことができる。１つ以上のプロセッサは、上述した種々の計算を実行できる。計算は、特定用途向け集積回路（ＡＳＩＣ）を使用して同様に実現可能である。用語「コンピュータ可読媒体」は、実行するためにプロセッサに命令を提供する際に関わる不揮発性媒体（例えば、光ディスク又は磁気ディスク）及び揮発性媒体（例えば、メモリ）及び伝送媒体を含むがそれらに限定されない媒体を示す。伝送媒体は、同軸ケーブル、銅線及び光ファイバを含むが、それらに限定されない。
【0269】
上述した特徴は、データ格納システムとの間でデータ及び命令を送受信するように接続された少なくとも１つのプログラマブルプロセッサと、少なくとも１つの入力装置と、少なくとも１つの出力装置とを含むプログラマブルシステム上で実行可能な１つ以上のコンピュータプログラムにおいて有利に実現可能である。コンピュータプログラムは、特定のアクティビティを実行するため又は特定の結果を得るためにコンピュータにおいて直接又は間接的に使用可能な命令セットである。コンピュータプログラムは、コンパイラ型言語又はインタープリタ型言語を含むあらゆる形態のプログラミング言語（例えば、Ｃ、Ｊａｖａ）で記述可能であり、スタンドアロンプログラム又はモジュール、構成要素、サブルーチン、ブラウザベースウェブアプリケーション、あるいはコンピュータ環境における使用に適した他のユニットとしての形態を含むあらゆる形態で配布可能である。
【0270】
命令プログラムを実行するのに適したプロセッサは、例えばあらゆる種類のコンピュータの汎用マイクロプロセッサ及び専用マイクロプロセッサの双方、デジタル信号プロセッサ及びソロプロセッサ、あるいは複数のプロセッサ又はコアのうちの１つを含む。一般に、プロセッサは読み出し専用メモリ又はランダムアクセスメモリ又はその双方から命令及びデータを受信する。コンピュータの不可欠な要素は、命令を実行するプロセッサと命令及びデータを格納する１つ以上のメモリである。一般にコンピュータは、データファイルを格納するための１つ以上の大容量記憶装置を更に含むか又はそれらと通信するために動作可能に接続される。そのような装置は、内部ハードディスク及び取外し可能ディスク、光磁気ディスク及び光ディスク等の磁気ディスクを含む。コンピュータプログラム命令及びデータを具体的に実現するのに適した記憶装置は、例えばＥＰＲＯＭ、ＥＥＰＲＯＭ及びフラッシュメモリ素子等の半導体メモリ素子、内部ハードディスク及び取外し可能ディスク等の磁気ディスク、光磁気ディスク、並びにＣＤ－ＲＯＭ及びＤＶＤ－ＲＯＭディスクを含む全ての形態の不揮発性メモリを含む。プロセッサ及びメモリは、ＡＳＩＣ（特定用途向け集積回路）により補足されるか又はＡＳＩＣに組み込まれることが可能である。
【0271】
本明細書は多くの特定の実現例の詳細を含むが、それらは、あらゆる発明又は特許請求されてもよいものの範囲を制限すると解釈されるべきでなく、特定の発明の特定の実施形態に固有の特徴の説明であると解釈されるべきである。個別の実施形態において本明細書中で説明される特定の特徴は、単一の実施形態において組み合わされても同様に実現可能である。反対に、単一の実施形態において説明される種々の特徴は、複数の実施形態において個別に又は何らかの適切な部分的組み合わせでも同様に実現可能である。
【0272】
同様に、動作は図中で特定の順序で示されるが、所望の結果を達成するためにはそのような動作が示される特定の順序で又は順次実行されるか、あるいは図示される全ての動作が実行される必要があると理解されるべきではない。特定の状況において、マルチタスク動作及び並行処理が有利な場合がある。更に、上述した実施形態における種々のシステム構成要素の分離は、そのような分離が全ての実施形態で必要であると理解されるべきでなく、説明されたプログラム構成要素及びシステムは一般に、単一のソフトウェアに一体化可能であるか又は複数のソフトウェアにパッケージ化可能であると理解されるべきである。
【0273】
従って、主題の特定の実施形態を説明した。他の実施形態は以下の請求項の範囲内である。いくつかの例において、請求項において示される動作は異なる順序で実行可能であり、所望の結果を依然として達成できる。更に、添付の図面に示される処理は、所望の結果を達成するために、示される特定の順序又は順番を必ずしも必要としない。特定の実現例において、マルチタスク動作及び並行処理が有利な場合がある。
【0274】
本発明の多くの実現例を説明した。しかし、本発明の主旨の範囲から逸脱することなく種々の変更を行うことができることが理解されるだろう。例えばチャネルプロービング信号のフレーム構成、プリアンブル１６２、バーカ符号１６４、ペイロード１６６及びトレーニング系列１６８のフォーマットは上述したものと異なってもよい。送信機１４２及び受信機１４４は、図示されていない更に多くの構成要素を含むことができる。例えばアナログデジタル変換器、デジタルアナログ変換器及び電力増幅器等の標準的な構成要素は図面から省略されている。
【0275】
上記の例では１つの送信機及び１つの受信機が示されるが、複数の送信機及び複数の受信器が存在してもよい。監視システムの送信機又は受信機は、ロケーションが固定されてもよく又は移動できてもよい。監視システムの送信機又は受信機は、単の一アンテナ又は複数のアンテナを有することができる。時間反転監視システムは、人間に加え、猫、犬、馬又は鶏等の動物も監視できる。監視システムにより特定のイベントが検出された場合、システムは例えば音声、ビデオ、ショートメッセージサービス又は電子メール等で警報信号を与えることができる。システムが警報信号を送出する場合、例えば３Ｇ、４Ｇ－ＬＴＥ、５Ｇ、ＷｉＦｉ等を使用できる。上記の例ではチャネル状態情報はＷｉＦｉを介して取得されると示したが、これはＬＴＥ／ＨＳＰＡマクロセル、ピコセル、フェムトセル、ホットスポット等による信号であってもよい。監視システムの帯域幅は、例えば２０ＭＨｚであってもよく、あるいは異なるアンテナからのＣＳＩ、異なる周波数サブキャリア／帯域／チャネルにおけるＣＳＩ、異なる時間インスタンスにおけるＣＳＩ等を連結することにより２０ＭＨｚより小さくてもよい。時間反転監視システムはスタンドアロンであってもよく、誤警報を減少するために他の監視システムを補助するために使用されてもよく、あるいは感度を向上するために他の監視システムを補助するために使用されてもよい。時間反転監視システムはプログラム可能であってもよい。
【0276】
いくつかの例において、カスタマイズされた送信機及び受信機が、チャネルインパルス応答（時間領域）を取得するためにチャネルサウンディングにおいて使用可能である。いくつかの例において、市販のＷｉＦｉ装置をチャネルサウンディングに使用でき、チャネル推定値出力はチャネル周波数応答として読み取り可能である。受信機はチャネルを推定する市販のＷｉＦｉ装置であってもよく、チャネル周波数応答は受信機装置から抽出される。市販のＷｉＦｉ装置のデータ処理能力が制限されるため、コンピュータが分類器のトレーニングに関連する計算及び式の計算を実行する。いくつかの例において、ＷｉＦｉ受信機がパワーデータプロセッサを有する場合、分類器のトレーニングに関連する計算は受信機により実行可能であり、全ての式は受信機において計算される。
【0277】
チャネル状態情報ｈ[１]、ｈ[２]、ｈ[Ｌ]等に対するサンプリング周波数は、システムの帯域幅に従って選択可能である。例えば、システムの帯域幅（両側）に少なくとも等しい。サンプリング周波数は高いほどよい。タップの数Ｌは大きいほどよい。いくつかの例において、収集されたチャネルのエネルギーが最大になるように、Ｌを徐々に増加することができる。例えばＬ＝３０とＬ＝５０の双方で同様のチャネル状態情報が得られる場合、Ｌ＝３０を使用できる。
【0278】
プロービング信号のフォーマットは図２７及び図２９に示したものと異なってもよく、ハンドシェイク処理は上述したものと異なってもよい。例えば無線信号を送受信するために市販のＷｉＦｉ装置が使用される場合、ハンドシェイク処理はＷｉＦｉプロトコルに従う。
【0279】
イベント検出システム又は人間動作検出システムは、多くの方法で実現可能である。例えば図３８を参照すると、イベント検出システム又は人間動作検出システム６３０は会場６３８内に配置された送信機６３２及び受信機６３６を含むことができ、送信機６３２は、システム６３０により検出されるべきイベント又は人間の動作により影響を受けるマルチパスチャネル６３４を介して受信機６３６に無線信号を送出する。この例において、受信機６３６は、上述した種々のデータ処理ステップを実行するためのデータプロセッサ６４０を含む。
【0280】
いくつかの実現例において、図３９を参照すると、イベント検出システム又は人間動作検出システム６５０は、送信機６３２、受信機６３６及びコンピュータ６５２を含むことができる。受信機６３６は基本的なデータ処理を実行するデータプロセッサ６４０を含み、複雑な計算の大部分は受信機６３６と通信するコンピュータ６５２により実行される。
【0281】
いくつかの実現例において、図４０を参照すると、イベント検出システム又は人間動作検出システム６６０は送信機６３２及び受信機６３６を含むことができ、受信機６３６は、基本的なデータ処理を実行するデータプロセッサ６４０を含む。受信機６３６はクラウドサーバ６６２と通信し、クラウドサーバ６６２が複雑な計算の大部分を実行する。
【0282】
従って、他の実施形態は以下の請求項の範囲内である。
【図1A】
【図1B】
【図1C】
【図2A】
【図2B】
【図3】
【図4A】
【図4B】
【図5A】
【図5B】
【図6A】
【図6B】
【図6C】
【図7A】
【図7B】
【図7C】
【図8A】
【図8B】
【図9A】
【図9B】
【図9C】
【図10A】
【図10B】
【図10C】
【図10D】
【図10E】
【図10F】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23A】
【図23B】
【図23C】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】
【図33】
【図34】
【図35】
【図36】
【図37】
【図38】
【図39】
【図40】