Matlab 数值计算----二分法求非线性方程组

bisect.m函数（二分法的实现算法）

function[x_star,index,it]=bisect(fun,a,b)
%求解非线性计算方程的二分法，其中，fun(x)为需要求根的函数;
%a,b为初始区间的端点;
%ep为精度，当(b-a)/2<ep时，算法能终止计算，
%缺省值为1e-5;
%当x_star迭代成功时，输出方程的根
%当x_start迭代失败时，输出两端点的值；
%index为指标变量，当index=1时，表明迭代成功，
%当index=0时，表明初始区间不是有根区间；
%it为迭代次数
if nargin<4ep=1e-5;
end
fa=feval(fun,a);fb=feval(fun,b);
if fa*fb>0x_star=[fa,fb];index=0;it=0;return;
end
k=1;
while abs(b-a)/2>=epx=(a+b)/2;fx=feval(fun,x);if fx*fa<0b=x;fb=fx;elsea=x;fa=fx;endk=k+1;
end
x_star=(a+b)/2;index=1;it=k;

test12.m函数（给定一个非线性方程组和区间进行二分法的计算测试）

fun=inline('x^3-x-1');
[x_star,index,it]=bisect(fun,1,2)

结果：

>> test12

x_star =

1.3247

index =

it =

Matlab 数值计算----二分法求非线性方程组相关推荐

二分法求非线性方程组Java_用C#编写二分法解一元非线性方程
在计算方法里有一种求一元非线性方程的解法,叫做二分法. 简单介绍如下: 函数F(x)在区间[a,b] 上连续,假定在区间内有唯一实根,计为x* 二分法的思想为:首先确定有根区间,将区间二等分,通过判断 ...
求非线性方程组的最小二乘解的广义逆法C实现
求非线性方程组的最小二乘解的广义逆法 #include "math.h"#include "stdlib.h"#include "6gmiv.c&qu ...
二分法求方程根matlab,matlab用二分法求方程的正根,要求误差小于0.0005
matlab编程题:用二分法求方程x^3-3*x-1=0的根先建立二分法的fun.m文件,代码如下:functionfun(a,b,e)%f是自定义的函数%a为隔根区间左端点,b为隔根区间右端点,e ...
蒙特卡洛法求非线性方程组
蒙特卡罗法的概念及应用蒙特卡洛法 (又称统计试验法)是描述装备运用过程中各种随机现象的基本方法,而且它特别适用于一些解析法难以求解甚至不可能求解的问题,因而在装备效能评估中具有重要地位. 用蒙特卡洛 ...
MATLAB符号数学解非线性方程组 -- 眼球反射模型
1.利用符号函数求解非线性方程组 1 定义符号变量 syms 2 方程求解 solve:求解出来的是个解析解,意味着其实不用matlab,自己一直算也能算出一个通式: vpasolve:数值解,可理解 ...
broyden matlab,Broyden方法求解非线性方程组的Matlab实现
<Broyden方法求解非线性方程组的Matlab实现>由会员分享,可在线阅读,更多相关<Broyden方法求解非线性方程组的Matlab实现(7页珍藏版)>请在人人文库网上搜 ...
matlab 非线性方程数值解法,非线性方程组的几种数值解法+matlab源代码
摘要很多领域都有涉及到非线性方程组,例如天气预报,石油地质勘探,电力系统计算等,甚至商业领域也有非线性优化问题,这些问题要从本质上解决就是求出非线性方程组的解.但是目前已知的数值解法并不完善,选择不同 ...
牛顿拉夫森法 matlab,【原创】牛顿-拉夫森迭代求非线性方程组
function [xf,fval,iter]=newtmullt(e,x0,es,maxit,varargin) % written by tubehu % www.MATLABsk ...
matlab怎么求解多元非线性方程组,MATLAB求解多元非线性方程组
我想求解一个非常复杂的方程组,原理很简单,就是解三个方程求出三个未知数a.b.c,但每一个方程都包含exp.log.sqrt函数,不知道能否用MATLAB求出解来,我试着写了程序,但是运行不出来,我第 ...

Matlab 数值计算----二分法求非线性方程组

Matlab 数值计算----二分法求非线性方程组相关推荐

最新文章

热门文章