北航数值分析大作业二纯高分版.docxVIP

目录 题目： 2 算法设计思路和方案 2 关于第一步矩阵 A 的拟上三角化 2 关于对矩阵 ????- ??进行带双步位移的 QR 分解迭代 3 关于求从属于矩阵 A 的实特征值 ????的特征向量 3 计算结果 4 发现的现象与问题 : 7 探究带双步位移的 QR 分解比一般 QR 分解节省的计算量 7 探究拟上三角化对 QR 分解迭代收敛速度的影响 7 关于直接单步 QR 分解计算发现的问题 9 源程序 10 数值分析上机实习作业二 王强 SY1413315 题目： 试求矩阵 ??= [?? ] × 的全部特征值，并对其中的每个实特征值求相应的特征向量， ????10 10 已知： sin ( 0.5??+ 0.2??) ??≠ ?? ??????= { 1.52 cos(??+ 1.2??) ??= ?? (???=, 1,2, ? ,10 ) 说明： 1. 在所用的算法中，凡是要给出精度水平的 ε，都取 ??= 10-12 。 打印以下内容： 采用带双步位移的 QR 分解法，说明算法设计方案和思路。 全部源程序。 (3) (??-1) 。 矩阵 A 经过拟上三角化后的矩阵 ?? (4) 对矩阵 (??-1) 实行 QR 方法迭代结束后所得的矩阵； ?? (5) 矩 阵 A 的全部实特征值??= ( ??, ) ( ??= 1,2, , ? ,10 ) ( ) ?? ???? ?? ?? ?? ????( ?) 。若 λ是实数，则令 ??= 0; ?? i ?? A 相应的实特征值的特征向量。 发现的现象与遇到的问题。 3. 采用 e 型数输出实型数，并且至少显示 12 位有效数字。 算法设计思路和方案 该问题的求解起始主要分为三个步骤， 第一步是对 A 拟上三角化得矩阵 (??-1) ；第二步 ?? (??-1) 进行带双步位移的 QR 分解迭代，并求出全部特征值；第三部是求出所有实 是对矩阵 ?? 特征值对应的特征向量。 关于第一步矩阵 A 的拟上三角化 对于此步书上已经有明确的思路和实现方法了， 只是在编程的过程中注意对下式在计算的时候应注意避免矩阵乘矩阵 (??+1) (??) ?? = ?????? ???? (??) 如果设 ?? ,(??= ??+ 1, ??+ 2, ? ,??不全为零，则有： ???? = (0, ? (??) ?? ?? ?? ,0, ?? ,? ,??????) ??+1,?? (??) ‖ ‖， (??) ‖ ‖ = -??????(??) ( 若?? = 0，择取 ????= ?? ?? ??+1,?? ??+1,?? ?? = ??- ???? ?? ?? = ???? | | ?? | ????| 2 ?? ????= ??- 2?????? 为了避免矩阵相乘采用如下步骤计算 (??+1) (??) ?? ′(??) ′(??) ?? = ???? ?? = ( ??- 2????) ?? = ?? - ?? ?? ?? ?? 其中  ?? ′(??) 2?????? ?? ?? ′(??) (??) (??) ?? ?? = (?? - 2?? ????) ?? ?? ； 计算中间向量 ?? = ?? ?? ?? 2. 从 ??= 1,2, ? ′(??) (??) ( ) ( ) ,??, ??= ??+ 1, ? ,??计算 ?? = ?? - 2???? ?????? ??, ???? ????  ′(??) ′(??) ?? = [?? ]； ???? 计算中间向量 从 ??= 1,2, ?  ′(??) = ???? ； ?? (??+1) ′(??) ,??,??= ??+ 1, ? ,??计算 ?? = ?? - ???? ????  ( ) ( ) (??+1) (??+1) ] ； 2???? ?????? ??, ?? = [?? ???? (??-??) 关于对矩阵 ?? 进行带双步位移的 QR 分解迭代 尽管书上给出了算法实现的 11 个步骤，然而在思路很混乱，在判断某步 QR 迭代后能 否对得到的矩阵 ?? 求部分特征值以达到降阶方面， 按照书上的判断方法， 几乎不可避免的 ??+1 得使用 goto 语句，这使得程序流程变得混乱。经过仔细分析后，可以发现能否求解 ?? 的 ??+1 部分特征值然后实现矩阵 ?? 降阶的关键在于判断一下拟上三角阵 ??+1 ??11 ??12 ? ? ??1?? ??21 ??22 ??23 ? ??2?? 0 ? ? ????-1??-2 ????-1??-1 ????-1?? [ 0 ? 0 ??????-1 ????