阻尼牛顿法vb编程书.docVIP

题目：利用阻尼牛顿法求的极小点。 利用VB6.0进行编程，求得极小值点和极小值。具体程序如下说明。 1、求解原理 1) 给定初始点，收敛精度，置。 2) 计算、、和 3) 求，其中为沿进行一维有哪些信誉好的足球投注网站的最佳步长。 4) 检查收敛精度。若，则，停机；否则置，返回步骤2，继续进行进行有哪些信誉好的足球投注网站。 改进后的阻尼牛顿法程序框图如下： 2、程序如下： Function minα(x1#, x2#, d1#, d2#, α0#, h#, a#, b#, c#, d#, e#, f#, g#, i#, l#, m#, ε1#, j#, k#, n0#) α1 = α0 Xk1 = x1 + α1 * d1: Xk2 = x2 + α1 * d2 Yα1 = a * Xk1 ^ 4 + b * Xk2 ^ 4 + c * Xk1 ^ 3 + d * Xk2 ^ 3 + e * Xk1 ^ 2 + f * Xk2 ^ 2 + g * Xk1 + i * Xk2 + l * Xk1 * Xk2 + m α2 = α0 + h Xk1 = x1 + α2 * d1: Xk2 = x2 + α2 * d2 Yα2 = a * Xk1 ^ 4 + b * Xk2 ^ 4 + c * Xk1 ^ 3 + d * Xk2 ^ 3 + e * Xk1 ^ 2 + f * Xk2 ^ 2 + g * Xk1 + i * Xk2 + l * Xk1 * Xk2 + m If Yα1 Yα2 Then h = 2 * h Else h = -0.25 * h α3 = α1: Yα3 = Yα1 α1 = α2: Yα1 = Yα2 α2 = α3: Yα2 = Yα3 End If α3 = α0 + h Xk1 = x1 + α3 * d1: Xk2 = x2 + α3 * d2 Yα3 = a * Xk1 ^ 4 + b * Xk2 ^ 4 + c * Xk1 ^ 3 + d * Xk2 ^ 3 + e * Xk1 ^ 2 + f * Xk2 ^ 2 + g * Xk1 + i * Xk2 + l * Xk1 * Xk2 + m Do While Yα2 = Yα3 h = 2 * h α1 = α2: Yα1 = Yα2 α2 = α3: Yα2 = Yα3 α3 = α0 + h Xk1 = x1 + α3 * d1: Xk2 = x2 + α3 * d2 Yα3 = a * Xk1 ^ 4 + b * Xk2 ^ 4 + c * Xk1 ^ 3 + d * Xk2 ^ 3 + e * Xk1 ^ 2 + f * Xk2 ^ 2 + g * Xk1 + i * Xk2 + l * Xk1 * Xk2 + m Loop If h 0 Then j = α1: Yαa = Yα1 k = α3: Yαb = Yα3 Else j = α3: Yαa = Yα3 k = α1: Yαb = Yα1 End If α1 = j + 0.382 * (k - j) α2 = j + 0.618 * (k - j) Xk1 = x1 + α1 * d1: Xk2 = x2 + α1 * d2 Yα1 = a * Xk1 ^ 4 + b * Xk2 ^ 4 + c * Xk1 ^ 3 + d * Xk2 ^ 3 + e * Xk1 ^ 2 + f * Xk2 ^ 2 + g * Xk1 + i * Xk2 + l * Xk1 * Xk2 + m Xk1 = x1 + α2 * d1: Xk2 = x2 + α2 * d2 Yα2 = a * Xk1 ^ 4 + b * Xk2 ^ 4 + c * Xk1 ^ 3 + d * Xk2 ^ 3 + e * Xk1 ^ 2 + f * Xk2 ^ 2 + g * Xk1 + i * Xk2 + l * Xk1 * Xk2 + m Do While Abs(k - j) ε1 If Yα1 Yα2 Then n0 = 0 k = α2 α2 = α1: Yα2 = Yα1