功能:牛顿法求解非线性方程组。

调用示例

命令行输入:

[r, n]=mul_Newton([0 0 0], 1e-5)

得到结果:

r =

0.4996 -0.0900 -0.5259

n =

5

源码

function [r, n] = mul_Newton(x0, eps)

% 牛顿法求解非线性方程组

% 输入参数说明:

% 初始迭代向量:x0

% 迭代精度:eps

% 输出参数说明:

% 解向量:r

% 迭代次数:n

r = x0- myfun(x0)/dmyfun(x0);

n = 1;

tol = 1;

while tol > eps

x0 = r;

r = x0 - myfun(x0)/dmyfun(x0);

tol = norm(r - x0);

n = n + 1;

if(n > 100000)

disp('迭代步数太多,可能不收敛!');

return;

end

end

function f = myfun(x)

x1 = x(1);

x2 = x(2);

x3 = x(3);

f(1) = 3*x1 - cos(x2*x3) -1/2;

f(2) = x1^2 - 81*(x2 + 0.1) + sin(x3) + 1.06;

f(3) = exp(-x1*x2) + 20*x3 + 1/3*(10*pi - 3);

function df =dmyfun(x)

x1 = x(1);

x2 = x(2);

x3 = x(3);

df=[3, x3*sin(x2*x3), x2*sin(x2*x3); 2*x1, -81, cos(x3); -x2*exp(-x1*x2), -x1*exp(-x1*x2), 20];

---------------------

作者:有个家伙喜欢代码

来源:CSDN

原文:blog.csdn.net/u012366767 …

版权声明:本文为博主原创文章,转载请附上博文链接!

matlab 牛顿法 非线性方程,Matlab学习手记——非线性方程组求解:牛顿法相关推荐

  1. Matlab学习手记——非线性方程组求解:牛顿下山法

    功能:牛顿下山法求解非线性方程组. 牛顿下山法 function [x, n] = NonLinearEquations_NewtonDown(x0, err) %{ 函数功能:牛顿下山法求解非线性方 ...

  2. fsolve:对非线性方程组求解

    fsolve 对非线性方程组求解 一:语法 x = fsolve (fun,x0) x = fsolve (fun,x0,options) x = fsolve (problem) [ x,fval ...

  3. matlab非线性方程组求解得到矩阵,非线性方程组求解——附Matlab原程序

    在科学与工程计算中,经常遇到求解非线性方程组的问题:非线性方程组在收敛速度及收敛性比线性方程组要差,特别对于非凸的非线性方程组,其求解更是困难.下面简要介绍非线性方程组的三种解法--牛顿法.拟牛顿法. ...

  4. matlab 非线性方程组画图,使用matlab进行非线性方程组求解

    今天,恰好遇到使用matlab进行了非线性方程组的求解. 情况是,一些无法手动表示,推导的方程组,比如像超越方程这种,可以用matlab进行直接求值,这在一些需要仿真数据的场景比较有用. 环境:mat ...

  5. 求解非线性方程组的牛顿法c语言,牛顿下山法求解非线性方程(组)(C实现)...

    1.算法描述 (1)符号说明与基本假设 对于非线性方程组:                                                        (1) 引入向量: 可将(1) ...

  6. matlab tolfun,非线性方程组求解问题(关于TolFun设置问题)

    function main clear all clc global n G dl M M1 b1 b2 C4 C6 C7 C8 C9 a T=273.2+112.13;                ...

  7. c语言割线法解非线性方程,用割线法解非线性方程组.doc

    用割线法解非线性方程组 自动化学院1011203050 陈晓祺 拟牛顿法解下列方程组 先将拟牛顿法的程序代码如下 Function[r,m]=mulVlineF,x0,A, eps) Format l ...

  8. 非线性方程组求解方法,神经网络的非线性函数

    1.rbf神经网络原理 rbf神经网络原理是用RBF作为隐单元的"基"构成隐含层空间,这样就可以将输入矢量直接映射到隐空间,而不需要通过权连接. 当RBF的中心点确定以后,这种映射 ...

  9. python fsolve说明_python用fsolve、leastsq对非线性方程组求解

    背景: 实现用python的optimize库的fsolve对非线性方程组进行求解.可以看到这一个问题实际上还是一个优化问题,也可以用之前拟合函数的leastsq求解.下面用这两个方法进行对比: 代码 ...

最新文章

  1. 你的团队能通过电梯测试吗?
  2. svn还原文件中去掉已经删除的文件
  3. Collections练习之按照字符串长度进行排序
  4. oracle数据库开多线程,学习笔记:Oracle表数据导入 DBA常用单线程插入 多线程插入 sql loader三种表数据导入案例...
  5. Bootstrap 可切换的标签导航
  6. js取thymeleaf值_26、模板引擎thymeleaf
  7. 随机读写工具,手写,百度云源码直接下载
  8. GitHub 使用教程图文详解(转)
  9. zend studio html插件安装,Zend Studio使用教程:将Zend Studio作为插件安装
  10. android dp转px的公式_dp与px的转换对照表
  11. Ubuntu设置静态IP地址的几种方法
  12. Java 服务器版超市管理系统,基于JAVA的超市进销存管理系统.doc
  13. zblog php建站教程_开源建站程序ZBlogPHP博客安装图文教程
  14. 3-3 uniapp、HTML5+、Native.js 功能代码汇总
  15. [攻防世界]getit
  16. MFC几种给对话框添加背景图的方法
  17. 【Python脚本抢红包】用Python写了一个自动抢红包工具,今年过年准备大干一场
  18. brew对redis的使用
  19. QQ浏览器X5内核问题汇总 转
  20. 【IPv6+燎原系列—第2期】IPv6+逐梦赛场,勇夺属于通信人的“金牌”

热门文章

  1. 创新性中小企业认定标准
  2. RS485与Modbus通信协议教程!
  3. 牛客 15898 黑猫的小老弟【法里数列】
  4. 2013年第十九届全国青少年信息学奥林匹克联赛初赛
  5. make clean和make clobber区别
  6. 企业级架构之LNMP
  7. python高效识别验证码 极简 dddd
  8. 中国聚酰亚胺泡沫市场预测与战略咨询报告(2022版)
  9. 最短路径 floyd最小环 洛谷2738 篱笆回路 网上题解
  10. python hook(钩子)函数