s03.rtf

机械优化设计中的惩罚函数法

                          常用优化方法——惩罚函数法 
                          ^^^^^^^^^^^^^^^^^^^^^^^^^^ 
                                                                               
一、初始数据
===============================================================================
   设计变量个数     N = 2
  -----------------------------------------------------------------------------
   不等式约束个数  KG = 3       等式约束个数      KH = 1
  -----------------------------------------------------------------------------
   惩罚因子         R = 1       惩罚因子降低系数   C = 0.2
  -----------------------------------------------------------------------------
   初始步长        T0 = 0.01       收敛精度         EPS = 0.0001
  -----------------------------------------------------------------------------
   无约束优化方法: POWELL法
  -----------------------------------------------------------------------------
   设计变量初始点 X0:
			X[1]=2
			X[2]=3
  -----------------------------------------------------------------------------
   设计变量下界 BL:
			BL[1]=0
			BL[2]=0
  -----------------------------------------------------------------------------
   设计变量上界 BU:
			BU[1]=6
			BU[2]=8
  -----------------------------------------------------------------------------
   初始点目标函数值 F(X0)= -13
  -----------------------------------------------------------------------------
   初始点处的不等约束函数值 G(X0):
				GX[1]= -3.000000E+00
				GX[2]= -2.000000E+00
				GX[3]= -3.000000E+00
  -----------------------------------------------------------------------------
   初始点处的等式约束函数值 H(X0):
				HX[1]= -1.200000E+01
-------------------------------------------------------------------------------
                                                              
二、计算过程__数据
===============================================================================
   IRC = 0	   R = 1.000000E+00	   PEN = 132.166666666667
  -----------------------------------------------------------------------------
   IRC = 1	   R = 2.000000E-01	   PEN = -5.60401469986731
  -----------------------------------------------------------------------------
   IRC = 2	   R = 4.000000E-02	   PEN = -28.4527570744455
  -----------------------------------------------------------------------------
   IRC = 3	   R = 8.000000E-03	   PEN = -28.6160466549713
  -----------------------------------------------------------------------------
   IRC = 4	   R = 1.600000E-03	   PEN = -31.309856602631
  -----------------------------------------------------------------------------
   IRC = 5	   R = 3.200000E-04	   PEN = -31.9021298869582
  -----------------------------------------------------------------------------
   IRC = 6	   R = 6.400000E-05	   PEN = -31.9625714995127
  -----------------------------------------------------------------------------
   IRC = 7	   R = 1.280000E-05	   PEN = -31.9800320324217
  -----------------------------------------------------------------------------
   IRC = 8	   R = 2.560000E-06	   PEN = -31.986356369792
  -----------------------------------------------------------------------------
   IRC = 9	   R = 5.120000E-07	   PEN = -31.9898563403618
  -----------------------------------------------------------------------------
                                                                               
三、优化结果__数据
===============================================================================
   罚函数构造次数          IRC = 10
  -----------------------------------------------------------------------------
   无约束优化方法调用次数  ITE = 10   一维搜索方法调用次数   ILI = 24
  -----------------------------------------------------------------------------
   惩罚函数值计算次数      NPE = 154   目标函数值计算次数     IFX = 0
  -----------------------------------------------------------------------------
   设计变量最优点 X*:
			X[1]= 1.001396E+00
			X[2]= 4.898729E+00
  -----------------------------------------------------------------------------
   最优值 F(X*)= -31.9919649346366
  -----------------------------------------------------------------------------
   最优点处的不等约束函数值 G(X*):
				GX[1]= -9.089605E-04
				GX[2]= -1.001396E+00
				GX[3]= -4.898729E+00
  -----------------------------------------------------------------------------
   最优点处的等式约束函数值 H(X*):
				HX[1]= 3.418881E-04
-------------------------------------------------------------------------------
--- STOP ---