《数学实验》实验报告——用Mathematica软件解微分方程.pdfVIP

《数学实验》实验报告

班级****学号****姓名****成绩

试验试验试验

用Mathematica软件解微分方程自选试验

内容类别时间

试验问题：

例1求解下列微分方程：

22

1)y(1y)(2y)

yz

2)

zy

x

3)

y3y3yy(x5)e

22

例2求常微分方程y′=x+y，满足初始条件y（0）=0的数值解

5t

例3求函数t和esint的拉氏变换

例4用拉氏变换解微分方程：

x

+3x″+3x′+x=1满足条件x″(0)=x′(0)=x(0)=0的解。

例5求微分方程组：

x2xy2y0

t

xy2x2e

′

满足条件x（0）=3，x（0）=2，y（0）=0的特解。

试验目的：

1.掌握用Mathematica求解微分方程（组）；（例1）

2.掌握用Mathematica求微分方程的数值解；(例2)

3.掌握用Mathematica进行拉氏变换与逆变换，并把这种变换应用到方程（组）中去。

（例3，例4，例5）

问题分析（可含问题的背景、相关知识、数学建模与求解的方法等）：

１

问题的背景：以上向个实例都是围绕常微分方程展开的，这些例子是来自《方程常微分》

第三版，作者：王高雄……。除了一些一般形式的方程（组）之外，还包含拉普拉斯变化

的方程组。

求解方法：如果用人工的方法来计算，只能根据方程的特征找到相应的解题方法，而且在

计算方面很容易出错。

如果用用Mathematica，很多种不同类型的题目可以用同一个函数来求解。很方便，效率很

高。

试验步骤（根据问题分析及试验目的所计划的试验步骤）：

1.在《方程常微分》找一些比较有代表性的题目

2.根据所选的题目，找到用Mathematica所需要的知识。

3．用Mathematica软件进行求解。

4．将求得的结果进行形式的变化，与书中一些结果进行对比。

试验过程（含详细试验步骤、程序清单及异常情况记录等）

例1求解下列微分方程：

22

1)y(1y)(2y)

In[1]：=DSolve[(y[x]^2)(1-y[x])(2-y[x])^2,y[x],x]

Out[1]=

书中结果为：yxc1/(xc)，其中c为任意常数。

yz

2)

zy

In[1]：=DSo