matlab求松弛迭代发,MATLAB实现迭代法最佳松弛因子的选取

1 迭代法最佳松弛因子的选取

一、问题提出:

针对矩阵430341014A ⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦

，b=[24;30;-24],用SOR 迭代求解。并选出最佳松弛因子。

理论分析 1.24ω=

=≈。做出()L ωρ关于ω函数的图像。 二、理论基础

选取分裂矩阵M 为带参数的下三角矩阵)(1wL D w

M -=

, 其中w>0为可选择的松弛因子.

于是,由⎪⎩⎪⎨⎧+=+f Bx x x k k )()1()0()(初始向量 (k=0,1,…,)可构造一个迭代法,其迭代矩阵为A wL D w I L w 1)(---≡

=).)1(()(1wU D w wL D +---

从而得到解Ax=b 的主次逐次超松弛迭代法.

解Ax=b 的SOR 方法为

⎪⎩⎪⎨⎧+=+f Bx x

x k k )()1()0()(初始向量 (k=0,1,…,) (1) 其中

w L =).)1(()(1wU D w wL D +---(2)

b wL D w f 1)(--=

下面给出解Ax=b 的SOR 迭代法的分量计算公式.记

,),...,,...,()()()(1)(T k n k i

k k x x x x = 由(1)式可得

,))1(()()()1(wb x wU D w x wL D k k ++-==-+

).()()()1()()1(k k k k k Dx Ux Lx b w Dx Dx -+++=++ (3) 由此,得到解Ax=b 的SOR 方法的计算公式