電磁波の伝搬路に前記電磁波を反射又は吸収する物体が無い解析空間において、所定の送信位置から電磁波が送信された際に、予め定められた受信位置における電界及び磁界の時間応答をシミュレーションする時間応答算出部と、電磁波の伝搬路に前記物体が存在する解析空間において、前記送信位置から電磁波が送信された際に、前記受信位置における電界及び磁界の時間応答をシミュレーションする伝搬路応答算出部と、前記時間応答算出部によって算出された電界の時間応答及び前記伝搬路応答算出部によって算出された電界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記時間応答算出部によって算出された磁界の時間応答及び前記伝搬路応答算出部によって算出された磁界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻判定部と、を備える到来時刻推定装置。前記到来時刻判定部は、前記伝搬路応答算出部によって算出された前記電界及び磁界の時間応答それぞれについて電界強度及び磁界強度のピークを検出するピーク検出部と、前記ピーク検出部によって検出された電界強度のピーク及び磁界強度のピークと、それぞれが検出された時刻とに基づいて、前記時間応答算出部によって算出された前記電界及び磁界の時間応答それぞれを前記ピーク部分の波形と近い形状となるように複数の態様で変形し、電界と磁界それぞれについて複数の変形応答を算出する変形応答算出部と、前記複数の変形応答それぞれについて、前記ピーク部分の波形との類似度を表す評価値を算出し、評価値に基づいて電界と磁界それぞれについて前記変形応答を選択する評価値算出部と、電界について選択された前記変形応答に基づいて電界の時間応答に基づいた前記時刻を推定し、磁界について選択された前記変形応答に基づいて磁界の時間応答に基づいた前記時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻推定部と、を備える、請求項1に記載の到来時刻推定装置。前記到来時刻判定部は、最初の処理においては、前記伝搬路応答算出部によって算出された時間応答から、前記選択された変形応答を減算することによって新たな時間応答を算出し、2回目以降の処理においては自身が直前に算出した時間応答から、前記選択された変形応答を減算することによってさらに新たな時間応答を算出する波形処理部と、前記波形処理部によって算出された新たな時間応答に基づいて繰り返し処理の終了を判定する終了判定部と、を備え、前記終了判定部が終了しないと判定した場合は、前記到来時刻判定部の各部は処理を繰り返し実行し、前記終了判定部が終了すると判定した場合は、前記到来時刻判定部の各部は処理を終了する、請求項2に記載の到来時刻推定装置。コンピュータを、請求項1~3のいずれか一項に記載の到来時刻推定装置として動作させるためのコンピュータプログラム。電磁波の伝搬路に前記電磁波を反射又は吸収する物体が無い解析空間において、所定の送信位置から電磁波が送信された際に、予め定められた受信位置における電界及び磁界の時間応答をシミュレーションする時間応答算出ステップと、電磁波の伝搬路に前記物体が存在する解析空間において、前記送信位置から電磁波が送信された際に、前記受信位置における電界及び磁界の時間応答をシミュレーションする伝搬路応答算出ステップと、前記時間応答算出ステップによって算出された電界の時間応答及び前記伝搬路応答算出ステップによって算出された電界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記時間応答算出ステップによって算出された磁界の時間応答及び前記伝搬路応答算出ステップによって算出された磁界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻判定ステップと、を有する到来時刻推定方法。前記到来時刻判定ステップは、前記伝搬路応答算出ステップによって算出された前記電界及び磁界の時間応答それぞれについて電界強度及び磁界強度のピークを検出するピーク検出ステップと、前記ピーク検出ステップによって検出された電界強度のピーク及び磁界強度のピークと、それぞれが検出された時刻とに基づいて、前記時間応答算出ステップによって算出された前記電界及び磁界の時間応答それぞれを前記ピーク部分の波形と近い形状となるように複数の態様で変形し、電界と磁界それぞれについて複数の変形応答を算出する変形応答算出ステップと、前記複数の変形応答それぞれについて、前記ピーク部分の波形との類似度を表す評価値を算出し、評価値に基づいて電界と磁界それぞれについて前記変形応答を選択する評価値算出ステップと、電界について選択された前記変形応答に基づいて電界の時間応答に基づいた前記時刻を推定し、磁界について選択された前記変形応答に基づいて磁界の時間応答に基づいた前記時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻推定ステップと、を有する、請求項5に記載の到来時刻推定方法。前記到来時刻判定ステップは、最初の処理においては、前記伝搬路応答算出ステップによって算出された時間応答から、前記選択された変形応答を減算することによって新たな時間応答を算出し、2回目以降の処理においては自身が直前に算出した時間応答から、前記選択された変形応答を減算することによってさらに新たな時間応答を算出する波形処理ステップと、前記波形処理ステップによって算出された新たな時間応答に基づいて繰り返し処理の終了を判定する終了判定ステップと、をさらに有し、前記終了判定ステップが終了しないと判定した場合は、前記到来時刻判定ステップの各ステップの処理を繰り返し実行し、前記終了判定ステップが終了すると判定した場合は、前記到来時刻判定ステップの各処理を終了する、請求項6に記載の到来時刻推定方法。FDTD法では、主にガウシアンパルスを用いて解析が実施される。ガウシアンパルスは波形が時間的な広がりを有しているため、到来時刻の推定精度が悪くなってしまう。さらに、複数の到来波が存在する環境においては、時間的な広がりを持った到来波が合成された状態で受信位置において観測される。そのため、このような原因によっても到来時刻の推定精度が悪くなってしまう。また、このように到来時刻の推定精度が悪くなってしまう問題は、ガウシアンパルスに限らず時間的な広がりを有した電磁波を用いた解析に共通する問題である。上記事情に鑑み、本発明は、時間的な広がりを有した電磁波の到来時刻を精度良く推定する技術を提供することを目的としている。本発明の一態様は、到来時刻推定装置であって、電磁波の伝搬路に前記電磁波を反射又は吸収する物体が無い解析空間において、所定の送信位置から電磁波が送信された際に、予め定められた受信位置における電界及び磁界の時間応答をシミュレーションする時間応答算出部と、電磁波の伝搬路に前記物体が存在する解析空間において、前記送信位置から電磁波が送信された際に、前記受信位置における電界及び磁界の時間応答をシミュレーションする伝搬路応答算出部と、前記時間応答算出部によって算出された電界の時間応答及び前記伝搬路応答算出部によって算出された電界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記時間応答算出部によって算出された磁界の時間応答及び前記伝搬路応答算出部によって算出された磁界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻判定部と、を備える。本発明の一態様は、上記の到来時刻推定装置であって、前記到来時刻判定部は、前記伝搬路応答算出部によって算出された前記電界及び磁界の時間応答それぞれについて電界強度及び磁界強度のピークを検出するピーク検出部と、前記ピーク検出部によって検出された電界強度のピーク及び磁界強度のピークと、それぞれが検出された時刻とに基づいて、前記時間応答算出部によって算出された前記電界及び磁界の時間応答それぞれを前記ピーク部分の波形と近い形状となるように複数の態様で変形し、電界と磁界それぞれについて複数の変形応答を算出する変形応答算出部と、前記複数の変形応答それぞれについて、前記ピーク部分の波形との類似度を表す評価値を算出し、評価値に基づいて電界と磁界それぞれについて前記変形応答を選択する評価値算出部と、電界について選択された前記変形応答に基づいて電界の時間応答に基づいた前記時刻を推定し、磁界について選択された前記変形応答に基づいて磁界の時間応答に基づいた前記時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻推定部と、を備える。本発明の一態様は、上記の到来時刻推定装置であって、前記到来時刻判定部は、最初の処理においては、前記伝搬路応答算出部によって算出された時間応答から、前記選択された変形応答を減算することによって新たな時間応答を算出し、2回目以降の処理においては自身が直前に算出した時間応答から、前記選択された変形応答を減算することによってさらに新たな時間応答を算出する波形処理部と、前記波形処理部によって算出された新たな時間応答に基づいて繰り返し処理の終了を判定する終了判定部と、を備え、前記終了判定部が終了しないと判定した場合は、前記到来時刻判定部の各部は処理を繰り返し実行し、前記終了判定部が終了すると判定した場合は、前記到来時刻判定部の各部は処理を終了する。本発明の一態様は、コンピュータを、上記の到来時刻推定装置として動作させるためのコンピュータプログラムである。本発明の一態様は、到来時刻判定方法であって、電磁波の伝搬路に前記電磁波を反射又は吸収する物体が無い解析空間において、所定の送信位置から電磁波が送信された際に、予め定められた受信位置における電界及び磁界の時間応答をシミュレーションする時間応答算出ステップと、電磁波の伝搬路に前記物体が存在する解析空間において、前記送信位置から電磁波が送信された際に、前記受信位置における電界及び磁界の時間応答をシミュレーションする伝搬路応答算出ステップと、前記時間応答算出ステップによって算出された電界の時間応答及び前記伝搬路応答算出ステップによって算出された電界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記時間応答算出ステップによって算出された磁界の時間応答及び前記伝搬路応答算出ステップによって算出された磁界の時間応答に基づいて、前記受信位置に前記電磁波が到来した時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻判定ステップと、を有する。本発明の一態様は、上記の到来時刻判定方法であって、前記到来時刻判定ステップは、前記伝搬路応答算出ステップによって算出された前記電界及び磁界の時間応答それぞれについて電界強度及び磁界強度のピークを検出するピーク検出ステップと、前記ピーク検出ステップによって検出された電界強度のピーク及び磁界強度のピークと、それぞれが検出された時刻とに基づいて、前記時間応答算出ステップによって算出された前記電界及び磁界の時間応答それぞれを前記ピーク部分の波形と近い形状となるように複数の態様で変形し、電界と磁界それぞれについて複数の変形応答を算出する変形応答算出ステップと、前記複数の変形応答それぞれについて、前記ピーク部分の波形との類似度を表す評価値を算出し、評価値に基づいて電界と磁界それぞれについて前記変形応答を選択する評価値算出ステップと、電界について選択された前記変形応答に基づいて電界の時間応答に基づいた前記時刻を推定し、磁界について選択された前記変形応答に基づいて磁界の時間応答に基づいた前記時刻を推定し、前記二つの推定結果に基づいて最終的に前記時刻を推定する到来時刻推定ステップと、を有する。本発明の一態様は、到来時刻判定方法であって、前記到来時刻判定ステップは、最初の処理においては、前記伝搬路応答算出ステップによって算出された時間応答から、前記選択された変形応答を減算することによって新たな時間応答を算出し、2回目以降の処理においては自身が直前に算出した時間応答から、前記選択された変形応答を減算することによってさらに新たな時間応答を算出する波形処理ステップと、前記波形処理ステップによって算出された新たな時間応答に基づいて繰り返し処理の終了を判定する終了判定ステップと、をさらに有し、前記終了判定ステップが終了しないと判定した場合は、前記到来時刻判定ステップの各ステップの処理を繰り返し実行し、前記終了判定ステップが終了すると判定した場合は、前記到来時刻判定ステップの各処理を終了する。本発明により、時間的な広がりを有した電磁波の到来時刻を精度良く推定することが可能となる。[概略]FDTD法等の電界および磁界についての解析が可能な電磁界解析を用いて、仮想の送信位置から送信した電磁波が、伝送路に物体が無い環境以下、理想環境という。において受信位置に到来する際の電界および磁界のそれぞれの時間応答をシミュレーションにより算出し保持する。以下、このような理想環境におけるシミュレーションによって算出される時間応答を理想時間応答という。次に、伝搬路に物体が配置された環境以下、反射環境という。において、同様にシミュレーションにより電界および磁界のそれぞれの時間応答を算出し保持する。以下、このような反射環境におけるシミュレーションによって算出される時間応答を反射時間応答という。なお、いずれのシミュレーションにおいても、送信される電磁波は、例えばガウシアンパルスのように時間的な広がりを有している。反射環境では、電磁波が伝送路において反射するため、受信位置に電磁波が到来するまでの時間に遅延が生じる。また、到来する電磁波の強度も理想環境の場合と異なる。この遅延及び強度の変化に基づいて理想時間応答の波形を変形し、反射時間応答の波形から連続して減算する。その電界および磁界それぞれの時間応答に基づいて、到来時刻を推定する。電界および磁界それぞれの時間応答を用いて到来時刻を推定することで、高精度な到来時刻の推定が可能になる。[実施形態]図1は、電磁波の到来時刻を推定する技術の概要を示す図である。仮想空間40内には、送信位置10と受信位置20とが配置される。図1の例では、仮想空間40は四方を壁で囲まれている。仮想空間内には1又は複数の物体30が配置される。物体30は、電磁波を反射したり吸収したりする。送信位置10からは所定の電磁波が送信され、壁や物体30上で電磁波が反射し、受信位置20に電磁波が到来する。受信位置20に電磁波が到来する時刻を推定する際には、受信位置20における電界の強さ電界強度と磁界の強さ磁界強度とがシミュレートされる。到来時刻推定装置は、受信位置20における電界強度及び磁界強度の時間変化に基づいて、送信位置10から送信された電磁波が受信位置20に到来する時刻以下、到来時刻という。を判定する。図2は、到来時刻推定装置500の構成を示す概略ブロック図である。到来時刻推定装置500は、バスで接続されたCPUCentralProcessingUnitやメモリや補助記憶装置などを備え、到来時刻推定プログラムを実行する。到来時刻推定プログラムの実行によって、到来時刻推定装置500は、パラメータ記憶部501、時間応答記憶部502、伝搬路応答記憶部503、到来時刻記憶部504、パラメータ設定部511、初期化部512、時間応答算出部520、伝搬路応答算出部530、到来時刻判定部540を備える装置として機能する。なお、到来時刻推定装置500の各機能の全て又は一部は、ASICApplicationSpecificIntegratedCircuitやPLDProgrammableLogicDeviceやFPGAFieldProgrammableGateArray等のハードウェアを用いて実現されても良い。到来時刻推定プログラムは、コンピュータ読み取り可能な記録媒体に記録されても良い。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM、CD-ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。パラメータ記憶部501は、電磁波の伝搬をシミュレートする際の各種のパラメータを記憶する。時間応答記憶部502は、理想環境においてシミュレーションを行った結果として、各時刻における電界強度を表すテーブル理想電界テーブルと、各時刻における磁界強度を表すテーブル理想磁界テーブルとを記憶する。伝搬路応答記憶部503は、反射環境においてシミュレーションを行った結果として、各時刻における電界強度を表すテーブル反射電界テーブルと、各時刻における磁界強度を表すテーブル反射磁界テーブルとを記憶する。到来時刻記憶部504は、到来時刻判定部540によって推定された到来時刻と、その到来時刻を表す識別情報とを対応付けたテーブル到来時刻テーブルを記憶する。パラメータ設定部511は、シミュレーションにおいて送信位置10から送信された電磁波が伝搬する空間、すなわちシミュレーション対象となる空間である解析空間のサイズ、図1に示すように解析空間内に設けられた閉空間の境界を構成する壁の位置、壁の電磁波の反射率、解析空間内に配置されている物体30の位置、サイズ、電磁波の反射率などのパラメータの入力を受け付け、受け付けたパラメータをパラメータ記憶部501に記録する。また、パラメータ設定部511は、解析空間における送信位置10の座標、受信位置20の座標、送信位置10から送信される電磁波を定義する情報、理想環境においてシミュレーションを行う際の離散時間間隔、反射環境においてシミュレーションを行う際の離散時間間隔、シミュレーションを行う期間などのパラメータの入力を受け付け、受け付けたパラメータをパラメータ記憶部501に記録する。初期化部512は、時間応答記憶部502、伝搬路応答記憶部503、到来時刻記憶部504の記憶領域に、予め定められた初期値を記憶させる。この処理によって、各記憶領域は初期化される。初期値は、例えば、零や無効値を示す値である。時間応答算出部520は、パラメータ記憶部501に記憶されている各パラメータを用いて、理想環境において送信位置10から電磁波が送信された場合における電磁波の伝搬のシミュレーションを行う。時間応答算出部520は、シミュレーションにおいて、受信位置20の電界強度を離散時間間隔ごとに算出し、算出した電界強度を時系列で理想電界テーブルに記憶させる。すなわち、時間応答算出部520は、シミュレーション対象となる解析空間における電磁波の伝搬を理想環境でシミュレートする。そして、時間応答算出部520は、受信位置20における電界の時間応答を表す理想電界テーブルを時間応答記憶部502に記憶させる。また、時間応答算出部520は、シミュレーションにおいて、受信位置20の磁界強度を離散時間間隔ごとに算出し、算出した磁界強度を時系列で理想磁界テーブルに記憶させる。すなわち、時間応答算出部520は、シミュレーション対象となる解析空間における電磁波の伝搬を理想環境でシミュレートする。そして、時間応答算出部520は、受信位置20における磁界の時間応答を表す理想磁界テーブルを時間応答記憶部502に記憶させる。伝搬路応答算出部530は、パラメータ記憶部501に記憶されている各パラメータを用いて、反射環境において送信位置10から電磁波が送信された場合における電磁波の伝搬のシミュレーションを行う。伝搬路応答算出部530は、シミュレーションにおいて、受信位置20の電界強度を離散時間間隔ごとに算出し、算出した電界強度を時系列で反射電界テーブルに記憶させる。すなわち、伝搬路応答算出部530は、シミュレーション対象となる解析空間における電磁波の伝搬を反射環境でシミュレートする。そして、伝搬路応答算出部530は、受信位置20における電界の時間応答を表す反射電界テーブルを伝搬路応答記憶部503に記憶させる。また、伝搬路応答算出部530は、シミュレーションにおいて、受信位置20の磁界強度を離散時間間隔ごとに算出し、算出した磁界強度を時系列で反射磁界テーブルに記憶させる。すなわち、伝搬路応答算出部530は、シミュレーション対象となる解析空間における電磁波の伝搬を反射環境でシミュレートする。そして、伝搬路応答算出部530は、受信位置20における磁界の時間応答を表す反射磁界テーブルを伝搬路応答記憶部503に記憶させる。図3は、時間応答算出部520及び伝搬路応答算出部530のシミュレーションによって得られる理想時間応答及び反射時間応答の具体例を表す図である。図3Aは、時間応答算出部520のシミュレーションによって得られる理想時間応答のグラフである。図3Bは、時間応答算出部520のシミュレーションによって得られる離散的な値を表すグラフである。図3Cは、伝搬路応答算出部530のシミュレーションによって得られる反射時間応答のグラフである。図3Dは、伝搬路応答算出部530のシミュレーションによって得られる離散的な値を表すグラフである。時間応答算出部520は、伝搬路に物体30が配置されていない仮想空間においてシミュレーションを行う。そのため、図3Aに示されるように理想時間応答の波形は乱れが少ない。一方、伝搬路応答算出部530は、伝搬路に物体30が配置されている仮想空間においてシミュレーションを行う。そのため、図3Cに示されるように反射時間応答の波形には乱れが生じる。したがって、反射時間応答の波形は、複数のピークを持った複雑な波形となる。図3B及び図3Dには、標本点60を示している。時間応答算出部520及び伝搬路応答算出部530は、所定の標本化周期サンプリング周期で電界の強度を測定し、その結果をテーブルとして記録する。そのため、テーブルには、各標本点60の電界強度Eとその測定時刻Tとを対応付けたレコード70が、標本点60毎に記録される。なお、図3は電界に関する時間応答を示すグラフであり、以上の説明も電界に関するものである。しかしながら、図3を用いて説明した事項は、理想時間応答と反射時間応答との違いに関するものであり、磁界に関しても同様のことがいえる。図4は、理想電界テーブル及び理想磁界テーブルの概略を示す図である。図4Aは、理想電界テーブルの概略を示す。図4Bは、理想磁界テーブルの概略を示す。理想電界テーブルは、時刻及び電界強度を対応付けたレコード70を複数有する。理想磁界テーブルは、時刻及び磁界強度を対応付けたレコード70を複数有する。レコード70は、図3B及び図3Dに示される標本点60に関する値を表す。なお、反射電界テーブル及び反射磁界テーブルの構成は、レコード70に記録される値が理想電界テーブル及び理想磁界テーブルと異なるものの、その構成の概略は同じであり図4に示す通りである。図2に戻って、到来時刻推定装置500の構成についての説明を続ける。到来時刻判定部540は、時間応答記憶部502に記憶されている情報と、伝搬路応答記憶部503に記憶されている情報とに基づいて、送信位置10から送信された電磁波が受信位置20に到達した時刻を表す情報を取得する。到来時刻判定部540は、ピーク検出部541、変形応答算出部542、評価値算出部543、到来時刻推定部544、波形処理部545、終了判定部546を備えている。ピーク検出部541は、各テーブルからピークを検出し、ピーク時の電界強度及び磁界強度とその時刻とを取得する。具体的には以下の通りである。ピーク検出部541は、時間応答記憶部502に記憶されている理想電界テーブルから、電界強度絶対値の最大値ei-maxと、このときの時刻tei-maxと、を取得する。また、ピーク検出部541は、時間応答記憶部502に記憶されている理想磁界テーブルから、磁界強度絶対値の最大値hi-maxと、このときの時刻thi-maxと、を取得する。さらに、ピーク検出部541は、伝搬路応答記憶部503に記憶されている反射電界テーブルから、電界強度絶対値の最大値er-maxと、このときの時刻ter-maxと、を取得する。また、ピーク検出部541は、伝搬路応答記憶部503に記憶されている反射磁界テーブルから、磁界強度絶対値の最大値hr-maxと、このときの時刻thr-maxと、を取得する。変形応答算出部542は、理想電界テーブルによって表される波形を、ピーク検出部541によって検出されたピーク部分周辺の波形と近い形状となるように変形する。具体的には、変形応答算出部542は、時間応答記憶部502に記憶されている理想電界テーブルによって表される波形を、ピーク検出部541によって検出された時刻ter-max及び電界強度er-maxに基づいて複数の態様に変形し、複数の電界変形応答を算出する。また、変形応答算出部542は、時間応答記憶部502に記憶されている理想磁界テーブルによって表される波形を、ピーク検出部541によって検出された時刻thr-max及び磁界強度hr-maxに基づいて複数の態様に変形し、複数の磁界変形応答を算出する。図5は、変形応答算出部542の処理の概略を示す図である。図5において、電界強度の軸縦軸に沿って伸びる複数の破線は、それぞれシミュレーションにおける離散時間間隔を表す。図5A~図5Fは、電界の反射時間応答のグラフのうち、ピーク検出部541によって検出された電界の最大値er-max付近を拡大して示した図である。変形応答算出部542は、まず、複数の候補標準点61の時刻ter-i:i=0~n:nは0以上の整数及び電界強度er-iを取得する。候補標準点61は、ピーク検出部541によって検出された電界強度の最大値を有する標準点60以下、ピーク標準点という。と、その周辺の標準点60とを含む。図5は、n=2の場合の例である。図5では、i=0の候補標準点61はピーク標準点であり、このときの時刻及び電界強度はそれぞれter-0及びer-0である。なおter-0=ter-maxであり、er-0=er-maxである。i=1の候補標準点61は、ピーク標準点よりも時刻が一つ前の標準点であり、このときの時刻及び電界強度はそれぞれter-1及びer-1である。i=2の候補標準点61は、ピーク標準点よりも時刻が一つ後の標準点であり、このときの時刻及び電界強度はそれぞれter-2及びer-2である。次に、変形応答算出部542は、各候補標準点61の時刻及び電界強度に応じて、電界の理想時間応答を変形する。具体的には、変形応答算出部542は以下の処理を候補標準点61毎に行う。まず、変形応答算出部542は時刻のオフセット量Oiを算出する。オフセット量Oiは、例えば式1を用いて算出できる。また、変形応答算出部542は、電界強度の係数Fiを算出する。係数Fiは、例えば式2を用いて算出できる。係数Fiは、電界の理想時間応答の波形を、その電界強度の最大値が候補標準点61の電界強度er-iとなるように電界強度軸に沿って拡大又は縮小するための係数である。次に、変形応答算出部542は、理想電界テーブルの各レコードの時刻に対してオフセット量Oiを加算し、各レコードの電界強度に対して係数Fiを乗算する。この処理の結果として得られるテーブルが、一つの電界変形応答のテーブルである。図6は、電界変形応答のテーブルの概略を示す図である。変形応答算出部542は、このような処理を候補標準点61毎に行うことによって、候補標準点61毎に電界変形応答のテーブルを算出する。また、変形応答算出部542は、同様の処理を磁界についても行うことによって、複数の磁界変形応答のテーブルを算出する。図2に戻って、到来時刻推定装置500の構成についての説明を続ける。評価値算出部543は、変形応答算出部542によって算出された複数の電界変形応答毎に評価値を算出する。そして、評価値算出部543は、算出された評価値に基づいて、複数の電界変形応答の中から一つの電界変形応答を選択する。また、評価値算出部543は、変形応答算出部542によって算出された複数の磁界変形応答毎に評価値を算出する。そして、評価値算出部543は、算出された評価値に基づいて、複数の磁界変形応答の中から一つの磁界変形応答を選択する。評価値算出部543が電界変形応答に関して算出する評価値は、各電界変形応答の波形が、伝搬路応答算出部530に記憶されている電界の反射時間応答の波形にどれだけ類似しているかを表す値である。同様に、評価値算出部543が磁界変形応答に関して算出する評価値は、各磁界変形応答の波形が、伝搬路応答算出部530に記憶されている磁界の反射時間応答の波形にどれだけ類似しているかを表す値である。そして、評価値算出部543は、類似度の高い例えば、最も類似している電界変形応答及び磁界変形応答を一つずつ選択する。例えば、評価値算出部543は以下のように評価値を算出しても良い。評価値算出部543は、電界変形応答のテーブル毎に、伝搬路応答記憶部503に記憶されている反射電界テーブルとの二乗誤差を算出し、算出した二乗誤差を評価値とする。ここで、二乗誤差とは、各時刻における電界強度の差の二乗和である。例えば、ある時刻における電界変形応答の電界強度から、反射電界テーブルにおいて同じ時刻に対応付けられている電界強度を減算し、減算結果を二乗した値の全時刻の総和が評価値である。評価値算出部543は、二乗誤差の値が最も小さい電界変形応答を選択する。磁界についても同様である。到来時刻推定部544は、変形応答算出部542によって算出された電界に関する候補標準点61の時刻に基づいて、電界に関する処理に基づいて得られる到来時刻以下、電界到来時刻という。を推定する。また、到来時刻推定部544は、変形応答算出部542によって算出された磁界に関する候補標準点61の時刻に基づいて、磁界に関する処理に基づいて得られる到来時刻以下、磁界到来時刻という。を推定する。そして、到来時刻推定部544は、電界到来時刻及び磁界到来時刻に基づいて、最終的に到来時刻を推定する。以下、到来時刻の推定方法の具体例について説明する。到来時刻推定部544は、電界に関する候補標準点61のうち最も評価値が小さかった最も類似していると評価された候補標準点61の時刻以下、最適電界時刻という。を、電界到来時刻として推定しても良い。この場合、到来時刻推定部544は、磁界に関する候補標準点61のうち最も評価値が小さかった最も類似していると評価された候補標準点61の時刻以下、最適磁界時刻という。を、磁界到来時刻として推定しても良い。そして、到来時刻推定部544は、電界到来時刻と磁界到来時刻の統計値例えば平均値を、最終的に到来時刻として推定しても良い。また、到来時刻推定部544は、電界到来時刻及び磁界到来時刻のうち、評価値がより小さかったより類似していると評価された時刻を、到来時刻として推定しても良い。到来時刻推定部544は、推定した到来時刻と、この到来時刻を表す識別情報とを対応付けて到来時刻記憶部504に記録する。図7は、到来時刻記憶部504に記録されるテーブルの具体例を示す図である。図7Aは、電界及び磁界の双方に基づいて推定された到来時刻のテーブルを示す図である。到来時刻推定部544は、上述したように推定した到来時刻と識別情報とを対応付けたテーブルを、図7Aのように到来時刻記憶部504に記録する。このとき、到来時刻推定部544は、識別情報として、到来時刻判定部540が後述するように処理を繰り返した回数を用いても良い。この場合、最初に到来時刻を推定した場合には、到来時刻推定部544は識別情報として1を用い、その後新たに到来時刻を推定する度に1ずつ増やした値を識別情報として用いても良い。図7Bは、電界到来時刻のテーブルを示す図である。到来時刻推定部544は、図7Aに示されるようなテーブルのみならず、電界到来時刻のテーブルや、磁界到来時刻のテーブルをさらに記録しても良い。図7Bにおける到来時刻及び電界強度は、例えば最適電界時刻及びその時の電界強度を示す。図7Bにおける識別情報は、図7Aにおける識別情報と同じ情報である。したがって、図7Aのテーブル及び図7Bのテーブルにおいて、同じ識別情報の到来時刻は、同じタイミングで時刻推定部544によって記録された値である。波形処理部545は、伝搬路応答記憶部503に記憶されている反射電界テーブルによって表される波形から、評価値算出部543によって選択された電界変形応答の波形を減算する。また、波形処理部545は、伝搬路応答記憶部503に記憶されている反射磁界テーブルによって表される波形から、評価値算出部543によって選択された磁界変形応答の波形を減算する。具体的には、波形処理部545は、反射電界テーブルの全ての時刻において、各時刻の電界強度から、電界変形応答のテーブルの同じ時刻の電界強度の値を減算する。また、波形処理部545は、反射磁界テーブルの全ての時刻において、各時刻の磁界強度から、磁界変形応答のテーブルの同じ時刻の磁界強度の値を減算する。図8は、波形処理部545の処理の概略を示す図である。図8Aは、反射磁界テーブルによって表される波形の具体例を示す図である。図8Bは、波形処理部545による減算処理が行われた結果の波形の具体例を示す図である。図8Bにおいて破線で示される部分は減算前の波形を表す。波形処理部545の減算処理によって、ピーク検出部541によって検出されたピーク部分標準点60-1が無くなり、それまで2番目に大きかった標準点60-2の電界強度が1番大きいピーク部分になる。波形処理部545は、減算結果のテーブルを、伝搬路応答記憶部503に記録する。図2に戻って、到来時刻推定装置500の構成についての説明を続ける。終了判定部546は、予め設定された終了条件を満たすか否かについて、波形処理部545の最新の処理の結果のテーブルに基づいて判定する。終了条件を満たす場合には、終了判定部546は、到来時刻判定部540の処理を終了させる。一方、終了条件を満たさない場合は、終了判定部546は、到来時刻判定部540の処理を終了させず、到来時刻判定部540は波形処理部545によって算出された結果を新たな反射時間応答として用いることによって、上述した処理を繰り返し実行する。この場合、例えば図8Bを例にとると、次の処理ではピーク検出部541によって標準点60-2がピークとして検出され、標準点60-2付近の時刻が到来時刻として新たに推定され、波形処理部545によって標準点60-2付近の波形が減算される。終了条件が満たされるまで、次々と大きなピークに関するこのような処理が繰り返し実行される。終了条件は、例えば最新の処理の結果のテーブルにおける電界強度又は磁界強度の最大値が、予め設定された電界閾値又は磁界閾値を下回ることであっても良い。また、終了条件は、例えば最新の処理の結果のテーブルに基づいて算出される電界強度の積分値又は磁界強度の積分値が、予め設定された電界閾値又は磁界閾値を下回ることであっても良い。また、終了条件は、到来時刻判定部540が所定の回数処理を繰り返したことであっても良い。終了条件はその他の条件であっても良い。図9は、到来時刻推定装置500の処理を示すフローチャートである。到来時刻の推定処理が開始されると、パラメータ設定部511が、理想環境に関する各種パラメータをパラメータ記憶部501に記憶させるステップS101。次に、初期化部512が、時間応答記憶部502に記憶される内容を初期化するステップS102。次に、時間応答算出部520は、パラメータ記憶部501に記憶されているパラメータに基づいて、送信位置10から電磁波が送信された時刻を基準例えば、時刻0として、シミュレーション対象に定められている解析空間内における電界成分及び磁界成分を算出する。そして、時間応答算出部520は、受信位置20における電界及び磁界の時間応答理想時間応答を、時間応答記憶部502に記憶させるステップS103。また、時間応答算出部520は、吸収境界条件例えば、Murの吸収境界条件や、BerengerのPML吸収境界条件を用いて、解析空間の境界に達した電磁波が反射しないようにする演算を行うステップS104。時間応答算出部520は、シミュレーションにおける所定の期間に亘って、受信位置20における電界強度及び磁界強度が終了閾値以下であるか否かを判定するステップS105。受信位置20における電界強度及び磁界強度が終了閾値以下である場合ステップS105-Yes、時間応答算出部520は、受信位置20に新たな電磁波が到来しないとみなして、電磁波の伝搬のシミュレーション時間応答計算処理を終了する。一方、受信位置20における電界強度及び磁界強度が終了閾値を超えている場合ステップS105-No、時間応答算出部520は、処理をステップS103に戻してシミュレーション時間応答計算処理を続行する。時間応答計算処理が終了すると、パラメータ設定部511が、反射環境に関する各種パラメータをパラメータ記憶部501に記憶させるステップS106。次に、初期化部512が、伝搬路応答記憶部503、到来時刻記憶部504に記憶される内容を初期化するステップS107。なお、伝搬路応答記憶部503、到来時刻記憶部504の初期化の処理は、ステップS102の際にまとめて行われても良い。次に、伝搬路応答算出部530は、パラメータ記憶部501に記憶されているパラメータに基づいて、送信位置10から電磁波が送信された時刻を基準例えば、時刻0として、シミュレーション対象に定められている解析空間内における電界成分及び磁界成分を算出する。そして、伝搬路応答算出部530は、受信位置20における電界及び磁界の時間応答反射時間応答を、伝搬路応答記憶部503に記憶させるステップS108。また、伝搬路応答算出部530は、吸収境界条件例えば、Murの吸収境界条件や、BerengerのPML吸収境界条件を用いて、解析空間の境界に達した電磁波が反射しないようにする演算を行うステップS109。伝搬路応答算出部530は、シミュレーションにおける所定の期間に亘って、受信位置20における電界強度及び磁界強度が終了閾値以下であるか否かを判定するステップS110。受信位置20における電界強度及び磁界強度が終了閾値以下である場合ステップS110-Yes、伝搬路応答算出部530は、受信位置20に新たな電磁波が到来しないとみなして、電磁波の伝搬のシミュレーションを終了する。一方、受信位置20における電界強度及び磁界強度が終了閾値を超えている場合ステップS110-No、伝搬路応答算出部530は、処理をステップS108に戻してシミュレーションを続行する。次に、ピーク検出部541がピークを検出し、ピーク時の電界強度及び磁界強度とその時刻とを取得するステップS111。次に、変形応答算出部542が、複数の電界変形応答と複数の磁界変形応答とを算出するステップS112。次に、評価値算出部543が、複数の電界変形応答及び複数の磁界変形応答のそれぞれについて評価値を算出し、電界変形応答及び磁界変形応答を一つずつ選択するステップS113。次に、到来時刻推定部544が、電界変形応答及び磁界変形応答に基づいて、到来時刻を推定する。そして、到来時刻推定部544が、推定結果の到来時刻とその識別情報とを対応付けて到来時刻記憶部504に記憶させるステップS114。次に、波形処理部545が、評価値算出部543によって選択された電界変形応答及び磁界変形応答に基づいて波形処理を行うステップS115。そして、終了判定部546が終了判定を行うステップS116。終了判定部546が終了と判定した場合ステップS116-YES、到来時刻判定部540の処理が終了する。一方、処理判定部546が終了と判定しない場合ステップS116-NO、到来時刻判定部540は波形処理部545によって算出された結果直前のステップS115の結果を新たな反射時間応答として用いることによって、ステップS111以降の処理を繰り返し実行する。以上の処理によって、図7Aに示されるように、識別情報と到来時刻とが対応付けられた複数のレコードを有するテーブルが到来時刻記憶部504に記録される。ユーザは、このテーブルをディスプレイやプリンタに出力して参照することによって、受信位置20における電磁波の到来時刻についてさらに解析を進めることが可能となる。例えば、最も大きいピークが到来した時刻が電磁波の到来時刻であると推測するならば、識別情報1に対応する到来時刻を電磁波の真の到来時刻であると推測することが可能である。このように構成された到来時刻推定装置500によれば、時間的な広がりを有した電磁波の到来時刻を精度良く推定することが可能となる。以下、この効果について詳細に説明する。電界のみに基づいて得られるシミュレーション結果と磁界のみに基づいて得られるシミュレーション結果とは異なることが多い。そのため、このずれが上述したような精度の悪化に寄与している場合があった。これに対し、到来時刻推定装置500は、電界に対して得られたシミュレーション結果と磁界に対して得られたシミュレーション結果との双方に基づいて到来時刻を推定する。すなわち、到来時刻推定装置500は、理想電界テーブル及び反射電界テーブルの結果に基づいて得られた候補標準点61の時刻と、理想磁界テーブル及び反射磁界テーブルの結果に基づいて得られた候補標準点61の時刻と、に基づいて到来時刻を推定する。そのため、上記のような電界と磁界とのずれを鑑みた上で推定結果が得られるため、到来時刻の推定精度を向上させることが可能となる。また、単にピーク時の電界強度や磁界強度に基づいて処理を行うのではなく、ピークに基づいて複数の候補標準点61を設定し、その中から評価値に基づいて最適電界時刻及び最適磁界時刻を選択し、これらの時刻に基づいて到来時刻を推定することによって、さらに推定精度を向上させることが可能となる。すなわち、離散時間間隔の間に実際のピークが存在する場合には、単にピーク検出部541によって検出されたピークのみによって処理を行うと誤差を生じさせるおそれがある。このような問題に対し、上述したような処理を行うことで、誤差を小さくすることが可能となる。<変形例>変形応答算出部542は、以下のように変形応答のテーブルを算出しても良い。上述したように、理想環境においてシミュレーションを行う際の離散時間間隔と、反射環境においてシミュレーションを行う際の離散時間間隔とは異なるように設定されても良い。理想環境における離散時間間隔が、反射環境における離散時間間隔よりも短い場合には、理想環境における離散時間間隔に合わせて、反射環境における離散時間間隔よりも細かい間隔で候補標準点61を設定しても良い。図10は、このように変形して構成された変形応答算出部542の処理の概略を示す図である。図10において、電界強度の軸縦軸に沿って伸びる複数の破線は、理想環境及び反射環境のシミュレーションに共通した離散時間を表す。電界強度の軸縦軸に沿って伸びる複数の一点鎖線は、理想環境のシミュレーションのみにおける離散時間を表す。したがって、理想電界テーブル及び理想磁界テーブルには、破線の離散時間及び一点鎖線の離散時間の両方の時刻に関するレコードが記録されている。一方、反射電界テーブル及び反射磁界テーブルには、破線の離散時間のみの時刻に関するレコードが記録されている。図10Aは、図5Aと同じ候補標準点61i=0を表す。図10Bは、図5Bと同じ候補標準点61i=1を表す。図10Cは、変形して構成された変形応答算出部542の処理によって新たに設定される候補標準点i=3を表す。図10Cに示されるように、変形応答算出部542は、理想電界テーブル及び理想磁界テーブルのみに記録されている時刻について候補標準点を設定しても良い。ただし、この場合この時刻に関する電界強度及び磁界強度は反射電界テーブル及び反射磁界テーブルには記録されていない。そのため、変形応答算出部542は、例えばその前後i=3の場合はi=0及びi=1の候補標準点における電界強度及び磁界強度に基づいて、新たに設定した候補標準点の電界強度及び磁界強度を取得しても良い。例えば、前後のいずれか一方の候補標準点61の値を用いても良いし、前後の候補標準点61の値に基づいて補間処理によって算出しても良い。このような処理を行うことによって、上述したような離散時間間隔の間に実際のピークが存在することによる誤差の問題を、より精度良く解決することが可能となる。以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。