(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022076378
(43)【公開日】2022-05-19
(54)【発明の名称】距離計測装置、距離計測方法、及びプログラム
(51)【国際特許分類】
G01S 13/34 20060101AFI20220512BHJP
【FI】
G01S13/34
【審査請求】有
【請求項の数】7
【出願形態】OL
(21)【出願番号】P 2020186767
(22)【出願日】2020-11-09
(71)【出願人】
【識別番号】503361400
【氏名又は名称】国立研究開発法人宇宙航空研究開発機構
(74)【代理人】
【識別番号】100106909
【弁理士】
【氏名又は名称】棚井 澄雄
(74)【代理人】
【識別番号】100161702
【弁理士】
【氏名又は名称】橋本 宏之
(74)【代理人】
【識別番号】100188592
【弁理士】
【氏名又は名称】山口 洋
(74)【代理人】
【識別番号】100181124
【弁理士】
【氏名又は名称】沖田 壮男
(74)【代理人】
【識別番号】100163496
【弁理士】
【氏名又は名称】荒 則彦
(72)【発明者】
【氏名】瀬在 俊浩
【テーマコード(参考)】
5J070
【Fターム(参考)】
5J070AB08
5J070AB24
5J070AC02
5J070AD01
5J070AH04
5J070AK22
(57)【要約】
【課題】計測精度、分解能を送信電波の周波数変化幅から決まる値より向上できる距離計測装置、距離計測方法、及びプログラムを提供することを目的とする。
【解決手段】距離計測装置は、周波数変調させたパルス波をターゲットに向けて送信し、ターゲットで反射されて戻って来たエコー波を送信したパルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部、を備える。
【選択図】
図1
【特許請求の範囲】
【請求項1】
周波数変調したパルス波をターゲットに向けて送信し、前記ターゲットで反射されて戻って来たエコー波を送信した前記パルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部、
を備える距離計測装置。
【請求項2】
前記信号処理部は、前記フィッティングにより得られた指数関数のうち、前記実数部の絶対値が所定値以下の指数関数を抽出する、
請求項1に記載の距離計測装置。
【請求項3】
前記信号処理部は、前記抽出した指数関数から振幅の値、及び前記虚数部の値を取得する、
請求項2に記載の距離計測装置。
【請求項4】
前記信号処理部は、前記取得した振幅の値、及び前記虚数部の値を含む指数関数情報から、前記ターゲットの相対振幅の値、及び距離の値を取得する、
請求項3に記載の距離計測装置。
【請求項5】
前記取得した相対振幅の値、及び距離の値を含むターゲット情報を表示する表示部を更に備える、
請求項4に記載の距離計測装置。
【請求項6】
コンピュータが、
周波数変調したパルス波をターゲットに向けて送信し、
前記ターゲットで反射されて戻って来たエコー波を送信した前記パルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする、
距離計測方法。
【請求項7】
コンピュータに、
周波数変調したパルス波をターゲットに向けて送信させ、
前記ターゲットで反射されて戻って来たエコー波を送信した前記パルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングさせる、
プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、距離計測装置、距離計測方法、及びプログラムに関する。
【背景技術】
【0002】
パルス圧縮レーダは、周波数が時間に対して直線的に変化するように変調させたパルス電波をターゲットに向けて送信する。パルス圧縮レーダは、ターゲットで反射されてレーダに戻って来たエコー電波をパルス時間幅毎に送信したパルス電波と逆相関処理を行うと、レーダとターゲットの距離の2倍に対応した時刻がピークとなる信号が得られる。パルス圧縮レーダは、ピークとなる時刻を検出すると、ターゲットの距離を計測することができる。以上がパルス圧縮レーダでの距離計測の原理である。この時、周波数変化幅をΔfとすると、距離分解能Δdは、{c/(2Δf)}で与えられる。ここでcは光速である(例えば非特許文献1参照)。なお、パルス圧縮レーダの距離分解能は、非特許文献1の式(1.14)Rres=cτ/2のパルス幅τに、パルス圧縮後の波形のパルス幅の1/Δfを代入すれば求めることができる。
【先行技術文献】
【非特許文献】
【0003】
【非特許文献1】吉田孝、“改訂 レーダ技術”、社団法人電子情報通信学会編、コロナ社、平成11年5月25日、p10-13、p274-277
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、従来技術では、ターゲットの距離が、非特許文献1も記載のように距離分解能単位での計測となる。このため、従来技術では、距離計測精度、分解能がパルス圧縮レーダの周波数変化幅から定まる値に制限されるという課題があった。
【0005】
本発明は、上記の問題点に鑑みてなされたものであって、計測精度、分解能を送信電波の周波数変化幅から決まる値より向上できる距離計測装置、距離計測方法、及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【0006】
上記目的を達成するため、本発明の一態様に係る距離計測装置は、周波数変調したパルス波をターゲットに向けて送信し、前記ターゲットで反射されて戻って来たエコー波を送信した前記パルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部、を備える。
【0007】
また、本発明の一態様に係る距離計測装置において、前記信号処理部は、前記フィッティングにより得られた指数関数のうち、前記実数部の絶対値が所定値以下の指数関数を抽出するようにしてもよい。
【0008】
また、本発明の一態様に係る距離計測装置において、前記信号処理部は、前記抽出した指数関数から振幅の値、及び前記虚数部の値を取得するようにしてもよい。
【0009】
また、本発明の一態様に係る距離計測装置において、前記信号処理部は、前記取得した振幅の値、及び前記虚数部の値を含む指数関数情報から、前記ターゲットの相対振幅の値、及び距離の値を取得するようにしてもよい。
【0010】
また、本発明の一態様に係る距離計測装置において、前記取得した相対振幅の値、及び距離の値を含むターゲット情報を表示する表示部を更に備えるようにしてもよい。
【0011】
上記目的を達成するため、本発明の一態様に係る距離計測方法は、コンピュータが、周波数変調したパルス波をターゲットに向けて送信し、前記ターゲットで反射されて戻って来たエコー波を送信した前記パルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする。
【0012】
上記目的を達成するため、本発明の一態様に係るプログラムは、コンピュータに、周波数変調したパルス波をターゲットに向けて送信させ、前記ターゲットで反射されて戻って来たエコー波を送信した前記パルス波と周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングさせる。
【発明の効果】
【0013】
本発明によれば、計測精度、分解能を送信電波の周波数変化幅から決まる値より向上できる。
【図面の簡単な説明】
【0014】
【
図1】実施形態に係る距離計測装置の構成例を示すブロック図である。
【
図2】実施形態に係る信号処理部の構成例を示すブロック図である。
【
図3】実施形態に係る距離測定の処理手順のフローチャートである。
【
図4】シミュレーションによって得られたターゲット情報の内、相対振幅が上位5までの情報を示した図である。
【
図5】本実施形態におけるシミュレーションによって得られた距離と相対強度の関係を示す図である。
【
図6】従来の手法によりシミュレーションを行った結果を示す図である。
【発明を実施するための形態】
【0015】
以下、本発明の実施の形態について図面を参照しながら説明する。なお、以下の例では、距離計測装置をレーダ装置に適用する例を説明する。
【0016】
<距離計測装置の構成例>
図1は、本実施形態に係る距離計測装置の構成例を示すブロック図である。
図1のように、距離計測装置1は、発振部11、遅延回路12、サーキュレータ部13、アンテナ14、逆相関部15、信号処理部16、および表示部17を備える。
【0017】
距離計測装置1は、例えばパルス圧縮レーダ装置であり、ターゲットに向けてパルス電波(「パルス波」ともいう)を送信し、ターゲットで反射されてレーダに戻って来たエコー電波(「エコー波」ともいう)を送信したパルス電波と周波数領域で逆相関して得られる信号を処理することにより、ターゲットの距離を計測する。
【0018】
発振部11は、周波数が時間に対して直線的に変化するように変調させたパルス電波を、遅延回路12とサーキュレータ部13へ出力する。
【0019】
遅延回路12は、発振部11から入力されたパルス電波をパルス時間の0倍を含む整数倍の時間を遅らせ、逆相関部15へ出力する。
【0020】
サーキュレータ部13は、発振部11から入力された電波をアンテナ14へ出力する。また、サーキュレータ部13は、アンテナ14から入力された電波を逆相関部15へ出力する。
【0021】
アンテナ14は、サーキュレータ部13から入力された電波をターゲットへ向けて放射するとともに、ターゲットで散乱されて戻ってきたエコー電波を受信し、サーキュレータ部13へ出力する。アンテナ14は、例えば、ダイポール系のアンテナ、ホーンアンテナ、パラボラアンテナ、アレイアンテナなど如何なるアンテナでも用いることができる。
【0022】
逆相関部15は、遅延回路12から入力されたパルス電波とサーキュレータ部13から入力された電波を周波数領域で逆相関して得られる信号を信号処理部16へ出力する。
【0023】
信号処理部16は、所定の信号処理を施し、ターゲットの相対振幅の値、及び距離の値を表示部17へ出力する。
【0024】
表示部17は、例えば液晶表示装置や有機EL(Electro Luminescence)表示装置等であり、信号処理部16からの出力に基づきターゲットの相対強度及び距離に関する情報を表示する。
【0025】
<信号処理部16の構成例>
次に、信号処理部16の構成例を説明する。
図2は、本実施形態に係る信号処理部の構成例を示すブロック図である。
図2のように、信号処理部16は、指数関数フィッティング部161、指数関数抽出部162、指数関数情報取得部163、およびターゲット情報取得部164を備える。
【0026】
指数関数フィッティング部161は、入力された信号を引数に実数部と虚数部とを有する指数関数でプロニー(Prony)法によりフィッティングする。なお、プロニー法は、データ系列を指数関数でフィッティングする手法である。
【0027】
指数関数抽出部162は、フィッティングにより得られた指数関数のうち、前記実数部の絶対値が所定値以下の指数関数を抽出する。
【0028】
指数関数情報取得部163は、抽出した指数関数から振幅の値、前記虚数部の値を取得する。
【0029】
ターゲット情報取得部164は、取得した指数関数情報から、ターゲットの相対振幅の値、及び距離の値を取得する。
【0030】
<距離測定の方法>
次に、距離測定方法例を説明する。
パルス圧縮レーダでは、送信した電波の周波数が時間に対して直線的に変化するように変調させたパルス電波をターゲットに向けて送信する。送信電波信号eT(t)は以下で表現される。
【0031】
【0032】
式(1)において、jは虚数単位、tは時間、f0はキャリア周波数、cは光速であり、Δfは周波数変化幅であり、TCはパルス時間幅である。ここでは説明を簡略化するために、ターゲットはレーダからcTC/2の距離内に存在する場合を考える。
レーダからターゲットまでの距離をRkとすると、ターゲットのエコー信号は、2Rk/cの時間後にレーダに到達するので、エコー信号eRk(t)は、次式(2)で表現される。
【0033】
【0034】
式(2)において、σkはターゲットの散乱係数である。また、式(2)はターゲットが1個存在している場合のエコー信号である。これを一般化し、p個のターゲットが存在している場合のエコー信号eRk(t)は、次式(3)となる。
【0035】
【0036】
エコー信号と送信信号の逆相関処理は、これらの信号をベースバンドに変換した信号で実施する。送信信号eT(t)をベースバンド信号bT(t)に変換した式は、次式(4)となる。
【0037】
【0038】
エコー信号eRk(t)をベースバンド信号bRT(t)に変換した式は、次式(5)となる。
【0039】
【0040】
逆相関処理はスペクトル領域では除算となるため、効率的に実施できる。エコー信号bRk(t)は、ベースバンド信号bT(t)の時間を2Rk/cずらし、係数を乗じた信号である。このため、ベースバンド信号bT(t)のスペクトル信号をBT(ω)とすると、bRk(t)のスペクトル信号BRk(ω)はBT(ω)を使用して次式(6)で表現される。
【0041】
【0042】
逆相関処理は、スペクトル領域でBRk(ω)をBT(ω)で割ることで実施される。BRk(ω)をBT(ω)で割った信号は、次式(7)で表現される。
【0043】
【0044】
この信号をフーリエ変換で時間領域に再変換すると逆相関処理が完了する。これが従来の逆相関処理である。
【0045】
式(7)から分かるように、スペクトル領域で逆相関処理をして得られる信号は、指数関数の和となる。実際に得られる信号は数値で構成されるデータ系列である。データ系列を指数関数にフィッティングする方法としてプロニー法が知られている。このため、本実施形態では、スペクトル領域で逆相関処理をして得られるデータ系列をプロニー法で指数関数にフィッティングすることで、指数関数の振幅σk・exp(-j4πf0Rk/c)と引数(-j2Rk)/cを求める。これにより、本実施形態では、ターゲットの距離Rkを求めることができる。
【0046】
この方法では、従来の方法のように帯域Δfの制約を受けることなく、ターゲットの距離Rkを数値で直接求めるため、精確に求めることができる。
式(7)は理想的な場合の式であるが、実際に得られる信号にはノイズが含まれるため、それをプロニー法でフィッティングすると、BRk(ω)をBT(ω)で割った信号{・^}は、次式(8)で表現される信号が得られる。
【0047】
【0048】
実際の信号をプロニー法でフィッティングすると、αk・exp(-j4πf0Rk/c)と(-j2Rk)/cの他にαkも求められる。式(7)に示したように、理想的な状態ではターゲットに対するαkは0であるので、ノイズが含まれる信号をプロニー法でフィッティングした場合、ターゲットに対するαkは0ではないとしても、小さい値であると考えられる。
【0049】
ノイズに対するαkは、小さい値であるとは限らず、任意の値であると考えられる。αkの大きさのこの違いを利用すると、得られた指数関数の中から、ノイズに対する指数関数の一部を取り除くことで、ターゲットに対する指数関数を効率よく抽出することが可能となる。
【0050】
<距離計測方法>
次に、距離測定の処理手順例を説明する。
図3は、本実施形態に係る距離測定の処理手順のフローチャートである。
【0051】
(ステップS1)アンテナ14は、サーキュレータ部13から入力された電波をターゲットへ向けて放射する。
【0052】
(ステップS2)アンテナ14は、ターゲットで散乱されて戻ってきたエコー電波を受信する。
【0053】
(ステップS3)遅延回路12は、発振部11から入力されたパルス電波をパルス時間の0倍を含む整数倍の時間を遅らせる。続けて、逆相関部15は、遅延回路12から入力されたパルス電波とサーキュレータ部13から入力された電波に対して周波数領域で逆相関を行う。
【0054】
(ステップS4)指数関数フィッティング部161は、逆相関処理で得られる信号をプロニー法によりフィッティングする。
【0055】
(ステップS5)指数関数抽出部162は、フィッティングにより得られた指数関数のうち、引数の実数部の絶対値が所定値以下の指数関数を抽出する。
【0056】
(ステップS6)指数関数情報取得部163は、抽出された指数関数から振幅の値、引数の虚数部の値を取得する。
【0057】
(ステップS7)ターゲット情報取得部164は、取得された指数関数情報から、ターゲットの相対振幅の値、及び距離の値を取得する。
【0058】
なお、上述した例では電波を送受信する例を説明したが、これに限らない。距離計測装置が送受信する対象は電波に限らず、光、音波等であってもよい。このように、光や音波を用いて送受信する装置であっても、パルス圧縮を行う装置に、本実施形態の手法を適用することができる。すなわち、「周波数変調したパルス波」は、電波、光および音波等によるものである。
【0059】
<評価>
次に、本実施形態の手法の効果を確認するために行った数値シミュレーション結果例を説明する。
シミュレーション条件は、キャリア周波数が76[GHz]で、1[μ秒]の間に周波数が直線的に500[MHz]増加するように変調させた電波をパルス圧縮レーダから送信する。そして、シミュレーションでは、レーダから20.0[m]と20.2[m]の距離に存在しているレーダ散乱断面積が等しい2つのターゲットから散乱され、レーダに戻ってくるエコー電波を信号対雑音比が20[dB]の環境で取得した際に得られるターゲット情報を検討した。
【0060】
中間周波数信号を引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングして得られた指数関数の内、引数の実数部の絶対値が0.003以下となる指数関数を抽出し、指数関数の振幅値からターゲットの相対振幅を取得し、指数関数の引数の虚数部からターゲットの距離を取得した結果を
図4に示す。
図4は、シミュレーションによって得られたターゲット情報の内、相対振幅が上位5までの情報を示した図である。
【0061】
シミュレーションで得られたターゲット情報の内、距離が10~30[m]までのターゲット情報を図示した結果を
図5に示す。
図5は、本実施形態におけるシミュレーションによって得られた距離と相対強度の関係を示す図である。
図5において、横軸は距離[m]であり、縦軸は相対強度[dB]である。
【0062】
図4より、相対振幅が上位2までと上位3以下とでは、相対振幅が20[dB]以上異なることより、上位2までがターゲットの情報で、上位3以下は雑音による情報であることが分かる。従って、ターゲットは2つ存在していることが分かる。2つのターゲット情報の内の距離は、実際の距離との誤差が4[cm]、9[cm]と非常に小さい。また、2つのターゲットの相対振幅の差は0.48[dB]と非常に小さい。
【0063】
シミュレーションで使用した周波数変化幅は500[MHz]であるので、従来の手法で得られる距離計測の精度、分解能は30[cm]である。これに対して、本実施形態によれば、
図4、
図5のように30[cm]以下に接近している2つのターゲットを識別することができている。また、得られた2つのターゲットの距離の誤差は最大で9[cm]であり、従来の手法の距離計測精度の3分の1未満である。これらより、本実施形態により、パルス圧縮レーダシステムの周波数変化幅から定まる値に制限されることなく、ターゲットの数、距離、及び相対振幅を求められる。
【0064】
ここで、比較例を説明する。
図6は、従来の手法によりシミュレーションを行った結果を示す図である。
図6において、横軸は距離[m]であり、縦軸は相対強度[dB]である。なお、シミュレーション条件は、本実施形態の手法の効果を確認するために行った数値シミュレーションと同じである。
図6のように、従来の手法では、ターゲットは20.0[m]付近に存在することが分かるのみで、ターゲットが2つあることを認識することができない。
【0065】
以上のように、本実施形態では、エコー電波と送信電波を周波数領域で逆相関して得られる信号を引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングし、フィッティングにより得られた指数関数の引数の実数部の絶対値が所定値以下の指数関数を抽出する。そして、本実施形態では、抽出された指数関数から振幅の値、引数の虚数部の値を取得し、取得された指数関数情報から、ターゲットの相対振幅の値、及び距離の値を取得するようにした。
【0066】
これにより、本実施形態によれば、計測精度、分解能を送信電波の周波数変化幅から決まる値より向上できる。
【0067】
なお、本発明における距離計測装置1の機能の全てまたは一部を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより距離計測装置1が行う処理の全てまたは一部を行ってもよい。なお、ここでいう「コンピュータシステム」とは、OSや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ホームページ提供環境(あるいは表示環境)を備えたWWWシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ROM、CD-ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ(RAM)のように、一定時間プログラムを保持しているものも含むものとする。
【0068】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク(通信網)や電話回線等の通信回線(通信線)のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であってもよい。
【0069】
以上、本発明を実施するための形態について実施形態を用いて説明したが、本発明はこうした実施形態に何等限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々の変形および置換を加えることができる。
【符号の説明】
【0070】
1…距離計測装置、11…発振部、12…遅延回路、13…サーキュレータ部、14…アンテナ、15…逆相関部、16…信号処理部、17…表示部、161…指数関数フィッティング部、162…指数関数抽出部、163…指数関数情報取得部、164…ターゲット情報取得部
【手続補正書】
【提出日】2022-03-16
【手続補正1】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項1】
周波数変調したパルス波をターゲットに向けて送信し、前記ターゲットで反射されて戻って来たエコー波と送信した前記パルス波それぞれをベースバンドに変換した信号に対して周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部、
を備える距離計測装置。
【請求項2】
前記信号処理部は、前記フィッティングにより得られた指数関数のうち、前記実数部の絶対値が所定値以下の指数関数を抽出する、
請求項1に記載の距離計測装置。
【請求項3】
前記信号処理部は、前記抽出した指数関数から振幅の値、及び前記虚数部の値を取得する、
請求項2に記載の距離計測装置。
【請求項4】
前記信号処理部は、前記取得した振幅の値、及び前記虚数部の値を含む指数関数情報から、前記ターゲットの相対振幅の値、及び距離の値を取得する、
請求項3に記載の距離計測装置。
【請求項5】
前記取得した相対振幅の値、及び距離の値を含むターゲット情報を表示する表示部を更に備える、
請求項4に記載の距離計測装置。
【請求項6】
コンピュータが、
周波数変調したパルス波をターゲットに向けて送信し、
前記ターゲットで反射されて戻って来たエコー波と送信した前記パルス波それぞれをベースバンドに変換した信号に対して周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする、
距離計測方法。
【請求項7】
コンピュータに、
周波数変調したパルス波をターゲットに向けて送信させ、
前記ターゲットで反射されて戻って来たエコー波と送信した前記パルス波それぞれをベースバンドに変換した信号に対して周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングさせる、
プログラム。
【手続補正2】
【補正対象書類名】明細書
【補正対象項目名】0006
【補正方法】変更
【補正の内容】
【0006】
上記目的を達成するため、本発明の一態様に係る距離計測装置は、周波数変調したパルス波をターゲットに向けて送信し、前記ターゲットで反射されて戻って来たエコー波と送信した前記パルス波それぞれをベースバンドに変換した信号に対して周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする信号処理部、を備える。
【手続補正3】
【補正対象書類名】明細書
【補正対象項目名】0011
【補正方法】変更
【補正の内容】
【0011】
上記目的を達成するため、本発明の一態様に係る距離計測方法は、コンピュータが、周波数変調したパルス波をターゲットに向けて送信し、前記ターゲットで反射されて戻って来たエコー波と送信した前記パルス波それぞれをベースバンドに変換した信号に対して周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングする。
【手続補正4】
【補正対象書類名】明細書
【補正対象項目名】0012
【補正方法】変更
【補正の内容】
【0012】
上記目的を達成するため、本発明の一態様に係るプログラムは、コンピュータに、周波数変調したパルス波をターゲットに向けて送信させ、前記ターゲットで反射されて戻って来たエコー波と送信した前記パルス波それぞれをベースバンドに変換した信号に対して周波数領域で逆相関して得られる信号を、引数に実数部と虚数部とを有する指数関数でプロニー法によりフィッティングさせる。