(数值分析 反幂法 反幂法用来求矩阵A的按模最小特征值及其相应的 特征向量 只要求出A1的按模最大的特征值,也就求出了A的按 模最小的特征值,及其相应的特征向量 §2 Jacobi方法 Jacobi方法用来求实对称矩阵的所有特征值和相应 的特征向量 经典的 Jacobi方法 经典的 Jacobi方法的特点是每次变换将绝对值最大 的非对角元素化为零 36 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 3636 三、 反幂法 反幂法用来求矩阵A的按模最小特征值及其相应的 特征向量. 只要求出A-1的按模最大的特征值, 也就求出了A的按 模最小的特征值, 及其相应的特征向量. §2 Jacobi方法 Jacobi方法用来求实对称矩阵的所有特征值和相应 的特征向量. 经典的Jacobi方法 经典的Jacobi方法的特点是每次变换将绝对值最大 的非对角元素化为零。 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 计算步骤为 1)首先在4中选择绝对值最大的非对角线元素,设 maxa.≠0 1≠ 选择平面旋转矩阵S1=S(i,i)使T1=S1AS1的非对角线 元素 )=0 (2)再在T1中选择绝对值最大的非对角线元素,设 (1) max ≠0 12J2 ≠ 选择平面旋转矩阵S2=S(2i)使T2=S2T1S2的非对角线 元素 (2)二f( 0 12/2 此时的t1(2)=51(2又可能变换成了非零元素 (3)重复以上过程,直到满足预定精度为止 37 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 3737 计算步骤为 (1)首先在A中选择绝对值最大的非对角线元素, 设 1 1 max 0 i j ij i j t a ¹ = ¹ 选择平面旋转矩阵S1 =S(i1 , j1 )使T1 =S1 TAS1的非对角线 元素 1 1 1 1 (1) (1) 0 i j j i t t = = (2)再在T1中选择绝对值最大的非对角线元素, 设 2 2 (1) (1) max 0 i j ij i j t t ¹ = ¹ 2 2 2 2 (2) (2) 0 i j j i t t = = 此时的t i1j1 (2)=t j1i1 (2)又可能变换成了非零元素. 选择平面旋转矩阵S2 =S(i2 , j2 )使T2 =S2 TT1S2的非对角线 元素 (3)重复以上过程, 直到满足预定精度为止. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 Jacobi过关法 经典的 Jacobi方法每次选取的是矩阵中绝对值最大 的非对角元素作为消去对象,需要在所有的非对角线元 素中进行比较选择,计算工作量相当大. Jacobi过关法也 称为“阈” Jacobi方法,是一种改进方法 计算 ∑∑ =lj=1,j≠ (1)设置阈值,也称为“关” ∑∑ i=1j=1,j 扫描矩阵所有非对角线元素,对绝对值小于阈值v1的 元素,就让其过“关”,暂不作处理.对绝对值 31函 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 Jacobi过关法 经典的Jacobi方法每次选取的是矩阵中绝对值最大 的非对角元素作为消去对象, 需要在所有的非对角线元 素中进行比较选择, 计算工作量相当大. Jacobi过关法也 称为“阈”Jacobi方法, 是一种改进方法. 2 0 1 1, n n ij i j j i v a = = ¹ 计算 =å å (1)设置阈值, 也称为“关” 2 1 1, 1 n n ij i j j i a v n = = ¹ = å å 扫描矩阵所有非对角线元素, 对绝对值小于阈值v1的 元素, 就让其过“关” , 暂不作处理. 对绝对值 3838 PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 大于等于阈值v1的元素,就构造平面旋转矩阵,并利用 旋转变换将其变为零.多次扫描非对角元素,直到所有 的非对角线元素的绝对值都小于阈值v为止 (2)缩小阈值 ∑∑ j=l,j≠i 多次扫描矩阵从(1)所得矩阵的非对角线元素,并做 相应的平面旋转变换.直到所有的非对角线元素的绝 对值都小于阈值v2为止 (k)设置阈值"k n少x≤E8为给定的误差精度 39 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 3939 大于等于阈值v1的元素, 就构造平面旋转矩阵, 并利用 旋转变换将其变为零. 多次扫描非对角元素, 直到所有 的非对角线元素的绝对值都小于阈值v1为止. (2)缩小阈值 2 1 1 1, 2 2 n n ij i j j i a v v n n = = ¹ = = å å 多次扫描矩阵从(1)所得矩阵的非对角线元素, 并做 相应的平面旋转变换. 直到所有的非对角线元素的绝 对值都小于阈值v2为止. . . . . . . . . . . . . . . . . . . . . . . . . . . . . (k)设置阈值 k 1 0 k k v v v n n e - = = £ e为给定的误差精度. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com
(数值分析 §3QR方法 QR方法可以用来求一般矩阵所有特征值 构造一个矩阵序列{4} A=FR k=1.2 RF 4是相似矩阵序列,因此它们具有相同特征值 (1)当取F为下三角矩阵L时,R为上三角矩阵时,就得 到矩阵的LR分解. 包)当取F为正交矩阵Q时,R为上三角矩阵时,就得 矩阵的QR分解 PDF檔案使用"pdfFactoryPro"試用版本建立www.pdffactory.com
数值分析 4040 §3 QR方法 QR方法可以用来求一般矩阵所有特征值. 1 1 1, 2, k k k k k k A A A F R k A + R F ì = ï í = = ï î = L 构造一个矩阵序列{Ak } {Ak }是相似矩阵序列, 因此它们具有相同特征值. (1)当取Fk为下三角矩阵L时, Rk为上三角矩阵时,就得 到矩阵的LR分解. (2)当取Fk为正交矩阵Q时, Rk为上三角矩阵时, 就得 到矩阵的QR分解. PDF 檔案使用 "pdfFactory Pro" 試用版本建立 www.pdffactory.com