精品文档
下面是基于 MATLAB 的迭代法程序，用迭代格式
p
n1
g
(
x
n
)
，求解方程
x g(x)
，其
中初始值为
p
0
。
**************************************************************************
function[p,k,err,P]=fixpt(f1021,p0,tol,max1)
% f1021 是给定的迭代函数。
% p0 是给定的初始值。
% tol 是给定的误差界。
% max1 是所允许的最大迭代次数。
% k 是所进行的迭代次数加 1。
% p 是不动点的近似值。
% err 是误差。
% P = {p1,p2,…,pn}
P(1) = p0;
for k = 2:max1
P(k) = feval('f1021', P(k-1));
k, err = abs(P(k) - P(k-1))
p = P(k);
if(err<tol),
break;
end
if k == max1
disp('maximum number of iterations exceeded');
end
end
P=P;
****************************************************************************
例 2.1 用上述程序求方程
sin x x
2
0
的一个近似解，给定初始值
x
0
0.5
，误差界为
10
5
。
解：先用 m 文件先定义一个名为 f1021.m 的函数文件。
function y = f1021(x)
收集于网络，如有侵权请联系管理员删除
