平方根法乔累斯基分解Cholesky_解线性方程组的直接解法

标签：计算方法实验

#include <stdio.h>
#include <math.h>const int maxn = 15;int main(){double a[maxn][maxn], b[maxn], y[maxn], x[maxn], l[maxn][maxn];int n, sum;freopen("sqrt.txt", "r", stdin);scanf("%d", &n);for(int i = 1; i <= n; i++){for(int j = 1; j <= n; j++)  scanf("%lf", &a[i][j]);scanf("%lf", &b[i]);}/*for(int i = 1; i <= n; i++){for(int j = 1; j <= n; j++)  printf("%10f", a[i][j]);printf("%10f\n", b[i]);}*/for(int i = 1; i <= n; i++){  //算L: 对正定矩阵A分解，A = LL^Tsum = 0;for(int k = 1; k <= i - 1; k++)  sum += (l[i][k] * l[i][k]);l[i][i] = sqrt(a[i][i] - sum);sum = 0;for(int j = i + 1; j <= n; j++){for(int k = 1; k <= i - 1; k++)  sum += (l[j][k] * l[i][k]);l[j][i] = (a[j][i] - sum) / l[i][i];}}for(int i = 1; i <= n; i++){  //求y: A = LL^T -> LL^Tx = b -> Ly = bsum = 0;for(int k = 1; k <= i - 1; k++)  sum += (l[i][k] * y[k]);y[i] = (b[i] - sum) / l[i][i];}for(int i = n; i >= 1; i--){  //求x: L^Tx = ysum = 0;for(int k = i + 1; k <= n; k++)  sum += (l[k][i] * x[k]);x[i] = (y[i] - sum) / l[i][i];}for(int i = 1; i <= n; i++)  printf("%10f\n", x[i]);return 0;
}

数据文件

实验结果

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法相关推荐

乔利斯基三角分解_解线性方程组的直接法4.1-2.ppt
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp高等教育&nbsp>&nbsp微积分解线性方程组的直接法4.1-2.ppt24页本文档一 ...
[矩阵的三角分解系列四] 乔累斯基(Cholesky)分解公式
乔累斯基分解公式简介 LLT分解证明具体解法稳定性 LDLT分解证明具体解法例子 LLT分解 LDLT分解引用矩阵的三角分解是求解线性方程组常用的方法,包括LU分解,LDU分解,杜利 ...
数值计算（一）之解线性方程组（高斯消去法，列选主元消去法，全选主元消去法，杜立特尔分解，克洛特分解，乔里斯基分解）
解线性方程组即解一个多元一次方程组,例如目录消去法分解法消去法原理没有学过高级的解法也没关系,凭借我们初高中的知识足以解决这个问题这是一个多元一次方程组,拥有n个未知量,也有n方程我们 ...
乔列斯基(Cholesky)法解方程(python,数值积分)
第四课乔列斯基(Cholesky)法解方程首先要清楚二次型和正定矩阵 "二次型"可以定义为n个变量的二次表达式如果这个二次型的所有变量X的值都等于或大于零,那么这个二次型就是 ...
乔里斯基分解（Cholesky decomposition）
在线性代数中,乔里斯基分解是将一个正定Hermite矩阵分解成为一个下三角阵和它的共轭转置阵的乘积. 如果矩阵A是正定Hermite阵,那么矩阵A可以做如下分解:其中L是一个下三角矩阵且主对角线元素严 ...
newmat乔里斯基分解
乔里斯基(cholesky)分解是将一个正定的厄米特(hermite)矩阵分解为一个下三角矩阵和它的共轭转置矩阵. 厄米特矩阵(共轭矩阵),矩阵对称的两个数字互相共轭.实数的共轭是自己本身,复数的共轭 ...
Cholesky分解、乔列斯基分解
一.简介 1.1 定理 Cholesky分解法又叫平方根法,是一种分解正定Hermite矩阵 (即 A=AH\boldsymbol A = \boldsymbol A^\mathrm HA=AH ...
参数随机场，随机参数生成python代码，基于乔列斯基分解中点法分解
本人用python较多,但是matlib使用很少,没装在这个软件,故采用python语言重现编写,里面没有真实的单元坐标,只能假设三个单元进行测试,如存在错误请及时和我联系,平附上对应matlib代码 ...
乔利斯基三角分解_二维二相模型的不完全乔利斯基——共轭斜量解法
前健罗厂刀二维二相模型是我国目前油田开发中使用最广泛的一种数值模型.这不仅因为它比较经济实用,对电子计算机的容量和速度要求不高,适宜于大规模普及推广应用,而且也因为现在我国油田的开发多采取注水保持压力 ...

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法

标签：计算方法实验

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法相关推荐

最新文章

热门文章

平方根法 乔累斯基分解Cholesky_解线性方程组的直接解法

平方根法 乔累斯基分解Cholesky_解线性方程组的直接解法

标签：计算方法实验

平方根法 乔累斯基分解Cholesky_解线性方程组的直接解法相关推荐

最新文章

热门文章

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法

平方根法乔累斯基分解Cholesky_解线性方程组的直接解法相关推荐