首页 > 代码库 > 解方程——不动点迭代

解方程——不动点迭代

定义:如果g(r)=r,那么实数r是函数g的一个不动点。

如果我们有方程f(x)=0,表示为不动点问题时,有:g(x)=x   (注:f(x)=g(x)-x)

不动点r是方程f(x)=0的一个根,几何表示为y=g(x)与y=x的交点就是不动点r

算法分析:x0 = 初始设定值

              x1 = g(x0)

              x2 = g(x1)

              x3 = g(x2)

               ...    ...

             x(k+1) = g(x(k))  

直到收敛至g(r) = g(lim x(i)) = lim x(i+1) = r

matlab代码实现:

%使用本函数需要一个事先建立的内联函数gfunction xc = fpi (g,x0,k)x(1) = x0;for i = 1:k    x(i+1) = g(x(i));endx‘ %x用列显示xc = x(k+1);

 

解方程——不动点迭代