微分方程相关问题的求解程序的分析_第1页
微分方程相关问题的求解程序的分析_第2页
微分方程相关问题的求解程序的分析_第3页
微分方程相关问题的求解程序的分析_第4页
微分方程相关问题的求解程序的分析_第5页
已阅读5页，还剩7页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、微分方程相关问题的求解程序的分析求解常微分方程的常用MATLAB方法有三类：一、有精确解命令格式：dsolve(odefun)例1：求微分方程 dsolve(Dy=x+x*y,x) ans = -1+exp(1/2*x2)*C1 dsolve(Dy=x+x*y) ans = -1+exp(x*t)*C1 dsolve(Dy=t+t*y)ans =-1+exp(1/2*t2)*C1的通解（注意括号里的变量）例2：求 r=dsolve(D2y=-a2*y,y(0)=1,Dy(pi/a)=0)r =cos(a*t)的解例3：求微分方程组 x,y=dsolve(Dx=y,Dy=-x)x =C1*sin

2、(t)+C2*cos(t) y =C1*cos(t)-C2*sin(t)的解练习：求解 dsolve(Dy=-2*y+2*t2+2*t,y(0)=1) ans = t2+exp(-2*t)二、常微分方程的初值问题的数值解t,y=ode45(odefun,tspan,y0)参数说明如下：odefun:用以表示f(x,y)函数，tspan表示t的取值。如果ode45()解法速度很慢，可以采用ode23s()。还有很多常微分数值解的算法，如果对于微分方程感兴趣的同学可以查阅常微分方程的一些资料及专门介绍MATLAB的比较厚一点的书籍。例：解微分方程组建立M函数文件eg6_3fun.mfunction

3、 f=eg6_3fun(t,x)f(1)=-x(1)3-x(2);f(2)=x(1)-x(2)3;f=f(:);在命令窗口执行：t,x=ode45(eg6_3fun,0:0.1:30,1;0.5)即可。练习：求解微分方程（数值解）求解范围为0,0.5建立m文件function f=edgt,xf=-2*x+2*t2+2*t;命令窗口输入t,x=ode45(edg,0:0.01:0.5,1)即可。三、高阶常微分方程 如果遇到高阶常微分方程则可以先把高阶常微分方程转换为一阶常微分方程组。一个任意的高阶常微分方程可以表述如下：上式可以化为如下一阶常微分方程组：设y1=y,y2=y,yn=y(n-1)例：求解微分方程组已知当时，解：首先引入辅助变量化为一阶方程组则得到建立M函数文件eg6_4fun.mfunction f=eg6_4fun(t,x)f(1)=x(2);f(2)=x(3);f(3)=-3*x(1)*x(3)+2*x(2)2-x(4);f(4)=x(5);f(5)=-2.1*x(1)*x(5);f=f(

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论