《数学实验》上机内容

1
《数学实验》上机内容
实验题目
实验一 解方程和方程组与极限运算
(一)．实验类型：验证型
(二)．实验类别：基础实验
(三)．每组人数：1
(四)．实验要求：选修
(五). 实验学时：4个学时
(六)．实验目的：（1）掌握Mathematica软件的计算器功能；（2）学会使用Mathematica软
件求各种类型方程（或方程组）的数值解和符号解；（3）通过本实验深刻理解极限概念；（4）
学习并掌握利用Mathematica求极限的基本方法。

(七)【预备知识】
（1）方程（或方程组）代数解法的基本理论，函数的零点，方程（或方程组）的解及
数值解；
（2）本实验所用命令：
●
用“= =”连接两个代数表达式构成一个方程
●
求方程（组）的代数解：
Solve[方程或方程组，变量或变量组]
●
求方程（组）的数值解：
NSolve[方程或方程组，变量或变量组]
●
从初始值开始搜索方程或方程组的解：
FindRoot[方程或方程组，变量或变量组初值]
●
在界定范围内搜索方程或方程组的解：
FindRoot[方程或方程组，变量或变量组范围]
●
绘图命令：
Plot[表达式，{变量，上限，下限}，可选项]
●
微分方程求解命令：
DSolve[微分方程, y[x], x]
（3）极限、左极限、右极限的概念；
2

（4）本实验所用Mathematica有关命令：
●
Limit[expr, x->x0] 求表达式在0xx时的极限
●
Limit[expr,x->x
0
,Direction -> 1] 求左极限

●
Limit[expr,x->x
0
,Direction ->-1] 求右极限

(八)．实验内容

（1）计算54564546；4567646545。
（2）对于方程0342234xxx，试用Solve和NSolve分别对它进行求解，并比较得
到的结果，体会代数解即精确解与数值解的差别。
（3）先观察函数CosxSinxxf)(的图形，然后选择一个初始点求解，并且根据图形
确定在某个区间中搜索它的零点。

（4）求方程组222111cybxacybxa的解，然后代入系数和常数项的一组初值，并求解。

（5）求微分方程xxyxyxye)(2)(3)(的通解。
（6）用 Mathematica软件计算下列极限：
（1）1233limnnnn； （2）xπxtanlim2； （3）xπxtanlim2；

（4）xxxxx3333lim； （5）nnznzn22lim； （6）210)sin(limxxxx；
（7）xxax1)1(lim0；（8）222limlimyxyxxy；（9）yxyyxxy3252223limlim；
（10）yxyyxyx3252232limlim；（11）222limlimyxyxyx；（12）)1sin(lim0xx。
3

实验二 积分运算与微分基本运算及函数的幂级数展开
(一)．实验类型：验证型
(二)．实验类别：基础实验
(三)．每组人数：1
(四)．实验要求：选修
(五). 实验学时：4个学时

(三)．实验目的：（1）通过本实验加深理解积分理论中分割、近似、求和、取极限的思想方
法；（2）学习并掌握二重积分及线性积分的计算方法；（3）学习常用积分命令。（4）掌握求
函数的导函数和偏导数方法；（5）学会使用Mathematica软件进行函数的幂级数展开。

(四)．【预备知识】
（1）定积分的概念、几何意义，二重积分的概念、二重积分化为定积分的过程及其计算方
法；
（2）本实验所用Mathematica有关命令：
●
无限积分：Integrate[f,x]
●
定积分：Integrate[f,{x，上限，下限}]
（3）函数的导函数、偏导数以及函数的幂级数展开式；
（4）本实验所用的Mathematica函数提示：
（a）求导数（或偏导数）
●
D[表达式F,x] 求F对于变量x的导数；
●
D[表达式F,x1,x2,...] 按顺序求F关于x1，x2,„的偏导数；
●
D[表达式F,{x,n}] 求F对x的n阶导数。
（b）幂级数展开
●
Series[表达式F,{x,x0,n}] 求F关于变量x在x0的n阶泰勒展式。
(五)．实验内容：

（1）求函数32)sin(xxaf的原函数；
（2）求xaxnd；
（3）求10dxaxn；
（4）求10122ddxxyxyx；
（5）求xyyxx00dcosd。
4

（6）求出被积函数F(x)=5312xxx的原函数和导函数，并画出被积函数、原函数和导
函数的图形，试分辨出哪一条曲线属于哪个函数。
（7）对函数sinx在0点展开10阶和20阶，并以图形方式对比展开的结果和sinx的差
别，并分析阶数高的展式对于原来函数的逼近程度是否优于阶数低的展式。
5
实验三 放射性废料的处理问题
(一)．实验类型：综合型
(二)．实验类别：基础实验
(三)．每组人数：1
(四)．实验要求：选修
(五). 实验学时：4个学时
(六)．实验目的：巩固和理解微分方程理论及其应用。

(七)．预备知识：常微分方程理论和Mathematica解方程的命令。
(八)．【实验内容与要求】

美国原子能委员会以往处理浓缩放射性废料的方法，一直是把它
们装入密封的圆桶里，然后扔到水深90多米的海底。生态学家和科
学家们表示担心，怕圆桶下沉到海底时与海底碰撞而发生破裂，从而
造成核污染。原子能委员会分辩说这是不可能的。为此工程师们进行
了碰撞实验，发现当圆桶下沉到海底时的速度超过12.2 m/s，圆桶与
海底碰撞会发生破裂。为避免圆桶碰裂，需要计算圆桶沉到海底时的
速度是多少？这时已知圆桶重为239.46 kg，体积为0.2058 m3,海水密
度为1035.71 kg/m3。如果圆桶下沉到海底时的速度小于12.2 m/s，就
说明这种方法是可靠的；否则就要禁止用这种方法来处理放射性废
料。假设水的阻力与速度大小成正比，其正比例常数为0.6。

（1）根据问题建立数学模型。
（2）根据数学模型求解的结果，判断这种处理废料的方法是否合理？
6

实验四 路程估计问题
(一)．实验类型
：综合型

(二)．实验类别
：基础实验

(三)．每组人数：1

(四)．实验要求：选修
(五). 实验学时：4个学时
(六)．实验目的：能用数学软件进行数据拟合。

(七). 预备知识：多元函数的极值求法；线性拟合的最小二乘法原理。

(八)【实验内容与要求】
外出旅行或行军作战等，都可能涉及到两地路程的估计问题。当
身边带有地图时，这似乎是件很容易的事。然而，从地图上量出的距
离却是两地的直线距离d，你能由此估计出两地的实际路程s吗？建
立s关于d的模型：)(dfs。
（1）要确定s与d的近似函数关系，必须收集若干s及与之相应的d的具体数据，通过
分析找出规律。这里将《中国地图》中量得四川省彭州市到其他几个城市的直线距离，并按
比例尺（1cm为20km）进行转换，以及从到汽车站了解到的对应的实际路程的有关数据列
于表1。
表1 城市间直线距离和实际路程

彭
州
市

成都 郫县 都江堰 什邡 德阳 新繁 广汉 温江 崇

庆
→

地图直线
距离(cm)
1.8 1.08 1.55 1.32 2.3 0.75 1.64 1.7 2.38

地图转换
距离
d(km)
36 21.6 31 26.4 46 15 32.8 34 47.6

实际路程
s(km)
42 30 58 43 68 16 43 50 65

（2）启动数学软件，将上表中d与s两组数据，按拟合时所需形式输入。
7

（3）画出数据散布图，观察它们是否大致在一条直线附近。
（4）进行直线拟合，并在同一图中显示拟合直线与数据点。观测拟合情况，并记下所
得到的模型（称为经验模型）。
（5）在只作粗略估计的情况下，为便于计算，若将上面得到的模型修改成bds5.1（简
单模型）行吗？根据表中数据，取b=3，试画出简单模型与样本数据点的图形，并与（4）
所得到的图形相对照。
（6）试计算由两个模型得到的估计值与实际值的差（残差），以大致观测一下两个模型
的差异。在只作粗略估计的前提下，你愿意用哪个模型？