Matlab實現——to find the fixed point

歡迎前往個人博客 駑馬點滴 和視頻空間 嗶哩嗶哩-《挨踢日誌》

fixpt.m

 

function [k,p,err,P]=fixpt(g,p0,tol,max1)

% Input - g is the iteration function input as a string ??g??

%           - p0 is the initial guess for the fixed-point

%           - tol is the tolerance

%           - max1 is the maximum number of iterations 

% Output - k is the number of iterations that were carried out

%         - p is the approximation to the fixed-point

%         - err is the error in the approximation

%         - P contains the sequence {pn}

P(1)=p0;

for k=2:max1

    P(k)=feval(g,P(k-1));

    err=abs(P(k)-P(k-1));

    relerr=err/(abs(P(k))+eps);

    p=P(k);

    if (err<tol) | (relerr<tol),break;end

end

if k==max1

    disp('maximum number of iterations exceeded')

end

P=P';

程序示例：

fun.m

 

function g=fun(x)

g=x^5-3*x^3-2*x^2+2;

Untitled3.m

p0=-1.6;

tol=10e-12;

max1=7;

g='fun';

[k,p,err,P]=fixpt(g,p0,tol,max1)

運行的結果說明在點(-1.6,-1.8)附近的不動點是排斥不動點

 

歡迎前往個人博客 駑馬點滴 和視頻空間 嗶哩嗶哩-《挨踢日誌》