目标
原函数: f(x) = 2x1x1 + x2x2 -2x1*x2 - 4x1 + 4
起始点 为(0,0)
精度 :误差要求小于 10**(-3)
过程:
直接上代码看注释吧()
import numpy as np
import math
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
def f(x):
return 2*x[0]**2 + x[1]**2 - 2*x[0]*x[1]-4*x[0]+4
def step(x): # 计算每次迭代的步长,H为该函数的hessian矩阵
H = np.matrix([[4, -