2-2 LU crout 分解

#include <math.h>
#include "..\CommonFiles\nrutil.h"
#define TINY 1.0e-20;/* Crout算法
** indx为输出向量，保存部分主元法而改变了行的行排列顺序
** 输出向量d为±1，表示行交换次数为偶数还是为奇数
*/
void ludcmp(float **a, int n, int *indx, float *d) {int i,imax,j,k;float big,dum,sum,temp;float *vv;              /* 用于保存每行的内含比例因子*/vv=vector(1,n);*d=1.0;for (i=1; i<=n; i++) {  /* 按行循环，求内含比例因子(每一行最大元素的倒数)*/big=0.0;for (j=1; j<=n; j++)if ((temp=fabs(a[i][j])) > big) big=temp;if (big == 0.0) nrerror("Singular matrix in routine ludcmp");vv[i]=1.0/big;}for (j=1; j<=n; j++) {            /* 外层列循环*/for (i=1; i<j; i++) {        /* 上三角部分*/sum=a[i][j];for (k=1; k<i; k++) sum -= a[i][k]*a[k][j];a[i][j]=sum;}big=0.0;                       /* 初始化最大主元*/for (i=j; i<=n; i++) {      /* i >= j 下三角部分*/sum=a[i][j];for (k=1; k<j; k++)sum -= a[i][k]*a[k][j];a[i][j]=sum;if ( (dum=vv[i]*fabs(sum)) >= big) {big=dum;imax=i;}}if (j != imax) {            /* 是否进行交换？*/for (k=1; k<=n; k++) {dum=a[imax][k];a[imax][k]=a[j][k];a[j][k]=dum;}*d = -(*d);             /* 改变d的奇偶性*/vv[imax]=vv[j];}indx[j]=imax;if (a[j][j] == 0.0) a[j][j]=TINY;if (j != n) {dum=1.0/(a[j][j]);for (i=j+1; i<=n; i++) a[i][j] *= dum;}}free_vector(vv,1,n);
}
#undef TINY

转载于:https://www.cnblogs.com/dLong/p/4477945.html

2-2 LU crout 分解相关推荐

crout分解计算例题_矩阵与数值计算（2）——矩阵三角分解LU、PALU、Cholesky三角分解、QR分解...
前言矩阵分解是设计算法的主要技巧,通过分解可以将复杂问题转化为几个简单问题求解,通常完成这一转化任务的主要技巧就是矩阵分解.例如,我们知道上三角矩阵和下三角矩阵是容易求解的,或者对角矩阵是最理想的求 ...
线性代数——LU(LR)分解
文章目录定义为什么要LU分解为什么能做到LU分解利用LU分解求行列式值利用LU分解求解线性方程组利用LU分解求逆矩阵对角线上元素有0的情况定义定义:给定矩阵A,将A表示成下三角矩阵L ...
python --- 实现LU=PA分解（部分主元的高斯消去）/带行交换的LU分解
题设: Pascal矩阵进行 LU=PA分解 import numpy as np import random from scipy.linalg import lu#Вычисление знач ...
Cholesky和LU矩阵分解
1.Cholesky分解在线性代数中,矩阵分解是将矩阵分解为矩阵的乘积.有许多不同的矩阵分解.其中之一就是Cholesky分解. Cholesky分解是将Hermitian正定矩阵分解为下三角矩阵及 ...
运筹系列70：pinv/inv/LU/SVD分解
1. LU分解 LU分解(LU Factorization)可以将一个矩阵分解为一个单位下三角矩阵和一个上三角矩阵的乘积.LU分解非常简单,下图可以直观的看出LU分解的求解方法. LU分解主要应用在数 ...
crout分解计算例题_专题：化学方程式计算
一. 根据化学方程式的简单计算 (1) 根据化学方程式计算的依据化学方程式表达的信息之一是反应物与生成物之间在"遵循固定质量比"的前提下的质量守恒.根据这一信息,可以利用化学方程 ...
crout分解计算例题_如何计算有理函数的不定积分
本文中所有不定积分都省略常数项.不同小节的相同字母没有关系. 部分分式分解一般微积分书上都会讲有理函数的部分分式分解,以及如何用这种方法算有理函数的不定积分.方法如下:假设是给定的有理函数,其中 ...
crout分解计算例题_初中数学因式分解习题大全（建议父母给孩子收藏直接打印）...
一．填空题(共10小题) 1．已知x+y=10,xy=16,则x2y+xy2的值为． 2．两位同学将一个二次三项式分解因式,一位同学因看错了一次项系数而分解成2(x﹣1)(x﹣9):另一位同学因看错 ...
矩阵的三角分解(LU分解)
矩阵的三角分解将矩阵分解成一个下三角矩阵L和一个上三角矩阵U的乘积. 定义:如果n阶矩阵A能够分解成一个下三角矩阵L和一个上三角矩阵U的乘积,则称这种分解为三角分解或LU分解,如果n阶矩阵A能够分解为 ...

2-2 LU crout 分解

2-2 LU crout 分解相关推荐

最新文章

热门文章