欢迎来到天天文库
浏览记录
ID：56998384
大小：108.50 KB
页数：5页
时间：2020-07-30
《追赶法求解三对角线性方程组.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、作业五：编写用追赶法求解三对角线性方程组的标准程序，并求下列方程组的解。matlabfunction[x,flag]=zgf(A,b)%A为方程组的系数矩阵%b为方程组的右端项%x为方程组的解%flag为指标向量，flag=‘failure’表示失败，flag=‘OK’表示成功[n,m]=size(A);n=length(b);[L,U,flag]=LU_decom(A);fori=1:nifi==1y(i)=b(i)/L(i,i);elsey(i)=(b(i)-L(i,i-1)*y(i-1))/L(i,i);endendx(n)=
2、y(n)/U(n,n);fork=n-1:-1:1y(k)=y(k)-U(k,k+1)*x(k+1);x(k)=y(k)/U(k,k);end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function[L,U,flag]=LU_decom(A)[n,m]=size(A);ifn~=merror('TherowsandcolumnsofmatrixAmustbeequal!');return;endL=eye(n);U=zeros(n);flag='OK';fork=1:nforj=k:nz=0;forq=1:k-
3、1z=z+L(k,q)*U(q,j);endU(k,j)=A(k,j)-z;endifabs(U(k,k))
此文档下载收益归作者所有
