MATLAB学习笔记之牛顿迭代法

牛顿迭代法的公式如下
在这里插入图片描述
给定函数f(x)和近似根值x0,由上式逐步逼近根值;直到误差小于给定值;
代码如下:

function x1=newton(fun,x0,delta)
syms x;%定义符号变量x,用于求导和代入
df=diff(fun(x));%对f(x)进行求导
x1=x0-fun(x0)/subs(df,x0);
while abs(x1-x0)>delta
x0=x1;
x1=x0-fun(x0)/subs(df,x0);
end

这里提一下subs函数,subs函数用于替换表达式中的符号变量,在第四行代码中,

x1=x0-fun(x0)/subs(df,x0);

即用x0的值来替换df表达式中的变量,否则编译器会报错,提示“ 未定义的变量值 ”。

你可能感兴趣的:(MATLAB学习笔记之牛顿迭代法)