东南大学C语言课程设计报告PPT.pptVIP

题目 n阶方阵求逆(B级) 功能要求： 输入一个n（n256）阶方阵A，方程系数矩阵与向量均从文本文件读入； 输出A的逆矩阵至文本文件； 将得到的逆矩阵与矩阵A相乘，验证其结果是否为单位矩阵。 提示：具体算法可参考相关文献。 题目分析 逆矩阵定义： 设A为n阶方阵，如果存在n阶方阵B，使得 AB=BA=E， 则称A是可逆的，并称B是A的逆矩阵。否 则，便说A是不可逆的。 题目分析 矩阵乘法： 若A是一个m*n阶矩阵，B是一个n*p阶矩 阵，则AB=C是一个m*p阶矩阵，而C中的每一 个（i,j）元都等于A的第i行中的各元和B的 第j列的各对应元之乘积的和。 算法设计 在数学问题中，求解方阵逆矩阵的方法有很多。其中比较常见的一种方法是用初等变换求方阵的逆矩阵，即高斯消元法求逆矩阵。高斯消元法可以用来找出一个可逆矩阵的逆矩阵。设A 为一个N * N的矩阵，其逆矩阵可被两个分块矩阵表示出来。将一个N * N单位矩阵 放在A 的右手边，形成一个N * 2N的分块矩阵B = [A,I] 。经过高斯消元法的计算程序后，矩阵B 的左手边会变成一个单位矩阵I ，而逆矩阵A - 1 会出现在B 的右手边。 //采用高斯消元法 #include stdio.h #include string.h #define maxn 1000 #define err 1e-8//宏定义保证计算中准确性 double a[maxn][maxn*3];//maxn是代码最大元素的下标 double abs(double a)//返回传递给它的参数的绝对值函数 { if (a0) return a; else return -a; }//定义绝对值函数 运行结果 不存在逆矩阵：不存在逆矩阵的情形有很多种，可能是初始矩阵的某行全为零，也有可能是初始矩阵经过变化之后其某行全部变为零。此时，无答案输出至文本文件。 2.存在逆矩阵的情形，这里选取了一6阶方阵，如图所示为一6阶方阵，在VC中运行运行后，我们可以得到如下的结果 * 姓 名 汤佳敏 学 号 一卡通号 213112153 ——n阶方阵求逆矩阵 开始 从文本读入系数矩阵及向量 将矩阵变成增广矩阵 找到一行中不为零的数 是 否 跳出当前循环 该矩阵为奇异矩阵 单位化矩阵 消元 直到循环玩每一行 输出逆矩阵至文本文件 逆矩阵与原矩阵相乘 结束 如图所示是程序的算法流程图。从中我们可以看出该程序主要分为七部分，分别为：文件的读入部分、判断是否为逆矩阵部分、单位化矩阵部分、消元部分、打印步骤部分、检验乘积部分及输出部分。其中单位化和消元部分是本程序设计的核心，即高斯消元法的主要步骤。具体程序如下： 在运算的过程中，我们要判断 一些数据是否为零，这时候我们采 用宏定义#define err 1e-8方式，通 过将得出数据在绝对值与e-8进行比 较，可以保证我们计算的准确性。 int main() { printf(Now read the matrix from m.txt..\n); int n,i,j,k; double l; FILE * f1 =fopen(m.txt,r); fscanf(f1,%d,n);//从文本文件读入方程系数矩阵 memset(a,0,sizeof(a));//数组a初始化 for (i=0;in;i++) for (j=0;jn;j++) { fscanf(f1,%lf,a[i][j]);//从文本读入矩阵 a[i][j+2*n]=a[i][j];//在n*2n阶矩阵后再加上原矩阵形成n*3n阶矩阵 } fclose(f1); printf(Done reading.\n); 从文本文件分别读入方程系数矩阵与向量 for (i=0;in;i++) a[i][i+n]=1;//建立单位矩阵，即将单位矩阵变成伴随矩阵 for (i=0;in;i++) { bool got=false;//定义布尔型变量 for (j=i;jn;j++) if (abs(a[j][i])err) { got=true; for (k=0;kn*2;k++)