追赶法求解三对角线性方程组

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

作业五:编写用追赶法求解三对角线性方程组的标准程序,并求下列方程组的解。matlab1233212417259xxxfunction[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)=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-1z=z+L(k,q)*U(q,j);endU(k,j)=A(k,j)-z;endifabs(U(k,k))epsflag='failure';return;endfori=k+1:nz=0;forq=1:k-1z=z+L(i,q)*U(q,k);endL(i,k)=(A(i,k)-z)/U(k,k);endend流程图开始赋值LU=A分解由Ly=b求y由Ux=y求出x输出x运行结果

1 / 4
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功