习题
一．编写 JACOBI 法求解实对称矩阵特征值和特征向量的程序：
1.从文本文件 mat.dat 中读取数据并存入变量 A；
2.如果 A 的非对角线绝对值最大元素的绝对值大于 eps 则需执行以下步骤：
1）搜索 A 的非对角线绝对值最大元素，得到此元素的行号 p 和列号 q。此步骤要求写成下
列函数：
function [p,q,s] = maxelement(A)
其中， A：矩阵 A；
p：行号；
q：列号；
s：A 中绝对值最大的元素的绝对值。
2）形成 Givens 矩阵。此步骤要求写成下列函数：
function [R] = givens(p,q,A)
其中， A：矩阵 A；
p：行号；
q：列号；
R：Givens 矩阵。
3）执行相似变换并计算特征向量，直到 s 小于 eps；
4）将特征值和特征向量按从小到大排列。步骤 1）至 4）写成下列函数：
function [V,D] = jacobi(A)
其中， A：矩阵；
V：由特征向量组成的矩阵，按列存放；
评论0