原创力文档- 1
- 0
- 约2.18千字
- 约 4页
- 2017-06-08 发布于重庆
- 举报
五线性方程组迭代解法
一、编写用谱半径判别法产生的迭代序列的敛散性的MATLAB函数，并调用之。
function B ddpbj H
B eig H ;mB norm B,inf ;
if mB 1 disp 谱半径不小于1，迭代序列发散，谱半径mB和H的所有特征值B如下：
else disp 谱半径小于1，迭代序列收敛，谱半径mB和H的所有特征值B如下：
end
mB
H [10 -1 -2;-1 10 -2;-1 -1 5]
H 10 -1 -2 -1 10 -2 -1 -1 5 B ddpbj H
谱半径不小于1，迭代序列发散，谱半径mB和H的所有特征值B如下：
mB 11.0000
B 4.1716 11.0000
9.8284
二、用jacobi迭代法解线性方程组的的函数。
function X jacdd A,b,X0,P,wucha,max1
[n m] size A ;
for j 1:m a j sum abs A :,j -2* abs A j,j ;
end
for i 1:n if a i 0 disp 系数矩阵A非严格对角占优，jacobi迭代不一定收敛 return end
end
if a i 0 disp 系数矩阵A严格对角占优，方程组有唯一解，且jacobi迭代收敛
end
for k 1:max1 k fo
文档评论（0）