用vrho()求谱半径。例如输入矩阵A=[1 0; 0 3],用vrho(A)求得谱半径为3
A=[10 3 1;2 -10 3;1 3 10];
b=[14 -5 14]';
D=diag(diag(A));
L=-tril(A,-1);U=-triu(A,1);
B=D\(L+U);f=D\b;
x=[0;0;0];
for k=1:9
x=B*x+f;
x'
end 其中B矩阵的矩阵半径:R=max(abs(eig(B)))=0.3873<1 收敛
高斯-赛德尔迭代
matlab2013中没法用 vrho(),怎么回事呢
R=max(abs(eig(B))) 其中的上、下三角矩阵 对角矩阵 自己应该知道吧?
还是看看这个百科 http://baike.baidu.com/view/2012163.htm 有例子的