共轭梯度优化方法MATLAB代码

应用该代码进行优化时，首先应将matlab的路径设置至该文件夹的子目录下，然后输入初值，输入函数，点击回车，即可求解

function [xvalk]=frcg（fungfunx0）
% 功能: 用FR共轭梯度法求解无约束问题:  min f（x）
%输入:  x0是初始点 fun gfun分别是目标函数和梯度
%输出:  x val分别是近似最优点和最优值  k是迭代次数.
maxk=5000;   %最大迭代次数
rho=0.6;sigma=0.4;
k=0;  epsilon=1e-4; 
n=length（x0）;
while（k    g=feval（gfunx0）;  %计算梯度
    itern=k-（n+1）*floor（k/（n+1））;
    itern=itern+1;
    %计算搜索方向
    if（itern==1）  
        d=-g;  
    else
        beta=（g‘*g）/（g0‘*g0）;
        d=-g+beta*d0;  gd=g‘*d;
        if（gd>=0.0）
            d=-g;  
        end
    end
    if（norm（g）    m=0; mk=0;
    while（m<20）   %Armijo搜索
        if（feval（funx0+rho^m*d）            mk=m; break;
        end
        m=m+1;
    end
    x0=x0+rho^mk*d;
    val=feval（funx0）;
    g0=g;  d0=d; 
    k=k+1;
end
x=x0;
val=feval（funx）; 

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

     文件        914  2009-09-22 05:34  共轭梯度法\frcg.m

     文件         52  2009-08-30 04:40  共轭梯度法\fun.m

     文件         81  2009-08-30 04:54  共轭梯度法\gfun.m

     目录          0  2012-06-23 13:54  共轭梯度法

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

                 1047                    4