%程序:nton.m
%功能:牛顿迭代法(MATLAB R2009a)
%说明:fx是char型的MATLAB的有效表达式,以x为自变量,并且可求导数
%     x0迭代初值
%     tol迭代误差
%     N迭代的上限
%     z返回值
%     n返回迭代的次数
%实例:用牛顿迭代法求方程x*exp(x)=1的正根
%     >>[z,n]=nton('x*exp(x)-1',0,10^(-5))
function [z,n]=nton(fx,x0,tol,N)
if nargin<2
    error('too less parameters...');
elseif nargin<3
    N=500;tol=10^(-3);
elseif nargin<4
    N=500; 
else 
    error('too many parameters...')
end
syms x;
dfx=char(diff(fx));  %sym型转换为char型
fx=inline(fx);        %内联函数,方便用feval求解
dfx=inline(dfx);  
x1=x0-feval(fx,x0)./feval(dfx,x0);
k=1;
while (abs(x1-x0)>tol)&(k<N)
    x0=x1;
    x1=x0-feval(fx,x0)/feval(dfx,x0);
    k=k+1;
end
if nargout==2
    z=x1;
    n=k;
elseif nargout==1
    z=x1;
elseif nargout==0
    disp('the answer is:'); disp(x1);
end

MATLAB 牛顿迭代算法相关推荐

  1. 秦九邵算法,牛顿迭代算法,盛金公式

    基础数学算法 1.牛顿迭代法 (1)为什么要牛顿迭代 (2)实现牛顿跌代 2.秦九邵算法 (1).历史由来 (2).算法介绍 3.盛金公式 1.牛顿迭代法 (1)为什么要牛顿迭代 啥是牛顿迭代法?大大 ...

  2. 用matlab画牛顿迭代图形,matlab牛顿迭代动画演示

    % 牛顿环演示的MATLAB程序 close all;clear;clc; figure('Position',[90 164 873 483]); L=632.8;R=5;H=5; a1=axes( ...

  3. python迭代算法含义_迭代算法原理及其Python编程实现

    黄旭 摘 要:迭代算法是数学算法在计算机中应用的一个热点,也是计算机解决问题的一般思路,本文结合数学中二分法求根的原理,阐述了数学迭代算法的一般原理,并采用了Python加以实现,为进一步对数学算法理 ...

  4. 修正牛顿法matlab,牛顿算法及其改进【阻尼牛顿法、修正牛顿法】

    牛顿算法 对于优化函数\(f(x)\),\(x=(x_1;x_2;...;x_n)\),二阶连续可导 在\(x_k\)处泰勒展开,取前三项,即对于优化函数二阶拟合 \[f(x)=f(x_k)+g_k( ...

  5. 【图像重建】基于matlab布雷格曼迭代算法集合ART算法CT图像重建【含Matlab源码 1905期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[图像重建]基于matlab布雷格曼迭代算法集合ART算法CT图像重建[含Matlab源码 1905期] 获取代码方式2: 通过订阅紫极神光 ...

  6. 有限元学习笔记-虚功平衡方程的线性化及牛顿-拉普森迭代算法

    有限元学习笔记-虚功平衡方程的线性化及牛顿-拉普森迭代算法 张量微分是解连续介质物理问题的利器.

  7. 求平方根的算法 牛顿迭代法和二分法

    牛顿法: public double sqr(double n){double x=n,y=0.0;while (Math.abs(x-y)>0.0001){y=x;x=(x+n/x)/2;}r ...

  8. BCH码和m序列参数估计(梅西迭代算法求多项式的MATLAB实现)

    关注公众号[逆向通信猿]试读更多内容!!! 伯利坎普迭代算法解决了BCH码译码中求错位多项式 σ ( x ) \sigma (x) σ(x)的问题,即用迭代算法求解方程组,极大地加快了译码速度.求最短 ...

  9. python牛顿迭代法_Jacobi迭代算法的Python实现详解

    import numpy as np import time 1.1 Jacobi迭代算法 def Jacobi_tensor_V2(A,b,Delta,m,n,M): start=time.perf ...

最新文章

  1. SpringBoot2.0 基础案例(05):多个拦截器配置和使用场景
  2. [BZOJ 1014][JSOI2008]火星人prefix(Splay+二分+hash)
  3. 检测数据类型的几种方式
  4. boost::posix_time模块实现计算时间段的测试程序
  5. python测试开发django-44.xadmin上传图片和文件
  6. Invalid nested tag br found, expected closing tag li 错误解决
  7. MySQL亿级数据数据库优化方案测试-银行交易流水记录的查询
  8. CSS 之 样式优先级机制
  9. NeurIPS 2021揭榜,接收率创九年新高,论文列表已公布,你的文章中了吗?
  10. iOS 刘海屏 隐藏横条 home键
  11. 轻松决绝疑难杂症:nodejs使用http报错: Request path contains unescaped characters
  12. 对西方国家的智能制造研究综述:过去现在和未来
  13. mysql触发器trigger 实例详解
  14. mac 误删引导分区
  15. 使用MySQL存储过程给表批量造数据
  16. 【第七篇】Flowable核心内容之任务分配
  17. 57个你没有听过的Google产品
  18. win10自带图片出现文件系统错误 (-2147219196)的修复方法
  19. jshop测试分析总览
  20. 结合NBA主题,Niantic再寻LBS AR游戏破局之道

热门文章

  1. [转贴]ATOM和RSS的区别
  2. 用asp.net 2.0实现网站二级域名(转)
  3. 公网mysql_(看图)为什么无法通过公网登录MySQL?
  4. 如何计算k段流水线执行n条指令的执行时间
  5. XML Schema全接触 (这里主要介绍W3C的Schema标准语法)
  6. python加载图片并显示_python OpenCV 读取并显示图像
  7. 学php5还是php7,学习猿地-php5.*与php7.* 区别
  8. python水平_python水平
  9. yml文件配置mysql表大小写_springboot常用配置(yml文件)
  10. XML--XML从入门到精通 Part 1 认识XML