解方程——不動點迭代

定義：若是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代碼實現：

%使用本函數須要一個事先創建的內聯函數g
function xc = fpi (g,x0,k)
x(1) = x0;
for i = 1:k
    x(i+1) = g(x(i));
end
x' %x用列顯示
xc = x(k+1);