(57)【要約】【課題】本発明が解決する課題は、従来のＳＴＦＴによる信号処理よりも格段に高速処理が可能で、かつＣＰＵへの負担が軽い、ＯＦＤＲ方式の計測装置が観測した干渉信号のあらたな高効率の信号処理手法を提供することにある。【解決手段】本発明のＯＦＤＲ方式光ファイバ計測の信号処理方法は、ブラッグ波長の初期信号とシフトされた信号から位相微分方式やＺ変換方式等の手法により群遅延を演算し、ブラッグ波長シフトの分布を求めることを特徴とする。【選択図】図１
2015105909.tif 000002
【特許請求の範囲】
【請求項１】
ＦＢＧを１本の光ファイバに複数個配置したＯＦＤＲ方式の光ファイバ計測装置において、観測されたブラッグ波長の初期信号とシフトされた信号から群遅延を演算してブラッグ波長シフトの分布を求めることを特徴とするＯＦＤＲ方式の光ファイバ計測装置の信号処理方法。
【請求項２】
前記群遅延の演算は、初期信号とシフトされた信号の位相を検出し、位相のアンラッピングを行い、それぞれの位相誘導計算をしてその差から時間群遅延計算を行う位相微分方式を用いるものである請求項１に記載のＯＦＤＲ方式光ファイバ計測の信号処理方法。
【請求項３】
前記群遅延の演算は、初期信号とシフトされた信号の位相を検出し、それぞれにフーリエ変換を施し、周波数応答を求め、逆フーリエ変換によりインパルス応答を算出し、該インパルス応答を用いて群遅延を計算するＺ変換方式を用いるものである請求項１に記載のＯＦＤＲ方式光ファイバ計測の信号処理方法。
【請求項４】
前記群遅延の演算は、初期信号とシフトされた信号の位相を検出し、それぞれにフーリエ変換を施し、周波数応答を求め、畳み込み積分によりインパルス応答を算出し、該インパルス応答を用いて群遅延を計算するＺ変換方式を用いるものである請求項１に記載のＯＦＤＲ方式光ファイバ計測の信号処理方法。
【請求項５】
インパルス応答ｈ(ｋ)およびｋｈ(ｋ)に窓関数をかけて特定の区間における次の関数を取り出し、その抜き出す位置を変えてＦＦＴを繰り返し計算し、それらの総和で群遅延を計算する請求項４に記載のＯＦＤＲ方式光ファイバ計測の信号処理方法。
【請求項６】
ＦＢＧを１本の光ファイバに複数個配置して、光の干渉強度周期的変化を利用するＯＦＤＲ方式の計測装置において、ブラッグ波長の初期信号とシフトされた信号から群遅延を演算し、ブラッグ波長シフトの分布を求める信号処理手段を備えたことを特徴とするＯＦＤＲ方式光ファイバ計測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光周波数領域リフレクトメトリを用いた光ファイバセンシング技術における信号処理に関する。
【背景技術】
【０００２】
光ファイバのコアに紫外線を照射して、規則的な縞（回折格子）を作ると、光ファイバの中を透過していく光の様々な波長の中でその縞の周期（屈折率）に合う波長だけが跳ね返（反射）される。この縞の周期は歪みや温度によって伸縮するため、このような光ファイバはセンサ（Fiber Bragg Grating：ＦＢＧ）になる。このＦＢＧを１本の光ファイバに複数個配置した周波数多重化（Optical Frequency Domain Reflectometry：ＯＦＤＲ）方式は、ファイバ全体を使って静歪み、変位、圧力、温度を計測するのに優れたものとして知られている。
【０００３】
ＦＢＧは通常、反射スペクトルのシフトを観測することでひずみ（あるいは温度）を計測する。ＦＢＧの反射スペクトルを観測する最も代表的な機器がスペクトラムアナライザである。代表例としてアンリツ社のスペクトラムアナライザは掃引速度が数Hzである。この場合、掃引速度がすなわち測定速度（周波数）となる。同社のＦＢＧセンサモニタはＦＢＧ計測用に特化した使用となっており、これは1.2kHz以上の測定速度を有する。
【０００４】
ＦＢＧの反射スペクトルそのものを直接観測するのではなく、反射信号の強度変調を追跡することで被測定物のひずみ・振動を計測する例もある。図１２のようにＦＢＧからの反射信号を受光器で観測する。受光器には図１１に示すような波長依存の光起電力特性を持たせておく。これにより被測定物のひずみ変化により生じるブラッグ波長シフトが起電力の変化として観測される。この手法により最大250Hzでの測定が実証された。
このようにＦＢＧを利用した高速測定の例は見られるが、これらはあくまでＦＢＧを単点測定用に利用したものである。ＦＢＧセンサモニタの場合はＦＢＧを多重化して計測することで準分布型の測定が可能である。
【０００５】
特許文献１に示された例で、ＯＦＤＲ方式の多点歪計測装置の概要を示す。図１３はその基本構成を示す図で、主要構成品は以下の通りである。
Ｄ1,Ｄ2,Ｄ3,Ｄ4：  光強度ディテクタ
Ｃ1：  ５％−９５％広帯域カップラ
Ｃ2-Ｃ5：  ５０％−５０％広帯域カップラ
Ｒ1-Ｒ4：  全反射が起きるような終端
ＦＣ／ＰＣ1-3：  光ファイバ結合用FC/PCコネクタ
ＢＮＣ1-3：  光強度電圧変換出力用BNCコネクタ
ＦＢＧ1-2：  ＦＢＧを５個配置した光ファイバ
計測は波長可変型光源よりレーザー光を供給し、波長をスイープさせながらディテクタ部で反射光の強度を測定する。データ処理には高速フーリエ変換やデジタルのフィルターを使用するので、一定間隔の光波数毎に計測を行う必要がある。ディテクタＤ1で計測される光強度は、次式のように表される光波数間隔△ｋ（＝π／ｎＬ）で周期的に変化する。光波数に対して周期的に変化するディテクタＤ1での光強度をトリガーとして使用して、センサ部である各ＦＢＧからの反射光の強度をディテクタＤ2、Ｄ3で計測する。
【０００６】
本発明者らが開発したＯＦＤＲ式分布測定システムは、光源にその波長掃引速度が最高100ｎｍ/secのものを用いており、当研究室のシステムの測定速度は10Hzとなっている。観測した干渉信号の信号処理効率も測定速度に大きく影響する。現在のＯＦＤＲ測定では波長掃引幅８ｎｍ，リファレンスアームであるＬ長（ディテクタＤ1信号用の二つの全反射端の行路長差）20ｍの条件で、一回の計測でバイナリ形式のファイルが0.76MBのサイズとなる。このサイズの大きさがデータ保存速度に影響し、測定速度に影響している。またこれは10Hzで５分間測定するとゆうに２GBを超える計算になり、長時間測定の場合は大きなストレージを用意しなくてはならない。実験室で測定データの取得を行った後に、解析作業を行いひずみ情報を算出する。ブラッグ波長分布を出力する解析は一つの測定データに対して１秒弱の時間がかかることになる。
現状では実験室でリアルタイムにひずみをモニタリングするのは難しい。というのは測定データの保存と解析、ひずみ情報の信号処理を並行して高速で行うことが出来ないためである。これらを並行して行うには、経験的に、一度の計測に一秒以上はかかってしまう。上述したように測定速度が10Hzといったときには測定データの保存だけを意味しており、ひずみ変動については実験後に解析を行ってからようやく把握している状況にある。
【０００７】
本発明者らのＯＦＤＲ式分布測定システムの信号処理は短時間フーリエ変換（STFT：Short time Fourier transform）に基づいた手法を採用している。短時間フーリエ変換は測定データの解析、ひずみ情報の信号処理を高速に行うための手法で、特許文献１に提示されている。ただし、特許文献１ではスペクトログラム解析との用語で記載されている。図１０にＳＴＦＴの概念図を示す。干渉信号に対して、ある波長を中心波長とした窓関数をかけることで干渉信号の一定波数区間を取り出す。取り出された窓関数内の干渉信号区間について高速フーリエ変換（FFT：fast Fourier transformation）解析を行うことで周波数スペクトルが算出される。周波数はファイバ位置に対応するので、ＦＦＴによりとある波長における反射光強度の空間分布を取得したことになる。すなわち、スペクトログラムにおける縦方向断面の一つである。そして窓関数をオーバーラップさせながら次々にずらして同様にＦＦＴ演算をしていくことで、すべての波数（波長）に亘って反射強度空間分布を算出できる。すなわちスペクトログラムを波長方向に掃引するように描画できる。なお、本手法では窓関数にハミングウィンドウを採用している。
表１に本研究で用いたことのあるＳＴＦＴ演算の条件例を示す。この場合はデータ点数7455個の区間に対するＦＦＴを1400回しなくてはならない。ＳＴＦＴはＦＦＴを数多く演算する、ＣＰＵ高負荷な演算である。
【表１】
2015105909.tif 000003
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００５−１４７９００号公報  「ＯＦＤＲ方式の歪連続分布計測装置」  平成１７年６月９日公開
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明が解決する課題は、従来のＳＴＦＴによる信号処理よりも格段に高速処理が可能で、かつＣＰＵへの負担が軽い、ＯＦＤＲ方式の計測装置におけるあらたな高効率の信号処理手法を提案することにある。
【課題を解決するための手段】
【００１０】
本発明のＯＦＤＲ方式光ファイバ計測の信号処理方法は、ＦＢＧを１本の光ファイバに複数個配置したＯＦＤＲ方式の光ファイバ計測装置において、観測されたブラッグ波長の初期信号とシフトされた信号から群遅延を演算してブラッグ波長シフトの分布を求めることを特徴とする。
前記群遅延の演算の第１の形態は、前記群遅延の演算は、初期信号とシフトされた信号の位相を検出し、位相のアンラッピングを行い、それぞれの位相誘導計算をしてその差から時間群遅延計算を行う位相微分方式を用いるものである請求項１に記載のＯＦＤＲ方式光ファイバ計測の信号処理方法。
前記群遅延の演算の第２の形態は、初期信号とシフトされた信号の位相を検出し、それぞれにフーリエ変換を施し、周波数応答を求め、逆フーリエ変換によりインパルス応答を算出し、該インパルス応答を用いて群遅延を計算するＺ変換方式を用いるものとした。
前記群遅延の演算の第３の形態は、初期信号とシフトされた信号の位相を検出し、それぞれにフーリエ変換を施し、周波数応答を求め、畳み込み積分によりインパルス応答を算出し、該インパルス応答を用いて群遅延を計算するＺ変換方式を用いるものとした。
上記形態において、ノイズ処理のためインパルス応答ｈ(ｋ)およびｋ ｈ(ｋ)に窓関数をかけて特定の区間における次の関数を取り出し、その抜き出す位置を変えてＦＦＴを繰り返し計算し、それらの総和で群遅延を計算するものとした。
【００１１】
本発明のＯＦＤＲ方式光ファイバ計測装置は、ＦＢＧを１本の光ファイバに複数個配置して、光の干渉強度周期的変化を利用するＯＦＤＲ方式の計測装置において、ブラッグ波長の初期信号とシフトされた信号から群遅延を演算し、ブラッグ波長シフトの分布を求める信号処理手段を備えたことを特徴とする。
【発明の効果】
【００１２】
本発明のＯＦＤＲ方式光ファイバ計測の信号処理方法は、従来のＳＴＦＴによる信号処理よりも格段に高速処理が可能で、かつＣＰＵへの負担が軽く、ＯＦＤＲ方式の計測装置における高効率の信号処理が実現できる。
Ｚ変換による手法では、波長掃引幅に依存しない傾向を示しており、従来のＳＴＦＴによる手法よりも高い空間分解能を示す。また、計測時間については位相微分による群遅延手法と比べてもこのＺ変換による手法の方が数倍速くなっている。
また、畳み込み積分によりインパルス応答を算出し、該インパルス応答を用いて群遅延を計算するＺ変換方式においては、従来のＳＴＦＴ手法の数十倍の計算効率が得られる。
【００１３】
本発明のＯＦＤＲ方式光ファイバ計測装置は、ブラッグ波長の初期信号とシフトされた信号から群遅延を演算し、ブラッグ波長シフトの分布を求める信号処理手段を備えたものであるから、従来のSTFTによる計測装置より、信号処理よりも格段に高速処理が可能で、かつＣＰＵへの負担が軽く、高効率の信号処理が実現可能なＯＦＤＲ方式の計測装置を提供することができる。
【図面の簡単な説明】
【００１４】
【図１】ブラッグスペクトルの群遅延を包絡線シフトで説明する図である。
【図２】ＳＴＦＴと群遅延計算による出力の違いを説明する図である。
【図３】周波数応答の位相分布を示す図である。
【図４】位相アンラッピングを説明する図である。
【図５】位相微分を用いた群遅延計算のフロー図である。
【図６】逆フーリエ変換によりインパルス応答を算出するＺ変換を用いた群遅延計算のフロー図である。
【図７】畳み込み積分によりインパルス応答を求めるＺ変換を用いた群遅延計算のフロー図である。
【図８】ノイズ処理を加えた畳み込み積分によりインパルス応答を求めるＺ変換を用いた群遅延計算のフロー図である。
【図９】ノイズ処理の過程で行う窓関数を用いた区間取り出しを説明する図である。
【図１０】ＳＴＦＴの概念図である。
【図１１】ブラッグスペクトルの光起電力と波長シフトの関係を説明する図である。
【図１２】ＦＢＧによるひずみ計測システムを説明する図である。
【図１３】ＯＦＤＲ方式の多点歪計測装置の概要を示す図である。
【発明を実施するための形態】
【００１５】
本発明が提案する新たなＯＦＤＲ信号処理である群遅延演算の原理・特徴について説明する。ＯＦＤＲにより生成される干渉信号は、各位置に準じた周波数を有する搬送波がのったブラッグスペクトルの重ね合わせである。とある位置（とある周波数）におけるブラッグスペクトルを取り出すと図１のようになる。搬送波の周波数は同じ位置を見ている限り一定であり、ブラッグ波長シフトは信号の包絡線（振幅変調）シフトで表される。この包絡線シフトが“群遅延”に対応する。
ＯＦＤＲが取得するＦＢＧ反射光の干渉信号は、ＦＢＧ長を微小区間に区切って考えると、各微小区間における反射スペクトルの重ね合わせである。各々の反射スペクトルは区間位置に線形相関した周波数を持つ搬送波がブラッグスペクトルに乗算されたものである。すなわち反射スペクトルの包絡線はブラッグスペクトルに対応する。ひずみ変動を経験する前後の包絡線シフトを算出すればひずみ変動を計測できるが、包絡線シフトは群遅延に相当する。従って周波数領域で群遅延を算出すれば、ＦＢＧ内各区間におけるひずみ変動を分布的に取得できる。
【００１６】
ひずみ変化前の干渉信号をｆ１(ｋ)、ひずみ変化後の干渉信号をｆ２(ｋ)とすると、これら二つはインパルス応答ｈ(ｋ)を用いて
【数１】
2015105909.tif 000004
と表現される。フーリエ変換によるｆ１(ｋ)，ｆ２(ｋ)の周波数スペクトルをＦ１(ω)，Ｆ２(ω)とすれば、それらの関係は周波数応答Ｈ(ｋ)を用いて
【数２】
2015105909.tif 000005
と表される。ここで周波数応答Ｈ(ｋ)の位相をφ(ω)とすれば、群遅延Ｔｇ(ω)は
【数３】
2015105909.tif 000006
となる。群遅延は周波数の関数として算出される。周波数が位置に対応していることを考えれば、上述のように群遅延を求めることでブラッグ波長シフトの分布を求めることができる。
【００１７】
ＳＴＦＴと群遅延計算では演算の意味合いでどのような差があるのかを説明すれば、ＳＴＦＴは反射光強度の空間分布を波長成分全域に亘って算出する。わかりやすくするため軸を違う方向からみれば、ブラッグ“スペクトル”の空間分布を算出する演算である。一方群遅延計算ではブラッグ“波長シフト”の空間分布を演算する。図２のように、ＳＴＦＴによる出力の方がスペクトル形状を観測できるために情報量が豊富だが、演算コストも高くなる。もちろんブラッグ波長シフトを出力するためにはシフト前と後の二つの信号に対してＳＴＦＴ演算を施さなくてはならない。それに対して群遅延計算ではブラッグ波長シフトしか算出されないためにひずみ変化といった情報しか引き出せないが、その分演算コストは低くなる。
【００１８】
本発明では従来のＳＴＦＴ演算による信号処理に代わる信号処理手法として群遅延計算に基づく新たなアプローチを提案する。本発明では群遅延計算に基づく信号処理手法として二種類のアプローチを提案する。一つは「位相微分を用いる方法」であり、これは直接的に群遅延を計算する手法である。もう一つは「Ｚ変換を用いる方法」である。これは微分演算をＺ領域に持ち込む手法である。
まず、位相微分を用いる方法であるが、この手法では式に示したように周波数応答の位相微分を行うことで群遅延を計算する。ここで位相アンラッピングの手間と不確定性について考えなければならない。周波数応答の位相は図３で例示されるように−π〜πの範囲で周期的に変動する。位相の微分を行うためには位相の“アンラッピング”を行わなくてはいけない。つまり位相が−πかπの境界をまたいだときに、図４に示す様に位相に±２πする必要が出てくる。ここでさらに問題となるのがアンラッピングの“不確定性”である。位相変動が−πかπの境界をまたいだのかどうかの判断には不確定性がつきまとう。例えば０であった位相が次点で２π/３になったとき、２π/３分の位相の増加によって２π/３となったのか、−４π/３分の位相の減少で２π/３となったのかはわからない。−４π/３分の位相減少の場合は境界をまたいでいるため、アンラッピングが必要である。点間における位相変化幅の許容基準を設けるなどしてアンラッピングする又はしないかの判断をしなくてはならない。
図５に本手法の演算フローを示す。本手法では二つの信号ｆ１，ｆ２それぞれについて位相を求め、アンラッピングを行う。位相の微分を行ったのちに二つの値の差を計算することで群遅延を求める。しかし、上記の手法は計算スピードにおいては従来のＳＴＦＴ法に比べ、格段の改善がなされるのであるが、アンラッピングの手間・不確定性の問題を伴いリアルタイムの信号処理と動作安定性において十分に満足できるものではなかった。
【００１９】
そこで、本発明らが推奨する手法は「Ｚ変換を用いる方法」である。これは微分演算をＺ領域に持ち込む手法である。
このＺ変換を用いる方法について以下に説明する。位相微分を用いる方法で述べたアンラッピングの手間・不確定性を気にしなくていいのがＺ変換を用いる方法である。一般に計算コストが高い微分演算を行わなくてもよいこともこの手法の特徴である。この手法では微分演算をＺ領域に持ち込むことで巧みに回避している。数式的な導出を以下に示す。
周波数応答Ｈを振幅Ａと位相πで表すと
【数４】
2015105909.tif 000007
となるため、位相は
【数５】
2015105909.tif 000008
と表される。このとき群遅延は
【数６】
2015105909.tif 000009
となるが、実関数であるために
【数７】
2015105909.tif 000010
となる。この計算をＺ領域で行うために、サンプル波数間隔をＴとして
【数８】
2015105909.tif 000011
を用いると群遅延は
【数９】
2015105909.tif 000012
と表される。Ｚ＝ｅｊ・Ｔより
【数１０】
2015105909.tif 000013
となる。ここでＺ変換をＺ[ ]で表すと
【数１１】
2015105909.tif 000014
となるので、Ｚ領域での群遅延は
【数１２】
2015105909.tif 000015
で表され、これを周波数領域にすると、フーリエ変換をfft[ ]で表して
【数１３】
2015105909.tif 000016
となる。
式(13)より、微分およびアンラッピングを行わずに群遅延を計算できる。しかしフーリエ変換や、周波数応答Ｈからインパルス応答ｈを算出するときに逆フーリエ変換を用いるため、演算過程でノイズ除去を行う必要が出てくる。
【００２０】
図６に上記手法の演算フローを示す。基本的には二つの信号ｆ１，ｆ２それぞれにフーリエ変換を施し、周波数応答を求め、逆フーリエ変換によりインパルス応答を算出する。インパルス応答を用いて群遅延を計算する。
ノイズ除去のために窓関数（ハン窓）を演算するタイミングが二度、ゼロパッディングを行うタイミングが一度ある。逆フーリエ変換によるノイズを除去するため、図６に示されるウィンドウ２のタイミングでインパルス応答ｈに行う演算は必須である。二つの干渉信号ｆ１，ｆ２それぞれをフーリエ変換する前に窓関数をかけるウィンドウ１の場合は、この演算の有無で群遅延計算結果の特性が変わってくる。ハン窓を用いるのは窓両端で振幅が０になるからである。本手法ではハン窓を採用した。
ゼロパッディングはｆ１信号の後ろに、ｆ２信号の前に同じ長さの０列を付け足す行程である。これにより周波数応答を算出するときにノイズが出にくくなる場合がある。しかしデータ量が二倍になるため計算コストもそれだけ高くなる。
【００２１】
前述した位相微分、Ｚ変換を用いた群遅延計算による信号処理手法の性能比較を実証データで表２に示す。Ｚ変換はウィンドウ１と２の演算を施している。信号処理性能はいずれの手法においてもＬ長には依らず、波長掃引幅（sweep幅）に依存していたため、波長掃引幅７ｎｍ，14ｎｍ，21ｎｍの３ケースを取り上げた。計測時間は同等の測定結果をもたらすＳＴＦＴ演算にかかる時間を１としている。
【表２】
2015105909.tif 000017
シミュレーションによる検討では、群遅延計算に基づく信号処理手法を採用することで、従来のＳＴＦＴによる手法よりも空間分解能、計測速度が大きく向上され、測定精度に関してもほぼ同様の性能を示すことがわかった。位相微分を用いる手法ではその性能は波長掃引幅に依存するような関係があることがわかった。Ｚ変換を用いた手法では特に空間分解がウィンドウのかけ方に影響されることがわかった。測定速度は位相微分を用いる手法よりも数倍速かった。
【００２２】
本発明者らは上記の信号処理手法に追加修正する形で更なる改良型のＺ変換を用いた信号処理方法を開発した。
この手法を以下に説明する。
前述したように、周波数応答Ｈ(ω)の位相をφ(ω)とすれば、群遅延Ｔｇ(ω)は式（3）で表される。これはインパルス応答ｈを用いれば
【数１４】
2015105909.tif 000018
と表せる。上式において、fft( )はフーリエ変換、Ｔは係数、Re( )は複素数関数の実部を示す。
図７を参照しながら信号処理の流れを説明する。
この手法ではインパルス応答ｈ(ｋ)は畳み込み積分により求める。
ｈ(ｋ)＝(ｆ２*ｆ１)(ｋ)            （15）
フーリエ変換後の関数を高速フーリエ変換fft ( )を用いて
Ｈ(ω)＝fft (ｈ(ｋ))               （16）
Ｈｋ(ω)＝fft (ｋ ｈ(ｋ))          （17）
とすれば、求める群遅延は
【数１５】
2015105909.tif 000019
となる。群遅延シフトをブラッグ波長シフトあるいはひずみ変動に変換するよう係数Ｔを設定し、角周波数をＦＢＧ位置に対応させることでひずみ分布を得る。
ＨｋやＨはｈｋやｈのフーリエ変換なので、
【数１６】
2015105909.tif 000020
とも書くことができる。
【数１７】
2015105909.tif 000021
となる。これは、ノイズ処理のない場合ではインパルス応答ｈを用いて行ったものである。
【００２３】
ノイズ処理を行う場合には図８に示すように、上記の信号処理手法に追加修正する形でノイズ処理手法となる。
インパルス応答ｈ(ｋ)およびｋ ｈ(ｋ)に窓関数をかけて特定の区間における次の関数を取り出す。
【数１８】
2015105909.tif 000022
上につくハットマークは、一部分であることを意味する。
その一部分について群遅延を求めると
【数１９】
2015105909.tif 000023
となる。ちなみに、抜きだす作業は窓関数をかけることで行うが、これは窓を通るもの以外は０にする作業である。
図９に示すように、この抜き出す位置を変えて（窓関数をずらして）
【数２０】
2015105909.tif 000024
を何度も計算し、それらの総和で群遅延を求めると
【数２１】
2015105909.tif 000025
が得られる。何度も計算した分母と分子それぞれの総和を割ることは、“重みづけして平均化する”ことに相当するので、この過程でノイズは抑えられることとなる。
この処理はＦＦＴ長の縮小により算出されるひずみ点数が減少する、すなわち空間分解能が低下する反面、ノイズを大きく低減する効果がある。空間分解能、ノイズ低減効果、信号処理時間は窓関数とそのスライドの長さにより調整される。この作業はＳＴＦＴ演算中の計算過程に似ているが、窓関数の長さ、窓関数をスライドする長さともにＳＴＦＴの場合よりはるかに大きく設定しても十分な効果が得られることが大きなメリットである。
【００２４】
本発明に係る改良型のＺ変換による群遅延演算を用いたＯＦＤＲ方式光ファイバ計測装置の効果をシミュレーションにより検証したデータを以下に示す。
この試験では図１３に示されるＯＦＤＲシステムを使用し、センシング領域を５ｍ（これはリファレンスアーム長Ｌ＝１０ｍに相当）とし、波長可変光源波長掃引を１０ｎｍとした。そして、群遅延演算におけるノイズ処理のための窓関数の幅を1.35ｎｍ、窓関数のスライド長さを0.7ｎｍとした。ＦＢＧ観測信号を生成したシミュレーション条件の詳細は表３に示す通りである。
【表３】
2015105909.tif 000026
生成された信号を従来のＳＴＦＴ、提案する改良型群遅延演算を用いて信号処理した結果を表４に示す。なお、計算ソフトには Matlab を用いた。
【表４】
2015105909.tif 000027
二つの信号処理手法を比較すると、出力されるブラッグ波長分布の波長測定精度、空間分解能は同等を維持して計算効率が２８倍であることを示している。
本発明の群遅延演算は信号処理速度において従来のＳＴＦＴの数十倍の性能が得られることが確認できた。
【符号の説明】
【００２５】
Ｄ1,Ｄ2,Ｄ3,Ｄ4：  光強度ディテクタ
Ｃ1：  ５％−９５％広帯域カップラ
Ｃ2-Ｃ5：  ５０％−５０％広帯域カップラ
Ｒ1-Ｒ4：  全反射が起きるような終端
【図１】
2015105909.tif 000028
【図２】
2015105909.tif 000029
【図３】
2015105909.tif 000030
【図４】
2015105909.tif 000031
【図５】
2015105909.tif 000032
【図６】
2015105909.tif 000033
【図７】
2015105909.tif 000034
【図８】
2015105909.tif 000035
【図９】
2015105909.tif 000036
【図１０】
2015105909.tif 000037
【図１１】
2015105909.tif 000038
【図１２】
2015105909.tif 000039
【図１３】
2015105909.tif 000040