追赶法(经典计算).doc

《追赶法(经典计算).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一、算法理论在一些实际问题中，例如解常微分方程边值问题，解热传导方程以及船体数学放样中建立三次样条函数等，都会要求解系数矩阵为对角占优的三对角线方程组，简记为.求解等价于解两个三角形方程组：；.从而得到解三对角线方程组的追赶法公式：(1)计算的递推公式(2)解(3)解我们将计算系数将计算方程组的解一、算法框图一、算法程序#include#include#include#defineN20doublea[N],b[N],c[N-1],f[N],r[N];intn;v

2、oidLUDecompose();//LU分解voidbackSubs();//回代voidmain(){printf("请输入方程的维数n＝");scanf("%d",&n);getchar();if(n>N

3、

4、n<=0){printf("由于该维数过于犀利,导致程序退出!");return;}printf("输入下三角元素");printf("输入%d个a值:",n-1);for(inti=1;i

5、n");printf("输入%d个b值:",n);for(i=0;i

6、printf("线性方程组的解为:");for(i=0;i

7、储空间f[0]=f[0]/b[0];for(inti=1;i=0;i--)f[i]=f[i]-c[i]*f[i+1];}一、算法实现例1．用该程序计算三对角线方程组，计算其方程组的解。解：运行程序（1）显示出请输入方程组的维数：n=5.，回车。（2）显示出请输入下三角元素输入4个a值：-1-1-1-1，回车。（3）显示出请输入主对角线元素输入5个b值：22222，回车。（4）显示出请输入上三角元

8、素输入4个c值：-1-1-1-1，回车。（5）显示出请输入5个方程组右端顶：10000，回车。其解为例2．用该程序计算三对角线方程组，计算其方程组的解。解：运行程序（1）显示出请输入方程组的维数：n=5.，回车。（2）显示出请输入下三角元素输入4个a值：1111，回车。（3）显示出请输入主对角线元素输入5个b值：33333，回车。（4）显示出请输入上三角元素输入4个c值：1111，回车。（5）显示出请输入5个方程组右端顶：11000，回车。其解为