特許 6844204 予測テーブル設計装置、方法及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6844204

(24)【登録日】2021年3月1日

(45)【発行日】2021年3月17日

(54)【発明の名称】予測テーブル設計装置、方法及びプログラム

(51)【国際特許分類】

   G06Q 10/04 20120101AFI20210308BHJP

   G06Q 50/04 20120101ALI20210308BHJP

   G06N 7/00 20060101ALI20210308BHJP

【ＦＩ】

   G06Q10/04

   G06Q50/04

   G06N7/00

【請求項の数】19

【全頁数】21

(21)【出願番号】特願2016-219111(P2016-219111)

(22)【出願日】2016年11月9日

(65)【公開番号】特開2017-174390(P2017-174390A)

(43)【公開日】2017年9月28日

【審査請求日】2019年7月3日

(31)【優先権主張番号】特願2016-56964(P2016-56964)

(32)【優先日】2016年3月22日

(33)【優先権主張国】JP

【新規性喪失の例外の表示】特許法第３０条第２項適用 平成２８年３月１日 一般社団法人日本鉄鋼協会「材料とプロセスＣＤ−ＲＯＭ（２０１６）Ｖｏｌ．２９ Ｎｏ．１」に発表

(73)【特許権者】

【識別番号】000006655

【氏名又は名称】日本製鉄株式会社

(74)【代理人】

【識別番号】100090273

【弁理士】

【氏名又は名称】國分 孝悦

(72)【発明者】

【氏名】塩谷 政典

【審査官】 梅岡 信幸

(56)【参考文献】

【文献】 特開２０１６−０４５７９９（ＪＰ，Ａ）

【文献】 特開２０１３−０８０４５８（ＪＰ，Ａ）

【文献】 特開２０１１−０８５９７０（ＪＰ，Ａ）

【文献】 特開２０１４−２３５４８１（ＪＰ，Ａ）

【文献】 米国特許第０８８２５３６１（ＵＳ，Ｂ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０６Ｑ １０／００−９９／００

Ｇ０６Ｎ ５／００− ７／０６

(57)【特許請求の範囲】

【請求項1】

複数の区分に分割された説明変数の各々の前記区分に対する予測対象の予測値がそれぞれ記述される予測テーブルを設計する予測テーブル設計装置であって、
前記説明変数から目的変数である前記予測対象の予測値を算出する予測モデルに、前記説明変数の各々の前記区分の代表値をそれぞれ入力して、各々の前記区分に対する前記予測対象の予測値を算出することで、前記説明変数の実績値に応じた前記予測対象の前記目的変数の実績値を予測するための前記予測テーブルを生成することを特徴とする予測テーブル設計装置。

【請求項2】

説明変数を離散化して形成される各区分に予測対象の予測値が記述される予測テーブルを設計する予測テーブル設計装置であって、
説明変数から目的変数である前記予測対象の予測値を算出する予測モデルを記憶する記憶手段と、
説明変数の実績値及び目的変数の実績値を格納する格納手段と、
説明変数を離散化する離散値配列の複数の候補について、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して生成される予測テーブルを、前記格納手段に格納する説明変数の実績値及び目的変数の実績値を用いて評価して、前記離散値配列の複数の候補から離散値配列を選択する評価手段とを備えたことを特徴とする予測テーブル設計装置。

【請求項3】

離散値配列の初期値を定め、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して、予測テーブルの初期値を生成する予測テーブル初期値生成手段と、
前記離散値配列の初期値から一又は複数の離散値を削除した離散値配列の複数の候補を作成する離散値配列候補作成手段とを備え、
前記評価手段は、前記離散値配列候補作成手段で作成した離散値配列の複数の候補について、評価を行い、離散値配列を選択することを特徴とする請求項２に記載の予測テーブル設計装置。

【請求項4】

前記評価手段で選択した離散値配列から生成される予測テーブルが、所定の条件を満たすか否かを判定する判定手段を備え、
前記判定手段で前記所定の条件を満たさないと判定した場合、
前記離散値配列候補作成手段は、前記評価手段で選択した離散値配列から一又は複数の離散値を削除した離散値配列の複数の候補を作成し、
前記評価手段は、前記離散値配列候補作成手段で作成した離散値配列の複数の候補について、評価を行い、離散値配列を選択し、
これを前記所定の条件を満たすまで繰り返すことを特徴とする請求項３に記載の予測テーブル設計装置。

【請求項5】

前記所定の条件として、予測テーブルのサイズに対する目標サイズ、又は、前記評価手段での評価結果に対する目標値が設定されることを特徴とする請求項４に記載の予測テーブル設計装置。

【請求項6】

前記評価手段で選択した離散値配列から一又は複数の離散値を削除し、かつ、前記評価手段で前記離散値配列を選択するまでに削除済みの離散値の中から、当該削除した離散値と同じ数の離散値を復活させて、改善用離散値配列候補を作成する改善用離散値配列候補作成手段と、
前記改善用離散値配列候補作成手段で作成した改善用離散値配列候補について、評価を行い、前記評価手段で選択した離散値配列よりも評価の良い改善用離散値配列候補があれば、それを前記評価手段で選択した離散値配列と置き換える改善用評価手段とを備えたことを特徴とする請求項４又は５に記載の予測テーブル設計装置。

【請求項7】

前記評価手段で離散値配列を選択した後、所定の実施条件を満たすときに、前記改善用離散値配列候補作成手段及び前記改善用評価手段による処理を実施することを特徴とする請求項６に記載の予測テーブル設計装置。

【請求項8】

前記所定の実施条件として、予測テーブルのサイズに対する条件、又は前記評価手段での評価結果に対する条件が設定されることを特徴とする請求項７に記載の予測テーブル設計装置。

【請求項9】

前記離散値配列候補作成手段は、削除する離散値を乱数を用いて決めることを特徴とする請求項３乃至８のいずれか１項に記載の予測テーブル設計装置。

【請求項10】

前記改善用離散値配列候補作成手段は、削除する離散値を乱数を用いて決めることを特徴とする請求項６乃至８のいずれか１項に記載の予測テーブル設計装置。

【請求項11】

説明変数の離散化形式は、離散値と離散値の間の予測値は一定値とする形式、及び、離散値と離散値の間の予測値はこれら離散値に対応する予測値を結ぶ直線で表わす形式のうち少なくともいずれか一方であることを特徴とする請求項２乃至１０のいずれか１項に記載の予測テーブル設計装置。

【請求項12】

前記評価手段は、予測テーブルの予測値と、目的変数の実績値との誤差を評価することを特徴とする請求項２乃至１１のいずれか１項に記載の予測テーブル設計装置。

【請求項13】

前記評価手段は、前記誤差の平方和を評価値とする、又は、前記誤差の絶対値和を評価値とすることを特徴とする請求項１２に記載の予測テーブル設計装置。

【請求項14】

前記改善用評価手段は、予測テーブルの予測値と、目的変数の実績値との誤差を評価することを特徴とする請求項６乃至８のいずれか１項に記載の予測テーブル設計装置。

【請求項15】

前記改善用評価手段は、前記誤差の平方和を評価値とする、又は、前記誤差の絶対値和を評価値とすることを特徴とする請求項１４に記載の予測テーブル設計装置。

【請求項16】

複数の区分に分割された説明変数の各々の前記区分に対する予測対象の予測値がそれぞれ記述される予測テーブルを設計する予測テーブル設計方法であって、
コンピュータが、前記説明変数から目的変数である前記予測対象の予測値を算出する予測モデルに、前記説明変数の各々の前記区分の代表値をそれぞれ入力して、各々の前記区分に対する前記予測対象の予測値を算出することで、前記説明変数の実績値に応じた前記予測対象の前記目的変数の実績値を予測するための前記予測テーブルを生成することを特徴とする予測テーブル設計方法。

【請求項17】

説明変数を離散化して形成される各区分に予測対象の予測値が記述される予測テーブルを設計する予測テーブル設計方法であって、
コンピュータが、
説明変数から目的変数である前記予測対象の予測値を算出する予測モデルを記憶する記憶手段と、説明変数の実績値及び目的変数の実績値を格納する格納手段とを用いて、
説明変数を離散化する離散値配列の複数の候補について、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して生成される予測テーブルを、前記格納手段に格納する説明変数の実績値及び目的変数の実績値を用いて評価して、前記離散値配列の複数の候補から離散値配列を選択することを特徴とする予測テーブル設計方法。

【請求項18】

複数の区分に分割された説明変数の各々の前記区分に対する予測対象の予測値がそれぞれ記述される予測テーブルを設計するためのプログラムであって、
前記説明変数から目的変数である前記予測対象の予測値を算出する予測モデルに、前記説明変数の各々の前記区分の代表値をそれぞれ入力して、各々の前記区分に対する前記予測対象の予測値を算出することで、前記説明変数の実績値に応じた前記予測対象の前記目的変数の実績値を予測するための前記予測テーブルを生成するためのプログラム。

【請求項19】

説明変数を離散化して形成される各区分に予測対象の予測値が記述される予測テーブルを設計するためのプログラムであって、
説明変数から目的変数である前記予測対象の予測値を算出する予測モデルを記憶する記憶手段と、
説明変数の実績値及び目的変数の実績値を格納する格納手段と、
説明変数を離散化する離散値配列の複数の候補について、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して生成される予測テーブルを、前記格納手段に格納する説明変数の実績値及び目的変数の実績値を用いて評価して、前記離散値配列の複数の候補から離散値配列を選択する評価手段としてコンピュータを機能させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、予測対象の予測値を求めるのに利用される予測テーブルを設計する予測テーブル設計装置、方法及びプログラムに関する。

【背景技術】

【0002】

予測対象の予測値を求めるために、様々な予測モデルが利用されている。例えば生産プロセスでは、所望の特性（品質や寸法等）を有する製品を製造する必要があるが、原材料の特性や温度等に依存して製品の特性も変わってしまう。そこで、製品の特性ｙを予測する式（１）の予測モデルを過去の操業データから作成し、製品の特性ｙが所望の仕様を満たすように生産プロセスの運転条件を変更することが行われている。式（１）において、変数ｙ’は製品の特性ｙの予測値であり、目的変数と呼ばれる。また、変数ｘ（＝［ｘ₁ ｘ₂・・・ｘ_M］^T）は製品の特性ｙを変化させる原因となる因子であり、説明変数と呼ばれる（Ｍは説明変数の個数である）。また、関数ｆ（ｘ）は、説明変数から目的変数を算出する関数であり、予測モデルと呼ばれる。
ｙ’＝ｆ（ｘ） ・・・（１）

【0003】

予測モデルの代表例として、式（２）のような重回帰モデルがある。これは、予測値ｙ’を説明変数ｘの線形式で表わす予測モデルであり、係数θ₀〜θ_Mは過去の操業データから算出することができる。
ｙ’＝θ₀＋θ₁ｘ₁＋θ₂ｘ₂＋・・・＋θ_Mｘ_M ・・・（２）
このような重回帰モデルは計算が容易であることから、幅広いプロセスで利用されている。

【0004】

目的変数と説明変数が線形関係にある場合は予測精度の良い予測モデルとなるが、目的変数と説明変数との間に強い非線形性がある場合には予測精度の悪い予測モデルとなってしまう。
そこで、様々な非線形予測モデルが提案されている。例えば特許文献１では、複数の領域毎に重回帰モデルを構築し、それぞれの予測値にファジー関数値を乗じた値を足し合わせ、最終的な予測値を算出する手法であり、云わば、区分線形型の予測モデルを過去操業データから自動的に構築する手法が開示されている。特許文献１以外にも、ニューラルネットワークや決定木、ランダムフォレスト等、様々な非線形予測モデルが提案されている。

【0005】

しかしながら、これらの非線形予測モデルは実産業界ではあまり普及しておらず、予測テーブルが多く利用されている。図３に予測テーブルの例を示す。図３の予測テーブルは、温度及び厚みを説明変数として、ある目的変数を予測するテーブルであり、温度及び厚みを離散化して形成される各区分（メッシュ）に予測値が記述される。図３（ａ）は温度及び厚みが共に階段形式で離散化される予測テーブルであり、例えば温度が７００〜７２０℃、厚みが５〜１０ｍｍの範囲に含まれる場合、予測値は１．３であることを表わす。また、図３（ｂ）は温度が折れ線形式、厚みが階段形式で離散化される予測テーブルであり、例えば温度ｘ₁が７００〜７２０℃、厚みが５〜１０ｍｍの範囲に含まれる場合、予測値は式（３）の線形式で表わされる。
ｙ’＝１．３＋（ｘ₁−７００）×（１．４−１．３）／（７２０−７００） ・・・（３）

【0006】

このような予測テーブルが広く普及している理由は、予測値を算出するための計算手順が単純であり、予測対象の挙動が理解し易く、予測精度の管理が容易なためである。
予測モデル（特に非線形予測モデル）では、入力変数を見ただけでは目的変数の予測値を人の頭の中で算出することは容易ではない。それに対して、予測テーブルでは、入力変数を見ただけで、予測値を容易に知ることができ、ユーザにとって非常に判り易いという利点がある。また、実操業で得られた目的変数の実績値とそれに対する予測値との誤差をメッシュ毎に集計することで、精度の悪いメッシュを容易に知ることができ、他のメッシュの精度に影響を与えることなく、精度の悪いメッシュのみ予測値を更新することでそのメッシュの精度を向上させることができる（折れ線形式の場合は隣のメッシュにも影響を与えるが、影響範囲は判り易く限定的である）。非線形予測モデルの場合、予測モデルの係数を変えると全ての製品の予測精度が変わってしまい、予測モデルの精度管理が難しいという問題がある。

【先行技術文献】

【特許文献】

【0007】

【特許文献1】特開２００９−１５１３８３号公報

【特許文献2】特開２０１２−９３３１４号公報

【特許文献3】特開平６−２５９１０７号公報

【発明の概要】

【発明が解決しようとする課題】

【0008】

上述したように実産業界では予測テーブルが広く普及しているが、予測精度の良い予測テーブルを設計するためには、メッシュの区分数や離散値、各メッシュの値（予測値）を適切に定める必要がある。離散値は、説明変数の離散化が階段形式である場合は、メッシュ間の境界値を表わし、折れ線形式である場合は、折れ線の頂点の値を表わす。例えば図３（ａ）、（ｂ）の厚みの離散値は５、１０、１５、・・・、４０となり、図３（ａ）、（ｂ）の温度の離散値は７００、７２０、・・・、８６０となる。

【0009】

従来、予測テーブルを設計するに際して、メッシュの値を定めるには、メッシュに属する実績データを取得して、予測対象の実績値の平均を用いることが一般的である。メッシュに属する実績データの個数が少ない場合には、特許文献３のように、他のメッシュの実績データを参考にメッシュの値を決定する手法も開示されている。
しかしながら、実績データからメッシュの値を算出する手法は簡便ではあるが、メッシュに属する（もしくはその周辺のメッシュに属する）少ない実績データからメッシュの値を算出することになり、過学習を引き起こしやすいという問題がある。過学習とは、過去の実績データに対しては予測精度が良いが、将来の未知のデータに対しては予測精度が悪くなる現象である。

【0010】

また、メッシュの区分数や離散値の定め方に特定の手法はなく、例えば特許文献２のように、説明変数の値域を等間隔で区切ったり、既に別の管理で用いている区分を使用したりする等している。
このように、メッシュの区分数や離散値に関しては、区分数を適当に決めた後、説明変数の値域を等間隔で区切ったり、予測値の変化が大きいところを細かく、変化が小さいところを粗く区切ったり、予測精度が悪ければ区分数を増やしたりする等、試行錯誤で定めることが一般的である。
しかしながら、試行錯誤でメッシュの区分数や離散値を定めるのでは、時間が掛るとともに、適切な数値を定めることは難しい。メッシュが粗いと、予測テーブルのサイズを小さくすることができるが、説明変数の値の細かな差異を反映させることができないため、予測精度が悪くなってしまう。一方、メッシュが細か過ぎると、予測テーブルのサイズが大きくなるばかりでなく、１つのメッシュに含まれる実績データの数が少なくなるため、メッシュの値を高精度に定めることが難しくなり、予測精度の管理も複雑化してしまう。

【0011】

本発明は上記のような点に鑑みてなされたものであり、過学習を引き起こしにくく、さらには適切な予測精度及びサイズの予測テーブルを設計できるようにすることを目的とする。

【課題を解決するための手段】

【0012】

上記の課題を解決するための本発明の要旨は、以下のとおりである。
［１］ 複数の区分に分割された説明変数の各々の前記区分に対する予測対象の予測値がそれぞれ記述される予測テーブルを設計する予測テーブル設計装置であって、
前記説明変数から目的変数である前記予測対象の予測値を算出する予測モデルに、前記説明変数の各々の前記区分の代表値をそれぞれ入力して、各々の前記区分に対する前記予測対象の予測値を算出することで、前記説明変数の実績値に応じた前記予測対象の前記目的変数の実績値を予測するための前記予測テーブルを生成することを特徴とする予測テーブル設計装置。
［２］ 説明変数を離散化して形成される各区分に予測対象の予測値が記述される予測テーブルを設計する予測テーブル設計装置であって、
説明変数から目的変数である前記予測対象の予測値を算出する予測モデルを記憶する記憶手段と、
説明変数の実績値及び目的変数の実績値を格納する格納手段と、
説明変数を離散化する離散値配列の複数の候補について、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して生成される予測テーブルを、前記格納手段に格納する説明変数の実績値及び目的変数の実績値を用いて評価して、前記離散値配列の複数の候補から離散値配列を選択する評価手段とを備えたことを特徴とする予測テーブル設計装置。
［３］ 離散値配列の初期値を定め、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して、予測テーブルの初期値を生成する予測テーブル初期値生成手段と、
前記離散値配列の初期値から一又は複数の離散値を削除した離散値配列の複数の候補を作成する離散値配列候補作成手段とを備え、
前記評価手段は、前記離散値配列候補作成手段で作成した離散値配列の複数の候補について、評価を行い、離散値配列を選択することを特徴とする［２］に記載の予測テーブル設計装置。
［４］ 前記評価手段で選択した離散値配列から生成される予測テーブルが、所定の条件を満たすか否かを判定する判定手段を備え、
前記判定手段で前記所定の条件を満たさないと判定した場合、
前記離散値配列候補作成手段は、前記評価手段で選択した離散値配列から一又は複数の離散値を削除した離散値配列の複数の候補を作成し、
前記評価手段は、前記離散値配列候補作成手段で作成した離散値配列の複数の候補について、評価を行い、離散値配列を選択し、
これを前記所定の条件を満たすまで繰り返すことを特徴とする［３］に記載の予測テーブル設計装置。
［５］ 前記所定の条件として、予測テーブルのサイズに対する目標サイズ、又は、前記評価手段での評価結果に対する目標値が設定されることを特徴とする［４］に記載の予測テーブル設計装置。
［６］ 前記評価手段で選択した離散値配列から一又は複数の離散値を削除し、かつ、前記評価手段で前記離散値配列を選択するまでに削除済みの離散値の中から、当該削除した離散値と同じ数の離散値を復活させて、改善用離散値配列候補を作成する改善用離散値配列候補作成手段と、
前記改善用離散値配列候補作成手段で作成した改善用離散値配列候補について、評価を行い、前記評価手段で選択した離散値配列よりも評価の良い改善用離散値配列候補があれば、それを前記評価手段で選択した離散値配列と置き換える改善用評価手段とを備えたことを特徴とする［４］又は［５］に記載の予測テーブル設計装置。
［７］ 前記評価手段で離散値配列を選択した後、所定の実施条件を満たすときに、前記改善用離散値配列候補作成手段及び前記改善用評価手段による処理を実施することを特徴とする［６］に記載の予測テーブル設計装置。
［８］ 前記所定の実施条件として、予測テーブルのサイズに対する条件、又は前記評価手段での評価結果に対する条件が設定されることを特徴とする［７］に記載の予測テーブル設計装置。
［９］ 前記離散値配列候補作成手段は、削除する離散値を乱数を用いて決めることを特徴とする［３］乃至［８］のいずれか一つに記載の予測テーブル設計装置。
［１０］ 前記改善用離散値配列候補作成手段は、削除する離散値を乱数を用いて決めることを特徴とする［６］乃至［８］のいずれか一つに記載の予測テーブル設計装置。
［１１］ 説明変数の離散化形式は、離散値と離散値の間の予測値は一定値とする形式、及び、離散値と離散値の間の予測値はこれら離散値に対応する予測値を結ぶ直線で表わす形式のうち少なくともいずれか一方であることを特徴とする［２］乃至［１０］のいずれか一つに記載の予測テーブル設計装置。
［１２］ 前記評価手段は、予測テーブルの予測値と、目的変数の実績値との誤差を評価することを特徴とする［２］乃至［１１］のいずれか一つに記載の予測テーブル設計装置。
［１３］ 前記評価手段は、前記誤差の平方和を評価値とする、又は、前記誤差の絶対値和を評価値とすることを特徴とする［１２］に記載の予測テーブル設計装置。
［１４］ 前記改善用評価手段は、予測テーブルの予測値と、目的変数の実績値との誤差を評価することを特徴とする［６］乃至［８］のいずれか一つに記載の予測テーブル設計装置。
［１５］ 前記改善用評価手段は、前記誤差の平方和を評価値とする、又は、前記誤差の絶対値和を評価値とすることを特徴とする［１４］に記載の予測テーブル設計装置。
［１６］ 複数の区分に分割された説明変数の各々の前記区分に対する予測対象の予測値がそれぞれ記述される予測テーブルを設計する予測テーブル設計方法であって、
コンピュータが、前記説明変数から目的変数である前記予測対象の予測値を算出する予測モデルに、前記説明変数の各々の前記区分の代表値をそれぞれ入力して、各々の前記区分に対する前記予測対象の予測値を算出することで、前記説明変数の実績値に応じた前記予測対象の前記目的変数の実績値を予測するための前記予測テーブルを生成することを特徴とする予測テーブル設計方法。
［１７］ 説明変数を離散化して形成される各区分に予測対象の予測値が記述される予測テーブルを設計する予測テーブル設計方法であって、
コンピュータが、
説明変数から目的変数である前記予測対象の予測値を算出する予測モデルを記憶する記憶手段と、説明変数の実績値及び目的変数の実績値を格納する格納手段とを用いて、
説明変数を離散化する離散値配列の複数の候補について、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して生成される予測テーブルを、前記格納手段に格納する説明変数の実績値及び目的変数の実績値を用いて評価して、前記離散値配列の複数の候補から離散値配列を選択することを特徴とする予測テーブル設計方法。
［１８］ 複数の区分に分割された説明変数の各々の前記区分に対する予測対象の予測値がそれぞれ記述される予測テーブルを設計するためのプログラムであって、
前記説明変数から目的変数である前記予測対象の予測値を算出する予測モデルに、前記説明変数の各々の前記区分の代表値をそれぞれ入力して、各々の前記区分に対する前記予測対象の予測値を算出することで、前記説明変数の実績値に応じた前記予測対象の前記目的変数の実績値を予測するための前記予測テーブルを生成するためのプログラム。
［１９］ 説明変数を離散化して形成される各区分に予測対象の予測値が記述される予測テーブルを設計するためのプログラムであって、
説明変数から目的変数である前記予測対象の予測値を算出する予測モデルを記憶する記憶手段と、
説明変数の実績値及び目的変数の実績値を格納する格納手段と、
説明変数を離散化する離散値配列の複数の候補について、前記予測モデルに各区分の説明変数の代表値を入力して予測値を算出して生成される予測テーブルを、前記格納手段に格納する説明変数の実績値及び目的変数の実績値を用いて評価して、前記離散値配列の複数の候補から離散値配列を選択する評価手段としてコンピュータを機能させるためのプログラム。

【発明の効果】

【0013】

本発明によれば、予測モデルから予測値を算出して予測テーブルを生成するようにしたので、過学習を引き起こしにくい予測テーブルを設計することができる。さらに、離散値配列の複数の候補について、評価を行い、離散値配列を選択するようにしたので、適切な予測精度及びサイズの予測テーブルを設計することができる。

【図面の簡単な説明】

【0014】

【図1】第１の実施形態に係る予測テーブル設計装置の機能構成を示す図である。

【図2】第１の実施形態に係る予測テーブル設計装置による予測テーブル設計方法を示すフローチャートである。

【図3】予測テーブルの例を示す図であり、（ａ）は温度及び厚みが共に階段形式で離散化された予測テーブルを示す図、（ｂ）は温度が折れ線形式、厚みが階段形式で離散化された予測テーブルを示す図である。

【図4】説明変数の離散値配列を説明するための図である。

【図5】図３の予測テーブルにおける説明変数の離散値配列を示す図である。

【図6】説明変数の代表値配列を説明するための図である。

【図7】図３の予測テーブルにおける説明変数の代表値配列を示す図であり、（ａ）は図３（ａ）の予測テーブルにおける説明変数の代表値配列を示す図、（ｂ）は図３（ｂ）の予測テーブルにおける説明変数の代表値配列を示す図である。

【図8】メッシュの区分数及び離散値を定め、メッシュの値を予測モデルから計算して設計する予測テーブルの例を示す図であり、（ａ）は図３（ａ）の予測テーブルに対応する予測テーブルを示す図、（ｂ）は図３（ｂ）の予測テーブルに対応する予測テーブルを示す図である。

【図9】予測テーブルの例を示す図である。

【図10】説明変数の離散化形式を説明するための図であり、（ａ）は階段形式を説明するための図、（ｂ）は折れ線形式を説明するための図である。

【図11】階段形式及び折れ線形式の説明変数が混在する場合を説明するための図である。

【図12】第２の実施形態に係る予測テーブル設計装置の機能構成を示す図である。

【図13】第２の実施形態に係る予測テーブル設計装置による予測テーブル設計方法を示すフローチャートである。

【図14】第２の実施形態における予測テーブル設計方法の概要を説明するための図である。

【発明を実施するための形態】

【0015】

以下、添付図面を参照して、本発明の好適な実施形態について説明する。
［第１の実施形態］
まず、本実施形態における予測テーブルのメッシュの区分数及び離散値、メッシュの値（予測値）の基本的な考え方について説明する。
予測テーブルを設計する際には、説明変数毎にメッシュの区分数及び離散値を定める必要がある。これは、図４に示すように、各説明変数の一次元配列（本願においては、離散値配列と呼ぶ）の要素数Ｔ_mと要素の値ｔ_m⁽ⁱ⁾を定めることに相当する。なお、ｍは説明変数の番号（ｍ：１〜Ｍ）、ｉは離散値の番号（ｉ：１〜Ｔ_m）である。

【0016】

説明変数の離散化形式には、階段形式と折れ線形式の２つがある。階段形式とは、図１０（ａ）に示すように、離散値と離散値の間の予測値は一定値とする形式である。また、折れ線形式とは、図１０（ｂ）に示すように、離散値と離散値の間の予測値は直線で表わされる形式である。なお、全ての説明変数が階段形式で離散化されていても、全ての説明変数が折れ線形式で離散化されていても良く、また、これらの形式が説明変数毎に混在していても構わない。
図４の説明変数の離散値配列の各要素の値ｔ_m⁽ⁱ⁾は、階段形式であればメッシュ間の境界値を意味し、折れ線形式であれば折れ線の頂点の値を意味する。この離散値配列が決まれば、予測テーブルは一義的に定まる。
図５には、図３（ａ）、（ｂ）の予測テーブルにおける説明変数の離散値配列を示す。

【0017】

また、本実施形態では、過学習を防ぐことを目的として、式（１）の予測モデルを構築した後、その予測モデルをテーブル形式に変換する。自由度調整済み決定係数や、ＡＩＣ（赤池情報量基準）、クロスバリデーション等の基準に基づき予測モデル選択することで、過学習を起こし難い予測モデルを構築することができる。
メッシュの区分数及び離散値を定めた後、各説明変数の代表値を予測モデルｆ（ｘ）に代入して、その結果得られる予測値を各メッシュの値にする。
図６には、図４の離散値配列から得られる説明変数の代表値配列を示す。階段形式の場合は、図１０（ａ）に示すように、各メッシュの代表値は離散値の中央値とし、式（４）により計算される。階段形式では、式（５）のように、代表値配列の数Ｒ_mは離散値配列の数Ｔ_mより１つ少なくなる。一方、折れ線形式の場合は、図１０（ｂ）に示すように、各メッシュの代表値配列は離散値配列と同じになる。
ｒ_m⁽ⁱ⁾＝（ｔ_m⁽ⁱ⁾＋ｔ_m⁽ⁱ⁺¹⁾）／２ ・・・（４）
Ｒ_m＝Ｔ_m−１ ・・・（５）
例えば図３（ａ）の予測テーブルでいえば、左上のメッシュは厚みが５〜１０ｍｍ、温度が７００〜７２０℃の範囲であるので、代表値は厚み７．５ｍｍ、温度７１０℃となる。また、図３（ｂ）の予測テーブルでいえば、折れ線の頂点（離散値）が代表値となり、左上のメッシュの代表値は厚み７．５ｍｍ、温度７００℃となる。図７（ａ）は、図３（ａ）の予測テーブルにおける説明変数の代表値配列を示し、図７（ｂ）は、図３（ｂ）の予測テーブルにおける説明変数の代表値配列を示す。

【0018】

以上のように、メッシュの区分数及び離散値が定まれば、各説明変数の代表値を用いて予測モデルにより各メッシュの値を計算することができる。図３（ａ）の予測テーブルを例にすれば、図７（ａ）の説明変数の代表値配列から、図８（ａ）のように予測テーブルが設計される。また、図３（ｂ）の予測テーブルを例にすれば、図７（ｂ）の説明変数の代表値配列から、図８（ｂ）のように予測テーブルが設計される。

【0019】

以下、本実施形態における予測テーブル設計の詳細について説明する。本実施形態では、過学習を引き起こしにくく、さらには適切な予測精度及びサイズの予測テーブルを設計することを目的として、以下に詳述するようにして予測テーブルを設計する。
図１に、第１の実施形態に係る予測テーブル設計装置１２００の機能構成を示す。
１０１は記憶手段として機能する予測モデル記憶部であり、式（１）のように、既に構築されている予測モデルｆ（ｘ）を記憶する。
１０２は格納手段として機能する実績データ格納部であり、過去の操業で得られた実績データ｛ｘ⁽ⁿ⁾，ｙ⁽ⁿ⁾｝を格納する。ｘ⁽ⁿ⁾＝［ｘ₁⁽ⁿ⁾，ｘ₂⁽ⁿ⁾，・・・ｘ_M⁽ⁿ⁾］^Tはｎ番目の実績データの説明変数、ｙ⁽ⁿ⁾はｎ番目の実績データの目的変数であり、実績データはＮ個の説明変数と目的変数の組み合せから構成される。

【0020】

１０３は予測テーブル初期値生成手段として機能する予測テーブル初期値生成部であり、予め設定された区分数で各説明変数を分割することにより離散値配列の初期値を定め、各説明変数の代表値を用いて予測モデルｆ（ｘ）により各メッシュの値を計算して予測テーブルの初期値を生成する。
１０４は離散値配列候補作成手段として機能する離散値配列候補作成部であり、対象の離散値配列から離散値を削除した離散値配列の複数の候補（以下、複数の離散値配列候補と記す）を作成する。対象の離散値配列とは、第一段階（後述の繰り返し計算の１回目）では、予測テーブル初期値生成部１０３で定めた離散値配列の初期値が、それ以降の段階では、評価部１０５で選択した離散値配列候補が該当する。

【0021】

１０５は評価手段として機能する評価部であり、実績データ格納部１０２に格納する実績データを用いて、離散値配列候補作成部１０４で作成した各離散値配列候補から生成される予測テーブルを評価し、その評価に基づいて離散値配列候補を選択する。
１０６は判定手段として機能する判定部であり、評価部１０５で選択した離散値配列候補から生成される予測テーブルが所定の条件である終了条件を満たすか否かを判定する。終了条件は、後述するように予測テーブルの予測精度及びサイズの観点から予め設定されている。終了条件を満たしていなければ、評価部１０５で選択した離散値配列候補を対象の離散値配列として、離散値配列候補作成部１０４及び評価部１０５による処理を行い、これを終了条件を満たすまで繰り返す。

【0022】

１０７は出力部であり、判定部１０６で終了条件を満たしていれば、予測テーブルを出力する。例えば予測テーブルそのもの、又は予測テーブルに関する情報を表示装置に表示したり、ネットワークを介して外部機器に送信したりすることをいう。

【0023】

図２は、第１の実施形態に係る予測テーブル設計装置１００による予測テーブル設計方法を示すフローチャートである。
ステップＳ１で、予測テーブル初期値生成部１０３は、予め設定された十分に大きな区分数で各説明変数を分割することにより離散値配列の初期値を定め、各説明変数の代表値を用いて予測モデルｆ（ｘ）により各メッシュの値を計算して予測テーブルの初期値を生成する。
ここで、十分に大きな区分数で各説明変数を分割するとは、予測モデルｆ（ｘ）をテーブル化したときに予測精度が落ちないと考えられる細かさで各説明変数を分割するという意味である。例えば技術的・経験的見地から、説明変数は０．５％の精度で計算すれば十分な予測精度が得られるということがわかっている場合、各説明変数の最小値と最大値を等間隔に２００個に分割して離散値を設定すれば良い。また、説明変数の変化に対して、予測値の変化が小さい領域がある場合、その変化の傾きに応じて離散値と離散値の間隔を調整しても良い。すなわち、傾きの大きい部分では区分を細かくし、傾きの小さい部分では区分を粗くする。
離散値配列の初期値を定めた後、前記（４）式又は（５）式を用いてメッシュ毎に各説明変数の代表値を求め、それら代表値用いて予測モデルｆ（ｘ）により各メッシュの値を計算して予測テーブルの初期値を生成する。図９に、予測テーブルの例を示す。

【0024】

ステップＳ２で、離散値配列候補作成部１０４は、対象の離散値配列から１つの離散値（ｍ番目の説明変数のｉ番目の離散値）を削除した複数の離散値配列候補を作成する。対象の離散値配列は、初回のループでは、ステップＳ１で定めた離散値配列の初期値が該当する。
この場合に、説明変数の個数と区分数が少ない場合、説明変数１の１番目の離散値から順番に全数探索（離散値配列候補数＝Ｒ₁＋Ｒ₂＋・・・＋Ｒ_M）を行っても良いが、全数探索に時間が掛る場合には、削除する離散値を乱数で決めて複数の離散値配列候補を作成するようにしても良い。
また、本実施形態では、離散値配列から１つの離散値を削除するようにしたが、複数の離散値を一度に削除するようにしても良い。

【0025】

ステップＳ３で、評価部１０５は、実績データ格納部１０２に格納する実績データを用いて、ステップＳ２で作成した各離散値配列候補から生成される予測テーブルを評価し、その評価に基づいて離散値配列候補を選択する。
予測テーブルの評価は、次のように行う。
実績データ格納部１０２に格納する実績データ｛ｘ⁽ⁿ⁾，ｙ⁽ⁿ⁾｝のそれぞれの説明変数ｘ_m⁽ⁿ⁾に対して、式（６）を満たす離散値の番号ｉ_mを求める。
ｔ_m^(im)≦ｘ_m⁽ⁿ⁾＜ｔ_m^(im+1) ・・・（６）
簡単のため、説明変数が１つのみの場合で説明すると、階段形式であれば、図１０（ａ）に示すように、離散値番号ｉ_m（代表値がｒ_m^(im)）の予測値を図９の予測テーブルから求め、目的変数の予測値ｙ’⁽ⁿ⁾とする。また、折れ線形式であれば、代表値ｒ_m^(im)とｒ_m^(im+1)の予測値を図９の予測テーブルからそれぞれ求め（それら予測値をｙ₁’⁽ⁿ⁾，ｙ₂’⁽ⁿ⁾とする）、図１０（ｂ）に示すように、それらを案分して予測値ｙ’⁽ⁿ⁾を求める。この計算式は、式（７）となる。
ｙ’⁽ⁿ⁾＝ｙ₁’⁽ⁿ⁾＋（ｘ_m⁽ⁿ⁾−ｒ_m^(im)）×（ｙ₂’⁽ⁿ⁾−ｙ₁’⁽ⁿ⁾）／（ｒ_m^(im+1)−ｒ_m^(im)） ・・・（７）

【0026】

説明変数が複数の場合も同様であり、例えば図１１は、説明変数１が階段形式、説明変数２が折れ線形式の場合に予測値ｙ’⁽ⁿ⁾を算出する例を説明するためのである。式（６）により、説明変数１と説明変数２が属する離散値番号ｉ₁、ｉ₂を求め、ｙ₁’⁽ⁿ⁾＝ｆ（ｒ₁⁽ⁱ¹⁾，ｒ₂⁽ⁱ²⁾），ｙ₂’⁽ⁿ⁾＝ｆ（ｒ₁⁽ⁱ¹⁾，ｒ₂⁽ⁱ²⁺¹⁾）を案分する点を式（７）により求めて予測値とすれば良い。２つの説明変数が共に折れ線形式である場合には、４つの点ｆ（ｒ₁⁽ⁱ¹⁾，ｒ₂⁽ⁱ²⁾），ｆ（ｒ₁⁽ⁱ¹⁺¹⁾，ｒ₂⁽ⁱ²⁾），ｆ（ｒ₁⁽ⁱ¹⁾，ｒ₂⁽ⁱ²⁺¹⁾），ｆ（ｒ₁⁽ⁱ¹⁺¹⁾，ｒ₂⁽ⁱ²⁺¹⁾）を案分する点が予測値となる。

【0027】

このようにして、実績データ｛ｘ⁽ⁿ⁾，ｙ⁽ⁿ⁾｝に対応する予測値ｙ₁’⁽ⁿ⁾を求め、式（８）のように、予測テーブルの予測精度を評価する評価値Ｊを計算する。なお、Σ^N_n=1の表記はΣの上にＮが、Σの下にｎ＝１が付されているものとする。
Ｊ＝Σ^N_n=1（ｙ’⁽ⁿ⁾−ｙ⁽ⁿ⁾）² ・・・（８）
なお、式（８）は予測値と実績値の誤差平方和を評価値とする式であるが、式（９）のように予測値と実績値の誤差絶対値和を評価値とする式としても良い。
Ｊ＝Σ^N_n=1｜ｙ’⁽ⁿ⁾−ｙ⁽ⁿ⁾｜ ・・・（９）

【0028】

以上のように、複数の離散値配列候補から生成される予測テーブルを評価し、評価値Ｊが最小となる、すなわち予測テーブルの予測精度が最も良い離散値配列候補を選択する。
なお、ステップＳ２で複数（例えば２つ）の離散値を一度に削除する場合、全ての説明変数の全ての離散値のうち、２つの離散値を削除した離散値配列候補を、全ての離散値の組み合せに対して作成し、最も予測精度の良い離散値配列候補を求めるようにしても良い。また、離散値の組み合わせの数が多い時には、乱数で選んだ２つの離散値を削除するようにしても良い。

【0029】

ステップＳ４で、判定部１０６は、ステップＳ３で選択した離散値配列候補から生成される予測テーブルが終了条件を満たすか否かを判定する。終了条件を満たしていれば、ステップＳ５に進み、得られた離散値配列と、その離散値配列から生成される予測テーブルを出力して終了する。一方、終了条件を満たしていなければ、ステップＳ３で選択した離散値配列候補を対象の離散値配列として、ステップＳ２及びＳ３を行い、これを終了条件を満たすまで繰り返す。

【0030】

ここで、予測テーブルは、主として次の２つの考え方に基づいて設計される。予測テーブルのサイズの上限値（以下、目標サイズと記載する）がメモリ制約や管理上の制約等で予め決まっている場合は、目標サイズに収まる予測テーブルの中で、最も予測精度の良い予測テーブルを設計することが望まれる。また、テーブルサイズに制限がなく、目標の予測精度が与えられている場合は、目標の予測精度を満足する予測テーブルの中で、最もサイズの小さい予測テーブルを設計することが望まれる。テーブルは小さい方が管理し易いからである。
このように、終了条件には、目標サイズを設定する場合と、目標精度を設定する場合とがある。

【0031】

目標サイズを設定する場合、この目標サイズ以下になるまで、ステップＳ２及びＳ３を繰り返す。予測テーブルの初期値を出発点として離散値を削除していくことで、予測テーブルのメッシュが粗くなり、予測テーブルのサイズは小さくなる。そこで、ステップＳ３で選択した離散値配列候補から生成される予測テーブルのサイズが目標サイズ以下になったならば、そこからさらにサイズを小さくする必要はないので、ステップＳ２及びＳ３を終了する。目標サイズは、説明変数毎の区分数の目標値でも、それらを乗じた値の目標値でも良く、それらの組合せでも構わない。
例えば目標サイズとして、条件１〜条件３が設定されたとする。
（条件１）説明変数１区分数≦２０、説明変数２区分数≦１０、説明変数３区分数≦１０
（条件２）説明変数１区分数≦２０、説明変数２区分数×説明変数３区分数≦５０
（条件３）説明変数１区分数×説明変数２区分数×説明変数３区分数≦２００
条件１では、全ての説明変数の区分数が目標値となったら終了する。この場合に、繰り返しの途中で目標値に到達した説明変数の離散値は削除しないようにする。また、条件２では、説明変数１の区分数が２０、説明変数２と説明変数３の区分数を乗じた値が５０以下になったら終了する。また、条件３では、３つの説明変数の区分数を乗じた値が２００以下になったら終了する。
このように目標サイズを設定することにより、目標サイズ以内で最も予測精度の良い予測テーブルを設計することができる。

【0032】

一方、目標精度を設定する場合、評価値Ｊに対する目標値を設定し、予測テーブルの評価値Ｊが目標値を上回るまで、ステップＳ２及びＳ３を繰り返す。予測テーブルの初期値を出発点として離散値を削除していくことで、予測テーブルのメッシュが粗くなり、予測値と実績値の誤差平方和或いは誤差絶対値和等、予測値と実績値との誤差で表わされる評価値Ｊは大きく、すなわち予測精度は悪くなる傾向となる。そこで、ステップＳ３で選択した離散値配列候補から生成される予測テーブルの評価値Ｊが目標値を上回ると、そこからさらに予測テーブルのメッシュを粗くすると目標精度から大きく劣ることになるので、ステップＳ２及びＳ３を終了する。この場合、最後の離散値配列から生成される予測テーブルは目標精度より劣るため、１回前の離散値配列及び予測テーブルを出力する。
このように目標精度を設定することにより、目標精度を満足する最もサイズの小さい予測テーブルを設計することができる。

【実施例1】

【0033】

厚板の製造プロセスでは、γ線厚み計を用いて鋼板の板厚を測定する。γ線厚み計は、鋼板を透過するγ線の減衰率を板厚に換算する方式であるが、鋼板の温度や材質等に依存して同じ厚みでも減衰率が異なることから、γ線の標準減衰率（温度や材質等に依存しない固定値）から計算した標準的な板厚であるγ線標準板厚を測定板厚にすることはできず、式（１０）のように、鋼板の温度や材質等から計算される補正係数Ｃｔを乗じて、最終的な測定値であるγ線測定板厚を算出する。
γ線測定板厚＝（１＋Ｃｔ）×γ線標準板厚 ・・・（１０）

【0034】

この補正係数Ｃｔを、検査工程にあるレーザ式板厚計の測定値（板厚の真値とみなす）とγ線測定板厚が一致するように、特許文献１の手法を用いて非線形予測モデルを設計して計算したところ、予測モデルの予測精度は表１のようになった。なお、説明変数は鋼板温度、変態温度、狙い厚の３つとし、２９０５８個の実績データを用いて非線形モデルを設計、得られた予測モデルの予測精度はそれとは異なる１９４９７個の実績データを用いて評価した。表１は既設テーブルからの改善度を示している。この既設テーブルは、鋼板温度２０区分（ほぼ等間隔）、鋼種１００区分のテーブルであり、テーブル管理者は定期的に予測精度をメッシュ毎に確認し、精度の悪いメッシュの値の修正や、鋼種の区分を見直しを手作業で行っていた。表１のＣｔ精度改善度の誤差、σ、ＲＭＳＥ（Root mean square error; 平均２乗誤差）はそれぞれ、レーザ板厚計の測定値から計算したＣｔの真値と既設テーブルに設定されているＣｔとの誤差の平均値，標準偏差、ＲＭＳＥと、前記Ｃｔの真値と本発明を適用して計算したＣｔとの誤差の平均値、標準偏差、ＲＭＳＥとの差異（既設の値−本発明を適用した値）を表わしている。板厚精度改善度の誤差、σ、ＲＭＳＥも同様である。また、以降で説明する表２〜表６についても表１と同様である。本発明を適用した予測テーブル設計手法により、補正係数ＣｔのＲＭＳＥは既設テーブルと比べて０．００９％改善し、板厚測定精度もＲＭＳＥが３５μｍ改善した。

【0035】

【表1】

【0036】

この非線形予測モデルを手動でテーブル化した場合の予測精度を表２に示す。予測テーブルは、鋼板温度は折れ線形式、変態温度及び狙い厚は階段形式とした。手動での予測テーブル設計においては、各説明変数の区分数を予測テーブルの予測誤差を見ながら調整し、離散値配列は各説明変数の最小と最大の範囲を等間隔に区切った。試行錯誤の結果、鋼板温度は２０区分、変態温度は１１区分、狙い厚は６区分に分割する予測テーブルが得られた。表２のように、予測テーブルを手動で設計すると、補正係数ＣｔのＲＭＳＥ改善度は０．００９％から０．００４％に、板厚測定精度のＲＭＳＥ改善度も３５μｍから２９μｍに減少してしまった。

【0037】

【表2】

【0038】

そこで、本発明を適用した予測テーブル設計手法を利用して、目標サイズを「鋼板温度区分数≦２０、変態温度区分数×狙い厚区分数≦６８」の条件で予測テーブルを設計した。これは、手動での設計とほぼ同じサイズの予測テーブルを求める条件である。その結果、鋼板温度を２０区分、変態温度を７区分、狙い厚を９区分する予測テーブルが得られ、予測テーブルの予測精度（既設テーブルからの改善度）は表３のようになった。表３のように、補正係数ＣｔのＲＭＳＥ改善度は０．０１３％、板厚測定精度のＲＭＳＥ改善度は３８μｍとなり、手動でテーブル化した場合と比較して大きな改善が得られたのみならず、もともとの非線形予測モデルより高精度の予測テーブルを設計することができた。元々の非線形予測モデルは、実績データの目的変数の値を完全に予測できるモデルではないため、このようにテーブル化することによってたまたま予測精度が向上する場合もある。

【0039】

【表3】

【0040】

また、目標精度を設定して、サイズの小さい予測テーブルを求めた。目標精度は、補正係数ＣｔのＲＭＳＥ改善度が０．００７％となるように設定し、予測テーブルを設計した（非線形予測モデルの補正係数ＣｔのＲＭＳＥよりも０．００２％悪いＲＭＳＥを目標精度とした）。その結果、鋼板温度を４区分、変態温度を４区分、狙い厚を４区分する予測テーブルが得られ、その予測テーブルの予測精度は表４のようになった。目標精度は学習データに対する条件のため、評価データの補正係数ＣｔのＲＭＳＥ改善度は目標精度と一致しないが、表１と比べると、板厚測定精度のＲＭＳＥ改善度は非線形予測モデルよりも若干悪い程度(改善度３５μｍ→３３μｍ)であり、手動で設計したときよりも圧倒的に小さい予測テーブルで（テーブルサイズ２０×１１×６＝１３２０→４×４×４＝６４)、高精度（板厚ＲＭＳＥ改善度２９μｍ→３３μｍ）の予測テーブルを設計することができた。

【0041】

【表4】

【0042】

［第２の実施形態］
次に、第２の実施形態を説明する。以下では、第１の実施形態との相違点を中心に説明し、第１の実施形態との共通点についてはその詳細な説明を省略する。
第１の実施形態では、離散値配列の初期値を出発点として、離散値を削除して複数の離散値配列候補を作成し、その中から離散値配列を選択することを繰り返す。このように離散値を段階的に少なくしていき、目標サイズ又は目標精度に達した離散値配列を求めるのであるが、そのときに、離散値配列に含まれる離散値を、削除済みの離散値と置換することにより、離散値配列の要素数は変えずに、予測精度を向上させられる可能性があることを見出した。
そこで、第２の実施形態では、離散値配列に含まれる離散値を、削除済みの離散値と置換するアルゴリズムを取り入れた形態を説明する。

【0043】

図１２に、第２の実施形態に係る予測テーブル設計装置１２００の機能構成を示す。なお、第１の実施形態に係る予測テーブル設計装置１００と同様の構成要素には同一の符号を付し、その説明を省略する。
予測テーブル設計装置１２００において、予測モデル記憶部１０１、実績データ格納部１０２、予測テーブル初期値生成部１０３、離散値配列候補作成部１０４、評価部１０５、判定部１０６、出力部１０７は第１の実施形態と同様であり、精度改善部１０８を追加している。
精度改善部１０８は、改善用離散値配列候補作成手段として機能する改善用離散値配列候補作成部１０９と、改善用評価手段として機能する改善用評価部１１０とを備える。

【0044】

改善用離散値配列候補作成部１０９は、評価部１０５で選択した離散値配列候補（以下、現在の離散値配列候補と呼ぶ）に含まれる離散値を、削除済みの離散値と置換する。詳しくは、改善用離散値配列候補作成部１０９は、現在の離散値配列候補から離散値を削除し、かつ、現在の離散値配列候補を選択するまでに削除済みの離散値の中から、当該削除した離散値と同数の離散値を復活させて、改善用離散値配列候補を作成する。同数の離散値を置換するので、改善用離散値配列候補は、現在の離散値配列候補と要素数が同じで、予測テーブルのサイズが変わらない離散値配列となる。

【0045】

改善用評価部１１０は、改善用離散値配列候補作成部１０９で作成した各改善用離散値配列候補について、評価部１０５で行った評価と同様の評価を行う。詳しくは、改善用評価部１１０は、実績データ格納部１０２に格納する実績データを用いて、改善用離散値配列候補作成部１０９で作成した各改善用離散値配列候補から生成される予測テーブルを評価する。その結果、現在の離散値配列候補よりも評価の良い改善用離散値配列候補がなければ、現在の離散値配列候補のままとする。それに対して、現在の離散値配列候補よりも評価の良い改善用離散値配列候補があれば、それを現在の離散値配列候補と置き換える。すなわち、現在の離散値配列候補ではなく、現在の離散値配列候補よりも評価の良い改善用離散値配列候補を採用する。なお、現在の離散値配列候補よりも評価の良い改善用離散値配列候補が複数あれば、その中で最も評価の良いものを選んで、それを現在の離散値配列候補と置き換えれば良い。
本実施形態では、改善用評価部１１０が評価部１０５で行った評価と同様の評価を行うとしたが、改善用評価部１１０と評価部１０５とで評価の方式を異ならせても良い。例えば、評価部１０５で使用した式（式（８）や式（９）等）を用いて計算した評価値Ｊに所定の定数値を加えた値を改善用評価部１１０の評価値としても良い。これにより、一度削除した離散値を復活させるためには、現在の離散値配列候補の評価値よりも、前記所定の定数値の分だけ評価値が良くならなくてはならず、一度削除した離散値は復活し難くなり、ある離散値が削除と復活を繰り返すことを防止することができる。

【0046】

図１３は、第２の実施形態に係る予測テーブル設計装置１００による予測テーブル設計方法を示すフローチャートである。
ステップＳ１〜Ｓ５は第１の実施形態と同様であり、本実施形態では、ステップＳ３とステップＳ４との間にステップＳ６〜Ｓ８を追加している。

【0047】

ステップＳ６で、精度改善部１０８は、所定の実施条件を満たすか否かを判定することにより、ステップＳ７及びＳ８の精度改善処理を実施するか否かを判定する。精度改善処理を実施する場合、ステップＳ７に進み、実施しない場合、ステップＳ４に進む。
所定の実施条件としては、例えばステップＳ３で選択した離散値配列候補（現在の離散値配列候補）から生成される予測テーブルのサイズが所定のサイズ（改善実施サイズと呼ぶ）や所定の予測精度（改善実施精度と呼ぶ）を満たすときに精度改善処理を実施するようにしても良い。改善実施サイズは、既述した目標サイズと同じでも良い。この場合、現在の離散値配列候補から生成される予測テーブルのサイズが目標サイズに達したときに、精度改善処理が実施される。改善実施サイズが目標サイズよりも大きい場合には、予測テーブルのサイズが目標サイズになるまで、ステップＳ４を経て、ステップＳ２〜Ｓ８が繰り返されることになる。同様に、改善実施精度は、既述した目標精度と同じでも良いし、目標精度よりも高い予測精度でも良い。
なお、ステップＳ６の判定を行わずに、毎回のループにおいて常に精度改善処理を実施するようにしても良い。

【0048】

ステップＳ７で、改善用離散値配列作成部１０９は、現在の離散値配列候補から１つの離散値（ｍ番目の説明変数のｉ番目の離散値）を削除し、かつ、削除済みの離散値の中から、１つの離散値を復活させて、改善用離散値配列候補を作成する。
ここで、削除済みの離散値とは、評価部１０５で現在の離散値配列候補を選択するまでに削除された離散値のことである。換言すれば、ステップＳ１で定めた離散値配列の初期値のうち、今回のループのステップＳ３で選択した離散値配列候補に含まれていない離散値である。
また、予測テーブルのサイズが変わることを防ぐため、削除する離散値の説明変数と復活させる離散値の説明変数とは同じ（説明変数の番号ｍが同じ）とした方が良い。例えば厚みの離散値を削除するのであれば、厚みの離散値を復活させ、温度の離散値を削除するのであれば、温度の離散値を復活させる。
また、改善用離散値配列候補を作成する際、説明変数の個数と区分数が少ない場合、説明変数１の現在の離散値配列候補の１番目の離散値から順番に削除する離散値を選択し、説明変数１の離散値配列の初期値の１番目の離散値から順番に、評価部１０５で現在の離散値配列候補を選択するまでに削除された離散値を探索するような全数探索を行っても良い。また、削除する離散値を乱数で決めて改善用離散値配列候補を作成するようにしても良い。

【0049】

ステップＳ８で、改善用評価部１１０は、実績データ格納部１０２に格納する実績データを用いて、ステップＳ７で作成した各改善用離散値配列候補から生成される予測テーブルを評価する。ここでの評価方法はステップＳ３と同じであり、各改善用離散値配列候補から生成される予測テーブルを用いて実績データに対応する予測値を求め、式（８）や式（９）で得られる評価値Ｊを計算する。そして、現在の離散値配列候補よりも評価の良い改善用離散値配列候補があれば、それを現在の離散値配列候補と置き換える。

【0050】

なお、離散値配列の要素数をＴ_m（ｍは説明変数の番号）、離散値配列の初期値の要素数をＴ_mSとすると、全数探索で求められる改善用離散値配列候補の数は、Ｔ₁×（Ｔ_1S−Ｔ₁）＋Ｔ₂×（Ｔ_2S−Ｔ₂）＋・・・＋Ｔ_M×（Ｔ_MS−Ｔ_M）となる。このため、Ｔ_m（もしくはＴ_mS−Ｔ_m）が小さいときには、全数探索を行っても改善用離散値配列候補の数はそれほど多くならない。したがって、Ｔ_m（もしくはＴ_mS−Ｔ_m）が小さいときのみステップＳ７及びＳ８の精度改善処理を全数探索によって実施することで、常に全数探索を行う場合と比較すると、予測テーブルの予測精度が多少悪くなる可能性があるものの、演算時間及び必要メモリを小さくすることができる。経験上、離散値配列の要素数Ｔ_mが小さいときのみ精度改善処理を実施することが効果的であるが、Ｔ_mS−Ｔ_mが小さいとき（離散値配列の要素数と初期の要素数が大きく変わらないとき）にも精度改善処理を実施しても良い。
また、１つの離散値を削除して、１つの離散値を復活させたが、削除する個数と復活させる個数が同じであれば良く、複数の離散値を一度に削除して、同数の離散値を復活させた改善用離散値配列候補を作成するようにしても良い。

【0051】

図１４を参照して、第２の実施形態における予測テーブル設計方法の概要を説明する。ここでは、説明の簡略化のため、単一の説明変数で、１つの離散値を削除していく例とする。
ステップＳ１で定めた離散値配列の初期値が、図１４の上段に示すように、［ｔ⁽¹⁾ ｔ⁽²⁾ ｔ⁽³⁾ ｔ⁽⁴⁾］であるとする。
初回のループのステップＳ２で、図１４の中段に示すように、４つの離散値配列候補が作成され、ステップＳ３で、星印で示すように、離散値ｔ⁽³⁾を削除した離散値配列候補［ｔ⁽¹⁾ ｔ⁽²⁾ ｔ⁽⁴⁾］が選択されたとする。
次に、ステップＳ６及びＳ４を経て、ステップＳ２に戻ったとする。この２回目のループのステップＳ２で、図１４の下段に示すように、３つの離散値配列候補が作成され、ステップＳ３で、星印で示すように、離散値ｔ⁽²⁾を削除した離散値配列候補［ｔ⁽¹⁾ ｔ⁽⁴⁾］が選択されたとする。
次に、ステップＳ６を経て精度改善処理を実施するとしてステップＳ７に進んだとする。ステップＳ７で、現在の離散値配列候補［ｔ⁽¹⁾ ｔ⁽⁴⁾］から１つの離散値を削除し、かつ、削除済みの離散値の中から、１つの離散値を復活させて、改善用離散値配列候補を作成する。この場合に、削除済みの離散値とは、現在の離散値配列候補を選択するまでに削除した離散値であり、離散値配列初期値［ｔ⁽¹⁾ ｔ⁽²⁾ ｔ⁽³⁾ ｔ⁽⁴⁾］のうち、ステップＳ３で選択した離散値配列候補［ｔ⁽¹⁾ ｔ⁽⁴⁾］に含まれていない離散値ｔ⁽²⁾、ｔ⁽³⁾となる。そして、図１４に示す４つの改善用離散値配列候補について評価を行い、現在の離散値配列候補［ｔ⁽¹⁾ ｔ⁽⁴⁾］よりも評価の良い改善用離散値配列候補があれば、それを現在の離散値配列候補［ｔ⁽¹⁾ ｔ⁽⁴⁾］と置き換える。
図１４に示す４つの改善用離散値配列候補について評価を行っても良いが、図１４から分るように、改善用離散値配列候補の［ｔ⁽¹⁾ ｔ⁽²⁾］［ｔ⁽²⁾ ｔ⁽⁴⁾］は、今回のループのステップＳ３で既に評価した離散値配列候補と同じである。従って、前記２つの改善用離散値配列候補は、現在の離散値配列候補［ｔ⁽¹⁾ ｔ⁽⁴⁾］よりも評価は良くないことが明らかであるため、これらは改善用離散値配列候補から除外しても良い。すなわち、改善用離散値配列候補を作成する際、ステップＳ２で求めた離散値配列候補と重複するものは、改善用離散値配列候補から除外しても良い。

【実施例2】

【0052】

表４で説明したケースのように、目標精度を設定して、サイズの小さい予測テーブルを求めるケースに対して、第２の実施形態の手法で全数探索により予測テーブルを作成した。目標精度は、表４で説明したケースと同じく、補正係数ＣｔのＲＭＳＥ改善度が０．００７％となるように設定した。
まず、常に精度改善処理を実施する予測テーブルを作成した。その結果、鋼板温度を４区分、変態温度を３区分、狙い厚を３区分する予測テーブルが得られ、その予測テーブルの予測精度は表５のようになった。表４と比べると、予測テーブルのサイズが６４から３６（４×３×３＝３６）に小さくなったにも関わらず、板厚測定精度のＲＭＳＥ改善度は３３μｍであり、表４と同じ値となった。

【0053】

【表5】

【0054】

ただし、ステップＳ６の判定を行わずに常に精度改善処理を実施する場合、計算時間が約１２時間も掛ってしまった。そこで、ステップＳ６において所定の実施条件として改善実施精度を目標精度と同じ値に設定し、予測テーブルの予測精度が改善実施精度よりも悪化すれば、精度改善処理を実施するようにした。その結果、表５の予測テーブルと変わらないサイズの予測テーブル（４×３×３＝３６）が得られ、その予測精度は表６のようになった。表６の予測テーブルは３０分で計算でき、板厚測定精度のＲＭＳＥ改善度もほぼ変わらない良好な結果が得られた。

【0055】

【表6】

【0056】

以上、本発明を種々の実施形態と共に説明したが、本発明はこれらの実施形態にのみ限定されるものではなく、発明の範囲内で変更等が可能である。例えば図１、図１２では、予測モデル記憶部１０１及び実績データ格納部１０２を予測テーブル設計装置１００、１２００の外部機器として示したが、内部に含むようにしても良い。また、予測テーブル設計装置１００、１２００が、予測モデルや実績データを、ネットワークを介して外部の予測モデル記憶部及び実績データ格納部から取得するような構成としても良い。
本発明を適用した予測テーブル設計装置は、例えばＣＰＵ、ＲＯＭ、ＲＡＭ等を備えたコンピュータ装置により実現される。
また、本発明は、本発明の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータがプログラムを読み出して実行することによっても実現可能である。

【符号の説明】

【0057】

１００、１２００：予測テーブル設計装置、１０１：予測モデル記憶部、１０２：実績データ格納部、１０３：予測テーブル初期値生成部、１０４：離散値配列候補作成部、１０５：評価部、１０６：判定部、１０７：出力部、１０８：精度改善部、１０９：改善用離散値配列候補作成部、１１０：改善用評価部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】