牛顿法迭代法 | matlab实现

function Newt(file_name, x0, xmin, xmax, n_points)
% file_name 是要进行迭代的函数
% x0 是选取的初始值
% xmin，xmax分别为图形横坐标的最小和最大值
% n_points 是自变量x的采样数clf, hold offwid_x = xmax - xmin;
dx = (xmax - xmin) / n_points;
xp = xmin: dx: xmax;
yp = feval(file_name, xp);
plot(xp, yp, 'r');
xlabel('x');
ylabel('f(x)');
title('Newton Iteration');
hold onymin = min(yp);
ymax = max(yp);
wid_y = ymax - ymin;
yp = 0.0 * xp;
plot (xp, yp);
x = x0;
xb = x + 999;
n=0;
del_x = 0.001;
while abs(x - xb) > 0.000001;
if n>300break;
endy = feval(file_name, x);plot([x, x], [y, 0], 'black');plot(x, 0, 'o')fprintf('n = %3.0f,  x = %12.5e, y = %12.5e\n', n, x, y);
if n<4text(x, -wid_y/10, [num2str(n)]);
endy_driv = (feval(file_name, x+del_x)-y)/del_x;xb = x;x = xb-y/y_driv;n = n+1plot([xb, x], [y, 0], 'g')
end
plot([x, x], [0.05 * wid_y, 0.2 * wid_y], 'r');
text(x, 0.25 * wid_y, 'Final solution');
plot([x, (x - wid_x * 0.004)], [0.01 * wid_y, 0.09 * wid_y], 'r');
plot([x, (x + wid_x * 0.004)], [0.01 * wid_y, 0.09 * wid_y], 'r');
disp("牛顿迭代法结束")
end

牛顿法迭代法 | matlab实现相关推荐

牛顿迭代法例题 matlab,牛顿迭代法-matlab程序(解线性方程组)
牛顿迭代法 matlab程序(解线性方程组) 作者:佚名来源:转载发布时间:2009-3-7 16:55:53 减小字体增大字体 1．功能本程序采用牛顿法,求实系数高次代数方程 f(x)=a0xn ...
牛顿迭代法解线性方程matlab程序,牛顿迭代法matlab程序(解线性方程组)
<牛顿迭代法matlab程序(解线性方程组)>由会员分享,可在线阅读,更多相关<牛顿迭代法matlab程序(解线性方程组)(4页珍藏版)>请在金锄头文库上搜索. 1.牛顿迭代法 ...
牛顿迭代法-matlab实现
牛顿迭代法-matlab实现牛顿迭代法简介: 牛顿迭代法又称为切线法,简单来说就是不断求切线与x轴的交点,来逐渐接近解的迭代过程.方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的 ...
matlab编jacobi方法,jacobi迭代法matlab
h Matlab 线性方程组的迭代解法(Jacobi 迭代法 Gauss-Seidel 迭代法)实验报告 2008年11月09日星期日 12:49 1.熟悉 Jacobi 迭代法,并编写 Matla ...
matlab 高斯迭代法求解,高斯迭代法matlab算例
Matlab 线性方程组的迭代解法 Gauss-Seidel 迭代法 Matlab 线性方程组的迭代解法 Gauss-Seidel 迭代法实验报告 1.熟悉 Gauss-Seidel 迭代法,并编写 ...
matlab牛顿法原理,牛顿法的matlab实现
牛顿法的matlab实现牛顿法的 MATLAB 实现摘要:性能学习是神经网络中的一类很重要的学习规则,其旨在找到一个最优点来提高网络的性能.牛顿法是一种基于二阶泰勒级数的算法,逐步迭代来实现进一步的 ...
非线性方程组牛顿迭代法matlab,matlab实现牛顿迭代法求解非线性方程组
<matlab实现牛顿迭代法求解非线性方程组>由会员分享,可在线阅读,更多相关<matlab实现牛顿迭代法求解非线性方程组(5页珍藏版)>请在人人文库网上搜索. 1.matla ...

牛顿法迭代法 | matlab实现

牛顿法迭代法 | matlab实现相关推荐

最新文章

热门文章