MATLAB 追赶法求解三对角方程组的算法原理例题与程序

MATLAB 追赶法求解三对角方程组的算法原理例题与程序

ID:41140169

大小:40.01 KB

页数:3页

时间:2019-08-17

MATLAB 追赶法求解三对角方程组的算法原理例题与程序_第1页
MATLAB 追赶法求解三对角方程组的算法原理例题与程序_第2页
MATLAB 追赶法求解三对角方程组的算法原理例题与程序_第3页
资源描述:

《MATLAB 追赶法求解三对角方程组的算法原理例题与程序》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、3)三对角形线性方程组二、数学原理设系数矩阵为三对角矩阵则方程组Ax=f称为三对角方程组。设矩阵A非奇异,A有Crout分解A=LU,其中L为下三角矩阵,U为单位上三角矩阵,记   可先依次求出L,U中的元素后,令Ux=y,先求解下三角方程组Ly=f得出y,再求解上三角方程组Ux=y。事实上,求解三对角方程组的2追赶法将矩阵三角分解的计算与求解两个三角方程组的计算放在一起,使算法更为紧凑。其计算公式为:(*)三、程序设计functionx=chase(a,b,c,f)%求解线性方程组Ax=f,其中A是三对角阵%a是矩阵A的下对角线元素a(1)=0%b是矩阵

2、A的对角线元素%c是矩阵A的上对角线元素c(n)=0%f是方程组的右端向量n=length(f);x=zeros(1,n);y=zeros(1,n);d=zeros(1,n);u=zeros(1,n);%预处理d(1)=b(1);fori=1:n-1u(i)=c(i)/d(i);d(i+1)=b(i+1)-a(i+1)*u(i);end%追的过程y(1)=f(1)/d(1);fori=2:ny(i)=(f(i)-a(i)*y(i-1))/d(i);end%赶的过程x(n)=y(n);fori=n-1:-1:1x(i)=y(i)-u(i)*x(i+1);en

3、d>>a=[0,-1,-1,-1,-1,-1,-1,-1,-1,-1];>>b=[4,4,4,4,4,4,4,4,4,4];>>c=[-1,-1,-1,-1,-1,-1,-1,-1,-1,0];>>f=[7,5,-13,2,6,-12,14,-4,5,-5];>>x=chase(a,b,c,f)x=2.00001.0000-3.00000.00001.0000-2.00003.0000-0.00001.0000-1.0000四、结果分析和讨论追赶法求解的结果为x=(2,1,-3,0,1,-2,3,0,1,-1)T。求解结果与精确解一样,这表明追赶法对于求解

4、三对角方程组具有非常高的精度,误差非常小。算法次数也较少,不选主元也可以有效的算出精确结果,是一种计算量少而数值稳定的方法。五、完成题目的体会与收获通过本题的求解,深刻的理解了追赶法求解三对角方程组的算法原理。学会了追赶法的matlab编程,学会了又一种求解特殊方程组的方法。追赶法在计算量方面有着巨大的优势,因此在可能的情况下应优先使用追赶法。加深了对数值计算教材知识的理解,收获非常大。

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。