JP5316433B2 - 最適化処理プログラム、方法及び装置 - Google Patents
最適化処理プログラム、方法及び装置 Download PDFInfo
- Publication number
- JP5316433B2 JP5316433B2 JP2010014062A JP2010014062A JP5316433B2 JP 5316433 B2 JP5316433 B2 JP 5316433B2 JP 2010014062 A JP2010014062 A JP 2010014062A JP 2010014062 A JP2010014062 A JP 2010014062A JP 5316433 B2 JP5316433 B2 JP 5316433B2
- Authority
- JP
- Japan
- Prior art keywords
- design parameter
- value
- storage unit
- model formula
- limit value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 97
- 238000005457 optimization Methods 0.000 title claims description 46
- 238000013461 design Methods 0.000 claims abstract description 288
- 230000014509 gene expression Effects 0.000 claims abstract description 106
- 238000011156 evaluation Methods 0.000 claims abstract description 42
- 238000012545 processing Methods 0.000 claims description 146
- 238000003860 storage Methods 0.000 claims description 84
- 238000013500 data storage Methods 0.000 claims description 67
- 238000004088 simulation Methods 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 16
- 238000012986 modification Methods 0.000 claims description 7
- 230000004048 modification Effects 0.000 claims description 7
- 238000003672 processing method Methods 0.000 claims description 7
- 230000001131 transforming effect Effects 0.000 claims description 5
- 230000008030 elimination Effects 0.000 abstract description 4
- 238000003379 elimination reaction Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 42
- 230000006870 function Effects 0.000 description 37
- 238000010586 diagram Methods 0.000 description 29
- 238000004364 calculation method Methods 0.000 description 27
- 238000005516 engineering process Methods 0.000 description 11
- 230000035945 sensitivity Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
雑誌FUJITSU2009-9月号,穴井 宏和, 金児 純司, 屋並 仁史, 岩根 秀直,「数式処理を用いた設計技術」(http://img.jp.fujitsu.com/downloads/jp/jmag/vol60-5/paper24.pdfから2009年10月取得可能)
Mats Jirstrand: Cylindrical Algebraic Decomposition - an Introduction,1995−10−18(http://www.control.isy.liu.se/research/reports/1995/1807.ps.Zから2009年10月取得可能)
∃x1∃x2∃x3∃x4∃x5∃x6∃x7(0≦x1≦1∧0≦x2≦1∧0≦x3≦1∧0≦x4≦1∧0≦x5≦1∧0≦x6≦1∧0≦x7≦1∧z≦F11∧z≦F12)
1274223392000000000*z−10521725891745634157≦0
1)x7を消去したQE問題を作成し、x7の値を求める。
∃x1∃x2∃x3∃x4∃x5∃x6∃z(0≦x1≦1∧0≦x2≦1∧0≦x3≦1∧0≦x4≦1∧0≦x5≦1∧0≦x6≦1∧0≦x7≦1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
2)x7の結果を条件に反映し、x6を消去したQE問題を作成し、x6の値を求める。
∃x1∃x2∃x3∃x4∃x5∃x7∃z(0≦x1≦1∧0≦x2≦1∧0≦x3≦1∧0≦x4≦1∧0≦x5≦1∧0≦x6≦1∧x7=1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
3)x6,x7の結果を条件に反映し、x5を消去したQE問題を作成し、x5の値を求める。
∃x1∃x2∃x3∃x4∃x6∃x7∃z(0≦x1≦1∧0≦x2≦1∧0≦x3≦1∧0≦x4≦1∧0≦x5≦1∧x6=1∧x7=1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
4)x5,x6,x7の結果を条件に反映し、x4を消去したQE問題を作成し、x4の値を求める。
∃x1∃x2∃x3∃x5∃x6∃x7∃z(0≦x1≦1∧0≦x2≦1∧0≦x3≦1∧0≦x4≦1∧x5=1∧x6=1∧x7=1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
5)x4,x5,x6,x7の結果を条件に反映し、x3を消去したQE問題を作成し、x3の値を求める。
∃x1∃x2∃x4∃x5∃x6∃x7∃z(0≦x1≦1∧0≦x2≦1∧0≦x3≦1∧x4=1∧x5=1∧x6=1∧x7=1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
6)x3,x4,x5,x6,x7の結果を条件に反映し、x2を消去したQE問題を作成し、x2の値を求める。
∃x1∃x3∃x4∃x5∃x6∃x7∃z(0≦x1≦1∧0≦x2≦1∧x3=1∧x4=1∧x5=1∧x6=1∧x7=1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
7)x2,x3,x4,x5,x6,x7の結果を条件に反映し、x1を消去したQE問題を作成し、x1の値を求める。
∃x2∃x3∃x4∃x5∃x6∃x7∃z(0≦x1≦1∧x2=1∧x3=1∧x4=1∧x5=1∧x6=1∧x7=1∧z≦F11∧z≦F12∧1274223392000000000*z−10521725891745634157=0)
x1=7013720457/14016457312=0.5003918,
x2=1,
x3=1,
x4=1,
x5=1,
x6=1,
x7=1,
z=10521725891745634157/1274223392000000000
F21b=8.500756371−0.8407133034*x1+0.1922736364*x1 2
F22b=7.953391257+1.351922165*x1+0.3143480967*x1 2
∃x1(0≦x1≦1∧z≦F21b∧z≦F22b)
238618787270928516302097379198561373969109372528*z2−3948898513981542665325272158255992498641993331585*z+16337501656369998932875768265187359652159276835066≦0 or
61093111211285133307320208872464185567392*z−505514462704138254311579799314432635828595<0
最終的には、変数zの実行可能領域は、z≦8.29074…となる。従って、変数zの最大値は、8.29074…となる。なお、本願で小数点以下の数字を「…」を含む形で表した場合、実際にはルート√を含む形で出力されるが、表記の都合上近似値を用いていることを表している。
∃x1∃z(0≦x1≦1∧z≦F21b∧z≦F22b∧238618787270928516302097379198561373969109372528*z2−3948898513981542665325272158255992498641993331585*z+16337501656369998932875768265187359652159276835066=0 or
61093111211285133307320208872464185567392*z−505514462704138254311579799314432635828595>0)
なお、この式においても表記の都合上近似値が用いられている。
147369062865673278681489897354299172230888310212951405753448549687*x1 4−1267585231880821011490655336656612299953463686527769426803057971300*x1 3+3026832164187229605677018935740423336948896624044750813252776489960*x1 2−1294814954233025589235990288953847635054037925349792811665028977000*x1+153374672525527731144963940745339046776486276154608383375827336700≧0 and
13783646935592814416188923068045742305401334162815422297006892092*x1 4−120537537699802114369719508264578689747806603680459929386166019204*x1 3+295964917461532064463036717507200901977755721564516665233454151547*x1 2−141846512456883797861995180769761248346815031482078133627654338810*x1+18989383858381045985849754233804905380729776996135581477310597575≧0 and
1835405414309906452852689890486911373094821422542*x1 2−8025279897097359340954146906679274179148864319177*x1+2159873695219589023157307374252815016052936539755≧0 and
545582478146275943640962608991635504887081954057*x1 2−2346395772026885507191646225202102006739498467150*x1+557302446970863503723288521839833928729371673530≦0 and
x1−1≦0 and
x1≧0
x1−1<0 and
184630436769134077176252*x1 2−799072005181775017876755*x1+199478730104451626416796=0
F21c=-2.464051045*x7+10.92366518−0.7848741460*x1+0.1922736364*x1 2+0.4114223556×10-1*x7 2−0.5583915742×10-1*x7*x1
F22c=4.739436358*x7+3.052324427+1.279031039*x1−0.3143480967*x1 2+0.1616304724*x7 2+0.7289112590×10-1*x7*x1
∃x1∃x7(0≦x1≦1∧0≦x7≦1∧z≦F21c∧z≦F22c)
76521692628371604200094665888089031213570348775719997908381149434130894497838804645345122280727232499000550400*z2−1266356274402535285733613562992528809428335713014481227199628488088920626515945478646176290112631818545299010496*z+5239207254170329463799573892494718766285329013399382293501296580397159499617737745543801847759123880991116365945≦0 and
463443635747785967425472797218181077323632916075879549238739474499290774293298730400*z−3834760678983826041395849039351295208012618621640274564860057389386590727945878723723<0
∃x1∃x7∃z(0≦x1≦1∧0≦x7≦1∧z≦F21c∧z≦F22c∧z=8.290742434)
x1=0.2659847199,
x7=0.9999999988
以下、QEツール300により実施される実行可能領域の算出処理について具体例を用いて説明する。
(1−1)関数F及びGを以下のように設定する。
F(X,Y)=X2+Y2−1
G(X,Y)=X3−Y2
(1−2)判別式dircriminant(F,Y)=−4x2+4=0を計算する。
X=−1,1
(1−3)判別式dircriminant(G,Y)=4x3=0を計算する。
X=0
(1−4)終結式resultant(F,G,Y)=(x3+x2−1)=0を計算する。
この式を満たすXの値をAとする。
X=A
(2−1)計算手順1で計算したXの値を整理する。
X={−1,0,A,1}
(2−2)Xの最小値より小さい値と、Xの計算値と計算値の間の値、Xの最大値より大きい値を追加する。
X={X1,X2=−1,X3,X4=0,X5,X6=A,X7,X8=1,X9}
Xの各々の値の場合について、F(X,Y)=0又はG(X,Y)=0となるYを計算する。
X=X1・・・該当無し
X=X2・・・Y={Y21=0}
X=X3・・・Y={Y31,Y32}
X=X4・・・Y={Y41=−1,Y42=0,Y43=1}
X=X5・・・Y={Y51,Y52,Y53,Y54}
X=X6・・・Y={Y61=−A3/2,Y62=A3/2}
X=X7・・・Y={Y71,Y72,Y73,Y74}
X=X8・・・Y={Y81=−1,Y82=0,Y83=1}
X=X9・・・Y={Y91,Y92}
Yの最小値より小さい値、Yの計算値と計算値の間の値、Yの最大値より大きい値を追加する。該当無しの場合には0を設定する。
X=X1・・・Y={YY11=0}
X=X2・・・Y={YY21,YY22=Y21,YY23}
X=X3・・・Y={YY31,YY32=Y31,YY33,YY34=Y32,YY35}
X=X4・・・Y={YY41,YY42=Y41,YY43,YY44=Y42,YY45,YY46=Y43,YY47}
X=X5・・・Y={YY51,YY52=Y51,YY53,YY54=Y52,YY55,YY56=Y53,YY57,YY58=Y54,YY59}
X=X6・・・Y={YY61,YY62=Y61,YY63,YY64=Y62,YY65}
X=X7・・・Y={YY71,YY72=Y71,YY73,YY74=Y72,YY75,YY76=Y73,YY77,YY78=Y74,YY79}
X=X8・・・Y={YY81,YY82=Y81,YY83,YY84=Y82,YY85,YY86=Y83,YY87}
X=X9・・・Y={YY91,YY92=Y91,YY93,YY94=Y92,YY95}
(5−1)計算手順4で算出された、XとYの組み合わせについて、F(X,Y)とG(X,Y)の符号を計算する。
(X,Y)=(X1,YY11)−>(F,G)=(+,−)
というような計算を全ての組み合わせについて実施する。
図16に示すように、XY平面において、F(X,Y)=0は原点を中心とした半径1の円であり、G(X,Y)=0はGで表す曲線である。X1は−1より小さい値であり、X=X1は点線で表す直線となる。YY11=0であるので、X=X1とY=0の交点が(X1,YY11)である。F及びGの定義から、(F,G)=(+,−)となる。
(X,Y)=(X2,YY23)−>(F,G)=(+,−)
(X,Y)=(X2,YY22)−>(F,G)=(0,−)
(X,Y)=(X2,YY21)−>(F,G)=(+,−)
(X,Y)=(X3,YY35)−>(F,G)=(+,−)
(X,Y)=(X3,YY34)−>(F,G)=(0,−)
(X,Y)=(X3,YY33)−>(F,G)=(−,−)
(X,Y)=(X3,YY32)−>(F,G)=(0,−)
(X,Y)=(X3,YY31)−>(F,G)=(+,−)
(X,Y)=(X4,YY47)−>(F,G)=(+,−)
(X,Y)=(X4,YY46)−>(F,G)=(0,−)
(X,Y)=(X4,YY45)−>(F,G)=(−,−)
(X,Y)=(X4,YY44)−>(F,G)=(−,0)
(X,Y)=(X4,YY43)−>(F,G)=(−,−)
(X,Y)=(X4,YY42)−>(F,G)=(0,−)
(X,Y)=(X4,YY41)−>(F,G)=(+,−)
(X,Y)=(X5,YY59)−>(F,G)=(+,−)
(X,Y)=(X5,YY58)−>(F,G)=(0,−)
(X,Y)=(X5,YY57)−>(F,G)=(−,−)
(X,Y)=(X5,YY56)−>(F,G)=(−,0)
(X,Y)=(X5,YY55)−>(F,G)=(−,+)
(X,Y)=(X5,YY54)−>(F,G)=(−,0)
(X,Y)=(X5,YY53)−>(F,G)=(−,−)
(X,Y)=(X5,YY52)−>(F,G)=(0,−)
(X,Y)=(X5,YY51)−>(F,G)=(+,−)
(X,Y)=(X6,YY65)−>(F,G)=(+,−)
(X,Y)=(X6,YY64)−>(F,G)=(0,0)
(X,Y)=(X6,YY63)−>(F,G)=(−,+)
(X,Y)=(X6,YY62)−>(F,G)=(0,0)
(X,Y)=(X6,YY61)−>(F,G)=(+,−)
(X,Y)=(X7,YY79)−>(F,G)=(+,−)
(X,Y)=(X7,YY78)−>(F,G)=(+,0)
(X,Y)=(X7,YY77)−>(F,G)=(+,+)
(X,Y)=(X7,YY76)−>(F,G)=(0,+)
(X,Y)=(X7,YY75)−>(F,G)=(−,+)
(X,Y)=(X7,YY74)−>(F,G)=(0,+)
(X,Y)=(X7,YY73)−>(F,G)=(+,+)
(X,Y)=(X7,YY72)−>(F,G)=(+,0)
(X,Y)=(X7,YY71)−>(F,G)=(+,−)
(X,Y)=(X8,YY87)−>(F,G)=(+,−)
(X,Y)=(X8,YY86)−>(F,G)=(+,0)
(X,Y)=(X8,YY85)−>(F,G)=(+,+)
(X,Y)=(X8,YY84)−>(F,G)=(0,+)
(X,Y)=(X8,YY83)−>(F,G)=(+,+)
(X,Y)=(X8,YY82)−>(F,G)=(+,0)
(X,Y)=(X8,YY81)−>(F,G)=(+,−)
(X,Y)=(X9,YY95)−>(F,G)=(+,−)
(X,Y)=(X9,YY94)−>(F,G)=(+,0)
(X,Y)=(X9,YY93)−>(F,G)=(+,+)
(X,Y)=(X9,YY92)−>(F,G)=(+,0)
(X,Y)=(X9,YY91)−>(F,G)=(+,−)
これは、制約条件としてG(X,Y)=0であり、F(X,Y)=Z<0であるためである。
(X,Y)=(X4,YY44),(X5,YY54),(X5,YY56)
(6−1)(X,Y)=(X4,YY44)を満足する条件を計算する。
図19及び上で述べた説明から分かるように、(X,Y)=(X4,YY44)=(0,0)である。
(6−2)(X,Y)=(X5,YY54)を満足する条件を計算する。
図20及び上で述べた説明から分かるように、X5は0<X<Aを満たす。また、Y<YY55=0である。さらに、G(X,Y)=0上の点であるから、X3−Y2=0である。
(6−3)(X,Y)=(X5,YY56)を満足する条件を計算する。
図20及び上で述べた説明から分かるように、X5は0<X<Aを満たす。また、Y>YY55=0である。さらに、G(X,Y)=0上の点であるから、X3−Y2=0である。
(X,Y)=(X4,YY44)、(X,Y)=(X5,YY54)又は(X,Y)=(X5,YY56)を満足し且つ計算手順6で算出された条件を整理する。
この例では、(X,Y)=(X4,Y44)があるので(X,Y)=(0,0)が含まれる。従って、X3−Y2=0且つ0≦X<Aとなる。すなわち、図25に示すように、G(X,Y)=0曲線上の一部分が実行可能領域となる。
予め定められた設計パラメータの値を複数セット生成して各セットについてシミュレーションを実行させることにより、予め定められた評価指標の値を取得し、前記評価指標と前記設計パラメータとの関係を表し且つ第1の次数についての第1のモデル式と、前記評価指標と前記設計パラメータとの関係を表し且つ前記第1の次数より高次の第2のモデル式とを生成し、モデル式データ格納部に格納するステップと、
前記モデル式データ格納部に格納されている前記第1のモデル式と予め定められた制約条件とから限定子除去法により前記第1のモデル式の最適解を実現する前記設計パラメータの値を算出させ、設計パラメータ値格納部に格納するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータを特定し、前記モデル式データ格納部に格納されている前記第2のモデル式に対して、特定された前記設計パラメータの値として前記上限値又は下限値を代入して前記第2のモデル式を変形する変形ステップと、
変形後の第2のモデル式と前記制約条件の該当部分とから前記限定子除去法により前記変形後の第2のモデル式の最適解を実現する、前記変形後の第2のモデル式内の設計パラメータの値を算出させ、前記設計パラメータ値格納部に格納するステップと、
を、コンピュータにより実行させるための最適化処理プログラム。
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータの中から、前記モデル式データ格納部に格納されている前記第1のモデル式に対する影響度の条件を満たす設計パラメータを抽出するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータのうち、抽出された前記設計パラメータ以外の設計パラメータの値として前記上限値又は下限値を、前記モデル式データ格納部に格納されている前記第2のモデル式に代入して前記第2のモデル式の第2の変形式を生成するステップと、
前記第2のモデル式の第2の変形式と前記制約条件の該当部分とから前記限定子除去法により前記第2の変形式の最適解を実現する、前記第2の変形式内の設計パラメータの値を算出し、前記設計パラメータ値格納部に格納するステップと、
をさらに前記コンピュータに実行させるための付記1記載の最適化処理プログラム。
予め定められた設計パラメータの値を複数セット生成して各セットについてシミュレーションを実行させることにより、予め定められた評価指標の値を取得し、前記評価指標と前記設計パラメータとの関係を表し且つ第1の次数についての第1のモデル式と、前記評価指標と前記設計パラメータとの関係を表し且つ前記第1の次数より高次の第2のモデル式とを生成し、モデル式データ格納部に格納するステップと、
前記モデル式データ格納部に格納されている前記第1のモデル式と予め定められた制約条件とから限定子除去法により前記第1のモデル式の最適解を実現する前記設計パラメータの値を算出させ、設計パラメータ値格納部に格納するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータを特定し、前記モデル式データ格納部に格納されている前記第2のモデル式に対して、特定された前記設計パラメータの値として前記上限値又は下限値を代入して前記第2のモデル式を変形する変形ステップと、
変形後の第2のモデル式と前記制約条件の該当部分とから前記限定子除去法により前記変形後の第2のモデル式の最適解を実現する、前記変形後の第2のモデル式内の設計パラメータの値を算出させ、前記設計パラメータ値格納部に格納するステップと、
を含み、コンピュータに実行される最適化処理方法。
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータの中から、前記モデル式データ格納部に格納されている前記第1のモデル式に対する影響度の条件を満たす設計パラメータを抽出するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータのうち、抽出された前記設計パラメータ以外の設計パラメータの値として前記上限値又は下限値を、前記モデル式データ格納部に格納されている前記第2のモデル式に代入して前記第2のモデル式の第2の変形式を生成するステップと、
前記第2のモデル式の第2の変形式と前記制約条件の該当部分とから前記限定子除去法により前記第2の変形式の最適解を実現する、前記第2の変形式内の設計パラメータの値を算出し、前記設計パラメータ値格納部に格納するステップと、
をさらに含む付記3記載の最適化処理方法。
予め定められた設計パラメータの値を複数セット生成して各セットについてシミュレーションを実行させることにより、予め定められた評価指標の値を取得し、前記評価指標と前記設計パラメータとの関係を表し且つ第1の次数についての第1のモデル式と、前記評価指標と前記設計パラメータとの関係を表し且つ前記第1の次数より高次の第2のモデル式とを生成し、モデル式データ格納部に格納するモデル式生成部と、
前記モデル式データ格納部に格納されている前記第1のモデル式と予め定められた制約条件とから限定子除去法により前記第1のモデル式の最適解を実現する前記設計パラメータの値を算出させ、設計パラメータ値格納部に格納する第1設計パラメータ値生成部と、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータを特定し、前記モデル式データ格納部に格納されている前記第2のモデル式に対して、特定された前記設計パラメータの値として前記上限値又は下限値を代入して前記第2のモデル式を変形し、変形後の第2のモデル式と前記制約条件の該当部分とから前記限定子除去法により前記変形後の第2のモデル式の最適解を実現する、前記変形後の第2のモデル式内の設計パラメータの値を算出させ、前記設計パラメータ値格納部に格納する第2設計パラメータ値生成部と、
を有する最適化処理装置。
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータの中から、前記モデル式データ格納部に格納されている前記第1のモデル式に対する影響度の条件を満たす設計パラメータを抽出し、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータのうち、抽出された前記設計パラメータ以外の設計パラメータの値として前記上限値又は下限値を、前記モデル式データ格納部に格納されている前記第2のモデル式に代入して前記第2のモデル式の第2の変形式を生成し、
前記第2のモデル式の第2の変形式と前記制約条件の該当部分とから前記限定子除去法により前記第2の変形式の最適解を実現する、前記第2の変形式内の設計パラメータの値を算出し、前記設計パラメータ値格納部に格納する
第3設計パラメータ値生成部
をさらに有する付記5記載の最適化処理装置。
3 モデル式生成部 4 モデル式データ格納部
5 第1QE処理部 6 第2QE処理部 7 第3QE処理部
8 処理結果格納部 9 出力部
Claims (6)
- 予め定められた設計パラメータの値を複数セット生成して各セットについてシミュレーションを実行させることにより、予め定められた評価指標の値を取得し、前記評価指標と前記設計パラメータとの関係を表し且つ第1の次数についての第1のモデル式と、前記評価指標と前記設計パラメータとの関係を表し且つ前記第1の次数より高次の第2のモデル式とを生成し、モデル式データ格納部に格納するステップと、
前記モデル式データ格納部に格納されている前記第1のモデル式と予め定められた制約条件とから限定子除去法により前記第1のモデル式の最適解を実現する前記設計パラメータの値を算出させ、設計パラメータ値格納部に格納するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータを特定し、前記モデル式データ格納部に格納されている前記第2のモデル式に対して、特定された前記設計パラメータの値として前記上限値又は下限値を代入して前記第2のモデル式を変形する変形ステップと、
変形後の第2のモデル式と前記制約条件の該当部分とから前記限定子除去法により前記変形後の第2のモデル式の最適解を実現する、前記変形後の第2のモデル式内の設計パラメータの値を算出させ、前記設計パラメータ値格納部に格納するステップと、
を、コンピュータにより実行させるための最適化処理プログラム。 - 前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータの中から、前記モデル式データ格納部に格納されている前記第1のモデル式に対する影響度の条件を満たす設計パラメータを抽出するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータのうち、抽出された前記設計パラメータ以外の設計パラメータの値として前記上限値又は下限値を、前記モデル式データ格納部に格納されている前記第2のモデル式に代入して前記第2のモデル式の第2の変形式を生成するステップと、
前記第2のモデル式の第2の変形式と前記制約条件の該当部分とから前記限定子除去法により前記第2の変形式の最適解を実現する、前記第2の変形式内の設計パラメータの値を算出し、前記設計パラメータ値格納部に格納するステップと、
をさらに前記コンピュータに実行させるための請求項1記載の最適化処理プログラム。 - 予め定められた設計パラメータの値を複数セット生成して各セットについてシミュレーションを実行させることにより、予め定められた評価指標の値を取得し、前記評価指標と前記設計パラメータとの関係を表し且つ第1の次数についての第1のモデル式と、前記評価指標と前記設計パラメータとの関係を表し且つ前記第1の次数より高次の第2のモデル式とを生成し、モデル式データ格納部に格納するステップと、
前記モデル式データ格納部に格納されている前記第1のモデル式と予め定められた制約条件とから限定子除去法により前記第1のモデル式の最適解を実現する前記設計パラメータの値を算出させ、設計パラメータ値格納部に格納するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータを特定し、前記モデル式データ格納部に格納されている前記第2のモデル式に対して、特定された前記設計パラメータの値として前記上限値又は下限値を代入して前記第2のモデル式を変形する変形ステップと、
変形後の第2のモデル式と前記制約条件の該当部分とから前記限定子除去法により前記変形後の第2のモデル式の最適解を実現する、前記変形後の第2のモデル式内の設計パラメータの値を算出させ、前記設計パラメータ値格納部に格納するステップと、
を含み、コンピュータに実行される最適化処理方法。 - 前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータの中から、前記モデル式データ格納部に格納されている前記第1のモデル式に対する影響度の条件を満たす設計パラメータを抽出するステップと、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータのうち、抽出された前記設計パラメータ以外の設計パラメータの値として前記上限値又は下限値を、前記モデル式データ格納部に格納されている前記第2のモデル式に代入して前記第2のモデル式の第2の変形式を生成するステップと、
前記第2のモデル式の第2の変形式と前記制約条件の該当部分とから前記限定子除去法により前記第2の変形式の最適解を実現する、前記第2の変形式内の設計パラメータの値を算出し、前記設計パラメータ値格納部に格納するステップと、
をさらに含む請求項3記載の最適化処理方法。 - 予め定められた設計パラメータの値を複数セット生成して各セットについてシミュレーションを実行させることにより、予め定められた評価指標の値を取得し、前記評価指標と前記設計パラメータとの関係を表し且つ第1の次数についての第1のモデル式と、前記評価指標と前記設計パラメータとの関係を表し且つ前記第1の次数より高次の第2のモデル式とを生成し、モデル式データ格納部に格納するモデル式生成部と、
前記モデル式データ格納部に格納されている前記第1のモデル式と予め定められた制約条件とから限定子除去法により前記第1のモデル式の最適解を実現する前記設計パラメータの値を算出させ、設計パラメータ値格納部に格納する第1設計パラメータ値生成部と、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータを特定し、前記モデル式データ格納部に格納されている前記第2のモデル式に対して、特定された前記設計パラメータの値として前記上限値又は下限値を代入して前記第2のモデル式を変形し、変形後の第2のモデル式と前記制約条件の該当部分とから前記限定子除去法により前記変形後の第2のモデル式の最適解を実現する、前記変形後の第2のモデル式内の設計パラメータの値を算出させ、前記設計パラメータ値格納部に格納する第2設計パラメータ値生成部と、
を有する最適化処理装置。 - 前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータの中から、前記モデル式データ格納部に格納されている前記第1のモデル式に対する影響度の条件を満たす設計パラメータを抽出し、
前記設計パラメータ値格納部に格納されている前記設計パラメータの値が前記制約条件に含まれる上限値又は下限値となった前記設計パラメータのうち、抽出された前記設計パラメータ以外の設計パラメータの値として前記上限値又は下限値を、前記モデル式データ格納部に格納されている前記第2のモデル式に代入して前記第2のモデル式の第2の変形式を生成し、
前記第2のモデル式の第2の変形式と前記制約条件の該当部分とから前記限定子除去法により前記第2の変形式の最適解を実現する、前記第2の変形式内の設計パラメータの値を算出し、前記設計パラメータ値格納部に格納する
第3設計パラメータ値生成部
をさらに有する請求項5記載の最適化処理装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010014062A JP5316433B2 (ja) | 2010-01-26 | 2010-01-26 | 最適化処理プログラム、方法及び装置 |
US13/009,041 US8498844B2 (en) | 2010-01-26 | 2011-01-19 | Optimization processing method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010014062A JP5316433B2 (ja) | 2010-01-26 | 2010-01-26 | 最適化処理プログラム、方法及び装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011154439A JP2011154439A (ja) | 2011-08-11 |
JP5316433B2 true JP5316433B2 (ja) | 2013-10-16 |
Family
ID=44309614
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010014062A Active JP5316433B2 (ja) | 2010-01-26 | 2010-01-26 | 最適化処理プログラム、方法及び装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8498844B2 (ja) |
JP (1) | JP5316433B2 (ja) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5402351B2 (ja) * | 2009-07-24 | 2014-01-29 | 富士通株式会社 | 多目的最適化設計支援装置、方法、及びプログラム |
JP5418409B2 (ja) * | 2010-06-01 | 2014-02-19 | 富士通株式会社 | モデル式生成方法、装置及びプログラム |
JP5477242B2 (ja) * | 2010-09-22 | 2014-04-23 | 富士通株式会社 | 最適化処理プログラム、方法及び装置 |
JPWO2013073125A1 (ja) * | 2011-11-15 | 2015-04-02 | 日本電気株式会社 | 端末、通信システムおよび操作共有方法 |
JP5765266B2 (ja) * | 2012-02-24 | 2015-08-19 | 富士通株式会社 | 性能評価方法、情報処理装置およびプログラム |
CN109918718B (zh) * | 2019-01-25 | 2023-01-10 | 沈阳无距科技有限公司 | 无人机飞行质量评估方法、装置、存储介质及电子设备 |
US12001779B2 (en) * | 2022-01-14 | 2024-06-04 | Huawei Cloud Computing Technologies Co., Ltd. | Method and system for automatically formulating an optimization problem using machine learning |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0754530B2 (ja) * | 1987-10-14 | 1995-06-07 | 沖電気工業株式会社 | アナログ回路特性合せ込み方法 |
JP3484344B2 (ja) * | 1998-05-13 | 2004-01-06 | 富士通株式会社 | 制御系解析・設計装置および制御系解析・設計処理方法 |
JP2003016114A (ja) * | 2001-06-29 | 2003-01-17 | Hitachi Ltd | パラメータ最適化装置 |
JP4101046B2 (ja) * | 2002-12-16 | 2008-06-11 | キヤノン株式会社 | 最適設計方法 |
JP2005129015A (ja) * | 2003-09-30 | 2005-05-19 | Fujitsu Ltd | シミュレーションにおけるモデルパラメータ決定プログラム、および決定装置 |
US20050257178A1 (en) | 2004-05-14 | 2005-11-17 | Daems Walter Pol M | Method and apparatus for designing electronic circuits |
US7516423B2 (en) | 2004-07-13 | 2009-04-07 | Kimotion Technologies | Method and apparatus for designing electronic circuits using optimization |
JP2006106817A (ja) * | 2004-09-30 | 2006-04-20 | Fujitsu Ltd | モデルパラメータ計算精度判定プログラム、およびそのプログラムを格納した記憶媒体 |
JP4802789B2 (ja) * | 2006-03-20 | 2011-10-26 | トヨタ自動車株式会社 | 設計値最適化方法及び設計値最適化システム |
JP2007272478A (ja) * | 2006-03-30 | 2007-10-18 | Fujitsu Ltd | モデルパラメータ決定プログラム、および決定装置 |
JP5018487B2 (ja) * | 2008-01-14 | 2012-09-05 | 富士通株式会社 | 製造バラつきを考慮した多目的最適化設計支援装置、方法、及びプログラム |
JP5458814B2 (ja) * | 2009-11-11 | 2014-04-02 | 富士通株式会社 | 数値処理プログラム、方法及び装置 |
-
2010
- 2010-01-26 JP JP2010014062A patent/JP5316433B2/ja active Active
-
2011
- 2011-01-19 US US13/009,041 patent/US8498844B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US8498844B2 (en) | 2013-07-30 |
JP2011154439A (ja) | 2011-08-11 |
US20110184706A1 (en) | 2011-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5316433B2 (ja) | 最適化処理プログラム、方法及び装置 | |
Papalambros et al. | Principles of optimal design: modeling and computation | |
US20150149128A1 (en) | Systems and methods for analyzing model parameters of electrical power systems using trajectory sensitivities | |
Mackman et al. | Investigation of an adaptive sampling method for data interpolation using radial basis functions | |
Parnianifard et al. | Recent developments in metamodel based robust black-box simulation optimization: An overview | |
US20170338802A1 (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
JP5402351B2 (ja) | 多目的最適化設計支援装置、方法、及びプログラム | |
WO2017131953A1 (en) | Design optimization for a wind turbine tower | |
JP2019537079A (ja) | 大規模再生可能エネルギーのデータについて確率モデルを構築する方法 | |
CN101852871A (zh) | 基于经验模态分解和数值集合预报的短期气候预测方法 | |
Varga et al. | Symbolic and numerical software tools for LFT-based low order uncertainty modeling | |
CN114564787B (zh) | 用于目标相关翼型设计的贝叶斯优化方法、装置及存储介质 | |
Dong et al. | Optimal layout of internal stiffeners for three-dimensional box structures based on natural branching phenomena | |
CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
JP5375429B2 (ja) | 設計支援装置、方法、及びプログラム | |
Sueyoshi et al. | Pitfalls and remedies in DEA applications: how to handle an occurrence of zero in multipliers by strong complementary slackness conditions | |
Yamazaki | Efficient multi-objective shape optimization using proper orthogonal decomposition with variable fidelity concept | |
Liu et al. | A sequential sampling generation method for multi-fidelity model based on Voronoi region and sample density | |
CN112434451B (zh) | 一种基于分块并行计算的有限元分析方法 | |
JP5458814B2 (ja) | 数値処理プログラム、方法及び装置 | |
JP2013235512A (ja) | 数理計画問題を解くための装置、プログラムおよび方法 | |
Johnson et al. | A Coupled Karhunen–Loève and anisotropic sparse grid interpolation method for the probabilistic load flow problem | |
Su et al. | Probabilistic load flow analysis based on sparse polynomial chaotic expansion | |
US8483999B2 (en) | Method and system for simplifying models | |
EP3989085A1 (en) | Optimization device, optimization method, and optimization program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20121005 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130531 |
|
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: 20130611 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130624 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5316433 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |