JP7330712B2 - Material property prediction device and material property prediction method - Google Patents
Material property prediction device and material property prediction method Download PDFInfo
- Publication number
- JP7330712B2 JP7330712B2 JP2019022666A JP2019022666A JP7330712B2 JP 7330712 B2 JP7330712 B2 JP 7330712B2 JP 2019022666 A JP2019022666 A JP 2019022666A JP 2019022666 A JP2019022666 A JP 2019022666A JP 7330712 B2 JP7330712 B2 JP 7330712B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- parameter
- prediction
- property prediction
- material property
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、材料特性予測装置および材料特性予測方法に係り、特に、材料開発時の実験データに基づいて材料の特性を予測して、適切な材料を選別する場合に人の知見を活かした評価を行うのに好適な材料特性予測装置および材料特性予測方法に関する。 The present invention relates to a material property prediction device and a material property prediction method, and in particular, predicts material properties based on experimental data at the time of material development, and evaluates using human knowledge when selecting an appropriate material. It relates to a material property prediction device and a material property prediction method suitable for performing.
データ分析を行うための統計処理技術の発展に伴い、材料科学においてもデータ分析を行う需要が高まっている。特に、材料科学分野では、新素材の開発を効率的に行うために、次の実験の候補の選定を既知のデータを元にして行うスクリーニングとよばれる手法が知られている。 With the development of statistical processing technology for data analysis, there is a growing demand for data analysis in materials science as well. In particular, in the field of material science, a method called screening is known, in which candidates for the next experiment are selected based on known data in order to efficiently develop new materials.
スクリーニングの手法としては、各種実験のデータを情報システムに入力して機械学習にかけることによって、実験結果の予測のモデルを構築し、その予測にもとづき、材料設計に関する様々なパラメータを引数にとり材料の特性を返す関数を回帰分析によって求める方法がよく知られている。なお、以下の説明においては、そのような関数によって表現される材料特性を表すモデルを材料特性予測モデルと呼ぶこととする。 As a screening method, data from various experiments are input into an information system and subjected to machine learning to construct a model for predicting experimental results. A well-known method is to find a function that returns a characteristic by regression analysis. In the following description, a model representing material properties represented by such a function is called a material property prediction model.
化合物に関するデータをコンピュータにより分析するシステムは、例えば、特許文献1に記載がある。特許文献1には、化学データを2次元空間上のマップウィンドウ(図6)にマッピングして可視化するとともに、可視化されたデータに対して利用者が編集を行うことにより、2次元への圧縮結果が見やすくなるように調整してスクリーニングする方法が開示されている。
A system for analyzing data on compounds by a computer is described, for example, in
また、薬剤を発見するためのスクリーニングを高速に行える技術は、例えば、特許文献2に記載がある。特許文献2では、比較的小さいスクリーニングデータから構造的特徴と生物活性を結びつける法則を導くために、パラメータを部分的に用いた空間である部分空間上で、類似するデータによい材料特性を持つものが多い候補を優先的に実験するスクリーニングの方法が開示されている。
Also, a technique for high-speed screening for drug discovery is described in
材料開発においては、材料特性予測モデルの精度を高めることより、より的確に新材料の候補に対する有望性を見極めることができ、不要な実験を省略することで効率的な材料開発ができるようになると期待される。 In material development, by increasing the accuracy of the material property prediction model, it will be possible to more accurately determine the promise of new material candidates, and by omitting unnecessary experiments, it will be possible to develop materials efficiently. Be expected.
しかしながら、実際には、材料特性予測の精度を高めることは容易ではない。その大きな理由の一つが、高い精度で予測を行うのに必要なデータ数が多いという点である。その原因の一つが、材料のパラメータの多さにある。どのパラメータが材料特性に寄与するかがわからないと、すべてのパラメータを予測に用いることになる。そのような場合には、機械学習において一般に過学習とよばれる現象を引き起こし、精度低下を招く。これを補うために実験を増やすことも考えられるが、予測に要するデータ数は、材料のパラメータ数に対して指数関数的に増大するため、かなり多くの実験をせねばならない。するとその実験の労力を新材料の探索に振り向けたほうが早かったというような事態もありうる。 However, in practice, it is not easy to improve the accuracy of material property prediction. One of the major reasons for this is the large amount of data required to make highly accurate predictions. One of the causes is the large number of material parameters. Without knowing which parameters contribute to the material properties, all parameters are used in the prediction. In such a case, a phenomenon generally called over-learning occurs in machine learning, resulting in a decrease in accuracy. It is conceivable to increase the number of experiments to compensate for this, but the number of data required for prediction increases exponentially with the number of material parameters, so a considerable number of experiments must be performed. Then, there may be a situation where it would have been quicker to direct the effort of the experiment to the search for new materials.
そこで、予測のみに頼るのではなく、多くのパラメータをもつデータを適切に2次元空間に投影圧縮して可視化することによってスクリーニングをしやすくする技術がある。 Therefore, instead of relying only on prediction, there is a technique that facilitates screening by appropriately projecting and compressing data with many parameters into a two-dimensional space and visualizing it.
上記特許文献1に記載された分析システムは、化学データを2次元空間上のマッピングウィンドウに可視化するとともに、可視化されたデータに対して利用者が編集を行うことにより、2次元への圧縮結果が見やすくなるように調整してスクリーニングするものである。この方法により、目視でのスクリーニングはしやすくなる。しかしながら、これに加えて材料特性予測の結果を反映するには利用者がすべてを勘案することが必要であり、労力を要する。
The analysis system described in
また、上記特許文献2に記載された分析方法は、パラメータを部分的に用いた空間である部分空間上で、類似するデータによい材料特性を持つものが多い候補を優先的に実験する方法であった。この方法によれば、パラメータが多くてもスクリーニングができるが、データがあまりなくパラメータが多すぎる状況では、実際には相関がないのに相関があるかのように判断してしまうことがある。これは上述の過学習と同様の現象であり、パラメータの取捨選択なしに、精度向上するのは容易ではない点は依然として解決していない。
In addition, the analysis method described in
本発明の目的は、材料開発時の実験データに基づいて材料の特性を予測して、適切な材料を選別する場合に、実験計画のスクリーニングのための予測値評価において、人の知見を活かした評価を行うことのできる材料特性予測装置および材料特性予測方法を提供することにある。 The purpose of the present invention is to predict the characteristics of materials based on experimental data at the time of material development, and when selecting appropriate materials, in the prediction value evaluation for screening of experimental plans, human knowledge is utilized. An object of the present invention is to provide a material property prediction device and a material property prediction method that can be evaluated.
本発明の材料特性予測装置は、好ましくは、材料特性を予測して表示する材料特性予測装置であって、材料特性とその材料と関連付けられた材料パラメータとを格納する材料データテーブルと、材料パラメータの次元圧縮した結果を二次元座標として格納する材料データ可視化テーブルと、材料パラメータの各々に対応付けられたパラメータ重要度を格納するパラメータ重要度テーブルとを保持し、材料パラメータから次元圧縮により二次元座標として、材料データテーブルに格納する教師データなし次元圧縮部と、材料データ可視化テーブルに格納された二次元座標を材料データ点として、二次元表示し、材料データ点の位置変更を受け付ける2D可視化部と、受け付けた材料データ点の位置変更の情報に基づいて、回帰分析を行って、パラメータ重要度テーブルに格納されたパラメータ重要度を更新し、位置変更された材料データ点と位置変更されていない材料データ点の材料パラメータを次元圧縮して、材料データ可視化テーブルに格納される二次元座標として更新する半教師データあり次元圧縮部と、パラメータ重要度を各々の材料パラメータに加重した材料パラメータの類似度に基づいて、回帰分析を行って、材料特性を予測する予測評価部とを有するようにしたものである。 The material property prediction device of the present invention is preferably a material property prediction device for predicting and displaying material properties, comprising: a material data table storing material properties and material parameters associated with the materials; A material data visualization table that stores the result of dimensional compression as two-dimensional coordinates, and a parameter importance table that stores the parameter importance associated with each material parameter. A dimensional compression unit without teacher data that stores the coordinates in the material data table, and a 2D visualization unit that two-dimensionally displays the two-dimensional coordinates stored in the material data visualization table as material data points and accepts position changes of the material data points. and, based on the received repositioned material data point information, perform a regression analysis to update the parameter importance stored in the parameter importance table and compare the repositioned material data points with the unrepositioned material data points. A dimensional compression part with semi-supervised data that dimensionally compresses the material parameters of the material data points and updates them as two-dimensional coordinates stored in the material data visualization table, and a material parameter similarity that weights the parameter importance to each material parameter. and a predictive evaluator for predicting material properties by performing regression analysis based on the degree.
本発明によれば、材料開発時の実験データに基づいて材料の特性を予測して、適切な材料を選別する場合に、実験計画のスクリーニングのための予測値評価において、人の知見を活かした評価を行うことのできる材料特性予測装置および材料特性予測方法を提供することができる。 According to the present invention, when predicting the characteristics of materials based on experimental data during material development and selecting appropriate materials, human knowledge is utilized in predictive value evaluation for screening experimental plans. It is possible to provide a material property prediction device and a material property prediction method that can be evaluated.
以下、本発明に係る各実施形態を、図1ないし図17を用いて説明する。 Hereinafter, each embodiment according to the present invention will be described with reference to FIGS. 1 to 17. FIG.
〔実施形態1〕
以下、本発明に係る実施形態1を、図1ないし図13を用いて説明する。
[Embodiment 1]
A first embodiment according to the present invention will be described below with reference to FIGS. 1 to 13. FIG.
先ず、図1および図2を用いて実施形態1の材料特性予測装置の構成について説明する。
材料特性予測装置101は、利用者の操作に従い、材料に関する情報から材料特性を予測するための必要な情報を提示する装置である。
First, the configuration of the material property prediction device of
The material
材料特性予測装置101は、図1に示されるように、機能構成として、予測評価部112、教師データなし次元圧縮部113、2D(Dimension)可視化部115、半教師データあり次元圧縮部116、材料特性予測提示部117からなる。
As shown in FIG. 1 , the material
予測評価部112は、材料のパラメータに加重して予測精度を算出する機能部である。教師データなし次元圧縮部113は、多次元の材料のパラメータを与え、それを2次元の座標へと変換をする機能部である。ここで、「教師データなし」としたのは、後に出てくる「半教師データあり」と対比したものであり、与えられた原データが格納された材料データテーブルの材料パラメータの値に対して、そのまま次元圧縮することを意味する。2D可視化部115は、材料のパラメータの情報が変換された2次元の座標を可視化する機能部である。半教師データあり次元圧縮部116は、材料のパラメータとそれを次元圧縮した2次元の座標データのユーザ編集結果に基づいて、材料のパラメータを再び次元圧縮して、次元の座標を算出する機能部である。「半教師データあり」というのは、利用者のユーザインタフェースによって、材料データに関する情報を編集して、それにより生成されるデータを用いて新たに次元圧縮をすることを意味する。材料特性予測提示部117は、利用者に最終的に材料特性の予測を提示する機能部である。
The
また、材料特性予測装置101は、データとして、材料データテーブル240、材料データ可視化テーブル250、パラメータ重要度テーブル260、材料・座標マッピングテーブル270を保持する。
The material
材料データテーブル240は、材料の特性と材料に与えるパラメータに関する情報を格納するテーブルである。材料データ可視化テーブル250、材料データの情報変換した2次元の座標と、それに関連する情報を格納するテーブルである。パラメータ重要度テーブル260は、材料予測の算出過程で用いる各々のパラメータのパラメータ重要度(詳細は、後述)を格納するテーブルである。材料・座標マッピングテーブル270は、利用者が画面に表示された材料の2次元表示の点を編集したとき(詳細は、後述)、材料と2次元座標の対応付けを格納するテーブルである。 The material data table 240 is a table that stores information on material properties and parameters given to materials. A material data visualization table 250 is a table that stores information-converted two-dimensional coordinates of material data and information related thereto. The parameter importance table 260 is a table that stores the parameter importance (details will be described later) of each parameter used in the material prediction calculation process. The material/coordinate mapping table 270 is a table that stores the correspondence between the material and the two-dimensional coordinates when the user edits the two-dimensional display point of the material displayed on the screen (details will be described later).
次に、図2を用いて材料特性予測装置のハードウェア・ソフトウェア構成について説明する。 Next, the hardware/software configuration of the material property prediction device will be described with reference to FIG.
材料特性予測装置101は、パーソナルコンピュータ、ワークステーションのような一般的な情報処理装置により実現することができる。材料特性予測装置101は、例えば、図2に示されるように、CPU(Central Processing Unit)201、メインメモリ202、ネットワークI/F206、表示I/F203、入出力I/F204、補助記憶I/F205が、バスにより結合された形態になっている。
The material
CPU201は、材料特性予測装置101の各部を制御し、メインメモリ202に必要なプログラムをロードして実行する。
The
メインメモリ202は、通常、RAMなどの揮発メモリで構成され、CPU201が実行するプログラム、参照するデータが記憶される。
The
ネットワークI/F206は、材料特性予測装置101に接続されたネットワークを介して他の装置と通信するためのインタフェースである。
A network I/
表示I/F203は、LCD(Liquid Crystal Display)などの表示装置210を接続するためのインタフェースである。
A display I/
入出力I/F204は、入出力装置を接続するためのインタフェースである。図2の例では、キーボード221とポインティングデバイスのマウス222が接続されている。
The input/output I/
補助記憶I/F205は、HDD(Hard Disk Drive)230やSSD(Solid State Drive)などの補助記憶装置を接続するためのインタフェースである。
Auxiliary storage I/
HDD230は、大容量の記憶容量を有しており、本実施形態を実行するためのプログラムが格納されている。材料特性予測装置101には、予測評価プログラム231、教師データなし次元圧縮プログラム232、2D可視化プログラム233、半教師データあり次元圧縮プログラム234、材料特性予測提示プログラム235がインストールされている。
The
予測評価プログラム231、教師データなし次元圧縮プログラム232、2D可視化プログラム233、半教師データあり次元圧縮プログラム234、材料特性予測提示プログラム235は、それぞれ予測評価部112、教師データなし次元圧縮部113、2D可視化部115、半教師データあり次元圧縮部116、材料特性予測提示部117の機能を実行するプログラムである。
The
また、HDD230には、材料データテーブル240、材料データ可視化テーブル250、パラメータ重要度テーブル260、材料・座標マッピングテーブル270が格納されている。
The
次に、図3ないし図6を用いて材料特性予測装置の取り扱うデータ構造について説明する。 Next, the data structure handled by the material property prediction device will be described with reference to FIGS. 3 to 6. FIG.
材料データテーブル240は、材料に関する情報を格納するテーブルであり、図3に示されるように、材料ID501、材料名称502、材料特性503、材料パラメータ504の項目からなるレコードにより構成されるテーブルである。 The material data table 240 is a table that stores information about materials, and as shown in FIG. 3, it is a table composed of records consisting of items of material ID 501, material name 502, material properties 503, and material parameters 504. .
材料ID501には、材料を一意に識別する識別子が格納される。材料名称502には、材料の名称が格納される。材料特性503は、実験などによって判明している材料の特性を表す数値が格納される。材料特性503については、実験などによって判明している数値が格納されるが、その材料について実験が行われていないときには、空白であってもよい。材料パラメータ504には、材料の構造情報、製造条件などの材料の特性を算出するためのパラメータが格納される。材料パラメータ504は、例えば、複数のパラメータIDと、その値のペアとして、(パラメータ1,値1),(パラメータ2,値2),…のような形式で格納される。材料パラメータ504には、化学構造式など材料の構造や生成過程を示すあらゆる情報を含んでよい。
The material ID 501 stores an identifier that uniquely identifies the material. The material name 502 stores the name of the material. The material property 503 stores a numerical value representing the property of the material that has been found by experiment or the like. The material properties 503 store numerical values that have been found through experiments, but may be blank if no experiments have been conducted on the material. The material parameter 504 stores parameters for calculating material characteristics such as material structure information and manufacturing conditions. The material parameters 504 are stored in a format such as (
材料データ可視化テーブル250は、材料データを可視化するにあたっての材料に関する作業データを格納するテーブルであり、図4に示されるように、材料ID601、材料特性602、材料パラメータ603、編集有無フラグ604、予測精度評価値605、X座標606、Y座標607、材料特性予測値608の項目からなるレコードにより構成されるテーブルである。 The material data visualization table 250 is a table that stores work data related to materials for visualizing material data. As shown in FIG. The table is composed of records having items of accuracy evaluation value 605 , X coordinate 606 , Y coordinate 607 , and material property prediction value 608 .
材料ID601、材料特性602、材料パラメータ603には、それぞれ、材料データテーブル240の材料ID501、材料特性503、材料パラメータ504の値を複製した値が格納される。編集有無フラグ604には、その材料ID501で示す材料に対して、ユーザが材料データ2D表示画面(後述)により編集を行ったか否かを示すフラグが格納される。この編集有無フラグ604の初期値は、FALSE値である。予測精度評価値605には、予測評価部112の交差検証による予測精度評価値(詳細は、後述)が格納される。なお、材料特性602が空白のものについての予測は行われず、予測精度評価値605は空白としたまま格納される。X座標606、Y座標607には、材料パラメータ504の値を次元圧縮(詳細は後述)して、二次元にしたときの値が可能される。材料特性予測値608には、材料特性602がある場合はそれが格納され、そうでない場合は予測評価部112による材料特性の予測値を格納する。
A material ID 601, a material property 602, and a material parameter 603 store duplicated values of the material ID 501, material property 503, and material parameter 504 in the material data table 240, respectively. The edit presence/absence flag 604 stores a flag indicating whether or not the material indicated by the material ID 501 has been edited by the user using a
材料・座標マッピングテーブル270は、ユーザが材料データ2D表示画面(後述)により編集を行ったときの結果を格納するテーブルであり、図5に示されるように、材料ID901、X座標901、Y座標902の項目からなるレコードにより構成されるテーブルである。材料・座標マッピングテーブル270では、ユーザが材料データ2D表示画面によりデータ編集を行ったときに、編集が行われたデータに対応する編集後のX座標901とY座標902が、材料ID601に記載されたIDと対応づく形で記録される。
The material/coordinate mapping table 270 is a table that stores the result of editing by the user using a
パラメータ重要度テーブル260は、各パラメータに対するパラメータ重要度を格納するテーブルであり、パラメータID1101、パラメータ重要度1102の項目からなるレコードにより構成されるテーブルである。 The parameter importance level table 260 is a table that stores the parameter importance level for each parameter, and is a table composed of records consisting of items of parameter ID 1101 and parameter importance level 1102 .
パラメータ重要度1102には、パラメータID1101に格納されたパラメータID(図3の材料パラメータ504に記述されたパラメータID)に対するパラメータ重要度の値が格納される。なお、パラメータIDは、全ての材料に対して一意的な識別子をつけるものとする。パラメータ重要度は、詳細は後述するが、ユーザが材料データ2D表示画面より編集されたデータも考慮にいれて、新たに材料パラメータの次元圧縮をして、X座標、Y座標を求めるときに、回帰分析により予測関数を構築するときのカーネル関数の各材料パラメータに対するある種の重み付けとして評価される。
The parameter importance level 1102 stores the parameter importance level value for the parameter ID stored in the parameter ID 1101 (the parameter ID described in the material parameter 504 in FIG. 3). It should be noted that the parameter ID shall be a unique identifier for all materials. Although the details of the parameter importance will be described later, when the X coordinate and the Y coordinate are obtained by newly compressing the dimension of the material parameter in consideration of the data edited by the user from the
次に、図7ないし図13を用いて材料特性予測装置の処理について説明する。 Next, the processing of the material property prediction device will be described with reference to FIGS. 7 to 13. FIG.
先ず、図7を用いてユーザが材料特性予測装置を利用して材料予測をするときの処理の概要について説明する。
材料特性予測装置101は、材料データ可視化処理を実行し、利用者に材料データ2D表示画面312を表示する((I)初期マップ提示)。材料データ2D表示画面は、後に詳細に説明するが、材料データ可視化テーブル250に示した材料パラメータの情報を2次元空間上に圧縮した結果とそれと関連した情報を表示する画面である
次に、利用者に材料データ2D表示画面312から表示画面上でデータ操作を行い、それに伴って、材料特性予測装置101は、パラメータ重要度編集処理を行う((II)マップ編集・更新)。これによってパラメータ重要度が利用者の知見に沿った形に修正されていくことになる。この処理により、パラメータ重要度が更新される。そして、利用者は、満足するまでこのデータ操作を繰り返し実行し、その記録は材料・座標マッピングテーブル270として蓄積され、最終的には、材料データ可視化テーブルに反映される。
First, with reference to FIG. 7, an outline of processing when a user makes a material prediction using a material property prediction device will be described.
The material
最後に、材料特性予測装置101は、マップ編集・更新処理が終了すると、材料特性予測提示処理331を実行し、材料特性予測表示画面を利用者に提示する(((III)予測結果閲覧)。これにより、利用者は、上記の(II)マップ編集・更新のフェーズのデータ操作が反映された材料特性を把握することができる。
Finally, when the map editing/updating process is finished, the material
次に、図8を用いて材料データ可視化処理について説明する。
材料データ可視化処理では、先ず、予測評価部112が、材料データテーブル240に基づいて交差検証(詳細は、後述)によって予測精度評価値(詳細は、後述)を算出する(S100)。
Next, material data visualization processing will be described with reference to FIG.
In the material data visualization process, first, the
この処理を詳細に説明すれば、以下のようになる。予測評価部112は、材料データ可視化テーブル250の材料特性602が記述されているもののみを取得して、交差検証に基づいて、材料特性の予測精度評価値を算出する。材料特性602が記述されていないものについては、材料特性602が記述されているデータを学習データとして用いて材料特性予測関数を構築して(詳細は、後述)、材料特性の予測値を算出する。その結果は、材料データ可視化テーブル250の予測精度評価値605、材料特性予測値608に格納される。
A detailed description of this process is as follows. The
予測評価部112は、材料データテーブル240の材料パラメータ504に記述された材料パラメータの値を引数に与えると材料特性の予測値を返す関数(以下、「材料特性予測関数」とよぶ)を構成する機能を有する。この機能は、二つのデータの間の類似度に基づいて、材料特性を予測するものであり、公知のカーネル回帰分析法やサポートベクタ回帰などにより実現できる。予測評価部112は、データ間の類似度の計算において、材料パラメータの各変数にパラメータ重要度に基づき加重して用いることができる点に特徴があるが(詳細は後述)、この材料データ可視化処理においては、パラメータ重要度は初期値に相当するため、任意の方法で決定してよく、例えば、所定の値に固定してもよいし、ランダムに決定してもよい。
The
ここで、交差検証(Cross-validation)とは、統計学において、標本データを分割し、標本データの一部を学習用データ、残りを評価用データとして、その学習用データでまず解析して、評価用データでその解析のテストを行い、解析自身の妥当性の検証・確認に当てる手法である。本実施形態の交差検証においては、材料データの学習用データにより、材料特性予測関数を構築し、評価用データの材料パラメータを入力することによって、材料特性の予測値を算出する。その上で、評価用データに元々ある材料特性(以下、「真値」と呼ぶ)と当該予測値の差を誤差として求める。この誤差の絶対値が小さいほど、精度がよい基準の予測精度評価値が求まる。なお、この予測値から予測精度評価値を算出する方法としては、利用者に精度として認識可能な方法であればどのような方法を用いてもよい。例えば、誤差を真値で除算する値を用いると、真値の数が多いほど誤差が大きくなる効果を緩和して評価することができる。あるいは、誤差の絶対値の逆数を用いると、大きいほど精度がよいという直感にあった指標とすることができる。あるいは、真値と予測値の比の対数を求めることで、精度がすでに高いときは更なる向上は難しい点に配慮した(例えば、元々80%正しいものを89%正しい状態にすることに対して、元々99%正しいものを99.9%正しい状態にすることは十倍簡単とはいえない)予測精度評価値となる。 Here, cross-validation means that in statistics, sample data is divided, part of the sample data is used as learning data, the rest is used as evaluation data, and the learning data is first analyzed, This is a method of testing the analysis using the evaluation data to verify and confirm the validity of the analysis itself. In the cross-validation of this embodiment, a material property prediction function is constructed from the learning data of the material data, and the material parameter of the evaluation data is input to calculate the predicted value of the material property. Then, the difference between the material properties originally included in the evaluation data (hereinafter referred to as "true value") and the predicted value is obtained as an error. The smaller the absolute value of this error, the more accurate the reference prediction accuracy evaluation value is obtained. As a method for calculating the prediction accuracy evaluation value from this prediction value, any method may be used as long as the method is recognizable as accuracy by the user. For example, if a value obtained by dividing the error by the true value is used, the effect of increasing the error as the number of true values increases can be mitigated for evaluation. Alternatively, if the reciprocal of the absolute value of the error is used, it can be used as an intuitive index that the larger the value, the better the accuracy. Alternatively, by obtaining the logarithm of the ratio of the true value to the predicted value, we took into consideration the fact that it is difficult to improve further when the accuracy is already high (for example, to make the original 80% correct state 89% correct , it cannot be said that it is ten times easier to change a 99% correct state from the original to a 99.9% correct state).
次に、教師データなし次元圧縮部113は、材料パラメータを、多変量に対する次元削減手法により、X座標、Y座標を生成し、材料データ可視化テーブル250のX座標606、Y座標607に格納する(S101)。ここで、多変量に対する次元削減手法とは、多変量のデータを多次元空間上の点データとみなしたときに、その点データ同士の関係をなるべく損なわないように(統計分析に必要となる特徴をできるだけ保持するように)、低次元空間の点データとしてマッピングする方法であり、主成分分析法、部分空間法、独立成分分析法、非負値行列分解法、特異値行列分解法、自己組織化マップ法、t-SNE法などが知られている。
Next, the non-teaching data
最後に、以上の結果に基づいて、2D可視化部115は、材料データ2D表示画面312の表示情報を生成して、表示装置219により利用者に提示する(S102)。
Finally, based on the above results, the
ここで、図12を用いて材料データ2D表示画面312の詳細について説明する。
材料データ2D表示画面312は、図12に示されるように、材料データ表示フレーム701、材料詳細情報表示フレーム707、色分け指定ラジオボタン708、更新ボタン709、完了ボタン710の表示要素からなる。
Details of the
The
材料データ表示フレーム701は、材料パラメータを次元圧縮して、X座標、Y座標にマッピングした点データ(以下、「材料データ点」という)をプロットする領域である。材料データ表示フレーム701には、材料データがその特徴を把握しやすいように、材料データ点がプロットされる。この図12に示した例では、材料特性が空白すなわち未計測のものは、白い円702で表示され、また、予測精度評価値605が高いものは黒いひし形703、ついで高いものがグレーのひし形705、低いものが白いひし形704のように異なる形式で表現される。
A material
材料詳細情報表示フレーム707は、材料データの詳細な情報を表示する領域である。利用者がマウス222を用いてカーソル706を操作し、材料データ表示フレーム701内の材料データ点をクリックすると、その材料に関する詳細な情報が左下の材料詳細情報表示フレーム707に表示される。表示される情報は、材料データ可視化テーブル250の材料パラメータ603の情報に基づいて、作成されるものであり、その材料の特徴がわかりやすく表示される。
A material detailed
色分け指定ラジオボタン708は、色分けの指定を利用者が切り換えられるときに設定するボタンである。ここでは、その材料の予測精度、材料特性で切り換える例が示されている。
A color coding
更新ボタン709は、利用者の変更をデータに反映されるときに指定するボタンであり、完了710ボタンは、利用者が作業を終了し、この材料データ2D表示画面312を閉じるためのボタンである。
The
次に、図9を用いてパラメータ重要度編集処理について説明する。
パラメータ重要度編集処理では、材料特性予測装置は、先ず、図12に示した材料データ2D表示画面312を提示して利用者のデータ移動の入力を受け付ける(S200)。利用者は、マウス220のようなポインティングデバイスにより、カーソル706を操作して材料データ表示フレーム701内の各点をドラッグすることにより各点の位置を移動させることができる。
Next, parameter importance editing processing will be described with reference to FIG.
In the parameter importance editing process, the material property prediction device first presents the
上記の次元圧縮手法の特徴から明らかなように、基本的には材料データ表示フレーム701内の各点は、材料パラメータが比較的似た材料が近くに配置されることが期待される。また、もし材料パラメータが似ていると材料特性も近いという傾向があるなら、材料パラメータの類似性を用いて予測する予測評価部112による予測精度も高いことが期待される。逆に、図12に示した予測精度評価値の低い点(白いひし形704)があるとき、材料特性の近さは、材料パラメータの類似性が整合しないことになる。このことは、材料データ表示フレーム701内での、その点の配置が適切でないことにつながる。
As is clear from the characteristics of the dimension compression method described above, basically, it is expected that materials having relatively similar material parameters are placed close to each point in the material
本発明の材料特性の予測の手法は、この予測精度と材料データ表示フレーム701内での配置の妥当性の関連性に着目したものである。すなわち、予測精度を高めることと、材料データ表示フレーム701内での配置の妥当性を高めることは類似しているので、パラメータ重要度により両者を結びつけることによって、材料データ表示フレーム701内での配置の修正という簡便な方法で予測精度が高くなるようなパラメータ重要度の調整ができるのである。
The material property prediction method of the present invention focuses on the relationship between this prediction accuracy and the validity of the arrangement within the material
本実施形態では、利用者は、予測精度評価値の低いデータを中心に、材料特性の観点でより似ていると思われる別の点の周囲に移動するという修正を施すことができる。しかるのち、更新ボタン709を押下することにより、データ移動をした場合に(S201:YES)、材料特性予測装置101は、その移動結果を受け付ける(S202)。このとき、編集の結果は、ワークテーブルとして、材料・座標マッピングテーブル270に記録され、最終的に、材料データ可視化テーブル250のX座標606、Y座標607に格納される。そして、編集有無フラグ604がFALSE(偽)からTRUE(真)へと変更される。
In this embodiment, the user can make corrections by moving data with low prediction accuracy evaluation values to other points that are considered to be more similar in terms of material properties. After that, when the
次に、材料特性予測装置101は、このデータ編集情報に基づき、次元圧縮を再実行するための、半教師データあり次元圧縮実行処理が実行される(S203)。半教師データあり次元圧縮実行処理の詳細は、後に図10を用いて説明する。
Next, the material
次は、予測評価部112が、パラメータ重要度テーブル260を用いて予測精度を評価する処理が実行される(S204)。
Next,
ここで本実施形態における、次元圧縮関数を構成するための回帰分析の手法とパラメータ重要度の意義について詳しく述べる。本実施形態では、公知のカーネルガウシアンプロセス回帰法(KGPR:Kernel Gaussian Process Regression)を用いる。一般に、KGPRは二つのデータの間の類似度を示す関数であるカーネル関数を用いて予測関数を構築するが、本実施形態においては二つの材料パラメータの加重ユークリッド距離を引数にとる関数を用いる。この加重ユークリッド距離とは、二つの材料パラメータ{mi}、{ni}に対して、以下の(式1)で定義されるD({mi},{ni})を意味する。 Here, the method of regression analysis for constructing the dimension reduction function and the significance of the parameter importance in this embodiment will be described in detail. In this embodiment, a known Kernel Gaussian Process Regression (KGPR) method is used. In general, KGPR constructs a prediction function using a kernel function, which is a function that indicates the degree of similarity between two pieces of data. This weighted Euclidean distance means D({m i }, {n i } ) defined by the following (Equation 1) for two material parameters {m i }, {n i }.
D({mi},{ni})=Σiγi(mi-ni)2 …(式1) D({m i }, {n i })=Σ i γ i (m i −n i ) 2 (equation 1)
ここで、γiは材料パラメータのパラメータ重要度を示すパラメータであり、Σiはすべての材料パラメータにわたって和をとることを意味する。本実施形態では、この加重ユークリッド距離を用いた以下の(式2)で表されるガウス関数を以ってカーネル関数とする。
k({mi},{ni})=exp(-D({mi},{ni})) …(式2)
where γ i is a parameter indicating the parameter importance of material parameters, and Σ i means sum over all material parameters. In this embodiment, a Gaussian function represented by the following (Equation 2) using this weighted Euclidean distance is used as a kernel function.
k({m i }, {n i })=exp(−D({m i }, {n i })) (Formula 2)
公知のKGPRの回帰分析では、学習に用いるデータすべてについてカーネル関数の値を求め、それを用いて予測関数を構築する。それと同時に、学習に用いるデータと予測関数のあてはまり度合いを示す尤度を計算することができる。ここでカーネル関数のγiは、公知の最適化法(例えば、勾配降下法)によって尤度が最大になるように決定することができる。 In the well-known regression analysis of KGPR, kernel function values are obtained for all data used for learning, and prediction functions are constructed using the values. At the same time, it is possible to calculate the likelihood that indicates the degree of fit between the data used for learning and the prediction function. Here, γ i of the kernel function can be determined by a known optimization method (for example, gradient descent method) so as to maximize the likelihood.
予測評価部112が、パラメータ重要度テーブル260を用いて予測精度を評価する処理では、このKGPRによって予測を行い、その予測精度を評価する。すなわち、二つの材料パラメータ{mi}、{ni}に対して、上記の(式1)で定義されるD({mi},{ni})を用いた(式2)で表されるガウス関数をカーネル関数として材料特性の予測関数を求める。
この際、γiは、パラメータ重要度テーブル260のパラメータ重要度1102の値を用いる。
In the process of evaluating the prediction accuracy using the parameter importance table 260, the
At this time, the value of parameter importance 1102 of parameter importance table 260 is used for γ i .
後に、詳細に説明するように、S203の半教師データあり次元圧縮処理で、利用者の材料データ2D表示画面の材料データ点の編集結果が、パラメータ重要度テーブル260のパラメータ重要度1102に反映される。そのため、これによって、利用者の知見を反映した予測ができるようになる。なお、このパラメータ重要度テーブル260のパラメータ重要度1102をそのまま用いるのではなく、パラメータ重要度テーブル260のパラメータ重要度1102の値から大きく異ならないようにしてもよい。例えば、パラメータ重要度テーブル260のパラメータ重要度1102が大きい上位のいくつかはそのままにし、残りを上記の次元圧縮関数を構成する際と同様にして、勾配降下法で最適化する方法が考えられる。他にも、パラメータ重要度テーブル260のパラメータ重要度1102から大きく異ならない範囲で最適化するように制約つきの最適化を行う方法や、パラメータ重要度テーブル260のパラメータ重要度1102からずれた度合いをペナルティとして与える方法などがあるが、パラメータ重要度テーブル260のパラメータ重要度1102を反映できる方法であればどのような方法でもよい。
As will be described in detail later, in the dimension compression processing with semi-supervised data in S203, the editing result of the material data points on the
次に、図10を用いて半教師データあり次元圧縮処理の詳細について説明する。
これは、図9に示したパラメータ重要度編集処理のS203に示した処理である。
パラメータ重要度編集処理のS203では、先ず、材料データ可視化テーブル250の編集有無フラグ604がTRUEのレコードが取得される(S300)。
Next, details of the dimension compression processing with semi-supervised data will be described with reference to FIG.
This is the process shown in S203 of the parameter importance editing process shown in FIG.
In S203 of the parameter importance editing process, first, a record in which the editing presence/absence flag 604 of the material data visualization table 250 is TRUE is obtained (S300).
次に、それらのレコードのデータに対して、回帰分析を実行する(S301)。この回帰分析は、材料パラメータを引数に与えると、編集後のX座標901とY座標902を与える以下の(式3)の次元圧縮関数gを構築する処理である。
(X,Y)=g({mi},{γi}) …(式3)
ここで、(X,Y)は、それぞれ、編集後のX座標901とY座標902の値、{mi}は、材料データ可視化テーブル250の材料パラメータ603の値、{γi}は、材料パラメータ603に示される材料パラメータに対応するパラメータ重要度テーブル260のパラメータ重要度1102の値である。この手法は、上で説明した予測評価部112の材料特性予測関数を構築したときと同様の手法を用い、ただし、目的変数のみが異なる形となる。
Next, regression analysis is performed on the data of those records (S301). This regression analysis is a process of constructing a dimensional compression function g of the following (Equation 3) that gives an X coordinate 901 and a Y coordinate 902 after editing when material parameters are given as arguments.
(X, Y)=g({m i }, {γ i }) (Formula 3)
Here, (X, Y) are the values of the X coordinate 901 and the Y coordinate 902 after editing respectively, {m i } is the value of the material parameter 603 of the material data visualization table 250, and {γ i } is the material It is the value of the parameter importance level 1102 of the parameter importance level table 260 corresponding to the material parameter indicated by the parameter 603 . This method uses the same method as when constructing the material property prediction function of the
すなわち、材料パラメータmiに対応するパラメータ重要度をγiとして、二つの材料パラメータ{mi}、{ni}に対して、(式1)で定義されるD({mi},{ni})を用いた(式2)で表されるガウス関数をカーネル関数として材料特性の予測関数を求める(KGPRによる予測)。 That is, let γ i be the parameter importance corresponding to the material parameter m i , and D({m i } , { n i }) using the Gaussian function represented by (Equation 2) as the kernel function to determine the prediction function of the material properties (prediction by KGPR).
本実施形態にこの方法を適用すると、材料パラメータを編集後のX座標901とY座標902に投影するような次元圧縮関数を求めると同時に、そのときのパラメータ重要度γiを決定できる。 When this method is applied to the present embodiment, it is possible to obtain a dimensional compression function that projects the material parameters to the edited X-coordinate 901 and Y-coordinate 902, and at the same time determine the parameter importance γ i at that time.
上記の手順により次元圧縮関数およびパラメータ重要度が求まったら、材料データ可視化テーブル250の可視化向け材料データから、編集有無フラグ604がFALSEのものを取得し(S302)、編集有無フラグ604がTRUEのものとあわせて、次元圧縮関数に材料パラメータを入力することにより、新たなX座標606、Y座標607を求めて材料データ可視化テーブル250に設定する(S303)。この際、編集有無フラグ604がTRUEのものについては、計算の簡略化のために、編集後のX座標901とY座標902をそのまま用いるようにしてもよい。 When the dimension compression function and the parameter importance are obtained by the above procedure, from the material data for visualization in the material data visualization table 250, those with the edit flag 604 set to FALSE are acquired (S302), and those with the edit flag 604 set to TRUE are obtained. In addition, by inputting the material parameters into the dimensional compression function, new X coordinates 606 and Y coordinates 607 are obtained and set in the material data visualization table 250 (S303). At this time, if the editing flag 604 is TRUE, the X coordinate 901 and the Y coordinate 902 after editing may be used as they are to simplify the calculation.
この次元圧縮関数を求める過程で得られたパラメータ重要度は、パラメータ重要度テーブル260のパラメータ重要度1102に記録される(S304)。 The parameter importance obtained in the process of obtaining this dimension compression function is recorded in the parameter importance 1102 of the parameter importance table 260 (S304).
次に、図11を用いて材料特性予測提示処理の詳細について説明する。
材料特性予測提示処理では、材料データ可視化テーブル250から、材料データ(材料ID601、材料特性602、材料パラメータ603、材料特性予測値608)を取得し(S400)、パラメータ重要度テーブル260から、パラメータID1101、パラメータ重要度1102を取得し(S401)、パラメータ重要度に基づいて、材料特性を予測した結果を材料特性予測提示画面332に表示する(S402)。
Next, details of the material property prediction presentation processing will be described with reference to FIG. 11 .
In the material property prediction presentation process, the material data (material ID 601, material property 602, material parameter 603, material property prediction value 608) is acquired from the material data visualization table 250 (S400), and from the parameter importance table 260, the parameter ID 1101 , the parameter importance 1102 is obtained (S401), and the result of predicting the material properties is displayed on the material property
ここで、図13を用いて材料特性予測提示画面について説明する。
材料特性予測提示画面332は、図13に示されるように、材料特性予測値表示フレーム1301、パラメータ重要度表示フレーム1310、材料情報詳細表示フレーム1320、完了ボタン710の表示要素からなる。
Here, the material property prediction presentation screen will be described with reference to FIG. 13 .
The material property
材料特性予測値表示フレーム1301には、材料データ可視化テーブル250の材料ID601、材料特性602、材料特性予測値608、材料データテーブル240の材料名称502の値が表示される。材料特性予測値表示フレーム1301に表示される材料特性表示欄1302については、材料特性602の値か、それがない場合には、材料特性予測値608の値を用いて表示する。この際、材料特性予測値608が表示されているときにはそれと判別がつくように、例えば下線をつけるなどして表示することができる。また、利用者が表示されている材料のうちの一つを選択すると、その詳細を材料情報詳細表示フレーム1320に表示するようにもできる。
In the material property predicted
パラメータ重要度表示フレーム1310には、材料IDの材料に対応するパラメータ重要度が、材料データ可視化テーブル250の材料パラメータ603、パラメータ重要度テーブル260のパラメータ重要度1102を参照して表示される。
The parameter
以上説明してきたように、本実施形態の材料特性予測装置により、実験計画のスクリーニングのために、材料特性予測するときに、利用者の持つ知見を容易に予測に組み込んで評価できるようになる。これにより実験計画がたてやすくなり、ひいては少ない実験回数でよい材料を開発可能となる。 As described above, the material property prediction device of the present embodiment makes it possible to easily incorporate the user's knowledge into the prediction when predicting material properties for screening of experimental plans. This makes it easier to plan experiments, which in turn makes it possible to develop good materials with fewer experiments.
〔実施形態2〕
以下、本発明に係る実施形態2を、図14ないし図17を用いて説明する。
実施形態1では、利用者が材料データ2D表示画面312に表示される材料データ点を編集することにより、その知見を活かした材料予測が可能になる材料特性予測装置に関するものであった。実施形態1では、二次元空間上に材料データが可視化されているが、その距離感が把握できないため、編集がしにくいことがある。本実施形態では既存の何らかの画像と二次元座標を対応付けることによって距離感を把握できるように表示できる。
[Embodiment 2]
A second embodiment according to the present invention will be described below with reference to FIGS. 14 to 17. FIG.
The first embodiment relates to a material property prediction device that enables material prediction that makes use of the user's knowledge by editing the material data points displayed on the
そこで、本実施形態の材料特性予測装置は、実施形態1と同様の思想の材料特性予測装置で、材料データ2D表示画面312の画面に背景図面を表示することにより、データ間の関係を把握しやすくするようにしたものである。
Therefore, the material property prediction device of this embodiment is a material property prediction device having the same idea as that of the first embodiment, and by displaying a background drawing on the screen of the
以下の本実施形態の材料特性予測装置の説明では、実施形態1と異なる所を中心に説明する。
先ず、図14を用い実施形態2の材料特性予測装置の構成について説明する。
本実施形態の材料特性予測装置101は、実施形態1の図1に示した材料特性予測装置101に加えて、材料可視化背景テーブル1401が加わっている。そして、2D可視化部115がこの材料可視化背景テーブル1401を用いて画面を構成する。
In the following description of the material property prediction device of the present embodiment, differences from the first embodiment will be mainly described.
First, the configuration of the material property prediction device of
The material
次に、図15を用いて材料可視化背景テーブル1401のデータ構造について説明する。
材料可視化背景テーブル1401には、材料データ2D表示画面312の材料データ表示フレーム701の背景画像の画像データを格納するテーブルであり、図15に示されるように、左下X座標1601、左下Y座標1602、右上X座標1603、右上Y座標1604、背景画像データ1605のからなるレコードにより構成されるテーブルである。
Next, the data structure of the material visualization background table 1401 will be described using FIG.
The material visualization background table 1401 is a table that stores the image data of the background image of the material
左下X座標1601、左下Y座標1602、右上X座標1603、右上Y座標1604は、それぞれ対応する背景画像データ1605の座標点を表している。 A lower left X coordinate 1601, a lower left Y coordinate 1602, an upper right X coordinate 1603, and an upper right Y coordinate 1604 represent coordinate points of the corresponding background image data 1605, respectively.
次に、図16を用いて実施形態2の材料データ可視化処理について説明する。
実施形態1との差異は、教師データなし次元圧縮部113によって材料パラメータの次元圧縮によって、材料データ点のX,Y座標が定まった後、材料可視化背景テーブル1401のデータを参照して座標原点と縮尺を整合させる処理(S110)がある点である。
Next, the material data visualization processing of the second embodiment will be described with reference to FIG. 16 .
The difference from the first embodiment is that after the X and Y coordinates of the material data point are determined by the dimension compression of the material parameter by the
2D可視化部115は、材料可視化背景テーブル1401に画像データが格納された背景画像の座標の中央を求め、また、材料データ可視化テーブル250のX座標606、Y座標607の重心を求め、両者の差を計算する。この差を材料データ可視化テーブル250のX座標606、Y座標607から減算することで、背景画像の座標と材料データ可視化テーブル250が整合するようにできる。この補正の方法については、背景画像に適切に重なるような方法であれば、どのような方法を用いてもよい。
The
次に、図17を用いて実施形態2の材料データ2D表示画面について説明する。
本実施形態の材料データ2D表示画面の材料データ表示フレーム701には、背景画像として、日本地図が重ねて表示される。このように、利用者になじみのある直感的に捉えやすい背景画像を表示することにより、材料データの間の距離を直感的に把握しやすくできる。
Next, the
A map of Japan is superimposed and displayed as a background image in the material
101…材料特性予測装置
112…予測評価部
113…教師データなし次元圧縮部
115…2D可視化部
116…半教師データあり次元圧縮部
117…材料特性予測提示部
240…材料データテーブル
250…材料データ可視化テーブル
260…パラメータ重要度テーブル
270…材料・座標マッピングテーブル
Reference numerals 101: material property prediction device 112: prediction evaluation unit 113: no teacher data dimension compression unit 115: 2D visualization unit 116: semi-supervised data presence dimension compression unit 117: material property prediction presentation unit 240: material data table 250: material data visualization Table 260: Parameter importance table 270: Material/coordinate mapping table
Claims (12)
材料特性とその材料と関連付けられた材料パラメータとを格納する材料データテーブルと、
前記材料パラメータの次元圧縮した情報を二次元座標として格納する材料データ可視化テーブルと、
前記材料パラメータの各々に対応付けられたパラメータ重要度を格納するパラメータ重要度テーブルとを保持し、
前記材料パラメータから次元圧縮により二次元座標として、前記材料データ可視化テーブルに格納する教師データなし次元圧縮部と、
前記材料データ可視化テーブルに格納された二次元座標を材料データ点として、二次元表示し、前記材料データ点の位置変更を受け付ける2D可視化部と、
受け付けた材料データ点の位置変更の情報に基づいて、変更された材料データ点を目的変数とし、変更された材料データ点に対応する前記材料データ可視化テーブルの材料パラメータと材料パラメータに対するパラメータ重要度を説明変数とする回帰分析によって次元圧縮関数と前記次元圧縮関数の材料パラメータに対するパラメータ重要度を算出し、算出したパラメータ重要度により、前記パラメータ重要度テーブルに格納されたパラメータ重要度を更新し、位置変更された材料データ点と位置変更されていない材料データ点の材料パラメータを入力して前記次元圧縮関数より次元圧縮した二次元座標として出力し、出力された二次元座標を前記材料データ可視化テーブルに格納される二次元座標として更新する半教師データあり次元圧縮部と、
更新された前記パラメータ重要度に基づいて、材料特性を予測する材料特性予測提示部とを有することを特徴とする材料特性予測装置。 A material property prediction device for predicting and displaying material properties,
a material data table storing material properties and material parameters associated with the material;
a material data visualization table storing dimensionally compressed information of the material parameters as two-dimensional coordinates;
holding a parameter importance table that stores parameter importance associated with each of the material parameters;
A dimension compression unit without teacher data that stores the material parameters as two-dimensional coordinates by dimension compression in the material data visualization table ;
a 2D visualization unit that two-dimensionally displays the two-dimensional coordinates stored in the material data visualization table as material data points and receives position changes of the material data points;
Based on the received position change information of the material data point, the changed material data point is set as an objective variable, and the material parameter of the material data visualization table corresponding to the changed material data point and the parameter importance for the material parameter are determined. Calculate the dimension compression function and the parameter importance of the dimension compression function with respect to the material parameter by regression analysis using explanatory variables, update the parameter importance stored in the parameter importance table with the calculated parameter importance, and position Input material parameters of changed material data points and unchanged material data points , output as two-dimensional coordinates dimensionally compressed by the dimension compression function, and store the output two-dimensional coordinates in the material data visualization table. a dimensional compression unit with semi-supervised data that updates as stored two-dimensional coordinates;
and a material property prediction presenting unit for predicting material properties based on the updated parameter importance .
前記2D可視化部は、前記予測精度に関する情報を表示することを特徴とする請求項2記載の材料特性予測装置。 The prediction evaluation unit calculates the prediction accuracy of material properties,
3. The material property prediction device according to claim 2 , wherein the 2D visualization unit displays information about the prediction accuracy.
前記2D可視化部は、前記材料データ点の背景に前記背景画像を表示することを特徴とする請求項1記載の材料特性予測装置。 further maintaining a background image for evaluating distances between said material data points ;
2. The material property prediction device according to claim 1, wherein the 2D visualization unit displays the background image on the background of the material data points.
前記材料特性予測装置は、
材料特性とその材料と関連付けられた材料パラメータとを格納する材料データテーブルと、
前記材料パラメータの次元圧縮した情報を二次元座標として格納する材料データ可視化テーブルと、
前記材料パラメータの各々に対応付けられたパラメータ重要度を格納するパラメータ重要度テーブルとを保持し、
前記材料特性予測装置が、前記材料パラメータから次元圧縮により二次元座標として、前記材料データテーブルに格納する教師データなし次元圧縮ステップと、
前記材料データ可視化テーブルに格納された二次元座標を材料データ点として、二次元表示し、前記材料データ点の位置変更を受け付ける2D可視化ステップと、
受け付けた材料データ点の位置変更の情報に基づいて、変更された材料データ点を目的変数とし、変更された材料データ点に対応する前記材料データ可視化テーブルの材料パラメータと材料パラメータに対するパラメータ重要度を説明変数とする回帰分析によって次元圧縮関数と前記次元圧縮関数の材料パラメータに対するパラメータ重要度を算出し、算出したパラメータ重要度により、前記パラメータ重要度テーブルに格納されたパラメータ重要度を更新し、位置変更された材料データ点と位置変更されていない材料データ点の材料パラメータを入力した前記次元圧縮関数より次元圧縮した二次元座標として出力し、出力された二次元座標を前記材料データ可視化テーブルに格納される二次元座標として更新する半教師データあり次元圧縮ステップと、
更新された前記パラメータに基づいて、材料特性を予測し、提示する材料特性予測提示ステップとを有することを特徴とする材料特性予測方法。 A material property prediction method for predicting and displaying material properties by a material property prediction device ,
The material property prediction device is
a material data table storing material properties and material parameters associated with the material;
a material data visualization table storing dimensionally compressed information of the material parameters as two-dimensional coordinates;
holding a parameter importance table that stores parameter importance associated with each of the material parameters;
A dimensional compression step without teacher data in which the material property prediction device stores the material parameters as two-dimensional coordinates by dimensional compression in the material data table;
a 2D visualization step of two-dimensionally displaying the two-dimensional coordinates stored in the material data visualization table as material data points and accepting position changes of the material data points;
Based on the received position change information of the material data point, the changed material data point is set as an objective variable, and the material parameter of the material data visualization table corresponding to the changed material data point and the parameter importance for the material parameter are determined. Calculate the dimension compression function and the parameter importance of the dimension compression function with respect to the material parameter by regression analysis using explanatory variables, update the parameter importance stored in the parameter importance table with the calculated parameter importance, and position The material parameters of the material data points whose positions have been changed and the material data points whose positions have not been changed are output as two-dimensional coordinates dimensionally compressed by the input dimension compression function , and the output two-dimensional coordinates are stored in the material data visualization table. a dimensionality reduction step with semi-supervised data that updates as two-dimensional coordinates of
and a material property prediction presenting step of predicting and presenting material properties based on the updated parameters .
前記2D可視化ステップにおいて、前記予測精度に関する情報を表示することを特徴とする請求項8記載の材料特性予測方法。 In the prediction evaluation step, calculating the prediction accuracy of the material properties,
9. The method of predicting material properties according to claim 8 , wherein information about the prediction accuracy is displayed in the 2D visualization step.
前記材料特性予測提示ステップにおいて、各々の材料特性に関する情報を表示することを特徴とする請求項7記載の材料特性予測方法。 Furthermore, having a material property prediction presentation step,
8. The material property prediction method according to claim 7 , wherein in said material property prediction presentation step, information on each material property is displayed.
前記2D可視化ステップにおいて、前記材料データ点の背景に前記背景画像を表示することを特徴とする請求項7記載の材料特性予測方法。 further maintaining a background image for evaluating distances between said material data points;
8. The method of claim 7, wherein the 2D visualization step displays the background image behind the material data points.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019022666A JP7330712B2 (en) | 2019-02-12 | 2019-02-12 | Material property prediction device and material property prediction method |
PCT/JP2020/002303 WO2020166299A1 (en) | 2019-02-12 | 2020-01-23 | Material characteristics prediction device and material characteristics prediction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019022666A JP7330712B2 (en) | 2019-02-12 | 2019-02-12 | Material property prediction device and material property prediction method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2020128962A JP2020128962A (en) | 2020-08-27 |
JP7330712B2 true JP7330712B2 (en) | 2023-08-22 |
Family
ID=72045024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019022666A Active JP7330712B2 (en) | 2019-02-12 | 2019-02-12 | Material property prediction device and material property prediction method |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7330712B2 (en) |
WO (1) | WO2020166299A1 (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4343578A4 (en) * | 2021-05-18 | 2024-11-13 | Panasonic Holdings Corporation | INFORMATION PROVIDING DEVICE, METHOD AND PROGRAM |
JP7574145B2 (en) | 2021-06-23 | 2024-10-28 | 株式会社Uacj | Data dimension compression method, computer program and data dimension compression device |
JPWO2023008447A1 (en) | 2021-07-28 | 2023-02-02 | ||
WO2023008449A1 (en) | 2021-07-28 | 2023-02-02 | パナソニックIpマネジメント株式会社 | Information display method, information display device, and program |
EP4379731A4 (en) | 2021-07-28 | 2024-11-06 | Panasonic Intellectual Property Management Co., Ltd. | INFORMATION DISPLAY METHOD, INFORMATION DISPLAY DEVICE AND PROGRAM |
WO2023048009A1 (en) * | 2021-09-22 | 2023-03-30 | 株式会社レゾナック | Program, condition search device, and condition search method |
JP7315124B1 (en) * | 2021-10-04 | 2023-07-26 | 株式会社レゾナック | Composition search method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001503546A (en) | 1996-11-04 | 2001-03-13 | 3―ディメンショナル ファーマシューティカルズ インコーポレイテッド | Systems, methods, and computer program products for visualizing and interactively processing and analyzing chemical data |
JP2005506511A (en) | 2000-08-09 | 2005-03-03 | グラクソ グループ リミテッド | Cell-based analysis of high-throughput screening data for drug discovery |
WO2018062398A1 (en) | 2016-09-30 | 2018-04-05 | 株式会社Uacj | Device for predicting aluminum product properties, method for predicting aluminum product properties, control program, and storage medium |
WO2018168580A1 (en) | 2017-03-13 | 2018-09-20 | 日本電気株式会社 | Relation search system, information processing device, method, and program |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3764252A4 (en) * | 2018-03-09 | 2021-12-01 | Showa Denko K.K. | DEVICE FOR PREDICTING THE PHYSICAL PROPERTIES OF POLYMER, STORAGE MEDIUM, AND METHOD OF PREDICTING THE PHYSICAL PROPERTIES OF POLYMER |
JP7122699B2 (en) * | 2018-08-23 | 2022-08-22 | パナソニックIpマネジメント株式会社 | Material information output method, material information output device, material information output system, and program |
-
2019
- 2019-02-12 JP JP2019022666A patent/JP7330712B2/en active Active
-
2020
- 2020-01-23 WO PCT/JP2020/002303 patent/WO2020166299A1/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001503546A (en) | 1996-11-04 | 2001-03-13 | 3―ディメンショナル ファーマシューティカルズ インコーポレイテッド | Systems, methods, and computer program products for visualizing and interactively processing and analyzing chemical data |
JP2005506511A (en) | 2000-08-09 | 2005-03-03 | グラクソ グループ リミテッド | Cell-based analysis of high-throughput screening data for drug discovery |
WO2018062398A1 (en) | 2016-09-30 | 2018-04-05 | 株式会社Uacj | Device for predicting aluminum product properties, method for predicting aluminum product properties, control program, and storage medium |
WO2018168580A1 (en) | 2017-03-13 | 2018-09-20 | 日本電気株式会社 | Relation search system, information processing device, method, and program |
Also Published As
Publication number | Publication date |
---|---|
WO2020166299A1 (en) | 2020-08-20 |
JP2020128962A (en) | 2020-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7330712B2 (en) | Material property prediction device and material property prediction method | |
JP7061536B2 (en) | Optimization device, simulation system and optimization method | |
JPWO2018079225A1 (en) | Automatic prediction system, automatic prediction method, and automatic prediction program | |
JP4230975B2 (en) | Data prediction apparatus, data prediction method, and program | |
US20240054385A1 (en) | Experiment point recommendation device, experiment point recommendation method, and semiconductor device manufacturing device | |
US20210097884A1 (en) | Work support apparatus, work support system, and work support method | |
JP2022090206A (en) | Optimization support device and optimization support method | |
CN111626789A (en) | House price prediction method, device, equipment and storage medium | |
US20250037381A1 (en) | Body model processing methods and apparatuses, electronic devices and storage media | |
Dergachev et al. | Evaluating the ergonomics of online store user interfaces based on visual analytics | |
JP7617820B2 (en) | Constraint condition optimization device and constraint condition optimization method | |
WO2023175921A1 (en) | Model analysis device, model analysis method, and recording medium | |
CN111340276B (en) | Method and system for generating prediction data | |
JP2020024542A (en) | Data analysis device and data analysis method | |
JP2024059400A (en) | Analysis device, analysis method, and analysis program | |
Packham et al. | Visualisation methods for supporting the exploration of high dimensional problem spaces in engineering design | |
JP4922644B2 (en) | Time series analysis program, time series analysis system, and time series analysis apparatus used therefor | |
JP7095744B2 (en) | Prediction status visualization device, forecast status visualization method and forecast status visualization program | |
JP7270454B2 (en) | Design support system, design support method and design support program | |
JP2006259940A (en) | Statistical information presentation device | |
JP2022014618A (en) | Prediction device and prediction method | |
JP2021179668A (en) | Data analysis system, data analysis method, and data analysis program | |
US9275484B2 (en) | Goodness of fit based on error calculation and fit type | |
D’Angelo et al. | Minimum contrast for the first-order intensity estimation of spatial and spatio-temporal point processes | |
JP7453895B2 (en) | Search condition presentation device, search condition presentation method, and search condition presentation program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220120 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230404 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230602 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230718 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230809 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7330712 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |