- 1、本文档共28页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Part 6——微分和差分方程
这次教程有点多,望大家耐心看哦!
第一章
1.1 代数方程求解
下面主要介绍几类常用的方程求根方法。
(1)多项式——roots
调用示例:
r = roots(p)
其中,p为多项式的系数,从高次到低次顺序,r为对应的多个根。
如多项式:,求解程序如下:
p = [1 -6 -72 -27];
r = roots(p)
r =
12.1229
-5.7345
-0.3884
注意:通过其它函数也可以多项式的根,但是不能求出所有的根,而roots可以求出多项式所有的根。下面将会看到。
参考函数:poly, residue
(2)一元函数——fzero
调用示例:
x = fzero(fun,x0)
x = fzero(fun,x0,options)
[x,fval] = fzero(...)
其中,fun为待求函数的名称,x0为初始值,x为解,fval为最优解对应的值(对于求根来说,理想情况下为0)。
如: ,求根程序如下:
f = @(x)x.^3-2*x-5;
z = fzero(f,2)
再看看(1)中的示例,程序
f = @(x)x.^3-6*x.^2-72*x-27;
z1 = fzero(f,10)
z2 = fzero(f,-4)
z3 = fzero(f,0)
很容易发现,对于不同的初始值,最后的根也不一样,最后的结果分别为:12.1229、-5.7345、-0.3884,而这三个分别为(1)中roots得到的解。实际上fzero函数的实现方式是以某个初始点来迭代得到方程解的,因为初始值的选取对最终结果影响很大。下面说到的牛顿迭代之类的算法也是类似情况,对初始值的选取很敏感,而智能算法对初始值的选取并不敏感。
注意:f = @(x)...是匿名函数的表达方式,即该函数没有名字,但是通过@来指定变量,或者待求变量。
参考函数:fminbnd, optimset, function_handle (@), Anonymous Functions
小技巧:快捷注释与注释消除键,选取待注释语句,然后ctrl+R为注释,ctrl+T为消除注释。
(3)非线性方程组——fsolve
调用示例:
x = fsolve(fun,x0) x = fsolve(fun,x0,options) [x,fval] = fsolve(fun,x0)
其中,fun为函数名,x0为初始值,options为参数设置(迭代算法、步长、精度等),x为最优解,fval为最优解对应的最优值。
示例如下:
x0 = [-5; -5]; % 随机初始值
options=optimset(Display,iter); % 显示每次迭代的过程
[x,fval] = fsolve(@myfun,x0,options) % 调用优化函数
其中,myfun如下:
function F = myfun(x)
F = [2*x(1) - x(2) - exp(-x(1));
-x(1) + 2*x(2) - exp(-x(2))];
迭代过程:
最后的解为:
x =
0.5671
0.5671
fval =
1.0e-006 *
-0.4059
-0.4059
注意:该函数求解性能很强,大家可以随意设置初始值x0看看结果如何,理论上结果是不变的,也就是该函数使用的迭代算法对初始值不敏感。
(4)符号方程——solve
调用示例:
solve(eq) solve(eq,var) solve(eq1,eq2,...,eqn) g = solve(eq1,eq2,...,eqn,var1,var2,...,varn)
其中,eq为待求解的符号方程,默认的求解符号是x;通过var设置带求解的符号变量;求解符号方程组的话,需要依次添加每个方程,以及带求解的所有符号变量。g是存放结果的结构体变量,通过g.x调用显示结果。
示例如下:
%% 求解 a*x^2 + b*x + c = 0 关于x的符号解
solve(a*x^2 + b*x + c)
%% 求解 a*x^2
文档评论(0)