移動平均値算出装置および移動平均値算出方法
【課題】入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値の出力を行うことができる移動平均値算出装置およびその移動平均値算出方法を提供する。
【解決手段】逐次入力される複数の入力値（T1〜Tn）の平均値を最新の入力値（Tn）から遡ってm個の入力値（但し、m≦n）を用いて、逐次算出する移動平均値算出部６は、最新の入力値（Tn）と過去の入力値（T1〜Tn-1のいずれか）との差（｜ΔTn｜）を算出し、該差（｜ΔTn｜）と、入力値に含まれ得るノイズの大きさ(N)から定められる閾値（Th）との大小関係により、前記平均値を算出する際に用いる入力値の個数（M）を設定する個数設定手段６２を備える。
【解決手段】逐次入力される複数の入力値（T1〜Tn）の平均値を最新の入力値（Tn）から遡ってm個の入力値（但し、m≦n）を用いて、逐次算出する移動平均値算出部６は、最新の入力値（Tn）と過去の入力値（T1〜Tn-1のいずれか）との差（｜ΔTn｜）を算出し、該差（｜ΔTn｜）と、入力値に含まれ得るノイズの大きさ(N)から定められる閾値（Th）との大小関係により、前記平均値を算出する際に用いる入力値の個数（M）を設定する個数設定手段６２を備える。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出装置および移動平均値算出方法に関する。
【背景技術】
【０００２】
従来、移動平均値算出装置としては、特許文献１に示すように、時間平均を算出する区間ｓを、ｓ（＝β／α）として時間平均値を逐次算出するものが知られている。ここで、βは、時間的に変化する変化量をデジタル化する際の最小分解能であり、αは、時間平均値のしきい値である。
【０００３】
かかる移動平均値算出装置によれば、例えば、温度検出値のデジタル変換値そのものではなく、その移動平均値に基づいて温度の変化率を求めるので、該変化率の落ち込みは生じ難くなり、また落ち込みが生じたとしてもその程度が小さくなる。そのため、変化率のデジタル的なばらつきの程度が小さくなり、変化率がしきい値αを超えているか否かの判定精度が向上する。
【０００４】
しかし、従来の移動平均値算出装置では、例えば、加熱炉のように急激な温度上昇・下降が起きた後に安定した状態が続くような測定環境では、急激な温度変動に追随できず、オーバーシュートしたかどうかの検知が困難であった。
【０００５】
これに対して、特許文献２に示す移動平均値算出装置では、圧力計における急激な圧力変化に対して、蓄積した検出値の平均化を行わず、検出値をそのまま出力する。この移動平均値算出装置によれば、急激な圧力変化が生じた場合にも、検出値がそのまま出力されるため、その変化に出力を追従させることができる。
【特許文献１】特開２０００−１２４８０６号公報
【特許文献２】特開平１１−１１８６５１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の移動平均値算出装置では、検出値の変動が、センサ故障や劣化に起因するノイズや一般的に生じる電気ノイズ等（以下、ノイズという）によるものである場合にも、これがそのまま出力されてしまうという問題がある。
【０００７】
これは、検出値を蓄積してその平均値を出力することにより、センサの出力がノイズによって急変することを回避して出力の安定化を図るという、移動平均を用いる本来の目的に反することとなる。
【０００８】
上記の事情に鑑みて、本発明は、入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値の出力を行うことができる移動平均値算出装置および移動平均値算出方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
第１発明の移動平均値算出装置は、逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出装置であって、前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定手段を備えることを特徴とする。
【００１０】
第１発明の移動平均値算出装置によれば、最新の入力値と過去の入力値との差に基づいて移動平均値を算出するために用いる個数を設定する。そのため、最新の入力値の変化が大きい場合には、例えば、入力値の個数を減少させて、この変化に平均値の出力を追従させ易くすることができる。一方、最新の入力値の変化が小さい場合には、例えば、入力値の個数を増加させて、ノイズの影響を低減させることにより安定した平均値の出力を行うことができる。
【００１１】
ここで、移動平均値を算出するために用いる個数は、前記差とノイズの大きさから定められる閾値との大小関係により設定される。そのため、例えば、閾値をノイズの大きさよりも大きくしておくことで、その閾値を超える入力値の変化にのみ平均値の出力を追従させ易くすることが可能となる。これにより、ノイズの影響を確実に低減して、入力値の変化に対応した適切な平均値の出力を行うことができる。
【００１２】
また、閾値は上述したノイズの大きさに基づいて定めるもののほか、ＳＮの比に基づき定めるもの、あるいは分解能／精度／再現性を考慮して定めるものなどが考えられる。いずれにせよ、ノイズの影響を低減して、入力値の変化に対応した適切な平均値の出力を行えるような値を閾値として設定すればよい。
【００１３】
なお、差を算出する際に用いる過去の入力値は、最新の入力値の直前（１周期前）の入力値でもよいが、該直前の入力値以外の値であってもよい。例えば、直前の入力値の更に直前の入力値（最新の入力値の２周期前に入力された入力値）や更にその直前の入力値（最新の入力値の３周期前に入力された入力値）を用いてもよい。これにより、最新の入力値の直前（１周期前）の入力値がノイズにより変化している場合にも、一定の間隔を空けることによりその影響を相対的に低減することができる。
【００１４】
第２発明の移動平均値算出装置は、第１発明の移動平均値算出装置において、前記個数設定手段は、前記所定の閾値である第１閾値に対して、前記差が該第１閾値より大きい場合には、前記入力値の個数を第１個数に設定し、該第１閾値以下である場合には、該個数を第１個数よりも多い第２個数に設定することを特徴とする。
【００１５】
第２発明の移動平均値算出装置によれば、第１閾値を、例えばノイズの大きさよりも大きい値に設定しておくことで、前記差が、第１閾値を超える入力値の変化である場合に、移動平均値を算出するために用いる個数が少ない第１個数で平均値を算出する。そのため、ノイズ成分に平均値の出力を追従させることなく、入力値の変化に平均値の出力を追従させることができる。一方、前記差が、第１閾値を以下の入力値の変化である場合に、該入力値の変化はノイズ成分に起因するものが含まれ得る。この場合、移動平均値を算出するために用いる個数が第１個数に比して多い第２個数で平均値を算出することにより、ノイズの影響を低減させることにより安定した平均値の出力を行うことができる。これにより、入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値の出力を行うことができる。
【００１６】
第３発明の移動平均値算出装置は、第２発明の移動平均値算出装置において、前記個数設定手段は、前記第１閾値よりも大きな値に設定された第２閾値に対して、前記差が該第２閾値より大きい場合には、前記第１個数を１とすることを特徴とする。
【００１７】
第３発明の移動平均値算出装置によれば、第１閾値よりも大きい第２閾値を設定し、前記差が、第２閾値を超える入力値の変化である場合に、移動平均値を算出するために用いる個数を１として、入力値をそのまま出力する。これにより、特に入力値の変化が著しい場合に、ノイズ成分に出力が追従することを回避しつつ、入力値の変化に出力を確実に追従させることができる。
【００１８】
第４発明の移動平均値算出装置は、第１〜第３発明のいずれかの移動平均値算出装置において、前記差は、前記最新の入力値と過去の入力値の平均値との差であることを特徴とする。
【００１９】
第４発明の移動平均値算出装置によれば、差を算出する際に用いる過去の入力値として過去の入力値の平均値を用いることで、ノイズにより個々変化し得る過去の入力値に依らず、差を算出することができる。そのため、ノイズの影響を低減して、入力値の変化を確実に捉え、これに対応した適切な平均値の出力を行うことができる。
【００２０】
第５発明の蛍光温度センサは、蛍光材料に投光する発光素子と、前記蛍光材料が発する蛍光を受光する受光素子と、前記受光素子の出力信号から温度検出値を生成する信号処理部とを含む蛍光温度センサにおいて、前記温度検出値を入力としてその移動平均値を出力する移動平均値算出部が請求項１乃至４のうちいずれか１項記載の移動平均値算出装置により構成されてなることを特徴とする。
【００２１】
第５発明によれば、蛍光温度センサは、一定の処理周期で取得された温度検出値に対して、移動平均値を算出して出力するため、第１〜第４発明の移動平均値算出装置を適用する対象として好適である。
【００２２】
第６発明の移動平均値算出方法は、逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出方法であって、前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定処理と、前記個数設定処理により設定された個数の入力値を用いて、これらの平均値を算出する平均値算出処理とを逐次実行することを特徴とする。
【００２３】
第６発明の移動平均値算出方法によれば、最新の入力値と過去の入力値との差に基づいて移動平均値を算出するために用いる個数を設定し、設定した個数の入力値を用いて平均値を逐次算出する。そのため、最新の入力値の変化が大きい場合には、入力値の個数を減少させて、この変化に平均値の出力を追従させ易くすることや、最新の入力値の変化が小さい場合には、入力値の個数を増加させて、ノイズの影響を低減させることにより安定した出力を行うことができる。
【００２４】
さらに、移動平均値を算出するために用いる個数は、前記差とノイズの大きさから定められる閾値との大小関係により設定されるため、閾値を例えばノイズの大きさよりも大きくしておくことで、その閾値を超える入力値の変化にのみ平均値の出力を追従させ易くることが可能となる。なお、閾値について種々の設定方法が考えられることは前述したとおりである。
【００２５】
このように、入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値を算出して出力することができる。
【発明を実施するための最良の形態】
【００２６】
本発明の一実施形態として、蛍光温度センサの移動平均値算出部に移動平均値算出装置を用いた例について、図１〜図３を参照して説明する。
【００２７】
まず、図１を参照して、本実施形態の蛍光温度センサの全体的な構成について説明する。蛍光温度センサは、温度によって異なる蛍光特性を示す蛍光材料１と、蛍光材料１に投光するＬＥＤ２と、ＬＥＤ２を駆動する電気信号を出力する駆動部３と、蛍光材料１が発する蛍光を受光するフォトダイオード４と、フォトダイオード４からの出力信号から対応する温度検出値を算出する信号処理部５と、信号処理部５に接続されて温度検出値の移動平均値を算出し、温度を表すデータとして出力する移動平均値算出部６（本発明の移動平均値算出装置に相当する）とを備える。
【００２８】
また、蛍光温度センサは、蛍光材料１への投光および蛍光材料１の蛍光の受光を行う光ファイバ７と備える。光ファイバ７は他端側が分岐して、ＬＥＤ２からの光を蛍光材料１に伝達する投光用光ファイバ７ａと、蛍光材料１の蛍光をフォトダイオード４に伝達する受光用光ファイバ７ｂとなっている。
【００２９】
蛍光材料１は、光ファイバ７の一端部を覆うように設けられた温度計測部１ａの中に、光ファイバ７のコア部に対向するように配置される。
【００３０】
ＬＥＤ２は、ＬＥＤモジュール２ａ内に配置された、例えば青色系の波長を発光色とする発光ダイオードである。ＬＥＤモジュール２ａは、投光用光ファイバ７ａが接続されるコネクタ部２ｂを有し、コネクタ部２ｂを介して接続された投光用光ファイバ７ａがＬＥＤ２の発光部と対向している。
【００３１】
駆動部３は、ＬＥＤ２の発光に必要な駆動電流の大きさおよび発光時間を規定したパルス電流を一定の処理周期でＬＥＤ２に印加する制御回路を備える。これにより、駆動部３は、例えば、蛍光材料１に対応して、一回の計測におけるＬＥＤ２の発光時間を１ｍｓ〜５００ｍｓの間のいずれかの時間とする所定の大きさのパルス電流をＬＥＤ２に印加する。
【００３２】
フォトダイオード４は、フォトダイオードモジュール４ａ内に配置されて、照射された光の光量（輝度）を測定する。フォトダイオードモジュール４ａは、受光用光ファイバ７ｂに接続されるコネクタ部４ｂを有し、コネクタ部４ｂを介して接続された受光用光ファイバ７ｂがフォトダイオード４の受光部と対向している。
【００３３】
信号処理部５は、前記処理周期（ＬＥＤ２の投光処理周期に対応した処理周期）で、フォトダイオード４によって測定された蛍光材料１の蛍光の減衰特性、特に蛍光緩和時間を計測する。具体的には、信号処理部５は、これがあらかじめ備える蛍光緩和時間と蛍光材料１との関係式（データテーブルやマップ等を含む）から、蛍光材料１が存在する温度測定環境の温度を算出して出力する。
【００３４】
移動平均値算出部６は、記憶保持手段６１と、個数設定手段６２と、平均値算出手段６３とを備え、信号処理部５から逐次入力される温度検出値の平均値を前記処理周期で逐次算出して出力する。
【００３５】
記憶保持手段６１は、図示しない内部メモリ等により構成され、信号処理部５から逐次入力された温度検出値を最新のものから一定の個数だけ記憶保持する。
【００３６】
個数設定手段６２は、記憶保持手段６１に記憶保持された温度検出値を読出し、最新の温度検出値と過去の温度検出値との差を算出し、該差と予め設定された閾値との大小関係により、平均値を算出する際に用いる温度検出値の個数を設定する。ここで、閾値は、温度検出値に含まれ得るノイズの大きさから定められる（図２のＳＴＥＰ２１および図３のＳＴＥＰ２２参照）。
【００３７】
平均値算出手段６３は、最新の温度検出値を含めて、個数設定手段６２により設定された個数の温度検出値を用いて、これらの平均値を算出する処理を逐次実行する。
【００３８】
なお、本実施形態において、各処理部３，５，６は、マイクロコンピュータ等により構成されたコントローラ１０に備えられてなる。コントローラ１０を構成する各処理部３，５，６は、ＣＰＵ、ＲＯＭ、ＲＡＭ等のハードウェアにより構成され、これらの各処理部３，５，６が共通のハードウェアによって構成されていてもよく、これらの各処理部３，５，６の一部又は全部が異なるハードウェアによって構成されていてもよい。
【００３９】
次に、図２に示すフローチャートを参照して、移動平均値算出部６により移動平均値を算出する処理について説明する。
【００４０】
まず、移動平均値算出部６は、前記処理周期毎に信号処理部５から出力される最新の温度検出値（Tn）を逐次取得する（ＳＴＥＰ１０）。ここで、温度検出値（T）の添え字n（n：自然数）は、信号処理部５から取得した温度検出値の順序を表し、nが小さいほど処理周期が過去のものを表す。例えば、１処理周期前の温度検出値は、Tn-1となる。
【００４１】
そして、移動平均値算出部６は、取得した温度検出値（Tn）を新たに記憶保持手段６１に記憶保持する（ＳＴＥＰ１１）。
【００４２】
次いで、移動平均値算出部６の個数設定手段６２は、記憶保持手段６１に記憶保持されている温度検出値（T1〜Tn）を読出し、最新の温度検出値（Tn）とその直前（１処理周期前）の温度検出値（Tn-1）との差ΔTn=(Tn)−(Tn-1)の絶対値（｜ΔTn｜）を算出する（ＳＴＥＰ１２）。
【００４３】
次に、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)より大きいか否かを判定する（ＳＴＥＰ２１）。ここで、第１閾値(Th1)は、温度検出値に含まれる可能性のあるノイズ(N)の大きさよりも大きな値に設定されている。このように、第１閾値(Th1)を設定することにより、温度検出値の変化のうち、ノイズ成分に起因する変化に拘らず、温度検出値自体の変化を捉えることができる。
【００４４】
そして、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)より大きい場合には（ＳＴＥＰ２１でＹＥＳ）、平均値を算出する際に用いる温度検出値の個数を第１個数（M1）とする（ＳＴＥＰ３１）。一方、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)以下の場合には（ＳＴＥＰ２１でＮＯ）、平均値を算出する際に用いる温度検出値の個数を第２個数（M2）とする（ＳＴＥＰ３２）。ここで、第１個数（M1）と第２個数（M2）では、M1＜＜M2の関係にある。第２個数（M2）は、当該蛍光温度センサにおいて移動平均値を算出する際に用いる通常の個数（例えば５００個）である。これに対して、第１個数（M1）は、第２個数（M2）の約２５％の個数（例えば１２５個）である。
【００４５】
このように、ノイズの大きさより大きな温度検出値の変化が生じた場合に、平均値を算出するために用いる個数を通常よりも少ない第１個数（M1）とすることで、ノイズ成分に平均値の出力を追従させることを回避しつつ、温度検出値の変化に平均値の出力を追従させ易くすることができる。一方で、ノイズに依る温度検出値の変化に対しては、第１個数（M1）より多い通常の第２個数（M2）で平均値を算出することで、ノイズの影響を低減させ、安定した温度の検出を行うことができる。
【００４６】
次いで、移動平均値算出部６の平均値算出手段６３は、最新の温度検出値（Tn）を含む、ＳＴＥＰ３１またはＳＴＥＰ３２で設定した個数の温度検出値で平均値を算出する（ＳＴＥＰ４０）。そして、移動平均値算出部６は、算出した平均値を当該温度センサの出力値（温度検出値）として出力する（ＳＴＥＰ４１）。
【００４７】
以上が、移動平均値算出部６の一連の処理であり、移動平均値算出部６は、これら一連の処理をＬＥＤ２の投光処理周期に対応した前記処理周期で逐次実行することにより、移動平均値を算出する。尚、上記実施形態において、ＳＴＥＰ２１の判定結果に基づいて、ＳＴＥＰ３１，３２で使用個数を設定する処理が、本発明の個数設定処理に相当し、ＳＴＥＰ４０の処理が、本発明の平均値算出処理に相当する。
【００４８】
次に、図３に示すフローチャートを参照して、図２に示すフローチャートの処理の変更例について説明する。なお、図３のフローチャートの処理は、図２のフローチャートの処理の一部を変更したものであるので、図２と同一処理については同一の参照符号を付して説明を省略する。
【００４９】
本実施例では、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)より大きい場合に（ＳＴＥＰ２１でＹＥＳ）、該差の絶対値が、第２閾値(Th2)より大きいか否かを判定する（ＳＴＥＰ２２）。ここで、第２閾値(Th2)は、第１閾値(Th1)よりも大きな値に設定されている。このように、第２閾値(Th2)を設定することにより、温度検出値の変化が著しい場合にもその変化を捉えることができる。
【００５０】
そして、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第２閾値(Th2)より大きい場合には（ＳＴＥＰ２２でＹＥＳ）、平均値を算出する際に用いる温度検出値の個数を１として、最新の温度検出値（Tn）を当該温度センサの出力値としてそのまま出力する（ＳＴＥＰ４２）。これにより、特に温度検出値の変化が著しい場合に、入力値の変化に出力を確実に追従させることができる。
【００５１】
一方、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第２閾値(Th2)以下の場合には（ＳＴＥＰ２２でＮＯ）、平均値を算出する際に用いる温度検出値の個数を第１個数（M1）として（ＳＴＥＰ３１）、以下、図２のフローチャートと同様に、第１個数（M1）の温度検出値で平均値を算出して（ＳＴＥＰ４０）、出力する（ＳＴＥＰ４１）。
【００５２】
以上が、移動平均値算出部６による処理の変更例である。尚、上記実施形態において、ＳＴＥＰ２１の判定結果に基づいてＳＴＥＰ３２で使用個数を設定する処理、およびＳＴＥＰ２２の判定結果に基づいてＳＴＥＰ３１，４２で使用個数を設定する処理が、本発明の個数設定処理に相当する。
【００５３】
このように、本実施形態の蛍光温度センサでは、温度検出値にノイズが含まれる場合にも、温度検出値の変化に対応した適切な平均値を逐次出力することができる。
【００５４】
尚、本実施形態では、差の絶対値（｜ΔTn｜）に対して第１および第２閾値(Th1,Th2)を設定し該閾値を超える場合に、使用個数を通常の第２個数（M2）を少ない第１個数（M1）等に設定しているが、これに限らず、閾値以内の場合に、使用個数を増加させてもよい。
【００５５】
また、本実施形態では、第１および第２の２つの閾値(Th1,Th2)を設けて、使用個数（M1,M2）を設定しているが、これに限らず、３つ以上の複数の閾値を設けて、段階的に個数を減少または増加させてもよい。これにより、温度検出値の変化に応じて最適な使用個数で平均値を逐次算出することができる。
【００５６】
さらに、本実施形態では、差（｜ΔTn｜）を最新の温度検出値（Tn）とその直前（１処理周期前）の温度度検出値（Tn-1）との差として算出しているが、これに限らず、１処理周期前の温度検出値（Tn-1）の代わりに、２処理周期前の温度検出値（Tn-2）や３処理周期前の温度検出値（Tn-3）を用いてもよい。これにより、１処理周期前の温度検出値（Tn-1）がノイズにより変化している場合にも、一定の間隔を空けることによりその影響を相対的に低減させることができる。
【００５７】
また、温度検出値の過去値（Tn-1，Tn-2）の代わりに、その平均値を用いてもよい。これにより、ノイズにより個々変化し得る過去の温度検出値に依らず、差を算出することができ、ノイズの影響を低減して、温度検出値の変化を確実に捉え、これに対応した適切な平均値の出力を行い、ひいては安定した温度出力を実現することができる。
【００５８】
以上、蛍光温度センサの温度検出値を入力値とする実施例について説明したが、本発明はこれに限らず、逐次入力される入力の平均を算出するものに適用できるものである。
【図面の簡単な説明】
【００５９】
【図１】本実施形態の蛍光温度センサの全体構成図。
【図２】移動平均値算出部における移動平均値を算出する処理を示すフローチャート。
【図３】図２に示すフローチャートの変更例を示す説明図。
【符号の説明】
【００６０】
１…蛍光材料、２…ＬＥＤ、３…駆動部、４…フォトダイオード、５…信号処理部、６…移動平均値算出部（移動平均値算出装置）、７…光ファイバ、１０…コントローラ、６１…記憶保持手段、６２…個数設定手段、６３…平均値算出手段。
【技術分野】
【０００１】
本発明は、逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出装置および移動平均値算出方法に関する。
【背景技術】
【０００２】
従来、移動平均値算出装置としては、特許文献１に示すように、時間平均を算出する区間ｓを、ｓ（＝β／α）として時間平均値を逐次算出するものが知られている。ここで、βは、時間的に変化する変化量をデジタル化する際の最小分解能であり、αは、時間平均値のしきい値である。
【０００３】
かかる移動平均値算出装置によれば、例えば、温度検出値のデジタル変換値そのものではなく、その移動平均値に基づいて温度の変化率を求めるので、該変化率の落ち込みは生じ難くなり、また落ち込みが生じたとしてもその程度が小さくなる。そのため、変化率のデジタル的なばらつきの程度が小さくなり、変化率がしきい値αを超えているか否かの判定精度が向上する。
【０００４】
しかし、従来の移動平均値算出装置では、例えば、加熱炉のように急激な温度上昇・下降が起きた後に安定した状態が続くような測定環境では、急激な温度変動に追随できず、オーバーシュートしたかどうかの検知が困難であった。
【０００５】
これに対して、特許文献２に示す移動平均値算出装置では、圧力計における急激な圧力変化に対して、蓄積した検出値の平均化を行わず、検出値をそのまま出力する。この移動平均値算出装置によれば、急激な圧力変化が生じた場合にも、検出値がそのまま出力されるため、その変化に出力を追従させることができる。
【特許文献１】特開２０００−１２４８０６号公報
【特許文献２】特開平１１−１１８６５１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の移動平均値算出装置では、検出値の変動が、センサ故障や劣化に起因するノイズや一般的に生じる電気ノイズ等（以下、ノイズという）によるものである場合にも、これがそのまま出力されてしまうという問題がある。
【０００７】
これは、検出値を蓄積してその平均値を出力することにより、センサの出力がノイズによって急変することを回避して出力の安定化を図るという、移動平均を用いる本来の目的に反することとなる。
【０００８】
上記の事情に鑑みて、本発明は、入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値の出力を行うことができる移動平均値算出装置および移動平均値算出方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
第１発明の移動平均値算出装置は、逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出装置であって、前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定手段を備えることを特徴とする。
【００１０】
第１発明の移動平均値算出装置によれば、最新の入力値と過去の入力値との差に基づいて移動平均値を算出するために用いる個数を設定する。そのため、最新の入力値の変化が大きい場合には、例えば、入力値の個数を減少させて、この変化に平均値の出力を追従させ易くすることができる。一方、最新の入力値の変化が小さい場合には、例えば、入力値の個数を増加させて、ノイズの影響を低減させることにより安定した平均値の出力を行うことができる。
【００１１】
ここで、移動平均値を算出するために用いる個数は、前記差とノイズの大きさから定められる閾値との大小関係により設定される。そのため、例えば、閾値をノイズの大きさよりも大きくしておくことで、その閾値を超える入力値の変化にのみ平均値の出力を追従させ易くすることが可能となる。これにより、ノイズの影響を確実に低減して、入力値の変化に対応した適切な平均値の出力を行うことができる。
【００１２】
また、閾値は上述したノイズの大きさに基づいて定めるもののほか、ＳＮの比に基づき定めるもの、あるいは分解能／精度／再現性を考慮して定めるものなどが考えられる。いずれにせよ、ノイズの影響を低減して、入力値の変化に対応した適切な平均値の出力を行えるような値を閾値として設定すればよい。
【００１３】
なお、差を算出する際に用いる過去の入力値は、最新の入力値の直前（１周期前）の入力値でもよいが、該直前の入力値以外の値であってもよい。例えば、直前の入力値の更に直前の入力値（最新の入力値の２周期前に入力された入力値）や更にその直前の入力値（最新の入力値の３周期前に入力された入力値）を用いてもよい。これにより、最新の入力値の直前（１周期前）の入力値がノイズにより変化している場合にも、一定の間隔を空けることによりその影響を相対的に低減することができる。
【００１４】
第２発明の移動平均値算出装置は、第１発明の移動平均値算出装置において、前記個数設定手段は、前記所定の閾値である第１閾値に対して、前記差が該第１閾値より大きい場合には、前記入力値の個数を第１個数に設定し、該第１閾値以下である場合には、該個数を第１個数よりも多い第２個数に設定することを特徴とする。
【００１５】
第２発明の移動平均値算出装置によれば、第１閾値を、例えばノイズの大きさよりも大きい値に設定しておくことで、前記差が、第１閾値を超える入力値の変化である場合に、移動平均値を算出するために用いる個数が少ない第１個数で平均値を算出する。そのため、ノイズ成分に平均値の出力を追従させることなく、入力値の変化に平均値の出力を追従させることができる。一方、前記差が、第１閾値を以下の入力値の変化である場合に、該入力値の変化はノイズ成分に起因するものが含まれ得る。この場合、移動平均値を算出するために用いる個数が第１個数に比して多い第２個数で平均値を算出することにより、ノイズの影響を低減させることにより安定した平均値の出力を行うことができる。これにより、入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値の出力を行うことができる。
【００１６】
第３発明の移動平均値算出装置は、第２発明の移動平均値算出装置において、前記個数設定手段は、前記第１閾値よりも大きな値に設定された第２閾値に対して、前記差が該第２閾値より大きい場合には、前記第１個数を１とすることを特徴とする。
【００１７】
第３発明の移動平均値算出装置によれば、第１閾値よりも大きい第２閾値を設定し、前記差が、第２閾値を超える入力値の変化である場合に、移動平均値を算出するために用いる個数を１として、入力値をそのまま出力する。これにより、特に入力値の変化が著しい場合に、ノイズ成分に出力が追従することを回避しつつ、入力値の変化に出力を確実に追従させることができる。
【００１８】
第４発明の移動平均値算出装置は、第１〜第３発明のいずれかの移動平均値算出装置において、前記差は、前記最新の入力値と過去の入力値の平均値との差であることを特徴とする。
【００１９】
第４発明の移動平均値算出装置によれば、差を算出する際に用いる過去の入力値として過去の入力値の平均値を用いることで、ノイズにより個々変化し得る過去の入力値に依らず、差を算出することができる。そのため、ノイズの影響を低減して、入力値の変化を確実に捉え、これに対応した適切な平均値の出力を行うことができる。
【００２０】
第５発明の蛍光温度センサは、蛍光材料に投光する発光素子と、前記蛍光材料が発する蛍光を受光する受光素子と、前記受光素子の出力信号から温度検出値を生成する信号処理部とを含む蛍光温度センサにおいて、前記温度検出値を入力としてその移動平均値を出力する移動平均値算出部が請求項１乃至４のうちいずれか１項記載の移動平均値算出装置により構成されてなることを特徴とする。
【００２１】
第５発明によれば、蛍光温度センサは、一定の処理周期で取得された温度検出値に対して、移動平均値を算出して出力するため、第１〜第４発明の移動平均値算出装置を適用する対象として好適である。
【００２２】
第６発明の移動平均値算出方法は、逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出方法であって、前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定処理と、前記個数設定処理により設定された個数の入力値を用いて、これらの平均値を算出する平均値算出処理とを逐次実行することを特徴とする。
【００２３】
第６発明の移動平均値算出方法によれば、最新の入力値と過去の入力値との差に基づいて移動平均値を算出するために用いる個数を設定し、設定した個数の入力値を用いて平均値を逐次算出する。そのため、最新の入力値の変化が大きい場合には、入力値の個数を減少させて、この変化に平均値の出力を追従させ易くすることや、最新の入力値の変化が小さい場合には、入力値の個数を増加させて、ノイズの影響を低減させることにより安定した出力を行うことができる。
【００２４】
さらに、移動平均値を算出するために用いる個数は、前記差とノイズの大きさから定められる閾値との大小関係により設定されるため、閾値を例えばノイズの大きさよりも大きくしておくことで、その閾値を超える入力値の変化にのみ平均値の出力を追従させ易くることが可能となる。なお、閾値について種々の設定方法が考えられることは前述したとおりである。
【００２５】
このように、入力値にノイズが含まれる場合にも、入力値の変化に対応した適切な平均値を算出して出力することができる。
【発明を実施するための最良の形態】
【００２６】
本発明の一実施形態として、蛍光温度センサの移動平均値算出部に移動平均値算出装置を用いた例について、図１〜図３を参照して説明する。
【００２７】
まず、図１を参照して、本実施形態の蛍光温度センサの全体的な構成について説明する。蛍光温度センサは、温度によって異なる蛍光特性を示す蛍光材料１と、蛍光材料１に投光するＬＥＤ２と、ＬＥＤ２を駆動する電気信号を出力する駆動部３と、蛍光材料１が発する蛍光を受光するフォトダイオード４と、フォトダイオード４からの出力信号から対応する温度検出値を算出する信号処理部５と、信号処理部５に接続されて温度検出値の移動平均値を算出し、温度を表すデータとして出力する移動平均値算出部６（本発明の移動平均値算出装置に相当する）とを備える。
【００２８】
また、蛍光温度センサは、蛍光材料１への投光および蛍光材料１の蛍光の受光を行う光ファイバ７と備える。光ファイバ７は他端側が分岐して、ＬＥＤ２からの光を蛍光材料１に伝達する投光用光ファイバ７ａと、蛍光材料１の蛍光をフォトダイオード４に伝達する受光用光ファイバ７ｂとなっている。
【００２９】
蛍光材料１は、光ファイバ７の一端部を覆うように設けられた温度計測部１ａの中に、光ファイバ７のコア部に対向するように配置される。
【００３０】
ＬＥＤ２は、ＬＥＤモジュール２ａ内に配置された、例えば青色系の波長を発光色とする発光ダイオードである。ＬＥＤモジュール２ａは、投光用光ファイバ７ａが接続されるコネクタ部２ｂを有し、コネクタ部２ｂを介して接続された投光用光ファイバ７ａがＬＥＤ２の発光部と対向している。
【００３１】
駆動部３は、ＬＥＤ２の発光に必要な駆動電流の大きさおよび発光時間を規定したパルス電流を一定の処理周期でＬＥＤ２に印加する制御回路を備える。これにより、駆動部３は、例えば、蛍光材料１に対応して、一回の計測におけるＬＥＤ２の発光時間を１ｍｓ〜５００ｍｓの間のいずれかの時間とする所定の大きさのパルス電流をＬＥＤ２に印加する。
【００３２】
フォトダイオード４は、フォトダイオードモジュール４ａ内に配置されて、照射された光の光量（輝度）を測定する。フォトダイオードモジュール４ａは、受光用光ファイバ７ｂに接続されるコネクタ部４ｂを有し、コネクタ部４ｂを介して接続された受光用光ファイバ７ｂがフォトダイオード４の受光部と対向している。
【００３３】
信号処理部５は、前記処理周期（ＬＥＤ２の投光処理周期に対応した処理周期）で、フォトダイオード４によって測定された蛍光材料１の蛍光の減衰特性、特に蛍光緩和時間を計測する。具体的には、信号処理部５は、これがあらかじめ備える蛍光緩和時間と蛍光材料１との関係式（データテーブルやマップ等を含む）から、蛍光材料１が存在する温度測定環境の温度を算出して出力する。
【００３４】
移動平均値算出部６は、記憶保持手段６１と、個数設定手段６２と、平均値算出手段６３とを備え、信号処理部５から逐次入力される温度検出値の平均値を前記処理周期で逐次算出して出力する。
【００３５】
記憶保持手段６１は、図示しない内部メモリ等により構成され、信号処理部５から逐次入力された温度検出値を最新のものから一定の個数だけ記憶保持する。
【００３６】
個数設定手段６２は、記憶保持手段６１に記憶保持された温度検出値を読出し、最新の温度検出値と過去の温度検出値との差を算出し、該差と予め設定された閾値との大小関係により、平均値を算出する際に用いる温度検出値の個数を設定する。ここで、閾値は、温度検出値に含まれ得るノイズの大きさから定められる（図２のＳＴＥＰ２１および図３のＳＴＥＰ２２参照）。
【００３７】
平均値算出手段６３は、最新の温度検出値を含めて、個数設定手段６２により設定された個数の温度検出値を用いて、これらの平均値を算出する処理を逐次実行する。
【００３８】
なお、本実施形態において、各処理部３，５，６は、マイクロコンピュータ等により構成されたコントローラ１０に備えられてなる。コントローラ１０を構成する各処理部３，５，６は、ＣＰＵ、ＲＯＭ、ＲＡＭ等のハードウェアにより構成され、これらの各処理部３，５，６が共通のハードウェアによって構成されていてもよく、これらの各処理部３，５，６の一部又は全部が異なるハードウェアによって構成されていてもよい。
【００３９】
次に、図２に示すフローチャートを参照して、移動平均値算出部６により移動平均値を算出する処理について説明する。
【００４０】
まず、移動平均値算出部６は、前記処理周期毎に信号処理部５から出力される最新の温度検出値（Tn）を逐次取得する（ＳＴＥＰ１０）。ここで、温度検出値（T）の添え字n（n：自然数）は、信号処理部５から取得した温度検出値の順序を表し、nが小さいほど処理周期が過去のものを表す。例えば、１処理周期前の温度検出値は、Tn-1となる。
【００４１】
そして、移動平均値算出部６は、取得した温度検出値（Tn）を新たに記憶保持手段６１に記憶保持する（ＳＴＥＰ１１）。
【００４２】
次いで、移動平均値算出部６の個数設定手段６２は、記憶保持手段６１に記憶保持されている温度検出値（T1〜Tn）を読出し、最新の温度検出値（Tn）とその直前（１処理周期前）の温度検出値（Tn-1）との差ΔTn=(Tn)−(Tn-1)の絶対値（｜ΔTn｜）を算出する（ＳＴＥＰ１２）。
【００４３】
次に、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)より大きいか否かを判定する（ＳＴＥＰ２１）。ここで、第１閾値(Th1)は、温度検出値に含まれる可能性のあるノイズ(N)の大きさよりも大きな値に設定されている。このように、第１閾値(Th1)を設定することにより、温度検出値の変化のうち、ノイズ成分に起因する変化に拘らず、温度検出値自体の変化を捉えることができる。
【００４４】
そして、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)より大きい場合には（ＳＴＥＰ２１でＹＥＳ）、平均値を算出する際に用いる温度検出値の個数を第１個数（M1）とする（ＳＴＥＰ３１）。一方、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)以下の場合には（ＳＴＥＰ２１でＮＯ）、平均値を算出する際に用いる温度検出値の個数を第２個数（M2）とする（ＳＴＥＰ３２）。ここで、第１個数（M1）と第２個数（M2）では、M1＜＜M2の関係にある。第２個数（M2）は、当該蛍光温度センサにおいて移動平均値を算出する際に用いる通常の個数（例えば５００個）である。これに対して、第１個数（M1）は、第２個数（M2）の約２５％の個数（例えば１２５個）である。
【００４５】
このように、ノイズの大きさより大きな温度検出値の変化が生じた場合に、平均値を算出するために用いる個数を通常よりも少ない第１個数（M1）とすることで、ノイズ成分に平均値の出力を追従させることを回避しつつ、温度検出値の変化に平均値の出力を追従させ易くすることができる。一方で、ノイズに依る温度検出値の変化に対しては、第１個数（M1）より多い通常の第２個数（M2）で平均値を算出することで、ノイズの影響を低減させ、安定した温度の検出を行うことができる。
【００４６】
次いで、移動平均値算出部６の平均値算出手段６３は、最新の温度検出値（Tn）を含む、ＳＴＥＰ３１またはＳＴＥＰ３２で設定した個数の温度検出値で平均値を算出する（ＳＴＥＰ４０）。そして、移動平均値算出部６は、算出した平均値を当該温度センサの出力値（温度検出値）として出力する（ＳＴＥＰ４１）。
【００４７】
以上が、移動平均値算出部６の一連の処理であり、移動平均値算出部６は、これら一連の処理をＬＥＤ２の投光処理周期に対応した前記処理周期で逐次実行することにより、移動平均値を算出する。尚、上記実施形態において、ＳＴＥＰ２１の判定結果に基づいて、ＳＴＥＰ３１，３２で使用個数を設定する処理が、本発明の個数設定処理に相当し、ＳＴＥＰ４０の処理が、本発明の平均値算出処理に相当する。
【００４８】
次に、図３に示すフローチャートを参照して、図２に示すフローチャートの処理の変更例について説明する。なお、図３のフローチャートの処理は、図２のフローチャートの処理の一部を変更したものであるので、図２と同一処理については同一の参照符号を付して説明を省略する。
【００４９】
本実施例では、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第１閾値(Th1)より大きい場合に（ＳＴＥＰ２１でＹＥＳ）、該差の絶対値が、第２閾値(Th2)より大きいか否かを判定する（ＳＴＥＰ２２）。ここで、第２閾値(Th2)は、第１閾値(Th1)よりも大きな値に設定されている。このように、第２閾値(Th2)を設定することにより、温度検出値の変化が著しい場合にもその変化を捉えることができる。
【００５０】
そして、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第２閾値(Th2)より大きい場合には（ＳＴＥＰ２２でＹＥＳ）、平均値を算出する際に用いる温度検出値の個数を１として、最新の温度検出値（Tn）を当該温度センサの出力値としてそのまま出力する（ＳＴＥＰ４２）。これにより、特に温度検出値の変化が著しい場合に、入力値の変化に出力を確実に追従させることができる。
【００５１】
一方、個数設定手段６２は、差の絶対値（｜ΔTn｜）が、第２閾値(Th2)以下の場合には（ＳＴＥＰ２２でＮＯ）、平均値を算出する際に用いる温度検出値の個数を第１個数（M1）として（ＳＴＥＰ３１）、以下、図２のフローチャートと同様に、第１個数（M1）の温度検出値で平均値を算出して（ＳＴＥＰ４０）、出力する（ＳＴＥＰ４１）。
【００５２】
以上が、移動平均値算出部６による処理の変更例である。尚、上記実施形態において、ＳＴＥＰ２１の判定結果に基づいてＳＴＥＰ３２で使用個数を設定する処理、およびＳＴＥＰ２２の判定結果に基づいてＳＴＥＰ３１，４２で使用個数を設定する処理が、本発明の個数設定処理に相当する。
【００５３】
このように、本実施形態の蛍光温度センサでは、温度検出値にノイズが含まれる場合にも、温度検出値の変化に対応した適切な平均値を逐次出力することができる。
【００５４】
尚、本実施形態では、差の絶対値（｜ΔTn｜）に対して第１および第２閾値(Th1,Th2)を設定し該閾値を超える場合に、使用個数を通常の第２個数（M2）を少ない第１個数（M1）等に設定しているが、これに限らず、閾値以内の場合に、使用個数を増加させてもよい。
【００５５】
また、本実施形態では、第１および第２の２つの閾値(Th1,Th2)を設けて、使用個数（M1,M2）を設定しているが、これに限らず、３つ以上の複数の閾値を設けて、段階的に個数を減少または増加させてもよい。これにより、温度検出値の変化に応じて最適な使用個数で平均値を逐次算出することができる。
【００５６】
さらに、本実施形態では、差（｜ΔTn｜）を最新の温度検出値（Tn）とその直前（１処理周期前）の温度度検出値（Tn-1）との差として算出しているが、これに限らず、１処理周期前の温度検出値（Tn-1）の代わりに、２処理周期前の温度検出値（Tn-2）や３処理周期前の温度検出値（Tn-3）を用いてもよい。これにより、１処理周期前の温度検出値（Tn-1）がノイズにより変化している場合にも、一定の間隔を空けることによりその影響を相対的に低減させることができる。
【００５７】
また、温度検出値の過去値（Tn-1，Tn-2）の代わりに、その平均値を用いてもよい。これにより、ノイズにより個々変化し得る過去の温度検出値に依らず、差を算出することができ、ノイズの影響を低減して、温度検出値の変化を確実に捉え、これに対応した適切な平均値の出力を行い、ひいては安定した温度出力を実現することができる。
【００５８】
以上、蛍光温度センサの温度検出値を入力値とする実施例について説明したが、本発明はこれに限らず、逐次入力される入力の平均を算出するものに適用できるものである。
【図面の簡単な説明】
【００５９】
【図１】本実施形態の蛍光温度センサの全体構成図。
【図２】移動平均値算出部における移動平均値を算出する処理を示すフローチャート。
【図３】図２に示すフローチャートの変更例を示す説明図。
【符号の説明】
【００６０】
１…蛍光材料、２…ＬＥＤ、３…駆動部、４…フォトダイオード、５…信号処理部、６…移動平均値算出部（移動平均値算出装置）、７…光ファイバ、１０…コントローラ、６１…記憶保持手段、６２…個数設定手段、６３…平均値算出手段。
【特許請求の範囲】
【請求項１】
逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出装置であって、
前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定手段を備えることを特徴とする移動平均値算出装置。
【請求項２】
請求項１記載の移動平均値算出装置において、
前記個数設定手段は、前記所定の閾値である第１閾値に対して、前記差が該第１閾値より大きい場合には、前記入力値の個数を第１個数に設定し、該第１閾値以下である場合には、該個数を第１個数よりも多い第２個数に設定することを特徴とする移動平均値算出装置。
【請求項３】
請求項２記載の移動平均値算出装置において、
前記個数設定手段は、前記第１閾値よりも大きな値に設定された第２閾値に対して、前記差が該第２閾値より大きい場合には、前記第１個数を１とすることを特徴とする移動平均値算出装置。
【請求項４】
請求項１乃至３のうちいずれか１項記載の移動平均値算出装置において、
前記差は、前記最新の入力値と過去の入力値の平均値との差であることを特徴とする移動平均値算出装置。
【請求項５】
蛍光材料に投光する発光素子と、前記蛍光材料が発する蛍光を受光する受光素子と、前記受光素子の出力信号から温度検出値を生成する信号処理部とを含む蛍光温度センサにおいて、
前記温度検出値を入力としてその移動平均値を出力する移動平均値算出部が請求項１乃至４のうちいずれか１項記載の移動平均値算出装置により構成されてなることを特徴とする蛍光温度センサ。
【請求項６】
逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出方法であって、
前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定処理と、
前記個数設定処理により設定された個数の入力値を用いて、これらの平均値を算出する平均値算出処理と
を逐次実行することを特徴とする移動平均値算出方法。
【請求項１】
逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出装置であって、
前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定手段を備えることを特徴とする移動平均値算出装置。
【請求項２】
請求項１記載の移動平均値算出装置において、
前記個数設定手段は、前記所定の閾値である第１閾値に対して、前記差が該第１閾値より大きい場合には、前記入力値の個数を第１個数に設定し、該第１閾値以下である場合には、該個数を第１個数よりも多い第２個数に設定することを特徴とする移動平均値算出装置。
【請求項３】
請求項２記載の移動平均値算出装置において、
前記個数設定手段は、前記第１閾値よりも大きな値に設定された第２閾値に対して、前記差が該第２閾値より大きい場合には、前記第１個数を１とすることを特徴とする移動平均値算出装置。
【請求項４】
請求項１乃至３のうちいずれか１項記載の移動平均値算出装置において、
前記差は、前記最新の入力値と過去の入力値の平均値との差であることを特徴とする移動平均値算出装置。
【請求項５】
蛍光材料に投光する発光素子と、前記蛍光材料が発する蛍光を受光する受光素子と、前記受光素子の出力信号から温度検出値を生成する信号処理部とを含む蛍光温度センサにおいて、
前記温度検出値を入力としてその移動平均値を出力する移動平均値算出部が請求項１乃至４のうちいずれか１項記載の移動平均値算出装置により構成されてなることを特徴とする蛍光温度センサ。
【請求項６】
逐次入力される入力値の平均値を最新の入力値を用いて逐次算出する移動平均値算出方法であって、
前記最新の入力値と過去の入力値との差を算出し、該差と所定の閾値との大小関係により、前記平均値を算出する際に用いる入力値の個数を設定する個数設定処理と、
前記個数設定処理により設定された個数の入力値を用いて、これらの平均値を算出する平均値算出処理と
を逐次実行することを特徴とする移動平均値算出方法。
【図１】
【図２】
【図３】
【図２】
【図３】
【公開番号】特開２００９−２９４１８０（Ｐ２００９−２９４１８０Ａ）
【公開日】平成２１年１２月１７日（２００９．１２．１７）
【国際特許分類】
【出願番号】特願２００８−１５０６８７（Ｐ２００８−１５０６８７）
【出願日】平成２０年６月９日（２００８．６．９）
【出願人】（０００００６６６６）株式会社山武 (1,808)
【Ｆターム（参考）】
【公開日】平成２１年１２月１７日（２００９．１２．１７）
【国際特許分類】
【出願日】平成２０年６月９日（２００８．６．９）
【出願人】（０００００６６６６）株式会社山武 (1,808)
【Ｆターム（参考）】
[ Back to top ]