matlab 牛顿法 非线性方程,Matlab学习手记——非线性方程组求解:牛顿法
功能:牛顿法求解非线性方程组。
调用示例
命令行输入:
[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学习手记——非线性方程组求解:牛顿法相关推荐
- Matlab学习手记——非线性方程组求解:牛顿下山法
功能:牛顿下山法求解非线性方程组. 牛顿下山法 function [x, n] = NonLinearEquations_NewtonDown(x0, err) %{ 函数功能:牛顿下山法求解非线性方 ...
- fsolve:对非线性方程组求解
fsolve 对非线性方程组求解 一:语法 x = fsolve (fun,x0) x = fsolve (fun,x0,options) x = fsolve (problem) [ x,fval ...
- matlab非线性方程组求解得到矩阵,非线性方程组求解——附Matlab原程序
在科学与工程计算中,经常遇到求解非线性方程组的问题:非线性方程组在收敛速度及收敛性比线性方程组要差,特别对于非凸的非线性方程组,其求解更是困难.下面简要介绍非线性方程组的三种解法--牛顿法.拟牛顿法. ...
- matlab 非线性方程组画图,使用matlab进行非线性方程组求解
今天,恰好遇到使用matlab进行了非线性方程组的求解. 情况是,一些无法手动表示,推导的方程组,比如像超越方程这种,可以用matlab进行直接求值,这在一些需要仿真数据的场景比较有用. 环境:mat ...
- 求解非线性方程组的牛顿法c语言,牛顿下山法求解非线性方程(组)(C实现)...
1.算法描述 (1)符号说明与基本假设 对于非线性方程组: (1) 引入向量: 可将(1) ...
- 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; ...
- c语言割线法解非线性方程,用割线法解非线性方程组.doc
用割线法解非线性方程组 自动化学院1011203050 陈晓祺 拟牛顿法解下列方程组 先将拟牛顿法的程序代码如下 Function[r,m]=mulVlineF,x0,A, eps) Format l ...
- 非线性方程组求解方法,神经网络的非线性函数
1.rbf神经网络原理 rbf神经网络原理是用RBF作为隐单元的"基"构成隐含层空间,这样就可以将输入矢量直接映射到隐空间,而不需要通过权连接. 当RBF的中心点确定以后,这种映射 ...
- python fsolve说明_python用fsolve、leastsq对非线性方程组求解
背景: 实现用python的optimize库的fsolve对非线性方程组进行求解.可以看到这一个问题实际上还是一个优化问题,也可以用之前拟合函数的leastsq求解.下面用这两个方法进行对比: 代码 ...
最新文章
- 你的团队能通过电梯测试吗?
- svn还原文件中去掉已经删除的文件
- Collections练习之按照字符串长度进行排序
- oracle数据库开多线程,学习笔记:Oracle表数据导入 DBA常用单线程插入 多线程插入 sql loader三种表数据导入案例...
- Bootstrap 可切换的标签导航
- js取thymeleaf值_26、模板引擎thymeleaf
- 随机读写工具,手写,百度云源码直接下载
- GitHub 使用教程图文详解(转)
- zend studio html插件安装,Zend Studio使用教程:将Zend Studio作为插件安装
- android dp转px的公式_dp与px的转换对照表
- Ubuntu设置静态IP地址的几种方法
- Java 服务器版超市管理系统,基于JAVA的超市进销存管理系统.doc
- zblog php建站教程_开源建站程序ZBlogPHP博客安装图文教程
- 3-3 uniapp、HTML5+、Native.js 功能代码汇总
- [攻防世界]getit
- MFC几种给对话框添加背景图的方法
- 【Python脚本抢红包】用Python写了一个自动抢红包工具,今年过年准备大干一场
- brew对redis的使用
- QQ浏览器X5内核问题汇总 转
- 【IPv6+燎原系列—第2期】IPv6+逐梦赛场,勇夺属于通信人的“金牌”