特許 7676912 設計支援装置、設計支援方法及び設計支援プログラム

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

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

(11)【特許番号】

(24)【登録日】2025-05-07

(45)【発行日】2025-05-15

(54)【発明の名称】設計支援装置、設計支援方法及び設計支援プログラム

(51)【国際特許分類】

   G06F 30/20 20200101AFI20250508BHJP

   G06F 30/27 20200101ALI20250508BHJP

   G06F 111/06 20200101ALN20250508BHJP

【ＦＩ】

G06F30/20

G06F30/27

G06F111:06

【請求項の数】 10

(21)【出願番号】P 2021075151

(22)【出願日】2021-04-27

(65)【公開番号】P2022169242

(43)【公開日】2022-11-09

【審査請求日】2024-02-29

(73)【特許権者】

【識別番号】000004455

【氏名又は名称】株式会社レゾナック

(74)【代理人】

【識別番号】100088155

【弁理士】

【氏名又は名称】長谷川 芳樹

(74)【代理人】

【識別番号】100128381

【弁理士】

【氏名又は名称】清水 義憲

(74)【代理人】

【識別番号】100169454

【弁理士】

【氏名又は名称】平野 裕之

(74)【代理人】

【識別番号】100153040

【弁理士】

【氏名又は名称】川井 夏樹

(72)【発明者】

【氏名】花岡 恭平

【審査官】松浦 功

(56)【参考文献】

【文献】特開２０２０－０５２７３７（ＪＰ，Ａ）

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

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

【文献】特開２０１６－２００９０２（ＪＰ，Ａ）

【文献】特表２０１６－５２３４０２（ＪＰ，Ａ）

【文献】国際公開第２０１９／０８８１８５（ＷＯ，Ａ１）

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

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

【文献】特開２０２０－０３０６８３（ＪＰ，Ａ）

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

【文献】目的変数が複数のときに実験計画法のベイズ最適化(Bayesian Optimization, BO)が対応！ [online]，2019年03月25日，[検索日 2024.11.01]，インターネット，ＵＲＬ：https://datachemeng.com/bayesian_optimization_for_multiple_y/

【文献】GONZALEZ, J. et al.，Batch Bayesian Optimization via Local Penalization [online]，2015年10月15日，[検索日 2024.11.01]，インターネット，ＵＲＬ：https://arxiv.org/abs/1505.08052

【文献】米津智弘 外４名，コスト考慮型ベイズ最適化による複数目的関数最適化とその材料分野への応用，電子情報通信学会技術研究報告，一般社団法人電子情報通信学会，2017年06月16日，Vol. 117, No. 110，pp. 207-213

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

Ｇ０６Ｆ ３０／００ －３０／３９８

Ｇ０６Ｑ １０／０４ －１０／０４７

Ｇ０６Ｑ ５０／０４

Ｇｏｏｇｌｅ Ｓｃｈｏｌａｒ

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

【請求項1】

複数の設計パラメータからなる設計パラメータ群に基づいて作製される製品、仕掛品、半製品、部品又は試作品の設計において、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより設計パラメータの最適化を図る手法に適用するために、製品、仕掛品、半製品、部品又は試作品の特性を示す複数の特性項目が向上するような、前記複数の設計パラメータを求める設計支援装置であって、
作製済みの前記製品、前記仕掛品、前記半製品、前記部品又は前記試作品に関しての、前記設計パラメータ群と前記複数の特性項目のそれぞれの観測値とからなる実績データを複数取得するデータ取得部と、
前記設計パラメータ群に基づいて、目的変数としての前記特性項目の観測値を確率分布若しくはその近似又は代替指標として予測する予測モデルを、前記実績データに基づいて構築するモデル構築部と、
各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を一のサンプリング点として、前記目的変数群を所定点数サンプリングするサンプリング部と、
前記目的変数群に含まれる目的変数の数を次元数とし各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の目的変数群の評価値を算出する評価値算出部と、
各サンプリング点の前記評価値の分布に基づいて、前記設計パラメータ群を入力とし、前記評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する獲得関数評価部と、
前記獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する設計パラメータ群取得部と、
前記設計パラメータ群取得部により取得された前記設計パラメータ群を出力する出力部と、
を備える設計支援装置。

【請求項2】

前記評価値算出部は、前記目的変数群に含まれる目的変数の重み付け和を含む前記評価値を算出する、
請求項１に記載の設計支援装置。

【請求項3】

前記評価値算出部は、各目的変数に目標値が設定されている場合には、前記目的変数群に含まれる複数の目的変数のうち、目標値との差が最も大きい目的変数の該目標値に対する差を更に含む前記評価値を算出する、
請求項２に記載の設計支援装置。

【請求項4】

前記獲得関数評価部は、ＬＣＢ（Ｌｏｗｅｒ Ｃｏｎｆｉｄｅｎｃｅ Ｂｏｕｎｄ）、ＵＣＢ（Ｕｐｐｅｒ Ｃｏｎｆｉｄｅｎｃｅ Ｂｏｕｎｄ）、ＥＩ（Ｅｘｐｅｃｔｅｄ Ｉｍｐｒｏｖｅｍｅｎｔ）及びＰＩ（Ｐｒｏｂａｂｉｌｉｔｙ ｏｆ Ｉｍｐｒｏｖｅｍｅｎｔ）のうちのいずれかの獲得関数により前記獲得関数評価値を出力する、
請求項１～３のいずれか一項に記載の設計支援装置。

【請求項5】

前記設計パラメータ群取得部は、前記獲得関数評価値を最適化する一つの設計パラメータ群を取得する、
請求項１～４のいずれか一項に記載の設計支援装置。

【請求項6】

前記設計パラメータ群取得部は、複数の前記設計パラメータ群を所定のアルゴリズムにより取得する、
請求項１～４のいずれか一項に記載の設計支援装置。

【請求項7】

前記予測モデルは、前記設計パラメータ群を入力とし、前記観測値の確率分布を出力とする回帰モデルまたは分類モデルであり、
前記モデル構築部は、前記実績データを用いた機械学習により、前記予測モデルを構築する、
請求項１～６のいずれか一項に記載の設計支援装置。

【請求項8】

前記予測モデルは、ベイズ理論に基づく予測値の事後分布、アンサンブルを構成する予測器の予測値の分布、回帰モデルの予測区間及び信頼区間の理論式、モンテカルロドロップアウト、及び、異なる条件で複数個構築した予測器の予測の分布のうちのいずれか一つを用いて観測値の確率分布若しくはその近似又は代替指標を予測する機械学習モデルである、
請求項７に記載の設計支援装置。

【請求項9】

複数の設計パラメータからなる設計パラメータ群に基づいて作製される製品、仕掛品、半製品、部品又は試作品の設計において、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより設計パラメータの最適化を図る手法に適用するために、製品、仕掛品、半製品、部品又は試作品の特性を示す複数の特性項目が向上するような、前記複数の設計パラメータを求める設計支援装置における設計支援方法であって、
作製済みの前記製品、前記仕掛品、前記半製品、前記部品又は前記試作品に関しての、前記設計パラメータ群と前記複数の特性項目のそれぞれの観測値とからなる実績データを複数取得するデータ取得ステップと、
前記設計パラメータ群に基づいて、目的変数としての前記特性項目の観測値を確率分布若しくはその近似又は代替指標として予測する予測モデルを、前記実績データに基づいて構築するモデル構築ステップと、
各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を１のサンプリング点として、前記目的変数群を所定点数サンプリングするサンプリングステップと、
前記目的変数群に含まれる目的変数の数を次元数とし各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の目的変数群の評価値を算出する評価値算出ステップと、
各サンプリング点の前記評価値の分布に基づいて、前記設計パラメータ群を入力とし、前記評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する獲得関数評価ステップと、
前記獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する設計パラメータ群取得ステップと、
前記設計パラメータ群取得ステップにより取得された前記設計パラメータ群を出力する出力ステップと、
を有する設計支援方法。

【請求項10】

コンピュータを、複数の設計パラメータからなる設計パラメータ群に基づいて作製される製品、仕掛品、半製品、部品又は試作品の設計において、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより設計パラメータの最適化を図る手法に適用するために、製品、仕掛品、半製品、部品又は試作品の特性を示す複数の特性項目が向上するような、前記複数の設計パラメータを求める設計支援装置として機能させるための設計支援プログラムであって、
前記コンピュータに、
作製済みの前記製品、前記仕掛品、前記半製品、前記部品又は前記試作品に関しての、前記設計パラメータ群と前記複数の特性項目のそれぞれの観測値とからなる実績データを複数取得するデータ取得機能と、
前記設計パラメータ群に基づいて、目的変数としての前記特性項目の観測値を確率分布若しくはその近似又は代替指標として予測する予測モデルを、前記実績データに基づいて構築するモデル構築機能と、
各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を１のサンプリング点として、前記目的変数群を所定点数サンプリングするサンプリング機能と、
前記目的変数群に含まれる目的変数の数を次元数とし各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の目的変数群の評価値を算出する評価値算出機能と、
各サンプリング点の前記評価値の分布に基づいて、前記設計パラメータ群を入力とし、前記評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する獲得関数評価機能と、
前記獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する設計パラメータ群取得機能と、
前記設計パラメータ群取得機能により取得された前記設計パラメータ群を出力する出力機能と、
を実現させる設計支援プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示の一側面は、設計支援装置、設計支援方法及び設計支援プログラムに関する。

【背景技術】

【0002】

機械学習を活用した製品設計が研究されている。製品設計の一分野として、例えば、機能性材料の設計においては、例えば、実験及び作製済みの材料に関する原材料配合比と特性とのペアからなる学習データとを用いた機械学習により材料の特性を推定するモデルを構築し、未実験の原材料配合比に対する特性の予測が行われている。このような特性の予測により実験計画を立てることにより、効率的に材料の特性及び原材料配合比等のパラメータを最適化することが可能となり、開発効率の向上が図られている。また、このような最適化の手法として、ベイズ最適化が有効であることが知られており、ベイズ最適化を用いて設計値を出力する設計装置が知られている。

【先行技術文献】

【特許文献】

【0003】

【文献】特開２０２０－５２７３７号公報

【発明の概要】

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

【0004】

一方、材料等の製品開発においては、複数の目的変数（特性）が与えられた状況で、設計変数に応じて変化する複数の特性を向上させるために、複数の目的変数の最適化が行われる。これを多目的最適化という。目的変数間にトレードオフがある場合において、最適解（パレート解）は複数存在し、一つに定まらない。例えば、各目的変数に対して目標値が設定される場合に、最適なパレート解を得るために、多くのパレート解を求めて、設計目標に近いパレート解を選択するというアプローチを取ることが考えられる。しかしながら、このようなアプローチでは、多くの目的関数の評価が必要であり、その処理負荷が膨大となり現実的ではない。このような問題は、材料設計に限られず、製品設計全般に共通のものである。

【0005】

そこで、本発明は、上記問題点に鑑みてなされたものであり、製品、仕掛品、半製品、部品又は試作品の作製プロセスにおいて目的変数を構成する製品の特性及び設計変数の最適化を、より少ない実験回数により低負荷で可能とすることを目的とする。

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

【0006】

本開示の一側面に係る設計支援装置は、複数の設計パラメータからなる設計パラメータ群に基づいて作製される製品、仕掛品、半製品、部品又は試作品の設計において、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより設計パラメータの最適化を図る手法に適用するために、製品、仕掛品、半製品、部品又は試作品の特性を示す複数の特性項目が向上するような、複数の設計パラメータを求める設計支援装置であって、作製済みの製品、仕掛品、半製品、部品又は試作品に関しての、設計パラメータ群と複数の特性項目のそれぞれの観測値とからなる実績データを複数取得するデータ取得部と、設計パラメータ群に基づいて、目的変数としての特性項目の観測値を確率分布若しくはその近似又は代替指標として予測する予測モデルを、実績データに基づいて構築するモデル構築部と、各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を一のサンプリング点として、目的変数群を所定点数サンプリングするサンプリング部と、目的変数群に含まれる目的変数の数を次元数とし各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の目的変数群の評価値を算出する評価値算出部と、各サンプリング点の評価値の分布に基づいて、設計パラメータ群を入力とし、評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する獲得関数評価部と、獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する設計パラメータ群取得部と、設計パラメータ群取得部により取得された設計パラメータ群を出力する出力部と、を備える。

【0007】

本開示の一側面に係る設計支援方法は、複数の設計パラメータからなる設計パラメータ群に基づいて作製される製品、仕掛品、半製品、部品又は試作品の設計において、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより設計パラメータの最適化を図る手法に適用するために、製品、仕掛品、半製品、部品又は試作品の特性を示す複数の特性項目が向上するような、複数の設計パラメータを求める設計支援装置における設計支援方法であって、作製済みの製品、仕掛品、半製品、部品又は試作品に関しての、設計パラメータ群と複数の特性項目のそれぞれの観測値とからなる実績データを複数取得するデータ取得ステップと、設計パラメータ群に基づいて、目的変数としての特性項目の観測値を確率分布若しくはその近似又は代替指標として予測する予測モデルを、実績データに基づいて構築するモデル構築ステップと、各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を一のサンプリング点として、目的変数群を所定点数サンプリングするサンプリングステップと、目的変数群に含まれる目的変数の数を次元数とし各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の目的変数群の評価値を算出する評価値算出ステップと、各サンプリング点の評価値の分布に基づいて、設計パラメータ群を入力とし、評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する獲得関数評価ステップと、獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する設計パラメータ群取得ステップと、設計パラメータ群取得ステップにより取得された設計パラメータ群を出力する出力ステップと、を有する。

【0008】

本開示の一側面に係る設計支援プログラムは、コンピュータを、複数の設計パラメータからなる設計パラメータ群に基づいて作製される製品、仕掛品、半製品、部品又は試作品の設計において、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより設計パラメータの最適化を図る手法に適用するために、製品、仕掛品、半製品、部品又は試作品の特性を示す複数の特性項目が向上するような、複数の設計パラメータを求める設計支援装置として機能させるための設計支援プログラムであって、コンピュータに、作製済みの製品、仕掛品、半製品、部品又は試作品に関しての、設計パラメータ群と複数の特性項目のそれぞれの観測値とからなる実績データを複数取得するデータ取得機能と、設計パラメータ群に基づいて、目的変数としての特性項目の観測値を確率分布若しくはその近似又は代替指標として予測する予測モデルを、実績データに基づいて構築するモデル構築機能と、各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を一のサンプリング点として、目的変数群を所定点数サンプリングするサンプリング機能と、目的変数群に含まれる目的変数の数を次元数とし各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の目的変数群の評価値を算出する評価値算出機能と、各サンプリング点の評価値の分布に基づいて、設計パラメータ群を入力とし、評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する獲得関数評価機能と、獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する設計パラメータ群取得機能と、設計パラメータ群取得機能により取得された設計パラメータ群を出力する出力機能と、を実現させる。

【0009】

このような側面によれば、実績データに基づいて特性項目の観測値を予測する予測モデルが構築される。この予測モデルは、目的変数としての観測値を確率分布若しくはその近似又は代替指標として予測するので、各特性項目の予測モデルから得られる観測値の多次元分布に基づいて、任意の点数の目的変数群をサンプリングすることができる。各サンプリング点における目的変数群を要素とするベクトルに対して所定の演算を行うことにより、スカラー値により表現された各サンプリング点に関する評価値を得ることができる。そして、各サンプリング点の評価値の分布に基づいて、所定の獲得関数を用いて出力された獲得関数評価値の最適化により、次の実験等に好適な設計パラメータ群を得ることができる。従って、評価値を直接学習して獲得関数を構築する通常の方法と比べて、より高精度な機械学習モデルを得ることが可能になり、結果として実験等に採用される設計パラメータ群の好適度の向上により、実験回数の低減が可能となる。

【0010】

他の側面に係る設計支援装置では、評価値算出部は、目的変数群に含まれる目的変数の重み付け和を含む評価値を算出することとしてもよい。

【0011】

このような側面によれば、目的変数の重み付け和が評価値に含まれることにより、各サンプリング点の特徴が適切に表現された評価値を得ることができる。また、評価値算出部は、目的変数を入力として、所定の処理により算出される任意のスカラー値を評価値として算出してもよい。

【0012】

他の側面に係る設計支援装置では、評価値算出部は、各目的変数に目標値が設定されている場合には、目的変数群に含まれる複数の目的変数のうち、目標値との差が最も大きい目的変数の該目標値に対する差を更に含む評価値を算出することとしてもよい。

【0013】

このような側面によれば、最適化を最小化問題とする場合には、評価値算出部は、目的変数群に含まれる複数の目的変数のうち目標値を基準として最も数値が高い目的変数の該目標値に対する差を更に含む評価値を算出することとしてもよい。また、最適化を最大化問題とする場合には、評価値算出部は、目的変数群に含まれる複数の目的変数のうち目標値を基準として最も数値が低い目的変数の該目標値に対する差を更に含む評価値を算出することとしてもよい。このように算出された評価値に基づいて得られた獲得関数評価値の最適化により、最適化の一過程において目標値から最もかけ離れた目的変数を、目標値に近付けるために効果的な設計パラメータ群を得ることができる。

【0014】

他の側面に係る設計支援装置では、獲得関数評価部は、ＬＣＢ（Ｌｏｗｅｒ Ｃｏｎｆｉｄｅｎｃｅ Ｂｏｕｎｄ）、ＵＣＢ（Ｕｐｐｅｒ Ｃｏｎｆｉｄｅｎｃｅ Ｂｏｕｎｄ）、ＥＩ（Ｅｘｐｅｃｔｅｄ Ｉｍｐｒｏｖｅｍｅｎｔ）及びＰＩ（Ｐｒｏｂａｂｉｌｉｔｙ ｏｆ Ｉｍｐｒｏｖｅｍｅｎｔ）のうちのいずれかの獲得関数により、前記獲得関数評価値を出力することとしてもよい。

【0015】

このような側面によれば、評価値の向上に好適な設計パラメータ群を評価するのに好適な獲得関数評価値が出力される。

【0016】

他の側面に係る設計支援装置では、設計パラメータ群取得部は、獲得関数評価値を最適化する一つの設計パラメータ群を取得することとしてもよい。

【0017】

このような側面によれば、特性項目に関する評価値の向上が可能な設計パラメータ群を得ることができる。

【0018】

他の側面に係る設計支援装置では、設計パラメータ群取得部は、複数の設計パラメータ群を所定のアルゴリズムにより取得することとしてもよい。

【0019】

このような側面によれば、次の実験に供される複数の設計パラメータ群を容易に得ることができる。

【0020】

他の側面に係る設計支援装置では、予測モデルは、設計パラメータ群を入力とし、観測値の確率分布を出力とする回帰モデルまたは分類モデルであり、モデル構築部は、実績データを用いた機械学習により、予測モデルを構築することとしてもよい。

【0021】

このような側面によれば、予測モデルが所定の回帰モデルまたは分類モデルとして構築されるので、特性項目の観測値の確率分布若しくはその近似又は代替指標の取得が可能な予測モデル得られる。

【0022】

他の側面に係る設計支援装置では、予測モデルは、ベイズ理論に基づく予測値の事後分布、アンサンブルを構成する予測器の予測値の分布、回帰モデルの予測区間及び信頼区間の理論式、モンテカルロドロップアウト、及び、異なる条件で複数個構築した予測器の予測の分布のうちのいずれか一つを用いて観測値の確率分布若しくはその近似又は代替指標を予測する機械学習モデルであることとしてもよい。

【0023】

このような側面によれば、設計パラメータ群に基づく特性項目の観測値の確率分布若しくはその近似又は代替指標としての予測が可能な予測モデルが構築される。

【発明の効果】

【0024】

本開示の一側面によれば、製品、仕掛品、半製品、部品又は試作品の作製プロセスにおいて目的変数を構成する製品等の特性及び設計変数の最適化を、より少ない実験回数により低負荷で可能とする。

【図面の簡単な説明】

【0025】

【図1】実施形態に係る設計支援装置が適用される材料設計のプロセスの概要を示す図である。

【図2】実施形態に係る設計支援装置の機能構成の一例を示すブロック図である。

【図3】実施形態に係る設計支援装置のハードブロック図である。

【図4】作製済みの材料に関する設計パラメータ群の例を示す図である。

【図5】作製済みの材料に関する観測値の例を示す図である。

【図6】材料設計における特性項目及び設計パラメータの最適化のプロセスを示すフローチャートである。

【図7】実施形態に係る設計支援装置における設計支援方法の内容の一例を示すフローチャートである。

【図8】設計支援プログラムの構成を示す図である。

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

【0026】

以下、添付図面を参照しながら本発明の実施形態を詳細に説明する。なお、図面の説明において同一又は同等の要素には同一の符号を付し、重複する説明を省略する。

【0027】

図１は、実施形態に係る設計支援装置が適用される製品、仕掛品、半製品、部品又は試作品の設計のプロセスの一例である材料設計のプロセスの概要を示す図である。なお、以下において、「製品、仕掛品、半製品、部品又は試作品」を「製品等」と記載する。本実施形態の設計支援装置１０は、当該製品等の特性を示す複数の特性項目及び各特性項目の最適化を要するあらゆる製品等の設計のプロセスに適用できる。設計支援装置１０は、設計パラメータの決定と決定された設計パラメータに基づく製品、仕掛品、半製品、部品又は試作品の作製との繰り返しにより製品等の設計パラメータ及び目的変数の最適化を図る手法に適用されることができる。具体的には、設計支援装置１０は、材料の開発・設計の他に、例えば、自動車及び薬品等の製品の設計、薬品の分子構造の最適化等に適用できる。本実施形態では、上述のとおり、製品等の設計の一例としての材料設計の例により、設計支援装置１０による設計支援処理を説明する。

【0028】

図１に示されるように、設計支援装置１０による設計支援処理は、プラント及び実験室Ａ等における材料の作製及び実験に適用される。即ち、設定された設計パラメータ群ｘにより、プラント及び実験室Ａ等において材料が作製され、作製された材料に基づいて、材料の特性を示す複数の特性項目の観測値ｙが取得される。なお、プラント及び実験室Ａにおける材料作製及び実験は、シミュレーションであってもよい。この場合には、設計支援装置１０は、次のシミュレーションの実行のための設計パラメータ群ｘを提供する。

【0029】

設計支援装置１０は、設計パラメータ群ｘ及び設計パラメータ群ｘに基づいて作製された材料の複数の特性項目の観測値ｙからなる実績データに基づいて、複数の特性項目及び設計パラメータの最適化を行う。具体的には、設計支援装置１０は、作製済みの材料に関する設計パラメータ群ｘ及び観測値ｙに基づいて、次の作製及び実験のための、より好適な特性を得られる可能性がある設計パラメータ群ｘを出力する。

【0030】

例えば、本実施形態の設計支援装置１０は、材料製品の設計において、複数の設計変数をチューニングして、複数の特性を向上するという目的のために適用される。材料製品の設計の一例として、ある材料を複数のポリマー及び添加剤を混ぜて作製する場合において、設計支援装置１０は、各ポリマー及び添加剤の配合量等の設計パラメータ群を設計変数とし、特性項目である弾性率、熱膨張率の観測値を目的変数として、複数の特性に対して設定した一つの評価値を向上するような設計パラメータ群のチューニングに用いられる。

【0031】

図２は、実施形態に係る設計支援装置の機能構成の一例を示すブロック図である。設計支援装置１０は、複数の設計パラメータからなる設計パラメータ群に基づいて作製される材料の設計において、材料の特性を示す複数の特性に対して設定された一つの評価値を向上させるような、複数の設計パラメータを求める装置である。図２に示すように、設計支援装置１０は、プロセッサ１０１に構成された機能部、設計パラメータ記憶部２１及び観測値記憶部２２を含み得る。各機能部については後述する。

【0032】

図３は、実施形態に係る設計支援装置１０を構成するコンピュータ１００のハードウェア構成の一例を示す図である。なお、コンピュータ１００は、設計支援装置１０を構成しうる。

【0033】

一例として、コンピュータ１００はハードウェア構成要素として、プロセッサ１０１、主記憶装置１０２、補助記憶装置１０３、および通信制御装置１０４を備える。設計支援装置１０を構成するコンピュータ１００は、入力デバイスであるキーボード、タッチパネル、マウス等の入力装置１０５及びディスプレイ等の出力装置１０６をさらに含むこととしてもよい。

【0034】

プロセッサ１０１は、オペレーティングシステムおよびアプリケーションプログラムを実行する演算装置である。プロセッサの例としてＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）およびＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）が挙げられるが、プロセッサ１０１の種類はこれらに限定されない。例えば、プロセッサ１０１はセンサおよび専用回路の組合せでもよい。専用回路はＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）のようなプログラム可能な回路でもよいし、他の種類の回路でもよい。

【0035】

主記憶装置１０２は、設計支援装置１０等を実現するためのプログラム、プロセッサ１０１から出力された演算結果などを記憶する装置である。主記憶装置１０２は例えばＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）およびＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）のうちの少なくとも一つにより構成される。

【0036】

補助記憶装置１０３は、一般に主記憶装置１０２よりも大量のデータを記憶することが可能な装置である。補助記憶装置１０３は例えばハードディスク、フラッシュメモリなどの不揮発性記憶媒体によって構成される。補助記憶装置１０３は、コンピュータ１００を設計支援装置１０等として機能させるための設計支援プログラムＰ１と各種のデータとを記憶する。

【0037】

通信制御装置１０４は、通信ネットワークを介して他のコンピュータとの間でデータ通信を実行する装置である。通信制御装置１０４は例えばネットワークカードまたは無線通信モジュールにより構成される。

【0038】

設計支援装置１０の各機能要素は、プロセッサ１０１または主記憶装置１０２の上に、対応するプログラムＰ１を読み込ませてプロセッサ１０１にそのプログラムを実行させることで実現される。プログラムＰ１は、対応するサーバの各機能要素を実現するためのコードを含む。プロセッサ１０１はプログラムＰ１に従って通信制御装置１０４を動作させ、主記憶装置１０２または補助記憶装置１０３におけるデータの読み出しおよび書き込みを実行する。このような処理により、対応するサーバの各機能要素が実現される。

【0039】

プログラムＰ１は、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、半導体メモリなどの有形の記録媒体に固定的に記録された上で提供されてもよい。あるいは、これらのプログラムの少なくとも一つは、搬送波に重畳されたデータ信号として通信ネットワークを介して提供されてもよい。

【0040】

再び図２を参照して、設計支援装置１０は、データ取得部１１、モデル構築部１２、サンプリング部１３、評価値算出部１４、獲得関数評価部１５、設計パラメータ群取得部１６及び出力部１７を備える。設計パラメータ記憶部２１及び観測値記憶部２２は、図２に示されるように、設計支援装置１０に構成されてもよいし、設計支援装置１０からアクセス可能な他の装置として構成されてもよい。

【0041】

データ取得部１１は、作製済みの材料に関しての実績データを複数取得する。実績データは、設計パラメータ群と複数の特性項目のそれぞれの観測値とのペアからなる。設計パラメータ記憶部２１は、実績データにおける設計パラメータ群を記憶している記憶手段であって、例えば主記憶装置１０２及び補助記憶装置１０３等に構成されてもよい。観測値記憶部２２は、実績データにおける観測値を記憶している記憶手段である。

【0042】

図４は、設計パラメータ記憶部２１に記憶されている設計パラメータ群の例を示す図である。図４に示されるように、設計パラメータ記憶部２１は、１回目（ｔ＝１）から（Ｔ－１）回目（ｔ＝Ｔ－１）の材料作製における設計パラメータ群ｘ_ｔを記憶している。設計パラメータ群ｘは、一例として、原材料Ａの配合量、原材料Ｂの配合量及び設計パラメータＤ等を含んでもよく、設計パラメータの数Ｄに応じた次元数のベクトルデータを構成しうる。設計パラメータは、例示したものの他、例えば、分子構造及び画像等の非ベクトルデータ等であってもよい。また、複数の分子の種類から最適な分子を選ぶ問題を扱う場合には、設計パラメータは、複数の分子のうちの選択肢を示すデータであってもよい。

【0043】

図５は、観測値記憶部２２に記憶されている観測値ｙの例を示す図である。図５に示されるように、観測値記憶部２２は、１回目（ｔ＝１）から（Ｔ－１）回目（ｔ＝Ｔ－１）の材料作製において作製された材料の特性を示す複数の特性項目ｋ（ｋ＝１～Ｋ）の観測値ｙ_ｋ，ｔを記憶している。特性項目ｋは、一例として、ガラス転移温度、接着力及び特性項目Ｋを含んでもよい。設計パラメータ群ｘ_ｔと観測値ｙ_ｋ，ｔとのペアが実績データを構成する。

【0044】

設計支援装置１０は、１回目（ｔ＝１）から（Ｔ－１）回目（ｔ＝Ｔ－１）の材料作製における実績データに基づいて、各特性項目の観測値が向上するような、Ｔ回目の設計パラメータ群ｘ_Ｔを求める。

【0045】

モデル構築部１２は、実績データに基づいて予測モデルを構築する。予測モデルは、設計パラメータ群ｘに基づいて、目的変数としての特性項目ｋの観測値ｙ_ｋを確率分布若しくはその近似又は代替指標として予測するモデルである。予測モデルを構成するモデルは、観測値ｙ_ｋを確率分布若しくはその近似又は代替指標として予測可能なモデルであればよく、その種類は限定されない。

【0046】

例えば、予測モデルは、設計パラメータｘを入力とし、観測値ｙ_ｋの確率分布を出力とする回帰モデルであってもよい。予測モデルが回帰モデルである場合には、予測モデルは、例えば、線形回帰、ＰＬＳ回帰、ガウス過程回帰、ランダムフォレスト、及び、ニューラルネットワークといった回帰モデルのうちのいずれか一つにより構成されてもよい。モデル構築部１２は、実績データを用いた周知の機械学習の手法により、予測モデルを構築してもよい。モデル構築部１２は、実績データを予測モデルに適用して当該予測モデルのパラメータを更新する機械学習の手法により、予測モデルを構築してもよい。

【0047】

ガウス過程回帰として構築される予測モデルでは、教師データの説明変数を構成する実績データにおける設計パラメータ群ｘ及び目的変数を構成する観測値ｙ並びに予測対象の設計パラメータｘをモデルに入力することにより、観測値の確率分布が予測される。

【0048】

また、予測モデルは、ベイズ理論に基づく予測値の事後分布、アンサンブルを構成する予測器の予測値の分布、回帰モデルの予測区間及び信頼区間の理論式、モンテカルロドロップアウトにより得られる分布、及び、異なる条件で複数個構築した予測器の予測の分布のうちのいずれか一つを用いて観測値の確率分布若しくはその近似又は代替指標を予測する機械学習モデルであってもよい。

【0049】

観測値の確率分布、又はその代替指標の予測は、モデル固有の手法によって得ることができる。観測値の確率分布若しくはその近似又は代替指標は、ガウス過程回帰及びベイジアンニューラルネットワークであれば予測値の事後分布に基づいて、ランダムフォレストであれば、アンサンブルを構成する予測器の予測の分布に基づいて、線形回帰であれば予測区間及び信頼区間に基づいて、及び、ニューラルネットワークであればモンテカルロドロップアウトに基づいて得ることができる。但し、各機械学習モデルに対する観測値の分布またはその代替指標の算出方法は上記手法に限定されない。

【0050】

また、任意のモデルを、観測値の確率分布またはその代替指標を予測できるモデルに拡張することもできる。例えば、ブートストラップ法等で複数個のデータセットを構築し、それぞれに対して予測モデルを構築することで得られる、各モデルの予測値の分布を、観測値の確率分布の代替指標として用いるモデルが例として挙げることができる。但し、機械学習モデルを観測値の確率分布またはその代替指標を予測できるモデルに拡張する方法は上記手法に限定されない。

【0051】

また、モデル構築部１２は、予測モデルのハイパーパラメータを周知のハイパーパラメータチューニングの手法により、チューニングしてもよい。即ち、モデル構築部１２は、実績データにおける説明変数である設計パラメータ群ｘを表すベクトルと、目的変数である観測値ｙを用いた最尤推定により、予測モデルのハイパーパラメータを更新してもよい。

【0052】

また、予測モデルは、分類モデルにより構築されてもよい。予測モデルが分類モデルである場合には、モデル構築部１２は、実績データを用いた周知の確率分布の評価が可能な機械学習の手法により予測モデルを構築できる。

【0053】

このように、モデル構築部１２が所定の回帰モデルまたは分類モデルにより予測モデルを構築することにより、任意の設計パラメータ群ｘに基づいて、特性項目の観測値の確率分布の取得が可能となる。

【0054】

また、予測モデルは、一の特性項目の観測値を確率分布若しくはその近似又は代替指標として予測するシングルタスクモデル、または、複数の特性項目の観測値を確率分布若しくはその近似又は代替指標として予測するマルチタスクモデルであってもよい。このように、特性項目の性質に応じて適宜に構成されたマルチタスクモデルまたはシングルタスクモデルにより予測モデルを構築することにより、予測モデルによる観測値の予測の精度を向上できる。

【0055】

サンプリング部１３は、各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を一のサンプリング点として、所定点数の目的変数群のサンプリングを行う。

【0056】

具体的には、例えば目的変数としての観測値ｙ_ｋ（ｋ＝１～Ｋ）が、予測モデルに基づいて多次元正規分布に従い、かつ互いに無相関な場合に、以下のように表記される。
ｙ_ｋ～Ｎ（ｍ_ｋ（ｘ），σ_ｋ（ｘ）^２）
このような場合において、サンプリング部１３は、１つの設計パラメータ群ｘに基づいて、目的変数ｙ_ｋの確率分布から、複数のｙ_ｋ（ｋ＝１～Ｋ）をサンプリングする。

【0057】

より具体的には、サンプリング部１３は、第ｎ（ｎ：１～Ｎ）のサンプリング点として、目的変数群ｙ_ｋ，ｎ（ｋ＝１～Ｋ、ｎ＝１～Ｎ）をそれぞれサンプリングする。第ｎのサンプリング点の目的変数群ｙ_ｎは、以下に表記されるように、目的変数群に含まれる目的変数の数を次元数とし、各目的変数の値を要素とするベクトルを構成する。
ｙ_ｎ＝［ｙ_１，ｎ，ｙ_２，ｎ，．．，ｙ_ｋ，ｎ，．．，ｙ_Ｋ，ｎ］

【0058】

そして、サンプリング部１３は、所定の点数であるＮ点のサンプリング点に相当する目的変数群集合Ｙを取得する。
Ｙ＝［ｙ_１，ｙ_２，．．，ｙ_ｎ，．．，ｙ_Ｎ］
目的変数群ｙ_１，ｙ_２，．．，ｙ_ｎ，．．，ｙ_Ｎは、各々がベクトルを構成する。

【0059】

なお、以上の例示では、目的変数としての観測値ｙ_ｋが正規分布に従い、且つ互いに無関係であることを仮定したが、相関があってもよく、また、正規分布であることにも限定されず、その他の確率分布に従うこととしてもよい。また、サンプリング部１３によるサンプリングについて、目的変数ごとにサンプリングする例により説明したが、そのような例には限定されず、サンプリング部１３は、例えば、各特性項目の予測モデルにより規定される目的変数の多次元正規分布に基づいて目的変数群を一括でサンプリングしてもよい。

【0060】

また、サンプリング部１３は、事前にサンプリング及び保存した標準正規分布からのサンプリング点に基づいて、目的変数群ｙ_ｋ，ｎのサンプルを得ることとしてもよい。具体的には、サンプリング部１３は、予め、平均０，分散１の正規分布である標準正規分布からサンプリング点ｙ＿ｓｔｄ_ｋ，ｎ（ｎ＝１～Ｎ）をサンプリングしておき、以下の変換式により、目的変数群ｙ_ｋ，ｎのサンプルを得ることができる。
ｙ_ｋ，ｎ＝ｙ＿ｓｔｄ_ｋ，ｎ＊σ_ｋ（ｘ）+ｍ_ｋ（ｘ）

【0061】

なお、コンピュータに実装された設計支援装置１０においては、サンプリング部１３は、Ｎ点のサンプリングを一括で実施してもよいし、複数の設計パラメータ群ｘのそれぞれに対応するサンプリングを一括で実施してもよい。

【0062】

評価値算出部１４は、一の設計パラメータ群ｘ及び一のサンプリング点に対応する目的変数群の評価値を算出する。具体的には、上述のとおり、一のサンプリング点の目的変数群は、目的変数群に含まれる目的変数の数を次元数とし、各目的変数の値を要素とするベクトルを構成するので、評価値算出部１４は、目的変数群を表すベクトルを所定の演算によりスカラー化することにより、各サンプリング点の評価値を算出する。各サンプリング点の評価値は、より一般には、目的変数を入力として所定の処理により算出される任意のスカラー値であってもよい。

【0063】

評価値算出部１４は、一例として、以下の式（１）に示すように、スカラー化関数（Ｓｃａｌａｒｉｚｉｎｇ Ｆｕｎｃｔｉｏｎ：ＳＦ）により評価値ｖ_ｎを算出する。
ｖ_ｎ＝ＳＦ（ｙ_ｎ） ・・・（１）
評価値ｖ_ｎは、スカラー値を構成する。スカラー化関数ＳＦは、目的変数群ｙ_ｎに含まれる目的変数の重み付け和を算出する項を含んでもよい。また、各目的変数に目標値が設定されている場合において、最適化を最小化問題とする場合には、スカラー化関数ＳＦは、目的変数群に含まれる複数の目的変数のうち目標値を基準として最も数値が大きい目的変数の該目標値に対する差を含む項を更に含んでもよい。また、最適化を最大化問題とする場合には、スカラー化関数ＳＦは、目的変数群に含まれる複数の目的変数のうち目標値を基準として最も数値が小さい目的変数の該目標値に対する差を含む項を更に含んでもよい。

【0064】

評価値算出部１４は、例えば、各目的変数に目標値が設定され、最適化を最小化問題とする場合には、以下の式（２）により評価値ｖ_ｎを算出してもよい。

【数1】

式（２）のうちの、第１項は目的変数群ｙ_ｎに含まれる目的変数の重み付け和を算出する項であり、ｗ_ｋは目的変数ｙ_ｋに対する正の重みであって目的変数間のスケールの違いを調整するためのものであり、ρは第１項と第２項とのバランスを調整するための任意の定数である。また、式（２）のうちの第２項は、目的変数群ｙ_ｎに含まれる複数の目的変数のうちの目標値との差が最も大きい目的変数の、該目標値に対する差を含む項であり、ｇ_ｋは、目的変数ｙ_ｋの目標値である。

【0065】

式（２）の第１項に例示されるように、評価値ｖ_ｎが目的変数群ｙ_ｎに含まれる目的変数の重み付け和を含むことにより、このように算出された評価値ｖ_ｎを用いたベイズ最適化のプロセスにおいて、各目的変数がパレート解に近付くことを促進させることができる。また、式（２）の第２項に例示されるように、評価値ｖ_ｎが目的変数群ｙ_ｎに含まれる複数の目的変数のうちの目標値との差が最も大きい目的変数の、該目標値に対する差を含む項を含むことにより、ベイズ最適化のプロセスにおいて、目的変数が目標値に近付くことを促進させることができる。

【0066】

評価値算出部１４は、第１～第Ｎのサンプリング点の目的変数群ｙ_ｎのそれぞれの評価値ｖ_ｎを算出することにより、第１～第Ｎの評価値ｖ_ｎからなる評価値集合Ｖを得ることができる。
Ｖ＝［ｖ_１，ｖ_２，．．，ｖ_ｎ，．．，ｖ_Ｎ］
上述のとおり、評価値ｖ_１，ｖ_２，．．，ｖ_ｎ，．．，ｖ_Ｎは、それぞれスカラー値を構成する。

【0067】

なお、コンピュータに実装された設計支援装置１０においては、評価値算出部１４は、Ｎ点のサンプリング点のそれぞれに対応する評価値ｖ_１，ｖ_２，．．，ｖ_ｎ，．．，ｖ_Ｎの算出を一括で実施してもよい。

【0068】

獲得関数評価部１５は、各サンプリング点の評価値の分布に基づいて、設計パラメータ群を入力とし、評価値の向上に関する獲得関数評価値を所定の獲得関数により出力する。

【0069】

獲得関数評価部１５は、例えば、ＬＣＢ（Ｌｏｗｅｒ Ｃｏｎｆｉｄｅｎｃｅ Ｂｏｕｎｄ）といった周知の獲得関数を用いて獲得関数評価値を出力してもよい。ＬＣＢは、関数の出力を最小化する場合に用いられ、ＬＣＢの値を最小化することで好適な設計パラメータ群ｘが得られる。

【0070】

獲得関数をＬＣＢにより構築する場合、獲得関数評価部１５は、以下の式（３）のように獲得関数評価値Ａ（ｘ）を定義及び構築する。
Ａ（ｘ）＝ｍｖ（ｘ）－ａσｖ（ｘ） ・・・（３）
獲得関数評価部１５は、評価値集合Ｖに含まれる評価値ｖ_ｎの分布に基づいて、平均ｍｖ（ｘ）及び標準偏差σｖ（ｘ）を評価及び取得し、式（３）に示される獲得関数により獲得関数評価値を出力する。式（３）におけるａは、任意のパラメータである。上記獲得関数の式（３）は、設計パラメータ群ｘをパラメータとしたときの、次の実験におけるｖ_ｎの観測値が正規分布に従うと仮定した場合の信頼区間下限を表す。

【0071】

また、獲得関数評価部１５は、ＵＣＢ（Ｕｐｐｅｒ Ｃｏｎｆｉｄｅｎｃｅ Ｂｏｕｎｄ）、ＥＩ（Ｅｘｐｅｃｔｅｄ Ｉｍｐｒｏｖｅｍｅｎｔ）及びＰＩ（Ｐｒｏｂａｂｉｌｉｔｙ ｏｆ Ｉｍｐｒｏｖｅｍｅｎｔ）といった周知の関数により獲得関数評価値Ａ（ｘ）を出力してもよい。

【0072】

設計パラメータ群取得部１６は、獲得関数評価部１５により出力された獲得関数評価値Ａ（ｘ）の最適化により少なくとも一つの設計パラメータ群を取得する。

【0073】

一例として、設計パラメータ群取得部１６は、獲得関数の出力を最適化する少なくとも一つの設計パラメータ群ｘを取得してもよい。具体的には、設計パラメータ群取得部１６は、獲得関数評価部１５により出力された獲得関数評価値Ａ（ｘ）を目的変数とする最適化を実施し、最適解として設計パラメータ群ｘを取得する。

【0074】

また、一例として、設計パラメータ群取得部１６は、複数の設計パラメータ群を所定のアルゴリズムにより取得してもよい。具体的には、設計パラメータ群取得部１６は、獲得関数に対して、バッチベイズ最適化の手法を適用することにより、複数の設計パラメータ群を取得してもよい。バッチベイズ最適化の手法は、例えば、Ｌｏｃａｌ Ｐｅｎａｌｉｚａｔｉｏｎ等の手法であってもよいが、その手法は限定されない。

【0075】

出力部１７は、設計パラメータ群取得部１６により取得された設計パラメータ群を出力する。即ち、出力部１７は、１回目（ｔ＝１）から（Ｔ－１）回目（ｔ＝Ｔ－１）の材料作製における実績データに基づいて得られた設計パラメータ群を、Ｔ回目の材料の作製のための設計パラメータ群ｘ_Ｔとして出力する。

【0076】

また、設計パラメータ群取得部１６により複数の設計パラメータ群が取得される場合には、出力部１７は、（Ｔ－１）回目の次回以降のＮ回分の材料作製のための設計パラメータ群として、取得された設計パラメータ群を出力する。複数回分の材料作製のための設計パラメータ群は、同時の実験及び材料作製に供されてもよい。

【0077】

出力の態様は限定されないが、出力部１７は、例えば、所定の表示装置に表示させたり所定の記憶手段に記憶させたりすることにより、設計パラメータ群候補を出力する。

【0078】

図６は、材料設計における特性項目及び設計パラメータ群の最適化のプロセスを示すフローチャートである。

【0079】

ステップＳ１において、設計パラメータ群が取得される。ここで取得される設計パラメータ群は、初期の材料作製（実験）のためのものであって、任意に設定された設計パラメータ群であってもよいし、既に行われた実験等に基づいて設定された設計パラメータ群であってもよい。

【0080】

ステップＳ２において、材料作製が行われる。ステップＳ３において、作製された材料の特性項目の観測値が取得される。ステップＳ２における作製条件としての設計パラメータ群とステップＳ３において取得された各特性項目の観測値とのペアは、実績データを構成する。

【0081】

ステップＳ４において、所定の終了条件が充足されたか否かが判定される。所定の終了条件は、設計パラメータ群及び特性項目の観測値の最適化のための条件であって任意に設定されてもよい。最適化のための終了条件は、例えば、作製（実験）及び観測値の取得の所定回数への到達、観測値の目標値への到達及び最適化の収束等であってもよい。所定の終了条件が充足されたと判定された場合には、最適化のプロセスが終了される。所定の終了条件が充足されたと判定されなかった場合には、プロセスは、ステップＳ５に進む。

【0082】

ステップＳ５において、設計支援装置１０による設計支援処理が行われる。設計支援処理は、次の材料作製のための設計パラメータ群を出力する処理である。そして、プロセスは、再びステップＳ１に戻る。

【0083】

なお、ステップＳ１～Ｓ５により構成される処理サイクルの１サイクル目において、設計パラメータ群及び特性項目の観測値のペアが初期データとして複数得られる場合には、ステップＳ１～Ｓ４の処理は省略される。初期データが得られない場合には、ステップＳ１において、例えば実験計画法及びランダムサーチ等の任意の方法で得られた設計パラメータ群が取得される。処理サイクルの２サイクル目以降では、ステップＳ１において、ステップＳ５において出力された設計パラメータ群が取得される。

【0084】

図７は、実施形態に係る設計支援装置１０における設計支援方法の内容の一例を示すフローチャートであって、図６におけるステップＳ５の処理を示す。設計支援方法は、プロセッサ１０１に設計支援プログラムＰ１が読み込まれて、そのプログラムが実行されることにより、各機能部１１～１７が実現されることにより実行される。

【0085】

ステップＳ１１において、データ取得部１１は、作製済みの材料に関しての実績データを複数取得する。実績データは、設計パラメータ群と特性項目のそれぞれの観測値とのペアからなる。

【0086】

ステップＳ１２において、モデル構築部１２は、実績データに基づいて、予測モデルを構築する。

【0087】

ステップＳ１３において、サンプリング部１３は、予測モデルに基づいて、一つの設計パラメータ群ｘに基づいて各予測モデルから得られる観測値の多次元確率分布からサンプリングした複数の目的変数群を一のサンプリング点として、所定点数の目的変数群のサンプリングを行う。

【0088】

ステップＳ１４において、評価値算出部１４は、目的変数群の各目的変数の値を要素とするベクトルを所定の演算によりスカラー化することにより、各サンプリング点の評価値を算出する。

【0089】

ステップＳ１５において、獲得関数評価部１５は、各サンプリング点の評価値の分布に基づいて、所定の獲得関数評価値を出力する。

【0090】

ステップＳ１６において、設計パラメータ群取得部１６は、ステップＳ１５において獲得関数評価部１５により得られた獲得関数評価値の最適化により少なくとも一つの設計パラメータ群を取得する。

【0091】

ステップＳ１７において、出力部１７は、ステップＳ１６において設計パラメータ群取得部１６により取得された設計パラメータ群を、次の材料作製（ステップＳ１）のための設計パラメータ群として出力する。

【0092】

次に、コンピュータを、本実施形態の設計支援装置１０として機能させるための設計支援プログラムについて説明する。図８は、設計支援プログラムの構成を示す図である。

【0093】

設計支援プログラムＰ１は、設計支援装置１０における設計支援処理を統括的に制御するメインモジュールｍ１０、データ取得モジュールｍ１１、モデル構築モジュールｍ１２、サンプリングモジュールｍ１３、評価値算出モジュールｍ１４、獲得関数評価モジュールｍ１５、設計パラメータ群取得モジュールｍ１６及び出力モジュールｍ１７を備えて構成される。そして、各モジュールｍ１１～ｍ１７により、データ取得部１１、モデル構築部１２、サンプリング部１３、評価値算出部１４、獲得関数評価部１５、設計パラメータ群取得部１６及び出力部１７のための各機能が実現される。

【0094】

なお、設計支援プログラムＰ１は、通信回線等の伝送媒体を介して伝送される態様であってもよいし、図８に示されるように、記録媒体Ｍ１に記憶される態様であってもよい。

【0095】

以上説明した本実施形態の設計支援装置１０、設計支援方法及び設計支援プログラムＰ１によれば、実績データに基づいて特性項目の観測値を予測する予測モデルが構築される。この予測モデルは、目的変数としての観測値を確率分布若しくはその近似又は代替指標として予測するので、各特性項目の予測モデルから得られる観測値の多次元分布に基づいて、任意の点数の目的変数群をサンプリングすることができる。各サンプリング点における目的変数群を要素とするベクトルに対して所定の演算を行うことにより、スカラー値により表現された各サンプリング点に関する評価値を得ることができる。そして、各サンプリング点の評価値の分布に基づいて、所定の獲得関数を用いて出力された獲得関数評価値の最適化により、次の実験等に好適な設計パラメータ群を得ることができる。従って、評価値を直接学習して獲得関数を構築する通常の方法と比べて、より高精度な機械学習モデルを得ることが可能になり、結果として実験等に採用される設計パラメータ群の好適度の向上により、実験回数の低減が可能となる。

【0096】

以上、本発明をその実施形態に基づいて詳細に説明した。しかし、本発明は上記実施形態に限定されるものではない。本発明は、その要旨を逸脱しない範囲で様々な変形が可能である。

【符号の説明】

【0097】

Ｐ１…設計支援プログラム、ｍ１０…メインモジュール、ｍ１１…データ取得モジュール、ｍ１２…モデル構築モジュール、ｍ１３…サンプリングモジュール、ｍ１４…評価値算出モジュール、ｍ１５…獲得関数評価モジュール、ｍ１６…設計パラメータ群取得モジュール、ｍ１７…出力モジュール、１０…設計支援装置、１１…データ取得部、１２…モデル構築部、１３…サンプリング部、１４…評価値算出部、１５…獲得関数評価部、１６…設計パラメータ群取得部、１７…出力部、２１…設計パラメータ記憶部、２２…観測値記憶部。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】