非线性方程和常微分方程的解法.pdfVIP

  1. 1、本文档共10页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多

非线性方程和常微分方程的解法

非线性方程和常微分方程的解法

实验8非线性方程和常微分方程的解法

一、实验目的

学会用MATLAB软件求解非线性方程和常微分方程。

二、实验内容与要求

1.非线性方程的整值解

(1)最小二乘法

格式:fsolve(’fun’,x0)%求方程fun=0在估计值x附近的近似解。

[例1.72]求方程xe0的解。

fc=inline(‘x-exp(-x)’);

xl=fsolve(fc,0)

xl=

0.5671

问题1.28:求解方程5xsinx-e0,观察知有多解,如何求之?

先用命令fplot(’5__^2*sin(x)-exp(-x),0]’,[0,10])作图1.13,注意

5__^2*sin(x)-exp(-x),0]中的“[,0]”是作y=0直线,即x轴。方程在[0,

10]区间从图中可看出有4个解,分别在0,3,6,9附近,所以用命令:

fun=inline(’5__^2*sin(x)-exp(-x)’);

fsolve(fun,[0,3,6,9],le-6)

得出结果:

ans=

0.59183.14076.28329.4248

【例1.73】求解方程组x-0.7sinx-0.2cosy=02-x-x

y0.7cosx0.2siny0

先编制函数文件fu.m:

functiony=fu(x)

y(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2));

y(2)=x(1)-0.7*cos(x(1))+0.2*sin(x(2));

y=[y(1),y(2)];

在命令窗口调用fu运算:

x1=fsolve(‘fu’,[0.5,0.5])

x1=

0.52650.5079

(2)零点法

格式:fzero(fun,x0)%求函数fun在x0附近的零点。

说明:估计值x0若为标量时,则在x0附近查找零点,x0=[x1,x2]向

量时,则首先要求函数值fun(x1)fun(x2)0,然后将严格在[x1,x2]区间内零

点,若找不到,系统将给出提示。

非线性方程和常微分方程的解法

【例1.74】求函数f(x)sinx2/xxex4的零点。

fn=inline(sin(x^2)/x+__exp(x)-4);

x=fzero(fn,[1,2])%这里的fn不要加单引号

x=

1.0748

注意:方程解的估计值可用fplot作图看出;用function建立函数文

件fn,求解调用时

fn两边要加单引号,而用inline时fn两边不要加单引号;这两种方

法也可解线性方程组。

⒉代数方程的符号解

格式:gsolve(eq)%求解方程eq0,输入参量eq可是符号表达式或

字符表达式。

gsolve(eq,var)%对eq中指定的变量var求解方程eq(var)0。

gsolve(eq1,eq2,,eqn)%求解方程组eq10,eq20,eqn0。

gsolve(eq1,eq2,eqn,var1,var2,,varn)%对方程组eq1,eq2,eqn中指定

的n

个变量加var1,var2,varn求解

【例1.75】

solve(a__^2+b__+c)

solve(a__^2+b__+c,b)

[x,y]=solve(x+y=1,x-11*y=5)

[a,u,v]=solve(a*u^2+v^2,u-v=1,a^2-5*a+6)

计算结果为:

ans=

文档评论(0)

130****4895 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档