最优化理论变尺度法Matlab实现

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

function [mHdaXg1fv] = DFP（GbcXe）
%X：为初始点，列向量表示
%G：对称正定矩阵；2*2
%b和c：函数参数，列向量表示
%e：精度
%m：最后导数值（精度）；H：参数H用于求方向；g1：导数值数组
%a：各阶段步长；d：迭代方向；X：各阶段迭代点；fv：函数极小值
if nargin<5
    error（‘输入参数必须为5‘）;
end
n=length（G）;
if n==2
format long
syms x1 x2 lema
f=1/2*[x1x2]*G*[x1;x2]+b‘*[x1;x2]+c;
f=f（11）;
g=[diff（fx1）;diff（fx2）];

g1=subs（subs（gx1X（11））x2X（21））;%代入x0
d=-g1;%第一个搜索方向
%a=-（d‘*g1）/（d‘*G*d）;%第一步步长

%----------中点法求零点----------------------------------
    df=[g1（:1）+lema*G*d（:1）]‘*d（:1）%输出代入具体函数值的关于步长的一阶导数，方便输入区间
    df=inline（df）;
    lf=-1;
    rt=1;
    while abs（rt-lf）>e
        mid=（lf+rt）/2;
        if df（mid）==0
            break;
        else if df（mid）*df（rt）<0
                lf=mid;
            else
                rt=mid;
          

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2438  2013-11-18 23:47  变尺度法刘敏\DFP.m

     文件       1376  2013-11-17 22:19  变尺度法刘敏\DFP2.m

     文件     814080  2013-11-14 23:21  变尺度法刘敏\例题是唐老师教材_1025_拟牛顿法（DFP算法）.ppt

     文件     309738  2013-11-19 00:41  变尺度法刘敏\变尺度法课程设计.doc

    I.A....     77824  2013-11-17 23:40  变尺度法刘敏\最优化课程设计要求.doc

     目录          0  2013-11-19 01:19  变尺度法刘敏

----------- ---------  ---------- -----  ----

              1205456                    6