追赶法求解三对角线性方程组Matlab编程.doc

《追赶法求解三对角线性方程组Matlab编程.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、追赶法求解三对角线性方程组Matlab编程对于n元线性方程组Ax=b，其系数矩阵A为非奇异的三对角矩阵（上下半带宽都是1的带状矩阵），此方程组可在Matlab中利用追赶法思想求解，Matlab编程如下：clear,clctemp=textread('exercise.txt','%n');n=temp(1);%读取系数矩阵的阶数[d,a,c,b]=textread('exercise.txt','%f%f%f%f','headerlines',1);fprintf('dacb')fori=1:1:n

2、fprintf('%9.1f%9.1f%9.1f%9.1f',d(i),a(i),c(i),b(i))endp(1)=a(1);fori=1:1:n-1q(i)=c(i)/p(i);p(i+1)=a(i+1)-d(i+1)*q(i);endy(1)=b(1)/p(1);fori=2:1:ny(i)=(b(i)-d(i)*y(i-1))/p(i);endx(n)=y(n);fori=n-1:-1:1x(i)=y(i)-q(i)*x(i+1);endx