JP2019185591A - Experiment assisting device and experiment assisting method - Google Patents

Experiment assisting device and experiment assisting method Download PDF

Info

Publication number
JP2019185591A
JP2019185591A JP2018078363A JP2018078363A JP2019185591A JP 2019185591 A JP2019185591 A JP 2019185591A JP 2018078363 A JP2018078363 A JP 2018078363A JP 2018078363 A JP2018078363 A JP 2018078363A JP 2019185591 A JP2019185591 A JP 2019185591A
Authority
JP
Japan
Prior art keywords
experiment
optimization
design
plan
optimization method
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.)
Pending
Application number
JP2018078363A
Other languages
Japanese (ja)
Inventor
達也 長谷部
Tatsuya Hasebe
達也 長谷部
一朗 片岡
Ichiro Kataoka
一朗 片岡
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2018078363A priority Critical patent/JP2019185591A/en
Publication of JP2019185591A publication Critical patent/JP2019185591A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

To efficiently determine an appropriate experiment plan.SOLUTION: An experiment assisting device of the present invention is characterized by comprising: an input processing unit for accepting track record data including a product design value and a performance value, an optimization goal, and constraints; a feature quantity calculation unit for calculating feature quantity of the accepted track record data; an experiment plan/optimization technique selection unit for selecting an experiment plan or an optimization technique in accordance with the calculated feature quantity; an experiment plan/optimization execution unit for executing the selected experiment plan or the optimization technique; an experiment plan/optimization technique effect calculation unit for calculating an effect of the experiment plan or the optimization technique on the basis of an execution result of the experiment plan or the optimization technique; and an output processing unit for outputting the effect calculated in association with the experiment plan or the optimization technique.SELECTED DRAWING: Figure 10

Description

本発明は、実験支援装置及び実験支援方法に関する。   The present invention relates to an experiment support apparatus and an experiment support method.

製品の設計条件は、その設計条件によって製品が所望の性能やコストを満たすように決定される。そのためには、様々な設計条件に応じた試作品による実験、あるいは設計条件に対応する解析条件でのシミュレーションが多数繰り返される。実験とシミュレーションの中にはコスト及び時間を要するものもあり、無駄な試行錯誤を最低限に抑えることが重要である。そこで、設計条件と性能やコストの関係や所望の性能・コストを満たすような設計条件を最小の実験またはシミュレーション回数で明らかにするような実験・シミュレーション条件の計画を行う必要がある。これを実験計画と呼ぶ。設計の初期段階では、まず設計条件と性能やコストの関係を把握するための実験計画を行う。十分関係が把握できれば、性能目標・コスト目標に近づくための実験計画を行うか、把握できている傾向に基づいて最も目標に近い設計案を得ることができる。これら、性能目標やコスト目標に最も近い設計条件を得ること、あるいは目標に最も近い設計条件得るための実験計画を行うことを最適化と呼ぶ。これら実験計画及び最適化をあわせて実験計画・最適化と呼称する。また、以後、実験およびシミュレーションをまとめて実験と呼称する。   The design conditions of the product are determined according to the design conditions so that the product satisfies desired performance and cost. For this purpose, many experiments with prototypes corresponding to various design conditions or simulations under analysis conditions corresponding to the design conditions are repeated. Some experiments and simulations are costly and time consuming, and it is important to minimize wasteful trial and error. Therefore, it is necessary to plan the experiment / simulation conditions so that the relationship between the design conditions and the performance / cost and the design conditions satisfying the desired performance / cost are clarified with the minimum number of experiments or simulations. This is called an experimental design. In the initial stage of design, an experiment plan is first designed to grasp the relationship between design conditions, performance, and cost. If a sufficient relationship can be grasped, an experiment plan for approaching the performance goal / cost goal can be performed, or a design plan closest to the goal can be obtained based on the grasped tendency. Obtaining the design conditions closest to the performance goal and cost goal, or performing an experimental design to obtain the design conditions closest to the goal is called optimization. These experimental design and optimization are collectively referred to as experimental design / optimization. Hereinafter, experiments and simulations are collectively referred to as experiments.

近年多くの実験計画・最適化の手法が活用されている。遺伝的アルゴリズムに代表される大域的最適化を行うヒューリスティックな手法、設計条件と性能との関係を示す応答曲面上で最適な設計条件を決定するベイズ最適化に代表される統計的手法、直交表等を使用した古典的手法等がその例である。これらの手法のそれぞれは、それによって得られる実験計画(最適解)の特徴及び数が異なる。また一般的に、製品の設計は、まず設計条件と性能との関連性の把握、次に得られた関連性に基づく実験による試行錯誤、さらに所望の性能を出力する設計条件の最適化という段階を順に辿る。前記した各手法は、どの段階においても好ましく適用され得るわけではない。各段階において、手法を取捨選択することが必要である。   In recent years, many experimental design and optimization techniques have been used. A heuristic method for global optimization represented by genetic algorithms, a statistical method represented by Bayesian optimization for determining optimal design conditions on the response surface indicating the relationship between design conditions and performance, orthogonal table An example is a classical method using the above. Each of these methods differs in the characteristics and number of experimental designs (optimal solutions) obtained thereby. In general, product design involves the steps of grasping the relationship between design conditions and performance, then trial and error through experiments based on the obtained relationships, and further optimizing the design conditions to output the desired performance. Trace in order. Each method described above cannot be preferably applied at any stage. At each stage, it is necessary to select a method.

特許文献1の実験計画支援方法は、ユーザが実験数、設計項目(因子)の水準数等を入力するのを受け付け、それに適した実験計画を出力する。特許文献2の割り付け方法は、設計項目の水準数及び設計項目間の交互作用をユーザが入力するのを受け付け、それに適した直交表のタイプを出力する。特許文献3のテスト計画作成方法は、テストケースに優先順位を設定したうえでテストケースの工数を予測し、優先順位にしたがって、工数に必要な人的リソースを割り当てる。さらに、当該テスト計画作成方法は、テストケースが所定の期限までに終了しない場合、人的リソースの割り付けを見直す。   The experiment planning support method of Patent Document 1 accepts that the user inputs the number of experiments, the number of levels of design items (factors), and the like, and outputs an experiment plan suitable for it. The allocation method of Patent Document 2 accepts the user's input of the number of design item levels and the interaction between design items, and outputs an orthogonal table type suitable for the input. The test plan creation method of Patent Document 3 predicts the man-hours of a test case after setting priorities to the test cases, and allocates human resources necessary for the man-hours according to the priorities. Furthermore, the test plan creation method reviews the allocation of human resources when the test case does not end by a predetermined time limit.

米国特許第6448983号明細書US Pat. No. 6,448,983 特開2001−265830号公報JP 2001-265830 A 特開2011−257947号公報JP 2011-257947 A

製品の設計者は、取り組んでいる問題に適した実験計画・最適化の手法を、その問題の段階に応じて、多くの候補のなかから的確に選択しなければならない。特に経験の少ない設計者にとって、このような選択は難しい。また、ノウハウが豊富な熟練設計者であっても、逆にそのノウハウが災いし客観的に適切な実験計画・最適化の手法に辿りつかない場合もある。特許文献1の実験計画支援方法及び特許文献2の割り付け方法が出力(提案)する実験計画は、直交表又は多元配置表を使用することを前提としている。これらは、実験初期において設計条件と性能との傾向を知る場合には有効である。しかしながら、既にある程度の実験の実績がある場合、別の手法(例えば統計的手法)が必要となることが多い。特許文献3のテスト計画作成方法は、ソフトウエア開発のテストケースに焦点を絞っており一般的な機械等の設計において実験計画・最適化の手法を決定するには、別途工夫が必要である。
そこで、本発明は、適切な実験計画を効率的に決定することを目的とする。
Product designers must select an appropriate experimental design / optimization method for the problem they are addressing from among many candidates, depending on the stage of the problem. Such a choice is difficult, especially for designers with little experience. In addition, even a skilled designer with abundant know-how may have a problem that the know-how is damaged and the objective experimental design / optimization method cannot be reached. The experimental design output (suggested) by the experimental design support method of Patent Literature 1 and the allocation method of Patent Literature 2 is premised on using an orthogonal table or a multi-way allocation table. These are effective when knowing the trend between design conditions and performance at the beginning of the experiment. However, if there is already some experience in experiments, another method (for example, a statistical method) is often required. The test plan creation method of Patent Document 3 focuses on a test case for software development, and a separate device is required to determine an experiment plan / optimization method in designing a general machine or the like.
Therefore, an object of the present invention is to efficiently determine an appropriate experiment plan.

本発明の実験支援装置は、製品の設計値及び性能値を含む実績データ、最適化目標、及び、制約を受け付ける入力処理部と、受け付けた実績データの特徴量を算出する特徴量算出部と、算出した特徴量に応じて実験計画又は最適化手法を選択する実験計画・最適化手法選択部と、選択した実験計画又は最適化手法を実行する実験計画・最適化実行部と、実験計画又は最適化手法の実行結果に基づき実験計画又は最適化手法の効果を算出する実験計画・最適化効果算出部と、実験計画又は最適化手法に関連付けて算出した効果を出力する出力処理部と、を備えることを特徴とする。
その他の手段については、発明を実施するための形態のなかで説明する。
The experiment support apparatus of the present invention includes an input processing unit that receives actual data including design values and performance values of products, an optimization target, and a constraint, and a characteristic amount calculation unit that calculates a characteristic amount of the received actual data, An experiment plan / optimization method selection unit that selects an experiment plan or optimization method according to the calculated feature value, an experiment plan / optimization execution unit that executes the selected experiment plan or optimization method, and an experiment plan or optimization An experiment plan / optimization effect calculation unit that calculates the effect of the experiment plan or the optimization method based on the execution result of the optimization method, and an output processing unit that outputs the effect calculated in association with the experiment plan or the optimization method It is characterized by that.
Other means will be described in the embodiment for carrying out the invention.

本発明によれば、適切な実験計画を効率的に決定することができる。   According to the present invention, an appropriate experimental design can be determined efficiently.

実験支援装置の構成を説明する図である。It is a figure explaining the structure of an experiment assistance apparatus. 製品を説明する図である。It is a figure explaining a product. 設計値管理情報を説明する図である。It is a figure explaining design value management information. 特徴量情報を説明する図である。It is a figure explaining feature-value information. 全体処理手順のフローチャートである。It is a flowchart of the whole processing procedure. 実験計画・最適化決定処理手順のフローチャートである。It is a flowchart of an experiment plan and optimization determination processing procedure. 削除試行処理手順のフローチャートである。It is a flowchart of a deletion trial processing procedure. 追加試行処理手順のフローチャートである。It is a flowchart of an additional trial process procedure. 入力画面の一例を示す図である。It is a figure which shows an example of an input screen. 出力画面の一例を示す図である。It is a figure which shows an example of an output screen.

以降、本発明を実施するための形態(“本実施形態”という)を、図等を参照しながら詳細に説明する。本実施形態は、ある程度蓄積された実験の実績データを有する設計者(本実施形態の実験支援装置のユーザ)が、次の段階の実験計画を作成する例である。   Hereinafter, a mode for carrying out the present invention (referred to as “the present embodiment”) will be described in detail with reference to the drawings. This embodiment is an example in which a designer (experiment support apparatus user of the present embodiment) having experiment result data accumulated to some extent creates an experiment plan for the next stage.

(実験支援装置の構成)
図1に沿って、実験支援装置1の構成を説明する。実験支援装置1は、一般的なコンピュータであり、中央制御装置11、キーボード、タッチパネル等の入力装置12、ディスプレイ等の出力装置13、主記憶装置14及び補助記憶装置15を備える。これらは、バスで相互に接続されている。補助記憶装置15は、設計値管理情報31、特徴量情報32及び実験計画・最適化手法群33(いずれも詳細後記)を格納している。
(Configuration of experiment support device)
The configuration of the experiment support apparatus 1 will be described with reference to FIG. The experiment support apparatus 1 is a general computer, and includes a central control device 11, an input device 12 such as a keyboard and a touch panel, an output device 13 such as a display, a main storage device 14, and an auxiliary storage device 15. These are connected to each other by a bus. The auxiliary storage device 15 stores design value management information 31, feature amount information 32, and an experiment plan / optimization method group 33 (all of which will be described later in detail).

主記憶装置14における、入力処理部21、特徴量算出部22、実験計画・最適化手法選択部23、実験計画・最適化実行部24、実験計画・最適化効果算出部25及び出力処理部26は、プログラムである。以降において、“○○部は”と動作主体を記した場合、それは、中央制御部11が補助記憶装置15から“○○部”を読み出して主記憶装置14にロードしたうえで、後記する処理を実行することを意味する。   In the main storage device 14, an input processing unit 21, a feature amount calculation unit 22, an experiment plan / optimization method selection unit 23, an experiment plan / optimization execution unit 24, an experiment plan / optimization effect calculation unit 25, and an output processing unit 26. Is a program. In the following, when the operation subject is described as “XX section”, it means that the central control section 11 reads “XX section” from the auxiliary storage device 15 and loads it into the main storage device 14, and then the processing described later. Means to execute.

(製品)
図2に沿って、製品を説明する。本実施形態における製品は、例えば、回転機械41である。回転機械41は、回転軸42、軸受43及び軸受44を有する。電動機(図示せず)は、回転軸42に対して回転力を与える。いま、ユーザは、回転軸42の振動、電動機の消費電力、軸受43及び軸受43における回転軸42の温度等を最低限に抑えるような設計条件を決定しようとしている。
(Product)
The product will be described with reference to FIG. The product in the present embodiment is, for example, a rotary machine 41. The rotating machine 41 includes a rotating shaft 42, a bearing 43, and a bearing 44. The electric motor (not shown) gives a rotational force to the rotary shaft 42. Now, the user is trying to determine design conditions that minimize the vibration of the rotating shaft 42, the power consumption of the motor, the bearing 43, the temperature of the rotating shaft 42 in the bearing 43, and the like.

(実験、シミュレーション等)
実験とは、製品の試作品を運転環境に置き、設計値(直ちに後記)で試作品を運転し、試作品の性能値(直ちに後記)を様々な種類のセンサ等で取得することである。シミュレーションとは、コンピュータ上でモデルに対して設計値を入力し、モデルに性能値を出力させることである。モデルとは、複数の変数を含む数式の集合である。実験計画・最適化手法とは、所望の性能を出力するような設計条件を、モデル上あるいは応答曲面(詳細後記)上で探索する最適化手法、及び、所望の性能を出力する設計条件を得るための実験を計画する実験計画手法のことである。実験計画・最適化手法を実現する異なる複数のソフトウエアアプリケーションが、図1における実験計画・最適化手法群33である。本実施形態の特徴は、適切な実験計画・最適化手法を選択することにある。
(Experiment, simulation, etc.)
An experiment is to place a prototype of a product in an operating environment, operate the prototype with design values (immediately described later), and obtain performance values (immediately later) of the prototype with various types of sensors. Simulation refers to inputting design values for a model on a computer and outputting performance values to the model. A model is a set of mathematical expressions including a plurality of variables. The experimental design / optimization method refers to an optimization method for searching for a design condition that outputs desired performance on a model or a response surface (detailed later), and a design condition for outputting desired performance. It is an experimental design method for planning an experiment. A plurality of different software applications for realizing the experiment design / optimization technique is the experiment design / optimization technique group 33 in FIG. The feature of this embodiment is that an appropriate experiment design / optimization method is selected.

(設計値)
回転軸42の種類及び用途によって、様々な設計値が想定され得る。ここでは、以下の項目が設計値であるとする。
・回転軸42の軸方向の長さ“a”
・軸受43の位置、すなわち、回転軸42の左端からの距離“b”
・軸受44の位置、すなわち、回転軸42の左端からの距離“c”
・回転軸42の回転速度
・回転軸42の任意の箇所で軸の中心に向かって掛けられる力。この力は、想定される負荷の値(タービン羽の重力等)であってもよいし、実験用の加振力であってもよい。
(Design value)
Various design values can be assumed depending on the type and application of the rotating shaft 42. Here, it is assumed that the following items are design values.
・ Axial length “a” of the rotating shaft 42
The position of the bearing 43, that is, the distance “b” from the left end of the rotating shaft 42
The position of the bearing 44, that is, the distance “c” from the left end of the rotating shaft 42
The rotational speed of the rotating shaft 42 The force applied to the center of the rotating shaft 42 at an arbitrary position. This force may be an assumed load value (such as the gravity of a turbine blade) or an experimental excitation force.

(性能値)
回転軸42の種類及び用途によって、様々な性能値が想定され得る。ここでは、以下の項目が性能値であるとする。
・回転軸42の一次モード振動周波数
・回転軸42の二次モード振動周波数
・回転軸42に回転力を与える電動機の消費電力
・軸受43の温度(軸受43における回転軸42の温度に等しい。)
・軸受44の温度(軸受44における回転軸42の温度に等しい。)
(Performance value)
Various performance values can be assumed depending on the type and application of the rotating shaft 42. Here, it is assumed that the following items are performance values.
The primary mode vibration frequency of the rotating shaft 42, the secondary mode vibration frequency of the rotating shaft 42, the electric power consumption of the motor that applies a rotational force to the rotating shaft 42, and the temperature of the bearing 43 (equal to the temperature of the rotating shaft 42 in the bearing 43).
The temperature of the bearing 44 (equal to the temperature of the rotating shaft 42 in the bearing 44)

前記よりわかるように、本実施形態の設計値は5次元のベクトルであり、性能値もまた5次元のベクトルである。両者の次数が一致しているのは偶然である。設計値が因果関係の“原因”であり、性能値が因果関係の“結果”であると解釈することも可能である。つまり、設計値Xをある関数Fの入力(説明変数)と解釈し、性能値Yを当該関数Fの出力(目的変数)と解釈することも可能である。F、X及びYの関係を以下に数式で表す。   As can be seen from the above, the design value of this embodiment is a five-dimensional vector, and the performance value is also a five-dimensional vector. It is a coincidence that the orders of both coincide. It is also possible to interpret that the design value is the “cause” of the causal relationship and the performance value is the “result” of the causal relationship. That is, the design value X can be interpreted as an input (explanatory variable) of a certain function F, and the performance value Y can be interpreted as an output (object variable) of the function F. The relationship among F, X, and Y is expressed by mathematical formulas below.

Y=F(X)
X=(x,x,x,x,x
Y=(y,y,y,y,y
Y = F (X)
X = (x 1 , x 2 , x 3 , x 4 , x 5 )
Y = (y 1 , y 2 , y 3 , y 4 , y 5 )

ここで、x、x、・・・、xは、前記した5種類の設計値である。y、y、・・・、yは、前記した5種類の性能値である。関数Fが未知であっても、設計値X及び性能値Yの組合せ(実験の実績データ)が充分多く存在していれば、実績データを学習データとして、関数Fを推定することが可能である。 Here, x 1 , x 2 ,..., X 5 are the above-described five types of design values. y 1 , y 2 ,..., y 5 are the five types of performance values described above. Even if the function F is unknown, if there are a sufficient number of combinations of the design value X and the performance value Y (experimental result data), the function F can be estimated using the actual data as learning data. .

(応答曲面)
実験やシミュレーションによって得られた設計値X及び性能値Yの離散的なデータを連続的な曲面で近似するとき、この曲面を応答曲面と呼ぶ。この応答曲面を作成するのに用いたデータXおよびYを学習データと呼ぶ。応答曲面の例として、関数Fを以下のように1次関数Gで近似する例を挙げる。
Y=F(X)=G(X)+誤差=a+a+a+・・・+a+誤差
ここで、a、a、a、・・・は学習によって推定されるパラメータであり、以下の関数Fとの2乗誤差が最小になるように決定される。
(F(X)―G(X))
このa、a、a、・・・の決定方法は最小2乗法と呼ばれる。
(Response surface)
When discrete data of the design value X and the performance value Y obtained by experiments and simulations are approximated by a continuous curved surface, this curved surface is called a response curved surface. Data X and Y used to create this response surface are called learning data. As an example of the response surface, an example in which the function F is approximated by a linear function G as follows is given.
Y = F (X) = G (X) + error = a 0 + a 1 x 1 + a 2 x 2 +... + A 5 x 5 + error where a 1 , a 2 , a 3 ,. It is a parameter estimated by learning, and is determined so that the square error with the following function F is minimized.
(F (X) -G (X)) 2
This determination method of a 1 , a 2 , a 3 ,... Is called a least square method.

(応答曲面のパラメータ)
学習データの他に必要な応答曲面の位置及び形状を決定するための値を“ハイパーパラメータ”と呼ぶ。ハイパーパラメータは、学習によって決められないパラメータであり、学習によって推定されるパラメータa、a、a、・・・とは区別される。
ハイパーパラメータの例として、最小2乗法を発展させたRidge回帰を挙げる。Ridge回帰では、a、a、a、・・・が過度に大きくならないように、以下の値が最小になるようにa、a、a、・・・を決定する。
(F(X)―G(X))+λ(a +a +・・・+a
ここで、新しく導入されたλはハイパーパラメータであり、λの大きさを変えることでa、a、a・・・の大きさをコントロールすることができ、過学習を防ぐことができる。λは学習によって決めることができないため、ハイパーパラメータである。学習されたパラメータ及びハイパーパラメータを含めて応答曲面のパラメータと呼ぶ。
(Response surface parameters)
In addition to the learning data, a value for determining the position and shape of the required response surface is called a “hyper parameter”. Hyper parameters are parameters that are not determined by learning, and are distinguished from parameters a 1 , a 2 , a 3 ,... Estimated by learning.
As an example of a hyperparameter, Ridge regression that is an extension of the least square method is given. In Ridge regression, a 1 , a 2 , a 3 ,... Are determined so that the following values are minimized so that a 1 , a 2 , a 3 ,.
(F (X) −G (X)) 2 + λ (a 0 2 + a 1 2 +... + A 5 2 )
Here, newly introduced λ is a hyper parameter, and by changing the size of λ, the sizes of a 1 , a 2 , a 3 ... Can be controlled, and overlearning can be prevented. . Since λ cannot be determined by learning, it is a hyperparameter. The learned parameters and hyperparameters are called response surface parameters.

本実施形態の実験計画・最適化手法には、以下2つの類型が存在する。   The experiment design / optimization method of the present embodiment has the following two types.

〈類型1:設計値の最適解を求める実験計画・最適化手法〉
類型1の実験計画・最適化手法は、最適化目標として選択された1又は複数の性能値yが所望の値を取るように、設計値Xを最適化する。したがって、ユーザは、最適化された設計値Xを使用して、試作品による実験を行い、実際に試作品が所望の性能値を出力することを確認することになる。類型1の具体例は、以下の通りである。
<Type 1: Experiment design / optimization method for finding optimal solution of design value>
The type 1 experiment design / optimization method optimizes the design value X so that one or a plurality of performance values y i selected as the optimization target take a desired value. Therefore, the user performs an experiment with a prototype using the optimized design value X, and confirms that the prototype actually outputs a desired performance value. Specific examples of Type 1 are as follows.

・導関数を使用する“勾配法”
・シンプレックス法等の“線形計画法”
・成績のよい2つの実験点を掛けあわせる“遺伝的アルゴリズム”
・ "Gradient method" using derivatives
・ "Linear programming" such as simplex method
・ "Genetic algorithm" that multiplies two experimental points with good results

〈類型2:設計値の最適解の候補を複数求める実験計画・最適化手法〉
類型2の実験計画・最適化手法は、最適解を1つに決定しないまでも、最適解となりそうな設計値Xの候補を複数決定する。したがって、ユーザは、複数の設計値Xの候補を試行錯誤的に使用して試作品による実験を行い、試作品が所望の性能値を出力するような設計値Xを自ら決定することになる。ユーザに試行錯誤を課する点で、類型2は、類型1に比してやや迂遠である。しかしながら、実験回数を削減できるということが、類型2の存在意義である。類型2の具体例は、以下の通りである。
<Type 2: Experiment design / optimization method for finding multiple candidates for the optimal solution of the design value>
The type 2 experiment design / optimization method determines a plurality of design value X candidates that are likely to be optimal solutions, even if the optimal solution is not determined as one. Therefore, the user uses a plurality of design value X candidates in a trial and error manner, performs an experiment with a prototype, and determines the design value X such that the prototype outputs a desired performance value. Type 2 is slightly more detoured than Type 1 in that it imposes trial and error on the user. However, the existence of Type 2 is that the number of experiments can be reduced. A specific example of type 2 is as follows.

・与えられた水準数、因子数(設計値の種類の数)、因子間の交互作用等に基づいて、直交表を作成する“直交表”
・与えられた水準数、因子数等に基づいて、設計値の総当たりの組合せを作成する“多元配置表”
・D−最適計画法、I−最適計画法等の、既存実験点を活用する “最適計画法”
・既存実験点の密度が低い部分に、新たなサンプルとして実験点を追加する“密度による計画法”
・クリギング、ブートストラップサンプリング及び機械学習を使用して応答曲面の予測値及び分散を算出し、EI(Expected Improvement:期待改善量)等の獲得関数を使用する“ベイズ最適化”
・ "Orthogonal table" that creates an orthogonal table based on the given number of levels, number of factors (number of types of design values), interaction between factors, etc.
-“Multiple layout table” that creates brute force combinations of design values based on the number of levels, factors, etc.
・ "Optimum planning" using existing experimental points such as D-optimal programming and I-optimal programming
・ "Design by density" that adds experimental points as a new sample to the part where the density of existing experimental points is low
“Bayesian optimization” that uses kriging, bootstrap sampling and machine learning to calculate the predicted value and variance of the response surface and uses an acquisition function such as EI (Expected Improvement)

例えば、設計値の種類が前記のように5つ存在し、そのそれぞれについて3つの水準が存在する場合、理論的には、3=243回の実験が必要になる。しかしながら直交表を作成する実験計画・最適化手法は、この243回を、例えば27回に削減する。 For example, if there are five types of design values as described above and there are three levels for each of them, theoretically, 3 5 = 243 experiments are required. However, the experimental design / optimization method for creating the orthogonal table reduces this 243 times to 27 times, for example.

(実験計画・最適化手法の効果)
実験支援装置1は、実験計画・最適化手法を実行した結果に対する効果を算出する。効果として使用される具体的な指標の例は、以下に定義する、性能改善度A(式(1))、予測精度改善度I(式(2))、スクリーニング達成度D(式(3))、信頼度C(式(4))及び必要実験回数の5つである。
(Effects of experimental design and optimization methods)
The experiment support apparatus 1 calculates the effect on the result of executing the experiment design / optimization technique. Examples of specific indices used as effects include the performance improvement degree A (formula (1)), the prediction accuracy improvement degree I (formula (2)), and the screening achievement degree D (formula (3)) defined below. ), Reliability C (formula (4)) and required number of experiments.

性能改善度Aは、式(1)によって定義される。ここで、“x”は、実験計画・最適化手法で提案される設計値である。白抜き太字のEは、確率変数の期待値である。“〜”は、確率変数が確率分布に従うことを示す。性能改善度Aは、追加する実験点におけるyの期待改善量(EI値)である。 The performance improvement degree A is defined by equation (1). Here, “x c ” is a design value proposed by the experiment design / optimization method. E in open bold is the expected value of the random variable. “˜” indicates that the random variable follows a probability distribution. The performance improvement degree A is an expected improvement amount (EI value) of y at the experimental point to be added.

Figure 2019185591
Figure 2019185591

予測精度改善度Iは、式(2)によって定義される。ここで、“L1”は、L1ノルムである。“μ”及び“σ”は、それぞれ線形回帰、クリギング等によって得られる性能値yの予想値の平均及び分散である。“μ’”及び“σ’”は、応答曲面を学習する際の学習データに(x,μ(x))を追加した場合のyの予測値の平均及び分散である。予測精度改善度Iが減少することは、探索領域(xの定義域)全体で性能値yの予測精度が改善されていることを示す。 The prediction accuracy improvement degree I is defined by equation (2). Here, “L1” is the L1 norm. “Μ” and “σ” are the mean and variance of the expected values of the performance value y obtained by linear regression, kriging or the like, respectively. “Μ ′” and “σ ′” are the mean and variance of the predicted value of y when (x c , μ (x c )) is added to the learning data when learning the response surface. The decrease in the prediction accuracy improvement degree I indicates that the prediction accuracy of the performance value y is improved in the entire search region (x definition region).

Figure 2019185591
Figure 2019185591

スクリーニング達成度Dは、式(3)によって定義される。ここで、“L∞”は、L∞ノルムである。スクリーニング達成度Dが減少することは、線形回帰等を行った場合の応答曲面のパラメータがよく予測できていることを示す。   Screening achievement degree D is defined by equation (3). Here, “L∞” is an L∞ norm. The decrease in the screening achievement degree D indicates that the parameters of the response surface can be well predicted when linear regression or the like is performed.

Figure 2019185591
Figure 2019185591

信頼度Cは、式(4)によって定義される。ここで、“p(x)”は、実験点密度である。“α”は、定数である。   The reliability C is defined by equation (4). Here, “p (x)” is the experimental point density. “Α” is a constant.

Figure 2019185591
Figure 2019185591

必要実験回数とは、例えば直交表を用いて実験をする場合の、直交表の行数である。   The required number of experiments is, for example, the number of rows in the orthogonal table when an experiment is performed using the orthogonal table.

(設計値管理情報)
図3に沿って、設計値管理情報31を説明する。設計値管理情報31においては、実験ID欄101に記憶された実験IDに関連付けて、性能値欄102には性能値が、設計値欄103には設計値が、取得時点欄104には取得時点が、実績/計画欄105には実績/計画フラグが、実験計画・最適化手法欄106には実験計画・最適化手法が、効果欄107には効果が、提案欄108には提案フラグが記憶されている。
(Design value management information)
The design value management information 31 will be described with reference to FIG. In the design value management information 31, in association with the experiment ID stored in the experiment ID column 101, the performance value column 102 has a performance value, the design value column 103 has a design value, and the acquisition time column 104 has an acquisition time point. However, the actual / plan column 105 stores the actual / plan flag, the experiment plan / optimization method column 106 stores the experiment plan / optimization method, the effect column 107 stores the effect, and the proposal column 108 stores the proposal flag. Has been.

実験ID欄101の実験IDは、実績データ、又は、実験計画・最適化手法による実験計画案若しくは最適化案を一意に特定する識別子である。
性能値欄102の性能値は、前記した性能値である。ここでは、前記した5種類の性能値が小欄102a〜102eに記憶されている。
設計値欄103の設計値は、前記した設計値である。ここでは、前記した5種類の設計値が小欄103a〜103eに記憶されている。
取得時点欄104の取得時点は、実験によって性能値が実際に取得された時点の年月日である。
The experiment ID in the experiment ID column 101 is an identifier that uniquely identifies the result data or the experiment plan or optimization plan based on the experiment plan / optimization method.
The performance value in the performance value column 102 is the above-described performance value. Here, the five types of performance values described above are stored in the small columns 102a to 102e.
The design value in the design value column 103 is the above-described design value. Here, the above-described five types of design values are stored in the small columns 103a to 103e.
The acquisition time point in the acquisition time column 104 is the date when the performance value is actually acquired through an experiment.

実績/計画欄105の実績/計画フラグは、“実績”又は“計画”のうちのいずれかである。“実績”は、その設計値を使用した実験が行われ、実験結果としての性能値が既に取得されていることを示す。“計画”は、その設計値を使用した実験が未だ行われていないことを示す。
実験計画・最適化手法欄106の実験計画・最適化手法は、実験計画・最適化手法の名称である。
効果欄107の効果は、前記した実験計画・最適化手法の効果を示す指標値(A、I、D、C、又は、必要実験回数)に対して所定の閾値を適用した結果であり、ここでは、“大”又は“小”のいずれかである。前記した5種類の効果のうちには、その指標値が大きくなるほど、ここでの効果が“大”となるものもあり、その指標値が小さくなるほど、ここでの効果が“大”となるものもある。
The actual / plan flag in the actual / plan column 105 is either “actual” or “plan”. “Actual” indicates that an experiment using the design value is performed and a performance value as an experimental result has already been acquired. “Plan” indicates that an experiment using the design value has not yet been performed.
The experiment plan / optimization technique in the experiment plan / optimization technique column 106 is the name of the experiment plan / optimization technique.
The effect in the effect column 107 is a result of applying a predetermined threshold to an index value (A, I, D, C, or necessary number of experiments) indicating the effect of the above-described experiment design / optimization method. Then, it is either “large” or “small”. Among the above-mentioned five types of effects, there are those in which the effect here becomes “large” as the index value increases, and the effect here becomes “increase” as the index value decreases. There is also.

提案欄108の提案フラグは、“採用”、“棄却”又は“追加”のうちのいずれかである。“棄却”は、その実験計画・最適化手法がユーザに対して出力(提案)されるべきものとして一旦決定されたが、見直しの結果、最終的に出力されなかったことを示す。“採用”は、その実験計画・最適化手法がユーザに対して出力されるべきものとして一旦決定され、予定通り最終的に出力されたことを示す。“追加”は、その実験計画・最適化手法がユーザに対して出力されるべきではないものとして一旦決定されたが、見直しの結果、最終的に出力されたことを示す。
図3の設計値管理情報31において、“#”は異なる設計値を省略的に表している。“♭”は、異なる性能値を省略的に表している。
The proposal flag in the proposal column 108 is one of “adopted”, “rejected”, and “added”. “Rejected” indicates that the experimental design / optimization method was once determined to be output (proposed) to the user, but was not finally output as a result of review. “Adopted” indicates that the experimental design / optimization method is once determined to be output to the user and finally output as scheduled. “Addition” indicates that the experimental design / optimization method was once determined as not to be output to the user, but was finally output as a result of the review.
In the design value management information 31 of FIG. 3, “#” represents different design values in an abbreviated manner. “♭” represents different performance values in an abbreviated manner.

図3において、二重線109より上のレコードは、実績データに関するレコードである。したがって、実験ID欄101〜実績/計画欄105にデータが記憶されている一方、実験計画・最適化手法欄106〜提案欄108は空欄である。   In FIG. 3, the record above the double line 109 is a record related to the performance data. Accordingly, data is stored in the experiment ID column 101 to the actual / plan column 105, while the experiment plan / optimization method column 106 to the proposal column 108 are blank.

二重線109より下のレコードは、実験支援装置1がユーザに対して提案する内容に関するレコードである。したがって、取得時点欄104は空欄である。性能値欄102は、データを記憶していてもよいし、空欄であってもよい。実験計画・最適化手法が類型1に属するもの(例えば勾配法)である場合、性能値欄102a〜102eに5つの“♭”が記憶されている。なお、そのうちの1又は複数は、最適化目標として所望の値になっている。実験計画・最適化手法が類型2に属するもの(例えば直交表L4)である場合、性能値欄102は空欄でもよい。なぜならば、類型2は、所望の性能値を出力させることが期待される設計値の候補を提案するに過ぎないからである。   A record below the double line 109 is a record related to the content proposed by the experiment support apparatus 1 to the user. Therefore, the acquisition time column 104 is blank. The performance value column 102 may store data or may be blank. When the experiment design / optimization method belongs to type 1 (for example, the gradient method), five “♭” are stored in the performance value columns 102a to 102e. One or more of them have a desired value as an optimization target. When the experiment design / optimization method belongs to type 2 (for example, orthogonal table L4), the performance value column 102 may be blank. This is because Type 2 only proposes candidate design values expected to output a desired performance value.

(特徴量情報)
図4に沿って、特徴量情報32を説明する。特徴量情報32においては、実績データ特徴量欄111に記憶された実績データ特徴量に関連付けて、実験計画・最適化手法欄112には実験計画・最適化手法が記憶されている。実験支援装置1は、実績データの特徴量に応じて実験計画・最適化手法を選択するとき、図4の特徴量情報32を使用してもよいし、図6で後記する実験計画・最適化手法決定処理手順を使用してもよい。特徴量情報32によれば、実験支援装置1は、実績データの特徴量が与えられれば、特徴量に相応しい1又は複数の実験計画・最適化手法を決定することができる。
(Feature information)
The feature amount information 32 will be described with reference to FIG. In the feature amount information 32, an experiment plan / optimization method is stored in the experiment plan / optimization method column 112 in association with the actual data feature amount stored in the result data feature amount column 111. The experiment support apparatus 1 may use the feature amount information 32 shown in FIG. 4 when selecting an experiment plan / optimization method according to the feature amount of the actual data, or the experiment plan / optimization described later in FIG. A technique determination procedure may be used. According to the feature amount information 32, the experiment support apparatus 1 can determine one or a plurality of experiment plans / optimization methods suitable for the feature amount when the feature amount of the performance data is given.

実績データ特徴量欄111は、小欄111a〜111gを有する。これらの小欄には、実績データの複数の特徴量が記憶されている。特徴量とは、ユーザが既に取得している(図3の二重線109より上のレコードの)設計値及び性能値に基づき、実験支援装置1が任意の方法で演算できる値である。“・・・”は、異なる特徴量を省略的に表している。   The actual data feature amount column 111 includes small columns 111a to 111g. In these small columns, a plurality of feature amounts of the performance data are stored. The feature amount is a value that can be calculated by the experiment support apparatus 1 by an arbitrary method based on the design value and the performance value that the user has already acquired (for the record above the double line 109 in FIG. 3). “...” Represents different feature amounts in an abbreviated manner.

実験点分布(欄111e)は、例えば、自己符号化器等によって次元圧縮された変数などの分布に関する情報である。応答曲面に関する情報(欄111f)は、例えば、線形回帰、カーネル回帰、グリギング、又は、アンサンブル学習、深層学習等の機械学習によって得られる、応答曲面に関する情報である。応答曲面から得られる情報(欄111g)は、例えば、グリギング又はブートストラップサンプリングを使用して推定できる予測値等の、応答曲面から得られる情報である。応答曲面に関する情報及び応答曲面から得られる情報は、応答曲面の学習時に推定されたパラメータを含む。   The experimental point distribution (column 111e) is information on the distribution of variables or the like that have been dimensionally compressed by a self-encoder or the like, for example. The information on the response surface (column 111f) is information on the response surface obtained by machine learning such as linear regression, kernel regression, grigging, ensemble learning, or deep learning. The information (column 111g) obtained from the response surface is information obtained from the response surface, such as a predicted value that can be estimated using grigging or bootstrap sampling. The information regarding the response surface and the information obtained from the response surface include the parameters estimated when learning the response surface.

(全体処理手順)
図5に沿って、全体処理手順を説明する。その説明のなかで、適宜図9、図10等を参照する。
ステップS201において、実験支援装置1の入力処理部21は、入力画面51(図9)を表示する。具体的には、入力処理部21は、入力画面51を出力装置13に表示する。
(Overall procedure)
The overall processing procedure will be described with reference to FIG. In the description, FIGS. 9 and 10 are referred to as appropriate.
In step S201, the input processing unit 21 of the experiment support apparatus 1 displays the input screen 51 (FIG. 9). Specifically, the input processing unit 21 displays the input screen 51 on the output device 13.

ステップS202において、入力処理部21は、実績データ等を受け付ける。具体的には、第1に、入力処理部21は、ユーザが入力画面51の実績データ欄52のファイル選択ボタン52a及び52bを操作して、それぞれ、設計値の実績データ及び性能値の実績データを選択するのを受け付ける。ここでは、図3の設計値管理情報31の二重線109より上のレコードの設計値及び性能値が、任意の外部サーバ(図1では図示せず)等に事前に格納されていることが前提となっている。その後、入力処理部21は、外部サーバ等から設計値及び性能値を取得し、取得した設計値及び性能値を使用して、設計値管理情報31を作成し、補助記憶装置15に記憶する。   In step S202, the input processing unit 21 receives performance data and the like. Specifically, first, the input processing unit 21 operates the file selection buttons 52a and 52b in the result data column 52 of the input screen 51, and the design value result data and the performance value result data, respectively. Accept to select. Here, the design value and performance value of the record above the double line 109 of the design value management information 31 in FIG. 3 may be stored in advance in an arbitrary external server (not shown in FIG. 1) or the like. It is a premise. Thereafter, the input processing unit 21 acquires design values and performance values from an external server or the like, creates design value management information 31 using the acquired design values and performance values, and stores them in the auxiliary storage device 15.

第2に、入力処理部21は、ユーザが最適化目標欄53のドロップダウンメニュー53aを操作して最適化目標とする性能値を選択するのを受け付ける。選択すべき性能値の候補は、図3の一次モード振動周波数、二次モード振動周波数、・・・、軸受44の温度である。さらに、入力処理部21は、ユーザがドロップダウンメニュー53bを操作して、“最大化”、“最小化”又は“特定の値”を選択するのを受け付ける。ユーザが“特定の値”を選択した場合、ユーザは、その値(又は値が取る範囲)を入力するものとする。さらに、入力処理部21は、ユーザが必要に応じて、最適化目標を追加するボタン53d又は最適化目標を削除するボタン53cを押下するのを受け付ける。   Secondly, the input processing unit 21 accepts that the user operates the drop-down menu 53a in the optimization target column 53 to select a performance value as an optimization target. The candidate performance values to be selected are the primary mode vibration frequency, the secondary mode vibration frequency,... Furthermore, the input processing unit 21 accepts that the user operates the drop-down menu 53b to select “maximize”, “minimize”, or “specific value”. When the user selects “specific value”, the user inputs the value (or the range that the value takes). Furthermore, the input processing unit 21 accepts that the user presses a button 53d for adding an optimization target or a button 53c for deleting an optimization target as necessary.

第3に、入力処理部21は、ユーザが制約欄54に設計値の制約を入力するのを受け付ける。ユーザは、設計値ごとに、下限テキストボックス54aに設計値の下限を入力し、上限テキストボックス54bに設計値の上限を入力する。ユーザは、必要に応じて、下限及び上限を入力するべき設計値を追加するボタン54d又は削除するボタン54cを押下する。ユーザは、不等式テキストボックス54eに、複数の設計値間の制約条件を定義する数式を入力する。ユーザは、必要に応じて、制約条件を追加するボタン54g又は削除するボタン54fを押下する。   Third, the input processing unit 21 accepts that the user inputs a design value constraint in the constraint column 54. For each design value, the user inputs the lower limit of the design value in the lower limit text box 54a and the upper limit of the design value in the upper limit text box 54b. The user presses the button 54d for adding or deleting the design value for inputting the lower limit and the upper limit as necessary. The user inputs a mathematical expression that defines a constraint between a plurality of design values in the inequality text box 54e. The user presses the button 54g for adding or deleting the restriction condition as necessary.

第4に、入力処理部21は、ユーザが実験可能回数テキストボックス55に、実験可能回数を入力するのを受け付ける。ここでの実験可能回数は、例えば、出力画面61(図10)において表示される直交表63bの最大行数である。
第5に、入力処理部21は、ユーザが実行ボタン56を押下するのを受け付ける。
Fourth, the input processing unit 21 accepts that the user inputs the number of possible experiments in the number of possible experiments text box 55. The number of possible experiments here is, for example, the maximum number of rows of the orthogonal table 63b displayed on the output screen 61 (FIG. 10).
Fifth, the input processing unit 21 accepts that the user presses the execution button 56.

ステップS203において、実験支援装置1の特徴量算出部22は、特徴量を算出する。具体的には、特徴量算出部22は、ステップS202において受け付けたデータに基づいて任意の方法で任意の特徴量を算出する。図4で前記したように、特徴量の例は、応答曲面に関する情報及び応答曲面から得られる情報を含む。そこで、特徴量算出部22は、実績データを学習データとして関数Fを学習する。例えば、特徴量算出部22は、多変数空間において重回帰分析を行い、関数Fを近似する曲面(応答曲面)を学習する。   In step S203, the feature amount calculation unit 22 of the experiment support apparatus 1 calculates a feature amount. Specifically, the feature amount calculation unit 22 calculates an arbitrary feature amount by an arbitrary method based on the data received in step S202. As described above with reference to FIG. 4, the feature amount example includes information on the response surface and information obtained from the response surface. Therefore, the feature amount calculation unit 22 learns the function F using the record data as learning data. For example, the feature amount calculation unit 22 performs multiple regression analysis in a multivariable space and learns a curved surface (response curved surface) that approximates the function F.

ステップS204において、実験支援装置1の実験計画・最適化手法選択部23は、実験計画・最適化手法を選択する。具体的には、第1に、実験計画・最適化手法選択部23は、特徴量情報32(図4)を使用して、又は、図6の実験計画・最適化手法決定処理手順(詳細後記)にしたがって、1又は複数の実験計画・最適化手法を選択する。なお、請求項の“実験計画又は最適化手法を選択する”ことは、ここで実験計画・最適化手法を選択することに相当する。このうち、“実験計画を選択する”ことには、例えば直交表のように、設計値の組合せを複数示す実験計画・最適化手法を選択することに相当する。“最適化手法を選択する”ことには、例えば勾配法のように、最適解であることが理論的に保証された設計値の組合せを示す実験計画・最適化手法を選択することに相当する。   In step S204, the experiment plan / optimization method selection unit 23 of the experiment support apparatus 1 selects an experiment plan / optimization method. Specifically, first, the experiment plan / optimization method selection unit 23 uses the feature amount information 32 (FIG. 4) or the experiment plan / optimization method decision processing procedure (details described later) in FIG. ) To select one or more experimental design / optimization methods. Note that “selecting an experiment design or optimization method” in the claims corresponds to selecting an experiment design / optimization method. Among these, “selecting an experimental design” corresponds to selecting an experimental design / optimization method that shows a plurality of combinations of design values, such as an orthogonal table. “Selecting an optimization method” corresponds to selecting an experimental design / optimization method that shows a combination of design values that are theoretically guaranteed to be optimal solutions, such as the gradient method. .

第2に、実験計画・最適化手法選択部23は、設計値管理情報31(図3)の新たなレコードを作成する。そして、実験計画・最適化手法選択部23は、新たなレコードの実験/計画欄105、実験計画・最適化手法欄106及び提案欄108に、それぞれ“計画”、ステップS204の“第1”において選択した実験計画・最適化手法、及び、“採用”を記憶する。   Second, the experiment design / optimization technique selection unit 23 creates a new record of the design value management information 31 (FIG. 3). Then, the experiment plan / optimization method selection unit 23 puts the new record in the experiment / plan column 105, the experiment plan / optimization method column 106, and the proposal column 108 in the “plan” and “first” in step S204, respectively. The selected experiment design / optimization method and “adopted” are stored.

ステップS205において、実験支援装置1の実験計画・最適化実行部24は、実験計画・最適化手法を実行する。具体的には、実験計画・最適化実行部24は、ステップS204の“第1”において選択された実験計画・最適化手法を、実験計画・最適化手法群33に記憶されているプログラムで実行する。   In step S205, the experiment plan / optimization execution unit 24 of the experiment support apparatus 1 executes the experiment plan / optimization technique. Specifically, the experiment plan / optimization execution unit 24 executes the experiment plan / optimization method selected in “first” in step S <b> 204 with a program stored in the experiment plan / optimization method group 33. To do.

実験計画・最適化実行部24は、設計値管理情報31のうち対応するレコードの設計値欄103に、ステップS205の実験計画・最適化の実行結果である実験案・設計値案を記憶する。   The experiment plan / optimization execution unit 24 stores the experiment plan / design value plan that is the execution result of the experiment plan / optimization at step S205 in the design value column 103 of the corresponding record in the design value management information 31.

実験計画・最適化実行部24は、ステップS205の処理をステップS204の“第1”において選択されたすべての実験計画・最適化手法について繰り返す。   The experiment plan / optimization execution unit 24 repeats the process of step S205 for all the experiment plans / optimization methods selected in “first” of step S204.

ステップS206において、実験支援装置1の実験計画・最適化効果算出部25は、効果を算出する。具体的には、第1に、実験計画・最適化効果算出部25は、それぞれの実験計画・最適化手法が出力した結果に基づいて、前記した性能改善度A、予測精度改善度I、スクリーニング達成度D、信頼度C及び必要実験回数のうちの少なくとも1つを算出する。
第2に、実験計画・最適化効果算出部25は、ステップS206の“第1”において取得した値に対して所定の閾値を適用し、効果として“大”又は“小”を取得する。
第3に、実験計画・最適化効果算出部25は、ステップS204の“第2”において作成したレコードの効果欄107に、“大”又は“小”を記憶する。
In step S206, the experiment plan / optimization effect calculation unit 25 of the experiment support apparatus 1 calculates an effect. Specifically, first, the experiment plan / optimization effect calculation unit 25 determines the performance improvement degree A, the prediction accuracy improvement degree I, the screening based on the results output by the respective experiment plans / optimization methods. At least one of achievement level D, reliability level C, and required number of experiments is calculated.
Second, the experiment design / optimization effect calculation unit 25 applies a predetermined threshold to the value acquired in “first” in step S206, and acquires “large” or “small” as the effect.
Third, the experiment design / optimization effect calculation unit 25 stores “large” or “small” in the effect column 107 of the record created in “second” in step S204.

ステップS207において、実験計画・最適化効果算出部25は、実験計画・最適化手法の削除を試みる。ステップS207の詳細は、図7の削除試行処理手順として後記する。結果的に、実験計画・最適化効果算出部25は、ステップS204の“第1”で選択された実験計画・最適化手法のうち、効果が“小”であるものに対して、提案フラグ“棄却”を立てる。例えば、既存実験点との距離が短いもの等に対して、提案フラグ“棄却”が立つことになる。もちろん、提案フラグ“棄却”を立てるべき実験計画・最適化手法が存在しない場合もある。   In step S207, the experiment plan / optimization effect calculation unit 25 attempts to delete the experiment plan / optimization technique. Details of step S207 will be described later as a deletion trial processing procedure in FIG. As a result, the experiment plan / optimization effect calculation unit 25 determines that the proposal flag “of” is selected for the experiment plan / optimization method selected in “first” in step S204 that has the effect “small”. Make a “rejection”. For example, the proposal flag “rejected” is set for an object having a short distance from the existing experimental point. Of course, there may be no experimental design / optimization method for setting the proposal flag “reject”.

ステップS208において、実験計画・最適化効果算出部25は、実験計画・最適化手法の追加を試みる。ステップS208の詳細は、図8の追加試行処理手順として後記する。結果的に、実験計画・最適化効果算出部25は、ステップS204の“第1”で選択された実験計画・最適化手法のすべてに提案フラグ“棄却”が立てられるのを防ぐ。   In step S208, the experiment design / optimization effect calculation unit 25 attempts to add an experiment design / optimization technique. Details of step S208 will be described later as an additional trial processing procedure in FIG. As a result, the experiment plan / optimization effect calculation unit 25 prevents the proposal flag “rejection” from being set for all the experiment plans / optimization methods selected in “first” in step S204.

ステップS209において、実験計画・最適化効果算出部25は、追加された実験計画・最適化手法があったか否かを判断する。具体的には、実験計画・最適化効果算出部25は、ステップS208において追加された実験計画・最適化手法があった場合(ステップS209“Yes”)、ステップS205に戻る。実験計画・最適化効果算出部25は、それ以外の場合(ステップS209“No”)、ステップS210に進む。   In step S209, the experiment plan / optimization effect calculation unit 25 determines whether there is an added experiment plan / optimization technique. Specifically, the experiment plan / optimization effect calculation unit 25 returns to step S205 when there is an experiment plan / optimization method added in step S208 (step S209 “Yes”). In other cases (step S209 “No”), the experiment plan / optimization effect calculation unit 25 proceeds to step S210.

ステップS210において、実験支援装置1の出力処理部26は、出力画面61(図10)を表示する。具体的には、第1に、出力処理部26は、出力画面61を出力装置13に表示する。
第2に、出力処理部26は、出力画面61の実験計画・最適化手法欄62bに、ユーザに提案する実験計画・最適化手法を表示する。設計値管理情報31(図3)のレコード(二重線109より下)のうち、提案フラグが“採用”又は“追加”であるレコードの実験計画・最適化手法が、実験計画手法欄62bに表示されることになる。
In step S210, the output processing unit 26 of the experiment support apparatus 1 displays the output screen 61 (FIG. 10). Specifically, first, the output processing unit 26 displays the output screen 61 on the output device 13.
Secondly, the output processing unit 26 displays the experiment plan / optimization technique proposed to the user in the experiment plan / optimization technique column 62 b of the output screen 61. Of the records of the design value management information 31 (FIG. 3) (below the double line 109), the experiment plan / optimization method of the record whose proposal flag is “adopted” or “added” is displayed in the experiment plan method column 62b. Will be displayed.

第3に、出力処理部26は、効果欄62cに、ステップS206の“第1”において算出した値のうちの4つ(性能改善度A、予測精度改善度I、スクリーニング達成度D、信頼度C)を表示する。これらの4つの値は、相互に比較しやすいように正規化されることが好ましい。
なお、出力処理部26は、図10に示した他にも、主成分分析、t−SNE等によって次元圧縮されたグラフ(例えば、設計値の散布図)を表示してもよい。
Third, the output processing unit 26 includes four values (performance improvement degree A, prediction accuracy improvement degree I, screening achievement degree D, reliability degree) calculated in “first” in step S206 in the effect column 62c. C) is displayed. These four values are preferably normalized to facilitate comparison with each other.
Note that the output processing unit 26 may display a graph (for example, a scatter diagram of design values) that is dimensionally compressed by principal component analysis, t-SNE, or the like, in addition to those illustrated in FIG. 10.

第4に、出力処理部26は、実験条件欄62aに、設計値を表示する。ここで表示される設計値は、多くの場合、最適化目標である所望の性能値を出力させることが理論的に保証されている設計値である(符号63a)。しかしながら、実験計画・最適化手法が類型2に属する直交表等である場合、ここで表示される設計値は、所望の性能値を出力させることが期待される複数の候補である(符号63b)。符号63bを見ると、設計値1として2つの水準(“10”及び“30”)並びに設計値2として2つの水準(“0”及び“0.1”)が、各列に2回ずつ表示されている(直交している)。因みに、符号63bの行数“4”は、必要実験回数として、“直交表(L4)”の効果の判定に使用されている。その後、全体処理手順を終了する。   Fourth, the output processing unit 26 displays the design value in the experiment condition column 62a. In many cases, the displayed design value is a design value that is theoretically guaranteed to output a desired performance value that is an optimization target (reference numeral 63a). However, when the experimental design / optimization method is an orthogonal table or the like belonging to type 2, the design values displayed here are a plurality of candidates expected to output desired performance values (reference numeral 63b). . Looking at the code 63b, two levels (“10” and “30”) as the design value 1 and two levels (“0” and “0.1”) as the design value 2 are displayed twice in each column. Has been (orthogonal). Incidentally, the number of rows “4” of the code 63b is used for determining the effect of the “orthogonal table (L4)” as the required number of experiments. Thereafter, the entire processing procedure is terminated.

(実験計画・最適化手法決定処理手順)
図6に沿って、実験計画・最適化手法決定処理手順を説明する。実験計画・最適化手法決定処理手順は、ステップS204(図5)の詳細である。
ステップS301において、実験支援装置1の実験計画・最適化手法選択部23は、特徴量を取得する。具体的には、実験計画・最適化手法選択部23は、ステップS203において算出された特徴量を取得する。
(Experimental design / optimization method decision process)
The experimental design / optimization method decision processing procedure will be described with reference to FIG. The experimental design / optimization technique determination processing procedure is the details of step S204 (FIG. 5).
In step S301, the experiment plan / optimization method selection unit 23 of the experiment support apparatus 1 acquires a feature amount. Specifically, the experiment design / optimization method selection unit 23 acquires the feature amount calculated in step S203.

ステップS302において、実験計画・最適化手法選択部23は、実験点数が閾値以上であるか否かを判断する。実験点数とは、実績データの数、すなわち、図3の設計値管理情報31における二重線109より上のレコードの数である。実験計画・最適化手法選択部23は、実験点数が所定の閾値以上である場合(ステップS302“Yes”)、ステップS304に進み、それ以外の場合(ステップS302“No”)、ステップS303に進む。   In step S302, the experiment plan / optimization technique selection unit 23 determines whether or not the number of experiments is equal to or greater than a threshold value. The number of experimental points is the number of actual data, that is, the number of records above the double line 109 in the design value management information 31 of FIG. The experiment plan / optimization method selection unit 23 proceeds to step S304 if the number of experiment points is equal to or greater than a predetermined threshold (step S302 “Yes”), and proceeds to step S303 otherwise (step S302 “No”). .

ステップS303において、実験計画・最適化手法選択部23は、直交表、最適計画法及び多元配置表を候補に追加する。具体的には、実験計画・最適化手法選択部23は、適当な設計値(引数)を与えたうえで、直交表、最適計画法及び多元配置表を実験計画・最適化手法の候補に追加する。その後、実験計画・最適化手法決定処理手順を終了する。   In step S303, the experiment design / optimization technique selection unit 23 adds the orthogonal table, the optimum design method, and the multi-way arrangement table to the candidates. Specifically, the experiment design / optimization method selection unit 23 gives an appropriate design value (argument) and then adds the orthogonal table, the optimum design method, and the multi-way layout table to the experiment design / optimization method candidates. To do. Thereafter, the experimental design / optimization method decision processing procedure is terminated.

ステップS304において、実験計画・最適化手法選択部23は、実験点密度が低い領域が閾値以上あるか否かを判断する。実験点密度とは、設計値の所定の領域における、実績データとなっている設計値の数である。実験計画・最適化手法選択部23は、実験点密度が低い領域が所定の閾値以上ある場合(ステップS304“Yes”)、ステップS305に進み、それ以外の場合(ステップS304“No”)、ステップS306に進む。   In step S304, the experiment design / optimization technique selection unit 23 determines whether or not there is an area with a low experiment point density that is equal to or greater than a threshold. The experimental point density is the number of design values that are actual data in a predetermined region of design values. The experiment plan / optimization method selection unit 23 proceeds to step S305 if the region having a low experiment point density is equal to or higher than the predetermined threshold (step S304 “Yes”), otherwise (step S304 “No”). The process proceeds to S306.

ステップS305において、実験計画・最適化手法選択部23は、直交表、最適計画法及び多元配置表を候補に追加する。具体的には、実験計画・最適化手法選択部23は、実験点密度が低い領域に対して、直交表、最適計画法及び多元配置表を実験計画・最適化手法の候補に追加する。その後、実験計画・最適化手法選択部23は、ステップS306に進む。   In step S305, the experiment design / optimization technique selection unit 23 adds the orthogonal table, the optimum design method, and the multi-way arrangement table to the candidates. Specifically, the experiment design / optimization technique selection unit 23 adds an orthogonal table, an optimum design method, and a multi-way arrangement table to the experiment design / optimization technique candidates for the region where the experiment point density is low. Thereafter, the experiment design / optimization technique selection unit 23 proceeds to step S306.

ステップS306において、実験計画・最適化手法選択部23は、応答曲面のパラメータが正常であるか否かを判断する。具体的には、実験計画・最適化手法選択部23は、特徴量として算出された応答曲面のパラメータが所定の正常範囲にある場合(ステップS306“Yes”)、ステップS308に進み、それ以外の場合(ステップS306“No”)、ステップS307に進む。   In step S306, the experiment design / optimization method selection unit 23 determines whether or not the response surface parameter is normal. Specifically, when the parameter of the response surface calculated as the feature amount is within a predetermined normal range (step S306 “Yes”), the experiment design / optimization method selection unit 23 proceeds to step S308, and otherwise If so (step S306 “No”), the process proceeds to step S307.

ステップS307において、実験計画・最適化手法選択部23は、最適計画法及び密度による計画法を実験計画・最適化手法の候補に追加する。その後、実験計画・最適化手法決定処理手順を終了する。
なお、“密度による計画法”は、特徴量として算出された対数実験点密度(logp)が最も小さい点に設計値のサンプルを追加する手法である。対数実験点密度は、式(5)によって定義される。式(5)の右辺の“2”は、ユークリッドノルムを示し、“t”は、ユーザが任意にその値を設定する正の定数である。
In step S307, the experiment design / optimization technique selection unit 23 adds the optimum design method and the density design method to the experiment design / optimization technique candidates. Thereafter, the experimental design / optimization method decision processing procedure is terminated.
The “planning method by density” is a method of adding a sample of design values to the point where the logarithmic experimental point density (logp) calculated as the feature quantity is the smallest. The log experimental point density is defined by equation (5). “2” on the right side of Expression (5) represents the Euclidean norm, and “t” is a positive constant that is arbitrarily set by the user.

Figure 2019185591
Figure 2019185591

ステップS308において、実験計画・最適化手法選択部23は、密度による計画法及びベイズ最適化を実験計画・最適化手法の候補に追加する。   In step S308, the experiment design / optimization technique selection unit 23 adds the density design method and Bayesian optimization to the experiment design / optimization technique candidates.

ステップS309において、実験計画・最適化手法選択部23は、応答曲面の分散が閾値以下であるか否かを判断する。具体的には、実験計画・最適化手法選択部23は、応答曲面の性能値yの分散が所定の閾値以下である場合(ステップS309“Yes”)、ステップS310に進み、それ以外の場合(ステップS309“No”)、実験計画・最適化手法決定処理手順を終了する。   In step S309, the experiment design / optimization method selection unit 23 determines whether or not the variance of the response surface is equal to or less than a threshold value. Specifically, the experiment design / optimization method selection unit 23 proceeds to step S310 if the variance of the performance value y of the response surface is equal to or smaller than a predetermined threshold (step S309 “Yes”), and otherwise ( In step S309 “No”), the experimental design / optimization method decision processing procedure ends.

ステップS310において、実験計画・最適化手法選択部23は、最適化目標が1つであるか否かを判断する。具体的には、実験計画・最適化手法選択部23は、最適化目標の設計値が1種類である場合(ステップS310“Yes”)、ステップS311に進み、それ以外の場合(ステップS310“No”)、ステップS312に進む。   In step S310, the experiment design / optimization technique selection unit 23 determines whether there is one optimization target. Specifically, the experiment plan / optimization method selection unit 23 proceeds to step S311 when the design value of the optimization target is one type (step S310 “Yes”), and otherwise (step S310 “No”). "), The process proceeds to step S312.

ステップS311において、実験計画・最適化手法選択部23は、勾配法を実験計画・最適化手法の候補に追加する。その後、実験計画・最適化手法決定処理手順を終了する。   In step S311, the experiment design / optimization technique selection unit 23 adds the gradient method to the experiment design / optimization technique candidates. Thereafter, the experimental design / optimization method decision processing procedure is terminated.

ステップS312において、実験計画・最適化手法選択部23は、線形計画法及び遺伝的アルゴリズムを実験計画・最適化手法の候補に追加する。線形計画法の例として、シンプレックス法が存在する。その後、実験計画・最適化手法決定処理手順を終了する。   In step S312, the experiment design / optimization technique selection unit 23 adds the linear programming method and the genetic algorithm to the experiment design / optimization technique candidates. As an example of linear programming, there is a simplex method. Thereafter, the experimental design / optimization method decision processing procedure is terminated.

実験計画・最適化手法決定処理手順において、実験計画・最適化手法を選択する基準は、実績データの特徴量である。しかしながら、より適切な実験計画・最適化手法を選択するためには、選択された実験計画・最適化手法の効果を検証する必要がある。実験支援装置1は、後記する削除試行処理手順において、効果が所定の基準に達しない実験計画・最適化手法を削除し、後記する追加試行処理手順において、効果が増加する可能性のある実験計画・最適化手法を追加する。   In the experiment planning / optimization technique determination processing procedure, the criterion for selecting the experiment planning / optimization technique is the feature value of the actual data. However, in order to select a more appropriate experiment design / optimization technique, it is necessary to verify the effect of the selected experiment design / optimization technique. The experiment support apparatus 1 deletes an experiment plan / optimization method whose effect does not reach a predetermined standard in a deletion trial process procedure to be described later, and an experiment plan in which the effect may increase in an additional trial process procedure to be described later・ Add optimization methods.

(削除試行処理手順)
図7に沿って、削除試行処理手順を説明する。削除試行処理手順は、ステップS207(図5)の詳細である。
ステップS401において、実験支援装置1の実験計画・最適化効果算出部25は、処理対象となる実験計画・最適化手法を取得する。具体的には、実験計画・最適化効果算出部25は、ステップS204の“第1”において選択された実験計画・最適化手法のうち未処理の1つを取得する。
(Deletion trial processing procedure)
The deletion trial processing procedure will be described with reference to FIG. The deletion trial processing procedure is the details of step S207 (FIG. 5).
In step S401, the experiment plan / optimization effect calculation unit 25 of the experiment support apparatus 1 acquires an experiment plan / optimization technique to be processed. Specifically, the experiment plan / optimization effect calculation unit 25 acquires an unprocessed one of the experiment plans / optimization methods selected in “first” in step S204.

ステップS402において、実験計画・最適化効果算出部25は、追加された実験計画・最適化手法であるか否かを判断する。具体的には、第1に、実験計画・最適化効果算出部25は、処理対象となっている実験計画・最適化手法が、後記するステップS502(図8)によって追加されたものであるか否かを調べる。
第2に、実験計画・最適化効果算出部25は、実験計画・最適化手法が追加されたものでない場合(ステップS402“No”)、ステップS403に進み、それ以外の場合(ステップS402“Yes”)、削除試行処理手順を終了する。
In step S402, the experiment plan / optimization effect calculation unit 25 determines whether the added experiment plan / optimization method is used. Specifically, first, the experiment plan / optimization effect calculation unit 25 determines whether the experiment plan / optimization technique to be processed has been added in step S502 (FIG. 8) to be described later. Check for no.
Second, the experiment plan / optimization effect calculation unit 25 proceeds to step S403 if the experiment plan / optimization method is not added (step S402 “No”), and otherwise (step S402 “Yes”). ”), The deletion trial processing procedure is terminated.

ステップS403において、実験計画・最適化効果算出部25は、既存実験点との距離が閾値以下であるか否かを判断する。具体的には、実験計画・最適化効果算出部25は、処理対象となっている実験計画・最適化手法において設定した設計値の代表値(例えば、平均値)と実績データの設計値との距離を算出する。そして、実験計画・最適化効果算出部25は、距離が所定の閾値以下である場合(ステップS403“Yes”)、ステップS405に進み、それ以外の場合(ステップS403“No”)、ステップS404に進む。   In step S403, the experiment plan / optimization effect calculation unit 25 determines whether the distance from the existing experiment point is equal to or less than a threshold value. Specifically, the experiment plan / optimization effect calculation unit 25 calculates a representative value (for example, an average value) of design values set in the experiment plan / optimization method to be processed and a design value of actual data. Calculate the distance. If the distance is equal to or smaller than the predetermined threshold (step S403 “Yes”), the experiment plan / optimization effect calculation unit 25 proceeds to step S405, and otherwise (step S403 “No”) proceeds to step S404. move on.

ステップS404において、実験計画・最適化効果算出部25は、効果が閾値以下であるか否かを判断する。具体的には、実験計画・最適化効果算出部25は、処理対象となっている実験計画・最適化手法の効果が所定の閾値以下である場合(ステップS404“Yes”)、ステップS405に進み、それ以外の場合(ステップS404“No”)、削除試行処理手順を終了する。   In step S404, the experiment plan / optimization effect calculation unit 25 determines whether the effect is equal to or less than a threshold value. Specifically, the experiment plan / optimization effect calculation unit 25 proceeds to step S405 when the effect of the experiment plan / optimization method to be processed is equal to or less than a predetermined threshold (step S404 “Yes”). In other cases (step S404 “No”), the deletion trial processing procedure is terminated.

ステップS405において、実験計画・最適化効果算出部25は、実験計画・最適化手法を削除する。具体的には、実験計画・最適化効果算出部25は、処理対象の実験計画・最適化手法に提案フラグ“棄却”を立てる。その後、削除試行処理手順を終了する。なお、実験計画・最適化効果算出部25は、削除試行処理手順を、未処理の実験計画・最適化手法ごとに繰り返す。   In step S405, the experiment plan / optimization effect calculation unit 25 deletes the experiment plan / optimization technique. Specifically, the experiment plan / optimization effect calculation unit 25 sets the proposal flag “rejected” to the experiment plan / optimization method to be processed. Thereafter, the deletion trial processing procedure is terminated. The experiment plan / optimization effect calculation unit 25 repeats the deletion trial processing procedure for each unprocessed experiment plan / optimization technique.

(追加試行処理手順)
図8に沿って、追加試行処理手順を説明する。追加試行処理手順は、ステップS208(図5)の詳細である。
ステップS501において、実験支援装置1の実験計画・最適化効果算出部25は、実験計画・最適化手法の候補が空であるか否かを判断する。“空である”とは、ステップS204の“第1”において実験計画・最適化手法が全く選択されなかった、又は、選択されたすべての実験計画・最適化手法に提案フラグ“棄却”が立てられたことを意味する。実験計画・最適化効果算出部25は、空である場合(ステップS501“Yes”)、ステップS502に進み、それ以外の場合(ステップS501“No”)、ステップS503に進む。
(Additional trial processing procedure)
The additional trial processing procedure will be described with reference to FIG. The additional trial processing procedure is the details of step S208 (FIG. 5).
In step S501, the experiment plan / optimization effect calculation unit 25 of the experiment support apparatus 1 determines whether the candidate experiment plan / optimization method is empty. “Empty” means that no experimental design / optimization method has been selected in “first” of step S204, or the proposal flag “reject” is set for all selected experimental designs / optimization methods. Means that The experiment plan / optimization effect calculation unit 25 proceeds to step S502 if it is empty (step S501 “Yes”), and proceeds to step S503 otherwise (step S501 “No”).

ステップS502において、実験計画・最適化効果算出部25は、密度による計画法及びベイズ最適化を実験計画・最適化手法の候補に追加する。ここで、ベイズ最適化を実験計画・最適化手法の候補に追加するのは、応答曲面のパラメータが正常範囲内にある場合である。ここで、実験計画・最適化効果算出部25は、設計値管理情報31(図3)の末尾に、提案フラグが“追加”であるレコードを記憶することになる。   In step S502, the experiment plan / optimization effect calculation unit 25 adds the density design method and Bayesian optimization to the experiment plan / optimization method candidates. Here, the Bayesian optimization is added to the experimental design / optimization method candidates when the parameters of the response surface are within the normal range. Here, the experiment plan / optimization effect calculation unit 25 stores a record in which the proposal flag is “added” at the end of the design value management information 31 (FIG. 3).

ステップS503において、実験計画・最適化効果算出部25は、処理対象となる実験計画・最適化手法を取得する。具体的には、実験計画・最適化効果算出部25は、ステップS204の“第1”において選択された、又は、ステップS502において追加された実験計画・最適化手法のうち未処理の1つを取得する。   In step S503, the experiment plan / optimization effect calculation unit 25 acquires the experiment plan / optimization technique to be processed. Specifically, the experiment plan / optimization effect calculation unit 25 selects one of the unprocessed experiment plans / optimization methods selected in “first” in step S204 or added in step S502. get.

ステップS504において、実験計画・最適化効果算出部25は、設計値等を調整する必要があるか否かを判断する。具体的には、第1に、実験計画・最適化効果算出部25は、ステップS205において使用された実験計画・最適化手法と、ステップS206において算出された効果との関係を、大小比較又は線形回帰やクリギングをはじめとする回帰により分析する。
第2に、実験計画・最適化効果算出部25は、分析の結果に基づき、設計値及び/又はパラメータを変更することによって効果が改善される可能性があるか否かを判断する。
第3に、実験計画・最適化効果算出部25は、可能性がある場合(ステップS504“Yes”)、ステップS505に進み、それ以外の場合(ステップS504“No”)、追加試行処理手順を終了する。
In step S504, the experiment design / optimization effect calculation unit 25 determines whether it is necessary to adjust the design value or the like. Specifically, first, the experiment design / optimization effect calculation unit 25 compares the relationship between the experiment design / optimization technique used in step S205 and the effect calculated in step S206 by comparing magnitude or linearity. Analyze by regression including regression and kriging.
Second, the experiment design / optimization effect calculation unit 25 determines whether there is a possibility that the effect is improved by changing the design value and / or the parameter based on the result of the analysis.
Third, the experiment plan / optimization effect calculation unit 25 proceeds to step S505 if there is a possibility (step S504 “Yes”), and otherwise performs the additional trial process procedure (step S504 “No”). finish.

ステップS505において、実験計画・最適化効果算出部25は、設計値等を調整する。具体的には、実験計画・最適化効果算出部25は、ステップS502において記憶された設計値管理情報31のレコードの設計値欄103を、ステップS504の“第2”における変更後の設計値で更新する。さらに、実験計画・最適化効果算出部25は、関数F(応答曲面)のパラメータをステップS504の“第2”における変更後のパラメータで更新する。
実験計画・最適化効果算出部25は、ステップS503〜S505の処理を未処理の実験計画・最適化手法ごとに繰り返す。その後、追加試行処理手順を終了する。
In step S505, the experiment plan / optimization effect calculation unit 25 adjusts the design value and the like. Specifically, the experiment design / optimization effect calculation unit 25 sets the design value column 103 of the record of the design value management information 31 stored in step S502 to the design value after the change in “second” in step S504. Update. Furthermore, the experiment design / optimization effect calculation unit 25 updates the parameter of the function F (response surface) with the parameter after the change in “second” in step S504.
The experiment plan / optimization effect calculation unit 25 repeats the processes of steps S503 to S505 for each unprocessed experiment plan / optimization technique. Thereafter, the additional trial processing procedure is terminated.

(本実施形態の効果)
本実施形態の実験支援装置の効果は以下の通りである。
(1)実験支援装置は、実験の実績データに基づき、適切な実験計画・最適化手法を効率的に出力できる。
(2)実験支援装置は、出力するべき実験計画・最適化手法を、その効果に応じて、削除又は追加することができる。
(3)実験支援装置は、実績データと実験計画・最適化手法との関係を予め決定しておくことができる。
(4)実験支援装置は、条件が整えば、最適化された設計値を出力することもでき、そうでなくても、直交表のように設計者の負担を軽減する実験計画を出力することができる。
(Effect of this embodiment)
The effects of the experiment support apparatus of this embodiment are as follows.
(1) The experiment support apparatus can efficiently output an appropriate experiment plan / optimization method based on the experiment result data.
(2) The experiment support apparatus can delete or add the experiment plan / optimization technique to be output according to the effect.
(3) The experiment support apparatus can determine in advance the relationship between the actual data and the experiment plan / optimization technique.
(4) The experiment support device can output an optimized design value if the conditions are met, or otherwise output an experiment plan that reduces the burden on the designer like an orthogonal table. Can do.

なお、本発明は前記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、前記した実施例は、本発明を分かり易く説明するために詳細に説明したものであり、必ずしも説明したすべての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。   In addition, this invention is not limited to an above-described Example, Various modifications are included. For example, the above-described embodiments have been described in detail for easy understanding of the present invention, and are not necessarily limited to those having all the configurations described. Further, a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added to the configuration of one embodiment. Further, it is possible to add, delete, and replace other configurations for a part of the configuration of each embodiment.

また、前記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウエアで実現してもよい。また、前記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウエアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、SSD(Solid State Drive)等の記録装置、又は、ICカード、SDカード、DVD等の記録媒体に置くことができる。
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしもすべての制御線や情報線を示しているとは限らない。実際には殆どすべての構成が相互に接続されていると考えてもよい。
Each of the above-described configurations, functions, processing units, processing means, and the like may be realized by hardware by designing a part or all of them with, for example, an integrated circuit. Each of the above-described configurations, functions, and the like may be realized by software by interpreting and executing a program that realizes each function by the processor. Information such as programs, tables, and files that realize each function can be stored in a recording device such as a memory, a hard disk, or an SSD (Solid State Drive), or a recording medium such as an IC card, an SD card, or a DVD.
In addition, the control lines and information lines are those that are considered necessary for the explanation, and not all the control lines and information lines on the product are necessarily shown. In practice, it may be considered that almost all the components are connected to each other.

1 実験支援装置
11 中央制御装置
12 入力装置
13 出力装置
14 主記憶装置
15 補助記憶装置
21 入力処理部
22 特徴量算出部
23 実験計画・最適化手法選択部
24 実験計画・最適化実行部
25 実験計画・最適化効果算出部
26 出力処理部
31 設計値管理情報
32 特徴量情報
33 実験計画・最適化手法群
DESCRIPTION OF SYMBOLS 1 Experiment support apparatus 11 Central control apparatus 12 Input apparatus 13 Output apparatus 14 Main memory 15 Auxiliary memory 21 Input processing part 22 Feature-value calculation part 23 Experiment plan and optimization method selection part 24 Experiment plan and optimization execution part 25 Experiment Planning / optimization effect calculation unit 26 Output processing unit 31 Design value management information 32 Feature quantity information 33 Experiment plan / optimization method group

Claims (5)

製品の設計値及び性能値を含む実績データ、最適化目標、及び、制約を受け付ける入力処理部と、
前記受け付けた実績データの特徴量を算出する特徴量算出部と、
前記算出した特徴量に応じて実験計画又は最適化手法を選択する実験計画・最適化手法選択部と、
前記選択した実験計画又は前記最適化手法を実行する実験計画・最適化実行部と、
前記実験計画又は前記最適化手法の実行結果に基づき前記実験計画又は前記最適化手法の効果を算出する実験計画・最適化効果算出部と、
前記実験計画又は前記最適化手法に関連付けて前記算出した効果を出力する出力処理部と、
を備えることを特徴とする実験支援装置。
An input processing unit that accepts performance data including product design values and performance values, optimization targets, and constraints;
A feature amount calculation unit for calculating a feature amount of the received record data;
An experiment plan / optimization method selection unit for selecting an experiment plan or an optimization method according to the calculated feature amount;
An experiment plan / optimization execution unit for executing the selected experiment plan or the optimization method;
An experiment plan / optimization effect calculation unit for calculating an effect of the experiment plan or the optimization method based on an execution result of the experiment plan or the optimization method;
An output processing unit for outputting the calculated effect in association with the experimental design or the optimization method;
An experiment support apparatus comprising:
前記実験計画・最適化効果算出部は、
前記効果が所定の基準に達しない実験計画又は最適化手法を削除し、
前記効果が増加する可能性のある実験計画又は最適化手法を追加すること、
を特徴とする請求項1に記載の実験支援装置。
The experiment design / optimization effect calculation unit
Delete experimental designs or optimization methods where the effect does not reach a predetermined standard,
Adding experimental designs or optimization techniques that may increase the effect,
The experiment support apparatus according to claim 1.
前記実験計画・最適化手法選択部は、
前記実績データと前記実験計画又は前記最適化手法との関係を予め定義した情報にしたがって、前記実験計画又は前記最適化手法を選択すること、
を特徴とする請求項2に記載の実験支援装置。
The experiment design / optimization method selection unit
Selecting the experiment plan or the optimization method according to information that defines in advance the relationship between the actual data and the experiment plan or the optimization method;
The experiment support apparatus according to claim 2.
前記実験計画・最適化手法選択部は、
最適計画法、密度による計画法、ベイズ最適化、勾配法、線形計画法、遺伝的アルゴリズム、直交表及び多元配置表のうちから前記実験計画又は前記最適化手法を選択すること、
を特徴とする請求項3に記載の実験支援装置。
The experiment design / optimization method selection unit
Selecting the experimental design or the optimization method from among optimal programming, density programming, Bayesian optimization, gradient method, linear programming, genetic algorithm, orthogonal table and multi-way table,
The experiment support apparatus according to claim 3.
実験支援装置の入力処理部は、
製品の設計値及び性能値を含む実績データ、最適化目標、及び、制約を受け付け、
前記実験支援装置の特徴量算出部は、
前記受け付けた実績データの特徴量を算出し、
前記実験支援装置の実験計画・最適化手法選択部は、
前記算出した特徴量に応じて実験計画又は最適化手法を選択し、
前記実験支援装置の実験計画・最適化実行部は、
前記選択した実験計画又は前記最適化手法を実行し、
前記実験支援装置の実験計画・最適化効果算出部は、
前記実験計画又は前記最適化手法の実行結果に基づき前記実験計画又は前記最適化手法の効果を算出し、
前記実験支援装置の出力処理部は、
前記実験計画又は前記最適化手法に関連付けて前記算出した効果を出力すること、
を特徴とする実験支援装置の実験支援方法。
The input processor of the experiment support device
Accept actual data including product design and performance values, optimization targets, and constraints,
The feature amount calculation unit of the experiment support apparatus includes:
Calculating a feature amount of the received performance data;
The experiment planning / optimization method selection unit of the experiment support device is:
Select an experiment plan or an optimization method according to the calculated feature amount,
The experiment plan / optimization execution unit of the experiment support apparatus is:
Performing the selected experimental design or the optimization technique;
The experiment plan / optimization effect calculation unit of the experiment support device
Calculate the effect of the experiment plan or the optimization method based on the execution result of the experiment plan or the optimization method,
The output processing unit of the experiment support apparatus includes:
Outputting the calculated effect in association with the experimental design or the optimization method;
The experiment support method of the experiment support apparatus characterized by this.
JP2018078363A 2018-04-16 2018-04-16 Experiment assisting device and experiment assisting method Pending JP2019185591A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018078363A JP2019185591A (en) 2018-04-16 2018-04-16 Experiment assisting device and experiment assisting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018078363A JP2019185591A (en) 2018-04-16 2018-04-16 Experiment assisting device and experiment assisting method

Publications (1)

Publication Number Publication Date
JP2019185591A true JP2019185591A (en) 2019-10-24

Family

ID=68341451

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018078363A Pending JP2019185591A (en) 2018-04-16 2018-04-16 Experiment assisting device and experiment assisting method

Country Status (1)

Country Link
JP (1) JP2019185591A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111784285A (en) * 2020-06-19 2020-10-16 上海华力集成电路制造有限公司 Engineering experiment scheme management system
KR20220029793A (en) * 2020-08-26 2022-03-08 쿠팡 주식회사 System and method for predicting an optimal stop point during an experiment test
WO2023080179A1 (en) * 2021-11-08 2023-05-11 株式会社レゾナック Design assitance device, design assitance method, and design assitance program
WO2023162959A1 (en) * 2022-02-22 2023-08-31 日東電工株式会社 Condition optimization device, condition optimization method, and program
JP7496498B2 (en) 2020-08-27 2024-06-07 パナソニックIpマネジメント株式会社 Information processing method, program, and information processing device
JP7496497B2 (en) 2020-08-27 2024-06-07 パナソニックIpマネジメント株式会社 Information processing method, program, and information processing device

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111784285A (en) * 2020-06-19 2020-10-16 上海华力集成电路制造有限公司 Engineering experiment scheme management system
KR20220029793A (en) * 2020-08-26 2022-03-08 쿠팡 주식회사 System and method for predicting an optimal stop point during an experiment test
KR102385928B1 (en) 2020-08-26 2022-06-23 쿠팡 주식회사 System and method for predicting an optimal stop point during an experiment test
JP7496498B2 (en) 2020-08-27 2024-06-07 パナソニックIpマネジメント株式会社 Information processing method, program, and information processing device
JP7496497B2 (en) 2020-08-27 2024-06-07 パナソニックIpマネジメント株式会社 Information processing method, program, and information processing device
WO2023080179A1 (en) * 2021-11-08 2023-05-11 株式会社レゾナック Design assitance device, design assitance method, and design assitance program
WO2023162959A1 (en) * 2022-02-22 2023-08-31 日東電工株式会社 Condition optimization device, condition optimization method, and program

Similar Documents

Publication Publication Date Title
JP2019185591A (en) Experiment assisting device and experiment assisting method
CN109657805B (en) Hyper-parameter determination method, device, electronic equipment and computer readable medium
US20180075357A1 (en) Automated system for development and deployment of heterogeneous predictive models
JP2017228086A (en) Machine learning management program, machine learning management method, and machine learning management device
CN111260073A (en) Data processing method, device and computer readable storage medium
JP7232122B2 (en) Physical property prediction device and physical property prediction method
KR102605481B1 (en) Method and Apparatus for Automatic Predictive Modeling Based on Workflow
Bolón-Canedo et al. Feature selection: From the past to the future
JP2021064049A (en) Calculator system and mathematical model generation support method
CN114365079A (en) Performance based system configuration as a pre-processing of system performance simulation
JP2017146888A (en) Design support device and method and program
JP2020144493A (en) Learning model generation support device and learning model generation support method
Ahmed et al. Optimum design of PI λ D μ controller for an automatic voltage regulator system using combinatorial test design
JP7479251B2 (en) Computer system and information processing method
Huang et al. Rough set theory: a novel approach for extraction of robust decision rules based on incremental attributes
US20230196109A1 (en) Non-transitory computer-readable recording medium for storing model generation program, model generation method, and model generation device
Huda et al. A group evaluation based binary PSO algorithm for feature selection in high dimensional data
JP2018077822A (en) Production support system and production support method
Hamidi et al. Constructing adaptive configuration dialogs using crowd data
Ghorbel et al. Smart adaptive run parameterization (SArP): enhancement of user manual selection of running parameters in fluid dynamic simulations using bio-inspired and machine-learning techniques
JP2020201611A (en) Planning system and method thereof
JP7196148B2 (en) Optimization support device and optimization support method
US20220027799A1 (en) Content classification method and classification model generation method
CN114239842A (en) Information processing apparatus, information processing system, and information processing method
JPWO2011016281A1 (en) Information processing apparatus and program for Bayesian network structure learning