数值分析计算实习作业二(共20页).docx

数 值 分 析计算实习题二学号：姓名： 院系：2015年11月28日数值分析计算实习作业二一、算法采用拟上三角化的双步位移的QR分解法计算特征值，并利用列主元素Guass消去法求实数特征值对应的特征向量，定义全局变量n=10。1 主程序在定义一维和二维数组时将数组的大小定为n+1（大小为11），这样在使用单个元素时的角标可以与矩阵或向量元素的角标相对应。先输入矩阵A的各个元素；保留数组A，定义数组B的元素与A相同，即B=A，后面的矩阵变换用数组B来进行；对B进行拟上三角化；输出拟上三角化后的矩阵B；对拟上三角化后的B进行QR分解，输出矩阵Q、R、RQ；对拟上三角化的数组矩阵B进行QR分解；输出A的全部特征值；判断所有特征值，如果特征值为实数（即储存特征值虚部的数组中元素为0），利用列主元素高斯消去法求其特征向量；输出实特征值对应的特征向量。2 拟上三角化程序对于r=1,2，n-2循环执行（1）求第r列的第r+2个元素到第n个元素