实用范文网

首页

追赶法解线性方程组

追赶法解线性方程组(C语言版)

C 语言 2010-06-17 20:23:00 阅读6 评论0 字号：大中小

#include

#define N 3

main()

{

double A[N][N]={{4,-1,0},{-1,4,-1},{0,-1,4}},b[N]={2,4,10};

int i;

A[0][1]=A[0][1]/A[0][0];

for(i=1;i

A[i][i+1]=A[i][i+1]/(A[i][i]-A[i][i-1]*A[i-1][i]);

for(i=1;i

A[i][i]=A[i][i]-A[i][i-1]*A[i-1][i];

b[0]=b[0]/A[0][0];

for(i=1;i

b[i]=(b[i]-A[i][i-1]*b[i-1])/A[i][i];

for(i=N-2;i>=0;i--)

b[i]=b[i]-A[i][i+1]*b[i+1];

for(i=0;i

printf("%.6lf ",b[i]);

}

原题：

追赶法解线性方程组(C语言版)

C 语言 2010-06-17 20:23:00 阅读6 评论0 字号：大中小

#include

#define N 3

main()

{

double A[N][N]={{4,-1,0},{-1,4,-1},{0,-1,4}},b[N]={2,4,10};

int i;

A[0][1]=A[0][1]/A[0][0];

for(i=1;i

A[i][i+1]=A[i][i+1]/(A[i][i]-A[i][i-1]*A[i-1][i]);

for(i=1;i

A[i][i]=A[i][i]-A[i][i-1]*A[i-1][i];

b[0]=b[0]/A[0][0];

for(i=1;i

b[i]=(b[i]-A[i][i-1]*b[i-1])/A[i][i];

for(i=N-2;i>=0;i--)

b[i]=b[i]-A[i][i+1]*b[i+1];

for(i=0;i

printf("%.6lf ",b[i]);

}

原题：

相关内容

迭代法实验

实验五线性方程组的迭代法实验一. 实验目的 (1)深入理解线性方程组的迭代法的设计思想,学会利用系数矩阵的性质以保证迭代过程的收敛性,以及解决某些实际的线性方程组求解问题. (2)熟悉Matlab编程环境,利用Matlab解决具体的方程求根问题. 二. 实验要求建立Jacobi迭代公式.Ga ...

数值计算基础

数值计算基础实验指导书 2010年目录实验一直接法解线性方程组的 ................................ 1 实验二插值方法 ........................................... 10 实验三数值积分 ............. ...

第二章解线性代数方程组的直接法

第二章解线性方程组的直接法本章研究的对象是n 阶线性方程组 ⎧a 11x 1+a 12x 2+... +a 1n x n =b 1⎪a x +a x +... +a x =b ⎪2112222n n 2⎨⎪ ⎪⎩a n 1x 1+a n 2x 2+... +a nn x n =b n 其矩阵形式 ...

预处理共轭梯度法求解线性方程组

[摘要]针对共轭梯度法求解线性方程组,提出一种预处理思想.基于次思想,首先给出预处理矩阵,然后求解预处理线性方程组,再使用共轭梯度法求解.最后通过几个数值试验,与直接使用共轭梯度法求解线性方程组相比较,本文的方法提高了收敛速度. [关键词]线性方程组,预处理,共轭梯度法中图分类号:E911 文献 ...

三对角方程组的追赶法

2013-2014(1)专业课程实践论文题目:三对角方程组的追赶法一.算法理论在一些实际问题中,例如解常微分方程边值问题,解热传导方程以及船体数学放样中建立三次样条函数等,都会要求解系数矩阵为对角占优的三对角线方程组 ⎛b 1c 1⎫⎛x 1⎫⎛ ⎪⎪ a b c x 2 22⎪2⎪ ⎪ ...

科学和工程计算复习题及答案

科学和工程计算基础复习题一. 填空题: 1. :2. 计算机计费的主要依据有两项:一是使用要由算数运算的次数决定; 二是占据存储器的空间, 3. 用计算机进行数值计算时, 4. 5. 6. 7. 8. 9. 10. 11. 敛的充分必要条件是选代矩阵B 的谱半径ρ(B ) (n +1) (x ...

解线性方程组的直接法

第三章解线性方程组的直接法许多科学技术问题要归结为解含有多个未知量x1, x2, „, xn的线性方程组 a11x1a12x2a1nxnb1axaxaxb2112222nn2   an1x1an2x2annxnbn (3.1) 这里aij (i, ...

深圳大学信息工程学院

深圳大学信息工程学院 <数值理论与技术方法>课程教学大纲一.课程基本信息课程编号:2313100201, 2313100202 课程名称:数值理论与技术方法课程类别:综合选修课适用专业:信息工程学院电子信息工程先修课程:高等数学.线性代数开课学期:2011-2012学年度第一 ...

滑差速器非圆锥齿轮线切割加工技术

第10卷第4期 2008年7月军事交通学院学报 JounrnalofAcademyofMilitaryTransportation V01．10No．4 July2008 ●基础理论与应用BasicTheoryandItsApplication 自适应变速比防滑差速器非圆锥齿轮线切割加工技术贾巨 ...