数值代数实验题目一.doc

数值代数实验题目一：矩阵分解的实现 实验原理： 定义：设非奇异，若有正交矩阵及上三角矩阵满足，则称为的分解。 可以利用Givens变换和Householder变换求非奇异矩阵的分解，分别称为Givens变换和Householder变换。 此外还有用经典Gram-Schmidt正交化方法（CGS方法）和修正Gram-Schmidt正交化方法（MGS方法）求非奇异矩阵的分解。 实验的目的和意义： 1．了解和掌握经典Gram-Schmidt正交化方法（CGS方法）和修正Gram-Schmidt正交化方法（MGS方法）。 2．了解和掌握Givens变换和Householder变换。 3．对于经典Gram-Schmidt正交化方法（CGS方法）和修正Gram-Schmidt正交化方法（MGS方法）充分理解其数值稳定性方面存在的问题。 4．对于Givens变换和Householder变换充分理解其程序设计方面存在的问题（关键点的程序处理）。 Givens变换的实际计算原则： （1）如何计算的问题。 为了尽可能减少运算次数，避免可能出现溢出现象，并考虑到算法的稳定性，对给定的，将按下面的基本算法来计算，使得的第个分量为零。 （1）给定； （2）如果，则； （3）如果，则 ，； 否则 ，。 （2）如何计算（）的问题。 注意到只改变的第行的元素，其余不变。如果知道，计算结果就可以存放在的存储单元。 Householder变换的实际计算原则： （1）前面的符号如何选取的问题： 如果是一个很接近的向量（只考虑方向，不考虑长度），则就很接近于零向量，从而单位化时会产生较大的误差。因此，应选择前面的符号与的第一个分量的符号相同，即 （2）如何避免计算时误差的问题。 由于对任意非零常数，向量与的单位化向量是一致的，因此可以选择代替来构造。另外，实际计算时也可以不用进行单位化：由于 其中。只需要给出即可。 （3）如何计算（）的问题。 只需注意到即可。 算法的设计与程序： 请同学自行完成！！！ 计算问题： 1．验证课堂例题的计算结果的正确性。 2．Hilbert矩阵 定义：设，其中 根据计算机的性能，分别计算时矩阵的分解，注意四种方法给出的结果，试分析产生如此结果的原因。 完成实验报告（主要写出对算法的认识和对计算结果的分析）。