输出ry停机否则转5
5若k

N置k
k
1r0
ry

x

maxabsx
转3
否则输出失败信息停机
3Matlab程序代码
fu
ctio
ryrayleighAx0eps
r是特征值，y是特征向量
k1r00
yx0maxabsx0
xAy
迭代格式计算新的x
rdotyxdotyyReyleigh商
ifabsrr0eps
retur

e
d
whileabsrr0epsk

kk1
r0r
yxmaxabsx
xAy
rdotyxdotyy
e
d
e
d类似得计算按模最小特征值的Rayleigh商加速法，如下：
fu
ctio
ryi
vrayleighAx0eps

k1r00
yx0maxabsx0LUluA迭代格式不同
zLy
xUz
rdotyxdotyy
ifabsrr0eps
retur

e
d
whileabsrr0epsk

kk1
fr0ryxmaxabsxzLyxUzrdotyxdotyye
dr1re
d4计算Hilb矩阵特征值此处不再举例，而是直接应用于15阶Hilb矩阵，初始向量选为o
es151，
结果如下，并将结果与幂法和反幂法得到结果比较
这与幂法得到结果一致，表明算法和代码正确。
f同理，最小特征值如下：与反幂法得到结果一致，表明代码和算法正确。
fffr