1、無約束最優(yōu)化方法阻尼牛頓法Function Y( As Double) As DoubleDim x1 As Double, x2 As Double, s1 As Double, s2 As DoubleDim m1 As Double, m2 As Double, m3 As Double, m4 As Double, m5 As Double, m6 As Doublem1 = Val(Text1.Text): m2 = Val(Text2.Text)m3 = Val(Text3.Text): m4 = Val(Text4.Text)m5 = Val(Text5.Text): m6 =

2、Val(Text6.Text)Y = m1 * (x1 + * s1) 2 + m2 * (x2 + * s2) 2 + m3 * (x1 + * s1) * (x2 + * s2) + m4 * (x1 + * s1) + m5 * (x2 + * s2) + m6End FunctionFunction yiweisosuo( As Double) As DoubleDim x1#, x2#, s1#, s2#,#, h#,0#,1#,2#, A1#, A2#, A3#,1#, a#, b#, yA1#, yA2#, yA3#h = Val(Text11.Text): 1 = Val(Te

3、xt12.Text)A1 =0: yA1 = Y(A1)A2 =0 + h: yA2 = Y(A2)If yA2 = yA1 Then h = -1 * hA3 = A1: yA3 = yA1A1 = A2: yA1 = yA2A2 = A3: yA2 = yA3Else h = 2 * hEnd IfA3 =0 + h: yA3 = Y(A3)Do Until yA2 A2 Thena = A2: b = A1Elsea = A1: b = A2End If2 = a + 0.618 * (b - a): 1 = b - 0.618 * (b - a)y1 = Y(1): y2 = Y(2

4、)Do While Abs(b - a) 1If y1 s1 = -(haiseni(0, 0) * tidu(0, 0) + haiseni(0, 1) * tidu(1, 0)s2 = -(haiseni(1, 0) * tidu(0, 0) + haiseni(1, 1) * tidu(1, 0) = yiweisosuo()x1 = x1 + * s1: x2 = x2 + * s2f = m1 * x1 2 + m2 * x2 2 + m3 * x1 * x2 + m4 * x1 + m5 * x2 + m6tidu(0, 0) = 2 * m1 * x1 + m3 * x2 + m4tidu(1, 0) = 2 * m2 * x2 + m3 * x1 + m5LoopText13.Text = x1= & x1Text14.Text = x