Matlab实验报告五(微分方程求解Euler折线法).doc_第1页
Matlab实验报告五(微分方程求解Euler折线法).doc_第2页
Matlab实验报告五(微分方程求解Euler折线法).doc_第3页
Matlab实验报告五(微分方程求解Euler折线法).doc_第4页
全文预览已结束

下载本文档

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

文档简介

数学与信息科学系实验报告实验名称 微分方程求解 所属课程 数学软件与实验 实验类型 综合型实验 专 业 信息与计算科学 班 级 学 号 姓 名 指导教师 一、实验概述【实验目的】 熟悉在Matlab环境下求解常微分方程组和偏微分方程组的方法，掌握利用Matlab软件进行常微分方程组和偏微分方程组的求解。【实验原理】 1.dsolve(equ1,equ2,.):matlab求微分方程的解析解。2. simplify(s):对表达式S使用MAPLE的化简规则进行化简。3. x,y=dslove(方程1，方程2，.初始条件1初始条件2,.自变量)：用字符串方程表示，自变量缺省值为t.4. ezplot(x,y,tmin,tmax):符号函数的作图命令。【实验环境】 MatlabR2010b二、实验内容问题1. 求微分方程组在初始条件下的解，并画出函数的图像. 1.分析问题本题是根据初始条件求微分方程组的特解，并根据t的范围画出函数的图形。2.问题求解syms x y tx,y=dsolve(Dx+x+y=0,Dy+x-y=0,x(0)=1,y(0)=0,t)x=simple(x)y=simple(y)ezplot(x,y,0,0.5);axis auto3.结果x =exp(2(1/2)*t)/2 + 1/(2*exp(2(1/2)*t) - (2(1/2)*exp(2(1/2)*t)/4 + 2(1/2)/(4*exp(2(1/2)*t)y =2(1/2)/(4*exp(2(1/2)*t) - (2(1/2)*exp(2(1/2)*t)/4x =cosh(2(1/2)*t) - (2(1/2)*sinh(2(1/2)*t)/2y =-(2(1/2)*sinh(2(1/2)*t)/24.结论及分析通过实验，结论正确，证明分析无误。问题2. 用Euler折现法求解常微分方程的数值解（步长），求解范围，并作出去图像. 1.分析问题本题是用Euler折线法根据已知条件求解微分方程组的数值解，并作出它的图形。2.问题求解clearf=sym(y-12*x2/y3);a=0;b=2;h=0.001;n=(b-a)/h+1;x=0;y=1;szj=x,y;for i=1:n-1 y=y+h*subs(f,x,y,x,y); x=x+h; szj=szj;x,y;endSzj;plot(szj(:,1),szj(:,2)3.结果4.结论及分析通过实验，结论正确，证明分析无误。三、实验小结 这次实
人人文库> 全部分类> 应用文书 > 作业报告

温馨提示

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

最新文档

评论

0/150

提交评论