?? dimethod.m
字號:
function [l,v]=dimethod(A,x0,u,eps)
if nargin==3
eps = 1.0e-6;
end
N = size(A);
n = N(1,1); %n為A的階
v = x0;
M = 5000; %迭代步數限制
m = 0;
l = 0;
for(k=1:M)
y = (A-u*eye(n,n))\v;
m = max(y); %m為按模最大的分量
v = y/m;
if(abs(m - l)<eps)
l = 1/m + u; %到所需精度,退出
return;
else
if(k==M)
disp('迭代步數太多,收斂速度太慢!');
l = 1/m + u;
else
l = m;
end
end
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -