最优化实验之变尺度法 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 SD(double h[],double grads_x[],double s[])//定义搜索方向向量函数 { s[0]=-(h[0]*grads_x[0]+h[1]*grads_x[1]); s[1]=-(h[2]*grads_x[0]+h[3]*grads_x[1]); } double findlambda (double coe[], double x[], double s[]) //求解lambda { return(-(2*coe[0]*x[0]*s[0]+2*coe[1]*x[1]*s[1]+coe[2]*x[0]*s[1]+coe[2]*x[1]*s[0]+coe[3]*s[0]+coe[4]*s[1])/(2*coe[0]*s[0]*s[0]+2*coe[1]*s[1]*s[1]+2*coe[2]*s[0]*s[1])); } void H (double h[],double x[],double temp_x[],double grads_x[],double temp_g[] )//求解变化矩阵H

动摇白拟牛顿法(变尺度法)DFP算法的cc++源码.doc

拟牛顿法(变尺度法)DFP算法的cc++源码.txt

变尺度法的dfp算法，包括确定搜索区间和步长，适合二元二次的目标函数

C++编写的拟牛顿法的源程序，非常实用的东东哦~~

变尺度法MATLAB仿真，代码后面有详细说明，每一步都有注释。变尺度法是求函数极小值的一种方法

输入多项式函数的系数矩阵及迭代精度，得到最后的函数的最小值

包括变尺度法 方向加速法 步长加速 公轭梯度法等代码，有界面，matlab程序

DFP法不仅综合了梯度法和牛顿法的优点，而且 具有对二次函数有限步收敛的性质，是无约束优化方法中最 有效的方法之一。