求解线性方程组Ax=b。手动输入A,其中A为n*n的N阶矩阵,x为列向量,即需要求解的未知的变量,b为列向量,是非齐次线性的常数
求解N阶线性矩阵方程
一:编程语言: C语言
二:问题描述
求解线性方程组Ax=b。手动输入A,其中A为n*n的N阶矩阵,x为列向量,即需要求解的未知的变量,b为列向量,是非齐次线性的常数。 三:问题分析及解题思路
首先,用户输入矩阵A和列向量b,这些用C语言很容易实现,然后判断方程是否有解,及解的个数。运用数学知识知道矩阵A的行列式不等于0才有唯一解。然后进行初等变换,使系数转换为上三角矩阵,这些代码第一次作业做过了,很容易实现。再然后根据数学知识代入即可求得原方程唯一解。
四:代码实现
#include"stdio.h"
#include"stdlib.h"
#include"math.h"
#define dim 10
double a[dim+1][dim+1],b[dim+1],x[dim+1];
double temp;
double getarray(int n);