数值分析3-平方根法求解线性方程组

《数值分析3-平方根法求解线性方程组》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、作业四：编写用平方根法求解线性方程组的标准程序，并求下列方程组的解。621-12410114-1-10-13x1x2x3x4=-1213流程图开始计算下三角矩阵G的数据值计算中间变量y计算x的值结束程序A=[621-1;2410;114-1;-10-13];[n,n]=size(A);G=zeros(n,n);G(1,1)=sqrt(A(1,1));G(2:n,1)=A(2:n,1)/G(1,1);forj=2:n-1G(j,j)=sqrt(A(j,j)-sum(G(j,1:j-1).^2));fori=j+1:nG(i,j)=(A(i,j)-sum(G(i,1:j-1).*G(j,1:

2、j-1)))/G(j,j);endendG(n,n)=sqrt(A(n,n)-sum((G(n,1:n-1)).^2));b=[-1213];x=zeros(n,1);%未知向量y=zeros(n,1);%中间向量y(1)=b(1)/G(1,1);fori=2:ny(i)=(b(i)-sum(G(i,1:i-1)'.*y(1:i-1)))/G(i,i);endG=G';x(n)=y(n)/G(n,n);fori=n-1:-1:1x(i)=(y(i)-sum(G(i,i+1:n)'.*x(i+1)))/G(i,i);end程序结果