Question

您好！想请教一下Matlab 拟牛顿求根法

function f=point1(x,xp,a1,a2,a3,a4,angle0,lc0,lb,mb,E,I,P)
一条方程，只有x是未知数，如何使用牛顿迭代呢

Answer 1

拟牛顿法是求解线性（非线性）方程（组）的一种数值分析方法。拟牛顿法的求解方法与割线法相类似，其目的是减少由于计算导数而带来的大量的计算量，构造的秩为1的迭代法。其迭代格式为

实现主要代码

n = 1;

while (norm(x - x0) > tol) && (n < 1000)

x0 = x;

x = x0 - A \ funm(x0);

p = x - x0;

q = funm(x) - funm(x0);

A = A + (q - A*p)*p'^rm(p)^2;

n = n + 1;

end