s02.rtf

来自「机械优化设计中的惩罚函数法」· RTF 代码 · 共 87 行

                          常用优化方法——惩罚函数法 
                          ^^^^^^^^^^^^^^^^^^^^^^^^^^ 
                                                                               
一、初始数据
===============================================================================
   设计变量个数     N = 2
  -----------------------------------------------------------------------------
   不等式约束个数  KG = 5       等式约束个数      KH = 0
  -----------------------------------------------------------------------------
   惩罚因子         R = 1       惩罚因子降低系数   C = 0.2
  -----------------------------------------------------------------------------
   初始步长        T0 = 0.01       收敛精度         EPS = 1E-5
  -----------------------------------------------------------------------------
   无约束优化方法: POWELL法
  -----------------------------------------------------------------------------
   设计变量初始点 X0:
			X[1]=1
			X[2]=1
  -----------------------------------------------------------------------------
   设计变量下界 BL:
			BL[1]=0
			BL[2]=0
  -----------------------------------------------------------------------------
   设计变量上界 BU:
			BU[1]=8
			BU[2]=8
  -----------------------------------------------------------------------------
   初始点目标函数值 F(X0)= 47
  -----------------------------------------------------------------------------
   初始点处的不等约束函数值 G(X0):
				GX[1]= -1.000000E+00
				GX[2]= -1.000000E+00
				GX[3]= -5.000000E+00
				GX[4]= -7.000000E+00
				GX[5]= -9.000000E+00
-------------------------------------------------------------------------------
                                                              
二、计算过程__数据
===============================================================================
   IRC = 0	   R = 1.000000E+00	   PEN = 49.4539682539683
  -----------------------------------------------------------------------------
   IRC = 1	   R = 2.000000E-01	   PEN = 16.3468803349782
  -----------------------------------------------------------------------------
   IRC = 2	   R = 4.000000E-02	   PEN = 13.1104311933622
  -----------------------------------------------------------------------------
   IRC = 3	   R = 8.000000E-03	   PEN = 11.8747845410832
  -----------------------------------------------------------------------------
   IRC = 4	   R = 1.600000E-03	   PEN = 11.3718892964575
  -----------------------------------------------------------------------------
   IRC = 5	   R = 3.200000E-04	   PEN = 11.1602767196316
  -----------------------------------------------------------------------------
   IRC = 6	   R = 6.400000E-05	   PEN = 11.0697619502954
  -----------------------------------------------------------------------------
   IRC = 7	   R = 1.280000E-05	   PEN = 11.0304136462488
  -----------------------------------------------------------------------------
   IRC = 8	   R = 2.560000E-06	   PEN = 11.0133419926234
  -----------------------------------------------------------------------------
   IRC = 9	   R = 5.120000E-07	   PEN = 11.0058788624432
  -----------------------------------------------------------------------------
   IRC = 10	   R = 1.024000E-07	   PEN = 11.0026236818822
  -----------------------------------------------------------------------------
   IRC = 11	   R = 2.048000E-08	   PEN = 11.0011577505621
  -----------------------------------------------------------------------------
                                                                               
三、优化结果__数据
===============================================================================
   罚函数构造次数          IRC = 12
  -----------------------------------------------------------------------------
   无约束优化方法调用次数  ITE = 32   一维搜索方法调用次数   ILI = 72
  -----------------------------------------------------------------------------
   惩罚函数值计算次数      NPE = 337   目标函数值计算次数     IFX = 0
  -----------------------------------------------------------------------------
   设计变量最优点 X*:
			X[1]= 5.999894E+00
			X[2]= 4.995081E+00
  -----------------------------------------------------------------------------
   最优值 F(X*)= 11.0003413701401
  -----------------------------------------------------------------------------
   最优点处的不等约束函数值 G(X*):
				GX[1]= -5.999894E+00
				GX[2]= -4.995081E+00
				GX[3]= -1.058933E-04
				GX[4]= -3.004919E+00
				GX[5]= -5.025241E-03
-------------------------------------------------------------------------------
--- STOP ---