%共軛剃度法極小化函數(shù)的通用子函數(shù)，gonge.m
%輸入變量為初始的迭代點(diǎn)，輸出變量為極小值點(diǎn)
function x0=gonge(x)
%判斷梯度范數(shù)是否滿足計(jì)算精度的要求.是,標(biāo)志變量設(shè)為1，輸出結(jié)果;否,標(biāo)志變量設(shè)為0，繼續(xù)迭代
if sum(abs(tidu(x)).^2)<1e-4
    flag=1;
    x0=x;
else 
    flag=0;
end
%第一次的迭代方法為負(fù)梯度方向
d1=-tidu(x);
a=gold(x,d1);
x1=x+a*d1;
%循環(huán)求解函數(shù)的極小點(diǎn)
while flag==0
    g1=tidu(x);
    g2=tidu(x1);
    %利用FR公式求解系數(shù)bata
    bata=(g2*g2')/(g1*g1');
    d2=-g1+bata*d1;
    a=gold(x1,d2);
    x=x1;
    x1=x+a*d2;
    %判斷下一次和上一次迭代點(diǎn)之差是否滿足計(jì)算精度的要求.是,標(biāo)志變量設(shè)為1，輸出結(jié)果;否,標(biāo)志變量設(shè)為0，繼續(xù)迭代
    if sum(abs(x1-x).^2)<1e-6 
       flag=1;
       x0=x1;
    else 
       flag=0;
    end
end