【文件属性】:
文件名称:最优化方法实验之阻尼牛顿法
文件大小:2KB
文件格式:CPP
更新时间:2013-03-22 09:49:12
阻尼牛顿法
最优化方法实验之阻尼牛顿法
#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];
}