追赶法解三对角阵方程组

原创 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的最后一个元素也要空出来

追赶法求解三对角方程组

在这篇文章里,我们介绍追赶法的基本原理,以及用追赶法求解三对角方程组的算法....
  • u010450214
  • u010450214
  • 2015年12月04日 09:37
  • 7475

数值分析 追赶法求解三对角线性方程组 MATLAB实现

函数主体部分编程算法  参考 数值分析 第四版 颜庆津 P27 运行结果截图: %追赶法求解三对角线性方程组,Ax=b,A用一维数组a,c,d存储。 function [L,U,x]=cr...
  • jingmiaa
  • jingmiaa
  • 2015年11月05日 11:48
  • 4718

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

一,要解决的问题选用合适的算法,求解三种线性方程组:一般线性方程组,对称正定方程组,三对角线性方程组。 方程略。二,数值方法1,使用Guass列主元消去法求解一般线性方程组。Guass列主元是为了防...
  • qq_24421591
  • qq_24421591
  • 2016年03月28日 22:05
  • 1433

特殊矩阵——三对角矩阵(Tridiagonal Matrix)

特殊矩阵——三对角矩阵(Tridiagonal Matrix) 1. 三对角矩阵的概念 三对角矩阵就是对角线、邻近对角线的上下次对角线上有元素,其他位置均为0的矩阵。 三对角矩阵是一种特殊的上He...
  • cainv89
  • cainv89
  • 2016年06月27日 13:06
  • 1775

三对角矩阵压缩存储--注意对角元素的下标

三对角矩阵压缩存储--注意对角元素的下标 对角矩阵的压缩存储   对角矩阵是指所有非零元素全部集中在中心几条对角线上的矩阵。下面以三对角矩阵(所有非零元素集中在中心三条对角线上)为例描述对角矩阵的...
  • moyanmalin
  • moyanmalin
  • 2016年08月13日 09:40
  • 2299

对角阵、特征值与特征向量

如何将矩阵P拆分成带有对角阵Λ的形式?
  • constantin_ouc
  • constantin_ouc
  • 2017年12月05日 01:40
  • 252

利用牛顿迭代法求解非线性方程组

最近一个哥们,是用牛顿迭代法求解一个四变量方程组的最优解问题,从网上找了代码去改进,但是总会有点不如意的地方,迭代的次数过多,但是却没有提高精度,真是令人揪心!        经分析,发现是这个方程...
  • liuguiyangnwpu
  • liuguiyangnwpu
  • 2014年05月08日 22:09
  • 31795

正交相似变换

线性空间的度量 首先用Hα=[hij]n∗nH_{\alpha}=\left [h_{ij} \right ]_{n*n}来表示线性空间的度量,其中hij=G(αi,αj)h_{ij}=G(\al...
  • zhengwei223
  • zhengwei223
  • 2017年12月27日 21:53
  • 103

对角矩阵的性质(diagonal matrix)

对角矩阵(英语:diagonal matrix)是一个主对角线之外的元素皆为0的矩阵。对角线上的元素可以为0或其他值。因此n行n列的矩阵{\displaystyle \mathbf {D} } =...
  • taoqick
  • taoqick
  • 2017年11月15日 21:53
  • 1054

数值计算——线性方程组的迭代法

数值计算——线性方程组的迭代法 与上一篇中的直接法相比,迭代法是从解的一个初始估计值除法,逐步对他进行改进,知道到达所需的精度,理论上来说,经过无限次的迭代之后就可以得到真解,但实际上只需要达到所需要...
  • qq_29721419
  • qq_29721419
  • 2017年03月30日 16:53
  • 841
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:追赶法解三对角阵方程组
举报原因:
原因补充:

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