2、新建dfun.m,求出一阶微分方程 1function df=dfun(x);2f=fun(x);3df=[diff(f,'x1');diff(f,'x2')]; %雅克比矩阵 3、建立newton.m,执行牛顿迭代过程 1clear;clc2format;3x0=[0 0]; %迭代初始值4eps = 0.00001; %定位精度要求5fori = 1:106f = double(subs(fun(x0),{'x1''x2'},...
⽜顿迭代法解⾮线性⽅程组（MATLAB版）⽜顿迭代法，⼜名切线法，这⾥不详细介绍，简单说明每⼀次⽜顿迭代的运算：⾸先将各个⽅程式在⼀个根的估计值处线性化（泰勒展开式忽略⾼阶余项），然后求解线性化后的⽅程组，最后再更新根的估计值。下⾯以求解最简单的⾮线性⼆元⽅程组为例...
2、新建dfun.m,求出一阶微分方程 1function df=dfun(x);2f=fun(x);3df=[diff(f,'x1');diff(f,'x2')]; %雅克比矩阵 3、建立newton.m,执行牛顿迭代过程 1clear;clc2format;3x0=[0 0]; %迭代初始值4eps = 0.00001; %定位精度要求5fori = 1:106f = double(subs(fun(x0),{'x1''x2'},...