MATLAB 牛顿迭代算法
%程序: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)为什么要牛顿迭代 (2)实现牛顿跌代 2.秦九邵算法 (1).历史由来 (2).算法介绍 3.盛金公式 1.牛顿迭代法 (1)为什么要牛顿迭代 啥是牛顿迭代法?大大 ...
- 用matlab画牛顿迭代图形,matlab牛顿迭代动画演示
% 牛顿环演示的MATLAB程序 close all;clear;clc; figure('Position',[90 164 873 483]); L=632.8;R=5;H=5; a1=axes( ...
- python迭代算法含义_迭代算法原理及其Python编程实现
黄旭 摘 要:迭代算法是数学算法在计算机中应用的一个热点,也是计算机解决问题的一般思路,本文结合数学中二分法求根的原理,阐述了数学迭代算法的一般原理,并采用了Python加以实现,为进一步对数学算法理 ...
- 修正牛顿法matlab,牛顿算法及其改进【阻尼牛顿法、修正牛顿法】
牛顿算法 对于优化函数\(f(x)\),\(x=(x_1;x_2;...;x_n)\),二阶连续可导 在\(x_k\)处泰勒展开,取前三项,即对于优化函数二阶拟合 \[f(x)=f(x_k)+g_k( ...
- 【图像重建】基于matlab布雷格曼迭代算法集合ART算法CT图像重建【含Matlab源码 1905期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[图像重建]基于matlab布雷格曼迭代算法集合ART算法CT图像重建[含Matlab源码 1905期] 获取代码方式2: 通过订阅紫极神光 ...
- 有限元学习笔记-虚功平衡方程的线性化及牛顿-拉普森迭代算法
有限元学习笔记-虚功平衡方程的线性化及牛顿-拉普森迭代算法 张量微分是解连续介质物理问题的利器.
- 求平方根的算法 牛顿迭代法和二分法
牛顿法: 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 ...
- BCH码和m序列参数估计(梅西迭代算法求多项式的MATLAB实现)
关注公众号[逆向通信猿]试读更多内容!!! 伯利坎普迭代算法解决了BCH码译码中求错位多项式 σ ( x ) \sigma (x) σ(x)的问题,即用迭代算法求解方程组,极大地加快了译码速度.求最短 ...
- 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 ...
最新文章
- SpringBoot2.0 基础案例(05):多个拦截器配置和使用场景
- [BZOJ 1014][JSOI2008]火星人prefix(Splay+二分+hash)
- 检测数据类型的几种方式
- boost::posix_time模块实现计算时间段的测试程序
- python测试开发django-44.xadmin上传图片和文件
- Invalid nested tag br found, expected closing tag li 错误解决
- MySQL亿级数据数据库优化方案测试-银行交易流水记录的查询
- CSS 之 样式优先级机制
- NeurIPS 2021揭榜,接收率创九年新高,论文列表已公布,你的文章中了吗?
- iOS 刘海屏 隐藏横条 home键
- 轻松决绝疑难杂症:nodejs使用http报错: Request path contains unescaped characters
- 对西方国家的智能制造研究综述:过去现在和未来
- mysql触发器trigger 实例详解
- mac 误删引导分区
- 使用MySQL存储过程给表批量造数据
- 【第七篇】Flowable核心内容之任务分配
- 57个你没有听过的Google产品
- win10自带图片出现文件系统错误 (-2147219196)的修复方法
- jshop测试分析总览
- 结合NBA主题,Niantic再寻LBS AR游戏破局之道
热门文章
- [转贴]ATOM和RSS的区别
- 用asp.net 2.0实现网站二级域名(转)
- 公网mysql_(看图)为什么无法通过公网登录MySQL?
- 如何计算k段流水线执行n条指令的执行时间
- XML Schema全接触 (这里主要介绍W3C的Schema标准语法)
- python加载图片并显示_python OpenCV 读取并显示图像
- 学php5还是php7,学习猿地-php5.*与php7.* 区别
- python水平_python水平
- yml文件配置mysql表大小写_springboot常用配置(yml文件)
- XML--XML从入门到精通 Part 1 认识XML