最优化方法实验之阻尼牛顿法 #include #include # define eps 0.01 double f (double coe[], double x[]) { return coe[0]*pow(x[0],2)+coe[1]*pow(x[1],2)+coe[2]*x[0]*x[1]+coe[3]*x[0]+coe[4]*x[1]+coe[5]; } void grads (double coe[], double x[],double grads_x[]) { grads_x[0] = 2*coe[0]*x[0]+coe[2]*x[1]+coe[3]; grads_x[1] = 2*coe[1]*x[1]+coe[2]*x[0]+coe[4]; } void H (double coe[], double h[])//求解海赛矩阵 { h[0] = 2.0*coe[0]; h[1] = coe[2]; h[2] = coe[2]; h[3] = 2.0*coe[1]; }