"製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する前記製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測装置であって、前記製造プロセスにおける過去の実績データを記憶する実績データ記憶部から、実績データを抽出するデータ抽出部と、前記データ抽出部により抽出された前記実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成部と、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出し、抽出された前記実績データに基づいて、前記品質予測対象製品の品質を予測する予測部と、を備え、前記検索用テーブル作成部は、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割部と、前記分割候補点で分割された前記プロセス変数値空間の各局所領域について、前記品質データを確率変数とする確率密度関数を用いて、前記局所領域における前記品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出部と、前記各局所領域について算出された前記バラツキ指標が最小となる局所領域を生成する前記分割候補点を確定分割点として決定し、前記確定分割点により前記プロセス変数値空間を確定分割するモデル選択部と、前記モデル選択部により前記確定分割点にて確定分割された前記プロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定部と、前記プロセス変数値空間の局所領域の分割結果を前記検索用テーブルとして出力する分割結果出力部と、を備えて、前記再分割対象領域判定部により再分割が必要と判定された局所領域に対して、前記プロセス変数値空間分割部による分割候補点設定と、前記確率密度関数算出部によるバラツキ指標算出と、前記モデル選択部による確定分割と、前記再分割対象領域判定部による再分割の要否判定とを繰り返し、前記予測部は、前記品質予測対象製品のプロセス変数値と同一の前記局所領域に存在するプロセス変数値を類似のプロセス変数値として、前記品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出することを特徴とする、品質予測装置。前記予測部は、前記検索用テーブルから、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データが複数抽出されたとき、抽出された前記実績データの各品質データに対して処理時刻に応じた重み付けを行い、重み付けされた前記各品質データの平均値を品質予測対象製品の予測品質データとすることを特徴とする、請求項1に記載の品質予測装置。前記予測部は、前記検索用テーブルから抽出された前記実績データの各品質データに対する重み付けを、現在時刻と処理時刻との時刻差が小さいほど大きく設定することを特徴とする、請求項2に記載の品質予測装置。前記予測部は、前記検索用テーブルから、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データが複数抽出されたとき、抽出された複数の前記実績データのプロセス変数値および品質データに基づいて、未定係数を乗じたプロセス変数値の線形和として品質データを推定する重回帰モデルの未定係数値を算出し、前記重回帰モデルに品質予測対象製品のプロセス変数値を設定して取得される値を予測品質データとすることを特徴とする、請求項1に記載の品質予測装置。前記予測部は、前記検索用テーブルから、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データが複数抽出されたとき、抽出された複数の前記実績データのプロセス変数値および品質データに基づいて、予め設定された確率密度関数のパラメータを算出し、前記確率密度関数モデルに該算出されたパラメータと予め設定した確率値を設定して取得される値を予測品質データとすることを特徴とする、請求項1に記載の品質予測装置。前記予測部は、生成した前記予測品質データをユーザに提示する出力装置に出力することを特徴とする、請求項1~5のいずれか1項に記載の品質予測装置。前記バラツキ指標は、前記確率密度関数の累積確率が予め設定された値であるときの確率変数値と確率分布の最頻値との差分であり、前記モデル選択部は、前記各局所領域における、前記バラツキ指標の最大値が最小である場合に、前記バラツキ指標が最小であるとすることを特徴とする、請求項1~6のいずれか1項に記載の品質予測装置前記バラツキ指標は、前記局所領域内の前記品質データの点数と該品質データの分散との積であり、前記モデル選択部は、前記各局所領域における、前記品質データの点数と該品質データの分散との積の合計値が最小である場合に、前記バラツキ指標が最小であるとすることを特徴とする、請求項1~6のいずれか1項に記載の品質予測装置。前記バラツキ指標は、前記局所領域内の前記品質データの分散であり、前記モデル選択部は、前記各局所領域における、前記品質データの分散の最大値が最小である場合に、前記バラツキ指標が最小であるとすることを特徴とする、請求項1~6のいずれか1項に記載の品質予測装置。前記プロセス変数値空間分割部は、前記プロセス変数値空間を、前記分割候補点においてそれぞれ2つの局所領域に分割し、前記確率密度関数算出部は、前記分割候補点で分割された前記局所領域のうち少なくともいずれか一方の前記局所領域に属するデータ点数が基準値未満となるとき、当該分割候補点を候補から除外することを特徴とする、請求項1~9のいずれか1項に記載の品質予測装置。前記再分割対象領域判定部は、前記プロセス変数値空間の分割回数が予め設定された最大分割数に達したとき、分割処理を終了することを特徴とする、請求項1~10のいずれか1項に記載の品質予測装置。前記プロセス変数値空間分割部は、前記局所領域の再分割時において、前記プロセス変数値空間分割部により設定された前記分割候補点のうち分割しようとする前記局所領域の領域外となる前記分割候補点を候補から除外することを特徴とする、請求項1~11のいずれか1項に記載の品質予測装置。前記請求項1~12のいずれか1項に記載の品質予測装置を用いて、製造途中の製品の予測品質データを生成し、生成した前記予測品質データに基づき製造される製品が所望の品質データとなるように前記製造プロセスの操業条件を決定する操業条件決定方法であって、前記データ抽出部により、前記製造プロセスにおける過去の実績データを記憶する実績データ記憶部から、既に工程が終了し実績が確定したプロセス変数値のみを確定実績データとして抽出し、前記検索用テーブル作成部により、前記確定実績データに基づいて前記検索用テーブルを作成し、前記予測部により、品質予測対象製品の予測品質データを出力するために必要な、前記確定実績データに含まれる操業条件以外の、決定対象である操業条件の決定されるべき値の候補である代表値を複数作成し、前記各代表値と前記確定実績データとの組合せについて前記予測品質データを生成して、所望の品質データに最も近い前記予測品質データとなった前記代表値を前記製造プロセスの操業条件として決定することを特徴とする、操業条件決定方法。前記プロセス変数値の代表値は、前記各プロセス変数値に対して前記検索用テーブルの作成時に採用された分割候補点で区分された各区間における代表値であることを特徴とする、請求項13に記載の操業条件決定方法。前記プロセス変数値の代表値は、前記各プロセス変数値に対して前記検索用テーブルの作成時に採用された分割候補点で区分された各区間における中間値であることを特徴とする、請求項14に記載の操業条件決定方法。前記予測部により、前記確定実績データと、当該確定実績データに含まれるプロセス変数値以外のプロセス変数値が取り得る代表値とからなる複数の入力データを作成し、前記各入力データに基づいて、前記検索用テーブルから当該入力データに類似する実績データを取得し、取得した前記実績データを未確定のプロセス変数値の取り得る代表値毎に層別して、それぞれについて前記予測品質データを生成し、所望の品質データに最も近い前記予測品質データの前記代表値を前記製造プロセスの操業条件として決定することを特徴とする、請求項13~15のいずれか1項に記載の操業条件決定方法。製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する前記製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測方法であって、前記製造プロセスにおける過去の実績データを記憶する実績データ記憶部から、実績データを抽出するデータ抽出ステップと、前記データ抽出ステップにより抽出された前記実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成ステップと、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出し、抽出された前記実績データに基づいて、前記品質予測対象製品の品質を予測する予測ステップと、を含み、前記検索用テーブル作成ステップは、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割ステップと、前記分割候補点で分割された前記プロセス変数値空間の各局所領域について、前記品質データを確率変数とする確率密度関数を用いて、前記局所領域における前記品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出ステップと、前記各局所領域について算出された前記バラツキ指標が最小となる局所領域を生成する前記分割候補点を確定分割点として決定し、前記確定分割点により前記プロセス変数値空間を確定分割するモデル選択ステップと、前記モデル選択ステップにて前記確定分割点にて確定分割された前記プロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定ステップと、前記プロセス変数値空間の局所領域の分割結果を前記検索用テーブルとして出力する分割結果出力ステップと、を含んで、前記再分割対象領域判定ステップにより再分割が必要と判定された局所領域に対して、前記プロセス変数値空間分割ステップと、前記確率密度関数算出ステップと、前記モデル選択ステップと、前記再分割対象領域判定ステップとを繰り返し、前記予測ステップは、前記品質予測対象製品のプロセス変数値と同一の前記局所領域に存在するプロセス変数値を類似のプロセス変数値として、前記品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出することを特徴とする、品質予測方法。コンピュータを、製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する前記製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測装置として機能させるためのコンピュータプログラムであって、前記製造プロセスにおける過去の実績データを記憶する実績データ記憶手段から、実績データを抽出するデータ抽出手段と、前記データ抽出手段により抽出された前記実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成手段と、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出し、抽出された前記実績データに基づいて、前記品質予測対象製品の品質を予測する予測手段と、を備え、前記検索用テーブル作成手段は、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割手段と、前記分割候補点で分割された前記プロセス変数値空間の各局所領域について、前記品質データを確率変数とする確率密度関数を用いて、前記局所領域における前記品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出手段と、前記各局所領域について算出された前記バラツキ指標が最小となる局所領域を生成する前記分割候補点を確定分割点として決定し、前記確定分割点により前記プロセス変数値空間を確定分割するモデル選択手段と、前記モデル選択手段により前記確定分割点にて確定分割された前記プロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定手段と、前記プロセス変数値空間の局所領域の分割結果を前記検索用テーブルとして出力する分割結果出力手段と、を備えて、前記再分割対象領域判定手段により再分割が必要と判定された局所領域に対して、前記プロセス変数値空間分割手段による分割候補点設定と、前記確率密度関数算出手段によるバラツキ指標算出と、前記モデル選択手段による確定分割と、前記再分割対象領域判定手段による再分割の要否判定とを繰り返し、前記予測手段は、前記品質予測対象製品のプロセス変数値と同一の前記局所領域に存在するプロセス変数値を類似のプロセス変数値として、前記品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出することを特徴とする品質予測装置として機能させるためのコンピュータプログラム。コンピュータに、製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する前記製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測装置として機能させるためのプログラムを記憶したコンピュータ読み取り可能な記憶媒体であって、前記製造プロセスにおける過去の実績データを記憶する実績データ記憶手段から、実績データを抽出するデータ抽出手段と、前記データ抽出手段により抽出された前記実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成手段と、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出し、抽出された前記実績データに基づいて、前記品質予測対象製品の品質を予測する予測手段と、を備え、前記検索用テーブル作成手段は、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割手段と、前記分割候補点で分割された前記プロセス変数値空間の各局所領域について、前記品質データを確率変数とする確率密度関数を用いて、前記局所領域における前記品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出手段と、前記各局所領域について算出された前記バラツキ指標が最小となる局所領域を生成する前記分割候補点を確定分割点として決定し、前記確定分割点により前記プロセス変数値空間を確定分割するモデル選択手段と、前記モデル選択手段により前記確定分割点にて確定分割された前記プロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定手段と、前記プロセス変数値空間の局所領域の分割結果を前記検索用テーブルとして出力する分割結果出力手段と、を備えて、前記再分割対象領域判定手段により再分割が必要と判定された局所領域に対して、前記プロセス変数値空間分割手段による分割候補点設定と、前記確率密度関数算出手段によるバラツキ指標算出と、前記モデル選択手段による確定分割と、前記再分割対象領域判定手段による再分割の要否判定とを繰り返し、前記予測手段は、前記品質予測対象製品のプロセス変数値と同一の前記局所領域に存在するプロセス変数値を類似のプロセス変数値として、前記品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する前記実績データを前記検索用テーブルから抽出することを特徴とする品質予測装置として機能させるためのプログラムを記憶したコンピュータ読み取り可能な記憶媒体。しかし、上記特許文献1では、プロセスの材料値に影響を及ぼす入力変数の数が増加すると、量子化された領域の数が指数的に増加するため、各領域に含まれる過去データの数が減少してしまう。さらには、過去データが存在しない領域が増加してしまい、予測精度が悪化するという問題があった。また、領域数の増加に伴い検索処理に要する時間が増え、品質予測処理の高速性が失われてしまう。上記特許文献2では、入力される素材成分や操業条件の指示値と蓄積されたデータを用いて品質予測を行っているが、ルールベースで行っている入力データの削減が製造プロセスの設備更新や経年変化によって適切ではなくなっている状況に対応していない。また、材質推定モデルを線形モデルで近似しているため、大規模、複雑、非線形かつ非定常なプロセスで製造される製品の品質予測においては十分な精度で予測できないという問題があった。さらに、大規模、複雑、非線形かつ非定常なプロセスの製品を予測するためには大量のデータが必要であるが、データベースに蓄積された過去事例が膨大になるにつれ、事例の検索に要する時間が増加し、オンラインでの予測に必要な高速性を維持することができないという問題もある。そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、大規模、複雑、非線形かつ非定常なプロセスの実績データから類似例を高速に検索し、製品の予測品質データを高精度に生成することが可能な、新規かつ改良された品質予測装置、操業条件決定方法、品質予測方法、コンピュータプログラムおよびコンピュータ読み取り可能な記憶媒体を提供することにある。上記課題を解決するために、本発明のある観点によれば、製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測装置が提供される。品質予測装置は、製造プロセスにおける過去の実績データを記憶する実績データ記憶部から実績データを抽出するデータ抽出部と、データ抽出部により抽出された実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成部と、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出し、抽出された実績データに基づいて、品質予測対象製品の品質を予測する予測部と、を備える。検索用テーブル作成部は、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割部と、分割候補点で分割されたプロセス変数値空間の各局所領域について、品質データを確率変数とする確率密度関数を用いて、局所領域における品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出部と、各局所領域について算出されたバラツキ指標が最小となる局所領域を生成する分割候補点を確定分割点として決定し、確定分割点によりプロセス変数値空間を確定分割するモデル選択部と、モデル選択部により確定分割点にて確定分割されたプロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定部と、プロセス変数値空間の局所領域の分割結果を検索用テーブルとして出力する分割結果出力部と、を備えて、再分割対象領域判定部により再分割が必要と判定された局所領域に対して、プロセス変数値空間分割部による分割候補点設定と、確率密度関数算出部によるバラツキ指標算出と、モデル選択部による確定分割と、再分割対象領域判定部による再分割の要否判定とを繰り返し、予測部は、品質予測対象製品のプロセス変数値と同一の局所領域に存在するプロセス変数値を類似のプロセス変数値として、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出することを特徴とする。予測部は、検索用テーブルから、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データが複数抽出されたとき、抽出された実績データの各品質データに対して処理時刻に応じた重み付けを行い、重み付けされた各品質データの平均値を品質予測対象製品の予測品質データとしてもよい。また、予測部は、検索用テーブルから抽出された実績データの各品質データに対する重み付けを、現在時刻と処理時刻との時刻差が小さいほど大きく設定してもよい。あるいは、予測部は、検索用テーブルから、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データが複数抽出されたとき、抽出された複数の実績データのプロセス変数値および品質データに基づいて、未定係数を乗じたプロセス変数値の線形和として品質データを推定する重回帰モデルの未定係数値を算出し、重回帰モデルに品質予測対象製品のプロセス変数値を設定して取得される値を予測品質データとしてもよい。また、予測部は、検索用テーブルから、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データが複数抽出されたとき、抽出された複数の実績データのプロセス変数値および品質データに基づいて、予め設定された確率密度関数のパラメータを算出し、確率密度関数モデルに該算出されたパラメータと予め設定した確率値を設定して取得される値を予測品質データとしてもよい。さらに、予測部は、生成した予測品質データをユーザに提示する出力装置に出力してもよい。また、バラツキ指標は、確率密度関数の累積確率が予め設定された値であるときの確率変数値と確率分布の最頻値との差分であり、モデル選択部は、各局所領域における、バラツキ指標の最大値が最小である場合に、バラツキ指標が最小であるとしてもよい。あるいは、バラツキ指標は、局所領域内の品質データの点数と該品質データの分散との積であり、モデル選択部は、各局所領域における、品質データの点数と該品質データの分散との積の合計値が最小である場合に、バラツキ指標が最小であるとしてもよい。さらに、バラツキ指標は、局所領域内の品質データの分散であり、モデル選択部は、各局所領域における、品質データの分散の最大値が最小である場合に、バラツキ指標が最小であるとしてもよい。また、プロセス変数値空間分割部は、プロセス変数値空間を、分割候補点においてそれぞれ2つの局所領域に分割し、確率密度関数算出部は、分割候補点で分割された局所領域のうち少なくともいずれか一方の局所領域に属するデータ点数が基準値未満となるとき、当該分割候補点を候補から除外してもよい。さらに、再分割対象領域判定部は、プロセス変数値空間の分割回数が予め設定された最大分割数に達したとき、分割処理を終了してもよい。また、プロセス変数値空間分割部は、局所領域の再分割時において、プロセス変数値空間分割部により設定された分割候補点のうち分割しようとする局所領域の領域外となる分割候補点を候補から除外してもよい。また、上記課題を解決するために、本発明の別の観点によれば、上記品質予測装置を用いて、製造途中の製品の予測品質データを生成し、生成した予測品質データに基づき製造される製品が所望の品質データとなるように製造プロセスの操業条件を決定する操業条件決定方法が提供される。上記品質予測装置を用いた操業条件決定方法は、データ抽出部により、製造プロセスにおける過去の実績データを記憶する実績データ記憶部から、既に工程が終了し実績が確定したプロセス変数値のみを確定実績データとして抽出し、検索用テーブル作成部により、確定実績データに基づいて検索用テーブルを作成し、予測部により、品質予測対象製品の予測品質データを出力するために必要な、確定実績データに含まれる操業条件以外の、決定対象である操業条件の決定されるべき値の候補である代表値を複数作成し、各代表値と確定実績データとの組合せについて予測品質データを生成して、所望の品質データに最も近い予測品質データとなった代表値を製造プロセスの操業条件として決定することを特徴とする。ここで、プロセス変数値の代表値は、各プロセス変数値に対して検索用テーブルの作成時に採用された分割候補点で区分された各区間における代表値としてもよい。例えば、プロセス変数値の代表値は、各プロセス変数値に対して検索用テーブルの作成時に採用された分割候補点で区分された各区間における中間値としてもよい。さらに、予測部により、確定実績データと、当該確定実績データに含まれるプロセス変数値以外のプロセス変数値が取り得る代表値とからなる複数の入力データを作成し、各入力データに基づいて、検索用テーブルから当該入力データに類似する実績データを取得し、取得した実績データを未確定のプロセス変数値の取り得る代表値毎に層別して、それぞれについて予測品質データを生成し、所望の品質データに最も近い予測品質データの代表値を製造プロセスの操業条件として決定してもよい。また、上記課題を解決するために、本発明の別の観点によれば、製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測方法が提供される。品質予測方法は、製造プロセスにおける過去の実績データを記憶する実績データ記憶部から実績データを抽出するデータ抽出ステップと、データ抽出ステップにより抽出された実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成ステップと、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出し、抽出された実績データに基づいて、品質予測対象製品の品質を予測する予測ステップと、を含む。そして、検索用テーブル作成ステップは、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割ステップと、分割候補点で分割されたプロセス変数値空間の各局所領域について、品質データを確率変数とする確率密度関数を用いて、局所領域における品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出ステップと、各局所領域について算出されたバラツキ指標が最小となる局所領域を生成する分割候補点を確定分割点として決定し、確定分割点によりプロセス変数値空間を確定分割するモデル選択ステップと、モデル選択ステップにて確定分割点にて確定分割されたプロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定ステップと、プロセス変数値空間の局所領域の分割結果を検索用テーブルとして出力する分割結果出力ステップと、を含んで、再分割対象領域判定ステップにより再分割が必要と判定された局所領域に対して、プロセス変数値空間分割ステップと、確率密度関数算出ステップと、モデル選択ステップと、再分割対象領域判定ステップとを繰り返し、予測ステップは、品質予測対象製品のプロセス変数値と同一の局所領域に存在するプロセス変数値を類似のプロセス変数値として、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出することを特徴とする。さらに、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測装置として機能させるためのコンピュータプログラムが提供される。かかるコンピュータプログラムは、製造プロセスにおける過去の実績データを記憶する実績データ記憶手段から実績データを抽出するデータ抽出手段と、データ抽出手段により抽出された実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成手段と、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出し、抽出された実績データに基づいて、品質予測対象製品の品質を予測する予測手段と、を備える。そして、検索用テーブル作成手段は、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割手段と、分割候補点で分割されたプロセス変数値空間の各局所領域について、品質データを確率変数とする確率密度関数を用いて、局所領域における品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出手段と、各局所領域について算出されたバラツキ指標が最小となる局所領域を生成する分割候補点を確定分割点として決定し、確定分割点によりプロセス変数値空間を確定分割するモデル選択手段と、モデル選択手段により確定分割点にて確定分割されたプロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定手段と、プロセス変数値空間の局所領域の分割結果を検索用テーブルとして出力する分割結果出力手段と、を備えて、再分割対象領域判定手段により再分割が必要と判定された局所領域に対して、プロセス変数値空間分割手段による分割候補点設定と、確率密度関数算出手段によるバラツキ指標算出と、モデル選択手段による確定分割と、再分割対象領域判定手段による再分割の要否判定とを繰り返し、予測手段は、品質予測対象製品のプロセス変数値と同一の局所領域に存在するプロセス変数値を類似のプロセス変数値として、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出することを特徴とする。また、上記課題を解決するために、本発明の別の観点によれば、コンピュータに、製造プロセスの操業条件およびプロセスの状態量からなるプロセス変数値と類似するプロセス変数値を有する製造プロセスにおける過去の実績データに基づいて、製品の品質を予測する品質予測装置として機能させるためのプログラムを記憶したコンピュータ読み取り可能な記憶媒体が提供される。かかる記憶媒体は、製造プロセスにおける過去の実績データを記憶する実績データ記憶手段から実績データを抽出するデータ抽出手段と、データ抽出手段により抽出された実績データに基づいて、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データの検索に用いる検索用テーブルを作成する検索用テーブル作成手段と、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出し、抽出された実績データに基づいて、品質予測対象製品の品質を予測する予測手段と、を備えるプログラムを記憶する。そして、記憶媒体に記憶されたプログラムの検索用テーブル作成手段は、過去の実績データのプロセス変数値からなるプロセス変数値空間をそれぞれが複数の局所領域に分割する分割候補点を複数設定するプロセス変数値空間分割手段と、分割候補点で分割されたプロセス変数値空間の各局所領域について、品質データを確率変数とする確率密度関数を用いて、局所領域における品質データのバラツキの度合いを表すバラツキ指標を算出する確率密度関数算出手段と、各局所領域について算出されたバラツキ指標が最小となる局所領域を生成する分割候補点を確定分割点として決定し、確定分割点によりプロセス変数値空間を確定分割するモデル選択手段と、モデル選択手段により確定分割点にて確定分割されたプロセス変数値空間の局所領域について、再分割の要否を判定する再分割対象領域判定手段と、プロセス変数値空間の局所領域の分割結果を検索用テーブルとして出力する分割結果出力手段と、を備えて、再分割対象領域判定手段により再分割が必要と判定された局所領域に対して、プロセス変数値空間分割手段による分割候補点設定と、確率密度関数算出手段によるバラツキ指標算出と、モデル選択手段による確定分割と、再分割対象領域判定手段による再分割の要否判定とを繰り返し、予測手段は、品質予測対象製品のプロセス変数値と同一の局所領域に存在するプロセス変数値を類似のプロセス変数値として、品質予測対象製品のプロセス変数値と類似するプロセス変数値を有する実績データを検索用テーブルから抽出することを特徴とする。以上説明したように本発明によれば、大規模、複雑、非線形かつ非定常なプロセスの実績データから類似例を高速に検索し、製品の予測品質データを高精度に生成することが可能な、品質予測装置、操業条件決定方法、品質予測方法、コンピュータプログラムおよびコンピュータ読み取り可能な記憶媒体を提供することができる。以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。<1.第1の実施形態>[品質予測装置の概略構成]まず、図1および図2を参照して、本発明の第1の実施形態に係る品質予測装置100の概略構成について説明する。なお、図1は、本実施形態に係る品質予測装置100の概略構成を示すブロック図である。図2は、本実施形態に係る検索用テーブル作成部120の構成を示すブロック図である。本実施形態に係る品質予測装置100は、製造プロセスのプロセス変数値および製造された製品の品質データである材質値に関するデータベース検索用テーブルを作成し、作成したデータベースから材質値を予測したい製品のプロセス変数値に類似するプロセス変数値を持つ過去の類似事例を検索することにより、製造される製品の材質値を予測する。本実施形態では、製造プロセスの一例として薄板の製造プロセス10を取り上げ、製品の品質予測装置100の構成およびその機能について説明する。製造プロセス10において製造される製品の品質予測を行うため、製造プロセス10より、プロセス変数値および材質値が取得される。薄板の製造プロセス10におけるプロセス変数値としては、例えば、処理時刻、試片番号、採取部位、圧延コイル番号、巻番号、注文厚み、注文巾、連続焼鈍ライン入側板厚、連続焼鈍ライン入側板巾、冷延材質コード、連続焼鈍ライン中央速度、加熱ゾーン出側板温、第1均熱ゾーン出側板温、第2均熱ゾーン出側板温、第1冷却ゾーン出側板温、再加熱ゾーン出側板温、過時効ゾーン出側板温、第2冷却ゾーン出側板温、調質圧延工程伸び率、調質圧延工程圧延力、圧延ラインにおける粗圧延最終スタンド出側温度、仕上出側温度、巻取り温度等がある。また、精錬工程を終了した時点での溶鋼内のC量、Si量、Mn量、P量、S量、Cu量、Ni量、Cr量、Mo量、Nb量、V量、Ti量、B量、Al量、N量、O量、Ca量等もプロセス変数値として用いられる。製造プロセス10では、これらの値を計測する各種センサや入力するための端末が複数設置されている。また、製造プロセス10には、材質値として、降伏強度YP、引張強度TS、延性EL、ランクフォード値R値を計測する装置、あるいはこれらの結果を入力するための端末が複数設置されている。各種センサにより計測され、あるいは端末より入力されたプロセス変数値や、装置により計測され、あるいは端末より入力された材質値は、データ処理部20に入力される。データ処理部20は、入力された各種情報を操業オペレータに提示して製造プロセスの操業状態を通知する。また、データ処理部20は、入力された各種情報のうち少なくとも1つを、後述する実績データベース30に記録する。さらに、データ処理部20は、操業オペレータによって品質を予測する対象のプロセス変数値等の情報が入力される入力装置としても機能する。データ処理部20は、具体的には、薄板の製造プロセス10における計測や制御、データ収集を行うプロセスコンピュータプロコンやビジネスコンピュータビジコンにより構成される。実績データベース30は、製造プロセス10における実績データとしてプロセス変数値および製造された製品の材質値を記憶する記憶部である。実績データベース30は、一定時刻毎、一定長さ毎、あるいは薄板製品の製造単位であるコイル毎に製造プロセスでの処理時刻と関連付けて、データ処理部20より入力されるプロセス変数値および材質値を格納している。この際、実績データベース30に格納されるデータは、周期や長さが揃えられた同一単位の情報とされている。品質予測装置100は、実績データベース30に記憶された実績データ、すなわちプロセス変数値および材質値を取得し、製品の品質を予測する装置である。品質予測装置100は、図1に示すように、データ抽出部110と、検索用テーブル作成部120と、記憶部130と、予測部140とを備える。データ抽出部110は、実績データベース30より、プロセス変数値および材質値を抽出する。実績データベース30より抽出するプロセス変数値および材質値は、予め設定することができる。具体的には、データ抽出部110は、例えば、材質を予測する製品の処理時刻から予め設定された所定時間だけ遡った過去時刻から現在時刻までのプロセス変数値および材質値を抽出する。データ抽出部110は、抽出したプロセス変数値および材質値を検索用テーブル作成部120へ出力する。検索用テーブル作成部120は、データ抽出部110より抽出されたプロセス変数値および材質値をプロセス変数値空間の木構造の分割を用いて分類した検索用テーブル132を作成する。検索用テーブル作成部120は、図2に示すように、データ入力部121と、プロセス変数値空間分割部122と、確率密度関数算出部123と、モデル選択部124と、再分割対象領域判定部125と、分割結果出力部126とからなる。データ入力部121は、データ抽出部110が実績データベース30から抽出したプロセス変数値および材質値が入力されるインタフェースである。データ入力部121は、データ抽出部110より入力されたプロセス変数値および材質値をプロセス変数値空間分割部122へ出力する。プロセス変数値空間分割部122は、データ入力部121より入力されたプロセス変数値および材質値からなるプロセス変数値空間を分割する分割候補点を複数設定する。確率密度関数算出部123は、プロセス変数値空間分割部122により作成された分割候補点について、プロセス変数値空間の仮分割を行う。そして、確率密度関数算出部123は、仮分割されたプロセス変数値空間の各局所領域について、材質値を確率変数とする確率密度関数を算出し、局所領域における品質データのバラツキの度合いを表すバラツキ指標を算出する。モデル選択部124は、分割対象の局所領域に対する分割候補点の中から、各局所領域について算出されたバラツキ指標が最小となる局所領域を生成する分割候補点を確定分割点として決定する。モデル選択部124は、確定分割点によりプロセス変数値空間を確定分割する。再分割対象領域判定部125は、モデル選択部124により決定された確定分割点においてプロセス変数値空間を分割したときに、全体の領域分割数と予め設定された最大分割数とを比較して、領域の再分割の要否を判定する。再分割対象領域判定部125により領域の再分割が必要と判断されると、プロセス変数値空間分割部122、確率密度関数算出部123およびモデル選択部124が機能され、領域の再分割が行われる。分割結果出力部126は、再分割対象領域判定部125により領域の再分割が不要であると判定されたとき、分割処理の結果を品質予測装置100の記憶部130や表示装置40へ出力する。分割結果出力部126から出力される分割処理の結果は検索用テーブル132として利用される。なお、検索用テーブル作成部120による検索用テーブル作成処理の詳細については後述する。図1の説明に戻り、記憶部130は、品質予測装置100による品質予測処理に必要な情報を記憶する。記憶部130には、例えば、分割結果出力部126より出力された分割処理の結果が検索用テーブル132として記憶されている。検索用テーブル132は、後述の予測部140による品質予測処理にて用いられる。検索用テーブル132には、例えば、図6下に示す表のように、プロセス変数値空間の各局所領域を特定するプロセス変数値の範囲が規定されている。なお、図6に示す検索用テーブル132の構成は一例であって、他の形式で構成されていてもよい。予測部140は、検索用テーブル132を用いて、予測対象に類似する類似事例を検索し、類似事例のプロセス変数値および材質値を用いて、品質を予測する対象製品の材質値予測品質データを予測する。類似事例とは、予測対象のプロセス変数値が存在する局所領域内に存在する過去の操業データである。なお、予測部140による品質予測処理の詳細については後述する。予測部140は、予測品質データを表示装置40へ出力する。表示装置40は、品質予測装置100において作成された検索用テーブル132の内容や予測品質データを操業オペレータに提示する出力装置の一例である。表示装置40としては、例えば液晶ディスプレイや有機ELディスプレイ等を用いることができる。なお、図1では、表示装置40を品質予測装置100とは別体として構成した例を示したが、本発明はかかる例に限定されず、品質予測装置100が表示装置40等の情報を提示する出力手段を備えるように構成することもできる。以上、本実施形態に係る品質予測装置100について、製造プロセス10の実績データを取得し、予測品質データを出力するシステムとともに説明した。[品質予測装置による品質予測処理]次に、図3~図9に基づいて、本実施形態に係る品質予測装置100による品質予測処理について説明する。なお、図3は、本実施形態に係る品質予測装置100による検索用テーブル作成処理を示すフローチャートである。図4は、プロセス変数値空間の領域分割決定方法を示す説明図である。図5は、プロセス変数値空間の領域再分割方法を示す説明図である。図6は、検索用テーブル132の一構成例を示す説明図である。図7は、予測部140における製品の品質予測処理を示すフローチャートである。図8は、製品の製造時刻に基づく重み係数の一設定例を示すグラフである。図9は、製品の製造時刻に基づく重み係数の他の一設定例を示すグラフである。本実施形態に係る品質予測装置100は、操業条件およびプロセスの状態量からなるプロセス変数値に基づいて、そのプロセス変数値と類似した過去の類似事例を検索し、検索の結果得られた類似事例の材質値に基づいて製品の材質値を予測する。このため、品質予測装置100は、まず、品質の予測を行う対象のプロセス変数値から、当該プロセス変数値と類似する類似事例を検索するための検索用テーブル132を作成する。そして、品質予測装置100は、検索用テーブル132を用いて、製品の品質を予測する品質予測処理を行う。以下、品質予測装置100にて行われる検索用テーブル作成処理および品質予測処理について、詳細に説明する。検索用テーブル作成処理検索用テーブル作成部120による検索用テーブル作成処理は、図3に示すように、まず、データ入力部121に、データ抽出部110が実績データベース30から抽出したプロセス変数値および材質値が入力されることにより開始されるS100。プロセス変数値データは複数p個のデータu1,u2,,upからなり、それらは連続値、離散値のいずれであってもよい。p個のプロセス変数値のデータをまとめてp次元ベクトルu=u1u2upと表現する。全体でN個のプロセス変数値データu1,u2,,uNが得られるとする。一方、材質値のデータyは、N個のプロセス変数値データu1,u2,,uNに対応して得られるものであり、各プロセス変数値データに対応する材質データをy1,y2,,yNと表す。これらの材質値データは連続値、離散値のいずれであってもよい。以下の処理において、材質値のデータyを確率変数とみなし、その確率密度関数はプロセス変数値による空間の領域に依存して決定されるものとする。データ入力部121は、プロセス変数値データおよび材質値データをプロセス変数値空間分割部122へ出力する。次いで、プロセス変数値空間分割部122は、プロセス変数値データuを構成する各データu1,u2,,upの分割候補点を作成するS102。分割候補点は、例えば、操業知見に基づき決定することもでき、また、離散データの場合にはデータの取り得る値の中間値としたり、一定間隔でデータを分割する値に設定したり、一定データ数間隔でデータを分割する値に設定したりすることができる。あるいはこれらの組合せによっても分割候補点を決定することができる。プロセス変数値データの各データu1,u2,,upに対し分割候補点がそれぞれm1,m2,,mp個であったとすると、すべての分割候補点数は、下記数式1により表される。数式1一例として、2つのプロセス変数値データを構成するデータu1,u2からなるプロセス変数値空間において分割候補点を設定した場合を考える。操業知見や取り得る値の中間値とする等の方法を用いて分割候補点を設定したとき、例えば図4に示すように、プロセス変数値のデータu1について3つの分割候補点u11,u12,u13が、データu2については2つの分割候補点u21,u22が設定されたとする。すなわち、プロセス変数値空間には5つの分割候補点が設定されている。プロセス変数値空間分割部122により分割候補点が設定されると、確率密度関数算出部123は、プロセス変数値空間を各プロセス変数値の各分割候補点で2分割するS106。なお、ステップS104の処理は、プロセス変数値空間の分割が2回目以降である場合に行われるため、かかる処理については後述する。ステップS106で行われる分割は正式な分割点確定分割点を決定するための仮の分割である。図4の例では、確率密度関数算出部123は、プロセス変数値空間を分割候補点u11、u12、u13、u21、u22でそれぞれ2分割して、分割候補1-1~1-5を作成する。さらに、確率密度関数算出部123は、各分割候補1-1~1-5について、分割して得られた2つの局所領域のいずれかに属するデータ点数が基準値以上であるか否かを判定するS108。局所領域に属するデータ点数が基準値より少ない場合には、確率密度関数をあてはめても精度が低いため、ステップS108では、確率密度関数のあてはめを行わない分割候補を外す処理を行っている。なお、基準値は、例えば10~100程度の数に設定するのがよく、全体のデータ数が多い場合にはデータ数の数パーセント程度の値に設定してもよい。ステップS108で局所領域に属するデータ点数が基準値より少ないと判定した分割候補については、確率密度関数算出部123は、確率密度関数のあてはめは行わず、ステップS116の処理を実行する。一方、局所領域に属するデータ点数が基準値以上であると判定した分割候補については、確率密度関数算出部123は、確率密度関数のあてはめを行うステップS110。なお、どの分割候補点で分割しても得られた局所領域のいずれかに属するデータ点数が基準値未満の場合は、確率密度関数算出部123は、その局所領域での分割は行わず、次にバラツキ指標が最大である局所領域を分割対象とする。ステップS110では、確率密度関数算出部123は、各分割候補について、分割された2つの局所領域それぞれにプロセス変数値ベクトルu1u2upが属する品質データyに対して確率密度関数fy;βをあてはめるS110。ここで、βは確率密度関数fのパラメータ以下、確率パラメータとも称する。である。確率密度関数によって1または複数のパラメータからなるが、ここではこれらをまとめて確率パラメータβとして表現している。確率密度関数の例としては、正規分布、対数正規分布、指数分布等がある。正規分布の確率密度関数は、下記数式2で表される。この場合、確率密度関数fの確率パラメータは、平均μおよび標準偏差σの2つである。したがって、これらの2つの確率パラメータをまとめてβ=μ、σと表現する。また、expは指数関数を表す。数式2また、対数正規分布の確率密度関数は、下記数式3で表される。この場合、確率密度関数fの確率パラメータは、対数平均μ ́および対数標準偏差σ ́の2つである。したがって、これらの2つの確率パラメータをまとめてβ=μ ́、σ ́と表現する。また、lnは自然対数を表す。数式3さらに、指数分布の確率密度関数は、下記数式4で表される。この場合、確率密度関数fの確率パラメータはλである。確率パラメータは1つであるが、他の確率密度関数と同様に、β=λと表現される。数式4なお、本実施形態では、対数正規分布の確率密度関数を用いている。確率密度関数の確率パラメータβは、対象となる領域Mjに属するプロセス変数値ベクトルui=ui1,ui2,,uipに対応する品質データ{yi|ui∈Mj}を用いて、統計確率手法により決定することができる。統計確率手法には、例えば最尤法がある。確率パラメータβは、領域Mjごとに決定されるので、βjと表現することができる。局所領域に対して確率密度関数のあてはめがされると、確率密度関数算出部123により各局所領域についてのバラツキ指標が算出され、モデル選択部124により各局所領域についてバラツキ指標の評価がなされるステップS112。モデル選択部124は、分割対象の局所領域に対する分割候補点の中からバラツキの評価指標バラツキ指標を基準として分割点を選択する。バラツキ指標は、領域に含まれる品質データのバラツキの度合いを示す指標である。バラツキ指標の例としては、例えば、領域内のデータの分散×データ点数σ2×nや、領域内のデータの分散σ2が考えられる。また、図4に示すように、確率密度関数の累積確率90の値yAと確率分布の最頻値yBとの差分VであるF-M値累積確率90の値-確率分布の最頻値を用いることも考えられる。バラツキ指標として、領域内のデータの分散×データ点数σ2×nを用いる場合は、分割して得られた2つの局所領域内のデータの分散×データ点数σ2×nの合計が最小であるものを分割点として選択するのがよい。また、バラツキ指標として、領域内のデータの分散σ2や、F-M値累積確率90の値-確率分布の最頻値を用いる場合は、分割して得られた2つの局所領域のバラツキ指標の大きい方が最小となる分割候補点を分割点として選択するのがよい。確率密度関数算出部123は、このような手法により各局所領域のバラツキ指標を算出する。そして、モデル選択部124は、各局所領域について算出されたバラツキ指標のうち、バラツキ指標が大きい方の領域について、バラツキ指標が最小であった局所領域を生成する分割候補点を確定分割点として決定する。例えば、図4の例では、分割候補1-1~1-5において、仮分割された2つの局所領域のバラツキ指標がそれぞれ算出される。そして、確率密度関数算出部123は、各分割候補1-1~1-5において、バラツキ指標が大きい方の領域を特定する。分割候補1-1においては、左側の局所領域がバラツキ指標の大きい方の領域となる。各分割候補1-1~1-5についてバラツキ指標の大きい方の領域を特定すると、確率密度関数算出部123は、これらのうちバラツキ指標が最小である分割候補の分割候補点を確定分割点として決定する。なお、本実施形態では、ステップS104~S114の処理を、分割候補毎に実施している。このため、ステップS112では、1つの分割候補についてのバラツキ指標が算出されると、それ以前に算出された分割候補のバラツキ指標のうち最小のものと比較し、バラツキ指標が小さい方の分割候補点でのデータを分割候補点データとして保存しておく。このとき、分割候補点、確率パラメータ、F‐M値等が分割データとして保存される。そして、すべての分割候補点についての判定を終えたとき、バラツキ指標が最小であった分割候補点が確定分割点として採用される。バラツキ指標がそれ以前に算出されたバラツキ指標の最小値より大きい場合には、当該分割候補点における分割データは保存せず、ステップS116の処理へ進む。一方、バラツキ指標がそれ以前に算出されたバラツキ指標の最小値以下である場合には、当該分割データをメモリ図示せず。に保存するS114。そして、モデル選択部124は、全プロセス変数値の全分割候補点についてステップS104~S114の処理を実行したか否かを判定するS116。全プロセス変数値の全分割候補点についてステップS104~S114の処理が終了していなければ、ステップS104に戻り、処理を繰り返す。一方、全プロセス変数値の全分割候補点についてステップS104~S114の処理が終了している場合には、確定分割点において局所領域を確定分割する領域分割処理を実行するS118。ステップS118にて領域分割処理が実行されると、再分割対象領域判定部125は、全体の領域分割数を1つ増やし、現在の領域分割数の値を更新するS120。そして、再分割対象領域判定部125は、更新された領域分割数が最大分割数に達したか否かを判定するS122。かかる処理は、全体の領域分割数があらかじめ指定した最大分割数に満たない場合に、分割して得られた局所領域をさらに分割して、予測精度を上げることを指向するために行われる。ステップS122にて、全体の領域分割数が最大分割数に達していないと判定された場合、再分割対象領域判定部125は、分割して得られた局所領域のうち、バラツキの評価指標が最大である局所領域を分割対象として決定するS124。ただし、その領域に属するデータ点数が基準値未満の場合は分割対象とせず、次にバラツキ指標が最大である局所領域を分割対象とする。なお、すべての領域について、各領域に属するデータ点数が基準値未満となった場合には、全体の領域分割数が最大分割数に達していない場合でも当該処理は終了し、ステップS126の処理が行われる。再分割対象とする局所領域が決定されると、検索用テーブル作成部120は、ステップS104からの処理を繰り返す。最大分割数は初期値を10~20程度として、得られた結果と予測精度とを見ながら順次増加させるのが適当である。ステップS124にて分割対象とされた局所領域に対するステップS104~S120の分割処理は、プロセス変数値空間の2回目以降の分割となる。プロセス変数値空間の分割が2回目以降となる場合、すなわち一旦分割された局所領域を再分割する場合、プロセス変数値空間分割部122は、領域外の分割候補点は除外するS104。このため、初回の分割候補点はM個であるが、再分割時はそれより少なくなる。例えば、図5において、直前の回でプロセス変数値u1の値u13が確定分割点として決定されたとき、プロセス変数値空間は局所領域A1、A2の2つに確定分割される。そして、全体の領域分割数が1増加されるとともに、生成された2つの局所領域それぞれについて確率密度関数が算出され、バラツキ指標が算出される。次の領域の確定分割点を決定する際には、局所領域A1、A2のうちバラツキ指標の大きい方についてステップS104~S120の処理が実行される。なお、図5の例では、バラツキ指標の大きい局所領域A1が次に分割される。一方、ステップS122にて、全体の領域分割数が指定した最大分割数に達したと判定した場合は、再分割対対象領域判定部125は分割処理を終了する。そして、分割結果出力部126より、ステップS102~S120により実行された分割結果を検索用テーブル132として、品質予測装置100の記憶部130に出力し記録するS126。この検索用テーブル132は、プロセス変数値による空間を局所領域に分割する方法を定義した情報である。分割結果出力部126は、表示装置40に分割結果を出力して、検索用テーブル132の内容を操業オペレータに提示するようにしてもよい。以上、本実施形態に係る品質予測装置100の検索用テーブル作成部120による検索用テーブル作成処理について説明した。かかる検索用テーブル作成処理により、データ抽出部110により実績データベース30から抽出された実績データのプロセス変数値からなるプロセス変数値空間は、例えば図6上に示すように、複数の局所領域例えば局所領域1~6に分割される。記憶部130の検索用テーブル132には、例えば図6下の表に示すように、各局所領域について、プロセス変数値空間を分割する範囲が規定されている。このように検索用テーブル132を作成することにより、プロセス変数値空間を分割する過程で操業データが層別される。したがって、本実施形態に係る品質予測装置100を用いた検索用テーブル132の作成では、従来予めモデルの入力変数をプロセスとの相関が高いものだけに絞り込む必要があったのに対し、入力変数中にプロセスとの相関の低い変数が含まれている場合であってもそのプロセス変数値は空間分割に採用されないため、モデルの精度に影響することがなく、また、予め入力変数をプロセスと相関の高いものだけに絞り込む必要がない。さらに、本実施形態に係る品質予測装置100を用いた検索用テーブル132の作成では、プロセス変数値空間を分割する際に、モデルの予測精度がよくなるように、プロセス変数値とその分割位置とが逐次検索される。これにより、プロセス変数値の数が増加しても空間分割数を適切に抑えることができる。品質予測処理検索用テーブル作成部120により検索用テーブル132が作成されると、予測部140は、検索用テーブル132を参照して、品質を予測する対象のプロセス変数値と類似する過去の類似事例の検索を行う。品質予測対象製品のプロセス変数値は、例えば、データ処理部20より品質予測装置100の予測部140へ入力されるS200。そして、予測部140は、データ処理部20から入力されたプロセス変数値に基づき、記憶部130内の検索用テーブル132を検索するための入力値を作成するS202。検索用テーブル132への入力値が作成されると、この入力値を検索キーとして、検索用テーブル132より類似事例の検索が行われるS204。ステップS204では、プロセス変数値空間において、予測したい製品のプロセス変数値データと同じ局所領域に存在する類似事例のプロセス変数値データおよび製品の材質値が抽出される。そして、予測部140は、検索用テーブル132の検索の結果得られた類似事例のプロセス変数値、製品の材質値、製造プロセス10での処理時刻あるいは製造プロセス10での処理時刻に対応する実績データベースの格納番号を用いて、予測対象の製品の材質値を予測するS206。なお、格納番号は、処理時刻に対応して時系列で付されているものとする。ステップS206での予測対象の製品の材質値の予測は、例えば、ステップS204で取得された複数の類似事例の材質値の平均値を、予測対象の製品の材質値とすることにより行うことができる。また、製造プロセス10においては、設備の経時変化等によりプロセス変数値および製品の材質値の特性が変化することがある。この場合、ステップS204で取得された複数の類似事例に対して、製造プロセス10での処理時間に基づき重み付けをして材質値の平均値を算出するようにすることで、予測値の精度を高めるようにすることができる。具体的には、まず、予測部140は、検索で得られたすべての類似事例に対して、その処理時刻と現在時刻との差で定義される時刻差tiを算出する。次いで、予測部140は、時刻差tiに基づいて重み係数Wiを算出し、重み係数Wiを用いて重み付き平均値を算出する。ここで、予測部140は、例えば図8に示すように重み付けを設定することができる。すなわち、例えば図8に示す、下記数式5で表される重み係数の分布を用いて、現在時刻に近づくほど指数的に大きくなる重みを類似事例の材質値に設定することができる。数式5において、expは指数関数、λは予め設定された指数関数の正値パラメータ、Mは類似事例の個数を表す。指数関数の正値パラメータλは、時刻差に基づく重みの度合いを調整するためのパラメータであり、値が大きいほど最新の類似事例に高い重みが設定される。このパラメータは、製造プロセスに関する知識から経験的に決定してもよく、あるいは予測値が最も実績値に一致するよう調整して決定してもよい。このような重み付けを行うことにより、設備の経時変化による特性の変化を反映した予測値を算出することができる。そして、予測部140は、検索で得られたすべての類似事例の製品の材質値yiと、これに対する重み係数Wiとより、数式6から重み付き平均値ypを取得する。数式5数式6なお、時刻差tiに基づいて重み係数を算出する式は、上記数式5に限定されるものではなく、例えば、図9に示す線形関数等、現在時刻において最も高い重みを有する任意の単調関数の式を用いてもよい。また、重み係数は、図8、図9に示したように関数式より算出する方法以外にも、例えば時刻差の範囲と重み係数を対応表にした重み係数テーブルを予め作成し、時刻差の値から重み係数テーブルを参照して重み係数を決定することもできる。なお、過去事例に対する重みを設定して平均値を算出する方法としては、予測対象製品のプロセス変数値と過去事例のプロセス変数値に基づいて、重み係数を決定する方法を用いてもよい。具体的には、予測したい事例のプロセス変数値と類似事例のプロセス変数値の差分量から重み係数を算出し、材質値に乗じた上で加算することで得られる差分量重み付き平均値を製品の予測値とするようにしてもよい。上記の製造プロセスの処理時刻や、あるいはプロセス変数値に基づいて重み係数を算出し、重み付き平均値を算出して予測値とする方法は、単独で適用してもよく、またそれぞれの手法を組み合わせて時刻差とプロセス変数値の両者を考慮した重み係数を求めて、予測値を算出してもよい。また、具体的な予測値の他の算出方法としては、検索の結果得られた複数個の類似事例のプロセス変数、および材質値の実績データから重回帰モデルを作成し、予測対象のプロセス変数値を、この重回帰モデルに入力して材質の予測値を求める方法がある。重回帰モデルの未定係数は、類似事例のプロセス変数と材質の実績値より、最小二乗法等の最適化手法で算出し、この未定係数値に予測したい時点のプロセス変数値を乗じて和を取った線形和を算出して、予測値とする。さらに、他の具体的な予測値の算出方法としては、検索の結果得られた複数個の類似事例の材質値に基づいて、予め設定された関数形の確率密度関数のパラメータを算出し、この確率密度関数モデルに基づいて予測値を算出する方法がある。一般にプロセス変数値と材質値の間には、測定ノイズ等の外乱や、センサで測定できない等の理由で時系列データベースに収集されていない因子の影響により、ほぼ同一のプロセス変数値であるにもかかわらず材質値が互いに異なるという、いわゆるバラツキが存在する。特に鉄鋼プロセスのような、高温環境等の理由で測定センサの設置に制約が多く、かつ外乱ノイズの多い製造プロセスにおいては、相当量のバラツキが存在する。このようなバラツキの大きい材質値の予測においては、顧客との取り決めで設定された材質の保証範囲を予測対象製品が満足しているかが重要であり、平均的な材質値を予測するよりも、ある所定の確率における材質値の範囲を予測することが有用であり、過去事例に基づく材質値の確率密度分布による予測が適切である。予測値の算出に用いる確率密度関数は、予め関数式を設定しておく必要があり、対象の製造プロセスに関する知識や、あるいは材質値のデータの度数分布に基づいて、分布を適切に表現できる関数を選定すればよい。具体的には、正規分布関数、指数分布関数、ポアソン分布関数、対数正規分布関数等統計学の研究分野で提案された関数を用いるのが一般的であるが、本発明は、これらの関数に限定されるものではなく、任意の関数系を用いることができる。これら確率密度関数のパラメータ例えば正規分布関数の場合は、平均値と標準偏差は、検索の結果得られた複数個の類似事例の材質値より算出される。この算出されたパラメータと、予め設定した確率値例えば材質値の平均値を中心とする確率90を、上記設定した確率密度関数の累積確率密度逆関数に代入して演算する処理によって、製品の材質値の範囲を算出し、予測値とすることができる。図10に、確率密度関数に基づく材質の予測値の算出する処理を模式的に表した図を示す。以上、本発明の第1の実施形態に係る品質予測装置100の構成と、これによる品質予測処理について説明した。本実施形態によれば、大規模、複雑、非線形かつ非定常な製造プロセスの実施データより類似例を高速に検索することができ、大規模、複雑、非線形かつ非定常な要因が影響を与える製品の材質の予測値を常に高い精度に維持することができる。また、操業並びに材質の安定化や不良発生頻度の低減に大きく寄与することができる。<2.第2の実施形態>次に、図11および図12に基づいて、本発明の第2の実施形態に係る品質予測装置による品質予測方法について説明する。本実施形態に係る品質予測装置は、材質値を予測したい製品の製造の途中工程において、望ましい材質を得るための操業条件を決定する。本実施形態の品質予測装置も第1の実施形態に係る品質予測装置100と同様に構成することができるので、ここでは装置構成の詳細な説明は省略し、製品の製造の途中工程における操業条件の決定方法を説明する。なお、図11は、本実施形態に係る品質予測装置により予測する操業条件を説明するための説明図である。図12は、本実施形態に係る品質予測装置による操業条件決定処理を示すフローチャートである。また、装置の各機能部については、以下の実施形態においても第1の実施形態と同一の符号を用いて説明する。[操業条件を決定する状況の説明]本実施形態に係る品質予測装置100では、製造プロセス10において、順次行われる工程の途中で、望ましい材質の製品を得るための製造条件を決定することができる。例えば、図11に示すように、製造プロセス10で行われる各工程を、実施順に、プロセス1、プロセス2、プロセス3、と表す。ここで、プロセス1は操業条件x1、プロセス2は操業条件x2で操業され、プロセス2に続いて行われるプロセス3は未実施であり、プロセス3の操業条件は未確定であるとする。このとき、本実施形態に係る品質予測装置100は、所望の材質値を有する製品を製造するためのプロセス3の操業条件を決定する。[製造条件決定処理]品質予測装置100は、図12に示すように、まず、データ入力部110より材質値を予測する製品のプロセス変数値をデータ処理部20から取得するS300。この際、すでに製造プロセス10における処理が完了し、実績が確定したプロセス変数に関する値のみデータ処理部20より入力される。次いで、予測部140は、実績が確定したプロセス変数値について、検索用テーブル132を検索するための入力値を作成するS302。さらに、予測部140は、材質予測に必要なプロセス変数値の中で、未だ処理が確定していないためにステップS300では取得されなかったプロセス変数値について、該プロセス変数値が取り得るすべての操業範囲を網羅するようにプロセス変数値の代表値を作成するS304。この代表値は、望ましい材質を得るために決定される操業条件の候補となる値であり、予測部140は、これらの中から最も望ましい材質を得ると予測された操業条件を最終的な操業条件として決定する。予測部140は、各プロセス変数値について、検索用テーブル132の作成時に採用された分割候補点で区切られたプロセス変数値の各区間について1つずつ代表値を決定する。例えば、図6上に示すプロセス変数空間において、データu2が未確定であるとき、予測部140は、分割候補点u20、u21、u22、u23の各区間について1つずつ代表値u21、u22、u23を決定する。このとき、代表値は、例えば、各区間の中間値とすることができる。このとき、さらに取得されていないプロセス変数値が複数ある場合は、予測部140は、それぞれのプロセス変数値の代表値のすべての組み合わせを作成する。その後、予測部140は、ステップS302で作成したプロセス変数値の実績についての入力値と、ステップS304で作成した未確定のプロセス変数値が取り得る代表値とに基づいて、検索用テーブル132を参照するための入力データを作成するS306。ステップS306において、入力データは、ステップS304で作成されたプロセス変数値の代表値のすべての組み合わせの数だけ作成される。これらの入力データは、予測対象の製品が今後、取り得る操業条件をすべて網羅したものに相当する。次いで、ステップS306で作成した入力データのうちの1つについて、記憶部130内の検索用テーブル132の検索を行い、入力データに類似した類似事例を選択する処理を行うS308。ステップS308の処理は、具体的には、図7のステップS204における類似事例の選択方法と同様に行うことができる。ステップS308の処理はステップS306で作成された入力データのすべてについて行われ、各入力データに類似する類似事例が抽出され、抽出された類似事例のデータからなる集合を作成する。このとき、異なる入力データに対して、同じ類似事例が抽出されることもある。この場合、類似事例のデータからなる集合では同じ類似事例のデータは1つしか含まないようにする。さらに、予測部140は、ステップS308で作成された類似事例のデータからなる集合を、ステップS304で作成した未確定のプロセス変数値の取り得る代表値毎に層別し、その層別毎に予測したい製品の材質値を算出する予測演算を行うS310。ステップS310の具体的な処理としては、例えば図7のステップS206の処理とすることができる。ステップS310の処理は、ステップS304で作成した未確定のプロセス変数値の取り得る代表値のすべてに対して行われるS312。これは、予測対象製品が今後取り得る操業条件の全てに対して、各々の材質予測値を算出することに相当する。ステップS312ですべての入力データの予測演算が完了していないと判定された場合には、ステップS310の処理を繰り返し行う。一方、ステップS312ですべての入力データの予測演算が完了したと判定された場合には、後述のステップS314の処理を行う。すべての入力データの予測演算が完了すると、予測部140は、ステップS300~S312までの処理で算出した各操業条件に対する材質の予測値を、予め設定した評価関数に基づいて評価する。そして、予測部140は、最も良好な材質予測値となった入力データを選択し、その入力データに相当するプロセス変数値の代表値に換算して、換算した値を出力するS314。以上、本実施形態に係る品質予測装置100による操業条件決定処理について説明した。本実施形態によれば、大規模、複雑、非線形かつ非定常な製造プロセスの実施データより類似例を高速に検索することができ、製造の途中工程において、望ましい材質を得るための操業条件を決定することができる。これにより、操業並びに材質の安定化や不良発生頻度の低減に大きく寄与することもできる。[実施例]以下、図13~図15に、製造プロセスとして、鉄鋼の薄板製造プロセスを対象とし、材質値を予測した例を示す。本実施例では、品質予測装置を薄板の製造プロセスに対して適用した例について説明するが、本発明はかかる例に限定されない。品質予測装置100は、操業条件およびプロセスの状態量からなるプロセス変数値と、製造された製品の材質値を取得可能な製造プロセスに対して適用することができる。薄板の製造プロセスに適用する場合、プロセス変数値としては、例えば、処理時刻、試片番号、採取部位、圧延コイル番号、巻番号、注文厚み、注文巾、連続焼鈍ライン入側板厚、連続焼鈍ライン入側板巾、冷延材質コード、連続焼鈍ライン中央速度、加熱ゾーン出側板温、第1均熱ゾーン出側板温、第2均熱ゾーン出側板温、第1冷却ゾーン出側板温、再加熱ゾーン出側板温、過時効ゾーン出側板温、第2冷却ゾーン出側板温、調質圧延工程伸び率、調質圧延工程圧延力、圧延ラインにおける粗圧延最終スタンド出側温度、仕上出側温度、巻取り温度等とすることができる。また、精錬工程を終了した時点での溶鋼内のC量、Si量、Mn量、P量、S量、Cu量、Ni量、Cr量、Mo量、Nb量、V量、Ti量、B量、Al量、N量、O量、Ca量等もプロセス変数値として用いられる。また、薄板の材質値として、降伏強度YP、引張強度TS、延性EL、ランクフォード値R値等を用いることができる。これらの値は、製造プロセスに設けられた各種センサや計測装置、情報を入力するための端末から入力される。図13は、材料の成分値のうち材質値に影響の大きい2つの成分値と、下流工程の温度条件の合計3因子をプロセス変数値として、上記特許文献1の手法を用いて材質値予測を行ったときの予測値と実績値との関係を示すグラフである。図14は、図13で用いた3因子に、上流工程での操業条件を加えて、合計7因子をプロセス変数値として、上記特許文献1の手法を用いて材質値予測を行ったときの予測値と実績値との関係を示すグラフである。図15は、図14と同じ7因子をプロセス変数値として、本実施形態に係る品質予測装置100による品質予測方法を用いて材質値予測を行ったときの予測値と実績値との関係を示すグラフである。本例では、値が小さくなるほど材質が良好であるように材質値を正規化している。製造プロセスにおいては、材料の成分値のバラツキや、上流工程での製造条件のバラツキが、下流工程の材質値に影響するため、過去の操業実績データとこれから製造する材料の成分値や上流工程の製造条件をもとに、下流工程の操業条件を調整している。図13と図14とを比較すると、上流工程での操業条件を加えることで、予測精度の向上を狙ったものであるが、かえって精度が悪くなっていることがわかる。一方、図14と図15とを比較すると、プロセス変数値を増やすことによって、予測精度が向上していることがわかる。これより、本実施形態に係る品質予測装置100による品質予測の精度が従来の品質予測手法と比較して高いことが確認された。[最適操業条件ガイダンスへの応用]また、本発明の第2の実施形態に係る品質予測装置100を用いて、最適操業条件を操業オペレータにガイダンスすることも可能である。図16に、品質予測装置100を用いて最適操業条件ガイダンスを行う場合の一処理例を示す。図16では、品質予測装置100を用いてコイル毎の最適操業条件温度条件をガイダンスする処理を示し、温度条件を変化させたときの品質指標を最良とする温度条件がガイダンスされる。本例では、累積確率80となる材質値を品質指標としており、これが最小となる温度条件を最適操業条件とする。図16に示すように、まず、予測対象コイルのプロセス変数値のデータが入力されるとS400、品質予測装置100は、評価用データ1本1本について、温度条件を変えて、検索用テーブル132より該当する過去の実績データを抽出するS402。すべての温度条件について実績データの抽出を終えると、品質予測装置100は、抽出された過去の実績データを操業条件毎に、ここでは温度条件毎に層別するS404。図17に示す例では、実績データを温度条件について水準1~5に層別している。そして、品質予測装置100は、操業条件別に確率分布モデルを作成し、品質指標を算出するS406。本例では、温度条件の層別それぞれについて確率分布モデルが作成され、累積確率80となる材質値が求められる。さらに、品質予測装置100は、品質指標が最適となる操業条件を決定するS408。本例では、品質指標である累積確率80となる材質値が最小となる温度条件が最適操業条件として決定される。その後、品質予測装置100は、ステップS402で抽出されステップS404で層別された過去の類似事例や、ステップS408で決定された最適操業条件を提示するS410。類似事例や最適操業条件は、例えば図17に示すように複合グラフとして図示することができる。図17では、温度条件毎の材質値を散布図として示すことで、各層に分布されたコイルの品質を目視することができる。また、温度条件毎の過去の実績データ数も棒グラフにより表してもよい。品質予測装置100は、累積確率80となる材質値が最小となる温度条件図17では水準4を最適操業条件として操業オペレータにガイダンスする。このように、品質予測装置100を用いて、最適操業条件をガイダンスする機能を構成することも可能である。<3.品質予測装置のハードウェア構成例>次に、図18を参照しながら、本発明の実施形態に係る品質予測装置100のハードウェア構成について、詳細に説明する。図18は、本発明の実施形態に係る品質予測装置100のハードウェア構成例を示すブロック図である。品質予測装置100は、主に、CPU901と、ROM903と、RAM905と、を備える。また、品質予測装置100は、更に、バス907と、入力装置909と、出力装置911と、ストレージ装置913と、ドライブ915と、接続ポート917と、通信装置919とを備える。CPU901は、演算処理装置および制御装置として機能し、ROM903、RAM905、ストレージ装置913、またはリムーバブル記録媒体921に記録された各種プログラムに従って、品質予測装置100内の動作全般またはその一部を制御する。ROM903は、CPU901が使用するプログラムや演算パラメータ等を記憶する。RAM905は、CPU901が使用するプログラムや、プログラムの実行において適宜変化するパラメータ等を一次記憶する。これらはCPUバス等の内部バスにより構成されるバス907により相互に接続されている。バス907は、ブリッジを介して、PCIPeripheralComponentInterconnect/Interfaceバスなどの外部バスに接続されている。入力装置909は、例えば、マウス、キーボード、タッチパネル、ボタン、スイッチおよびレバーなどユーザが操作する操作手段である。また、入力装置909は、例えば、赤外線やその他の電波を利用したリモートコントロール手段いわゆる、リモコンであってもよいし、品質予測装置100の操作に対応したPDA等の外部接続機器923であってもよい。さらに、入力装置909は、例えば、上記の操作手段を用いてユーザにより入力された情報に基づいて入力信号を生成し、CPU901に出力する入力制御回路などから構成されている。品質予測装置100のユーザは、この入力装置909を操作することにより、品質予測装置100に対して各種のデータを入力したり処理動作を指示したりすることができる。出力装置911は、取得した情報をユーザに対して視覚的または聴覚的に通知することが可能な装置で構成される。このような装置として、CRTディスプレイ装置、液晶ディスプレイ装置、プラズマディスプレイ装置、ELディスプレイ装置およびランプなどの表示装置や、スピーカおよびヘッドホンなどの音声出力装置や、プリンタ装置、携帯電話、ファクシミリなどがある。出力装置911は、例えば、品質予測装置100が行った各種処理により得られた結果を出力する。具体的には、表示装置は、品質予測装置100が行った各種処理により得られた結果を、テキストまたはイメージで表示する。他方、音声出力装置は、再生された音声データや音響データ等からなるオーディオ信号をアナログ信号に変換して出力する。なお、本発明の実施形態では、出力装置911は、品質予測装置100とは別体の表示装置40として設けられている。ストレージ装置913は、品質予測装置100の記憶部の一例として構成されたデータ格納用の装置である。ストレージ装置913は、例えば、HDDHardDiskDrive等の磁気記憶部デバイス、半導体記憶デバイス、光記憶デバイス、または光磁気記憶デバイス等により構成される。このストレージ装置913は、CPU901が実行するプログラムや各種データ、および外部から取得した各種のデータなどを格納する。ドライブ915は、記録媒体用リーダライタであり、品質予測装置100に内蔵、あるいは外付けされる。ドライブ915は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体921に記録されている情報を読み出して、RAM905に出力する。また、ドライブ915は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記録媒体921に記録を書き込むことも可能である。リムーバブル記録媒体921は、例えば、CDメディア、DVDメディア、Blu-rayメディア等である。また、リムーバブル記録媒体921は、コンパクトフラッシュ登録商標CompactFlash:CF、フラッシュメモリ、または、SDメモリカードSecureDigitalmemorycard等であってもよい。また、リムーバブル記録媒体921は、例えば、非接触型ICチップを搭載したICカードIntegratedCircuitcardまたは電子機器等であってもよい。接続ポート917は、機器を品質予測装置100に直接接続するためのポートである。接続ポート917の一例として、USBUniversalSerialBusポート、IEEE1394ポート、SCSISmallComputerSystemInterfaceポート、RS-232Cポート等がある。この接続ポート917に外部接続機器923を接続することで、品質予測装置100は、外部接続機器923から直接各種のデータを取得したり、外部接続機器923に各種のデータを提供したりする。通信装置919は、例えば、通信網925に接続するための通信デバイス等で構成された通信インタフェースである。通信装置919は、例えば、有線または無線LANLocalAreaNetwork、Bluetooth登録商標、またはWUSBWirelessUSB用の通信カード等である。また、通信装置919は、光通信用のルータ、ADSLAsymmetricDigitalSubscriberLine用のルータ、または、各種通信用のモデム等であってもよい。この通信装置919は、例えば、インターネットや他の通信機器との間で、例えばTCP/IP等の所定のプロトコルに則して信号等を送受信することができる。また、通信装置919に接続される通信網925は、有線または無線によって接続されたネットワーク等により構成され、例えば、インターネット、家庭内LAN、赤外線通信、ラジオ波通信または衛星通信等であってもよい。以上、本発明の実施形態に係る品質予測装置100の機能を実現可能なハードウェア構成の一例を示した。上記の各構成要素は、汎用的な部材を用いて構成されていてもよいし、各構成要素の機能に特化したハードウェアにより構成されていてもよい。従って、本実施形態を実施する時々の技術レベルに応じて、適宜、利用するハードウェア構成を変更することが可能である。以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。"