共轭梯度法matlab实例.docx

《共轭梯度法matlab实例.docx》由会员分享，可在线阅读，更多相关《共轭梯度法matlab实例.docx（2页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上clear all;close all;clc;syms x1 x2 t;f=2*x12+x22-x1*x2;f_grad=diff(f,x1);diff(f,x2);X0=1;1;n=10;epsonal=0.01;fx=inline(f);fx_grad=inline(f_grad);X=X0;fx0=fx(X(1),X(2);fx0_grad=fx_grad(X(1),X(2);while 1 if fx0_grad.*fx0_grad=epsonal break; end P0=-fx0_grad; k=0; while 1 Xk=X+t*P0; fx1=fx

2、(Xk(1),Xk(2); tk,y=equation_extremum(fx1,t,-1,5,epsonal); Xk=X+tk*P0; fx0_k=fx(Xk(1),Xk(2); fx0_grad_k=fx_grad(Xk(1),Xk(2); if fx0_grad_k.*fx0_grad_k=epsonal fx0_grad=fx0_grad_k; break; end if k+1=n X=X; else lamdak=fx0_grad_k.*fx0_grad_k/(fx_grad(X(1),X(2).*fx_grad(X(1),X(2); P0=-fx0_grad_k+lamdak*P0; X=Xk; k=k+1; end endendXkfx0_k运行结果：Xk = 1.0e-03 * 0.0486 0.5350fx0_k = 2.6495e-07专心-专注-专业