matlab 三对角矩阵 追赶法 计算方法追赶法解三对角矩阵方程

计算方法实验报告

实验：追赶法解三对角矩阵方程

1.实验目的

①熟悉matlab 软件的使用

②掌握追赶法的基本方法

③编程实现追赶法解三对角矩阵方程

2.实验内容

題目：用追赶法解方程组

????????????????--------4114114114114????????????????54321x x x x x =???????

?????????100200200200100 3.求解：

1.追赶法代码

function [x]=ZhuiGan(a,b,c,f)

r=size(a);

m=r(2);

r=size(b);

n=r(2);

if size(a)~=size(b)|(m~=n-1)|size(b)~=size(f)

error('变量不匹配,检查变量输入情况!');

end

p=ones(1,m);

Y=ones(1,n);

x=Y;

p(1)=a(1)/b(1);

Y(1)=f(1)/b(1);

t=0;

for i=2:m

t=b(i)-a(i-1)*p(i-1);

p(i)=c(i)/t;

Y(i)=(f(i)-a(i-1)*Y(i-1))/t;