LUP分解的C++实现

简单粗暴上代码

在这里插入代码片#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream>
using namespace std;void Swap(double* a, double* b) {//交换函数 double temp = *a;*a = *b;*b = temp;
}
void Swap1(int* a, int* b) {//交换函数 int temp = *a;*a = *b;*b = temp;
}double Round(double dVal, short iPlaces) {//四舍五入保留小数位数 double dRetval;double dMod = 0.0000001;if (dVal < 0.0) dMod = -0.0000001;dRetval = dVal;dRetval += (5.0 / pow(10.0, iPlaces + 1.0));dRetval *= pow(10.0, iPlaces);dRetval = floor(dRetval + dMod);dRetval /= pow(10.0, iPlaces);return(dRetval);
}void Input(double** A,  int n) {int i, j;cout << "请输入A矩阵:" << endl;for (i = 1; i <= n; i++) {for (j = 1; j <= n; j++) {cin >> A[i][j];}}/*cout << "请输入b数组:" << endl;for (i = 1; i <= n; i++) {cin >> b[i];}*/
}void LUP_DECOMPOSITION(double** A, double** L, double** U, int* P, int n) {int i, j, k, k1;double temp, temp1;;for (i = 1; i <= n; i++) {P[i] = i;}for (k = 1; k <= n; k++) {int p = 0;for (i = k; i <= n; i++) {if (fabs(A[i][k]) > p) {p = fabs(A[i][k]);k1 = i;}}if (p == 0) {cout << "singular matrix" << endl;}else {Swap1(&P[k], &P[k1]);}for (i = 1; i <= n; i++) {Swap(&A[k][i], &A[k1][i]);}for (i = k + 1; i <= n; i++) {temp = A[i][k] / A[k][k];A[i][k] = Round(temp, 2);for (j = k + 1; j <= n; j++) {temp = A[i][k] * A[k][j];temp1 = Round(temp, 2);A[i][j] = A[i][j] - temp1;}}}for (i = 1; i <= n; i++) {for (j = 1; j <= n; j++) {if (i == j) {L[i][j] = 1;}else {if (i < j) {L[i][j] = 0;}else {L[i][j] = A[i][j];}}}}for (i = 1; i <= n; i++) {for (j = 1; j <= n; j++) {if (i <= j) {U[i][j] = A[i][j];}else {U[i][j] = 0;}}}for (i = 1; i <= n; i++) {P[i] -= 1;}}void Print(int* x, int n) {int i, j;for (i = 1; i <= n; i++) {cout << x[i] << "\t";}
}void PrintL(double** L, int n) {int i, j;for (i = 1; i <= n; i++) {for (j = 1; j <= n; j++) {cout << L[i][j] << "\t";}cout << endl;}
}int main() {int n;double* x = NULL;double* y = NULL;double** L = NULL;double** U = NULL;int* P = NULL;double** A = NULL;double* b = NULL;cin >> n;int i = 0;L = new double* [n + 1];for (i = 1; i <= n; i++) {L[i] = new double[n + 1];}U = new double* [n + 1];for (i = 1; i <= n; i++) {U[i] = new double[n + 1];}A = new double* [n + 1];for (i = 1; i <= n; i++) {A[i] = new double[n + 1];}P = new int[n ];x = new double[n + 1];y = new double[n + 1];b = new double[n + 1];Input(A,  n);LUP_DECOMPOSITION(A, L, U, P, n);//LUP_SOLVE(A, L, U, P, x, y, b, n);//Print(x, n);cout << "L:" << endl;PrintL(L, n);cout << "U:" << endl;PrintL(U, n);//PrintL(A, n);cout << "P:" << endl;Print(P, n);return 0;
}/*1 0 0 0 1 0 1 0
0 0 0 0 1 0 1 1
0 1 0 0 0 1 0 0
0 0 0 0 1 1 1 1
0 0 1 0 1 0 0 1
0 0 0 0 0 1 1 0
0 0 0 1 0 1 0 0
0 0 0 0 0 0 1 11 0 0 0 1 1 0 0
1 1 0 0 1 0 1 0
0 1 1 0 0 1 0 0
1 0 1 1 1 1 1 1
1 1 0 1 0 0 1 0
0 1 1 0 1 0 0 0
0 0 1 1 0 1 0 0
0 0 0 1 1 0 1 0
*/

如果对您有用请点赞哦，欢迎评论交流

LUP分解的C++实现相关推荐

矩阵分解及其Eigen实现
主要是用来记录自己的学习过程,内容也主要来自于网上的各种资料,然后自己总结而来,参考的资料都以注明,感谢这些作者的分享.如果内容有误,请大家指点. LU分解1 理论定义将矩阵等价为两 ...
matlab lud矩阵分解,MIT线性代数总结笔记——LU分解
MIT线性代数总结笔记--LU分解矩阵分解矩阵分解(Matrix Factorizations)就是将一个矩阵用两个以上的矩阵相乘的等式来表达.而矩阵乘法涉及到数据的合成(即将两个或多个线性变换的 ...
LU分解完成利用节点电压法的简易电路求解程序（二）
上一篇博客中我们回顾了LU分解,分析了一个简单的电路并观察到电导矩阵中的对角线项为节点相连的电导之和,非对角线项等于两个节点之间电导的相反数与电流项为独立电流源流入节点与流出节点电流代数和,其中 ...
ACM 全部算法总结
ACM 所有算法附带我学过的算法的博客链接数据结构栈队列链表哈希表哈希数组堆优先队列双端队列可并堆左偏堆二叉查找树 Treap 伸展树并查集集合计数问题二分图的识别平 ...
ACM基础知识及算法
ACM 算法难度数据结构栈栈 1 单调栈队列一般队列 1 优先队列/单调队列 1 循环队列 2 双端队列 2 链表一般链表 ...
你是如何坚持读完《算法导论》这本书的？(帖子收集)
你是如何坚持读完<算法导论>这本书的? <算法导论>不够猛,答者顺便补充 "你是如何坚持读完<计算机编程的艺术>这本书的?" 罗必成: CLRS ...
ACM所有算法大全（持续更新）
转载自: http://blog.sina.com.cn/s/blog_adb6743801019h29.html ACM 所有算法数据结构栈,队列,链表哈希表,哈希数组堆,优先队列双端队列 ...
【心随意动】20-SIST研一秋季上学期课程总结【补】
呀,自从来到了SIST,每天每夜都要么在写作业,要么就是在去写作业的路上,每天都是它们的deadline呐!好不容易有了一个早上可以睡觉养精神,我也是用来刷手机或发会儿呆.以至于上学期期间完全没有时间 ...
基于C++的简易RLC电路仿真器与滤波器仿真测试
电路建模分析与仿真算法对于任意一个电路系统,我们基本上都能使用RLC和受控源进行等效,因此RLC系统的求解是任意一个系统求解的核心. 我们熟知,对于任意一个复杂的电路系统,基尔霍夫定律(Kirchh ...
MATLAB处理矩阵的一些命令
矩阵的表示 A=[10 8 6 4 1;2 5 8 9 4;6 0 9 9 8] 矩阵中的元素用空格隔开,或用逗号隔开.不同行用分号隔开或用回车键换行. 2.命令 det(A) %计算行列式A的值 [ ...

LUP分解的C++实现

LUP分解的C++实现

LUP分解的C++实现相关推荐

最新文章

热门文章