数值优化中共轭梯度法程序的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

function [x,minf]=mingt(f,x0,var,eps)
format long;
syms s;
if nargin==3
    eps=1.0e-6;
end
g=gradient(f);
n=length(var);
g0=Funval(g,var,x0);
p0=-g0;
k=0;
while 1
    g0=Funval(g,var,x0);
    t0=norm(g0);
    if t0<=eps
        x=x0;
        break;
    else        
        y=x0+s*p0';
        z=Funval(f,var,y);
        h=minndf(z,0);
         x1=x0+h*p0';
        g1=Funval(g,var,x1);
        t1=norm(g1);      
    end
        if t1<=eps
            x=x1;
            break;
        else if k+1==n               
                x0=x1;
                g0=Funval(f,var,x0);
                p0=-g0;
            else  r=t1^2./t0^2;
                  p0=-g1+r*p0;
                  k=k+1;
                  x0=x1;
            end
       end
end

minf=Funval(f,var,x);
format short;

syms t s;
f=3*t^2+s^2-2*t*s+4*t+3*s;
[x,minf]=mingt(f,[0 0],[t s])

为啥老是出不来？麻烦大神给我看一下，到底哪里编错了？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：数值优化 function gradient format nargin function gradient