- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Mathematica__常微分方程、拉氏变换和级数实验.doc
常微分方程、拉氏变换与级数实验
[学习目标]
会用Mathematica求解微分方程(组);
能用Mathematica求微分方程(组)的数值解;
会利用Mathematica进行拉氏变换与逆变换;
能进行幂级数和傅里叶级数的展开。
常微分方程(组)
Mathematica能求常微分方程(组)的准确解,能求解的类型大致覆盖了人工求解的范围,功能很强。但不如人灵活(例如在隐函数和隐方程的处理方面),输出的结果与教材上的答案可能在形式上不同。另外,Mathematica求数值解也很方便,且有利于作出解的图形。在本节中,使用Laplace变换解常微分方程(组)的例子也是十分成功的,过去敬而远之的方法如今可以轻而易举的实现了。
求准确解的函数调用格式如下:
DSolve[eqn,y[x],x] 求方程eqn的通解y(x),其中自变量是x。
DSolve[{eqn,y[x0]= =y0},y[x],x] 求满足初始条件y(x0)= y0的特解y(x)。
DSolve[{eqn1,eqn2,…},{y1[x],y2[x],…},x] 求方程组的通解。
DSolve[{equ1,…,y1[x0]= =y10,…},{y1[x],y2[x],…},x] 求方程组的特解。
说明:应当特别注意,方程及各项参数的表述方式很严格,容易出现输入错误。微分方程的表示法只有通过例题才能说清楚。
解下列常微分方程(组):
(1),(2), (3) ,
(4)的通解及满足初始条件y(0)=0,z(0)=1的特解。
解:In[1]:=DSolve[y′[x]= =2y[x]/(x+1)+(x+1)^(5/2),
y[x],x]
Out[1]=
In[2]:=DSolve[y′[x]= =(1+y[x]^2)/((x+x^3)y[x]),y[x],x]
Out[2]={{}, {}}
In[3]:=DSolve[{y′[x]= =z[x],z′[x]= = -y[x]},
{y[x],z[x]},x]
Out[3]={{y[x]→C[1]Cos[x]+ C[2]Sin[x],
z[x]→C[2]Cos[x]- C[1]Sin[x]}}
In[4]:=DSolve[{y′[x]= =z[x],z′[x]= = -y[x],y[0]= =0,z[0]= =1},
{y[x],z[x]},x]
Out[4]={{y[x]→Sin[x],z[x]→Cos[x]}}
提示:认真观察上例,可以从中学习输入格式,未知函数总带有自变量,等号用连续键入两个等号表示,这两点由于不习惯会出错!导数符号用键盘上的撇号,连续两撇表示二阶导数,这与习惯相同。自变量、未知量、初始值的表示法与普通变量相同。
说明:输出结果总是尽量用显式解表出,有时反而会使表达式变得复杂,这与教科书的习惯不同。当求显式解遇到问题时,会给出提示。通解中的任意常数用C[1],C[2],…表示。
求解下列微分方程:
(1),(2),(3)。
解:In[1]:=DSolve[+3y″[x] +3y′[x] + y[x] = =(x - 5)Exp[-x],
y[x],x]
Out[1]={{
}}
In[2]:=Simplify[%]
Out[2]={{}}
In[3]:=DSolve[x^2 + y′[x]^2 = = 1,y[x],x]
Out[3]={{},
{}}
In[4]:=DSolve[Sqrt[y′[x]] = = x y[x],y[x],x]
Out[4]={{}}
说明:由以上可以看出对方程的类型并无限制,但是输出的答案未必符合习惯,例如第一个方程的答案需要化简,有时即使化简后也未必与教材上的答案一致。
求微分方程的通解。
解:In[1]:=DSolve[y′[x]+2x y[x]= = x E^(-x^2),y[x],x]
Out[1]={{y[x]→}}
这就是所给微分方程的通解。式中的C
您可能关注的文档
最近下载
- 新外研版高中英语必修三Unit 3 The world of science Period 1 Starting out & Understanding ideas说课稿.doc
- 项目回款登记表.xlsx VIP
- 《药品记录与数据管理要求》解读.pptx
- 购物广场商场超市消防安全知识培训PPT.pptx VIP
- 油田运输有限公司HSE管理的体系程序文件.doc
- 人教版2022-2023学年二年级数学第一学期期中测试卷(含答案).docx
- 全科医师培训课件.pptx
- 汉字听写大会试题库(附答案).pdf VIP
- 专科医生如何带教全科医生课件.pptx
- 2024小学生汉字听写大会竞赛试题库(附答案).docx VIP
文档评论(0)