数值分析第二题 北航 大作业.doc

《数值分析》计算实习题目二 PAGE 32 数值分析第二题 目 录 TOC \o 1-3 \h \z \u 数值分析第二题 1 1 引言： 2 1.1 矩阵的拟上三角化 2 1.2 矩阵的特征值求解 3 1.3 矩阵的特征向量求解 5 2 算法的程序实现 6 2.1 主程序 6 2.2 子程序的实现 8 3计算结果 9 3.2矩阵Q、R 以及乘积RQ 10 3.3 各实特征值及其相对应的特征向量 11 4 实验结论 12 附录 源代码 13 1 引言 1.1 矩阵的拟上三角化 为了减少计算量，对矩阵利用Householder矩阵进行相似变换，把化为上三角矩阵A(n-1)。 对拟上三角化，得到拟上三角矩阵，具体算法如下： 记，并记的第r列至第n列的元素为。 对于执行 若全为零，则令,转5；否则转2。 计算 令。 计算 继续。 1.2 矩阵的特征值求解 使用带双步位移的QR方法计算矩阵的全部特征值，也是的全部特征值，具体算法如下： 给定精度水平和迭代最大次数。 记，令。 如果，则得到的一个特征值，置（降阶），转4；否则转5。 如果，则得到的一个特征值，转11；如果，则转3。 求2阶子阵 的两个特征值和，即计算二次方程 的两个根和。 如果，则得到的两个特征值和，转11；否则转7。 如果，则得到的两个特征值和，置（降阶），转4；否则转8 如果，则计算终止，未得到的全部特征值；否则转9。 记，计算 置，转3。 的全部特征值已计算完毕，停止计算。 其中，的分解与的计算用下列算法实现： 记。 对于执行 若全为零，则令,转5；否则转2。 计算 令。 计算 继续。 此算法执行完后，就得到。 1.3 矩阵的特征向量求解 用列主元素Gauss消去法计算矩阵对应于实特征值的特征向量，具体算法如下： 记 消元过程 对于执行 选行号，使。 交换与所含的数值。 对于计算 回代过程 最终得到的向量的即为对应于实特征值的特征向量。 2 算法的程序实现 2.1 主程序 #include iostream #includematrix-tool.h #includesimi-trian.h #includeQR_decompose.h #includeGauss_by_max.h #define threhold (1e-12) using namespace std; int n = 10; extern double *a; extern double *Mk; extern double *a_lambda; extern double *b_lambda; double *spec_u; //高斯求解方程组的右端常数列向量 double *G_b; int main() { simi_Initial(); cout处理之前的A矩阵为：endl; Print_Matrix(a,n); cout拟上三角化后的A矩阵为：endl; pseudo_triangle(a,n); BI_STEP_QR(a,n); /******************特征向量求取代码段*****************/ //初始化spec_u[10] spec_u = new double [n]; G_b = new double [n]; char i; simi_Initial(); //为G_b 赋值 for(i=0; in; i++) { (*(G_b+i)) = 0; } for(i=0; in; i++) { if( (*(b_lambda+i)) == 0 ) { Gauss_root( a, G_b, n, (*(a_lambda+i)),spec_u); cout属于(*(a_lambda+i)) i (*(b_lambda+i))的特征向量为： endl; Print_Vector(spec_u,n); } } system(pause); return 0; } 2.2 子程序的实现 2.2.1 “matrix-tool.h” void Matrix_M_Vector(double *A,double *b,double *y,int n); 函数实现 N*N的矩阵乘以N*1的向量，得到N*1向量Y void Matrix_Converted_M_Vector(double *A,double *b,double *y,int n); 函数实现转置后N*N矩阵乘以N*1向量,得到N*1向量Y vo