Matrix说明文档矩阵说明文档.docVIP

Matrix说明文档 功能说明：解实数范围内的超定方程组。 使用说明：在与Matrix.exe的同一目录下建立文本文件, a.txt 和b.txt。a.txt中存放方程组的系数矩阵，b.txt中存放方程组的常数项矩阵。点击Matrix.exe运行结果如下： 在“矩阵行数”中输入方程组系数矩阵的行数，在“矩阵列数”中输入方程组系数矩阵的列数，然后点击“读入系数矩阵”和“读入常数项矩阵”按钮，最后点击“计算”按钮，即可得到结果。 用以下超定方程组作为测试，方程组如下：2X1- X2=1; 8X1+4X2=0; 2X1+X2=1; 7X1- X2=8; 4X1 =3。已知结果为：X1=0.793，X2=-1.464。与用Matrix.exe计算出的结果符合。 操作步骤及代码说明： 首先用MFC建立基于对话框的应用程序，取工程名为“Matrix”。 用“建立类向导”中为“矩阵行数”编辑框和“矩阵列数”编辑框分别添加成员变量int m_h和int m_l。 m_h中包含了系数矩阵的行数，m_l包含了系数矩阵的列数。 在MatrixDlg.cpp中添加代码： #include fstream using namespace std; const int Max=10000; bool control=0; float **a;//系数矩阵 float **b;//常数项矩阵 float **at;//a的转置 float **ata;//at与a的乘积 float **atb;//at与b的乘积 float **u; float* fOneDPointer(int num) //对一维矩阵初始化 { register int i; float* p = new float[num]; for(i=0; inum; i++){ p[i]=0.0f; } return(p); } void Free_fOneDPointer(float *p) //释放一维矩阵占用的空间 { delete [] p; } float** fTwoDPointer(int raw,int col) //对二维矩阵初始化 { register int i,j; float* arr = new float[raw*col]; float** p = new float*[raw]; for(i=0; iraw; i++){ p[i] = arr + i*col; for(j=0; jcol; j++) p[i][j]=(float)0.0; } return(p); } void Free_fTwoDPointer(float **p) //释放二维矩阵占用的空间 { delete [] *p; delete [] p; } void CMatrixDlg::InitMatrix(int m,int n)//对各个二维数组进行初始化 { UpdateData(true); if(control==0){ m=m_h; n=m_l; a=fTwoDPointer(m,n); b=fTwoDPointer(m,1); at=fTwoDPointer(n,m); ata=fTwoDPointer(n,n); atb=fTwoDPointer(n,1); u=fTwoDPointer(n,n+1); control=1; } } void CMatrixDlg::FreeMatrix()//释放各个二维数组的空间 { Free_fTwoDPointer(a); Free_fTwoDPointer(b); Free_fTwoDPointer(at); Free_fTwoDPointer(ata); Free_fTwoDPointer(atb); Free_fTwoDPointer(u); } void CMatrixDlg::MatrixInver(float **a,int m,int n) //求系数矩阵的转置 { UpdateData(true); int i,j; m=m_h; n=m_l; for(i=0;im;i++) for(j=0;jn;j++) { at[j][i]=a[i][j]; } } void CMatrixDlg::MatrixMultiple(float **at,float **a,int m,int n) //用于求系数矩阵的转置和系数矩阵的乘积 { UpdateData(true); int i,j,k; m=m_h