首页 > 其他 > 详细

fslove - Matlab求解多元多次方程组

时间:2019-09-25 22:17:59      阅读:230      评论:0      收藏:0      [点我收藏+]

fslove - Matlab求解多元多次方程组

简介: 之前看到网上的一些资料良莠不齐,各种转载之类的,根本无法解决实际问题,所以我打算把自己的学到的总结一下,以实例出发讲解fsolve。
示例如下:
\[ \begin{cases} 2x_1 - x_2 = e^{ax_1} \-x_1 + 2x_2 = e^{ax_2} \\end{cases} \]
具体的求解过程在后面 点击跳转

1. fsolve的基本使用

调用格式一:

X = fslove(FUN,X0)

功能:给定初值X0,求解方程组的解,X就是返回的解

调用格式二:

X = fsolve(FUN,X0,OPTIONS)

功能:同上,并解决默认参数优化为options指定值

调用格式三:

[X,FVAL] = fslove(FUN,X0,...)

功能:返回X处目标函数值

调用格式四:

[X,FVAL,EXITFLAG] = fslove(FUN,X0,...)

功能:返回EXITFLAG的值,用来描述计算退出的条件,其中EXITFLAG取值和相应的含义如下表。(主要作为判断条件来使用)

EXITFLAG 含义
1 函数fslove收敛于解X处
2 X的变化小于限制
3 残差变化小于限制
4 重要搜索方向小于限制
0 达到最大迭代次数或者评价标准
-1 算法由输出函数终止
-2 算法无法收敛到解的点
-3 信赖域半径太小
-4 线搜索在当前不能充分减少残差

调用格式五:

[X,FVAL,EXITFLAG,OUTPUT] = fslove(FUN,X0,...)

功能:包含OUTPUT的输出

调用格式六:

[X,FVAL,EXITFLAG,OUTPUT,JACOB] = fslove(FUN,X0,...)

功能:返回雅各比矩阵

2.方程求解

(1) 编制函数文件fun.m

编写函数主要用来书写函数的表达式。

function f = fun(x,a,b,c)  % b c可以是随意的参数
f1 = 2*x(1)-x(2)-exp(a*x(1));
f2 = -x(1)+2*x(2)-exp(a*x(2));
f = [f1;f2];

% 也可以写成下面的方式
% f = [2*x(1)-x(2)-exp(a*x(1));-x(1)+2*x(2)-exp(a*x(2))];

(2) 给定函数的参数值和初值(解在周围寻找)

调用求解函数 fslove

>> a = -1;
>> x0 = [-5,-4];
>> [x,FVAL,EXITFLAG,OUTPUT,JACOB] = fsolve(@(x)fun(x,a,1,1),x0);

@(x)fun(x,a,1,1)调用fun函数,函数的参数是a,1,1,求解x的值
执行后调用x返回,也就是X的解。

x =

??0.5671???0.5671

调用FVAL显示在目标解的函数值,可以看出,FVAL越小越接近真实解。

FVAL =

??1.0e-09 *

??-0.4242
??-0.3753

调用EXITFLAG 结合上面的表格可以知道,函数FSOLVE收敛于解X处。

EXITFLAG =

??1

fslove - Matlab求解多元多次方程组

原文:https://www.cnblogs.com/teaurn/p/11587502.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!