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

1、坐标轮换法程序框图Vb运行界面Vb编程程序Private Sub Command1_Click() Dim s110 As DoubleDim s100 As DoubleDim s200 As DoubleDim s210 As DoubleDim x100(10000000) As Double Dim x110(10000000) As Double Dim x200(10000000) As Double Dim x210(10000000) As Double Dim 1 As DoubleDim 2 As DoubleDim k As LongDim a1#, a2#, a3#,

2、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)h0 = 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

3、.Text): a5 =Val(Text5.Text): a6 = Val(Text6.Text)e = Val(Text9.Text)k = 0m1: s100 = 1: s110 = 0'一维搜索求最优解h = h01 = 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

4、* 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 * xx

5、2 + 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 -

6、 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) >

7、 e1If yy1 < yy2 Thenn0 = 0b = 22 = 1yy2 = yy1Elsen0 = 1a = 11 = 2yy1 = yy2End IfIf n0 = 0 Then1 = a + 0.382 * (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 + a6Else2 = a + 0.618 * (b - a)xx1 = x100(k) + s100 * 2: xx

8、2 = x110(k) + s110 * 2: yy2 = a1 * xx1 2 + a2 * xx2 2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a6End IfLoop1 = (a + b) / 2x200(k) = x100(k) + s100 * 1x210(k) = x110(k) + s110 * 1s200 = 0: s210 = 1'一维搜索求最优解h = h01 = 0xx1 = x200(k) + s200 * 1: xx2 = x210(k) + s210 * 1: yy1 = a1 * xx1 2 + a2 * xx2

9、2 + a3 * xx1 * xx2 + a4 * xx1 + a5 * xx2 + a62 = 0 + hxx1 = x200(k) + s200 * 2: xx2 = x210(k) + s210 * 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 = x200(k) + s20

10、0 * 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 + a6Lo

11、opIf 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 *

12、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

13、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):

14、 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 * xx2 + a6Elsex1