JP6863393B2 - Parameter value determination method, parameter value determination program and parameter value determination device - Google Patents
Parameter value determination method, parameter value determination program and parameter value determination device Download PDFInfo
- Publication number
- JP6863393B2 JP6863393B2 JP2018562775A JP2018562775A JP6863393B2 JP 6863393 B2 JP6863393 B2 JP 6863393B2 JP 2018562775 A JP2018562775 A JP 2018562775A JP 2018562775 A JP2018562775 A JP 2018562775A JP 6863393 B2 JP6863393 B2 JP 6863393B2
- Authority
- JP
- Japan
- Prior art keywords
- parameter value
- value
- signal output
- adjusted
- model
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 28
- 238000012937 correction Methods 0.000 claims description 80
- 238000012545 processing Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 description 8
- 230000006399 behavior Effects 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000012447 hatching Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
Description
本発明は、パラメータ値決定方法、パラメータ値決定プログラム及びパラメータ値決定装置に関する。 The present invention relates to a parameter value determination method, a parameter value determination program, and a parameter value determination device.
従来、光伝送装置などの信号を出力する製品の量産製造時においては、信号出力制御するファームウェアのパラメータ値を、信号出力が既定の許容範囲内に収まるような適切な値に設定する工程が行われている。 Conventionally, in mass production of products that output signals such as optical transmission devices, a process of setting the parameter value of the firmware that controls the signal output to an appropriate value so that the signal output falls within a predetermined allowable range is performed. It has been.
一方、特許文献1等には、調整対象(例えば電動機駆動装置)の調整時間を短縮するため、新規の調整対象の調整に必要な諸定数と調整仕様が同一又は類似する過去の調整対象を抽出し、抽出した過去の調整対象に基づいて調整の初期値を設定する技術が開示されている。 On the other hand, in Patent Document 1 and the like, in order to shorten the adjustment time of the adjustment target (for example, the electric motor drive device), the past adjustment targets having the same or similar adjustment specifications as the constants required for the adjustment of the new adjustment target are extracted. However, a technique for setting an initial value of adjustment based on the extracted past adjustment target is disclosed.
しかるに、信号を出力する製品の場合、出力信号は製品に搭載されているデバイスの素子値(個体特徴量)のばらつきの影響を大きく受けるため、過去に別の製品に設定したパラメータ値を流用することは難しい。このため、実際には、装置毎に適当なパラメータの初期値を設定して、出力信号を観察しながら、出力が許容範囲に入るようにパラメータ値を変更しつつ探索する必要がある。 However, in the case of a product that outputs a signal, the output signal is greatly affected by the variation in the element value (individual feature amount) of the device mounted on the product, so the parameter value set in another product in the past is diverted. It's difficult. Therefore, in practice, it is necessary to set an appropriate initial value of the parameter for each device, observe the output signal, and change the parameter value so that the output falls within the permissible range.
信号を出力する製品のパラメータ値を上述のように探索する場合、パラメータ値を変更すると信号出力が安定するまでに時間M(例えばM=1分)を要する。このため、信号出力が許容範囲に入るまでにパラメータ値をN回変更した場合には、M×N以上の調整時間を要することになる。現状では、調整時間は、20〜30分程度となっている。 When searching for the parameter value of a product that outputs a signal as described above, it takes time M (for example, M = 1 minute) for the signal output to stabilize when the parameter value is changed. Therefore, if the parameter value is changed N times before the signal output falls within the permissible range, an adjustment time of M × N or more is required. At present, the adjustment time is about 20 to 30 minutes.
1つの側面では、本発明は、信号出力を行う装置の調整時間を短縮することが可能なパラメータ値決定方法、パラメータ値決定プログラム及びパラメータ値決定装置を提供することを目的とする。 In one aspect, it is an object of the present invention to provide a parameter value determination method, a parameter value determination program, and a parameter value determination apparatus capable of shortening the adjustment time of the device that outputs a signal.
一つの態様では、パラメータ値決定方法は、パラメータ値に応じた信号出力を行う調整対象の装置に対して入力するパラメータ値を決定するパラメータ値決定方法であって、過去に調整対象となった装置の信号出力が予め定められた許容範囲の境界近傍の値を示したときに該装置に入力されていたパラメータ値と、該装置の個体特徴量と、の組み合わせのデータに基づいて、新たな調整対象の装置の信号出力が前記許容範囲の境界近傍の値を示すようなパラメータ値を該新たな調整対象の装置の個体特徴量から推定するためのモデルを設定し、前記新たな調整対象の装置の信号出力が前記許容範囲に入る確率を所定以上にするために前記モデルを用いて推定したパラメータ値を補正する補正値を、前記許容範囲内において予測されるパラメータ値に対する信号出力の変化に基づいて算出し、前記新たな調整対象の装置の個体特徴量と前記モデルとに基づいてパラメータ値を推定し、推定したパラメータ値を前記補正値で補正することで、前記新たな調整対象の装置に入力するパラメータ値を算出する、処理をコンピュータが実行するパラメータ値決定方法である。 In one embodiment, the parameter value determination method is a parameter value determination method for determining a parameter value to be input to an adjustment target device that outputs a signal according to the parameter value, and is a device that has been adjusted in the past. New adjustment based on the data of the combination of the parameter value input to the device when the signal output of the device indicates a value near the boundary of the predetermined allowable range and the individual characteristic amount of the device. A model for estimating a parameter value such that the signal output of the target device indicates a value near the boundary of the permissible range from the individual characteristic amount of the new adjustment target device is set, and the new adjustment target device is set. The correction value for correcting the parameter value estimated using the model in order to make the probability that the signal output of the above is within the allowable range exceeds a predetermined range is based on the change in the signal output with respect to the parameter value predicted within the allowable range. The parameter value is estimated based on the individual feature amount of the new adjustment target device and the model, and the estimated parameter value is corrected by the correction value to obtain the new adjustment target device. This is a parameter value determination method in which a computer executes a process for calculating a parameter value to be input.
信号出力を行う装置の調整時間を短縮することができる。 The adjustment time of the device that outputs the signal can be shortened.
以下、一実施形態に係る調整システムについて、図1〜図10に基づいて詳細に説明する。 Hereinafter, the adjustment system according to the embodiment will be described in detail with reference to FIGS. 1 to 10.
図1には、本実施形態に係る調整システム100の構成が概略的に示されている。図1の調整システム100は、調整対象の製品50と接続されるパラメータ値決定装置としての調整装置10と、製品50から出力される信号を検出する信号検出装置40と、を備える。調整対象の製品50は、例えば、信号を出力する光伝送装置などの装置である。
FIG. 1 schematically shows the configuration of the
製品50は、図3に示すように、信号発生器制御部52と、信号発生器54とを有する。信号発生器制御部52には、設定されたパラメータ値に応じて信号発生器54に信号を発生させるファームウェアが組み込まれている。信号発生器54は、信号発生器制御部52からの指示に基づいて信号を発生する。
As shown in FIG. 3, the
調整装置10は、製品50の信号発生器制御部52に設定するパラメータ値を探索する。具体的には、調整装置10は、信号検出装置40で検出される製品50の出力信号を許容範囲に含めることができるパラメータ値を探索する。
The adjusting
図2には、調整装置10のハードウェア構成が示されている。調整装置10は、図2に示すように、CPU(Central Processing Unit)90、ROM(Read Only Memory)92、RAM(Random Access Memory)94、記憶部(ここではHDD(Hard Disk Drive))96、入出力インタフェース97、及び可搬型記憶媒体用ドライブ99等を備えている。これら調整装置10の構成各部は、バス98に接続されている。調整装置10では、ROM92あるいはHDD96に格納されているプログラム(パラメータ値決定プログラムを含む)、或いは可搬型記憶媒体用ドライブ99が可搬型記憶媒体91から読み取ったプログラム(パラメータ値決定プログラムを含む)をCPU90が実行することにより、図3に示す各部の機能が実現される。なお、図3には、調整装置10のHDD96等に格納されている履歴DB(database)22、モデルパラメータDB24、補正値DB26、及び許容範囲DB28も図示されている。
FIG. 2 shows the hardware configuration of the
図3には、調整装置10及び製品50の機能ブロック図が示されている。図3に示すように、調整装置10は、CPU90がプログラムを実行することで、モデル設定部12、補正値算出部14、初期値決定部16、パラメータ値出力部18、としての機能が実現されている。
FIG. 3 shows a functional block diagram of the
モデル設定部12は、新たな調整対象の製品50に入力すべきパラメータ値の初期値(補正前)を決定するために用いるモデルのモデルパラメータを設定する。具体的には、モデル設定部12は、履歴DB22に格納されている調整済みの製品の調整データ(最終のパラメータ値(最終調整値)と個体特徴量)に基づいて、モデルパラメータを設定する。そして、モデル設定部12は、設定したモデルパラメータをモデルパラメータDB24に格納する。図4には、履歴DB22のデータ構造の一例が示されている。図4の履歴DB22には、調整済みの各製品の個体特徴量(X1,X2,…)が格納されているとともに、調整済みの各製品の最終調整値(y)が格納されている。なお、個体特徴量は、素子値とも呼ばれる。ここで、本実施形態では、過去の調整においては、所定のパラメータ値を初期値に定め、図5に示すように所定時間ごとにパラメータ値を変更しながら信号出力を繰り返す処理を行っている。そして、信号出力が許容範囲に入った段階で、そのときのパラメータ値を最終調整値として保存し、調整を終了することとしている。このように、信号出力が許容範囲に入った段階で調整を終了するのは、製品の製造時間短縮のためである。したがって、本実施形態では、許容範囲ぎりぎりの(許容範囲の境界に収束した)信号が出力されたときのパラメータ値が最終調整値として図4の履歴DB22に格納されている。なお、最終調整値yは、図4に示すように製品ごとにばらついている。なお、このばらつきは、個体特徴量の違いに起因するものである。
The
モデル設定部12は、履歴DB22に格納されている個体特徴量(X1,X2,…)と、調整済みの各製品の最終調整値(y)とに基づいて、次式(1)で示す線形モデルの定数a0及び係数a1、a2、…を算出し、算出結果をモデルパラメータとしてモデルパラメータDB24に格納する。
y=a0+a1×X1+a2×X2+… …(1)The
y = a 0 + a 1 x X1 + a 2 x X2 + ... (1)
なお、上式(1)に、新たな調整対象の製品50の個体特徴量(X1,X2,…)を代入することで、製品50の出力信号が許容範囲ぎりぎりの値を示す可能性が最も高いパラメータ値yを決定することができる。以下においては、パラメータ値yを「境界予測値」とも呼ぶものとする。
By substituting the individual feature quantities (X1, X2, ...) Of the
補正値算出部14は、履歴DB22のデータに基づいて、モデル設定部12が設定したモデルを用いて求められるパラメータ値(境界予測値)yを補正するための補正値を算出する。前述のように、境界予測値yは、製品50の出力信号が許容範囲ぎりぎりの値を示す可能性が最も高いパラメータ値である。したがって、補正値算出部14は、境界予測値yを、信号出力が許容範囲に含まれる可能性が最も高いパラメータ値に補正するための補正値gを算出する。ここで、補正後のパラメータ値y0は、y0=y+gと表すことができる。補正値算出部14は、算出した補正値gを補正値DB26に格納する。The correction
なお、補正値算出部14による補正値gの算出方法の詳細については、後述する。
The details of the method of calculating the correction value g by the correction
初期値決定部16は、モデル設定部12が設定したモデルのモデルパラメータをモデルパラメータDB24から読み出すとともに、補正値算出部14が算出した補正値を補正値DB26から読み出す。そして、初期値決定部16は、読み出したモデルパラメータと上式(1)とを用いて製品50の境界予測値yを算出するとともに、境界予測値yを補正値gで補正して、製品50に設定するのに適切なパラメータ値の初期値y0を算出する。初期値決定部16は、決定したパラメータ値の初期値y0をパラメータ値出力部18に送信する。なお、初期値決定部16は、製品50が調整装置10に接続された際に、製品50の個体特徴量を取得できるようになっている。なお、初期値決定部16は、モデルパラメータDB24からモデルを取得するモデル取得部、補正値DB26から補正値を取得する補正値取得部、モデルと補正値に基づいてパラメータ値の初期値を算出する算出部として機能している。The initial
パラメータ値出力部18は、初期値決定部16から受信したパラメータ値の初期値を製品50の信号発生器制御部52に対して出力する。また、信号検出装置40によって検出された製品50からの出力信号の値が、許容範囲DB28に格納されている許容範囲に入っているか否かを判断し、許容範囲に入っていない場合には、パラメータ値を変更して、信号発生器制御部52に対して再度出力する。なお、許容範囲DB28には、許容範囲を示す値(図7の例では、[z0,z1])が格納されている。The parameter
本実施形態では、上述のように、初期値決定部16において、モデル設定部12で設定されたモデルと補正値算出部14で算出された補正値gとを用いて、調整対象の製品50を調整する際のパラメータ値の初期値が決定される。この場合、決定されたパラメータ値の初期値は、許容範囲に含まれる可能性が最も高い。このため、パラメータ値出力部18が、決定されたパラメータ値の初期値を製品50に出力することで、許容範囲に含まれるパラメータ値を探索する時間を短くすることができる。これにより、製品50の調整を短時間で行うことが可能となる。
In the present embodiment, as described above, in the initial
次に、本実施形態の調整装置10の補正値算出部14による補正値の算出方法について、詳細に説明する。
Next, a method of calculating the correction value by the correction
本実施形態においては、許容範囲近傍での信号出力とパラメータ値の関係を示す曲線(収束曲線)は線形性を有すると仮定している。すなわち、許容範囲の境界に収束する際の挙動が、図6に示すように直線になると仮定することで、許容範囲内の挙動を予測する。図6において、許容範囲内の挙動は、次式(2)にて表すことができる。
z=c1×y+c2 …(2)In this embodiment, it is assumed that the curve (convergence curve) showing the relationship between the signal output and the parameter value in the vicinity of the permissible range has linearity. That is, the behavior within the permissible range is predicted by assuming that the behavior when converging to the boundary of the permissible range becomes a straight line as shown in FIG. In FIG. 6, the behavior within the permissible range can be expressed by the following equation (2).
z = c 1 × y + c 2 … (2)
なお、c1は、パラメータ値に対する感度(予測値)を意味し、c2は、定数(切片)を意味する。Note that c 1 means the sensitivity (predicted value) to the parameter value, and c 2 means the constant (intercept).
ここで、感度c1には、予測誤差があり、前述した式(1)を用いて算出される境界予測値yにも予測誤差がある。図7には、各予測誤差が図示されている。各予測誤差は、図7に示すように正規分布であるものとする。この場合、新たな調整対象の製品50の挙動を式(1)のモデルを用いて推定した場合、ハッチングを付して示す予測誤差が発生することになる。Here, the sensitivity c 1 has a prediction error, and the boundary prediction value y calculated by using the above-mentioned equation (1) also has a prediction error. FIG. 7 illustrates each prediction error. It is assumed that each prediction error has a normal distribution as shown in FIG. In this case, when the behavior of the
本実施形態では、補正値算出部14は、信号出力が許容範囲内に含まれる確率が最大となる値y0を算出するための補正値gを算出する。In the present embodiment, the correction
具体的には、許容範囲の境界へ収束するパラメータ値(境界予測値)yの予測誤差をσe 2、感度c1の平均をcave、感度の予測誤差をσc 2として、信号出力zに対する調整対象の製品50の存在確率を、補正値gの関数として表す。そして、許容範囲に入る確率が最大となるような補正値gを算出する。図8(b)は、図8(a)に示すようにパラメータ値を境界予測値とした場合において、信号出力が許容範囲に含まれる確率を、ハッチングを付した範囲の面積として示す図である。また、図9(b)は、図9(a)に示すようにパラメータ値を境界予測値から補正値gだけずらした場合において、信号出力が許容範囲に含まれる確率を、ハッチングを付した範囲の面積として示す図である。ここで、製品50の存在確率は、図8(b)、図9(b)に示すような正規分布の面積で表すことができる。本実施形態では、許容範囲内における製品50の存在確率を示す面積が最大となるような補正値gを算出する。ここで、許容範囲は[z0,z1]とする。Specifically, the signal output z is defined as the prediction error of the parameter value (boundary prediction value) y that converges to the boundary of the allowable range is σ e 2 , the average of the sensitivity c 1 is cave , and the sensitivity prediction error is σ c 2. The existence probability of the
この場合、c1〜N(cave,σc 2)とし、Δz=z1−z0とすると、Δzの平均E[Δz]及び分散Var[Δz]は、次式(3)、(4)にて表すことができる。In this case, if c 1 to N ( cave , σ c 2 ) and Δz = z 1 −z 0 , then the average E [Δz] and the variance Var [Δz] of Δz are given by the following equations (3) and (4). ) Can be expressed.
そして、次式(5)に上式(3)、(4)を代入することで、gを求めることができる。 Then, g can be obtained by substituting the above equations (3) and (4) into the following equation (5).
本実施形態では、補正値算出部14は、上式(5)を解くことで、信号出力が許容範囲に含まれる可能性を最大にすることができる補正値gを算出することができる。
In the present embodiment, the correction
(変形例)
上述した例では、許容範囲近傍での信号出力とパラメータ値との関係を示す曲線(収束曲線)が線形であると仮定して説明したが、非線形と仮定してもよい。非線形と仮定した場合、図10に示すように、新たに調整する製品50の境界予測値μyにおける収束曲線z=f(y)の傾きを用いて、近似的に、上述した線形と仮定した場合の手法を用いるようにすればよい。(Modification example)
In the above example, the curve (convergence curve) showing the relationship between the signal output and the parameter value in the vicinity of the permissible range is assumed to be linear, but it may be assumed to be non-linear. Assuming non-linearity, as shown in FIG. 10, using the slope of the convergence curve z = f (y) at the boundary prediction value μ y of the newly adjusted product 50, it was assumed to be approximately the above-mentioned linearity. The method of the case may be used.
具体的には、z=f(y)に対して、μyにおける次式(6)の近似を行う。Specifically, the following equation (6) in μ y is approximated to z = f (y).
そして、上式(6)のc1の平均caveと、分散σc 2を求め、上述した線形と仮定した場合の手法を用いて補正値gを求めるようにすればよい。Then, the average cave of c 1 and the variance σ c 2 of the above equation (6) may be obtained, and the correction value g may be obtained by using the above-mentioned method assuming linearity.
(平均caveと分散σc 2の求め方(収束曲線が2次多項式の場合))
ここで、収束曲線が2次多項式の場合には、c1の平均caveと、分散σc 2は、以下のようにして求めることができる。( How to find the mean cave and variance σ c 2 (when the convergence curve is a quadratic polynomial))
Here, when the convergence curve is a quadratic polynomial, the mean cave of c 1 and the variance σ c 2 can be obtained as follows.
2次元多項式をz=f(y)=(1/2)ay2+by+cとし、パラメータ(a,b,c)が既に推定されており、平均μa、μb、分散σa 2,σb 2は、既知であるとする。この場合、2次元多項式f(y)のパラメータ(a,b,c)をθと置き、f(y)=f(y;θ)とし、f(y;θ)をyで微分すると、f’(y;θ)=ay+bとなる。したがって、c1の平均caveと、分散σc 2は、次式(7)、(8)から求めることができる。なお、μyは、yの平均である。
Let the two-dimensional polynomial be z = f (y) = (1/2) ay 2 + by + c, and the parameters (a, b, c) have already been estimated, and the mean μ a , μ b , variance σ a 2 , σ b. 2 is known. In this case, if the parameters (a, b, c) of the two-dimensional polynomial f (y) are set as θ, f (y) = f (y; θ), and f (y; θ) is differentiated by y, then f '(Y ; θ ) = ay + b. Therefore, the average cave of c 1 and the variance σ c 2 can be obtained from the following equations (7) and (8). Note that μ y is the average of y.
(平均caveと分散σc 2の求め方(収束曲線が高次元多項式の場合))
収束曲線が高次元多項式の場合には、c1の平均caveと、分散σc 2は、以下のようにして求めることができる。( How to find the mean cave and variance σ c 2 (when the convergence curve is a high-dimensional polynomial))
When the convergence curve is a high-dimensional polynomial, the mean cave of c 1 and the variance σ c 2 can be obtained as follows.
高次元多項式の関数f(y)のパラメータをθとし、近傍のデータを用いてθの平均μθ、分散σθ 2を、次式(9)より推定する。
f(y;θ)=Σnθnyn …(9)Let θ be the parameter of the function f (y) of the high-dimensional polynomial, and estimate the mean μ θ and variance σ θ 2 of θ from the following equation (9) using nearby data.
f (y; θ) = Σ n θ n y n … (9)
また、f(y;θ)をyで微分し、f’(y;θ)を算出する。この場合、f’(y;θ)は、yの多項式Σnnθnyn-1と表現される。Further, f (y; θ) is differentiated by y to calculate f'(y; θ). In this case, f'(y; θ) is expressed as the polynomial Σ n n θ n y n-1 of y.
この場合、c1の平均caveは、次式(10)にて表すことができる。
cave=E[f’(μy;θ)]=ΣnnE[θnyn-1] …(10)In this case, the average cave of c 1 can be expressed by the following equation (10).
c ave = E [f'(μ y ; θ)] = Σ n nE [θ n y n-1 ]… (10)
この場合、θとyの独立性を仮定すると、
E[θnyn-1]=E[θn]E[yn-1]=μθnE[yn-1] …(11)
となる。In this case, assuming the independence of θ and y,
E [θ n y n-1 ] = E [θ n ] E [y n-1 ] = μ θ n E [y n-1 ]… (11)
Will be.
ここで、yn-1をμyで2次までTaylor展開することで、平均E[yn-1]をμy、σe 2で近似する(Delta Method)と、次式(12)のようになる。Here, by Taylor-expanding y n-1 to the second order with μ y , the average E [y n-1 ] is approximated by μ y and σ e 2 (Delta Method), and the following equation (12) Will be.
また、E[yn-1]を求めるために、yn-1をμyの周りでTaylor展開し、yn-1の分散(Var(yn-1))をμy、σe 2で近似すると、次式(13)のようになる。
Further, in order to obtain E [y n-1 ], y n-1 is expanded by T y lor around μ y , and the variance of y n-1 (Var (y n-1 )) is μ y , σ. Approximating with e 2 , it becomes as shown in the following equation (13).
上式(12)、(13)を用いることで、任意の確率変数のべき乗表現の平均、分散を求めることができる。 By using the above equations (12) and (13), the average and variance of the power expression of any random variable can be obtained.
また、σc 2についても、yの多項式として表せるので、上記と同様にμy、σe 2を用いて、次式(14)のように表現することができる。なお、θnとθk(n≠k)は独立であるものとする。Further, since σ c 2 can also be expressed as a polynomial of y, it can be expressed as in the following equation (14) by using μ y and σ e 2 in the same manner as described above. It is assumed that θ n and θ k (n ≠ k) are independent.
なお、θnのべき乗表現の平均、分散についても、上式(12)、(13)と同様に算出することができる。The mean and variance of the power expression of θ n can also be calculated in the same manner as in the above equations (12) and (13).
(平均caveと分散σc 2の求め方(収束曲線が多項式でない非線形曲線の場合))
収束曲線が多項式でない非線形曲線の場合にも、Taylor展開することで、近似的に多項式に展開し、上述した手法に適用することができる。例えば、収束曲線の式が、f(y)=Aeayの場合、次式(15)のようにTaylor展開して多項式として扱えばよい。( How to find the mean cave and variance σ c 2 (when the convergence curve is a non-linear curve that is not a polynomial))
Even when the convergence curve is a non-linear curve that is not a polynomial, it can be approximately expanded to a polynomial by Taylor expansion and applied to the above-mentioned method. For example, when the equation of the convergence curve is f (y) = Ae ay , it may be treated as a polynomial by Taylor expansion as in the following equation (15).
以上、詳細に説明したように、本実施形態によると、モデル設定部12は、過去に調整対象となった製品の信号出力が予め定められた許容範囲の境界近傍の値を示したときに該製品に入力されていたパラメータ値と、個体特徴量と、の組み合わせのデータ(履歴DB22のデータ)に基づいて、新たな調整対象の製品50の境界予測値yを個体特徴量から推定するためのモデルを設定し、補正値算出部14は、製品50の信号出力が許容範囲に入る確率を最大にするために境界予測値yを補正する補正値gを、許容範囲内において予測されるパラメータ値に対する信号出力の変化に基づいて算出する。そして、初期値決定部16は、製品50の個体特徴量とモデルとに基づいて境界予測値yを算出(推定)し、算出した境界予測値yを補正値gで補正することで、製品50に入力するパラメータ値の初期値y0を算出する。これにより、過去に調整された製品のデータとして、出力信号が許容範囲の境界近傍の値を示したときのパラメータ値しか存在しない場合であっても、新たな調整対象の製品50に入力するパラメータ値の初期値として、信号出力が許容範囲に入る確率が最大であるパラメータ値を算出することができる。したがって、製品50において信号出力が許容範囲に入るパラメータ値を探索する時間を短縮することができるので、製品50の調整時間を短縮することが可能となる。これにより、製造リードタイムを短縮することができ、生産効率向上を図ることができる。また、調整装置10の数を増やさずに、調整できる製品の数を増加させることができるので、コスト低減を図ることが可能である。As described in detail above, according to the present embodiment, when the signal output of the product to be adjusted in the past shows a value near the boundary of the predetermined allowable range, the
また、本実施形態では、許容範囲内において予測されるパラメータ値に対する信号出力の変化が、線形な変化又は非線形な変化であると仮定し、補正値算出部14は、信号出力の変化に応じた方法で補正値を算出するので、補正値として、許容範囲内の挙動に応じた適切な値を算出することができる。
Further, in the present embodiment, it is assumed that the change in the signal output with respect to the parameter value predicted within the permissible range is a linear change or a non-linear change, and the correction
また、本実施形態では、補正値算出部14は、許容範囲内においてパラメータ値に対する信号出力の変化が非線形な変化を示す場合に、線形近似して補正値を算出することとしている。これにより、簡単な計算により補正値を算出することができる。
Further, in the present embodiment, the correction
また、本実施形態では、補正値算出部14は、境界予測値の誤差(分散)σe 2と、感度の予測誤差(分散)σc 2に基づいて、補正値を算出するので、各誤差を考慮した適切な補正値を算出することができる。Further, in the present embodiment, the correction
なお、上記実施形態では、補正値算出部14は、補正値gとして、製品50の信号出力が許容範囲に入る確率を最大にするため値を算出する場合について説明したが、これに限られるものではない。例えば、補正値算出部14は、確率が所定値以上となる補正値gを算出するようにしてもよい。
In the above embodiment, the correction
なお、モデル設定部12や補正値算出部14は、モデル設定や補正値算出を所定期間ごとに実行(更新)することとしてもよいし、履歴DB22に所定数のデータが新規追加されるごとに実行することとしてもよい。
The
なお、上記実施形態では、調整装置10が、モデル設定部12や補正値算出部14を有する場合について説明したが、これに限らず、モデル設定部12や補正値算出部14は、別装置が有していてもよい。すなわち、モデルの設定や補正値の算出は、別装置が行い、調整装置10は、モデルや補正値を別装置から取得することとしてもよい。別装置は、調整装置10に直接接続された装置であってもよいし、ネットワークを介して接続されたサーバ(クラウドサーバ)であってもよい。
In the above embodiment, the case where the
なお、上記の処理機能は、コンピュータによって実現することができる。その場合、処理装置が有すべき機能の処理内容を記述したプログラムが提供される。そのプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体(ただし、搬送波は除く)に記録しておくことができる。 The above processing function can be realized by a computer. In that case, a program that describes the processing content of the function that the processing device should have is provided. By executing the program on a computer, the above processing function is realized on the computer. The program describing the processing content can be recorded on a computer-readable recording medium (however, the carrier wave is excluded).
プログラムを流通させる場合には、例えば、そのプログラムが記録されたDVD(Digital Versatile Disc)、CD−ROM(Compact Disc Read Only Memory)などの可搬型記録媒体の形態で販売される。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することもできる。 When a program is distributed, it is sold in the form of a portable recording medium such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory) on which the program is recorded. It is also possible to store the program in the storage device of the server computer and transfer the program from the server computer to another computer via the network.
プログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送されるごとに、逐次、受け取ったプログラムに従った処理を実行することもできる。 The computer that executes the program stores, for example, the program recorded on the portable recording medium or the program transferred from the server computer in its own storage device. Then, the computer reads the program from its own storage device and executes the processing according to the program. The computer can also read the program directly from the portable recording medium and execute the processing according to the program. In addition, the computer can sequentially execute processing according to the received program each time the program is transferred from the server computer.
上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。 The embodiments described above are examples of preferred embodiments of the present invention. However, the present invention is not limited to this, and various modifications can be made without departing from the gist of the present invention.
10 調整装置(パラメータ値決定装置)
12 モデル設定部
14 補正値算出部
16 初期値設定部(モデル取得部、補正値取得部、算出部)10 Adjustment device (parameter value determination device)
12
Claims (9)
過去に調整対象となった装置の信号出力が予め定められた許容範囲の境界近傍の値を示したときに該装置に入力されていたパラメータ値と、該装置の個体特徴量と、の組み合わせのデータに基づいて、新たな調整対象の装置の信号出力が前記許容範囲の境界近傍の値を示すようなパラメータ値を該新たな調整対象の装置の個体特徴量から推定するためのモデルを設定し、
前記新たな調整対象の装置の信号出力が前記許容範囲に入る確率を所定以上にするために前記モデルを用いて推定したパラメータ値を補正する補正値を、前記許容範囲内において予測されるパラメータ値に対する信号出力の変化に基づいて算出し、
前記新たな調整対象の装置の個体特徴量と前記モデルとに基づいてパラメータ値を推定し、推定したパラメータ値を前記補正値で補正することで、前記新たな調整対象の装置に入力するパラメータ値を算出する、
処理をコンピュータが実行することを特徴とするパラメータ値決定方法。 It is a parameter value determination method that determines the parameter value to be input to the device to be adjusted that outputs a signal according to the parameter value.
The combination of the parameter value input to the device when the signal output of the device to be adjusted in the past shows a value near the boundary of the predetermined allowable range and the individual characteristic amount of the device. Based on the data, a model is set for estimating the parameter value such that the signal output of the device to be adjusted newly indicates the value near the boundary of the allowable range from the individual characteristic amount of the device to be adjusted. ,
A correction value that corrects a parameter value estimated using the model in order to increase the probability that the signal output of the new device to be adjusted falls within the allowable range is a parameter value predicted within the allowable range. Calculated based on the change in signal output to
A parameter value is estimated based on the individual feature amount of the new adjustment target device and the model, and the estimated parameter value is corrected by the correction value to input the parameter value to the new adjustment target device. To calculate,
A method for determining parameter values, characterized in that the processing is performed by a computer.
前記補正値を算出する処理では、前記信号出力の変化に応じた方法で、前記補正値を算出することを特徴とする請求項1又は2に記載のパラメータ値決定方法。 Assuming that the change in signal output with respect to the parameter value predicted within the permissible range is a linear change or a non-linear change.
The parameter value determination method according to claim 1 or 2, wherein in the process of calculating the correction value, the correction value is calculated by a method according to a change in the signal output.
前記補正値を算出する処理では、前記非線形な変化を線形近似して、前記補正値を算出することを特徴とする請求項3に記載のパラメータ値決定方法。 When the change in the signal output with respect to the parameter value predicted within the permissible range is a non-linear change,
The parameter value determination method according to claim 3, wherein in the process of calculating the correction value, the non-linear change is linearly approximated to calculate the correction value.
過去に調整対象となった装置の信号出力が予め定められた許容範囲の境界近傍の値を示したときに該装置に入力されていたパラメータ値と、該装置の個体特徴量と、の組み合わせのデータに基づいて、新たな調整対象の装置の信号出力が前記許容範囲の境界近傍の値を示すようなパラメータ値を該新たな調整対象の装置の個体特徴量から推定するためのモデルを設定し、
前記新たな調整対象の装置の信号出力が前記許容範囲に入る確率を所定以上にするために前記モデルを用いて推定したパラメータ値を補正する補正値を、前記許容範囲内において予測されるパラメータ値に対する信号出力の変化に基づいて算出し、
前記新たな調整対象の装置の個体特徴量と前記モデルとに基づいてパラメータ値を推定し、推定したパラメータ値を前記補正値で補正することで、前記新たな調整対象の装置に入力するパラメータ値を算出する、
処理をコンピュータに実行させることを特徴とするパラメータ値決定プログラム。 A parameter value determination program that determines the parameter value to be input to the device to be adjusted that outputs a signal according to the parameter value.
The combination of the parameter value input to the device when the signal output of the device to be adjusted in the past shows a value near the boundary of the predetermined allowable range and the individual characteristic amount of the device. Based on the data, a model is set for estimating the parameter value such that the signal output of the device to be adjusted newly indicates the value near the boundary of the allowable range from the individual characteristic amount of the device to be adjusted. ,
A correction value that corrects a parameter value estimated using the model in order to increase the probability that the signal output of the new device to be adjusted falls within the allowable range is a parameter value predicted within the allowable range. Calculated based on the change in signal output to
A parameter value is estimated based on the individual feature amount of the new adjustment target device and the model, and the estimated parameter value is corrected by the correction value to input the parameter value to the new adjustment target device. To calculate,
A parameter value determination program characterized by having a computer execute processing.
過去に調整対象となった装置の信号出力が予め定められた許容範囲の境界近傍の値を示したときに該装置に入力されていたパラメータ値と、該装置の個体特徴量と、の組み合わせのデータに基づいて設定した、新たな調整対象の装置の信号出力が前記許容範囲の境界近傍の値を示すようなパラメータ値を該新たな調整対象の装置の個体特徴量から推定するためのモデルを取得するモデル取得部と、
前記許容範囲内において予測されるパラメータ値に対する信号出力の変化に基づいて算出した、前記新たな調整対象の装置の信号出力が前記許容範囲に入る確率を所定以上にするために前記モデルを用いて推定したパラメータ値を補正する補正値を取得する補正値取得部と、
前記新たな調整対象の装置の個体特徴量と前記モデルとに基づいてパラメータ値を推定し、推定したパラメータ値を前記補正値で補正することで、前記新たな調整対象の装置に入力するパラメータ値を算出する算出部と、
を備えるパラメータ値決定装置。 It is a parameter value determining device that determines the parameter value to be input to the device to be adjusted that outputs a signal according to the parameter value.
The combination of the parameter value input to the device when the signal output of the device to be adjusted in the past shows a value near the boundary of the predetermined allowable range and the individual characteristic amount of the device. A model for estimating a parameter value set based on the data such that the signal output of the new device to be adjusted indicates a value near the boundary of the permissible range from the individual characteristic amount of the device to be adjusted. Model acquisition department to acquire and
Using the model to make the probability that the signal output of the new device to be adjusted falls within the permissible range, which is calculated based on the change in the signal output with respect to the parameter value predicted within the permissible range. A correction value acquisition unit that acquires a correction value that corrects the estimated parameter value,
A parameter value is estimated based on the individual feature amount of the new adjustment target device and the model, and the estimated parameter value is corrected by the correction value to input the parameter value to the new adjustment target device. And the calculation unit that calculates
A parameter value determining device comprising.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2017/001486 WO2018134903A1 (en) | 2017-01-18 | 2017-01-18 | Parameter value determination method, parameter value determination program, and parameter value determination device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2018134903A1 JPWO2018134903A1 (en) | 2019-11-07 |
JP6863393B2 true JP6863393B2 (en) | 2021-04-21 |
Family
ID=62907918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018562775A Active JP6863393B2 (en) | 2017-01-18 | 2017-01-18 | Parameter value determination method, parameter value determination program and parameter value determination device |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP6863393B2 (en) |
WO (1) | WO2018134903A1 (en) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06339295A (en) * | 1993-05-26 | 1994-12-06 | Toshiba Corp | Regulator |
JPH11296561A (en) * | 1998-04-07 | 1999-10-29 | Toshiba Corp | Method and device for generating worst case model parameter |
JP2004334714A (en) * | 2003-05-09 | 2004-11-25 | Yamaha Motor Co Ltd | Parameter optimization method, parameter optimization device, parameter optimization program, and sailing control device |
-
2017
- 2017-01-18 WO PCT/JP2017/001486 patent/WO2018134903A1/en active Application Filing
- 2017-01-18 JP JP2018562775A patent/JP6863393B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
WO2018134903A1 (en) | 2018-07-26 |
JPWO2018134903A1 (en) | 2019-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110686352B (en) | Temperature detection value compensation method and device, storage medium and air conditioner | |
JPWO2018047655A1 (en) | Time-series data feature extraction device, time-series data feature extraction method, and time-series data feature extraction program | |
CN112468097B (en) | Temperature compensation method, radio frequency device and storage medium | |
WO2016147722A1 (en) | Estimating device, estimating method and program | |
CN112445357B (en) | Calibration method, calibration device and storage medium | |
US20180337837A1 (en) | Data center management method, management apparatus, and data center system | |
JP6863393B2 (en) | Parameter value determination method, parameter value determination program and parameter value determination device | |
CN112067166B (en) | Thermal resistor checking method | |
WO2023199588A1 (en) | Control device, control method, and program | |
JP2019039764A (en) | Impedance estimating device | |
CN109814085B (en) | Laser radar temperature compensation method, device, system and equipment | |
JP6862841B2 (en) | Processing device, adjustment parameter prediction model estimation method, and adjustment parameter prediction model estimation program | |
CN114296496A (en) | Temperature control method and device | |
JP6379709B2 (en) | Signal processing apparatus, signal processing method, and program | |
CN112152232B (en) | Load frequency control method and device for power system | |
US20210194071A1 (en) | Secondary battery management device, secondary battery management method, and non-transitory computer readable storage medium | |
JP6204318B2 (en) | Similarity Evaluation Device, Similarity Evaluation System, Similarity Evaluation Device Method, and Similarity Evaluation Program | |
JP2019008470A (en) | Management device, simulation system, and simulation method | |
US11726470B2 (en) | Method and device for determining a measure of quality continuously | |
CN112989632B (en) | Solar energy resource short-term forecasting method based on satellite radiation product | |
US11703411B2 (en) | Signal processing device and signal processing method | |
CN113485100B (en) | Robot leg length planning method and device, readable storage medium and robot | |
JP2004178189A (en) | Estimation method of quantization error, identification method of plant, control method, estimation device of quantization error, and program | |
CN113536213B (en) | Method and processor for determining seed vigor | |
CN118362037A (en) | Sensor calibration method, device, medium, terminal equipment and program product |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190527 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190527 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200714 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200828 |
|
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: 20210302 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210315 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6863393 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |