数值分析最速下降法MATLAB代码

%A为系数矩阵，b为等式右侧元素，x0为迭代初值，N代表最大迭代次数，p为误差允许范围
%其中需要输入x0、N、p
%需要输入nA，改变nA的值来改变方程组阶数
A=zeros(nA);
A(1,:)=-1;
A(:,1)=-1;
A(nA,:)=-1;
A(:,nA)=-1;
	for i=1:nA
	A(i,i)=30;
	end
A(1,nA)=0; 
A(nA,1)=0;
%以上构建方程组A的矩阵
nA=size(A,1);
b=zeros(nA,1);
for i=1:nA
    b(i,1)=2*i;
end
%以上构建方程组等式右侧矩阵b
x0=transpose(x0);  %转置
r=b-A*x0;
k=0;
while k<=N
    r=b-A*x0;
    q=dot(r,r)/dot(A*r,r);
    x=x0+q*r;%迭代公式
    k=k+1;
    h=norm((x-x0),inf);
    if h<=p
        break
    end    
    x0=x;
end
if k>N
    disp(['迭代次数= ',num2str(k),',算法超过最大迭代次数。']);    
else
    disp(['迭代次数= ',num2str(k)]);    
end
x=x0;