C实验报告高斯消元法.docxVIP

高斯肖元法 C++ 上机实验报告 学生姓名： 学 号： 专业班级： 实验类型： 综合 — 实验项目名称 全选主元高斯消去法解线性方程组二 实验原理 设有 n 元线性方程组(考虑便于 C++程序数组表示,方程的下标从 0 开始), ? a x ? 00 0 ? a x ? 00 0 ?? a x 01 1 a x 01 1 ? ? ? ? a ? ? ? ? a  0 , n ?1 0 , n ?1 x ? b n ?1 0 x ? b n ?1 1 ? a x ? a x ? ? ? ? a x ? b ? n1 0 n ?1,1 1 n ?1, n ?1 n ?1 n ?1 写为矩阵形式为 Ax=b,其中 A 为线性方程组的系数矩阵,x 为列向量,是方程组的解,b 也是列向量. 一般来讲，可以假定矩阵 A 是非奇异阵。（n 阶矩阵 A 的行列式不为零，即 |A|≠，0 则称 A 为非奇异矩 阵） ? a ? 00 a a ?01 0 , n ? ? ? x ??? ? 0 ? ? ? ? b ? ? ? 0 ? A ? ? a10 a a 11 1, n x ? ? x1 ? b ? ? b ? 1? ? ? ? 1 ? ? ， ? ? ? ， ? ? ? a? ? ? ? ? ? a ?? an ?1,0 ? a n ?1,1 n ?1, n ?1 ? ? xn ?1 ? ?b n ?1 ? 将系数矩阵 A 和向量 b 放在一起，形成增广矩阵 B： ? a ? 00 b ? ( A, b ) ? ? a10 a a b ? ?01 0 , n 0 ? ? ??a a b ? ? 11 1, n 1 ? ? ? ?a ? n ?1,0 ? ? ?a ? n ?1,1 ? a n ?1, n ?1 ? ? ? bn ?1 ? 全选主元消去就在矩阵 B 上进行，整个过程分为如下两个步骤： 第一步：消去过程。 对于 k 从 0 开始到 n-2 结束，进行以下三步。 首先，从系数矩阵 A 的k 行k 列开始的子矩阵中选取绝对值最大的元素作为主元素。例如： a i , j 1 1 ? max a ? 0 k ? i ? n ij k ? j ? n 然后交换 B 的第 k 行与第i 1 行，第 k 列与第 k 1 列，这样，这个子矩阵中具有最大绝对 值的元素被交换到 k 行 k 列的位置上. ?? a  ? a / a  , j ? k ? 1,? , n ? 1 ?其次，进行归一化计算。计算方法为： ? kj ? ? b k kj kk ? b / a k kk ?? a 最后进行消去计算： ? ij ? a ? a a ij ik kj , j , i ? k ? 1,? , n ? 1 ?? b ? b ? a b , i ? k ? 1,? , n ? 1 ? x ? n ?1 i ? b n ?1 i ik k / a n ?1, n ?1 ? ?第二步，回带过程 ? ? ?? x b ? n ?1 ? a x , i ? n ? 2, ? ,1, 0 三 代码的实现 i i ij j j ? i ?1 整个程序分为 5 个独立文件，Matrix.h 文件中包括矩阵类 Matrix 的定义，Matrix.cpp 文件中包括该类成员函数的实现，LinearEqu.h 文件中包括线性方程组类 LinearEqu 的定义，LinearEqu.cpp 文件中包括该类的成员函数实现文件；7-9.cpp 文件包括程序的主函数，主函数中定义了一个类 LinearEqu 的对象，通过这个对象求解一个四元线性方程组。 //Matrix.h 文件一，定义一个CMatrix 类 #ifndef _MATRIX_H #define _MATRIX_H class CMatrix //定义CMatrix基类 { public: //外部接口 CMatrix(int size=2); //构造函数 ~CMatrix(); //析构函数 void setMatrix(const double*values); //矩阵赋初值void printMatrix() const; //显示矩阵 int getsize() const{return size;} //得到矩阵大小 double element(int i,int j){return elements[i*size+j];} double element(int i,int j)const{return elements[i*size+j];} private: //保护数据成员 int size; //定