JP2004355189A - Result predicting device - Google Patents

Result predicting device Download PDF

Info

Publication number
JP2004355189A
JP2004355189A JP2003150345A JP2003150345A JP2004355189A JP 2004355189 A JP2004355189 A JP 2004355189A JP 2003150345 A JP2003150345 A JP 2003150345A JP 2003150345 A JP2003150345 A JP 2003150345A JP 2004355189 A JP2004355189 A JP 2004355189A
Authority
JP
Japan
Prior art keywords
condition
result
distance
calculating
similarity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2003150345A
Other languages
Japanese (ja)
Other versions
JP4474847B2 (en
Inventor
Hiroyasu Shigemori
弘靖 茂森
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JFE Steel Corp
Original Assignee
JFE Steel Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by JFE Steel Corp filed Critical JFE Steel Corp
Priority to JP2003150345A priority Critical patent/JP4474847B2/en
Publication of JP2004355189A publication Critical patent/JP2004355189A/en
Application granted granted Critical
Publication of JP4474847B2 publication Critical patent/JP4474847B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Investigating And Analyzing Materials By Characteristic Methods (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To highly quickly predict a result corresponding to a requested condition without inputting any special rule even when an object is complicated and non-linear. <P>SOLUTION: This result predicting device is provided with a result database 10 for storing the value of the past condition and a result acquired from the condition, a means 12 for calculating an influence coefficient for the result of each condition in the neighborhood of a requested condition whose results should be predicted in a condition space specified by the condition stored in the result database 10, a means 14 for converting the axis of the condition space based on the acquired influence coefficient, and for calculating a distance between the value of the past condition stored in the result database and the requested condition in the converted condition space, a means 16 for calculating similarity between the value of each condition and the requested condition based on the acquired distance, a means 18 for preparing a prediction formula in the neighborhood of the requested condition based on the acquired similarity and a means 20 for calculating a result corresponding to the requested condition based on the acquired prediction formula. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

【0001】
【発明が属する技術分野】
本発明は、結果予測装置に係り、特に過去に適用した条件及びその条件を適用して得られた結果を実績データとして格納したデータベースを使用し、任意の要求条件に対する結果を予測する際に適用して好適な結果予測装置に関する。
【0002】
【従来の技術】
任意の要求条件に対して結果を予測する方法としては、例えば、実績データベースに格納されている過去の各観測データの条件と、所望の要求条件との距離を計算し、その距離から観測データ(実績データ)の重みを計算し、その重みから要求条件の近傍をフィッティングする関数を作成し、その関数を用いて要求条件に基づく結果を予測する方法が知られている(例えば、非特許文献1、2、特許文献1を参照)。
【0003】
任意の要求条件と、過去の実績データの各条件との距離の計算には、各条件を変数とする軸に規定される条件空間におけるノルムが用いられる。非特許文献1、2には、ノルムの一種である、ユークリッド距離、正規化ユークリッド距離が、又、特許文献1には、マハラノビス距離を用いることが記載されている。
【0004】
特に、特許文献1には、正規化ユークリッド距離を使用すると、個々の条件空間の軸(条件軸)が正規分布になるように正規化されているため、ユークリッド距離を用いる場合に問題があった項目による分布の違いが影響しないようになることが記述されている。又、マハラノビス距離を用いる場合は、正規化ユークリッド距離に加え、項目間の相関を取り除くために新たに軸の方向を取り直し、互いに無相関な軸とすることにより、特別な操業を他の操業と区別することができ、より良くフィッティングできるようになることが記述されている。
【0005】
又、特許文献2には、条件空間の位置により結果に対して影響を与える要因が変化する対象に対する予測方法が開示されている。これは、条件空間の位置に応じて条件空間の軸をルールに従って決定するようにしたことにより、複雑・非線形な対象においても要求条件に対する結果を精度良く予測できるようにしたものである。
【0006】
【特許文献1】
特開2001−290508号公報
【特許文献2】
特開2002−236119号公報
【非特許文献1】
William S.Cleveland and Susan J.Devlin;
Locally Weighted Regression:An approach to
Regression Analysis by Local Fitting,Journal of the American Statistical Association,Vol.8 3, No.403,September 1988.
【非特許文献2】
Zheng,Q.,and H.Kimura;Locally Weighted
Regression Based on k Bipartite Neighbors,
Proc.42nd Japan Joint Automatic Control
Conference,pp143,Tokyo (1999).
【0007】
【発明が解決しようとする課題】
しかしながら、非特許文献1、2及び特許文献1に記載されているような方法、即ち実績データの条件空間を線形変換する方法(ユークリッド距離は実績データの条件軸をそのまま用いる。正規化ユークリッド距離は実績データの条件軸を各条件の分布(標準偏差)によりスケーリングする。マハラノビス距離は実績データの条件軸を互いに無相関になるように線形変換する。)には、条件空間の位置に応じて、各条件により得られる結果に対する影響が変化する複雑・非線形な対象に対する予測精度には限界がある。
【0008】
又、特許文献2に記載されている方法、即ち条件空間の位置に応じて条件空間の軸をルールに従って決定する手段を持つ方法は、各条件により得られる結果に対する影響が変化する複雑・非線形な対象に対する予測精度は良いが、熟練者がそのルールを入力する必要があり、予測モデルのメンテナンス性に問題がある。
【0009】
本発明は、前記従来の問題点を解決するべくなされたもので、条件空間の位置に応じて、各条件の結果に対する影響が変化する複雑・非線形な対象であっても、特別なルールを入力することなく、任意の要求条件に対する結果を高精度に予測することができる結果予測装置を提供することを課題とする。
【0010】
【課題を解決するための手段】
本発明は、過去の条件の値と、その条件によって得られた結果とを保存した実績データベースと、実績データベースに保存されている条件により規定される条件空間において、結果を予測したい要求条件の近傍における各条件の結果に対する影響係数を計算する手段と、得られた影響係数に基づいて条件空間の軸を変換し、変換された条件空間において、前記実績データベースに保存されている過去の条件の値と前記要求条件との距離を計算する手段と、得られた距離に基づいて、各条件の値と前記要求条件との類似度を計算する手段と、得られた類似度に基づいて、前記要求条件近傍の予測式を作成する手段と、得られた予測式に基づいて、要求条件に対する結果を計算する手段とを備えたことにより、前記課題を解決したものである。
【0011】
即ち、本発明においては、結果を予測したい要求条件から実績データとして保存されている各条件までの距離を、該要求条件の近傍における各条件の結果に対する影響係数に基づいて軸変換した条件空間において計算するようにしたので、該距離に基づいて特別なルールを入力することなく、複雑・非線形な対象についても、要求条件に対する結果を高精度に予測することができる。
【0012】
本発明は、又、前記影響係数が、前記要求条件の近傍における各条件と結果との関係を近似する線形式における各条件に対する係数であるようにすることができる。
【0013】
本発明は、又、予測計算する前記結果が、鋼材の材質であるようにすることができる。
【0014】
【発明の実施の形態】
以下、図面を参照して、本発明の実施の形態について詳細に説明する。
【0015】
図1は、本発明に係る一実施形態の結果予測装置の要部構成を示すブロック図である。
【0016】
本実施形態の結果予測装置は、過去に適用した条件の値と、その条件によって得られた結果とを保存した実績データベース10と、該実績データベース10に保存されている条件により規定される条件空間において、結果を予測したい要求条件の近傍における各条件が、結果に対して影響する程度を表わす影響係数を計算する影響係数計算部12と、得られた影響係数に基づいて、条件空間の軸を変換し、変換された条件空間において、前記実績データベースに保存されている過去の条件の観測値と前記要求条件との距離を計算する距離計算部14と、得られた距離に基づいて、各条件の観測値と前記要求条件との類似度を計算する類似度計算部16と、得られた類似度に基づいて、前記要求条件近傍の予測式を作成する予測式作成部18と、得られた予測式に基づいて、要求条件に対する結果を計算する予測計算部20とを備えている。以下、本実施形態の予測装置について詳述する。
【0017】
前記実績データベース10には、図2のテーブルで示すように、実績データとして、過去に適用された条件であるM個の入力変数と、これらの入力変数の組合せにより得られた(観測された)結果である出力変数からなるN個の観測データが、予め保存されているとする。このようにデータベースで与えられる観測データとしては、例えば製鉄の場合であれば、鉄鋼の素材成分や操業条件を入力変数(条件)とし、鉄鋼の材質、例えば強度を出力変数(結果)とする例を挙げることができる。
【0018】
ここでは、図示してあるように、出力変数の項目名称をY、M個の入力変数の項目名称をXm(m=1,2,…,M)とする。観測データはN個あり、n番目(n=1,2,…,N)の出力変数の値をyとし、入力変数の値をx と表記することにする。
【0019】
前記影響係数計数部12では、図2に示したM個の入力変数により規定される条件空間において、要求条件の近傍に位置する各条件について、過去に得られている結果に対する影響係数を計算により求める。ここでは、結果を予測したい要求条件を入力ベクトルとし、これを
=[x ,x ,…,x …(1)
で表記する。
【0020】
まず、大域的な回帰式のパラメータを推定する。即ち、図2で、与えられたN個の観測データを用いて、結果(Y)を予測するための回帰式モデルを作成し、該回帰式のパラメータを推定する。
【0021】
モデル式は次の線形式
Y=b+a・X+a・X+…+a・X …(2)
とし、この回帰式のパラメータ:b,a,a,…,aを最小2乗法により求める。
【0022】
このパラメータから定数bを除いて係数のみを抽出した次式の偏回帰係数ベクトル
α=[a,a,…,a …(3)
を、次に説明する距離計算に用いる影響係数とする。
【0023】
前記距離計算部14では、各観測データの入力変数について、前記要求条件からの距離計算を行なう。そのために、まず入力空間(条件空間)のある点x=[x,x,…,xに対する、前記(1)式の要求条件xからの距離Lを計算するための距離関数を、前記(3)式の影響係数を考慮した次式
【数1】

Figure 2004355189
により定義する。
【0024】
この(4)式では、各変数と要求条件との差の絶対値に、それぞれ影響係数aの絶対値を掛けたものを、全ての変数について足し合わせる処理を行なっている。
【0025】
前記(3)式で与えられる偏回帰係数(影響係数)aは、出力変数Yの変化量に対する各入力変数Xの寄与度と考えることができる。従って、上記(4)式の距離関数は、その寄与度を加味した重み付きの距離を表わしていることになる。
【0026】
又、この距離関数により距離を計算することは、同時にこの影響係数により条件空間における軸変換の操作を実行していることになる。これを、便宜上、X、Xの2次元に対する出力変数Yの場合の観測データの空間におけるデータ分布のイメージが、図3(A)に破線で囲んだ点で示すようであるとし、このデータ分布における要求条件近傍の回帰式が、
Y=b+a+a …(2´)
で表わされるとすると、同図(B)に示すように、影響係数a、aを用いて軸をX/|a|、X/|a|に変換し、この軸変換された空間におけるxとxとの距離Lを計算していることになる。因みに、正規化ユークリッド距離の場合は、各変数に対応する条件軸をそれぞれのデータ分布の標準偏差で割っているが、ここでは係数で割っている。
【0027】
次いで、前記(4)式で定義した距離関数を用いて、各観測データの要求条件からの距離を計算する。即ち、前記図2に示したN個の観測データのそれぞれについて、要求条件xからの距離を求める。
【0028】
具体的には、n番目(n=1,2,…,N)の観測データxの要求条件からの距離は、次の式
=L(x,x,α) …(5)
ここで、x=[x ,x ,…,x
n=1,2,…,N
から求めることができる。又、1〜N番目の観測データについて計算された要求条件からの各距離をまとめて次式
l=[L,L,…,L …(6)
のように表記する。
【0029】
前記類似度計算部16では、以上のように、対象とする全ての観測データについて、要求条件からの距離計算を実行した後、各観測データの要求条件からの類似度を計算する。そのために、まず要求条件からの近さを表わす類似度関数Wを、次式
W(L,p,l)=exp{−(L/(p・σ(l)))} …(7)
ここで、σ(l):正規化に使用するlの標準偏差
p:調整パラメータ(初期値:1.5)
のように定義する。
【0030】
図4には、この類似度関数の特徴を示す。即ち、前記(5)式により得られる各観測データの要求条件からの距離が短いほど類似度が高く、長いほど低い値をとる。なお、類似度関数はこれに限定されず、同様の特徴を持つ、例えば折れ線関数としても、あるいは、前記非特許文献1に記載されているトリキューブ関数を用いてもよい。
【0031】
次に、上記のように定義した類似度関数を用いて、各観測データの要求条件からの類似度を計算する。即ち、前記図2のN個の観測データそれぞれについて、前記(5)式により計算された距離を用いて要求条件からの類似度を求める。
【0032】
n番目(n=1,2,…,N)の観測データの要求条件からの類似度は、次の式
=W(L,p,l) …(8)
(n=1,2,…,N)
から求めることができる。又、ここでは、1〜N番目の観測データの要求条件からの類似度を求めて次式
w=[W,W,…,W …(9)
のように表記する。
【0033】
前記予測式作成部18では、以上のように全ての観測データについて要求条件からの類似度の計算が終了した後、局所回帰式のパラメータを推定計算し、与えられたN個の観測データと、それぞれの類似度wを用いて、回帰式モデルを作成する。
【0034】
モデル式は、次の線形式
Y=b+a・X+a・X+…+a・X …(10)
とする。この式が、要求条件の結果を予測するために使用する最終的な予測式である。
【0035】
便宜上、この予測式(10)は、前記(2)式の線形式と同一式で表わされているが、この(10)式では、パラメータθ=[b,a,a,…,aを、類似度wを重みとする重み付き最小2乗法により求める。
【0036】
このようにすることにより、類似度の大きい観測データ(要求点(条件)に近いデータ)は、重みが大きく、類似度の小さい観測データ(要求点から遠いデータ)は、重みが小さくなるような回帰式が得られ、要求条件の近傍のデータをより精度良くフィッティングする回帰式モデルができる。
【0037】
ここに、(10)式の局所回帰式と前述した(2)式の大域的な回帰式との差異を説明する。局所回帰式と大域的な回帰式は、いずれも実績データベース10に蓄積されているすべての観測データを用いて、パラメータを最小2乗法を用いて推定することにより求めるが、大域的回帰式(2)は、すべての観測データの重みを等しくして、最小2乗法によりパラメータを推定しているため、どの要求条件においても、パラメータは同じ値になり、製造条件空間すべてにおいて共通な、即ち大域的に使用できる回帰式である。
【0038】
これに対し、局所回帰式(10)は、要求条件に近い観測データの重みを大きくして、遠い観測データの重みを小さくして、最小2乗法によりパラメータを推定しているため、要求条件の値によって、パラメータの値は異なり、局所的にしか使用できない(有効でない)が、精度の高い回帰式である。
【0039】
前記予測計算部20では、要求条件の値を、上で得られた局所回帰式(10)の右辺に与えて、結果予測値を計算する。
【0040】
以上詳述した本実施形態の予測装置による予測処理の手順の概要を図5に示す。
【0041】
ステップ1では、前記影響係数計算部12により影響係数の計算を実行し、ステップ2では前記距離計算部14と類似度計算部16により距離と類似度の計算を実行し、ステップ3では予測式作成部18によりステップ3のパラメータを推定して予測式を作成し、ステップ4では前記予測計算部20によりステップ4の予測計算を実行し、出力している。
【0042】
本発明の方法を、薄鋼板の強度推定に適用したところ、図6に示すような実績値に対する予測値の結果が得られた。なお、この図の単位MPaはメガパスカルである。
【0043】
この図6に示した予測精度評価は、100個の熱延鋼板の観測データを基に評価を行った結果であり、100個の中の1つの観測データの例を次に示す。
【0044】
板厚:15.66[mm]、板幅:1257[mm]、
化学成分[%]
C:0.063、Si:0.19、Mn:1.44、P:0.018、S:0.0023、Al:0.021、
Nb:0.04、V:0.019、Ti:0.008、Cu:0.01、Ni:0.01、Cr:0.03、
Ca:0.0001、N:0.0038、O:0.0036、Mo:0.001、B:0.001,
加熱炉抽出温度:1191[℃]、仕上ミル前面温度:948[℃]、仕上ミル後面温度:824[℃]、巻取温度:519[℃]
【0045】
この図6の結果より、1σ=9.9[MPa]の精度が得られ、前記特許文献1に開示されている従来方法(条件空間のマハラノビス距離を用いる場合)に比べて、推定誤差の標準偏差を35%低減できた。
【0046】
又、特許文献2の方法のように、入力変数を限定するツールを熟練者が入力する必要もなく、実績データから自動的に要求条件近傍における各入力変数の重みを計算することができたので、予測装置のメンテナンス性を向上することもできた。
【0047】
なお、結果予測装置の具体的な構成は、前記図1に示したものに限定されない。又、本発明を適用する具体例として、鋼材の材質の一つである薄鋼板の強度を推定する場合を示したが、他の材質であってもよいことはいうまでもなく、更には複数の条件とその条件を適用したときの結果の間に因果関係があり、又、実績データ間に実質上の連続性が成立するような対象であれば、任意のものに適用可能である。
【0048】
【発明の効果】
以上説明したとおり、本発明によれば、条件空間の位置により、各条件の結果に対する影響が変化する複雑・非線形な対象であっても、特別なルールを入力することなく、要求条件に対する結果を高精度に予測することができる。
【図面の簡単な説明】
【図1】本発明に係る一実施形態の予測装置の要部を示すブロック図
【図2】実績データベースに保存されているデータのイメージを示す図表
【図3】影響係数による条件空間の軸変換のイメージを示す線図
【図4】類似度関数の一例の特徴を示す線図
【図5】実施形態の予測装置における処理手順の概要を示すフローチャート
【図6】本発明方法を薄鋼板の強度推定に適用した結果を示す線図
【符号の説明】
10…実績データベース
12…影響係数計算部
14…距離計算部
16…類似度計算部
18…予測式作成部
20…予測計算部[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a result predicting apparatus, and more particularly to a method for predicting a result for an arbitrary required condition using a database storing conditions applied in the past and results obtained by applying the condition as actual data. And a suitable result prediction device.
[0002]
[Prior art]
As a method of predicting a result with respect to an arbitrary request condition, for example, a distance between a condition of each past observation data stored in an actual result database and a desired request condition is calculated, and the observation data ( There is known a method of calculating a weight of actual data, creating a function for fitting a neighborhood of a requirement from the weight, and predicting a result based on the requirement using the function (for example, Non-Patent Document 1). , 2, Patent Document 1).
[0003]
In calculating the distance between an arbitrary required condition and each condition of past performance data, a norm in a condition space defined on an axis having each condition as a variable is used. Non-Patent Documents 1 and 2 describe Euclidean distance and normalized Euclidean distance, which are a kind of norm, and Patent Document 1 describes using Mahalanobis distance.
[0004]
In particular, Patent Literature 1 has a problem when using the Euclidean distance because the normalized Euclidean distance is normalized so that the axes of the individual condition spaces (condition axes) have a normal distribution. It is described that the difference in distribution depending on the item has no effect. When using the Mahalanobis distance, in addition to the normalized Euclidean distance, the direction of the axis is newly re-oriented in order to remove the correlation between items, and the axes are mutually uncorrelated. It is described that they can be distinguished and can be fitted better.
[0005]
Further, Patent Document 2 discloses a prediction method for an object in which a factor affecting a result changes depending on a position of a condition space. This is because the axis of the condition space is determined according to the rule according to the position of the condition space, so that the result with respect to the required condition can be accurately predicted even for a complex / non-linear object.
[0006]
[Patent Document 1]
JP 2001-290508 A [Patent Document 2]
Japanese Patent Application Laid-Open No. 2002-236119 [Non-Patent Document 1]
William S.M. Cleveland and Susan J. et al. Devlin;
Locally Weighted Regression: An Approach to
Regression Analysis by Local Fitting, Journal of the American Statistical Association, Vol. 83, No. 403, September 1988.
[Non-patent document 2]
Zheng, Q. et al. , And H .; Kimura; Locally Weighted
Regression Based on k Bipartite Neighbors,
Proc. 42nd Japan Joint Automatic Control
Conference, pp143, Tokyo (1999).
[0007]
[Problems to be solved by the invention]
However, the method described in Non-Patent Documents 1 and 2 and Patent Document 1, that is, a method of linearly converting the condition space of the actual data (the Euclidean distance uses the conditional axis of the actual data as it is. The normalized Euclidean distance is The condition axis of the actual data is scaled by the distribution (standard deviation) of each condition. The Mahalanobis distance is linearly transformed so that the condition axes of the actual data are uncorrelated with each other.) There is a limit to the prediction accuracy for a complex / non-linear object whose influence on the result obtained by each condition changes.
[0008]
Further, the method described in Patent Document 2, that is, a method having means for determining the axis of the condition space according to the rule according to the position of the condition space, has a complicated and nonlinear effect in which the effect on the result obtained by each condition changes. Although the prediction accuracy for the target is good, it is necessary for an expert to input the rule, and there is a problem in the maintainability of the prediction model.
[0009]
The present invention has been made to solve the above-mentioned conventional problems, and a special rule is input even for a complex / non-linear object whose influence on the result of each condition changes according to the position of the condition space. It is an object of the present invention to provide a result predicting apparatus capable of predicting a result with respect to any required condition with high accuracy without performing.
[0010]
[Means for Solving the Problems]
The present invention provides a performance database storing values of past conditions and results obtained under the conditions, and a condition space defined by the conditions stored in the performance database, in a condition space defined by the conditions, in the vicinity of a request condition whose result is to be predicted. Means for calculating an influence coefficient with respect to the result of each condition in, and converting the axis of the condition space based on the obtained influence coefficient, and in the converted condition space, the value of the past condition stored in the performance database Means for calculating a distance between the request and the required condition; means for calculating a similarity between the value of each condition and the required condition on the basis of the obtained distance; and The object has been achieved by providing means for creating a prediction formula near a condition and means for calculating a result for a required condition based on the obtained prediction formula.
[0011]
That is, in the present invention, the distance from the required condition for which the result is to be predicted to each condition stored as actual data is calculated in a condition space in which the axis is transformed based on the influence coefficient on the result of each condition in the vicinity of the required condition. Since the calculation is performed, the result with respect to the required condition can be predicted with high accuracy even for a complex / non-linear object without inputting a special rule based on the distance.
[0012]
According to the present invention, the influence coefficient may be a coefficient for each condition in a linear form that approximates a relation between each condition and a result near the required condition.
[0013]
According to the present invention, the result of the predictive calculation may be a material of a steel material.
[0014]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
[0015]
FIG. 1 is a block diagram illustrating a main configuration of a result prediction device according to an embodiment of the present invention.
[0016]
The result prediction apparatus according to the present embodiment includes a result database 10 storing values of conditions applied in the past and results obtained under the conditions, and a condition space defined by the conditions stored in the results database 10. In the above, an influence coefficient calculation unit 12 that calculates an influence coefficient indicating the degree to which each condition in the vicinity of a required condition for which a result is to be predicted influences the result, and an axis of a condition space based on the obtained influence coefficient In the converted condition space, the distance calculation unit 14 that calculates the distance between the observed value of the past condition stored in the performance database and the required condition, and, based on the obtained distance, each condition A similarity calculator 16 for calculating the similarity between the observed value of the above and the required condition, and a prediction formula creator 18 for creating a prediction formula near the required condition based on the obtained similarity. Based on the obtained prediction equation, and a prediction calculating unit 20 to calculate the result for the request condition. Hereinafter, the prediction device of the present embodiment will be described in detail.
[0017]
As shown in the table of FIG. 2, the result database 10 obtained (observed) as result data by M input variables which are conditions applied in the past and a combination of these input variables. It is assumed that N observation data including output variables as results are stored in advance. As the observation data provided in the database as described above, for example, in the case of steelmaking, an example in which the material composition and operating conditions of steel are used as input variables (conditions) and the material of steel, for example, strength is used as an output variable (result). Can be mentioned.
[0018]
Here, as shown, the item names of the output variables are Y, and the item names of the M input variables are Xm (m = 1, 2,..., M). Observation data exactly N, n-th (n = 1,2, ..., N ) the values of the output variables and y n, the value of the input variable to be referred to as x m n.
[0019]
The influence coefficient counting unit 12 calculates an influence coefficient for a result obtained in the past for each condition located in the vicinity of the required condition in the condition space defined by the M input variables shown in FIG. Ask. Here, the requirements to be predictive of results as input vector, which x r = [x 1 r, x 2 r, ..., x M r] T ... (1)
Notation.
[0020]
First, the parameters of the global regression equation are estimated. That is, in FIG. 2, a regression equation model for predicting the result (Y) is created using the given N pieces of observation data, and parameters of the regression equation are estimated.
[0021]
The model formula is represented by the following linear form Y = b + a 1 · X 1 + a 2 · X 2 +... + A M · X M (2)
The parameters b, a 1 , a 2 ,..., A M of the regression equation are obtained by the least square method.
[0022]
A partial regression coefficient vector α = [a 1 , a 2 ,..., A M ] T ...
Is an influence coefficient used for distance calculation described below.
[0023]
The distance calculator 14 calculates a distance from the required condition for the input variable of each observation data. Therefore, first, the input space x = [x 1, x 2 , ..., x M] point in the (condition-space) for the T, the distance to calculate the distance L from the requirements x r of the equation (1) The function is expressed by the following equation taking into account the influence coefficient of the above equation (3).
Figure 2004355189
Defined by
[0024]
In this equation (4), the absolute value of the difference between each variable and requirements, a multiplied by the absolute value of each influence coefficient a m, is performed a process of adding up all the variables.
[0025]
Wherein (3) the partial regression coefficient (influence coefficient) a m given by equation can be considered as the contribution of each input variable X m corresponding to the change in the output variable Y. Therefore, the distance function of the above equation (4) represents a weighted distance in consideration of the contribution.
[0026]
Calculating the distance with this distance function means that the axis conversion operation in the condition space is being executed by the influence coefficient at the same time. For convenience, it is assumed that the image of the data distribution in the space of the observation data in the case of the output variable Y with respect to the two dimensions of X 1 and X 2 is as shown by a point surrounded by a broken line in FIG. The regression equation near the requirement in the data distribution is
Y = b + a 1 X 1 + a 2 X 2 ( 2 ′ )
As shown in FIG. 4B, the axes are converted into X 1 / | a 1 | and X 2 / | a 2 | using the influence coefficients a 1 and a 2 , and this axis conversion is performed. it means that by calculating the distance L between x and x r in space. Incidentally, in the case of the normalized Euclidean distance, the condition axis corresponding to each variable is divided by the standard deviation of each data distribution. Here, the condition axis is divided by a coefficient.
[0027]
Next, the distance from the required condition of each observation data is calculated using the distance function defined by the above equation (4). That is, for each of the N observation data shown in FIG. 2, obtaining the distance from the requirements x r.
[0028]
Specifically, n-th (n = 1,2, ..., N ) Distance from requirements of observation data x n of the following expression L n = L (x n, x r, α) ... (5 )
Here, x n = [x 1 n , x 2 n, ..., x M n] T
n = 1, 2,..., N
Can be obtained from Also, the distances from the required conditions calculated for the first to Nth observation data are put together and the following equation l = [L 1 , L 2 ,..., L N ] T (6)
Notation as follows.
[0029]
As described above, the similarity calculation unit 16 calculates the similarity of each observation data from the required conditions after executing the distance calculation from the required conditions for all the target observation data. For this purpose, first, a similarity function W representing proximity from a required condition is expressed by the following equation: W (L, p, l) = exp {− (L / (p · σ (l))) 2 } (7)
Here, σ (l): standard deviation of l used for normalization p: adjustment parameter (initial value: 1.5)
Is defined as
[0030]
FIG. 4 shows the characteristics of this similarity function. That is, the similarity increases as the distance from the required condition of each observation data obtained by the equation (5) decreases, and the similarity decreases as the distance increases. Note that the similarity function is not limited to this, and may have a similar feature, for example, a polygonal function, or a tricube function described in Non-Patent Document 1.
[0031]
Next, using the similarity function defined as above, the similarity from the required condition of each observation data is calculated. That is, for each of the N pieces of observation data shown in FIG. 2, the similarity from the required condition is obtained using the distance calculated by the equation (5).
[0032]
n-th (n = 1,2, ..., N ) similarity from requirements of observation data, the following equation W n = W (L n, p, l) ... (8)
(N = 1, 2,..., N)
Can be obtained from Here, the similarity of the first to N-th observation data from the required conditions is obtained, and the following expression is obtained: w = [W 1 , W 2 ,..., W N ] T (9)
Notation as follows.
[0033]
After the calculation of the similarity from the required condition is completed for all the observation data as described above, the prediction formula creation unit 18 estimates and calculates the parameters of the local regression equation, and gives the given N observation data, A regression equation model is created using each similarity w.
[0034]
The model formula is represented by the following linear form Y = b + a 1 · X 1 + a 2 · X 2 +... + A M × X M (10)
And This is the final prediction formula used to predict the outcome of the requirement.
[0035]
For convenience, this prediction equation (10) is expressed by the same equation as the linear form of the above equation (2), but in this equation (10), the parameters θ = [b, a 1 , a 2 ,. a M ] T is obtained by a weighted least squares method using the similarity w as a weight.
[0036]
In this way, observation data with a high degree of similarity (data close to a request point (condition)) has a large weight, and observation data with a low degree of similarity (data far from the request point) has a low weight. A regression equation is obtained, and a regression equation model that fits data near the required conditions with higher accuracy can be created.
[0037]
Here, the difference between the local regression equation of equation (10) and the global regression equation of equation (2) will be described. The local regression equation and the global regression equation are both obtained by estimating parameters using the least squares method using all the observation data accumulated in the performance database 10, but the global regression equation (2 ) Estimates the parameters by the least-squares method with the weights of all the observation data equal, so that the parameters have the same value in any of the required conditions, and are common in all manufacturing condition spaces, that is, global. This is a regression equation that can be used for
[0038]
On the other hand, the local regression equation (10) estimates the parameters by the least squares method by increasing the weight of the observation data close to the required condition and decreasing the weight of the distant observation data. Depending on the value, the value of the parameter differs and can only be used locally (not valid), but is a highly accurate regression equation.
[0039]
The prediction calculation unit 20 calculates the result prediction value by giving the value of the required condition to the right side of the local regression equation (10) obtained above.
[0040]
FIG. 5 shows an outline of the procedure of the prediction process by the prediction device of the present embodiment described in detail above.
[0041]
In step 1, the influence coefficient is calculated by the influence coefficient calculation unit 12, in step 2, the distance and similarity are calculated by the distance calculation unit 14 and the similarity calculation unit 16, and in step 3, a prediction formula is created. The prediction equation is created by estimating the parameters in step 3 by the unit 18, and the prediction calculation in step 4 is executed and output by the prediction calculation unit 20 in step 4.
[0042]
When the method of the present invention was applied to the strength estimation of a thin steel sheet, a result of a predicted value with respect to an actual value as shown in FIG. 6 was obtained. The unit MPa in this figure is megapascal.
[0043]
The prediction accuracy evaluation shown in FIG. 6 is a result of evaluation based on observation data of 100 hot-rolled steel sheets, and an example of one observation data out of 100 is shown below.
[0044]
Board thickness: 15.66 [mm], board width: 1257 [mm],
Chemical composition[%]
C: 0.063, Si: 0.19, Mn: 1.44, P: 0.018, S: 0.0023, Al: 0.021,
Nb: 0.04, V: 0.019, Ti: 0.008, Cu: 0.01, Ni: 0.01, Cr: 0.03,
Ca: 0.0001, N: 0.0038, O: 0.0036, Mo: 0.001, B: 0.001,
Heating furnace extraction temperature: 1191 [° C], finishing mill front surface temperature: 948 [° C], finishing mill rear surface temperature: 824 [° C], winding temperature: 519 [° C]
[0045]
From the result of FIG. 6, an accuracy of 1σ = 9.9 [MPa] is obtained, and compared with the conventional method disclosed in Patent Document 1 (in the case of using the Mahalanobis distance in the condition space), the standard error of the estimation error is obtained. The deviation could be reduced by 35%.
[0046]
Further, unlike the method of Patent Document 2, it is not necessary for an expert to input a tool for limiting input variables, and the weight of each input variable in the vicinity of a required condition can be automatically calculated from actual data. In addition, the maintainability of the prediction device could be improved.
[0047]
The specific configuration of the result prediction device is not limited to that shown in FIG. In addition, as a specific example to which the present invention is applied, the case where the strength of a thin steel plate, which is one of the steel materials, is estimated, but it is needless to say that other materials may be used. There is a causal relationship between the above condition and the result when the condition is applied, and the present invention can be applied to any object as long as the actual data has substantial continuity.
[0048]
【The invention's effect】
As described above, according to the present invention, even for a complex / non-linear object whose influence on the result of each condition changes depending on the position of the condition space, the result for the required condition can be obtained without inputting a special rule. Prediction can be made with high accuracy.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a main part of a prediction apparatus according to an embodiment of the present invention. FIG. 2 is a table showing an image of data stored in an actual result database. FIG. FIG. 4 is a diagram showing characteristics of an example of a similarity function. FIG. 5 is a flowchart showing an outline of a processing procedure in a prediction device according to an embodiment. Diagram showing results applied to estimation [Explanation of symbols]
DESCRIPTION OF SYMBOLS 10 ... Result database 12 ... Influence coefficient calculation part 14 ... Distance calculation part 16 ... Similarity calculation part 18 ... Prediction formula preparation part 20 ... Prediction calculation part

Claims (3)

過去の条件の値と、その条件によって得られた結果とを保存した実績データベースと、
実績データベースに保存されている条件により規定される条件空間において、結果を予測したい要求条件の近傍における各条件の結果に対する影響係数を計算する手段と、
得られた影響係数に基づいて条件空間の軸を変換し、変換された条件空間において、前記実績データベースに保存されている過去の条件の値と前記要求条件との距離を計算する手段と、
得られた距離に基づいて、各条件の値と前記要求条件との類似度を計算する手段と、
得られた類似度に基づいて、前記要求条件近傍の予測式を作成する手段と、
得られた予測式に基づいて、要求条件に対する結果を計算する手段とを備えたことを特徴とする結果予測装置。
A performance database that stores the values of past conditions and the results obtained under those conditions,
Means for calculating an influence coefficient for the result of each condition in the condition space defined by the conditions stored in the performance database, in the vicinity of the required condition whose result is to be predicted;
Converting the axis of the condition space based on the obtained influence coefficient, means for calculating the distance between the value of the past condition stored in the performance database and the required condition in the converted condition space,
Means for calculating a similarity between the value of each condition and the required condition based on the obtained distance,
Means for creating a prediction formula near the required condition based on the obtained similarity;
Means for calculating a result for a required condition based on the obtained prediction formula.
前記影響係数が、前記要求条件の近傍における各条件と結果との関係を近似する線形式における各条件に対する係数であることを特徴とする請求項1に記載の結果予測装置。The result prediction apparatus according to claim 1, wherein the influence coefficient is a coefficient for each condition in a linear form that approximates a relation between each condition and a result near the required condition. 予測計算する前記結果が、鋼材の材質であることを特徴とする請求項1又は2に記載の結果予測装置。3. The result prediction apparatus according to claim 1, wherein the result of the prediction calculation is a material of a steel material.
JP2003150345A 2003-05-28 2003-05-28 Result prediction device Expired - Lifetime JP4474847B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003150345A JP4474847B2 (en) 2003-05-28 2003-05-28 Result prediction device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003150345A JP4474847B2 (en) 2003-05-28 2003-05-28 Result prediction device

Publications (2)

Publication Number Publication Date
JP2004355189A true JP2004355189A (en) 2004-12-16
JP4474847B2 JP4474847B2 (en) 2010-06-09

Family

ID=34046173

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003150345A Expired - Lifetime JP4474847B2 (en) 2003-05-28 2003-05-28 Result prediction device

Country Status (1)

Country Link
JP (1) JP4474847B2 (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006309709A (en) * 2005-03-30 2006-11-09 Jfe Steel Kk Result prediction device, control device and quality design device
WO2007080688A1 (en) 2006-01-13 2007-07-19 Jfe Steel Corporation Prediction formula making device and prediction formula making method
JP2008112288A (en) * 2006-10-30 2008-05-15 Jfe Steel Kk Prediction type creation device, result prediction device, quality design device, prediction type creation method and method for manufacturing product
JP2009167481A (en) * 2008-01-17 2009-07-30 Jfe Steel Corp Apparatus for controlling impeller desulfuration and method therefor
JP2009230412A (en) * 2008-03-21 2009-10-08 Jfe Steel Corp Result predicting device, and product quality prediction method using same
JP2009245141A (en) * 2008-03-31 2009-10-22 Kobe Steel Ltd Output value prediction method, device thereof, and program for the method
JP2010007992A (en) * 2008-06-27 2010-01-14 Jfe Steel Corp Sintering apparatus and sintering method
JP2011039763A (en) * 2009-08-11 2011-02-24 Kobe Steel Ltd Output value prediction method, output value prediction device, and program for the method
JP2015041321A (en) * 2013-08-23 2015-03-02 大日本印刷株式会社 Information processing apparatus, display system, and program
JP2020128975A (en) * 2019-02-07 2020-08-27 Jfeスチール株式会社 Method for generating prediction model, method for predicting corrosion amount of metallic material, program for generating prediction model, and device of generating prediction model
WO2021014804A1 (en) 2019-07-22 2021-01-28 Jfeスチール株式会社 Method for generating quality prediction model, quality prediction model, quality prediction method, method for manufacturing metal material, device for generating quality prediction model, and quality prediction device
WO2022149344A1 (en) 2021-01-06 2022-07-14 Jfeスチール株式会社 Quality abnormality analysis method, metal material manufacturing method, and quality abnormality analysis device
JP7483827B1 (en) 2022-10-27 2024-05-15 エヌ・ティ・ティ・コミュニケーションズ株式会社 Learning device, learning method, and learning program

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006309709A (en) * 2005-03-30 2006-11-09 Jfe Steel Kk Result prediction device, control device and quality design device
US8374981B2 (en) 2006-01-13 2013-02-12 Jfe Steel Corporation Apparatus and method for constructing prediction model
WO2007080688A1 (en) 2006-01-13 2007-07-19 Jfe Steel Corporation Prediction formula making device and prediction formula making method
KR101011546B1 (en) * 2006-01-13 2011-01-27 제이에프이 스틸 가부시키가이샤 Prediction formula making device and prediction formula making method
CN101361085B (en) * 2006-01-13 2012-07-04 杰富意钢铁株式会社 Result prediction device and method
JP2008112288A (en) * 2006-10-30 2008-05-15 Jfe Steel Kk Prediction type creation device, result prediction device, quality design device, prediction type creation method and method for manufacturing product
JP2009167481A (en) * 2008-01-17 2009-07-30 Jfe Steel Corp Apparatus for controlling impeller desulfuration and method therefor
JP2009230412A (en) * 2008-03-21 2009-10-08 Jfe Steel Corp Result predicting device, and product quality prediction method using same
JP2009245141A (en) * 2008-03-31 2009-10-22 Kobe Steel Ltd Output value prediction method, device thereof, and program for the method
JP2010007992A (en) * 2008-06-27 2010-01-14 Jfe Steel Corp Sintering apparatus and sintering method
JP2011039763A (en) * 2009-08-11 2011-02-24 Kobe Steel Ltd Output value prediction method, output value prediction device, and program for the method
JP2015041321A (en) * 2013-08-23 2015-03-02 大日本印刷株式会社 Information processing apparatus, display system, and program
JP2020128975A (en) * 2019-02-07 2020-08-27 Jfeスチール株式会社 Method for generating prediction model, method for predicting corrosion amount of metallic material, program for generating prediction model, and device of generating prediction model
JP7020500B2 (en) 2019-02-07 2022-02-16 Jfeスチール株式会社 Prediction model generation method, corrosion amount prediction method for metal materials, prediction model generation program and prediction model generation device
WO2021014804A1 (en) 2019-07-22 2021-01-28 Jfeスチール株式会社 Method for generating quality prediction model, quality prediction model, quality prediction method, method for manufacturing metal material, device for generating quality prediction model, and quality prediction device
KR20220007653A (en) 2019-07-22 2022-01-18 제이에프이 스틸 가부시키가이샤 Quality prediction model generation method, quality prediction model, quality prediction method, metal material manufacturing method, quality prediction model generation apparatus and quality prediction apparatus
WO2022149344A1 (en) 2021-01-06 2022-07-14 Jfeスチール株式会社 Quality abnormality analysis method, metal material manufacturing method, and quality abnormality analysis device
KR20230116870A (en) 2021-01-06 2023-08-04 제이에프이 스틸 가부시키가이샤 Quality anomaly analysis method, metal material manufacturing method, and quality anomaly analysis device
JP7483827B1 (en) 2022-10-27 2024-05-15 エヌ・ティ・ティ・コミュニケーションズ株式会社 Learning device, learning method, and learning program

Also Published As

Publication number Publication date
JP4474847B2 (en) 2010-06-09

Similar Documents

Publication Publication Date Title
JP2004355189A (en) Result predicting device
JP6439780B2 (en) Magnetic property prediction device and magnetic property control device for electrical steel sheet
JP5265355B2 (en) Method and apparatus for optimizing flatness control in strip rolling process
KR101011546B1 (en) Prediction formula making device and prediction formula making method
JP6508185B2 (en) Result prediction device and result prediction method
JP2019074969A (en) Quality predicting device and quality predicting method
JP5867349B2 (en) Quality prediction apparatus, operation condition determination method, quality prediction method, computer program, and computer-readable storage medium
JP2019087152A (en) Manufacturing process state prediction device, method, and program, and manufacturing process control system
JP2008112288A (en) Prediction type creation device, result prediction device, quality design device, prediction type creation method and method for manufacturing product
JP6662109B2 (en) Product state prediction device, product state control device, product state prediction method and program
JP5789958B2 (en) Cooling stop temperature control device and cooling stop temperature control method
JP5167891B2 (en) Result prediction apparatus and product quality prediction method using the same
JP2014018844A (en) Heat transfer coefficient predictor for steel material and cooling control method
JP2009072807A (en) Method and device for controlling planar shape of rolled stock and method of manufacturing thick steel plate
JP2006309709A (en) Result prediction device, control device and quality design device
JP5577568B2 (en) Impeller desulfurization control apparatus and method
JP6020263B2 (en) Rolling load learning control device and learning control method, and metal plate manufacturing method using the same
JP2007033067A (en) Method of estimating dent rigidity
Li et al. Online Monitor and Control of Cooling Temperature on Run‐out Table of Hot Strip Mill
Gao et al. Using transient measurements in iterative steady-state optimizing control
JP2007260729A (en) Method and device for controlling edge drop
JP2010207900A (en) Method of predicting deformation resistance in hot rolling
JP6414086B2 (en) Rolling time prediction method and heating furnace extraction time determination method
JP2020071493A (en) Result prediction device, result prediction method and program
JP2023007165A (en) Characteristic estimation device of ferrous material, characteristic estimation model generation device of ferrous material, characteristic estimation method of ferrous material, characteristic estimation model generation method of ferrous material and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060327

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090421

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090622

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20100216

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100301

R150 Certificate of patent or registration of utility model

Ref document number: 4474847

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130319

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130319

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140319

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term