追赶法解三对角阵方程组

原创 2007年09月17日 18:30:00
//追赶法解三对角阵方程组Ax=b。
 //n---方程组的阶数。
 //a1,a2,a3依次为三对角阵的三对角。
 //成功返回1,失败返回0。
 int RunEq(int n,double *a1,double *a2,double *a3,double *b,double *x)
 {
  int i;
  double l;
  for(i=1;i<n;i++)
  {
   if(fabs(a2[i-1])<0.00000000001) return 0;
   l=a1[i]/a2[i-1];
   a2[i]-=a3[i-1]*l;
   b[i]-=b[i-1]*l;
  }
  if(fabs(a2[n-1])<0.00000000001) return 0;
  x[n-1]=b[n-1]/a2[n-1];
  for(i=n-2;i>=0;i--)
  {
   x[i]=(b[i]-a3[i]*x[i+1])/a2[i];
  }
  return 1;
 }
 
注意:a1的第一个元素需要空出来,因为矩阵第一行不含第一各对角的元素,同样a3的最后一个元素也要空出来

相关文章推荐

追赶法求解三对角方程组

在这篇文章里,我们介绍追赶法的基本原理,以及用追赶法求解三对角方程组的算法....

追赶法解三对角方程组

  • 2010年07月23日 20:40
  • 2KB
  • 下载

解三对角矩阵的追赶法

  • 2014年01月09日 11:41
  • 635B
  • 下载

三对角阵的LU分解和三对角方程组的求解(C语言)

/*三对角阵的LU分解和三对角方程组的求解 -------------A=LU的分解算法------- 参考教材:《数值分析》李乃成,梅立泉,科学出版社     《计算方法教程》第二版 凌永...

MATLAB追赶法求解三对角方程

  • 2017年03月28日 22:13
  • 350B
  • 下载

Matlab实现——严格对角占优三对角方程组求解(高斯赛尔德Gauss-Seidel迭代、超松弛)

严格对角占优三对角方程组求解 对中等规模的n阶的(n)线性方程组,直接法的准确性和可靠性,所以常采用直接法 对于较高阶的方程组,特别是地于某些偏微分方程离散化后得到的大型稀疏方程组(系统矩 阵绝...

追赶法求解三对角矩阵

  • 2013年06月18日 21:12
  • 431B
  • 下载

追赶法解线性一般带型方程组

  • 2009年04月17日 20:54
  • 750B
  • 下载

Guass列主元、平方根法、追赶法求解方程组的C++实现

一,要解决的问题选用合适的算法,求解三种线性方程组:一般线性方程组,对称正定方程组,三对角线性方程组。 方程略。二,数值方法1,使用Guass列主元消去法求解一般线性方程组。Guass列主元是为了防...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:追赶法解三对角阵方程组
举报原因:
原因补充:

(最多只允许输入30个字)