matlab最速下降法和共轭梯度法

function [x,c]=ZSXJF(A,b,eps,N,x0)

c=0;x=x0;r0=b;

while(r0'*r0>=eps&c<=N)

r0=b-A*x;

t0=r0'*r0/(r0'*A*r0);

x=x+t0*r0;

c=c+1;

end

b0=p0'*A*r0/(p0'*A*p0);

p0=r0-b0*p0;

c=c+1;

end

输 入： A=[4 -1 0 -1 0 0;-1 4 -1 0 -1 0;0 -1 4 0 0 -1;-1 0 0 4 -1 0;0 -1 0 -1 4 -1;0 0 -1 0 -1

4];b=[0;5;0;6;-2;6];x0=zeros(size(b));eps=0.00001;N=25;[x,c]=ZSXJF(A,b,eps,N,x0);

或者[x,c]=GETDF(A,b,eps,N,x0);

输出：…