北航数值分析大作业-第二题-QR分解.docxVIP

北航数值分析大作业-第二题-QR分解 北航数值分析大作业-第二题-QR分解 PAGE 北航数值分析大作业-第二题-QR分解 《数值分析B》课计算实习第一题 设计文档与源程序 姓名：杨彦杰 学号：SY1017134  1 算法的设计方案 （1）运行平台 操作系统：Windows XP； 开发平台：++； 工程类型：文档视图类； 工程名：Numanalysis； （2）开发描述 首先新建类CMetrix，该类完成矩阵之间的相关运算，包括相乘、加减等，以主程序方便调用； 题目的解算过程在视图类CNumanalysisView中实现，解算结果在视图界面中显示； （3）运行流程 计算开始 计算开始 动态生成内存空间 矩阵A的求解 矩阵A拟上三角化 严格按照解算步骤求解矩阵A的全部特征值 求解实根的特征向量 求解结果显示 关闭内存空间 计算结束 （4）运行界面 2、全部源代码 （1）类CMetrix 文件： class CMetrix { public: double** MetrixMultiplyConst(double**A,int nRow,int nCol,double nConst); e, cr = %.6e, hr = %.6e,dr,cr,hr); for (int row = 0;row nRow;row++) ealnum = A[m-1][m-1]; eigenvalue[m-1].Imagnum = ; m = m - 1; } if (m == 1) { eigenvalue[m-1].Realnum = A[m-1][m-1]; eigenvalue[m-1].Imagnum = ; DisplayText(已求出A的全部特征值：); break; } b = -(A[m-2][m-2]+A[m-1][m-1]); e + j%.6e,eigenvalue[i].Realnum,eigenvalue[i].Imagnum); DisplayText(str); } magnum != { (特征值%.6e+j%.6e为虚数，不需要求特征向量。,eigenvalue[i].Realnum, eigenvalue[i].Imagnum); DisplayText(str); continue; } tempA = (I,m,m,eigenvalue[i].Realnum); for (row = 0;row m;row++) { for (col = 0;col m;col++)Ar[row][col] = tempA[row][col]; } tempA = (A,Ar,m,m); e 的特征向量为：,eigenvalue[i].Realnum); DisplayText(str); DisplayVector(tempV,nRow); } delete[] A,Ar,tempA,Cr,C,I; delete V,ur,vr,pr,qr,wr,tempV,eigenvalue; } void CNumanalysisView::DisplayText(CString str) { int Len = GetWindowTextLength(); if(Len) str += \r\n; GetEditCtrl().SetSel(Len,Len); GetEditCtrl().ReplaceSel(str); } void CNumanalysisView::DisplayMetrix(double **A, int nRow, int nCol) { CString str,strtemp; for (int i = 0;i nRow;i++) { for (int j = 0;j nCol;j++) { (%.6e, ,A[i][j]); strtemp += str; } strtemp += \r\n; } DisplayText(strtemp); } void CNumanalysisView::DisplayVector(double *V, int nV) { CString str,strtemp; for (int i = 0;i nV;i++) { (%.6e, ,V[i]); strtemp += str; } strtemp += \r\n; DisplayText(strtemp); } 3拟上三角化矩阵A（n-1） 矩阵A拟上三角化后所得的矩