/*************************************
* 用追赶法求线性方程组
*
* |- -| |- -| |- -|
* | 2 -1 | | x1 | | 3 |
* |-1 3 -2 | | x2 | | 1 |
* | -2 4 -2 |*| x3 | = | 0 |
* | -2 5 | | x4 | | -5 |
* |- -| |- -| |- -|
*
**************************************/
#include<stdio.h>
#include<math.h>
#include<conio.h>
#include<string.h>
#include<stdlib.h>
#define n 5
int main() {
float a[n] = { 0, 0, -1, -2, -2 };
float b[n] = { 0, 2, 3, 4, 5 };
float c[n] = { 0, -1, -2, -2, 0 };
float d[n] = { 0, 3, 1, 0, -5 };
float x[n] = { 0, 0, 0, 0, 0 };
float r[n] = { 0, 0, 0, 0, 0 };
float y[n] = { 0, 0, 0, 0, 0 };
float q;
int k;
r[1] = c[1] / b[1];
y[1] = d[1] / b[1];
for (k = 2; k < n - 1; k++) {
q = b[k] - r[k - 1] * a[k];
r[k] = c[k] / q;
y[k] = (d[k] - y[k - 1] * a[k]) / q;
}
y[n - 1] = (d[n - 1] - y[n - 2] * a[n - 1])
/ (b[n - 1] - r[n - 2] * a[n - 1]);
x[n - 1] = y[n - 1];
for (k = n - 2; k >= 1; k--)
x[k] = y[k] - r[k] * x[k + 1];
for (k = 1; k < n; k++)
printf("x[%d] = %f\n", k, x[k]);
return 0;
}
分享到:
相关推荐
MFC追赶法求线性方程组,个人编写小软件,结合了计算方法
解线性方程组的直接法是用有限次运算求出线性方程组 Ax=b 的解的方法。线性方程组的直接法主要有Gauss消元法及其变形、LU(如Doolittle、Crout方法等)分解法和一些求解特殊线性方程组的方法(如追赶法、LDLT法等)。...
结合一个题目演示了如何自编代码实现高斯消元法和追赶法,各函数文件独立,便于移植,题目附有解答,题目来自西工大数值计算方法作业。
○1高斯-赛德尔迭代法比雅克比的迭代次数少一半多,说明高斯-赛德尔迭代法优于雅克比迭代法,收敛速度更快,效率更高。但也不绝对,有时也会出现雅克比迭代法收敛而高斯-赛德尔迭代法不收敛的情况。 ○2SOR法的迭代...
C语言实现L U分解 追赶法 平方根法对线性方程组的求解 计算方法中算法的实现
一.试验目的:练习线性方程组的直接解法。 二.实验内容:求解如下三对角线方程组: 实验方法:追赶法
内容包括:线性方程组的直接解法:Gauss消去法、Gauss列主元消去法、Gauss全主元消去法、列主元消去法应用『列主元求逆矩阵、列主元求行列式、矩阵的三角分解』、LU分解法、平方根法、改进的平方根法、追赶法(解三对...
东北大学理学院数学系2017年12月13日计算数学教研室东北大学解线性方程组的直接方法东北大学解线性方程组的直接方法追赶法向量和矩阵的范数Gauss消去法平方根
非线性方程与方程组的求解195 4.1 求非线性方程实根的对分法195 4.2 求非...方法211 4.7 求实系数代数方程全部根的牛顿下山法216 4.8 求复系数代数方程全部根的牛顿下山法225 4.9 求非线性方程组一组实根的梯度法233 ...
Gauss消去法与矩阵三角分解法(Doolittle分解法相比Crout分解法更常用)及其选择列主元的改进方法、Doolittle分解法的延伸(实对称正定矩阵利用Cholesky分解得到的平方根法、三对角矩阵作为线性方程组系数矩阵的追赶...
实验二 线性方程组的直接解法 合理利用Gauss消元法、LU分解法、追赶法求解方程组: 实验三 线性方程组的迭代解法 线性方程组的迭代解法 使用雅可比迭代法或高斯-赛德尔迭代法对下列方程组进行求解。 实验四 矩阵...
现代数值计算方法pdf实验二 线性方程组数值求解 1 用列主元高斯消元法解方程组2.2(p26) 2 用追赶法解方程组2.4(p26) 3 分别用雅可比迭代法、高斯赛德尔迭代法求解线性方程组8.1(p131),自设迭代求解精度 4 用...
4.11 求非线性方程组最小二乘解的广义逆法 4.12 求非线性方程一个实根的蒙特卡洛法 4.13 求实函数或复函数方程的一个复根的蒙特卡洛法 4.14 求非线性方程组一组实根的蒙特卡洛法 第5章 插值 5.1 插值类...
实验一 直接法解线性方程组的 1 实验二 插值方法 4 实验三 数值积分 6 实验四 常微分方程的数值解 8 实验五 迭代法解线性方程组与非线性方程 10
计算方法-课件+实例程序源代码+误差集锦 ...【Lagrange 插值】,【Newton 插值】,【分段线性插值】,【复合梯形公式求定积分】,【列主元高斯】,【牛顿迭代】,【数据拟合】,【线性方程组迭代】,【追赶法】
线性方程组的求解,LU分解(高斯消元法)、正定矩阵的Cholesky分解、追赶法求解特殊结构的矩阵. 非线性方程的求解:牛顿迭代法、不动点迭代、二分法等等其他先进的算法。 函数逼近中学习的插值算法:拉格朗日插值、...
里面有十几个计算方法的C++程序实现,包括QR方法求矩阵特征值、列主元LU分解、全主元高斯消去法、列主元高斯消去法、消元形式的追赶法、一次求解多个线性方程组等。跟大家分享一下,希望能对大家有所帮助。
4.11 求非线性方程组最小二乘解的广义逆法246 4.12 求非线性方程一个实根的蒙特卡洛法262 4.13 求实函数或复函数方程一个复根的蒙特卡洛法265 4.14 求非线性方程组一组实根的蒙特卡洛法269 第5章 插值与逼近274 5.1 ...
7.10 求非线性方程组最小二乘解的广义逆法 7.11 求非线性方程一个实根的蒙特卡洛法 7.12 求实函数或复函数方程一个复根的蒙特卡洛法 7.13 求非线性方程组一组实根的蒙特卡洛法 第8章 插值与逼近 8.1 一元全区间插值 ...