无约束多维优化坐标轮换法vb编程.doc

坐标轮换法程序框图Vb运行界面Vb编程程序Private Sub Command1_Click()Dim s110 As DoubleDim s100 As DoubleDim s200 As DoubleDim s210 As DoubleDim x100(10000000) As DoubleDim x110(10000000) As DoubleDim x200(10000000) As DoubleDim x210(10000000) As DoubleDim 1 As DoubleDim 2 As DoubleDim k As LongDim a1#, a2#, a3#, a4#, a5#, a6#Dim h#, e1#, e#Dim 0#, 1#, 2#, 3#, a#, b#Dim xx1#, xx2#, yy1#, yy2#, yy3#Dim f2#, f1#, fxx#0 = Val(Text10.Text)h = Val(Text11.Text)e1 = Val(Text12.Text)x100(0) = Val(Text7.Text)x110(0) = Val(Text8.Text)a1 = Val(Text1.Text): a2 = Val(Text2.Text): a3 = Val(Text3.Text): a4 = Val(Text4.Text): a5 = Val(Text5.Text): a6 = Val(Text6.Text)e = Val(Text9.Text)k = 0m1: s100 = 1: s110 = 0一维搜索求最优解1 = 0xx1 = x100(k) + s100 * 1: xx2 = x110(k) + s110 * 1: yy1 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a62 = 0 + hxx1 = x100(k) + s100 * 2: xx2 = x110(k) + s110 * 2: yy2 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6If yy1 yy2 Thenh = 2 * hElseh = -0.25 * h3 = 11 = 22 = 3yy3 = yy1yy1 = yy2yy2 = yy3End If3 = 0 + hxx1 = x100(k) + s100 * 3: xx2 = x110(k) + s110 * 3: yy3 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6Do While yy2 yy3h = 2 * h1 = 2yy1 = yy22 = 3yy2 = yy33 = 0 + hxx1 = x100(k) + s100 * 3: xx2 = x110(k) + s110 * 3: yy3 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6LoopIf h 0 Thena = 1ya = yy1b = 3yb = yy3Elsea = 3b = 1ya = yy3yb = yy1End If1 = a + 0.382 * (b - a)2 = a + 0.618 * (b - a)xx1 = x100(k) + s100 * 1: xx2 = x110(k) + s110 * 1: yy1 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6xx1 = x100(k) + s100 * 2: xx2 = x110(k) + s110 * 2: yy2 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6Do While Abs(b - a) e1If yy1 yy2 Thenh = 2 * hElseh = -0.25 * h3 = 11 = 22 = 3yy3 = yy1yy1 = yy2yy2 = yy3End If3 = 0 + hxx1 = x200(k) + s200 * 3: xx2 = x210(k) + s210 * 3: yy3 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6Do While yy2 yy3h = 2 * h1 = 2yy1 = yy22 = 3yy2 = yy33 = 0 + hxx1 = x200(k) + s200 * 3: xx2 = x210(k) + s210 * 3: yy3 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6LoopIf h 0 Thena = 1ya = yy1b = 3yb = yy3Elsea = 3b = 1ya = yy3yb = yy1End If1 = a + 0.382 * (b - a)2 = a + 0.618 * (b - a)xx1 = x200(k) + s200 * 1: xx2 = x210(k) + s210 * 1: yy1 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6xx1 = x200(k) + s200 * 2: xx2 = x210(k) + s210 * 2: yy2 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6Do While Abs(b - a) e1If yy1 yy2 Thenn0 = 0b = 22 = 1yy2 = yy1Elsen0 = 1a = 11 = 2yy1 = yy2End IfIf n0 = 0 Then1 = a + 0.382 * (b - a)xx1 = x200(k) + s200 * 1: xx2 = x210(k) + s210 * 1: yy1 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6Else2 = a + 0.618 * (b - a)xx1 = x200(k) + s200 * 2: xx2 = x210(k) + s210 * 2: yy2 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6End IfLoop2 = (a + b) / 2x200(k + 1) = x200(k) + s200 * 2x210(k + 1) = x210(k) + s210 * 2xx1 = x200(k): xx2 = x210(k): f1 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6xx1 = x200(k + 1): xx2 = x210(k + 1): f2 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6If Abs(x200(k + 1) - x200(k) 2 + (x210(k + 1) - x210(k) 2) = e And Abs(f2 - f1) / f2) = e Thenxx1 = x200(k + 1): xx2 = x210(k + 1)fxx = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 *