MATLAB简单迭代法求方程解的例子例题
将方程f(x)=0,经过简单变形处理,转化为等价的方程x= ϕ (x),设置迭代格式xk+1 = ϕ (x k ),若迭代收敛,即方程收敛到不动点x*。
例题 设置迭代格式求解以下方程:
x-e -x =0。
解: 由x=e -x ,可以得到迭代格式:x k+1 =e -x k 。
程序如下:
function y=exam(n)
%简单迭代法求根,n表示迭代次数
y=ones(1,n);
for i=1:n
y(i+1)=exp(-y(i));
end
在命令窗口下执行:
y=exam(10)
结果如下:
ans =
1.0000 0.3679 0.6922 0.5005 0.6062 0.5454
0.5796 0.5601 0.5711 0.5649 0.5684