最优化方法实验之最速下降法
#include<math.h>
#include<stdio.h>
# 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];
}