(57)【要約】  周期的な繰り返し現象を観測するフーリエ変換型分光におけるスペクトル分解能とスペクトル確度の制限を解消し、理論上、無限小のスペクトル分解（無限大のスペクトル分解能）を実現し、スペクトル確度を向上できるフーリエ変換型分光法を提供する。現象の繰り返し周期を正確かつ十分に安定化した上で、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と繰り返し周期とを厳密に一致させて時間波形を取得し、フーリエ変換することにより、観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する。繰り返し周期を変化させながら計測を繰り返し、離散分布スペクトルの間隙を補完する。これにより現象の存在時間が繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度を向上できる。
2014034085.tif 000002
【特許請求の範囲】
【請求項１】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式１が成立し、
現象の存在時間が前記繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度が向上し得る、ことを特徴とするフーリエ変換型分光法。
（数１）
ｆｎ＝ｎ／Ｔ      ・・・（式１）
【請求項２】
前記観測時間窓Ｔにおいて、前記現象の繰り返し周期と一致させる替わりに、前記繰り返し周期よりも短い時間で現象を観測して時間波形を取得し、取得した時間波形データに対して繰り返し周期に一致させるように時間波形データを補間し、繰り返し周期と一致した時間波形を取得する、ことを特徴とする請求項１に記載のフーリエ変換型分光法。
【請求項３】
前記繰り返し周期を変化させるステップと、
繰り返し周期を変化させた後の時間波形を取得し、フーリエ変換することにより、離散分布スペクトルを取得するステップと、
異なる繰り返し周期の複数の離散分布スペクトルを重畳するステップと、
を更に備え、
上記の取得した離散分布スペクトルの各プロットの間の間隙を補完することを特徴とする請求項１又は２に記載のフーリエ変換型分光法。
【請求項４】
前記繰り返し周期を、周波数標準器を基準として安定化させた、請求項１〜３のいずれかに記載のフーリエ変換分光法。
【請求項５】
前記繰り返し周期の整数分の１の時間間隔で時間波形データを離散化して取得した、或いは、前記繰り返し現象の周期を時間波形データの離散化時間間隔の整数倍に設定して取得した、前記繰り返し周期１周期分の離散化した時間波形データを離散フーリエ変換することにより離散分布スペクトルを得る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換分光法。
【請求項６】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）に用いられ、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、
ことを特徴とする請求項１〜５のいずれかに記載のフーリエ変換型分光法。
【請求項７】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）に用いられ、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、
ことを特徴とする請求項１〜５のいずれかに記載のフーリエ変換型分光法。
【請求項８】
請求項１〜５のいずれかに記載のフーリエ変換型分光法を用いた、
核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析（ＦＴ−ＭＳ）装置，フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーのいずれかのフーリエ変換型分光装置。
【請求項９】
周期的な繰り返し現象を観測するフーリエ変換型分光装置において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段と、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段と、
を備え、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式２が成立し、
現象の存在時間が前記繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度が向上し得る、
ことを特徴とするフーリエ変換型分光装置。
（数２）
ｆｎ＝ｎ／Ｔ      ・・・（式２）
【請求項１０】
前記繰り返し周期を変化させる手段と、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段と、
を更に備え、
上記の取得した離散分布スペクトルの各プロットの間の間隙を補完することを特徴とする請求項９に記載のフーリエ変換型分光装置。
【請求項１１】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）を行うテラヘルツ時間領域分光装置において、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザー手段の各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザー手段をそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、
ことを特徴とする請求項９又は１０に記載のフーリエ変換型分光装置。
【請求項１２】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）を行うフーリエ変換型赤外分光装置において、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、
ことを特徴とする請求項９又は１０に記載のフーリエ変換型分光装置。
【請求項１３】
周期的な繰り返し現象を観測するフーリエ変換型分光計測プログラムであって、
コンピュータを、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段、
として機能させるプログラムであって、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式３が成立し、
現象の存在時間が前記繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度が向上し得る、
ことを特徴とするフーリエ変換型分光計測プログラム。
（数３）
ｆｎ＝ｎ／Ｔ      ・・・（式３）
【請求項１４】
コンピュータを、
前記繰り返し周期を変化させる手段、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段、
を更に備え、
上記の取得した離散分布スペクトルの各プロットの間の間隙を補完することを特徴とする請求項１３に記載のフーリエ変換型分光計測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スペクトル分解能とスペクトル確度を向上するフーリエ変換型分光法、分光装置および分光計測プログラムに関するものである。
【背景技術】
【０００２】
フーリエ変換型分光法（フーリエ変換型周波数解析法）とは、測定信号強度の時間変化を測定してからフーリエ変換することで周波数ごとのスペクトルを得る分光法である。
ある周期信号を考えた場合、時間領域で観測された時間波形信号ｈ(ｔ)と周波数領域で観測されたスペクトル信号Ｈ(ｆ)は、下記数式に示すようにフーリエ変換（数式（１））および逆フーリエ変換（数式（２））によって結ばれている。【０００３】
【数１】
2014034085.tif 000003
【０００４】
【数２】
2014034085.tif 000004
【０００５】
フーリエ変換と逆フーリエ変換で結ばれることで、時間領域と周波数領域で計測した信号は、それぞれ等価であると言える。コンピュータの登場以前は、フーリエ変換を行うことが容易ではなかったため、スペクトル信号を計測する場合は、一般に周波数領域の計測が行われていた。
例えば、回折分光器に代表される分散型分光計による光波の分光計測では、光波を分散素子（回折格子、プリズムなど）で波長毎に空間的に分離し、スリットによって特定の波長成分のみを抽出することにより、スペクトルを取得していた。
しかし、コンピュータの登場により、フーリエ変換を行うことが極めて容易になったため、フーリエ変換型分光法、すなわち、測定信号強度の時間変化を測定してからフーリエ変換することで周波数ごとのスペクトルを取得することが可能になった。【０００６】
フーリエ変換型分光法の代表的な物として、パルス・フーリエ変換核磁気共鳴分光（ＦＴ−ＮＭＲ）、核磁気共鳴画像法（ＭＲＩ）、フーリエ変換赤外分光法（ＦＴ−ＩＲ）、テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）などがあり、産業分野や医療分野で幅広く使われている。ＦＴ−ＮＭＲでは、パルス状高周波電波（ラジオ波）を、超伝導磁石の中央にセットされたサンプル管に照射することにより、ＦＩＤ（Free Induction Decay；自由誘導減衰）信号という時間領域信号を測定し、これをフーリエ変換することによりＮＭＲスペクトルを得る。核磁気共鳴画像法（ＭＲＩ）では、ＮＭＲスペクトルをコンピュータ断層撮影に応用する。ＦＴ−ＩＲやＴＨｚ−ＴＤＳでは、電磁波の干渉波形や電場波形を時間の関数として観測し、それをフーリエ変換することによってスペクトルを得る。【０００７】
フーリエ変換型分光法の特徴は、下記１）〜４）に示すとおりである。
１）全スペクトル領域を同時に観測するため、信号強度が強い。
フーリエ変換型分光法における干渉信号の測定時には、入力全スペクトル領域を同時に観測するため、信号強度が強く、結果として高い信号雑音比が得られる。
２）スリットが不要なので、分散型分光器よりも明るい光学系が使用可能である。
分散型分光計と異なりスリットが不要で、分散型分光計よりも明るい光学系が使用できるため、高い光学的スループットが実現でき高い信号雑音比が得られる。
３）連続スペクトルを計測でき、スペクトル精度が高い。
連続した波長の光を観測でき、波数精度が非常に高い。
４）様々な電磁波領域で利用可能である。【０００８】
上記のフーリエ変換型分光法の特徴から、フーリエ変換型分光法は分光計測の主流になりつつあり、様々な領域で幅広く使われている。【０００９】
【表１】
2014034085.tif 000005
【００１０】
分光計測において重要な性能は、スペクトル分解能とスペクトル確度である。
図１は、観測信号の時間波形と、フーリエ変換によって取得される振幅スペクトルを示している。フーリエ変換型分光法におけるスペクトル分解能は、時間波形を観測する時間幅（観測時間窓）の逆数によって決まる。
したがって、観測時間窓を大きくすればするほどスペクトル分解能は向上することになるが、周期的に繰り返す信号を観測する場合には、設定可能な観測時間窓は繰り返し周期になるため、得られるスペクトル分解能の上限は繰り返し周波数（理論限界スペクトル分解能）で制限される。しかし、測定信号の大部分は時間的に局在している場合、観測時間窓の拡大はノイズ信号成分を増大させ、測定ＳＮ比を低下させる。【００１１】
また、測定時間も、観測時間窓の拡大に比例して増大する。特に、ＴＨｚ−ＴＤＳでは、時間波形取得のために機械式時間遅延走査を行う必要があるが、時間遅延走査を行う機械式ステージの移動量によって観測時間窓が制限され、現実的なスペクトル分解能は繰り返し周波数で決定されるスペクトル分解能を大きく下回っているのが実情であり、十分に大きな観測時間窓を実現するのは容易ではない。このような点を考慮して実際の観測時間窓が設定され、スペクトル分解能が決まる。
一方、スペクトル確度は、時間波形観測の時間精度に依存する。
分光分析における物質識別能力を今まで以上に向上させるためには、更なるスペクトル分解能とスペクトル確度の向上が必要である。【００１２】
既に、本発明者の一部の発明者は、ＴＨｚ−ＴＤＳにおいて、ＴＨｚスペクトルを理論限界スペクトル分解（＝繰り返し周波数＝モード同期周波数）で高速に分光測定できるような計測装置を提案している（特許文献１を参照）。
本発明は、ＴＨｚ−ＴＤＳを含めたフーリエ変換型分光におけるスペクトル分解能の制限を解消し、理論上、無限小のスペクトル分解（スペクトル分解能が無限大）を実現すると共に、スペクトル確度を大幅に向上できるものである。
【先行技術文献】
【特許文献】
【００１３】
【特許文献１】特許第４５６５１９８号公報
【発明の概要】
【発明が解決しようとする課題】
【００１４】
上述の如く、分光分析における物質識別能力を今まで以上に向上させるためには、更なるスペクトル分解能とスペクトル確度の向上が必要である。
上記状況に鑑みて、本発明では、フーリエ変換型分光におけるスペクトル分解能の制限を解消し、理論上、無限小のスペクトル分解を実現すると共に、スペクトル確度を大幅に向上できるフーリエ変換型分光法、分光装置および分光計測プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１５】
上記目的を達成すべく、本発明のフーリエ変換型分光法は、周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する。
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式３が成立する。【００１６】
【数３】
2014034085.tif 000006
【００１７】
これにより、現象の存在時間が繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度を向上できる。【００１８】
フーリエ変換型分光において、現象の繰り返し周期を正確かつ十分に安定化した上で、観測時間窓Ｔを繰り返し周期と厳密に一致させることにより、フーリエ変換型分光法におけるスペクトル分解能の上限（現象の繰り返し周波数で決定される）を解消し、理論上、無限小のスペクトル分解を実現する。
上記の観測時間窓Ｔにおいて、繰り返し周期と一致させる替わりに、繰り返し周期よりも短い時間で現象を観測して時間波形を取得し、取得した時間波形データに対して繰り返し周期に一致させるように時間波形データを補間して（例えば、０データの補間）、繰り返し周期と一致した時間波形を取得することでも構わない。ただし、繰り返し周期と補間後のデータは、厳密に一致させる必要がある。この場合、繰り返し周期と観測時間窓Ｔを一致させて時間波形を取得する場合より、データ補間部分の時間波形の欠損があるためフーリエ変換スペクトルの確度が低下するデメリットが生ずる。【００１９】
上記の本発明のフーリエ変換型分光法において、繰り返し周期を変化させるステップと、繰り返し周期を変化させた後の時間波形を取得し、フーリエ変換することにより、離散分布スペクトルを取得するステップと、異なる繰り返し周期の複数の離散分布スペクトルを重畳するステップと、を更に備えることが好ましい。
繰り返し周波数を変化させて取得した離散分布スペクトルの各プロットの間の間隙を補完することにより、微細化スペクトルを取得する。また、繰り返し周期を正確かつ十分に安定化することにより、スペクトル分解能のみならず、スペクトル確度も大幅に向上できる。【００２０】
ここで、繰り返し周期を、周波数標準器を基準として安定化させたことが好ましい。周波数標準器を基準とすることで、繰り返し周期を正確かつ十分に安定化できる。
また、繰り返し周期の整数分の１の時間間隔で時間波形データを離散化して取得した、或いは、前記繰り返し現象の周期を時間波形データの離散化時間間隔の整数倍に設定して取得した、前記繰り返し周期１周期分の離散化した時間波形データを離散フーリエ変換することにより離散分布スペクトルを得ることがより好ましい。フーリエ変換スペクトルの確度をより高めることができる。【００２１】
上記の本発明のフーリエ変換型分光法は、テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）に好適に用いられる。その際、テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用いる。２台のフェムト秒レーザーの各々のモード同期周波数は周波数標準器を基準として正確かつ十分に安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御する。そして、一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用いる。
これにより、繰り返し周期であるパルス周期の安定化を図ることができる。【００２２】
また、上記の本発明のフーリエ変換型分光法は、フーリエ変換型赤外分光法（ＦＴ−ＩＲ）を行うフーリエ変換型赤外分光装置に好適に用いられる。その際、フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を用いる。２台のフェムト秒レーザーの各々のモード同期周波数は周波数標準器を基準として正確かつ十分に安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御する。そして、一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用いる。
これにより、繰り返し周期の安定化を図ることができる。【００２３】
下記文献に開示されている配置では、赤外インコヒーレント光源（セラミック光源、タングステンハロゲン光源、タングステン・ヨウ素ランプなど）を用いたホモダイン干渉でなく、２台のフェムト秒レーザーと非線形光学結晶で発生させた赤外コヒーレント光によるヘテロダイン干渉を用いており、上記の本発明のフーリエ変換型分光法に用いることができる。【００２４】
（文献）Ｆ．  Ｋｅｉｌｍａｎｎ，  Ｃ．  Ｇｏｈｌｅ，  ａｎｄ  Ｒ．  Ｈｏｌｚｗａｒｔｈ，  “Ｔｉｍｅ-ｄｏｍａｉｎ  ｍｉｄ-ｉｎｆｒａｒｅｄ  ｆｒｅｑｕｅｎｃｙ-ｃｏｍｂ  ｓｐｅｃｔｒｏｍｅｔｅｒ，”  Ｏｐｔｉｃｓ  Ｌｅｔｔｅｒｓ，  Ｖｏｌ．  ２９  Ｉｓｓｕｅ  １３，  ｐｐ．１５４２-１５４４  （２００４）．【００２５】
上記の本発明のフーリエ変換型分光法は、核磁気共鳴分光法（ＮＭＲ），核磁気共鳴画像法（ＭＲＩ），フーリエ変換質量分析法（ＦＴ−ＭＳ），フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーのいずれかのフーリエ変換型分光法に好適に用いることができる。核磁気共鳴分光法（ＮＭＲ），核磁気共鳴画像法（ＭＲＩ），フーリエ変換質量分析法（ＦＴ−ＭＳ），フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーは、繰り返し現象を観測して、フーリエ変換型分光法を用いて計測するからである。【００２６】
次に、本発明のフーリエ変換型分光装置について説明する。
本発明のフーリエ変換型分光装置は、周期的な繰り返し現象を観測するフーリエ変換型分光装置において、現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段と、取得した時間波形からフーリエ変換することにより、観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段と、を備える。
かかる構成によって、現象の存在時間が繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度を向上できる。
なお、上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、上述の数式３が成立する。【００２７】
ここで、上記の本発明のフーリエ変換型分光装置において、繰り返し周期を変化させる手段と、異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段と、を更に備えることが好ましい。
繰り返し周期を変化させて取得した離散分布スペクトルの各プロットの間の間隙を補完し、微細化スペクトルを取得することにより、スペクトル分解能が大幅に向上できる。また、繰り返し周期を正確かつ安定化することにより、スペクトル確度も大幅に向上できる。【００２８】
また、上記の本発明のフーリエ変換型分光装置は、テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）を行うテラヘルツ時間領域分光装置として好適に用いることができる。
この場合、テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備える。
２台のフェムト秒レーザー手段の各々のモード同期周波数は、周波数標準器を基準として正確かつ十分に安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザー手段をそれぞれ独立して同時に制御する。
そして、一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用いる。
これにより、繰り返し周期であるパルス周期の安定化を図ることができる。【００２９】
また、上記の本発明のフーリエ変換型分光装置は、フーリエ変換型赤外分光法（ＦＴ−ＩＲ）を行うフーリエ変換型赤外分光装置として好適に用いることができる。
この場合、フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備える。
２台のフェムト秒レーザーの各々のモード同期周波数は、周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御する。
そして、一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用いる。
これにより、繰り返し周期の安定化を図ることができる。【００３０】
次に、本発明のフーリエ変換型分光計測プログラムについて説明する。
本発明のフーリエ変換型分光計測プログラムは、周期的な繰り返し現象を観測するフーリエ変換型分光計測プログラムであって、コンピュータを、現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段、取得した時間波形からフーリエ変換することにより、観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段、として機能させるためのプログラムである。
このプログラムによって、現象の存在時間が繰り返し周期よりも長い観測対象の場合に、取得した離散分布スペクトルのスペクトル分解が無限小（スペクトル分解能が無限大）となり、スペクトル確度を向上できる。
なお、上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、上述の数式３が成立する。【００３１】
ここで、上記の本発明のフーリエ変換型分光計測プログラムにおいて、繰り返し周期を変化させる手段と、異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段と、を更に備え、コンピュータをこれらの手段として機能させることが好ましい。
繰り返し周期を変化させて取得した離散分布スペクトルの各プロットの間の間隙を補完し、微細化スペクトルを取得することにより、スペクトル分解能のみならず、スペクトル確度も大幅に向上できる。【００３２】
上記の本発明のフーリエ変換型分光計測プログラムは、核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析（ＦＴ−ＭＳ）装置，フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーなどのフーリエ変換型分光装置に好適に搭載することができる。
いずれの装置の場合も、繰り返し現象を観測するものであり、スペクトル分解能とスペクトル確度を向上することが可能である。
【発明の効果】
【００３３】
本発明によれば、フーリエ変換型分光におけるスペクトル分解能とスペクトル確度の制限を解消し、理論上、無限小のスペクトル分解を実現すると共に、スペクトル確度を大幅に向上できるといった効果を有する。
本発明を利用すると、核磁気共鳴分光法（ＮＭＲ），フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーなどのフーリエ変換型分光法のハード面を大幅に改造することなく、性能を大幅に向上できるようになる。フーリエ変換型分光法におけるスペクトル分解能とスペクトル確度の向上は、各種物質同定の識別能力を高め、各種機器の高性能化につながる。その結果、半導体や分析化学の分野における利用を更に促進する。また、ＦＴ−ＮＭＲをコンピュータ断層撮影法に応用した核磁気共鳴画像法（ＭＲＩ）においても、生体断層画像の高精細化に繋がる可能性がある。
【図面の簡単な説明】
【００３４】
【図１】フーリエ変換型分光法の説明図
【図２】測定時間波形のフーリエ変換の説明図
【図３】繰り返し周期よりも長い緩和信号を繰り返し現象として観測するイメージ図
【図４】異なるタイミングで誘起され時間的に異なった緩和現象を１周期の観測時間窓で観測するイメージ図
【図５】１周期で観測された繰り返し現象のフーリエ変換の説明図（１）
【図６】１周期で観測された繰り返し現象のフーリエ変換の説明図（２）
【図７】スペクトルを構成するプロット分布の様子（プロット間隔１／Ｔ）。(ａ)無限小スペクトル幅で切り出されたプロットの離散分布（本発明の場合），(ｂ)平均値で与えられたプロットの連続分布（従来方法の場合）
【図８】観測時間窓が１周期と等しくない場合の説明図
【図９】離散分布スペクトル間隙の補完の説明図
【図１０】従来のＴＨｚ−ＴＤＳ装置（機械式時間遅延走査ポンプ・プローブ法）の説明図
【図１１】実施例１のＴＨｚ−ＴＤＳの原理の説明図
【図１２】実施例１のＴＨｚ−ＴＤＳ装置の構成図
【図１３】水蒸気の吸収スペクトル（回転遷移１１０−＞１０１）。(ａ)本発明１（離散分布スペクトル間隙の補完無し），(ｂ)本発明２（離散分布スペクトル間隙の補完有り），(ｃ)従来方法（ゼロ・フィリング法）
【図１４】吸収緩和現象が繰り返し周期より長い場合の吸収スペクトル。(ａ) 観測窓＝１周期，(ｂ) 観測窓＝０．９９９９５周期
【図１５】吸収緩和現象が繰り返し周期より短い場合の吸収スペクトル。(ａ) 観測窓＝１周期，(ｂ) 観測窓＝０．９９９９５周期
【図１６】本発明（離散分布スペクトル間隙の補完無し），本発明（離散分布スペクトル間隙の補完有り）および従来方法（ゼロ・フィリング法）における圧力拡がり線幅と吸収スペクトルの相関グラフ
【図１７】フーリエ変換型分光計測プログラムの機能ブロック図
【図１８】フーリエ変換型分光計測プログラムの処理フロー図
【図１９】コンピュータ・ハードウェアの内部構成図
【図２０】アセトニトリルガスの微細構造を示す振幅スペクトル
【図２１】０．６〜０．７ＴＨｚの範囲の吸収スペクトル
【図２２】Ｊ＝３４〜３５遷移付近（０．６４ＴＨｚ付近）を拡大した吸収スペクトル
【図２３】モード同期周波数を変化させ得られた吸収スペクトルであって、０．６４３２ＴＨｚ付近を拡大した吸収スペクトル
【発明を実施するための最良の形態】
【００３５】
以下、本発明の実施形態について、図面を参照しながら詳細に説明していく。なお、本発明の範囲は、以下の実施例や図示例に限定されるものではなく、幾多の変更及び変形が可能である。
実施例を説明する前に、本発明の原理について図を参照しながら詳細に説明する。【００３６】
（本発明の原理について）
時間波形信号ｈ(ｔ)を取得し、そのスペクトル波形Ｈ(ｆ)をフーリエ変換によって取得する場合を想定する。測定信号ｈ(ｔ)のある周波数成分ｆを得るためには、前述のフーリエ変換（数式（１））から、測定信号ｈ(ｔ)と観測したい周波数信号ｃｏｓ２πｆｔの掛け算を行い、時間領域で積分することによって算出される（図２を参照、実数部のみについて記述。）。
さらに、各周波数成分に対して同様な処理を行うことにより、スペクトル波形を取得できる。ここで、時間積分領域の大きさ（観測時間窓に相当）が、スペクトル分解能を決定することになる。
したがって、スペクトル分解能を向上するためには、時間積分領域（観測時間窓）を大きく取る必要があるが、通常は必要なスペクトル分解能、測定ＳＮ比、測定時間、観測時間窓の技術的制限などを考慮して、有限の観測時間窓で計測が行われる。【００３７】
次に、繰り返し周期よりも長い緩和信号（例えば、ＮＭＲにおける自由誘導減衰など）を繰り返し現象として観測する場合を考える。
図３は、繰り返し周期よりも長い緩和信号を繰り返し現象として観測するイメージ図を示している。図３の場合、異なるタイミングの４つのパルスによって励起され、時間的に異なった４つの緩和現象信号（直近の信号１０，１つ前の信号１１，２つ前の信号１２，３つ前の信号１３）が時間的に重なった状態となる。【００３８】
このような時間的に重なった状態の信号を、１周期の観測時間窓で観測すると、図４に示すように、直近のパルス信号１０による緩和現象の時間領域（Ａ）、１つ前のパルス信号１１による緩和現象の時間領域（Ｂ）、２つ前のパルス信号１２による緩和現象の時間領域（Ｃ）、３つ前のパルス信号１３による緩和現象の時間領域（Ｄ）、が観測時間窓に含まれることになる。
図５のように、１周期で観測された繰り返し現象のフーリエ変換スペクトルを求めるためには、時間領域（Ａ），（Ｂ），（Ｃ），（Ｄ）に含まれる各緩和信号の和と観測したい周波数信号ｃｏｓ２πｆｔの掛け算を時間領域で積分すればよい。【００３９】
ここで、緩和現象を引き起こす光源の繰り返し周期を、正確にかつ十分に安定化した上で、その繰り返し周期と観測時間窓が厳密に一致するようにして測定すると、図６に示すように、時間領域（Ａ），（Ｂ），（Ｃ），（Ｄ）に含まれる各信号を時間的に連続な信号として連結することが可能になる。
すなわち、繰り返し周期に等しい観測時間窓であるにも関わらず、それよりも十分に長い緩和現象の信号を無限長の観測時間窓（時間積分領域）で測定することと等価になり、スペクトル分解は理論上無限小となり、スペクトル分解能が無限大となるのである。【００４０】
ここで、繰り返し周期（観測時間窓）をＴとすると、無限小のスペクトル分解（無限大のスペクトル分解能）でスペクトル情報を得られる周波数は、図７（ａ）に示すように、１／Ｔ間隔で離散分布し、各プロットの周波数は上述の数式（３）で与えられる。ここで、ｆｎは各プロットの周波数、ｎはプロットの次数（整数）である。
また、図７は、スペクトルを構成するプロット分布の様子（プロット間隔１／Ｔ）を示している。図７(ａ)は、無限小スペクトル幅で切り出されたプロットの離散分布（本発明の場合）であり、一方、図７(ｂ)は、平均値で与えられたプロットの連続分布（従来方法の場合）である。【００４１】
次に、サンプリングを行うことで、無限小スペクトル幅で切り出されたプロットの離散分布が得られることについて説明する。
単一現象ｈ（ｔ）を、ｈ（ｔ）の最大周波数の２倍以上のサンプリング周波数で取得すると、サンプリング定理より、その情報を失わずに離散化できる。
ｈ（ｔ）の最大周波数の２倍以上のサンプリング周波数で取得したスペクトルは、フーリエ級数展開を用いて、下記数式（４）で表すことができる。ここで、ｍは、整数でありτはサンプリング間隔である。ここで、サンプリング定理を満たすために、１／τは、ｈ（ｔ）の最大周波数の２倍以上でなければならない。【００４２】
【数４】
2014034085.tif 000007
【００４３】
一方、単一現象ｈ（ｔ）を、サンプリング周期のＮ倍の時間間隔Ｔ＝Ｎτで繰り返し、重ね合わせると下記数式（５）となる。【００４４】
【数５】
2014034085.tif 000008
【００４５】
上記数式（５）におけるｇ（ｍτ）において、繰り返し１周期分、すなわち、  ｍ＝０，・・・，Ｎ−１のデータを用いた離散フーリエ変換は、下記数式（６）で定義される。【００４６】
【数６】
2014034085.tif 000009
【００４７】
上記数式（６）において、ｆｎは、前述の数式（３）であり、ｎは整数である。
上記数式（６）は、下記数式（７）に展開することができる。【００４８】
【数７】
2014034085.tif 000010
【００４９】
すなわち、現象の繰り返し周期ＴをＮ等分して信号をサンプリング（ただし、サンプリング定理を満たすため、Ｎ／Ｔを現象の最大周波数の２倍以上とすることが望ましい。）、或いは、現象の繰り返し周期を時間波形データのサンプリング間隔のＮ倍に設定し、繰り返し周期１周期分の離散化データを取得して、これを離散フーリエ変換して離散スペクトルを得ると、その離散スペクトルの各点は、元となった単一現象を無限小のスペクトル分解（無限大のスペクトル分解能）で観測した値に等しいことになる。なお、Ｎが２のべき乗の場合、高速離散フーリエ変換を使用できるため、演算時間を大幅に小さくすることが可能となる。【００５０】
図７において、注目すべき点は、スペクトル分解が無限小（スペクトル分解能が無限大）として得られるのは、離散的になっている点である。これは、従来のスペクトル（例えば、分散型分光法）を構成する各プロットの場合、図７（ｂ）に示すように、プロット間隔（１／Ｔ）の平均値として連続分布になっているのに対して、本発明の場合、図７（ａ）に示すように、各プロットのスペクトル分解（切り出しスペクトル幅）が無限小の離散分布になっている。【００５１】
一方、仮に観測時間窓が厳密に１周期に等しくなければ、時間領域（Ａ），（Ｂ），（Ｃ），（Ｄ）に含まれる各信号を時間的に連続な信号として連結することができず、時間的不連続点を含む信号に変化してしまう（図８を参照）。その結果、緩和現象の信号を無限長の観測時間窓（時間積分領域）で測定することと等価にならない。すなわち、上記数式（３）を満たされなければ、正しいスペクトルは得られないのである。【００５２】
図７（ａ）に示したスペクトルの各プロットは、無限小のスペクトル分解（無限大のスペクトル分解能）を有しているものの、離散的に分布している。そのため、このまま利用すると、実用上のスペクトル分解能はプロット間隔（パルス周期の逆数１／Ｔ）になる。
各プロットが有する無限小のスペクトル分解（無限大のスペクトル分解能）を有効利用するためには、各プロットの間に存在する間隙部分を補完する必要がある。そのためには、パルス周期Ｔを変化させながら、図７（ａ）のスペクトル波形を逐次取得し、それらを最終的に重ねて補完することにより、より狭いプロット間隔を有するスペクトルを得ることが可能になる（図９を参照）。
これは、上記数式（３）において、パルス周期Ｔを変化させながら、離散的なスペクトル情報を計測することと等価である。その結果、実用上プロット点数の増大分だけ、スペクトル分解能が大幅に向上する。【００５３】
また、緩和現象の繰り返し周期が周波数標準器を基準として予め正確かつ十分に安定化され、その繰り返し周期と観測時間窓（Ｔ）を厳密に一致させて測定しているので、各プロットの間隔は常に一定で、スペクトル波長（周波数）の絶対精度は極めて正確になる。その結果、スペクトル分解能の向上により波長（周波数）の読み取り精度が上がり、スペクトル確度も向上することになる。【００５４】
以下の実施例では、ＴＨｚ−ＴＤＳにおけるフーリエ変換型分光法と分光装置を例に挙げて、本発明について説明する。また、低圧状態の水蒸気をサンプルに用いて、本発明のフーリエ変換型分光法の有用性を評価する。
【実施例１】
【００５５】
実施例１として、ＴＨｚ−ＴＤＳにおける本発明のフーリエ変換型分光法の実施例を示す。ＴＨｚ−ＴＤＳでは、自由空間を伝搬するパルス状ＴＨｚ波（ＴＨｚパルス）の電場時間波形を取得し、図１のように、フーリエ変換によって振幅スペクトルを取得する。しかし、検出エレクトロニクスの帯域不足から、サブピコ秒〜ピコ秒オーダーのＴＨｚパルスの電場時間波形を直接計測（ＴＨｚパルスを実時間測定）することは不可能であるため、機械式時間遅延走査に基づいたポンプ・プローブ法と呼ばれる手法がこれまで用いられてきた。【００５６】
図１０は、典型的なＴＨｚ−ＴＤＳ装置（機械式時間遅延走査ポンプ・プローブ法を用いた装置）の構成を示している。フェムト秒レーザー３１光は、ビームスプリッター３２によって、ＴＨｚ発生用ポンプ光（ｐｕｍｐ  ｌｉｇｈｔ）及びＴＨｚ検出用プローブ光（ｐｒｏｂｅ  ｌｉｇｈｔ）に分割される。ポンプ光をＴＨｚ発生素子３４（光伝導アンテナなど）に入射することにより、ＴＨｚパルスが発生する。自由空間を伝搬したＴＨｚパルスは、ＴＨｚ検出素子３５（光伝導アンテナなど）に入射される。プローブ光は、時間遅延走査ステージ３３によってタイミング（光路長）調整された後、ＴＨｚ検出素子３５に入射される。【００５７】
ここで、図１０（２）に示すように、ＴＨｚパルスとプローブ光が同じタイミングでＴＨｚ検出素子に入射されると、ＴＨｚパルスの電場時間波形の一部を、プローブ光によって時間的に切り出すことが可能になる。したがって、時間遅延走査ステージ３３によってプローブ光の入射タイミングを少しずつずらしながら、ＴＨｚパルスとプローブ光が重なるタイミングでＴＨｚパルスの電場時間波形の切り出しを行うことにより、ＴＨｚ電場時間波形を再構成することが可能になる。
しかし、従来の機械式時間遅延走査ポンプ・プローブ法を用いた装置構成では、繰り返し周期であるパルス周期に等しい観測時間窓を実現するのに技術的困難を伴う（数メートル以上の時間遅延走査ステージ３３が必要）。また、フェムト秒レーザー光によって発生させたＴＨｚパルスの周期を正確かつ十分に安定化することも困難である。このような理由から、従来の機械式時間遅延走査ポンプ・プローブ法を用いた装置で、本発明のフーリエ変換型分光法を実施するのは容易でない。【００５８】
そこで、繰り返し周波数がわずかに異なるように予め制御された２台のフェムト秒レーザー（各々の繰り返し周波数ｆ１，ｆ２）をそれぞれポンプ光とプローブ光に用い、ＴＨｚパルスの発生及び検出に用いることにより、ピコ秒オーダー（図１１（２）に示す１／ｆ１）のＴＨｚパルス波形をマイクロ秒オーダー（図１１（２）に示す１／Δ（Δ＝ｆ１−ｆ２））まで時間スケールを正確に拡大する。その結果、機械式時間遅延走査ステージを必要とせず、オシロスコープ等でパルスを実時間測定することが可能になる。これにより、機械式時間遅延走査に伴う観測時間窓の制限を解消し、繰り返し周期であるパルス周期まで観測時間窓を拡大する。その上、任意の観測時間窓を設定可能であるため、観測時間窓を繰り返し周期であるパルス周期と厳密に一致することが可能になる。また、マイクロ波原子時計を基準としたレーザー制御により、パルス周期を正確かつ十分に安定化制御することが可能である。
さらに、プロット間隔を走査して、スペクトルを重ね合わせ、プロット間隙の補完により、更なるスペクトル分解能の向上とスペクトル確度を高めることができる。【００５９】
図１２は、実施例１のＴＨｚ−ＴＤＳ装置の構成図を示している。
２台のフェムト秒ファイバーレーザー（５１，５２）（平均パワー５００ｍＷ、中心波長１５５０ｎｍ、パルス幅５０ｆｓ）のパルス周波数（ｆ１，ｆ２）をルビジウム原子時計５３に位相同期することにより、パルス周波数（ｆ１，ｆ２）およびパルス周波数差がわずかに異なるようにレーザー制御５４する（ｆ１＝250,000,000Ｈｚ,ｆ２＝250,000,050Ｈｚ、Δｆ＝ｆ２−ｆ１＝５０Ｈｚ）。ファイバーレーザー１（５１）をポンプレーザーとして用い、ファイバーレーザー２（５２）をプローブレーザーとして用いる。【００６０】
両レーザー光を非線形光学結晶で波長変換した後、ＴＨｚ発生３４及びＴＨｚ検出３５用の光伝導アンテナ（ＰＣＡ）を用いて、ＴＨｚ−ＴＤＳ光学系を構築する。ここで、ＴＨｚパルスとプローブ光のパルス周期はわずかに異なるので、ＴＨｚパルスとプローブ光がＴＨｚ検出用ＰＣＡで重なるタイミングがパルス毎に自動的にずれていき、図１１（２）のように高速サンプリングされる。その結果、サブピコ秒オーダーのＴＨｚパルス電場時間波形の時間スケールを、任意の時間スケール拡大率（＝ｆ１／Δｆ）で拡大できる。時間スケールがＲＦ領域まで拡大された電流信号がＴＨｚ検出用ＰＣＡから出力されることから、電流電圧変換後、増幅器５７で増幅された時間波形を高速デジタイザー５６で取得する。また、レーザー光の一部は、ＳＦＧ（和周波発生光）相互相関計５５に導かれ、発生したＳＦＧ光は時間原点信号として、デジタイザー５６のトリガー信号に用いられる。また、ルビジウム原子時計５３からの信号は、デジタイザー５６の基準クロック信号として使われる。【００６１】
（ＴＨｚ−ＴＤＳによる低圧ガス分光測定による有用性の評価）
次に、低圧状態の水蒸気をサンプルに用いて、ＴＨｚ−ＴＤＳにおける本発明のフーリエ変換型分光法の有用性を評価した結果を説明する。低圧状態の水蒸気は、回転遷移に伴う鋭い吸収線をＴＨｚ領域において示し、その吸収線幅は圧力によって調節できる（以下、“理論吸収線幅”と称する）。そこで、水蒸気と窒素の混合ガスを、ガスセル（長さ５００ｍｍ、直径４０ｍｍ）に封入し、０．５５７ＴＨｚ付近の理論吸収線幅（回転遷移１１０−＞１０１）が十分に狭線幅となるような低圧状態に保ちながら、スペクトル分解能の評価を行った。【００６２】
まず、理論吸収線幅が２５０ＭＨｚになるように設定された水蒸気サンプル（水蒸気１７０Ｐａ；窒素３２００Ｐａ）を測定した。図１３(ａ)は、繰り返し周期であるパルス周期と観測時間窓が厳密に一致させた状態で計測した吸収スペクトルである（観測された吸収スペクトルの線幅を、以下、“観測吸収線幅”と称する）。図１３（ａ）は、前述の図７（ａ）と等価である。吸収スペクトルは、パルス周波数（＝２５０ＭＨｚ）に等しい間隔のプロットから構成されており、観測吸収線幅は理論吸収線幅から若干拡がっている。これは、プロットの分布が理論吸収線幅に比べて粗であるため、詳細なスペクトル形状が得られていないからである。【００６３】
次に、プロット間の隙部を補完するようにスペクトルを逐次走査した後、スペクトルを重ねることによって得られた微細化スペクトルを図１３(ｂ)に示す。図１３（ｂ）は、前述の図９右側と等価である。２５０ＭＨｚ程度の観測吸収線幅が得られているが、より詳細なスペクトル形状が反映されている。【００６４】
比較データとして、ＦＴ−ＩＲ等で得られたスペクトルのデータ点数を増大させる方法、すなわち、従来法としてよく用いられるゼロ・フィリング法によって取得したスペクトルを図１３(ｃ)に示す。図１３（ｃ）ではプロット数は大幅に増えているが（図１３(ａ)の１０倍）、スペクトル分解能は本質的に向上しておらず、スペクトル形状は図１３(ａ)とほぼ同様であることがわかる。また、回転遷移１１０−＞１０１のＮＡＳＡ（National Aeronautics and Space Administration）のスペクトル・データベース値（０．５５６９４６ＴＨｚ）と比較すると、図１３(ｂ)のみ、スペクトル確度が向上していることが確認できる。【００６５】
次に、繰り返し周期と観測時間窓を厳密に一致させる必要性の確認を行った。
まず、水蒸気による吸収緩和現象が繰り返し周期（４ｎｓ）よりも長く継続するよう、理論吸収線幅が１０ＭＨｚとなるように圧力調節した。低圧状態の水蒸気のサンプルとして、水蒸気と窒素の混合ガス（水蒸気６Ｐａ；窒素１４０Ｐａ）を用いた。理論吸収線幅１０ＭＨｚの吸収現象は、時間領域において、理論吸収線幅の逆数である１００ｎｓの時間をかけて吸収が緩和していくことになる。緩和時間１００ｎｓは、繰り返し周期（４ｎｓ）の２５周期分に相当する。【００６６】
上記条件下で、観測時間窓＝１周期と観測時間窓＝０．９９９９５周期の２つの場合について、微細化吸収スペクトル（図９右側を参照）を求めた。測定結果を図１４に示す。観測時間窓＝１周期の場合、鋭く対称的な吸収スペクトルが得られ、観測吸収線幅は２９ＭＨｚであった（図１４(ａ)を参照）。理論吸収線幅（１０ＭＨｚ）からのズレは、サンプルのガス圧力の不安定性によるものである。一方、観測時間窓＝０．９９９９５周期の場合、スペクトル波形が歪んで観測された（図１４(ｂ)を参照）。これは、図８で説明したような時間的不連続接続が発生したためであると推察する。【００６７】
比較データとして、吸収緩和現象の継続時間を繰り返し周期より短くした場合についても、同様な計測を行った。ここでは、理論吸収線幅が５００ＭＨｚとなるように圧力を調節し（水蒸気１０００Ｐａ；窒素３５００Ｐａ）、緩和現象の継続時間が繰り返し周期（４ｎｓ）の半分となる２ｎｓとなるようにした。測定結果を図１５に示す。観測時間窓＝１周期の場合の測定結果（図１５（ａ）参照）と観測時間窓＝０．９９９９５周期の場合の測定結果（図１５（ｂ）参照）の両方において、理論吸収線幅とほぼ等しい観測吸収線幅を得た。これは、吸収緩和現象の継続時間が繰り返し周期よりも短いため、観測時間窓＝０．９９９９５周期の測定結果である図１５（ｂ）の場合でも、観測時間波形に時間的不連続接続が含まれないことを意味している。これらの結果から、本発明のフーリエ変換型分光法においては、緩和現象が繰り返し周期より長い場合において、繰り返し周期と観測時間窓を厳密に一致させることが極めて重要であることが分かる。【００６８】
次に、達成可能なスペクトル分解能を評価するために、圧力調節により理論吸収線幅を変化させながら吸収スペクトルを計測した。
本発明（離散分布スペクトル間隙の補完無し），本発明（離散分布スペクトル間隙の補完有り）および従来方法（ゼロ・フィリング法）という３つの場合について、理論吸収線幅および観測吸収線幅の圧力依存性を比較した結果を図１６のグラフに示している。ここで、図中の点線は、理論吸収線幅の理論曲線である。【００６９】
本発明（離散分布スペクトル間隙の補完無し）および従来方法（ゼロ・フィリング法）では、１ｋＰａ付近までは圧力の低下に伴い、観測吸収線幅は減少したが、１ｋＰａ以下の圧力では観測吸収線幅は約２５０ＭＨｚのままであった。これは、理論吸収線幅を反映しているのではなく、装置のスペクトル分解能の限界のため、観測吸収線幅が拡がって観測されていることを示している。すなわち、本発明（離散分布スペクトル間隙の補完無し）および従来方法（ゼロ・フィリング法）のスペクトル分解が、理論限界スペクトル分解であるプロット間隔（繰り返し周波数）の２５０ＭＨｚであることを示している。【００７０】
一方、本発明（離散分布スペクトル間隙の補完有り）における観測吸収線幅は、理論吸収線幅とほぼ同様な変化を示している。例えば、今回の最低ガス圧である５Ｐａの実験では、観測されたスペクトル幅は３ＭＨｚであり、理論限界スペクトル分解である繰り返し周波数（２５０ＭＨｚ）を大幅に上回っている。また、低圧領域における理論曲線からのズレは、サンプルのガス圧力の不安定性によるものであり、発明装置のスペクトル分解能の限界を示しているわけではない。本発明のフーリエ変換型分光法においては、理論上、無限小のスペクトル分解（無限大のスペクトル分解能）が達成可能であるが、実際には繰り返し周波数の安定性などによって、スペクトル分解能の限界が決まっている。
【実施例２】
【００７１】
実施例２では、フーリエ変換型分光計測プログラムについて説明する。
図１７は、フーリエ変換型分光計測プログラムの機能ブロック図を示している。
実施例２のフーリエ変換型分光計測プログラムは、周期的な繰り返し現象を観測するプログラムであり、コンピュータを下記のＭ１）〜Ｍ４）の手段として機能させる。
（Ｍ１）時間波形取得手段
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓）と、を一致させて時間波形を取得する。
（Ｍ２）フーリエ変換スペクトル算出手段
上記Ｍ１の時間波形取得手段により取得した時間波形からフーリエ変換スペクトルを算出する。
（Ｍ３）繰り返し周期走査手段
離散分布スペクトルの間隙を補完するために、繰り返し周波数を走査する（繰り返し周波数を変化させる）。
（Ｍ４）スペクトル重畳手段
上記Ｍ３の繰り返し周期走査手段により、繰り返し周期が変化した現象の時間波形を、上記Ｍ１の時間波形取得手段で取得し、上記Ｍ２のフーリエ変換スペクトル算出手段で得られた各スペクトル、すなわち、異なる繰り返し周期の各スペクトルを重畳して、離散分布スペクトルの間隙が補完された微細化スペクトルを取得する。【００７２】
また、図１８は、フーリエ変換型分光計測プログラムの処理フローを示している。
フーリエ変換型分光計測プログラムは、周期的な繰り返し現象を観測すべく、現象の繰り返し周期と観測時間窓を一致させて、現状の時間波形を取得する。そして、取得した時間波形からフーリエ変換スペクトルを算出する。そして、離散分布スペクトルの間隙を補完するために、繰り返し周波数を走査する。離散分布スペクトルの間隙を補完するために必要なだけ走査を行う。
そして、繰り返し周期の走査により、繰り返し周期が変化した現象の時間波形を取得し、それらをフーリエ変換して得られた各スペクトル（異なる繰り返し周期の各スペクトル）を重畳して、離散分布スペクトルの間隙が補完された微細化スペクトルを取得する。【００７３】
図１９は、フーリエ変換型分光計測プログラムを実行するコンピュータ・ハードウェアの内部構成を示す。図１９において、コンピュータ・ハードウェアの内部構成は、ＣＰＵ１１１、ＲＯＭ１１２、ハードディスク１１３、キーボード１１４、マウス１１５、ディスプレイ１１６、光学式ドライブ１１７、ＲＡＭ１１８を備えており、システムバス１１９に接続されている。ＲＯＭ１１２はコンピュータを起動するためのブートアッププログラム等のプログラムを記憶する。ＲＡＭ１１８はフーリエ変換型分光計測プログラムの命令を一時的に記憶するとともに一時記憶空間を提供する。ハードディスク１１３はフーリエ変換型分光計測プログラム、システムプログラム及びデータを記憶する。キーボード１１４、マウス１１５は、コンピュータの操作者からの命令を受け付ける。ディスプレイ１１６は、フーリエ変換型分光計測プログラムにおける出力データを表示する。なお、コンピュータは、ネットワークへの接続を行うためのネットワークインターフェース（図示しない）を含んでもよい。【００７４】
フーリエ変換型分光計測プログラムは、核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析（ＦＴ−ＭＳ）装置，フーリエ変換光スペクトラム・アナライザー装置，フーリエ変換スペクトラム・アナライザーに標準装備されているコンピュータ・ハードウェア（ＣＰＵ，ＲＯＭ，ハードディスク等）を共用して使用することができ、標準搭載プログラムの追加オプションとして搭載できる。
【実施例３】
【００７５】
上述の実施例１では、低圧状態の水蒸気をサンプルに用いて、ＴＨｚ−ＴＤＳにおける本発明のフーリエ変換型分光法の有用性を評価した結果を説明した。
実施例３では、アセトニトリルガス（ＣＨ３ＣＮ）の吸収線の計測を例に挙げ、ＴＨｚ−ＴＤＳにおける本発明のフーリエ変換型分光法の有用性を説明する。
アセトニトリルは、星間物質のひとつであり、宇宙天文学の分野でも宇宙の物理状態を知る目的で、アセトニトリルの吸収線を計測する研究は重要と考えられている。アセトニトリルは対称コマ型分子であり、回転遷移による吸収線の周波数ｖは、下記式式（８）で表される。【００７６】
【数８】
2014034085.tif 000011
【００７７】
ここで、Ｂは:対称軸周りの回転定数であり、ＤＪ，ＤＪＫは分子回転による遠心力ひずみ定数であり、Ｊは全角運動量回転量子数、Ｋは分子対称軸方向回転量子数である。
ＤＪ，ＤＪＫは、係数Ｂに比べて非常に小さいので、広い範囲で観察すると、吸収線は周波数２Ｂ（＝１８．４ＧＨｚ）でほぼ等間隔で現れることから、下記数式（９）で表される。【００７８】
【数９】
2014034085.tif 000012
【００７９】
さらに、狭い範囲で観察すると、量子数Ｋによる微細構造も現れるようになる。従来のテラヘルツ分光装置ではこれらを一度に観測することは困難であった。このため、量子数Ｊによる構造を観測する際には、広帯域のＴＨｚ−ＴＤＳ装置を用い、量子数Ｋによる微細構造を観測する際は、高分解能のＣｏｎｔｉｎｕｏｕｓ−Ｗａｖｅ（ＣＷ）−ＴＨｚ分光装置を用いる必要があった。周波数スケールのダイナミックレンジ（スペクトル分解とスペクトル範囲の比）が極めて高い本発明のフーリエ変換型分光法を用いれば、これらの両構造を１台のテラヘルツ分光装置で観測することが可能である。【００８０】
本実施例では、アセトニトリルガスの微細構造が圧力拡がりによって重畳してしまうことを避けるために、圧力を３０Ｐａに設定して計測を実施した。
Ｊ＝３４〜３５遷移付近の微細構造を観測するために、０．６５ＴＨｚ付近でスペクトルの各プロット間隔が１２．５ＭＨｚずつになるように、モード同期周波数を２０回変化させた。
図２０に、スペクトル全体の振幅スペクトルを示す。０．３〜１．０ＴＨｚ付近にアセトニトリル分子の周期的吸収線が確認できる。【００８１】
次に、吸収スペクトルを算出した。図２１に、０．６〜０．７ＴＨｚの範囲での吸収スペクトルを示す。アセトニトリル分子の６本の吸収線が、予測値とほぼ同じ約１８．４ＧＨｚの間隔で存在していることが確認できる。【００８２】
さらに、図２２に、Ｊ＝３４〜３５遷移付近（０．６４ＴＨｚ付近）を拡大した吸収スペクトルを示す。図２２では、各々のプロット（図中に○で示している）をローレンツ関数によってフィッティングした曲線を示している。また、Ｋ＝０〜１０の破線の延長上に描くバーは、ＮＡＳＡ（National Aeronautics and Space
Administration）データベースから引用したアセトニトリル分子の絶対周波数と積分強度を示している。
図２２から、量子数Ｋによる微細構造が確認することができ、Ｋ＝２〜１０を識別することができることがわかる。【００８３】
一方で、周波数間隔が１２．３ＭＨｚである二つの吸収線Ｋ＝０，１については、スペクトル・プロット間隔が１２．５ＭＨｚの実験条件では分離・識別することができていない。図２２では分離・識別することができない二つの吸収線Ｋ＝０，１を識別するために、さらに０．６４３２ＴＨｚ周辺でスペクトル・プロット間隔が１．２５ＭＨｚずつ２０回ずれるように、モード同期周波数を変化させた。【００８４】
図２３は、モード同期周波数を変化させて、複数の離散分布スペクトルを重畳して得られた吸収スペクトルであって、０．６４３２ＴＨｚ付近を拡大した吸収スペクトルを示している。図２３では、周波数間隔が１２．３ＭＨｚの二つの吸収線Ｋ＝０，１を、はっきり区別できていることがわかる。このように、スペクトル・プロット間隔を狭めれば狭めるほどスペクトル分解能が向上するという結果は、本発明において各スペクトル・プロットが無限小のスペクトル分解（無限大のスペクトル分解能）を達成していることを意味している。【００８５】
また、スペクトル確度においても周波数走査ステップ１．２５ＭＨｚの半分の範囲内（０．６２５ＭＨｚ）で一致しており、高い精度が得られている。
上記のことから、本発明のフーリエ変換型分光法であって、モード同期周波数を変化させて、複数の離散分布スペクトルを重畳して、取得した離散分布スペクトルの各プロットの間の間隙を補完する方法では、ＴＨｚ−ＴＤＳの広帯域特性とＣＷ−ＴＨｚ分光法の高いスペクトル分解能を併せ持ち、さらにスペクトル確度がマイクロ波周波数標準によって担保されていることがわかるであろう。【００８６】
（その他の実施例）
上記の実施例では、テラヘルツ時間領域分光装置（ＴＨｚ−ＴＤＳ）に適用する例を説明したが、その他、核磁気共鳴分光装置（ＮＭＲ），核磁気共鳴画像装置（ＭＲＩ），フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析装置（ＦＴ−ＭＳ），フーリエ変換光スペクトラム・アナライザー装置，フーリエ変換スペクトラム・アナライザーに、本発明のフーリエ変換型分光法を適用し、またフーリエ変換型分光計測プログラムを搭載できる。
【産業上の利用可能性】
【００８７】
本発明は、テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析装置（ＦＴ−ＭＳ），フーリエ変換光スペクトラム・アナライザー装置，フーリエ変換スペクトラム・アナライザーなどのフーリエ変換型分光装置に利用できる。
【符号の説明】
【００８８】
１０〜１３          時間波形
２０                周波数信号
３１，４１，４２    フェムト秒レーザー
３２                ビームスプリッター
３３                時間遅延走査ステージ
３４                ＴＨｚ発生素子
３５                ＴＨｚ検出素子
３６，３６ａ，３６ｂ      集光レンズ
５１，５２          ファイバーレーザー
５３                ルビジウム原子時計
５４                レーザー制御装置
５５                ＳＦＧ相互相関計
５６                デジタイザー
５７                増幅器
５８                光ファイバー
６０                ガスセル
６１                光ファイバー端
６２                集光レンズ
【図１】
2014034085.tif 000013
【図２】
2014034085.tif 000014
【図３】
2014034085.tif 000015
【図４】
2014034085.tif 000016
【図５】
2014034085.tif 000017
【図６】
2014034085.tif 000018
【図７】
2014034085.tif 000019
【図８】
2014034085.tif 000020
【図９】
2014034085.tif 000021
【図１０】
2014034085.tif 000022
【図１１】
2014034085.tif 000023
【図１２】
2014034085.tif 000024
【図１３】
2014034085.tif 000025
【図１４】
2014034085.tif 000026
【図１５】
2014034085.tif 000027
【図１６】
2014034085.tif 000028
【図１７】
2014034085.tif 000029
【図１８】
2014034085.tif 000030
【図１９】
2014034085.tif 000031
【図２０】
2014034085.tif 000032
【図２１】
2014034085.tif 000033
【図２２】
2014034085.tif 000034
【図２３】
2014034085.tif 000035
【手続補正書】
【提出日】平成26年4月23日(2014.4.23)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式１が成立し、
前記観測時間窓Ｔを前記現象の繰り返し周期と一致させる替わりに、前記繰り返し周期よりも短い時間で現象を観測して時間波形を取得し、取得した時間波形データに対して繰り返し周期に一致させるように時間波形データを補間し、繰り返し周期と一致した時間波形を取得する、
ことを特徴とするフーリエ変換型分光法。
【数１】
2014034085.tif 000036
【請求項２】
（削除）
【請求項３】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式２が成立し、
前記繰り返し周期を変化させ、繰り返し周期を変化させた後の時間波形を取得し離散分布スペクトルを算出し、異なる繰り返し周期の各離散分布スペクトルを重畳する、
ことを特徴とするフーリエ変換型分光法。
【数２】
2014034085.tif 000037
【請求項４】
前記繰り返し周期を、周波数標準器を基準として安定化させた、請求項１      に記載のフーリエ変換分光法。
【請求項５】
前記繰り返し周期の整数分の１の時間間隔で時間波形データを離散化して取得した、或いは、前記繰り返し現象の周期を時間波形データの離散化時間間隔の整数倍に設定して取得した、前記繰り返し周期１周期分の離散化した時間波形データを離散フーリエ変換することにより離散分布スペクトルを得る、ことを特徴とする請求項１〜３のいずれかに記載のフーリエ変換分光法。
【請求項６】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）に用いられ、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換型分光法。
【請求項７】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）に用いられ、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換型分光法。
【請求項８】
請求項１〜４のいずれかに記載のフーリエ変換型分光法を用いた、
核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析（ＦＴ−ＭＳ）装置，フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーのいずれかのフーリエ変換型分光装置。
【請求項９】
周期的な繰り返し現象を観測するフーリエ変換型分光装置において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段と、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段と、
前記繰り返し周期を変化させる手段と、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段と、
を備え、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式３が成立する、
ことを特徴とするフーリエ変換型分光装置。
【数３】
2014034085.tif 000038
【請求項１０】
（削除）
【請求項１１】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）を行うテラヘルツ時間領域分光装置において、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザー手段の各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザー手段をそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、
ことを特徴とする請求項９に記載のフーリエ変換型分光装置。
【請求項１２】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）を行うフーリエ変換型赤外分光装置において、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、
ことを特徴とする請求項９に記載のフーリエ変換型分光装置。
【請求項１３】
周期的な繰り返し現象を観測するフーリエ変換型分光計測プログラムであって、
コンピュータを、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段、
前記繰り返し周期を変化させる手段、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段、
として機能させるプログラムであって、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式４が成立する、
ことを特徴とするフーリエ変換型分光計測プログラム。
【数４】
2014034085.tif 000039
【手続補正書】
【提出日】平成27年2月18日(2015.2.18)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式１が成立し、
前記観測時間窓Ｔを前記現象の繰り返し周期と一致させる替わりに、前記繰り返し周期よりも短い時間で現象を観測して時間波形を取得し、取得した時間波形データに対して繰り返し周期に一致させるように時間波形データを補間し、繰り返し周期と一致した時間波形を取得する、ことを特徴とするフーリエ変換型分光法。
（数１）
ｆｎ＝ｎ／Ｔ      ・・・（式１）
【請求項２】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式２が成立し、
前記繰り返し周期を変化させ、繰り返し周期を変化させた後の時間波形を取得し離散分布スペクトルを算出し、異なる繰り返し周期の各離散分布スペクトルを重畳する、ことを特徴とするフーリエ変換型分光法。
（数２）
ｆｎ＝ｎ／Ｔ      ・・・（式２）
【請求項３】
前記繰り返し周期を、周波数標準器を基準として安定化させた、請求項１又は２に記載のフーリエ変換分光法。
【請求項４】
前記繰り返し周期の整数分の１の時間間隔で時間波形データを離散化して取得した、或いは、前記繰り返し現象の周期を時間波形データの離散化時間間隔の整数倍に設定して取得した、前記繰り返し周期１周期分の離散化した時間波形データを離散フーリエ変換することにより離散分布スペクトルを得る、ことを特徴とする請求項１〜３のいずれかに記載のフーリエ変換分光法。
【請求項５】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）に用いられ、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換型分光法。
【請求項６】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）に用いられ、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換型分光法。
【請求項７】
請求項１〜４のいずれかに記載のフーリエ変換型分光法を用いた、核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析（ＦＴ−ＭＳ）装置，フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーのいずれかのフーリエ変換型分光装置。
【請求項８】
周期的な繰り返し現象を観測するフーリエ変換型分光装置において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段と、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段と、
前記繰り返し周期を変化させる手段と、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段と、
を備え、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式３が成立する、ことを特徴とするフーリエ変換型分光装置。
（数３）
ｆｎ＝ｎ／Ｔ      ・・・（式３）
【請求項９】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）を行うテラヘルツ時間領域分光装置において、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザー手段の各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザー手段をそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項８に記載のフーリエ変換型分光装置。
【請求項１０】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）を行うフーリエ変換型赤外分光装置において、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項８に記載のフーリエ変換型分光装置。
【請求項１１】
周期的な繰り返し現象を観測するフーリエ変換型分光計測プログラムであって、
コンピュータを、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段、
前記繰り返し周期を変化させる手段、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段、
として機能させるプログラムであって、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式４が成立する、ことを特徴とするフーリエ変換型分光計測プログラム。
（数４）
ｆｎ＝ｎ／Ｔ      ・・・（式４）
【手続補正書】
【提出日】平成28年4月25日(2016.4.25)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得するか、
もしくは、前記観測時間窓Ｔを前記現象の繰り返し周期と一致させる替わりに、前記繰り返し周期よりも短い時間で現象を観測して時間波形を取得し、取得した時間波形データに対して繰り返し周期に一致させるように時間波形データを補間し、繰り返し周期と一致した時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットが、ｆｎを各プロットの周波数、ｎをプロットの次数（整数）として、下記数式１を満足することを特徴とするフーリエ変換型分光法。
（数１）
ｆｎ＝ｎ／Ｔ      ・・・（式１）
【請求項２】
周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得し、
上記の離散分布スペクトルの各プロットが、ｆｎを各プロットの周波数、ｎをプロットの次数（整数）として、下記数式２を満足し、
前記繰り返し周期を変化させ、繰り返し周期を変化させた後の時間波形を取得し離散分布スペクトルを算出し、異なる繰り返し周期の各離散分布スペクトルを重畳する、ことを特徴とするフーリエ変換型分光法。
（数２）
ｆｎ＝ｎ／Ｔ      ・・・（式２）
【請求項３】
前記繰り返し周期を、周波数標準器を基準として安定化させた、請求項１又は２に記載のフーリエ変換分光法。
【請求項４】
前記繰り返し周期の整数分の１の時間間隔で時間波形データを離散化して取得した、或いは、前記繰り返し現象の周期を時間波形データの離散化時間間隔の整数倍に設定して取得した、前記繰り返し周期１周期分の離散化した時間波形データを離散フーリエ変換することにより離散分布スペクトルを得る、ことを特徴とする請求項１〜３のいずれかに記載のフーリエ変換分光法。
【請求項５】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）に用いられ、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換型分光法。
【請求項６】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）に用いられ、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザーを用い、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項１〜４のいずれかに記載のフーリエ変換型分光法。
【請求項７】
請求項１〜４のいずれかに記載のフーリエ変換型分光法を用いた、核磁気共鳴分光（ＮＭＲ）装置，核磁気共鳴画像（ＭＲＩ）装置，テラヘルツ時間領域分光（ＴＨｚ−ＴＤＳ）装置，フーリエ変換型赤外分光（ＦＴ−ＩＲ）装置，フーリエ変換質量分析（ＦＴ−ＭＳ）装置，フーリエ変換光スペクトラム・アナライザー，フーリエ変換スペクトラム・アナライザーのいずれかのフーリエ変換型分光装置。
【請求項８】
周期的な繰り返し現象を観測するフーリエ変換型分光装置において、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段と、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段と、
前記繰り返し周期を変化させる手段と、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段と、
を備え、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式３が成立する、ことを特徴とするフーリエ変換型分光装置。
（数３）
ｆｎ＝ｎ／Ｔ      ・・・（式３）
【請求項９】
テラヘルツ時間領域分光法（ＴＨｚ−ＴＤＳ）を行うテラヘルツ時間領域分光装置において、
テラヘルツ時間領域分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザー手段の各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザー手段をそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をＴＨｚ発生用ポンプ光として用い、他方のフェムト秒レーザーの出力光をＴＨｚ検出用プローブパルス光として用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項８に記載のフーリエ変換型分光装置。
【請求項１０】
フーリエ変換型赤外分光法（ＦＴ−ＩＲ）を行うフーリエ変換型赤外分光装置において、
フーリエ変換型赤外分光法に用いる光源として、レーザーパルスの繰り返し周波数（モード同期周波数）が異なる２台のフェムト秒レーザー手段を備え、
前記２台のフェムト秒レーザーの各々のモード同期周波数が周波数標準器を基準として安定化され、かつ、モード同期周波数の差が所定の一定値に保持されるように２台のフェムト秒レーザーをそれぞれ独立して同時に制御し、
一方のフェムト秒レーザーの出力光をサンプル計測用赤外光発生のために用い、他方のフェムト秒レーザーの出力光をヘテロダイン干渉用局部発振器赤外光発生のために用い、
前記繰り返し周期の安定化を図る、ことを特徴とする請求項８に記載のフーリエ変換型分光装置。
【請求項１１】
周期的な繰り返し現象を観測するフーリエ変換型分光計測プログラムであって、
コンピュータを、
現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得する手段、
取得した時間波形からフーリエ変換することにより、前記観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する手段、
前記繰り返し周期を変化させる手段、
異なる繰り返し周期の複数の離散分布スペクトルを重畳する手段、
として機能させるプログラムであって、
上記の離散分布スペクトルの各プロットは、ｆｎが各プロットの周波数、ｎがプロットの次数（整数）として、下記数式４が成立する、ことを特徴とするフーリエ変換型分光計測プログラム。
（数４）
ｆｎ＝ｎ／Ｔ      ・・・（式４）
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１５
【補正方法】変更
【補正の内容】
【００１５】
上記目的を達成すべく、本発明のフーリエ変換型分光法は、周期的な繰り返し現象を観測するフーリエ変換型分光法（フーリエ変換型周波数解析法）において、現象の繰り返し周期と、現象の時間波形を観測する時間幅（観測時間窓Ｔ）と、を一致させて時間波形を取得し、フーリエ変換することにより、観測時間窓Ｔの逆数を周波数データ間隔とする離散分布スペクトルを取得する。
上記の離散分布スペクトルの各プロットが、ｆｎを各プロットの周波数、ｎをプロットの次数（整数）として、下記数式３を満足する。
【国際調査報告】
2014034085.tif 000040
2014034085.tif 000041
2014034085.tif 000042