共轭梯度法求解方程程序.doc

求的极小点，取初始点 = 。 用matlab程序如下： function mg=MG() % %共轭梯度法求解% clc; clear; n=2; x=[2 0]; max_k=100; count_k=1; trace(1,1)=x(1); trace(2,1)=x(2); trace(3,1)=f_fun(x); k=0; g1=f_dfun(x); s=-g1; while count_k=max_k if k==n g0=f_dfun(x); s=-g0; k=0; else r_min=fminbnd(@(t) f_fun(x+t*s),-100,100); x=x+r_min*s; g0=g1; g1=f_dfun(x); if norm(g1)10^(-6) break; end m=(norm(g1)^2)/(norm(g0)^2); s=-g1+m*s; count_k=count_k+1; trace(1,count_k)=x(1); trace(2,count_k)=x(2); trace(3,count_k)=f_fun(x); k=k+1; end end count_k x f=f_fun(x) function g=f_dfun(x) g(1,1)=48*x(1)^3-24*x(1)*x(2)+8*x(1)-4; g(2,1)=6*x(2)-12*x(1)^2; function f=f_fun(x) f=12*x(1)^4-12*x(1)^2*x(2)+4*x(1)^2-4*x(1)+1+3*x(2)^2; 运行结果如下： k = 9 = 0.5000 0.5000 = -1.1102e-016 从上述运行结果可以得??：最优解为=0.5 ，最小值约为= -1.1102e-016。