(57)【要約】【課題】  簡易な構成で測定精度を低下させずに時間応答を速く中性子を測定する。【解決手段】  中性子測定装置１は、中性子照射によってアナログ信号を発生させる中性子検出器１１と、アナログ信号を所定のサンプリング周期でサンプリングし、デジタル信号に変換するアナログデジタルコンバータ２と、デジタル信号を受信し、二乗演算データとして二乗演算する二乗演算部３と、Ｎ回サンプリング分の二乗演算データを受信しＮ回平均データとして平均する平均データ演算部４と、Ｎ回平均データを複数回分保持し、移動平均して出力する移動平均フィルタ５と、複数の平均データの前後の信号変化率を演算する変化率演算部６と、変化率演算部で演算された信号変化率が所定の値を上回ったときに、アナログデジタルコンバータ２のサンプリング周波数を高く設定する周波数変更部７とを備える。【選択図】  図１
2013200259.tif 000002
【特許請求の範囲】
【請求項１】
中性子照射によってアナログ信号を発生させる中性子検出器と、
この中性子検出器からのアナログ信号を所定のサンプリング周期でサンプリングし、デジタル信号に変換するアナログデジタルコンバータと、
このアナログデジタルコンバータからデジタル信号を受信し、二乗演算データとして二乗演算する二乗演算部と、
Ｎ回サンプリング分の前記二乗演算データを受信しＮ回平均データとして平均する平均データ演算部と、
このＮ回平均データを複数回分保持し、移動平均して出力する移動平均フィルタと、
複数の前記平均データの前後の信号変化率を演算する変化率演算部と、
この変化率演算部で演算された前記信号変化率を所定の上限閾値と比較し、この所定の上限閾値を前記信号変化率が上回ったときに、前記アナログデジタルコンバータの前記サンプリング周波数を高く設定する周波数変更部とを備えることを特徴とする中性子測定装置。
【請求項２】
前記周波数変更部は、前記信号変化率と前記所定の上限閾値との比較に代えて、前記信号変化率に対する適切な前記サンプリング周波数をあらかじめテーブル化して記憶し、受信した前記信号変化率に対応する前記サンプリング周波数を前記アナログデジタルコンバータの前記サンプリング周波数として設定することを特徴とする請求項１に記載の中性子測定装置。
【請求項３】
前記周波数変更部は、前記信号変化率が所定の下限閾値を下回ったときに、前記アナログデジタルコンバータの前記サンプリング周波数を低く設定することを特徴とする請求項１または請求項２に記載の中性子測定装置。
【請求項４】
初期周波数、初期平均データ数、初期移動平均数をそれぞれ前記アナログデジタルコンバータ、前記平均データ演算部、前記移動平均フィルタに送信して中性子測定開始時の初期値の設定を行う初期値設定部をさらに備えることを特徴とする請求項１乃至請求項３の何れか一項に記載の中性子測定装置。
【請求項５】
中性子の測定対象となる原子炉の種類を含む中性子測定環境に応じてあらかじめ求めた前記初期周波数、前記初期平均データ数、前記初期移動平均数を記憶し、前記初期値設定部に送信可能に接続される初期値データベースをさらに備えることを特徴とする請求項４に記載の中性子測定装置。
【請求項６】
中性子照射によってアナログ信号を発生させ、所定のサンプリング周期でサンプリングし、デジタル信号に変換するデジタル変換工程と、
前記デジタル信号を二乗演算データとして二乗演算する二乗演算工程と、
Ｎ回サンプリング分の前記二乗演算データをＮ回平均データとして平均する平均データ演算工程と、
このＮ回平均データを複数回分保持し、移動平均して出力する移動平均工程と、
複数の前記平均データの前後の信号変化率を演算する信号変化率演算工程と、
この信号変化率が所定の値を上回ったときに、前記デジタル変換工程における前記サンプリング周波数を高く設定する周波数変更工程とを備えることを特徴とする中性子測定方法。 
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、中性子を測定する中性子測定装置およびその中性子測定方法に関する。
【背景技術】
【０００２】
原子力発電所において原子炉で発生する中性子を測定する中性子測定装置には、中性子照射によって中性子検出器から信号を発生させ、原子炉出力が低い順にパルスモード、ＭＳＶ（Ｍｅａｎ  Ｓｑｕａｒｅ  Ｖａｌｕｅ：二乗平均値）モード、電流モードによって信号を測定するものが開発されている（例えば、特許文献１参照。）。
【０００３】
このうちＭＳＶモードでは、中性子検出器からの個々のパルス信号が重なり合いパルス計測ができないため、重なり合った信号の揺らぎ（脈動）の幅の二乗平均値が原子炉出力に比例する特性、いわゆるキャンベル特性を利用して、中性子検出器からの信号を二乗演算して中性子を測定する。
【０００４】
核融合実験施設においても、プラズマから放射される中性子発生量を測定するために同様の中性子測定装置が設けられ、ＭＳＶモードにおいて中性子ビーム入射（ＮＢＩ：Ｎｅｕｔｒａｌ  Ｂｅａｍ  Ｉｎｊｅｃｔｉｏｎ）時のプラズマ生成から消滅までの中性子発生量の変化を測定する。ここで、中性子ビーム入射における中性子発生量は急激な上昇を示すため、中性子測定装置には特に速い時間応答性（１〜１０ｍｓ程度）が要求される。
【０００５】
そこで、ＭＳＶモードにおいて中性子照射による信号変化を監視し、信号の過渡変化時には時定数を短くして時間応答性を高くするデジタルフィルタが開発されている（例えば、特許文献２参照。）。また、異なる仮の時定数で計数率を測定し、測定値の差に基づいて適切な時定数を決定する放射線測定装置が開発されている（例えば、特許文献３参照。）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開昭６４−８８１７９号公報
【特許文献２】特公平６−８３０１８号公報
【特許文献３】特開２００９−２７０８５０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ここで中性子測定装置は、中性子検出器からのアナログ信号を所定のサンプリング周期でサンプリングし、複数サンプリング分のデータを移動平均して出力している。そして中性子測定装置の時間応答性は、検出器からの信号をサンプリングし、移動平均等を行って出力するまでの時間を示す時定数によって決定される。しかしながら、上述した特許文献２に記載のデジタルフィルタにおいて、サンプリング周波数を一定として、移動平均数を少なくすることによって時定数を短くすると、サンプル数が不足して、信号の揺らぎやノイズによって統計的な測定精度が低下するおそれがある。また、上述した特許文献３に記載の技術は、異なる複数の時定数で計数率を測定しているため、設備が冗長化し、コストが増大するおそれがある。
【０００８】
そこで本発明では、簡易な構成で測定精度を低下させずに時間応答を速く中性子を測定することができる中性子測定装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するために、本発明の中性子測定装置は、中性子照射によってアナログ信号を発生させる中性子検出器と、中性子検出器からのアナログ信号を所定のサンプリング周期でサンプリングし、デジタル信号に変換するアナログデジタルコンバータと、アナログデジタルコンバータからデジタル信号を受信し、二乗演算データとして二乗演算する二乗演算部と、Ｎ回サンプリング分の二乗演算データを受信しＮ回平均データとして平均する平均データ演算部と、Ｎ回平均データを複数回分保持し、移動平均して出力する移動平均フィルタと、複数の平均データの前後の信号変化率を演算する変化率演算部と、変化率演算部で演算された信号変化率が所定の値を上回ったときに、アナログデジタルコンバータのサンプリング周波数を高く設定する周波数変更部とを備えることを特徴とする。
【００１０】
さらに上記目的を達成するために、本発明の中性子測定方法は、中性子照射によってアナログ信号を発生させ、所定のサンプリング周期でサンプリングし、デジタル信号に変換するデジタル変換工程と、デジタル信号を二乗演算データとして二乗演算する二乗演算工程と、Ｎ回サンプリング分の二乗演算データをＮ回平均データとして平均する平均データ演算工程と、Ｎ回平均データを複数回分保持し、移動平均して出力する移動平均工程と、複数の平均データの前後の信号変化率を演算する信号変化率演算工程と、信号変化率が所定の値を上回ったときに、デジタル変換工程におけるサンプリング周波数を高く設定する周波数変更工程とを備えることを特徴とする。
【発明の効果】
【００１１】
本発明によれば、簡易な構成で測定精度を低下させずに時間応答を速く中性子を測定することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の第１の実施形態に係る中性子測定装置の概略構成図。
【図２】本発明の第１の実施形態に係る中性子測定装置の模擬入力データに対する応答を示すグラフ。
【図３】核融合実験施設における中性子ビーム入射による中性子発生量の例を示すグラフ。
【図４】本発明の第２の実施形態に係る中性子測定装置の概略構成図。
【図５】ステップ信号に対する時定数別の応答シミュレーションを示すグラフ。
【発明を実施するための形態】
【００１３】
以下、本発明の実施形態を説明する。
【００１４】
（第１の実施形態）
（構成）
以下、本発明の第１の実施形態に係る中性子測定装置について図１乃至図３を参照して説明する。図１は、本発明の第１の実施形態に係る中性子測定装置の概略構成図である。
【００１５】
中性子測定装置１は、アナログデジタル（ＡＤ）コンバータ２と、二乗演算部３と、平均データ演算部４と、移動平均フィルタ５と、変化率演算部６と、周波数変更部７と、検出器１１と、表示装置１２とから構成される。また、図１に示すように二乗演算部３と、平均データ演算部４と、移動平均フィルタ５と、変化率演算部６と、周波数変更部７は、デジタル演算処理部としてＡＳＩＣやＦＰＧＡで構成することができる。
【００１６】
検出器１１は、中性子照射によって内部の電離物質を電離させてアナログ信号２１を発生させることができる装置である。検出器１１は、ＡＤコンバータ２にアナログ信号２１を送信することができるように接続される。
【００１７】
ＡＤコンバータ２は、アナログ信号２１を所定のサンプリング周波数ｆｓでサンプリングし、デジタル信号２２にデジタル変換する。ＡＤコンバータ２は、二乗演算部３にデジタル信号２２を送信することができるように接続される。
【００１８】
二乗演算部３は、デジタル信号２２を二乗演算して二乗演算データ２３に変換する。二乗演算部３は、二乗演算データ２３を平均データ演算部４に送信することができるように接続される。平均データ演算部４は、二乗演算データ２３をＮ回サンプリング分受信し平均化してＮ回平均データ２４を算出する。平均データ演算部４は、Ｎ回平均データ２４を移動平均フィルタ５に送信することができるように接続される。さらに平均データ演算部４は、受信し保持しているＮ回平均データ２４を今回平均データ２４ａとして変化率演算部６に送信することができるように接続される。
【００１９】
移動平均フィルタ５は、Ｎ回平均データ２４をＭ回平均分受信し移動平均して移動平均出力２７として出力する。移動平均フィルタ５は、移動平均出力２７を表示装置１２に送信することができるように接続される。さらに移動平均フィルタ５は、Ｎ回平均データ２４を前回平均データ２４ｂとして変化率演算部６に送信することができるように接続される。
【００２０】
変化率演算部６は、今回平均データ２４ａと前回平均データ２４ｂとから信号変化率を演算する。変化率演算部６は、演算した信号変化率を信号変化率情報２５として周波数変更部７に送信することができるように接続される。周波数変更部７は、受信した信号変化率情報２５に基づいて周波数変更指令２６をＡＤコンバータ２に送信することができるように接続される。
【００２１】
（作用）
以下、本発明の第１の実施形態の作用について説明する。中性子測定装置１の作用は、中性子照射によって検出器１１で発生した信号を測定し出力する作用と、信号過渡変化を監視してサンプリング周波数を変化させる作用に大別される。
【００２２】
まず、中性子照射による信号を測定する作用について説明する。検出器１１は、中性子照射によってアナログ信号２１を発生させＡＤコンバータ２に送信する。ＡＤコンバータ２は、アナログ信号２１を所定のサンプリング周波数ｆｓでサンプリングしてデジタル信号２２にデジタル変換する。サンプリング周波数ｆｓは、サンプリングするアナログ信号２１の周波数の５倍ないし１０倍以上とすることが望ましい。ＡＤコンバータ２は、サンプリング周波数ｆｓの１周期毎にアナログ信号２１をサンプリングしてデジタル信号２２に変換し、各サンプリング毎にデジタル信号２２を二乗演算部３に送信する。二乗演算部３は、デジタル信号２２を二乗演算して二乗演算データ２３に変換する。二乗演算部３は、二乗演算データ２３を平均データ演算部４に順次送信する。
【００２３】
平均データ演算部４は、二乗演算データ２３をＮ回サンプリング分受信し平均化してＮ回平均データ２４を算出する。なお、平均データ演算部４における平均化サンプリング数Ｎは測定対象となる検出器１１からのアナログ信号２１の一周期分以上となるように設定することが望ましい。
【００２４】
平均データ演算部４は、Ｎ回平均データ２４を移動平均フィルタ５に順次送信する。移動平均フィルタ５は、Ｎ回平均データ２４を受信して、少なくとも直近のＭ回分保持する。移動平均フィルタ５は、Ｎ回平均データ２４を順次受信するたびに直近のＭ回分のＮ回平均データ２４を移動平均して移動平均出力２７として表示装置１２に出力する。
【００２５】
移動平均フィルタ５は、平均データ演算部４からＭ＋１回目のＮ回平均データ２４を受信すると、２回目からＭ＋１回目までのＮ回平均データ２４を移動平均して移動平均出力２７を出力する。このとき、１回目のＮ回平均データ２４は移動平均に不要であるため消去することによって、Ｎ回平均データ２４を最小限のＭ回分保持することで移動平均フィルタ５内の記憶データを最小限にすることができる。
【００２６】
表示装置１２は、移動平均出力２７を受信して中性子測定結果として表示する。なお、移動平均フィルタ５には、表示装置１２に代えて中性子測定結果を印字する印字装置や中性子測定結果を閾値と比較して警報を発する警報装置を接続する構成としてもよい。
【００２７】
次に、信号過渡変化を監視してサンプリング周波数ｆｓを変化させる作用について説明する。平均データ演算部４は、Ｎ回平均データ２４を今回平均データ２４ａとして変化率演算部６に送信する。また移動平均フィルタ５は、Ｍ−１回目のＮ回平均データ２４を前回平均データ２４ｂとして変化率演算部６に送信する。
【００２８】
変化率演算部６は、今回平均データ２４ａを前回平均データ２４ｂによって除して信号変化率を演算する。変化率演算部６は、演算した信号変化率を信号変化率情報２５として周波数変更部７に送信する。周波数変更部７は、受信した信号変化率情報２５が示す信号変化率を所定の上限閾値と比較する。信号変化率が所定の閾値の閾値を上回った場合、周波数変更部７はＡＤコンバータ２のサンプリング周波数ｆｓを高く設定させる周波数変更指令２６をＡＤコンバータ２に送信する。ＡＤコンバータ２は、周波数変更指令２６を受信するとサンプリング周波数ｆｓを高く設定してサンプリングを継続する。
【００２９】
なお、信号変化率が所定の閾値を上回るときにサンプリング周波数ｆｓを高く設定するだけでなく、信号変化率に対する適切なサンプリング周波数ｆｓをあらかじめ求めてテーブル化して周波数変更部７において記憶し、周波数変更部７が受信した信号変化率に対応するサンプリング周波数ｆｓをテーブルから読み出して周波数変更指令２６として送信する構成とすることもできる。
【００３０】
以下、サンプリング周波数ｆｓと時定数τの関係について説明する。（式１）に移動平均時間Ｔ、平均化サンプリング数Ｎ、サンプリング周波数ｆｓ、移動平均数Ｍとの関係を示す。
【００３１】
Ｔ＝Ｎ・Ｍ／ｆｓ［ｓ］・・・（１）
移動平均フィルタ５における時定数τは、移動平均時間／２と同じとみなせるので、時定数τと移動平均時間Ｔは（式２）の関係がある。
【００３２】
τ≒Ｔ／２＝Ｎ・Ｍ／２ｆｓ［ｓ］・・・（２）
（式２）より、時定数τはサンプリング周波数ｆｓに反比例する関係にある。したがって、サンプリング周波数ｆｓを高くすることで、時定数τを小さくすることができる。また平均データ化数Ｎおよび移動平均数Ｍは固定としておくことで、時定数τを小さくしながらサンプル数を一定にして測定精度を保つことができる。
【００３３】
ここで、時定数と移動平均の関係について述べる。図５は、ステップ信号に対する時定数別の応答シミュレーションを示すグラフである。τ＝２ｍｓ、２５ｍｓ、Ｔ＝４ｍｓ時のステップ応答シミュレーション結果を示す。幅１０ｍｓのステップ信号に対して、τ＝２５ｍｓのフィルタでは追従できず、ステップ信号の立ち上がりを再現できていないが、τ＝２ｍｓのフィルタでは十分追従できることがわかる。また、τ＝２ｍｓとＴ＝４ｍｓの結果を比較すると、ステップ信号の立上りから２ｍｓ後に０．６３で交差し、（式２）において、時定数τの時間が一致することがわかる。また、移動平均フィルタはτ以降の時間では１次遅れフィルタを使用するよりも速く追従できる特性がある。
【００３４】
図２は、本発明の第１の実施形態に係る中性子測定装置の模擬入力データに対する応答を示すグラフである。中性子測定装置１において、例えば平均データ演算部４における平均化サンプリング数Ｎを８、移動平均フィルタ５における移動平均数Ｍを４０に固定する。
【００３５】
定常状態のサンプリング周波数ｆｓを１ｋＨｚとすると、時定数は１６ｍｓとなる。さらに移動平均数Ｍを一定保ちながらサンプリング周波数ｆｓを高くして、時定数を８ｍｓ、４ｍｓ、２ｍｓと順に小さくする。これらにＳＮ比を０．１としたランダムノイズを重畳した模擬入力データを入力したシミュレーション結果を示す。
【００３６】
定常状態において時定数τが１６ｍｓであるとき、図２に示すように信号の過渡状態では時間応答が遅く追従性が悪い。そこで時定数を８ｍｓ、４ｍｓ、２ｍｓと順に小さくすると時間応答が速くなり追従性が改善され、時定数が４ｍｓで十分な追従性が得られる。
【００３７】
信号変化率の演算と信号変化率に基づくサンプリング周波数ｆｓの変更は平均データ演算部４におけるＮ回平均データ２４の演算ごとに行われるため、１０ｍｓ以内に信号の過渡状態に追従できる。
【００３８】
図３は、核融合実験施設における中性子ビーム入射による中性子発生量の例を示すグラフである。核融合実験施設における中性粒子ビーム入射（ＮＢＩ：Ｎｅｕｔｒａｌ  Ｂｅａｍ  Ｉｎｊｅｃｔｉｏｎ）時のプラズマ生成から消滅までには、図３に示すような中性子の急激な立上りと減衰が生じる。中性子測定装置１は、この中性子の立上りと減衰を信号変化率によって監視し、信号変化率に応じてサンプリング周波数ｆｓを変化させることによって、中性子の増減に追従して測定を行うことができる。
【００３９】
（効果）
本発明の第１の実施形態によれば、前後のＮ回平均データ２４の信号変化率に応じてサンプリング周波数ｆｓを変化させることにより、簡易な構成で測定精度を低下させずに時間応答を高め、中性子が増減する過渡状態を測定することができる。
【００４０】
また、検出器１１の信号の信号変化率演算において、ＡＤコンバータ２におけるデジタル信号２２の前後を比較して信号変化率を演算すると、デジタル信号２２のノイズや揺らぎが信号変化率にも反映され、安定してサンプリング周波数を変更することができない。しかしながら本実施形態では、デジタル信号２２を二乗演算してＮ回分を平均化しノイズや揺らぎの影響を抑えたＮ回平均データ２４によって信号変化率を演算することによって、信号変化率への揺らぎやノイズの影響を抑え、安定してサンプリング周波数を変更することができる。
【００４１】
なお、周波数変更部７においてさらに所定の下限閾値を設定し、信号変化率が所定の下限閾値を下回ったときにはサンプリング周波数ｆｓを低く設定する周波数変更指令２６を送信することもできる。この場合、サンプリング周波数ｆｓを下げることによって時定数を下げ、単位時間当たりの移動平均出力２７の出力数を下げることによって演算の負荷を防ぐことができる。さらに移動平均出力２７を外部記憶装置によって記憶する場合、必要最小限の時間応答を満たす移動平均出力２７を記憶することができ、記録データの増大を防ぐことができる。
【００４２】
（第２の実施形態）
（構成）
以下、本発明の第２の実施形態に係る中性子測定装置について図４を参照して説明する。第１の実施形態に係る中性子測定装置の各部と同一部分には同一符号を付し、同一の構成についての説明は省略する。
【００４３】
図４は、本発明の第２の実施形態に係る中性子測定装置を示す概略構成図である。第２の実施形態が第１の実施形態と異なる点は、初期値入力部１３および初期値データベース（ＤＢ）１４をさらに設けた点である。初期値入力部１３は、ＡＤコンバータ２、平均データ演算部４、移動平均フィルタ５にそれぞれ初期周波数３２、初期平均データ数３３、初期移動平均数３４を送信することができるように接続される。初期値ＤＢ１４は、初期値情報３１を初期値入力部１３に送信することができるように接続される。
【００４４】
（作用）
以下、本発明の第２の実施形態の作用について説明する。初期値入力部１３は、キーボードやマウスといった入力インターフェイスによって、平均化サンプリング数Ｎ、移動平均データ数Ｍ、サンプリング周波数ｆｓの初期値を入力して時定数τを設定することができる。また、上述した式（２）より平均化サンプリング数Ｎ、移動平均データ数Ｍを固定として時定数τのみを入力してサンプリング周波数ｆｓを決定する構成としてもよい。
【００４５】
初期値入力部１３は、これら初期値をそれぞれ初期周波数３２、初期平均データ数３３、初期移動平均数３４としてＡＤコンバータ２、平均データ演算部４、移動平均フィルタ５に送信して初期値の設定を行う。測定開始後は、Ｍは固定され、上述した信号変化率に応じたサンプリング周波数ｆｓの変更が行われる。
【００４６】
さらに、初期値ＤＢ１４内に中性子測定装置１を設置する原子炉の種類や出力、運転計画といった中性子測定環境において初期周波数３２、初期平均データ数３３、初期移動平均数３４をあらかじめ初期値情報３１として求めて記憶させ、初期値情報３１を初期値入力部１３に読み込ませることによって初期値を設定することもできる。
【００４７】
（効果）
本発明の第２の実施形態によれば、任意の初期値を設定できることにより、中性子測定環境に対応した中性子の測定をすることができる。
【００４８】
なお、本発明の実施形態は上述した実施形態に限られないことは言うまでもない。例えば、直近の前後のＮ回平均データ２４を用いて信号変化率を演算するだけでなく、前後のＮ回平均データ２４であれば、移動平均フィルタ５に保持されるＮ回平均データ２４のうち一回おきや、２回おきのＮ回平均データ２４を用いて信号変化率を演算することができる。さらに信号変化率の演算は平均データ演算部４における平均化演算毎全てにおいて行うだけでなく、所定のタイミングの平均化演算において行ってもよい。
【００４９】
また、信号変化率演算において、移動平均フィルタ５内に保持されたＮ回平均データ２４を前回平均データ２４ｂとして用いるだけでなく、変化率演算部６において前後２つのＮ回平均データ２４を平均データ演算部４から受信して保持し、それぞれ今回平均データ２４ａ、前回平均データ２４ｂとして信号変化率を演算してもよい。
【００５０】
さらに周波数変更部７において、信号変化率が閾値を上回ったときにサンプリング周波数ｆｓを高く設定するだけでなく、信号変化率について適切なサンプリング周波数ｆｓをあらかじめテーブル化し、各々の変化率についてサンプリング周波数ｆｓを設定するようにして、信号変化率が高くなるにつれ、サンプリング周波数を速く設定することもできる。 
【符号の説明】
【００５１】
１・・・中性子測定装置
２・・・ＡＤ（アナログ／デジタル）コンバータ
３・・・二乗演算部
４・・・平均データ演算部
５・・・移動平均フィルタ
６・・・変化率演算部
７・・・周波数変更部
１１・・・検出器
１２・・・表示装置
１３・・・初期値入力部
１４・・・初期値ＤＢ
２１・・・アナログ信号
２２・・・デジタル信号
２３・・・二乗演算データ
２４・・・平均データ
２４ａ・・・今回平均データ
２４ａ・・・前回平均データ
２５・・・信号変化率
２６・・・周波数変更指令
２７・・・移動平均出力
３１・・・初期値情報
３２・・・初期周波数
３３・・・初期平均データ数
３４・・・初期移動平均数 
【図１】
2013200259.tif 000003
【図２】
2013200259.tif 000004
【図３】
2013200259.tif 000005
【図４】
2013200259.tif 000006
【図５】
2013200259.tif 000007