(57)【要約】【課題】光パルス計測装置及びその計測方法を提供する。【解決手段】光パルスの計測方法はパルス光源から発せられる光パルスの計測に用いられる。前記方法は、光パルスが分光された後に再度測定位置に集光され、分光後の光パルスが測定位置で自己相関効果を発生させると共にガスプラズマを発生させる工程と、ガスプラズマから発生する音声を収音させ、プラズマ音声信号が生成される工程と、プラズマ音声信号に基づいて光パルスの特性が計算される工程とを含む。その計測装置も提出する。【選択図】図１
2019203874.tif 000002
【特許請求の範囲】
【請求項１】
パルス光源から発せられる光パルスの計測に用いられる光パルス計測方法であって、
前記光パルスが分光された後に再度測定位置に集光され、前記分光後の光パルスが前記測定位置で自己相関効果を発生させると共にガスプラズマを発生させる工程と、
前記ガスプラズマから発生する音声が収音され、プラズマ音声信号が生成される工程と、
前記プラズマ音声信号に基づいて前記光パルスの特性が計算される工程とを含むことを特徴とする光パルス計測方法。
【請求項２】
前記プラズマ音声信号は測定チャートを含み、前記測定チャートは時間軸に沿って記録される前記音声の強度であり、且つ前記プラズマ音声信号に基づいて前記光パルスの特性が計算される工程は、
前記測定チャートに基づいて前記光パルスの特性が計算される工程を含むことを特徴とする請求項１に記載の光パルス計測方法。
【請求項３】
前記測定チャートに基づいて前記光パルスの特性が計算される工程は、
フィット曲線により前記測定チャートにフィットさせる工程と、
前記フィット曲線により前記光パルスの特性が計算される工程とを含むことを特徴とする請求項２に記載の光パルス計測方法。
【請求項４】
前記フィット曲線は、
【数１】
2019203874.tif 000016
を含み、αは前記測定チャートに対応させる振幅の定数であり、
ここでは、
【数２】
2019203874.tif 000017
であり、ｋ１は定数であり、
ここでは、
【数３】
2019203874.tif 000018
であり、ｋ２は定数であり、且つ前記測定チャート中の前記ガスプラズマの音声に対応する信号は（ａ−ｋ１）から（ｂ＋ｋ２）の間の範囲に位置されることを特徴とする請求項３に記載の光パルス計測方法。
【請求項５】
前記フィット曲線により前記光パルスの特性が計算される工程は、
前記フィット曲線の計測幅により前記光パルスの時間軸上における幅が計算され、前記計測幅は（ｂ−ａ）である工程を更に含むことを特徴とする請求項４に記載の光パルス計測方法。
【請求項６】
前記ｆＬ（ｔ）関数及び前記ｆＲ（ｔ）関数は、Error function、Step function、Heaviside step function、Fermi-Dirac distribution function、Maxwell-Boltzmann distribution、Boltzmann distribution、Sigmoid function、Smoothstep function、Generalized logistic function、Richards’ curve、Logistic function、Logistic curve、Cumulative distribution function、percent point function、Quantile function、Survival function、またはReliability functionであることを特徴とする請求項４に記載の光パルス計測方法。
【請求項７】
前記光パルスが前記測定位置で分散される場合、前記フィット曲線は、
−ｇ（ｔ）を更に含み、ここでは、−ｇ（ｔ）は前記プラズマ音声信号図中の分散現象が発生する部分図形に対応することを特徴とする請求項４に記載の光パルス計測方法。
【請求項８】
前記ｇ（ｔ）はガウス関数であることを特徴とする請求項７に記載の光パルス計測方法。
【請求項９】
前記測定チャートに基づいて前記光パルスの特性が計算される工程は、
前記測定チャートの幅を利用して前記光パルスの時間軸上における幅が計算される工程を更に含むことを特徴とする請求項２に記載の光パルス計測方法。
【請求項１０】
前記光パルスが前記測定位置で分散される場合、前記測定チャートに基づいて前記光パルスの特性が計算される工程は、分散補正値及び前記測定チャートを利用して前記光パルスの特性が計算される工程を更に含むことを特徴とする請求項2に記載の光パルス計測方法。
【請求項１１】
パルス光源から発せられる光パルスの計測に用いられる光パルス計測装置であって、
前記光パルスが計測位置で発生させる自己相関効果を利用してガスプラズマを発生させる光モジュールと、
ガスプラズマから発生する音声を収音させる収音ユニットと、
前記収音ユニットに接続される処理ユニットとを備え、
ここでは、前記収音ユニットはガスプラズマから発生する音声をプラズマ音声信号に変換させて前記処理ユニットに伝送させ、
前記処理ユニットは前記プラズマ音声信号に基づいて前記光パルスの特性の計算を行うことを特徴とする光パルス計測装置。
【請求項１２】
前記プラズマ音声信号は測定チャートを含み、前記測定チャートは時間軸に沿って記録される前記音声の強度であり、前記処理ユニットは前記測定チャートに基づいて前記光パルスの光学特性の計算を行うことを特徴とする請求項１１に記載の光パルス計測装置。
【請求項１３】
前記処理ユニットはフィット曲線により前記測定チャートにフィットさせ、且つ前記フィット曲線により前記光パルスの特性の計算を行うことを特徴とする請求項１２に記載の光パルス計測装置。
【請求項１４】
前記フィット曲線は、
【数４】
2019203874.tif 000019
を含み、αは前記測定チャートに対応させる振幅の定数であり、
ここでは、
【数５】
2019203874.tif 000020
であり、ｋ１は定数であり、
ここでは、
【数６】
2019203874.tif 000021
であり、ｋ２は定数であり、且つ前記測定チャート中の前記ガスプラズマの音声に対応する信号は（ａ−ｋ１）から（ｂ＋ｋ２）の間の範囲に位置されることを特徴とする請求項１３に記載の光パルス計測装置。
【請求項１５】
前記処理ユニットは前記フィット曲線の計測幅を利用して前記光パルスの時間軸上における幅の計算を行い、且つ前記計測幅は（ｂ−ａ）であることを特徴とする請求項１４に記載の光パルス計測装置。
【請求項１６】
前記ｆＬ（ｔ）関数及び前記ｆＲ（ｔ）関数は、Error function、Step function、Heaviside step function、Fermi-Dirac distribution function、Maxwell-Boltzmann distribution、Boltzmann distribution、Sigmoid function、Smoothstep function、Generalized logistic function、Richards’ curve、Logistic function、Logistic curve、Cumulative distribution function、percent point function、Quantile function、Survival function、またはReliability functionであることを特徴とする請求項１４に記載の光パルス計測装置。
【請求項１７】
前記光パルスが前記測定位置で分散される場合、前記フィット曲線は、
−ｇ（ｔ）を更に含み、ここでは、−ｇ（ｔ）は前記プラズマ音声信号図中の分散現象が発生する図形に対応することを特徴とする請求項１４に記載の光パルス計測装置。
【請求項１８】
前記ｇ（ｔ）はガウス関数であることを特徴とする請求項１７に記載の光パルス計測装置。
【請求項１９】
前記光モジュールは、
第一分光ユニットと、
第一反射ユニットと、
第二反射ユニットであって、前記第一分光ユニットにより前記光パルスが前記第一反射ユニット及び前記第二反射ユニットに分光されることと、
前記第一反射ユニット及び前記第二反射ユニットにより反射される光パルスを併合させる第二分光ユニットと、
前記第二分光ユニットからの光パルスを前記計測位置に集光させ、且つ前記計測位置でガスプラズマを発生させるフォーカスユニットとを備えることを特徴とする請求項１１に記載の光パルス計測装置。
【請求項２０】
前記第一分光ユニット及び前記第二分光ユニットはビームスプリッターであり、前記フォーカスユニットは軸外放物面鏡であり、前記収音ユニットはマイクロフォンであることを特徴とする請求項１９に記載の光パルス計測装置。
【請求項２１】
遮音空間を有し、少なくとも前記計測位置及び前記収音ユニットが前記遮音空間内に位置される遮音カバーとを更に備えることを特徴とする請求項１１に記載の光パルス計測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、計測装置または計測方法に関し、より詳しくは、フォトガスプラズマの計測装置または計測方法に関する。
【背景技術】
【０００２】
レーザーの発展以来、高出力レーザーは重要な発展課題であった。一般的な連続波レーザー（Continuous wave laser）より、パルスレーザー（Pulsed Laser）は高エネルギーを提供しやすい。パルスレーザーの発展過程において多くの人がより短いパルス幅を追い求めるようになっており、このため、パルス幅の計測を行う自己相関計（Autocorrelator）がこの技術分野では鍵となる道具となっている。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、上述の自己相関計は計測時に光源の強度を落とさなければならないため、高エネルギーのレーザーパルスを計測できず、もし計測したとしても自己相関計の非線形結晶が高エネルギーのレーザーパルスによって破壊されてしまう。また、非線形結晶自体のコストも高く、自己相関計の核心部分の部材が計測時にレーザービームにより破壊されて計測を行えなくなれば更なる損害を被る。このため、高エネルギーレーザーのパルスを直接計測可能な計測装置及び方法によってレーザーパルスの幅を計測することが求められた。
【０００４】
本発明は、低エネルギー及び高エネルギーのレーザー光のパルスの光学的性質を計測可能な光パルス計測装置及びその計測方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記課題を解決するために、本発明のある態様の光パルス計測方法は、パルス光源から発せられる光パルスの計測に用いられ、光パルスが分光された後に再度測定位置に集光され、分光後の光パルスが測定位置で自己相関効果を発生させると共にガスプラズマを発生させる工程と、ガスプラズマから発生する音声を収音させ、プラズマ音声信号が生成される工程と、プラズマ音声信号に基づいて光パルスの特性が計算される工程とを含む。
【０００６】
一実施形態において、上述のプラズマ音声信号は測定チャートを含む。測定チャートは時間軸に沿って記録される音声の強度であり、且つ上述の方法では測定チャートに基づいて光パルスの特性が計算される。
【０００７】
一実施形態において、上述の方法はフィット曲線により測定チャートにフィットさせ、且つフィット曲線により光パルスの特性が計算される。
【０００８】
一実施形態において、上述のフィット曲線は、
【数１】
2019203874.tif 000003
を含み、αは測定チャートに対応させる振幅の定数であり、
ここでは、
【数２】
2019203874.tif 000004
であり、ｋ１は定数であり、
ここでは、
【数３】
2019203874.tif 000005
であり、ｋ２は定数であり、且つ測定チャート中のガスプラズマの音声に対応する信号は（ａ−ｋ１）から（ｂ＋ｋ２）の間の範囲に位置される。
【０００９】
一実施形態において、上述の方法では、フィット曲線の計測幅を利用して光パルスの時間軸上における幅の計算を行い、計測幅は（ｂ−ａ）である。
【００１０】
一実施形態において、上述のｆＬ（ｔ）関数及びｆＲ（ｔ）関数は、Error function、Step function、Heaviside step function、Fermi-Dirac distribution function、Maxwell-Boltzmann distribution、Boltzmann distribution、Sigmoid function、Smoothstep function、Generalized logistic function、Richards’ curve、Logistic function、Logistic curve、Cumulative distribution function、percent point function、Quantile function、Survival function、またはReliability functionである。
【００１１】
一実施形態において、分光後の光パルスが測定位置で分散される場合、フィット曲線は−ｇ（ｔ）を更に含み、ここでは、−ｇ（ｔ）はプラズマ音声信号図中の分散現象が発生する図形に対応する。
【００１２】
一実施形態において、上述のｇ（ｔ）はガウス関数である。
【００１３】
一実施形態において、上述の方法では、測定チャートの幅を利用して光パルスの時間軸上における幅の計算が行われる。
【００１４】
一実施形態において、分光後の光パルスが測定位置で分散される場合、フィット曲線の計測幅を利用して光パルスの時間軸上における幅の計算が行われる工程では、
分散補正値及び計測幅を利用して光パルスの時間軸上における幅の計算が行われる工程を更に含む。
【００１５】
本発明に係る計測装置はパルス光源から発せられる光パルスの計測に用いられる。計測装置は、
光パルスが計測位置で発生させる自己相関効果を利用してガスプラズマを発生させる光モジュールと、
ガスプラズマから発生する音声を収音させる収音ユニットと、
収音ユニットに接続される処理ユニットとを備える。
ここでは、収音ユニットはガスプラズマから発生する音声をプラズマ音声信号に変換させて処理ユニットに伝送させる。
処理ユニットはプラズマ音声信号に基づいて光パルスの特性の計算を行う。
【００１６】
一実施形態において、上述のプラズマ音声信号は測定チャートを含み、測定チャートは時間軸に沿って記録される音声の強度である。処理ユニットは測定チャートに基づいて光パルスの光学的特性の計算を行う。
【００１７】
一実施形態において、上述の処理ユニットはフィット曲線により測定チャートにフィットさせ、且つフィット曲線により光パルスの特性の計算を行う。
【００１８】
一実施形態において、上述のフィット曲線は、
【数４】
2019203874.tif 000006
を含み、αは測定チャートに対応させる振幅の定数である。
ここでは、
【数５】
2019203874.tif 000007
であり、ｋ１は定数であり、
ここでは、
【数６】
2019203874.tif 000008
であり、ｋ２は定数であり、且つ測定チャート中のガスプラズマの音声に対応する信号は（ａ−ｋ１）から（ｂ＋ｋ２）の間の範囲に位置される。
【００１９】
一実施形態において、上述の処理ユニットはフィット曲線の計測幅を利用して光パルスの時間軸上における幅の計算を行い、且つ計測幅は（ｂ−ａ）である。
【００２０】
一実施形態において、上述のｆＬ（ｔ）関数及びｆＲ（ｔ）関数は、Error function、Step function、Heaviside step function、Fermi-Dirac distribution function、Maxwell-Boltzmann distribution、Boltzmann distribution、Sigmoid function、Smoothstep function、Generalized logistic function、Richards’ curve、Logistic function、Logistic curve、Cumulative distribution function、percent point function、Quantile function、Survival function、またはReliability functionである。
【００２１】
一実施形態において、分光後の光パルスが測定位置で分散される場合、フィット曲線は、
−ｇ（ｔ）を更に含み、ここでは、−ｇ（ｔ）はプラズマ音声信号図中の分散現象が発生する図形に対応する。
【００２２】
一実施形態において、上述のｇ（ｔ）はガウス関数である。
【００２３】
一実施形態において、上述の光モジュールは、
第一分光ユニットと、
第一反射ユニットと
第二反射ユニットであって、第一分光ユニットにより光パルスが第一反射ユニット及び第二反射ユニットに分光されることと、
第一反射ユニット及び第二反射ユニットにより反射される光パルスを併合させる第二分光ユニットと、
第二分光ユニットからの光パルスを計測位置に集光させ、且つ計測位置でガスプラズマを発生させるフォーカスユニットとを備える。
【００２４】
一実施形態において、上述の第一分光ユニット及び第二分光ユニットはビームスプリッターであり、フォーカスユニットは軸外放物面鏡であり、収音ユニットはマイクロフォンである。
【００２５】
一実施形態において、上述の計測装置は、遮音空間を有し、少なくとも計測位置及び収音ユニットが遮音空間内に位置される遮音カバーを更に備える。
【発明の効果】
【００２６】
上述のように、本発明に係る計測装置及び計測方法は、ガスプラズマを発生させることにより光パルスの特性の分析を行う。この方法では非線形結晶や他の計測媒体が不要であるため、コストを節約でき、部材が破壊される危険性も回避される。
【図面の簡単な説明】
【００２７】
【図１】本発明の測定位置の実施形態を示す概略図である。
【図２Ａ】プラズマ音声信号の実施形態を示す概略図である。
【図２Ｂ】プラズマ音声信号の実施形態を示す概略図である。
【図３Ａ】測定チャートのフィット曲線の実施形態を示す概略図である。
【図３Ｂ】測定チャートのフィット曲線の実施形態を示す概略図である。
【図４】従来の非線形結晶の計測結果及び計測方法の計測結果の比較実施形態を示す概略図である。
【発明を実施するための形態】
【実施例】
【００２８】
本発明に係る計測方法及び計測装置はパルス光源に応用可能であり、パルス光源から発せられるパルス光の計測を行うことにより光パルスの光学的特性を取得させる。パルス光源は好ましくはレーザーパルス光源であり、例えば、フェムト秒（Femtosecond、１０-１５ sec）のレーザー、自己モード同期（Self mode-locking）により発生する短パルスレーザー、またはチャープパルス増幅（Chirped pulse amplification、CPA）により発生する高出力パルスレーザーである。
【００２９】
上述のパルスレーザーの出力は１０１３から１０１６Ｗ／ｃｍ２に達し、レーザーパルスのエネルギー密度が伝達媒体中の材料の破壊閾値（Threshold）より大きい場合、材料が融解して蒸発する。例えば、空気の破壊エネルギー閾値は約４×１０１３Ｗ／ｃｍ２であり、上述のレーザーパルスが集光された後にこのような出力によりガスプラズマを発生させる。
【００３０】
以下、本発明の実施形態に関わる計測装置を参照しながら本発明に関わる計測方法を説明する。本発明の実施形態に係る計測装置はパルス光源５０から発せられる光パルスＬの計測を行い、計測装置１００は光モジュール１１０と、収音ユニット１２０と、処理ユニット１３０とを備える（図１参照）。
【００３１】
光モジュール１１０により光パルスＬが受光されて計測位置Ｐに集光され、且つ計測位置Ｐで自己相関効果が発生すると共にガスプラズマが発生する。本発明は光モジュール１１０の詳細な部材構成に限られず、本分野で通常知識を有する技術者ならば適切な部材で置換させることで光モジュール１１０に同じ効果を発揮させることができ、ガスプラズマを発生させるガスの種類も限定されない。以下では空気を例に説明する。以下に例を挙げて本実施形態に係る計測装置１００について具体的に説明する。本実施形態に係る光モジュール１１０は、第一分光ユニット１１１と、第一反射ユニット１１２と、第二反射ユニット１１３とを備え、第一分光ユニット１１１により光パルスＬがＬ１及びＬ２に分光されると共に第一反射ユニット１１２及び第二反射ユニット１１３にそれぞれ伝送される。第一分光ユニット１１１は例えばビームスプリッター（Beam splitter）であり、ビームの一部分が透過させると共に一部分が反射されることにより光パルスＬが異なる２つの方向に伝達される光パルスＬ１、Ｌ２となる。
【００３２】
光モジュール１１０は第二分光ユニット１１４を更に備え、且つ第一反射ユニット１１２及び第二反射ユニット１１３は例えばリトロリフレクター（Retroreflector）である。第一反射ユニット１１２により２回反射されることにより光パルスＬ１が第二分光ユニット１１４に伝達され、第二反射ユニット１１３により２回反射されることにより光パルスＬ２も第二分光ユニット１１４に伝達される。光パルスＬ１及び光パルスＬ２は第一反射ユニット１１２及び第二反射ユニット１１３により各々反射された後に第二分光ユニット１１４に伝達され、第二分光ユニット１１４は例えば第一分光ユニット１１１に類似するビームスプリッターであり、その方向は互いに反対であるが、第二分光ユニット１１４は光パルスＬ１及び光パルスＬ２を透過及び反射させて同じ光路に伝達させ、同時に分散の補償を完成させる。
【００３３】
光モジュール１１０のフォーカスユニット１１５は第二分光ユニット１１４からの光パルスＬ１、Ｌ２を受光させると共にこれら光パルスＬ１、Ｌ２を計測位置Ｐに集光させるために用いられる。フォーカスユニット１１５は例えば軸外放物面鏡であり、好ましくは金めっき軸外放物面鏡であり、凹面によりこれら光パルスＬ１、Ｌ２を反射させて集光させる。
【００３４】
上述したように、本実施形態は空気を例にし、パルスレーザーが集光された後、パルスエネルギー密度が空気の破壊エネルギー閾値を超えれば空気がガスプラズマを発生させる。大出力パルスレーザーと材料とが相互作用する際に材料の表面に超音速パルスが発生する。前記音波を発生させるメカニズムは、電気歪（Electrostriction）と、熱膨張（Thermal expansion）と、蒸発（Vaporization）と、光学破壊（Optical breakdown）とを含む。レーザーパルスが材料を融解させると共に蒸発させて高温、高圧のプラズマを発生させる。発生するプラズマの温度及び圧力は非常に高く、プラズマが高速膨張し、プラズマが持続的に爆発して強烈に環境中の空気を圧縮させ、衝撃波を発生させて超音速で伝播させた後、高速で減衰して一般的な音波となる。
【００３５】
本実施形態に係る計測装置１００の収音ユニット１２０により計測位置Ｐから発生するガスプラズマの音波が収音されると共にプラズマ音声信号Ｓに変換される。収音ユニット１２０は例えばマイクロフォンであり、好ましくは周波数が５ｋHｚの良好な音声の周波数応答を有するマイクロフォンである。例えば、本実施形態に係る収音ユニット１２０と計測位置Ｐとの間の距離ｄは約１５ミリメートルであり、収音ユニット１２０がガスプラズマから発生する音声を収音した際に収音ユニット１２０が収音可能な音量範囲に維持されるが、本発明はこれに限られない。他の実施形態では、本分野で通常知識を有する者ならばマイクロフォンの違いに基づいて適切な計測距離に調整させることでプラズマ音声信号を取得可能である。
【００３６】
本実施形態に係る計測装置１００は処理ユニット１３０を更に備える。処理ユニット１３０は例えばパソコンやノートパソコンの中央処理装置（Central Processing Unit、CPU）であるが、本発明はこれに限られない。処理ユニット１３０は収音ユニット１２０に接続されると共にガスプラズマの音声に対応するプラズマ音声信号Ｓを受信させる。
【００３７】
更に詳しくは、本発明に係る計測装置１００は別途信号処理素子を組み合わせて良好なプラズマ音声信号Ｓを提供してもよい。例えば、計測装置１００がボックスカー１２１（Boxcar）及びロックインアンプ１２２（Lock-in amplifier）を更に備える。収音ユニット１２０により収音された音声はボックスカー１２１を経て絶対値が取得された後、ロックインアンプ１２２に伝送されることにより上述のプラズマ音声信号Ｓが取得される。但し、本発明はこれに限られず、本分野で通常知識を有する技術者ならば他の同等の効果を有する部材により類似する標本化効果を達成させることが可能である。
【００３８】
本実施形態に係る処理ユニット１３０はプラズマ音声信号Ｓに基づいて光パルスＬの光学的特性の計算を行う。換言すれば、本実施形態に係る計測装置１００は別途非線形結晶を使用して光パルスの時間軸上における長さ等の光学的特性を計測する必要がなく、全体的な製造コストが更に減少する。また、計測装置１００は光パルス強度が高すぎるために非線形結晶が破壊されることによる材料の損耗及びこれにより計測不能となる事態を回避させる。
【００３９】
より詳しくは、本実施形態に係る計測装置１００は、計測装置１００を設置させるための遮音空間Ａを有する遮音カバー１４０を更に備える。これにより、外界の雑音が収音ユニット１２０に影響を与えることがなくなる。他の実施形態では、上述の遮音カバーは上述の収音ユニット１２０及び計測位置Ｓを収容させるための遮音容器でもよく、収音ユニット１２０がガスプラズマから発生する音声を明確に収音可能になる。
【００４０】
以下、更に詳細な例を挙げて本発明に係る計測装置及びその計測方法について説明する。しかしながら、本発明に係る計測装置は以下の計測方式に制限されず、本分野で通常知識を有する者ならば需要に応じて上述のプラズマ音声信号Ｓを利用して必要な光パルスに関する各種特性を計測することが可能である。
【００４１】
図１に示されるように、本実施形態に係るプラズマ音声信号Ｓは測定チャート２０１を含み、測定チャート２０１はガスプラズマの時間軸に沿った音声強度に対応する。上述の収音ユニット１２０により収音された音声により、計測装置１００が取得させたプラズマ音声信号Ｓが含まれる測定チャート２０１は光パルスＬの光学的特性の計測に用いられ、測定チャート２０１はプラズマ音声信号中のガスプラズマの音声に対応する信号である。例えば、測定チャート２０１の時間軸上における幅Ｗを計測することにより、処理ユニット１３０は光パルスＬの時間的長さを取得させ、これにより光パルスＬの計測が完成する。
【００４２】
本実施形態に係るプラズマ音信号Ｓは測定チャート２０１にフィットさせることによりパルス光Ｌの光学的特性が計算される。図１に示されるように、処理ユニット１２０が測定チャート２０１に対してフィット曲線Ｇ１でフィットさせ、フィット曲線Ｇ１は例えばガウス関数（Gauss function）であり、処理ユニット１２０はガウス関数に基づいて計測幅Ｗの計算を行う。この計測幅Ｗは例えばフィット曲線Ｇ１の半値幅であり、即ち、測定チャート２０１の両側であり、各々の最小値から最大値までの平均点の時間軸上における幅Ｗである。この計測幅Ｗを計測時間値に換算させ、且つ測定チャート１２０がパルス光Ｌが発生させる空気パルスの時間的長さに対応することにより、処理ユニット１２０がこの計測時間値に基づいてパルス光Ｌの時間的長さに換算させる。
【００４３】
つまり、上述のフィット曲線Ｇ１により、本実施形態に係る計測装置１００がパルス光Ｌの光学的特性の計測を有効的に行い、同時に例えば非線形結晶の材料の消耗も回避させる。
【００４４】
しかしながら、本発明に係る計測装置及びその計測方法は上述のフィット曲線に限定されず、以下では更に例を挙げて他の実施形態における計測方法がどのようにフィット曲線を利用して光パルスの光学的特性の計測を行うかについて説明する。
【００４５】
図２Ａ及び図２Ｂを参照し、図２Ａは本発明の他の実施形態におけるプラズマ音声信号を示す概略図であり、図２Ｂは本発明の他の実施形態におけるプラズマ音声信号を示す概略図である。横軸は時間を示し、縦軸は音声強度に対応する電圧値を示す。以下の実施形態の説明を明確にするため、以下の実施形態で参照する測定チャートはフィット曲線で直接示し、図３Ａは図２Ａに対応するフィット曲線図であり、図３Ｂは図２Ｂに対応するフィット曲線図であるが、本発明に係る計測方法及び計測装置を限定させるものではない。
【００４６】
図３Ａを参照し、本発明の他の実施形態におけるプラズマ音声信号を示す測定チャートであり、フィット曲線Ｇ２によりフィットさせ、且つフィット曲線Ｇ２は例えば
【数７】
2019203874.tif 000009
であり、
ここでは、αは前記測定チャートに対応させる振幅の定数であり、
ここでは、
【数８】
2019203874.tif 000010
であり、ｋ１は定数であり、
ここでは、
【数９】
2019203874.tif 000011
であり、ｋ２は定数である。
また、測定チャート中の前記ガスプラズマの音声に対応する信号は（ａ−ｋ１）から（ｂ＋ｋ２）の間の範囲に位置される。具体的には、プラズマ音声信号の測定チャートは（ａ−ｋ１）の位置から上昇を始め、測定チャートの上昇部分はｆＬ（ｔ）によりフィットさせる。また、プラズマ音声信号の測定チャートは（ｂ＋ｋ２）の位置で終結し、終結前の測定チャートの下降部分はｆＲ（ｔ）によりフィットさせる。αは測定チャートの最大振幅に対応することにより上述の
【数１０】
2019203874.tif 000012
を測定チャートに完全フィットさせるために用いられ、これにより光パルスの光学的特性をフィット曲線Ｇ２により計測可能になる。
【００４７】
例えば、フィット曲線Ｇ２の計測幅は（ｂ−ａ）であり、即ち、フィット曲線中のｆＬ（ｔ）が半分まで上昇した時間点とｆＲ（ｔ）が半分まで下降した時間点の間の距離である。この計測幅（ｂ−ａ）によりこの測定チャートに対応する光パルスの時間的長さが計算される。
【００４８】
一方、ｆＬ（ｔ）は例えば誤差関数（error function）であり、即ち、ガウス関数の積分関数であり、ｆＲ（ｔ）は他の反対の誤差関数である。本実施形態に係るフィット曲線Ｇ２は誤差関数により測定チャートの図形にフィットさせ、これによりプラズマ音声信号が計測される。しかしながら、本発明は上述の関数の種類に制限されず、上述の測定チャートの形にフィットさせる関数として、Step function、Heaviside step function、Fermi-Dirac distribution function、Maxwell-Boltzmann distribution、Boltzmann distribution、Sigmoid function、Smoothstep function、Generalized logistic function、Richards’ curve、Logistic function、Logistic curve、Cumulative distribution function、percent point function、Quantile function、Survival function、またはReliability functionを含む。
【００４９】
図３Ｂは本発明の他の実施形態におけるプラズマ音声信号を示す測定チャートであり、他のフィット曲線により更にフィットさせ、且つこのフィット曲線は光パルスが分散される場合の測定チャートに応用可能である。ここで述べる分散とは、例えば、光パルスが伝達媒体が発生させるパルス幅により変化することであり、例えば、
【数１１】
2019203874.tif 000013
で示す。ここでは、ｔｏｕｔは媒体を経た後のパルス幅であり、ｔｉｎは媒体に進入する前のパルス幅であり、gは伝達媒体の厚さを考慮した群遅延分散（Group Delay Dispersion）であり、βは修正に用られる比率係数（Coefficient of proportion、後述する）である。図３Ｂに示されるフィット曲線は例えばgが2454fsの場合に光パルスが発生させるガスプラズマの測定チャートのフィット曲線である。
【００５０】
図３Ｂを参照し、光パルスが分散される場合、上述の計測装置１００では、計測位置Ｓに発生するガスプラズマが発生させるプラズマ音声信号の測定チャート中には陥没する図形が生成される。このため、本実施形態ではフィット曲線Ｇ３を利用して上述の測定チャートにフィットさせる。フィット曲線Ｇ３は、
【数１２】
2019203874.tif 000014
を含み、ここでは、−ｇ（ｔ）は前記プラズマ音声信号図中の分散現象により発生する図形に対応する。即ち、上述の陥没部分に対応する。これにより、ガスプラズマのプラズマ音声信号が発生させる測定チャートの形に完全にフィットさせる。ｇ（ｔ）は例えばガウス関数であり、陥没する図形に対応するが、但し、本発明は関数の種類が制限されない。
【００５１】
また、上述のフィット曲線は測定チャートの両側に対してｆＬ（ｔ）、ｆＲ（ｔ）を使用してフィットさせる。このため、同様に２つの関数の間の時間差により幅が計測され、光パルスのパルスの時間的長さが更に算出される。特に説明する必要なのは、ここで述べる関数ｆＬ（ｔ）、ｆＲ（ｔ）は上述の誤差関数等と同じ種類のフィット関数を指し、異なる測定チャートにフィットさせる場合、異なる数値を使用して修正を行い、上述の図３Ａの内容で用いられるｆＬ（ｔ）、ｆＲ（ｔ）と完全に同じになるわけではない。
【００５２】
図４は従来の非線形結晶を利用して計測されたパルス幅Ｓ１（図中の円データ点に対応する）及び本発明の実施形態に係る計測装置に用いられる計測方法を利用して計測されたパルス幅Ｓ２（図中の三角データ点に対応する）の異なる分散状況（即ち、分散補正係数、上述の関連分散係数gに対応する）における数値を示し、第一幅の縦座標は本発明の実施形態に係る計測装置に用いられる計測方法によって計測された数値に対応し、第二幅の縦座標は従来の非線形結晶によって計測されたパルス幅の数値に対応する。図３から分かるように、本発明の実施形態に係る計測装置により計測されたパルス幅は全て従来の非線形結晶により計測されたパルス幅に線形対応し、よって本発明に係る計測装置及びこれを使用する計測方法は光パルスの幅を有効的に計測可能であり、更に異なる分散状況に応じて修正を加えることも可能である。換言すれば、図３の計測結果によって、各分散状況において対応する分散補正値を算出可能である。例えば、
【数１３】
2019203874.tif 000015
の比率係数βを１．８４に修正し、パルス幅Ｓ２をパルス幅Ｓ１に修正する。計測装置が上述の例えば図２Ｂに示されるフィット曲線により測定チャートの形にフィットさせると共に幅を計測させた後、分散補正値及び比率係数の修正により計測された幅の値の補正を行うことでより正確になる。
【００５３】
総合すると、本発明に係る計測装置は光パルスが空気中に発生させるガスプラズマにより光パルスの特性を直接計測させる。このため、計測装置に高価な非線形結晶が不要になり、全体的な製造コストが低下する。また、本発明に係る計測方法により、非線形結晶がなくても光パルスの光学的特性を計測可能になる。
【符号の説明】
【００５４】
ｄ  距離
Ｌ光パルス
Ｌ１  光パルス
Ｌ２  光パルス
Ｐ  計測位置
Ｓ  プラズマ音声信号
Ｗ  幅
５０  パルス光源
１００  計測装置
１１０  光モジュール
１１１  第一分光ユニット
１１２  第一反射ユニット
１１３  第二反射ユニット
１１４  第二分光ユニット
１１５  フォーカスユニット
１２０  収音ユニット
１２１  ボックスカー
１２２  ロックインアンプ
１３０  処理ユニット
１４０  遮音カバー
２０１  測定チャート
【図１】
2019203874.tif 000022
【図２Ａ】
2019203874.tif 000023
【図２Ｂ】
2019203874.tif 000024
【図３Ａ】
2019203874.tif 000025
【図３Ｂ】
2019203874.tif 000026
【図４】
2019203874.tif 000027