weixin_41639107的博客

function Newton() x0=[0.1;0.5]; x1=x0-inv(myJacobi(x0))*myfun(x0); while norm(x1-x0)>1e-3 x0=x1; x1=x0-inv(myJacobi(x0))*myfun(x0); end x1 end function f=myfun(x) syms x1 x2 f1=(15*x1+10...

%% 牛顿法主程序 clear;clc format; x0=[0 0]; % 迭代初始值 eps = 0.00001; % 定位精度要求 for i = 1:10 f = double(subs(fun(x0),{'x1' 'x2'},{x0(1) x0(2)})); df = double(subs(dfun(x0),{'x1' 'x2'},{x0(1) x0(2)})...

function Newton() format long x0=[0;0;20;0;0;0]; x1=x0-inv(myJacobi(x0))*myfun(x0); while norm(x1-x0)>1e-3 x0=x1; x1=x0-inv(myJacobi(x0))*myfun(x0); end x1 end function f=myfun(x) syms x...

在如下链接中，本人看到了牛顿法求解非线性方程的一个程序，给人耳目一新的感觉。 https://blog.youkuaiyun.com/weixin_43429936/article/details/83720382 但不知何故，本人未能运行该博主写的程序，于是我对原程序做了一定修改，程序得以运行，代码如下： 牛顿法程序代码： function gen=newton(f,x,tol) %f为函数，x...

另一个牛顿法求解非线性方程的程序，优点是速度特别快，缺点是需要建立两个子程序，分别是待求函数子程序和待求函数的导数子程序。 主函数如下： function gen=fnewton(func,dfunc,x,tol) % func是函数表达式 % dfunc是这个函数的导数 % tol为截断误差 % x是初始值 x0=x; d=feval(func,x0)/feval(dfunc,x0);将...