(57)【要約】【課題】簡易な構成で、構造物周辺における風に関する情報を精度よく計測すること。【解決手段】計測装置は、音波がスピーカーから出力されてから複数のマイクのそれぞれに到達するまでの反射時間を計測する音波反射時間計測部と、風速、風向、及び対象物との距離と、記反射時間との関係を示す音波反射時間モデルと、反射時間の計測値とに基づいて、風速、及び対象物との距離を推定する推定部と、を備える。【選択図】図１
2019095219.tif 000002
【特許請求の範囲】
【請求項１】
音波がスピーカーから出力されてから複数のマイクのそれぞれに到達するまでの反射時間を計測する音波反射時間計測部と、
風速、風向、及び対象物との距離と、前記反射時間との関係を示す音波反射時間モデルと、前記反射時間の計測値とに基づいて、前記風速、及び前記対象物との距離を推定する推定部と、
を備える計測装置。
【請求項２】
前記複数のマイクは、少なくとも４以上のマイクであり、
前記推定部は、
前記音波反射時間モデルと、前記音波反射時間計測部により計測された前記反射時間の計測値とに基づいて、前記風速、前記風向、前記対象物との距離、及び前記対象物に対する角度を推定する、
請求項１に記載の計測装置。
【請求項３】
前記音波反射時間計測部は、
予め測定した前記スピーカーから出力される音波を直接的に計測した直接音波形と、前記複数のマイクのそれぞれで受信した受信波形との相互相関関数を求めるとともに、該相互相関関数の複数のピーク値と前記音波反射時間モデルに基づいて、前記反射時間の計測値を決定する、
請求項１または請求項２に記載の計測装置。
【請求項４】
前記音波反射時間計測部は、
前記スピーカーから位相をずらした音波を出力させる、
を備える請求項１から請求項３のいずれか一項に記載の計測装置。
【請求項５】
前記スピーカーと、
前記複数のマイクと、
請求項１から請求項４のいずれか一項に記載の計測装置と、
を備える移動体。
【請求項６】
計測装置の計測方法であって、
音波がスピーカーから出力されてから複数のマイクのそれぞれに到達するまでの反射時間を計測する音波反射時間計測ステップと、
風速、風向、及び対象物との距離と、前記反射時間との関係を示す音波反射時間モデルと、前記反射時間の計測値とに基づいて、前記風速、及び前記対象物との距離を推定する推定ステップと、
を有する計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、計測装置、移動体、及び計測方法に関する。
【背景技術】
【０００２】
小型の無人飛行体（例えば、ドローン）の運用機会が増加してきている。例えば、運送、災害情報収集、構造物や設備の点検等での運用が期待されている。しかしながら、ドローンは、軽量であるため、風によって機体姿勢が乱れやすく、また、それによって操縦ミス等のヒューマンエラーが重複しやすいといった懸念がある。そのため、強風や突風等の悪環境化でのロバスト性が求められている。
【０００３】
主に大型の航空機を対象としたものでは、風速遠隔計測センサから得られる事前情報を制御に利用することが検討されている。例えば、風速遠隔計測の手法としては、レーザー光を発射し空気中のエアロゾルで反射する散乱光を解析することで風速を計測する技術（例えば、非引用文献１、２）や、音波を利用して空気中の密度が異なる不連続面で反射する音波の散乱を解析して風向や風速を測定する技術（例えば、非引用文献３、４）が開示されている。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】ドップラーライダーＧａｌｉｏｎの製品紹介（日本気象株式会社）  https://n-kishou.com/offshorewind/galionlidar.html
【非特許文献２】航空機搭載型ドップラーライダーの高高度飛行実証、日本航空宇宙学会論文集  Ｖｏｌ．６２（２０１４）Ｎｏ．６  ｐ．１９８−２０３
【非特許文献３】風速鉛直分布測定装置（ドップラーソーダー）  ＦＡＳシリーズ（英弘精機株式会社）  http://eko.co.jp/windpower/win_products/0322.html
【非特許文献４】ミニドップラーソーダを用いた相関法による風況予測、第２０回風工学シンポジウム論文集、２００８
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述した非特許文献１〜４に記載された技術は、主に大型の航空機を対象としたものであり、装置のサイズ、重量、消費電力が大きいためドローンへの搭載が困難である。また、ドローンは、構造物付近（近傍）を飛行することがあるが、構造物の周辺では風が急激に変化しやすいため、風に関する事前情報の推定が重要となる。しかしながら、従来技術では、レーザー光や音波の散乱よりも構造物等からの反射光や反射波のレベルの方が大きくなり、正確な測定ができない可能性がある。
【０００６】
本発明は、上記した事情に鑑みてなされたもので、簡易な構成で、風に関する情報を計測できる計測装置、移動体、及び計測方法を提供することを目的の一つとする。
【課題を解決するための手段】
【０００７】
本発明は上記の課題を解決するためになされたものであり、本発明の一態様は、音波がスピーカーから出力されてから複数のマイクのそれぞれに到達するまでの反射時間を計測する音波反射時間計測部と、風速、風向、及び対象物との距離と、前記反射時間との関係を示す音波反射時間モデルと、前記反射時間の計測値とに基づいて、前記風速、及び前記対象物との距離を推定する推定部と、を備える計測装置である。
【０００８】
また、本発明の他の態様は、前記スピーカーと、前記複数のマイクと、請求項１から請求項４のいずれか一項に記載の計測装置と、を備える移動体である。
【０００９】
また、本発明の他の態様は、計測装置の計測方法であって、音波がスピーカーから出力されてから複数のマイクのそれぞれに到達するまでの反射時間を計測する音波反射時間計測ステップと、風速、風向、及び対象物との距離と、前記反射時間との関係を示す音波反射時間モデルと、前記反射時間の計測値とに基づいて、前記風速、及び前記対象物との距離を推定する推定ステップと、を有する計測方法である。
【発明の効果】
【００１０】
本発明によれば、簡易な構成で、構造物周辺における風に関する情報を精度よく計測できる。
【図面の簡単な説明】
【００１１】
【図１】第１の実施形態に係る風速遠隔計測手法の概念を示す図。
【図２】第１の実施形態に係る音波反射モデルを示す図。
【図３】スピーカーへの入力（印加電圧波形）を示す図。
【図４】マイクの出力（受信波形）を示す図。
【図５】スピーカーへの入力（途中で位相を変更した印加電圧波形）を示す図。
【図６】マイクの出力（途中で位相を変更した音波の受信波形）を示す図。
【図７】直接音波形と受信波形との相互相関関数の値の算出結果の一例を示す図。
【図８】第１の実施形態に係るドローンの構成の一例を示すブロック図。
【図９】第１の実施形態に係る制御部の構成の一例を示すブロック図。
【図１０】第１の実施形態に係る風速遠隔計測処理の一例を示すフローチャート。
【図１１】実験の概要を説明する図。
【図１２】実験に用いたスピーカーとマイクの設置状態の写真。
【図１３】実験に用いた風洞の写真。
【図１４】実験条件を示す表。
【図１５】実験システムの構成を示すブロック図。
【図１６】実験結果を示す図。
【図１７】第２の実施形態に係る音波反射モデルの一例を示す図。
【図１８】第３の実施形態に係る音波反射モデルの一例を示す図。
【発明を実施するための形態】
【００１２】
以下、図面を参照しながら本発明の実施形態について説明する。なお、各図において同一部分には同一符号を付している。
［第１の実施形態］
本発明の第１の実施形態について説明する。
まず、本発明の実施形態に係る風速遠隔計測手法の概要について説明する。
図１は、本発明の実施形態に係る風速遠隔計測手法の概要を示す図である。本実施形態に係る風速遠隔計測装置は、構造物等に対する音波の反射時間が風によって変化することを利用して計測する。具体的には、風速遠隔計測装置は、スピーカー１５から音波（例えば、超音波）を出力させ、構造物等の壁面で反射した音波をマイク１６で受信し、受信波形（受信音波波形）を計測する。また、風速遠隔計測装置は、この受信波形に基づいて音波反射時間判別手法を用いて音波の反射時間の判別を行い、音波の反射時間の計測値（音波反射時間計測値）を決定する。次に、風速遠隔計測装置は、風速、風向、構造物等の壁面との距離及び角度等と、音波の反射時間との関係を示す音波の反射時間モデル式（音波反射時間モデル式）と上記音波反射時間計測値とに基づいて、推定アルゴリズムを用いて、風速、風向、構造物等の壁面との距離または角度等を推定する（いわゆる逆問題）。なお、反射時間モデル式は、数式でなくてもよく、例えば値の集合であってもよい。
以下、上述の音波反射時間モデル式、音波反射時間計測値の判別手法、及び音波反射時間モデル式と音波反射時間計測値とに基づく推定アルゴリズムについて詳しく説明する。
【００１３】
（音波反射時間モデル式）
まず、本実施形態に係る音波反射時間モデル式について説明する。図２は、本実施形態に係る音波反射モデルを示す図である。この図は、２次元空間における音波反射モデルの例であり、ｘｙ座標で定義されている。また、ここでは、マイクが２個の場合の音波反射モデルを示している。図示する音波反射モデルにおいて、風速をＶ、風向（壁面に対する風向の角度）をΦとする。
【００１４】
スピーカー１５の出力の中心をｘ軸とｙ軸との交点（原点）とし、２個のマイク１６（１６ａ、１６ｂ）がｘ軸上にスピーカー１５を挟んで相対する位置にあるとする。マイク１６ａがマイク１、マイク１６ｂがマイク２であり、マイクの後ろに付している「１」、「２」は、マイクを識別するためのマイク番号である。スピーカー１５とマイク１６との距離をｄｉとする（ｉは、マイク番号）と、マイク１６ａ（マイク１）のｘｙ座標が（ｄ１、０）、マイク１６ｂ（マイク２）のｘｙ座標が（−ｄ２、０）となる。また、スピーカー１５の出力面に対して垂直方向（ｙ軸方向）にスピーカー１５から壁面までの距離をＬとすると、壁面とｙ軸との交点のｘｙ座標が（０、Ｌ）となる。また、壁面に対するセンサ角度（壁面に対するスピーカー１５及びマイク１６の角度、即ち、壁面とｘ軸との角度）をθとする。
【００１５】
また、スピーカー１５から出力された音波がマイク番号ｉのマイク１６に到達するまでの間に壁面に反射する位置のｘｙ座標を（ｘｉ、ｙｉ）とし、スピーカー１５から壁面のｘｙ座標（ｘｉ、ｙｉ）の位置までの距離をｌｉｎｉ、壁面のｘｙ座標（ｘｉ、ｙｉ）の位置からマイク１６までの距離をｌｒｅｉとする。例えば、スピーカー１５から出力された音波が壁面のｘｙ座標（ｘ１、ｙ１）の位置で反射してマイク１６ａ（マイク１）に到達するとした場合、スピーカー１５から壁面のｘｙ座標（ｘ１、ｙ１）の位置までの距離をｌｉｎ１、壁面のｘｙ座標（ｘ１、ｙ１）の位置からマイク１６ａ（マイク１）までの距離をｌｒｅ１とする。また、スピーカー１５から壁面のｘｙ座標（ｘ１、ｙ１）への線分とｘ軸との角度をθｉｎ１とする。一方、スピーカー１５から出力された音波が壁面のｘｙ座標（ｘ２、ｙ２）の位置で反射してマイク１６ｂ（マイク２）に到達するとした場合、スピーカー１５から壁面のｘｙ座標（ｘ２、ｙ２）の位置までの距離をｌｉｎ２、壁面のｘｙ座標（ｘ２、ｙ２）の位置からマイク１６ｂ（マイク２）までの距離をｌｒｅ２とする。また、スピーカー１５から壁面のｘｙ座標（ｘ２、ｙ２）への線分とｘ軸との角度をθｉｎ２、とする。
【００１６】
上述した図２に示す音波反射モデルにおいて、音波反射時間計算値Ｔｉを算出する音波反射時間モデル式は、以下に示す数式１となる。なお、Ｃは音速である。
【００１７】
【数１】
2019095219.tif 000003
【００１８】
上記数式１に示す音波反射時間モデル式を用いて、音波反射時間計算値を求めることができる。
【００１９】
（推定アルゴリズム）
次に、上記音波反射時間モデル式と音波反射時間計測値とに基づいて、風速、風向、壁面との距離または角度等を推定する推定アルゴリズムについて説明する。
風速Ｖ、風向Φ、壁面との距離Ｌ、及び壁面に対するセンサ角度θを最適化パラメータｘとすると、音波反射時間モデル式は、以下の数式２のように表すことができる。
【００２０】
【数２】
2019095219.tif 000004
【００２１】
また、実際にスピーカー１５から出力された音波が壁面で反射した反射波をマイク１６で受信した受信結果に基づく音波反射時間計測値を以下のように表す。
【００２２】
【数３】
2019095219.tif 000005
【００２３】
そして、上述の音波反射時間モデル式及び音波反射時間計測値を用いて、風速、風向、壁面との距離や角度等を推定する推定アルゴリズムとして、評価関数Ｊを以下の数式３で表すことができる。
【００２４】
【数４】
2019095219.tif 000006
【００２５】
上記の数式３は、評価関数Ｊを最小化する最適化パラメータｘを求めることを示している。評価関数Ｊを最小化する最適化パラメータｘを求めることにより、風速Ｖ、風向Φ、壁面との距離Ｌ、及び壁面に対するセンサ角度θを推定することができる。
【００２６】
例えば、風向Φ、及び壁面に対するセンサ角度θを既知の値とし、図２に示す音波反射モデルで２個のマイク１６を用いて風速Ｖ及び壁面との距離Ｌを推定する場合には、以下の数式４に示すように最適化パラメータｘは、風速Ｖと、壁面との距離Ｌとの変数となる。
【００２７】
【数５】
2019095219.tif 000007
【００２８】
この場合、最適化パラメータｘを用いて、音波反射時間モデル式は、以下の数式５のように表すことができる。
【００２９】
【数６】
2019095219.tif 000008
【００３０】
そして、推定アルゴリズムとして、評価関数Ｊを以下の数式６で表すことができる。
【００３１】
【数７】
2019095219.tif 000009
【００３２】
この数式５において、評価関数Ｊを最小化する最適化パラメータｘを求めることにより、風速Ｖ及び壁面との距離Ｌを推定することができる。
【００３３】
（音波反射時間判別手法）
次に、スピーカー１５から出力した音波の反射波をマイク１６で受信した受信波形に基づいて音波反射時間計測値を決定する際の音波反射時間判別手法について説明する。
スピーカー１５は、例えば、超音波を出力する超音波スピーカーである。理想的には音波の発信及び停止を急激に行うことができることが望ましいが、機械的な共振現象を利用しているため、それが難しい。
【００３４】
図３及び図４に、スピーカー１５とマイク１６とを１０ｃｍの距離で対向させ、スピーカー１５への入力（印加電圧波形）と、該入力に応じてスピーカー１５から出力される音波（超音波）を直接的に受信したときのマイク１６の出力（受信波形）とを示す。図３は、スピーカー１５への入力（印加電圧波形）を示す図であり、縦軸がスピーカー１５への入力電圧、横軸が時間である。一方、図４は、マイク１６の出力（受信波形）を示す図であり、縦軸がマイク１６の出力電圧、横軸が音波の到達時刻からの時間である。図３及び図４を参照すると、スピーカー１５への入力に対してマイク１６の振幅の増加及び減少が遅いことが分かる。つまり、スピーカー１５は、音波の発信及び停止を急激に行うことができないことが分かる。また、図４を参照すると、このマイク１６の受信波形から音波の一波目を検出することは難しいことが分かる。よって、受信波形に基づいて音波反射時間計測値を決定する際においても、マイク１６の受信波形から音波の反射波の一波目を検出することが難しく、反射時間の判別が困難である。また、スピーカー１５への印加電圧波形とマイク１６の受信波形との相互相関関数を求めることによる反射時間の判別が困難である。
【００３５】
また、スピーカー１５の共振周波数と印加電圧波形の周波数とを一致させないと大きな音圧が得られないという問題もある。このため、周波数スイープした場合に、スピーカー１５への印加電圧波形とマイク１６の受信波形との相互相関関数を求めることが困難であり、反射時間の判別も困難である。
【００３６】
そこで、まずスピーカー１５とマイク１６とを対向させて、スピーカー１５から出力される音波を直接的に計測した直接音波形（例えば、図４に示す受信波形）を事前に予め測定しておく。そして、この直接音波形と、実際に壁面に当たって反射した音波の受信波形との相互相関関数を求めることにより、音波反射時間を判別する。これにより、スピーカー１５への入力（印加電圧波形）との相互相関関数を求めるよりは、正確に音波反射時間の計測値を決定することができる。
【００３７】
また、さらに特徴のある受信波形が得られるような印加電圧波形をスピーカー１５への入力（印加）することにより、音波反射時間の判別精度を向上させてもよい。例えば、図５に示すように、測定の途中で位相を変更した印加電圧波形をスピーカー１５への入力する。図５は、スピーカー１５への入力（印加電圧波形）を示す図であり、測定の途中で位相を変更した印加電圧波形の一例を示す図である。この図において、縦軸がスピーカー１５への入力電圧、横軸が時間である。図示する例では、途中で位相を１８０度変化させた４０ｋＨｚ矩形波信号がスピーカー１５への入力されている。図６は、スピーカー１５に対向するマイク１６の出力（途中で位相を変更した音波の受信波形）を示す図であり、上記図５に示す印加電圧波形のときにスピーカー１５から出力される音波（超音波）を直接的に受信したマイク１６の出力（受信波形）の一例を示す図である。この図６に示すように、音圧を維持しつつ、途中で位相が変更されたことにより、特徴のある受信波形が得られるため、相互相関関数による音波反射時間の判別の精度を向上させることができる。
【００３８】
なお、上述したように、スピーカー１５から出力される音波を直接的に計測した直接音波形と、実際に壁面に反射した音波の受信波形との相互相関関数を求めることにより、音波反射時間の判別精度を向上させることができるが、それでもなお、様々な計測誤差等により、直接音波形と比較して受信波形には反射の過程で変化が生じてしまうことがある。それにより、例えば、直接音波形と受信波形との相互相関関数による音波反射時間の判別結果が整数波長分ずれてしまうことがある。
【００３９】
図７は、直接音波形と受信波形との相互相関関数の算出結果の一例を示す図である。この図において、縦軸が相互相関関数の値であり、横軸が時間シフト量である。例えば、相互相関関数の値が最大となる時間シフトτ３を反射時間とするところ、上述したように、様々な計測誤差等により受信波形に変化が生じている場合、この時間シフトτ３が正しい反射時間とは限らない。
【００４０】
そこで、本実施形態では、相互相関関数の最大値を含む複数のピーク値に基づいて、音波反射時間計測値を決定する。例えば、図７に示すように、相互相関関数の最大値に加えて前後それぞれの２つのピーク値に対応する５つの時間シフトτ１〜τ５を反射時間の候補とする。そして、マイク毎に時間シフトτ１〜τ５を求める。例えば、マイク１６ａ（マイク１）とマイク１６ｂ（マイク２）の２個のマイクを使用した場合には、２個のマイクのそれぞれで５つの時間シフトτ１〜τ５を求める。そして、求めた全ての組み合わせに対して、前述した推定アルゴリズム（例えば、数式５参照）を適用し、解が得られた時間シフト量を音波反射時間計測値とする。これにより、様々な計測誤差等により、直接音波形と比較して受信波形に反射の過程で変化が生じてしまうことがあっても、相互相関関数による音波反射時間の判別精度をより向上させることができる。
【００４１】
次に、本実施形態に係る風速遠隔計測手法を適用した小型の無人飛行体の構成例を説明する。小型の無人飛行体とは、例えばドローンであり、遠隔操作や自動制御によって無人で飛行できる小型の航空機である。
図８は、本実施形態に係るドローンの構成の一例を示すブロック図である。図示するドローン１０は、通信部１１と、駆動部１２と、記憶部１３と、制御部１４と、スピーカー１５と、２個のマイク１６（１６ａ、１６ｂ）とを備えている。
【００４２】
通信部１１は、ドローン１０を遠隔から操縦するためのコントローラ（不図示）と無線通信を行う。例えば、通信部１１は、ユーザ操作に基づいてコントローラ（不図示）から操縦内容を示す情報を受信する。また、通信部１１は、必要に応じて、ドローン１０の状態を示す情報をコントローラ（不図示）へ送信する。
【００４３】
なお、通信部１１は、インターネット等の通信ネットワークを介してサーバ装置や他の機器と通信してもよい。
【００４４】
駆動部１２は、複数のモータと、各モータのモータ軸の回転にそれぞれ連動する複数のプロペラとを含んで構成されている。駆動部１２は、制御部１４の制御に応じて各モータの回転を駆動し、プロペラを回転させることで揚力や推力を発生させ、ドローン１０の飛行を可能にする。
【００４５】
なお、ドローン１０は、本図では図示していないが、例えば、機体の動きや姿勢を検出するための加速度センサ、ジャイロセンサ、または方位センサや、各部に給電するバッテリー等を備えている。また、ドローン１０は、上空から撮影が可能なカメラ等を備えていてもよい。
【００４６】
記憶部１３は、ドローン１０の各種の情報を記憶する。例えば、記憶部１３は、制御部１４による飛行制御に関する制御プログラムを記憶する。また、記憶部１３は、前述した風速遠隔計測手法に関する音波反射時間モデルに関するパラメータ、音波反射時間モデル式、推定アルゴリズム等を記憶する。
【００４７】
制御部１４は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を備え、記憶部１３に記憶されている制御プログラムを実行することにより、ドローン１０を制御するための各種の処理を実行する。例えば、制御部１４は、ユーザが操縦するコントローラ（不図示）から通信部１１を介して受信する操縦内容を示す情報、及びドローン１０の機体の動きや姿勢の検出結果等に基づいて、駆動部１２の駆動を制御することによりドローン１０の飛行を制御する。また、制御部１４は、スピーカー１５と、マイク１６（１６ａ、１６ｂ）とを用いて、前述した風速遠隔計測手法により、ドローン１０の周辺の風に関する情報を推定し、ドローン１０の飛行制御に利用する。
【００４８】
スピーカー１５は、例えば、超音波を出力可能な超音波スピーカーである。マイク１６（１６ａ、１６ｂ）は、例えば、超音波を受信可能な超音波マイクである。スピーカー１５及び２個のマイク１６（１６ａ、１６ｂ）は、ドローン１０の機体の側面に設けられている。具体的には、図２に示すように、スピーカー１５と２個のマイク１６（１６ａ、１６ｂ）とが直線上に、且つ２個のマイク１６（１６ａ、１６ｂ）がスピーカー１５を挟んで相対する位置になるように、機体の側面に設置されている。例えば、スピーカー１５と２個のマイク１６（１６ａ、１６ｂ）とは、飛行時の水平方向に対応する方向へ直線上に並べて機体の側面に設置されている。また、スピーカー１５は、機体面に対面する方向に音波（超音波）を出力するように設置されている。２個のマイク１６（１６ａ、１６ｂ）のそれぞれは、スピーカー１５から出力され音波（超音波）の反射波を受信可能なように、機体面に対面する方向に向けて設置されている。
【００４９】
なお、上述の説明では、スピーカー１５及び複数のマイク１６は、直線上に設置されている例を説明した。このように、スピーカー１５及び複数のマイク１６が一直線に並ぶことで、簡潔なモデル式とすることができ、それにより音波反射時間の算出を容易に行うことができる。なお、スピーカー１５及び複数のマイク１６は、直線上に設置されなくてもよく、直線上からずれた位置に設置されても音波反射時間の算出は可能である。また、スピーカー１５及び２個のマイク１６（１６ａ、１６ｂ）は、ドローン１０の機体の側面の一つに設置されてもよいし、複数の側面のそれぞれに設置されてもよい。例えば、機体の側面の一つに設置されている場合、該側面が進行方向を向くように操縦され、進行方向の風の情報を遠隔計測してもよい。また、機体の側面の全て（例えば、側面が４面ある場合には４面の全て）に設置することにより、機体の全周方向（例えば、４方向）の風の情報をまとめて遠隔計測してもよい。また、スピーカー１５及び２個のマイク１６（１６ａ、１６ｂ）は、ドローン１０の周辺に音波を出力し且つその音波の反射波を受信可能な位置あれば、機体の側面に限らず、機体の任意の場所に設置することができる。
【００５０】
次に、図９を参照して、制御部１４の構成について詳しく説明する。ここでは、制御部１４の構成のうち、風速遠隔計測手法によりドローン１０の周辺の風に関する情報を推定する機能構成について説明する。図９は、本実施形態に係る制御部１４の構成の一例を示すブロック図である。図示する制御部１４は、前述した風速遠隔計測手法により、ドローン１０の周辺の風に関する情報を推定する計測装置１４１と、その推定結果等に基づいて駆動部１２の駆動を制御する駆動制御部１４２とを備えている。
【００５１】
計測装置１４１は、音波反射時間計測部１４１１と、推定部１４１２とを備えている。音波反射時間計測部１４１１は、スピーカー１５から出力した音波の壁面からの反射波を２個のマイク１６（１６ａ、１６ｂ）で受信した受信結果に基づいて、音波がスピーカー１５から出力されてから２個のマイク１６（１６ａ、１６ｂ）のそれぞれに到達するまでの反射時間を計測する。例えば、音波反射時間計測部１４１１は、図５に示すように、スピーカー１５から所定の周波数の音波（例えば、超音波）を出力させる。ここで、音波反射時間計測部１４１１は、位相を変更させた（ずらした）音波を出力させる。例えば、音波反射時間計測部１４１１は、途中で位相を１８０度変更する。なお、変更する位相の量は、１８０度に限られるものではなく、例えば９０度等としてもよい。
【００５２】
また、音波反射時間計測部１４１１は、予め測定したスピーカー１５から出力される音波を直接的に計測した直接音波形と、２個のマイク１６（１６ａ、１６ｂ）のそれぞれで受信した受信波形との相互相関関数を求めるとともに、該相互相関関数の最大値を含む複数のピーク値に基づいて、音波反射時間計測値を決定してもよい。例えば、音波反射時間計測部１４１１は、図７に示すように、直接音波形と受信波形との相互相関関数最大値に加えて前後それぞれの２つのピーク値に対応する５つの時間シフトτ１〜τ５を反射時間の候補（音波反射時間計測値の候補）とし、マイク毎に時間シフトτ１〜τ５を求める。
【００５３】
推定部１４１２は、スピーカー１５または２個のマイク１６（１６ａ、１６ｂ）と壁面との距離、風速、及び風向と反射時間との関係を示す音波反射時間モデル式（図２及び数式１参照）と、音波反射時間計測部１４１１により計測された音波反射時間計測値とに基づいて、推定アルゴリズム（例えば、数式３、数式５参照）を用いて、風に関する情報を推定する。例えば、推定部１４１２は、風向Φ、及び壁面に対するセンサ角度θを既知の値として、５つの音波反射時間計測値の候補の全ての組み合わせに対して推定アルゴリズム（数式３、数式５参照）を適用し、評価関数Ｊを最小化する最適化パラメータｘを求めることにより、ドローン１０の周辺の風速及び壁面との距離を推定する。
【００５４】
駆動制御部１４３は、計測装置１４１による推定結果等に基づいて、駆動部１２の駆動を制御する。例えば、駆動制御部１４３は、ユーザが操縦するコントローラ（不図示）から通信部１１を介して受信する操縦内容を示す情報、及びドローン１０の機体の動きや姿勢の検出結果等に基づいて駆動部１２の駆動を制御するとともに、計測装置１４１による推定結果に基づいてドローン１０の周辺の風に関する情報を加味して駆動部１２の駆動を制御する。例えば、駆動制御部１４３は、計測装置１４１が遠隔計測した風の分布を利用して，風の影響を受けずに目的地まで移動する最適経路の生成や，風の影響を軽減する最適制御の生成を行い、それらに基づき駆動部１２の駆動を制御することもできる。これにより、強風や突風等の悪環境化におけるドローン１０の飛行のロバスト性を高め、安定的な飛行制御を可能とする。
【００５５】
次に、本実施形態に係る計測装置１４１が風速遠隔計測手法により風に関する情報を推定する風速遠隔計測処理について説明する。
図１０は、本実施形態に係る風速遠隔計測処理の一例を示すフローチャートである。
【００５６】
（ステップＳ１０）計測装置１４１は、スピーカー１５から位相を変更させた音波（例えば、超音波）を出力させるとともに、出力させる音波の位相を計測の途中で変更する。例えば、計測装置１４１は、途中で位相を１８０度変更する（図５参照）。そして、ステップＳ１２の処理に進む。
【００５７】
（ステップＳ１２）計測装置１４１は、スピーカー１５から出力した音波例えば、超音波）の壁面からの反射波を２個のマイク１６（１６ａ、１６ｂ）で受信した受信結果を取得する。そして、ステップＳ１４の処理に進む。
【００５８】
（ステップＳ１４）計測装置１４１は、予め測定したスピーカー１５から出力される音波を直接的に計測した直接音波形（図６参照）と、ステップＳ１２で取得した受信結果（２個のマイク１６（１６ａ、１６ｂ）のそれぞれで受信した受信波形）との相互相関関数を求めるとともに、該相互相関関数の最大値を含む複数のピーク値に基づいて、音波反射時間計測値を決定する。例えば、計測装置１４１は、直接音波形と受信波形との相互相関関数最大値に加えて前後それぞれの２つのピーク値に対応する５つの時間シフトτ１〜τ５を音波反射時間計測値の候補とし、マイク毎に時間シフトτ１〜τ５を求める（図７参照）。そして、ステップＳ１６の処理に進む。
【００５９】
（ステップＳ１６）計測装置１４１は、スピーカー１５または２個のマイク１６（１６ａ、１６ｂ）と壁面との距離、風速、及び風向と反射時間との関係を示す音波反射時間モデル式（図２及び数式１参照）と、ステップＳ１４において計測した音波反射時間計測値とに基づいて、推定アルゴリズム（数式３、数式５等参照）を用いて、風に関する情報を推定する。例えば、計測装置１４１は、風向Φ、及び壁面に対するセンサ角度θを既知の値として、５つの音波反射時間計測値の候補の全ての組み合わせに対して推定アルゴリズム（数式３、数式５参照）を適用し、評価関数Ｊを最小化する最適化パラメータｘを求めることにより、ドローン１０の周辺の風速Ｖ及び壁面との距離Ｌを推定する。
【００６０】
（風速遠隔計測手法の実験例）
次に、上記で説明した風速遠隔計測手法の妥当性を検証するための実験結果の一例を説明する。
図１１は、実験の概要を説明する図である。実験では、風洞内に１個のスピーカーと２個のマイク（マイク１、マイク２）が設置されたものを使用した。風洞の内壁が壁面に相当し、基本的には図２に示す音波反射モデルの条件を再現した。図１２は、実験に用いたスピーカーとマイクの設置状態の写真である。なお、この写真では、マイク１とマイク２のそれぞれに２個のマイクが設置されているが、本実験では、それぞれ１個のマイクのみ（合計２個のマイク）を実験に使用している。また、図１３は、実験に用いた風洞の写真である。風洞は、金沢大学航空宇宙システム研究室所蔵の、流速：２〜２５ｍ/s、試験区間：２００×２００×５００ｍｍのものを使用した。
【００６１】
図１４は、実験条件を示す表である。マイク１及びマイク２は、スピーカーを挟んで直線上に相対する位置にあり、スピーカーとの距離ｄを２０ｃｍとした。また、スピーカーと風洞の内壁との距離Ｌ０（内壁から内壁面に対して垂直方向の距離）を、１２．２ｃｍとした。また、風向は０ｄｅｇとし、内壁に対するセンサ角度θを、−５、０、５ｄｅｇの３種類に変化させて測定した。また、風速Ｖを、０、２、４、６、８、１０ｍ／ｓの６種類に変化させて測定した。なお、サンプリング周期を２μｓ、音波の発信間隔を０．０１ｓとした。
【００６２】
また、図１５は、実験システムの構成を示すブロック図である。ファンクションジェネレータから出力されるスピーカー駆動信号がスピーカーアンプとデータロガーに入力される。スピーカーアンプは、スピーカー駆動信号を増幅してスピーカーへ出力する。マイクアンプは、マイクの出力信号を増幅した増幅信号をデータロガーへ出力する。データロガーは、マイクアンプから出力される増幅信号を、ファンクションジェネレータから出力されるスピーカー駆動信号に基づいてサンプリングして記録し、パーソナルコンピュータ（パソコン）へ出力する。なお、マイクアンプの数は、マイクの数に対応して備えられている。ここでは、２個のマイクに対応する２個のマイクアンプが備えられている。各マイクアンプは、それぞれに接続されているマイクの出力信号を増幅してデータロガーへ出力する。なお、この実験では、ファンクションジェネレータは、エヌエフ回路設計ブロック製ＤＦ１９０６を使用した。また、データロガーは、ＫＥＹＥＮＣＥ製ＮＲ-６００を使用した。
【００６３】
図１６は、実験結果を示す図である。ここでは、実験結果として、風速Ｖ及び壁面との距離Ｌの推定結果をグラフ化して示している。この図において、上段は、風速Ｖの推定結果を示しており、内壁に対するセンサ角度θを、−５、０、５ｄｅｇの３種類に変化させた場合の３種類のグラフを示している。一方、下段は、壁面との距離Ｌの推定結果を示しており、同様に内壁に対するセンサ角度θを、−５、０、５ｄｅｇの３種類に変化させた場合の３種類のグラフを示している。いずれのグラフも横軸は設定した風速（設定風速）である。これらのグラフにおいて、設定風速０、２、４、６、８、１０ｍ／ｓのそれぞれについての推定値と真値とを示している。これらのグラフを参照すると、推定結果と真値とに大きな差は無く、風速Ｖ及び壁面との距離Ｌを精度よく推定できていることが分かる。
【００６４】
以上説明したように、本実施形態に係る計測装置１４１は、音波反射時間計測部１４１１と、推定部１４１２とを備えている。音波反射時間計測部１４１１は、スピーカー１５から出力した音波の壁面（対象物の一例）からの反射波を複数のマイクで受信した受信結果に基づいて、該音波がスピーカー１５から出力されてから複数のマイクのそれぞれに到達するまでの反射時間（音波反射時間計測値）を計測する。推定部１４１２は、風速、風向、及びスピーカー１５または複数のマイク１６と壁面（対象物の一例）との距離と、上記反射時間との関係を示す音波反射時間モデル式と、音波反射時間計測部１４１１により計測された音波反射時間計測値とに基づいて、風速及び壁面（対象物の一例）との距離を推定する。
【００６５】
これにより、計測装置１４１は、例えば１個のスピーカー１５と２個のマイク１６（１６ａ、１６ｂ）を用いた簡易な構成で、風に関する情報と壁面（対象物の一例）までの距離とをまとめて推定できるため、構造物周辺における風に関する情報（例えば、風速などの事前情報）を精度よく計測することができる。
【００６６】
例えば、音波反射時間計測部１４１１は、スピーカー１５から位相を変更させた音波を出力（例えば、１８０度変更）してもよい。これにより、計測装置１４１は、音圧を維持しつつ、途中で位相が変更されたことにより、特徴のある受信波形が得られるため、相互相関関数による音波反射時間の判別の精度を向上させることができる。
【００６７】
また、音波反射時間計測部１４１１は、予め測定したスピーカー１５から出力される音波を直接的に計測した直接音波形と、複数のマイク１６のそれぞれで受信した受信波形との相互相関関数を求めるとともに、該相互相関関数の複数のピーク値に基づいて、音波反射時間計測値を決定する。例えば、音波反射時間計測部１４１１は、相互相関関数の最大値を含む複数のピーク値に基づいて、音波反射時間計測値を決定する。
【００６８】
これにより、計測装置１４１は、直接音波形と比較して受信波形に反射の過程で変化が生じてしまうことがあっても、相互相関関数による音波反射時間の判別精度をより向上させることができる。
【００６９】
なお、本実施形態では、図７に示すように、相互相関関数の最大値に加えて前後それぞれの２つのピーク値に基づいて音波反射時間計測値を決定する例を説明したが、これに限られるものではない。例えば、相互相関関数の最大値の前後で用いるピーク値の数は任意に決めることができる。また、相互相関関数の最大値以外の複数のピーク値を用いてもよい。
【００７０】
また、計測装置１４１を備えたドローン１０は、ドローン１０の周辺の風に関する情報を加味してドローン１０の飛行を制御することができるため、強風や突風等の悪環境化におけるドローン１０の飛行のロバスト性を高め、安定的な飛行制御を可能とする。例えば、ドローン１０は、計測装置１４１が遠隔計測した風の分布を利用して，風の影響を受けずに目的地まで移動する最適経路の生成や，風の影響を軽減する最適制御の生成を行い、それらに基づき駆動部１２の駆動を制御することもできる。
【００７１】
（第２の実施形態）
次に、本発明の第２の実施形態について説明する。第１の実施形態では、２個のマイクを使用して風速Ｖ及び壁面との距離Ｌを推定する例を説明したが、マイクを増設することにより、風向Φ、及び壁面に対するセンサ角度θについても推定が可能となる。本実施形態では、マイクを４個備える構成について説明する。
【００７２】
図１７は、本実施形態に係るマイクが４個の場合の音波反射モデルの一例を示す図である。この図に示す音波反射モデルは、図２に示す音波反射モデルに対してマイクの個数が異なる。スピーカー１５を挟んでマイク１６ａ（マイク１）及びマイク１６ｃ（マイク３）とマイク１６ｂ（マイク２）及びマイク１６ｄ（マイク４）とが相対する位置にある。マイクが４個の場合も、音波の反射時間の計算値Ｔｉを算出する音波反射時間モデル式は、前述した数式１で表すことができる。
【００７３】
最適化パラメータｘは、以下の数式７に示すように風速Ｖ、風向Φ、壁面との距離Ｌ、及び壁面に対するセンサ角度θの変数となる。
【００７４】
【数８】
2019095219.tif 000010
【００７５】
また、音波反射時間モデル式は、最適化パラメータｘを用いて以下の数式８のように表すことができる。
【００７６】
【数９】
2019095219.tif 000011
【００７７】
また、音波反射時間計測値は、第１の実施形態と同様に以下のように表す。
【００７８】
【数１０】
2019095219.tif 000012
【００７９】
そして、風速Ｖ、風向Φ、壁面との距離Ｌ、及び壁面に対するセンサ角度θを推定する推定アルゴリズムとして、評価関数Ｊを以下の数式９のように表すことができる。
【００８０】
【数１１】
2019095219.tif 000013
【００８１】
この数式９において、評価関数Ｊを最小化する最適化パラメータｘを求めることにより、風速Ｖ、風向Φ、壁面との距離Ｌ、及び壁面に対するセンサ角度θを推定することができる。
【００８２】
このように、本実施形態では、計測装置１４１は、４個のマイク１６を備えている。そして、推定部１４１２は、音波反射時間モデル式と、音波反射時間計測部１４１１により計測された音波反射時間計測値とに基づいて、風速、風向、壁面（対象物の一例）との距離、及び壁面（対象物の一例）に対するスピーカー１５及びマイク１６の角度を推定する。
【００８３】
これにより、計測装置１４１は、例えば１個のスピーカー１５と４個のマイク１６（１６ａ、１６ｂ）を用いた簡易な構成で、風速と壁面（対象物の一例）までの距離とに加えて、風向、及び壁面に対するセンサ角度についても推定が可能となるため、構造物周辺におけるより多くの風に関する情報を計測することができる。
【００８４】
なお、設置するマイクの数は、少なくとも４個以上であればよく、例えば、５個または６個等であってもよい。
【００８５】
（第３の実施形態）
次に、本発明の第３の実施形態について説明する。第１及び第２の実施形態では、２次元空間における風速、風向、壁面との距離、または壁面に対するセンサ角度等を推定する例を説明したが、本実施形態では３次元空間に拡張した例を説明する。
【００８６】
図１８は、本実施形態に係る３次元空間での音波反射モデルの一例を示す図である。この図に示す音波反射モデルは、図２に示す音波反射モデルが２次元空間をｘｙ座標で定義していたのに対して、３次元空間をｘｙｚ座標で定義している。ここでは、ドローン１０等の機体（センシングシステム（測定装置）筐体）に固定した座標軸を、ｘａ軸、ｙａ軸、及びｚａ軸と定義している。スピーカー１５の出力の中心をｘａ軸、ｙａ軸、及びｚａ軸の交点（原点）とし、ドローン１０等の機体（センシングシステム（測定装置）筐体）に、スピーカー１５と、少なくとも３個以上のマイクが設けられているとする。図示する例では、スピーカー１５を中心として、放射方向にマイク１６ａ（マイク１）、マイク１６ｂ（マイク２）、及びマイク１６ｉ（マイクｉ）がそれぞれ配置されており、スピーカー１５から各マイク１６までの距離ベクトルをｄｉ（ｉはマイク番号）とする。
【００８７】
また、図２の例と同様に、スピーカー１５の出力面に対して垂直方向にスピーカー１５から壁面までの距離をＬとし、スピーカー１５から出力された音波が壁面で反射してマイク１６ａ（マイク１）に到達する経路において、スピーカー１５から壁面までの距離をｌｉｎ１、壁面からマイク１６ａ（マイク１）までの距離をｌｒｅ１とする。また、スピーカー１５から出力された音波が壁面で反射してマイク１６ｂ（マイク２）に到達する経路において、スピーカー１５から壁面までの距離をｌｉｎ２、壁面からマイク１６ｂ（マイク２）までの距離をｌｒｅ２とする。また、スピーカー１５から出力された音波が壁面で反射してマイク１６i（マイクｉ）に到達する経路において、スピーカー１５から壁面までの距離をｌｉｎｉ、壁面からマイク１６ｂ（マイク２）までの距離をｌｒｅｉとする。
【００８８】
また、壁面に対する機体又はセンサ姿勢角度（壁面に対するスピーカー１５及びマイク１６の角度）をθ１、θ２とし、風速ベクトルをＶ、音速をＣとすると、音波の反射時間の計算値Ｔｉを算出する音波反射時間モデル式は、以下に示す数式１０で表すことができる。
【００８９】
【数１２】
2019095219.tif 000014
【００９０】
最適化パラメータｘは、以下の数式１１に示すように風速ベクトルＶ、壁面との距離Ｌ、及び壁面に対するセンサ姿勢角度θ１、θ２の変数となる。
【００９１】
【数１３】
2019095219.tif 000015
【００９２】
また、音波反射時間モデル式は、最適化パラメータｘを用いて以下の数式１２のように表すことができる。
【００９３】
【数１４】
2019095219.tif 000016
【００９４】
また、音波反射時間計測値は、第１の実施形態と同様に以下のように表す。
【００９５】
【数１５】
2019095219.tif 000017
【００９６】
そして、推定アルゴリズムとして、評価関数Ｊを以下の数式１３のように表すことができる。
【００９７】
【数１６】
2019095219.tif 000018
【００９８】
このように、３次元空間に対応した音波反射モデルを用いて、上記数式１３において、評価関数Ｊを最小化する最適化パラメータｘを求めることにより、３次元空間における風速ベクトルＶ及び壁面との距離Ｌを推定することができる。また、スピーカーをさらに増設することにより、３次元空間における風速ベクトルＶ、壁面との距離Ｌに加えて、壁面に対するセンサ姿勢角度θ１、θ２の推定の可能となる。
【００９９】
以上、図面を参照してこの発明の実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。例えば、上述の第１から第３の実施形態において説明した各構成は、任意に組み合わせることができる。
【０１００】
また、上述した実施形態では、マイク１６の数が２個の場合、または少なくとも４個以上の場合等、複数個設置される例について説明したが、スピーカー１５の数も１個に限られるものではなく、複数個としてもよい。また、スピーカー１５は、フェーズドアレイスピーカーとしてもよい。
【０１０１】
また、上述した実施形態では、音波が反射する対象物が構造物等の壁面である例を説明したが、対象物は、音波が反射する物であれば、どのようなものであってもよい。
【０１０２】
また、計測装置１４１は、計測した風に関する情報（風速、風向、壁面との距離、または壁面に対するセンサ角度等）を、通信部１１を介して他の装置に送ってもよい。送信先は、通信ネットワークを介して接続されるサーバ装置であってもよいし、ドローン１０を遠隔から操縦するためのコントローラであってもよい。これにより、ドローン１０の周辺の風に関する情報を、ドローン１０から離れた場所でも確認または利用することができる。
【０１０３】
なお、上述した実施形態における制御部１４の一部または全部をコンピュータで実現するようにしても良い。その場合、この制御部１４の一部または全部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現しても良い。なお、ここでいう「コンピュータシステム」とは、ドローン１０に内蔵されたコンピュータシステムであって、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでも良い。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
【０１０４】
また、上述した実施形態における制御部１４の一部または全部を、ＬＳＩ（Ｌａｒｇｅ  Ｓｃａｌｅ  Ｉｎｔｅｇｒａｔｉｏｎ）等の集積回路として実現してもよい。また、例えば、制御部１４の一部または全部を集積してプロセッサ化してもよい。また、集積回路化の手法はＬＳＩに限らず専用回路、または汎用プロセッサで実現してもよい。また、半導体技術の進歩によりＬＳＩに代替する集積回路化の技術が出現した場合、当該技術による集積回路を用いてもよい。
【０１０５】
また、上述した実施形態では、計測装置１４１がドローン１０に備えられている構成を説明したが、計測装置１４１は、ドローン１０とは異なる装置に備えられてもよい。例えば、計測装置１４１は、ドローン１０を遠隔から操縦するためのコントローラに備えられてもよいし、或いは、通信ネットワークを介して接続されるサーバ装置に備えられてもよい。これらの場合、ドローン１０は、マイク１６が受信した受信結果を、上記コントローラ或いはサーバ装置へ送信することにより、上記コントローラ或いはサーバ装置で風に関する情報を推定し、推定結果がドローン１０へフィードバックされてもよい。なお、このとき、スピーカー１５の出力は、ドローン１０が自律的に制御してもよいし、上記コントローラ或いはサーバ装置から制御信号をドローン１０へ送信することにより、該制御信号に基づいてドローン１０が制御してもよい。
【０１０６】
また、上述した実施形態では、計測装置１４１をドローン１０に適用する例を説明したが、ドローン１０に限らず、いずれの飛行体に適用してもよいし、また飛行する物に限らず、車等の移動体に適用してもよい。
【符号の説明】
【０１０７】
１０  ドローン、１１  通信部、１２  駆動部、１３  記憶部、１４  制御部、１５  スピーカー、１６  マイク、１４１  計測装置、１４１１  音波反射時間計測部、１４１２  推定部、１４２  駆動制御部
【図１】
2019095219.tif 000019
【図２】
2019095219.tif 000020
【図３】
2019095219.tif 000021
【図４】
2019095219.tif 000022
【図５】
2019095219.tif 000023
【図６】
2019095219.tif 000024
【図７】
2019095219.tif 000025
【図８】
2019095219.tif 000026
【図９】
2019095219.tif 000027
【図１０】
2019095219.tif 000028
【図１１】
2019095219.tif 000029
【図１２】
2019095219.tif 000030
【図１３】
2019095219.tif 000031
【図１４】
2019095219.tif 000032
【図１５】
2019095219.tif 000033
【図１６】
2019095219.tif 000034
【図１７】
2019095219.tif 000035
【図１８】
2019095219.tif 000036