MATLAB当中一些简单的数值分析函数总结

求解函数的近似根：
二分法相关的MATLAB函数

function [x]=dichotomy(f,a,b,TOL) 
% f：原函数，a，b：求根区间，需要自定义给出，TOL：误差范围
N=1+fix(log2((b-a)/TOL));
%fix：取整数，并根据|X_N-X*|<=(b-a)/2^N<TOL,求N
n=1;
if (f(a)*f(b)>0)
    disp('根不在输入的区间里，请重新输入区间');
    return
else
    while n <= N
        x=(a+b)/2;
        if f(a)*f(x)>0
            a=x;
        else
            b=x;
        end
        n=n+1;
    end
    disp(x);
end
end

不动点迭代法相关的MATLAB函数
p.s. 请注意，不动点迭代法的原理当中，首先需要对原函数f(x)进行变化，再代入到下列函数中

function [x] = fixPoint(f,x0,TOL)
% f：原函数，x0：初始的迭代值，TOL：误差范围
x1 = f(x0);
while abs(x1 - x0)>TOL
    x0 = x1;
    x1 = f(x0);
end
x = x1;
disp(x)
end

Newton法相关的MATLAB函数

function [X]=Newton(f,p0,TOL,N)
format long;
n=1;%初始迭代次数
syms x;
while n<=N
    if abs(subs(diff(f(x)),x,p0))<TOL
        X=p0;
        break;
    else
        if subs(diff(f(x),2),x,p0)==0
            disp('Method failed');
            break;
        else
            p=p0-f(p0)/subs(diff(f(x)),x,p0);
            p=eval(p);%将exp的值转为小数值
            if(abs(p-p0)<TOL)
                X=p;
                break;
            else