收稿日期:2011-05-27.
基金项目:2011年河南省教育厅人文社会科学研究项目—高职数学教育的隐忧和对策的研究,项目编号:2011—
GH
—
075.
作者简介:王大鹿(1974—),男,河南开封人,硕士,讲师,研究方向:数学建模,数学学习心理与CAI应用.
·高职高专教学·
基于
Mathematica
软件的函数作图及性质分析
王大鹿
(黄河水利职业技术学院基础部,河南开封
475003
)
摘
要:利用Mathematica软件可以轻松作出各种函数的二维或三维图像。探讨通过
Mathematica软件强大的绘图功能学习高等数学如定积分概念
、多元函数的极限、重积分的计算等
内容,有助于提高学生数学分析能力和解决综合应用问题的能力。
关键词:数学软件;Mathematica5.0;函数作图
中图分类号:O313文献标识码:A
文章编号:1006-7353(2011)05-0074-03
0
引言
Mathematic
数学软件由美国
Wolfram
Research
公司开发的一个完全集成环境下的符
号运算系统。其卓越的数值处理、符号运算和作
图功能,已被用来帮助解决许多数学、物理和工程
计算等问题。该软件具有强大的图形输出功能,
将二维、三维图形的绘制和显示集成在
Plot
、
Plot3D
、
ParametricPlot3D
和
Show
等内部函数
中,并有图形修饰选项和多种函数用于特殊图形
的描绘和特殊效果的显示。其丰富的图形表现效
果给高等数学教学提供了很好的可视化素材,使
学生对抽象的高等数学概念获得感性认识,形象
地把握住函数的某些性质,对于加深知识理解,培
养学生的计算能力和分析、解决问题的能力具有
极大的帮助
[1]
。
1
定积分的概念及应用
定积分概念体现了微积分的思想和方法,其
应用几乎涵盖了所有的自然学科。定积分概念既
有其抽象性,也有其具体内容。用
Mathematica
可以将抽象的定积分概念与相关应用变得更具
体,计算也大为简化,降低了学生学习的难度,有
利于激发学生学习的积极性。
根据定积分定义,连续函数在某一闭区间上
的定积分等于该区间上划分的最大长度趋于零时
黎曼和的极限值
[2]
。以函数
y=x
2
在区间上的定
积分为例,编制程序如下:
Clear
[
f
,
x
,
a
,
b
]
f
[
x
_]
=x^2
;
a=0
;
b=1.5
;
m=0
g1=Plot
[
f
[
x
],{
x
,
a
,
b
},
PlotStyle->
{
RGBColor
[
1
,
0
,
0
]},
DisplayFunction->
Identity
]
For
[
j=3
,
j<=50
,
j+=5
,
m=j
;
tt1=
{};
tt2
=
{}
For
[
i=0
,
i<m
,
i++
,
x1=a+i*
(
b-a
)/
m
;
x2=x1+
(
b-a
)/
m
tt1=Append
[
tt1
,
Graphics
[{
RGBColor
[
1
,
0.8
,
1
]
Rectangle
[{
x1
,
0
},{
x2
,
f
[
x2
]}]}]]
tt2=Append
[
tt2
,
Graphics
[{
RGBColor
[
0
,
0
,
0
],
Rectangle
[{
x1
,
f
[
x1
]},{
x2
,
0
}]}]]];
Show
[
tt1
,
tt2
,
g1
,
DisplayFunction->
$DisplayFunction
,
PlotLabel->m
"
intervals
"]]
执行上述命令,可以得到系列图形(取其中一
幅,见图
1
),双击任意一幅图形即可形成动画,当
分割越来越细时,仔细观察小矩形面积之和与曲
47
第24卷第5期高等函授学报(自然科学版)Vol.24No.5
2011年10月JournalofHigherCorrespondenceEducation
(
NaturalSciences
)
2011
边梯形面积之间的关系,有助于理解定积分的概
念和几何意义。
图1
定积分的应用十分广泛,在高等数学中二维
图形面积和旋转体体积的计算为其代表,利用
Mathematic
软件的绘图和数值计算功能,可以帮
助学生快速找到积分的范围,列出算式,得到答
案。这里仅就平面图形的面积举例。
例
1
求由抛物线
y
2
=x
和直线
y=x-2
所
围成图形的面积。
解
首先画出函数图形(如图
2
所示):
图2
Plot
[{
Sqrt
[
x
],
-Sqrt
[
x
],
x-2
},{
x
,
0
,
5
}]
然后求出两条曲线的交点:
Solve
[{
y^2-x==0
,
y-x+2==0
},{
x
,
y
}],得{{
x→1
,
y→-1
},{
x→4
,
y→2
}}
再以
y
为积分变量求面积:
Integrate
[
y+2-y^2
,{
y
,
-1
,
2
}],运行后得
9
/
2
2
多元函数的极限
考察极限
lim
x→0
x→0
xy
x
2
+y
2
,当点
P
(
x
,
y
)沿着直线
y=kx
无限趋近于点(
0
,
0
)时,有
lim
x→0
x→0
xy
x
2
+y
2
=
lim
x→0
kx
2
x
2
+k
2
x
2
=
K
1+K
2
,它是随着
k
值的不同而
改变,所以
lim
x→0
x→0
f
(
x
,
y
)不存在。
图3
取
y=4x
和
y=x
。在
Mathematica
中绘制
出三个图形,通过利用菜单栏里选
Input>>3D
ViewPointSelector
,或者
Shift+Ctrl+V
。在弹
出的对话框里有一个立体图,用鼠标随意转动到
合适的角度,然后把光标点在绘图命令“]”的前
面,加个逗号,然后在
3DViewPointSelector
的
对话框里面点击
Paste
,这样逗号后面就有了
ViewPoint
选项的字样,然后
Shift+Enter
就可
以了直接观察或操作处理了。也可以通过调用一
个画图包来实现,在开始时输入
<<RealTime3D
`
,然后画
3D
图就行了,在生成的图上按住鼠标左
键,随意旋转至理想角度即可。从图形上进行观
察,由于点(
x
,
y
)沿着不同的曲线趋于点(
0
,
0
)时
(见图
3
),
z
趋于不同值,故函数的极限不存在。
借助几何图形,可以更好地理解二元函数极限不
存在的内涵,深刻地理解多元函数极限的实质,熟
悉通过特殊曲线判断多元函数极限不存在的
方法。
3
重积分的计算
在重积分和曲面积分的计算中,确定积分限
是关键。对于比较复杂的图形,要用手工描绘是
非常困难的,但通过
Mathematica
的图形输出功
能,可以清晰地构造各种形体图形和空间各种几
何形体的相互关系,有助于学生及时确定积
分限
[3]
。
例
2
求
Ω
z
2
dv
,其中
Ω
是由曲面
z=x
2
+
57
第24卷第5期高等函授学报(自然科学版)Vol.24No.5
2011年10月JournalofHigherCorrespondenceEducation
(
NaturalSciences
)
2011
y
2
,
y=x
2
,
y=1
,
z=0
围成
解
先分别作出各曲面的图形,以便作出积
分区域的图形,键入
<<RealTime3D
’
s=ParametricPlot3D
[{
u
,
v
,
u^2+v^2
},{
u
,
-2
,
2
},{
v
,
-2
,
2
},
PlotRange->
{
0
,
2
}]
t=ParametricPlot3D
[{
u
,
u^2
,
v
},{
u
,
-2
,
2
},{
v
,
0
,
2
}]
m=ParametricPlot3D
[{
u
,
1
,
v
},{
u
,
-2
,
2
},
{
v
,
0
,
2
}]
n=ParametricPlot3D
[{
u
,
v
,
0
},{
u
,
-2
,
2
},
{
v
,
0
,
2
}]
图4
Show
[
s
,
t
,
m
,
n
,
Boxed->False
,
Mesh→
False
]
运行后可得积分区域图,在生成的图上按住
鼠标左键随意转动到适当的角度(见图
4
),即可
观察到该区域是以
z=x
2
+y
2
为顶,
z=0
为底,
y=x
2
,
y=1
为侧面的一个空间区域,所以
Ω
z
2
dx=
∫
1
-1
dx
∫
1
x
2
dy
∫
x
2
+y
2
0
z
2
dz
因此键入
Integrate
[
z^2
,{
x
,
-1
,
1
},{
y
,
x^2
,
1
},{
z
,
0
,
x^2+y^2
}],运行后得
7712
/
27027
4
曲线拟合
在实践中,常常需要从一组测量数据即根据
给定的
N
个点(
i=1
,
2
,
3
…,
N
)求一条最接近该
组数据点的曲线,以显示这些点总的趋势,此过程
称之为曲线拟合,该曲线方程为回归方程。最小
二乘法是保证具有最佳拟合与回归的常用方法。
但在实际应用时,直接采用最小二乘法手工处理
数据非常繁琐,运用
Mathematic
可以简洁方便
地对数据进行曲线拟合。下面仅就直线拟合进行
举例简单说明。
例
3
设某次实验数据如下;
X1.361.491.721.811.932.15
Y14.0715.0916.8617.3518.4219.96
利用
Mathematica
编制程序如下:
data=
{{
1.36
,
14.07
},{
1.49
,
15.09
},
{
1.72
,
16.86
},{
1.81
,
17.35
},{
1.93
,
18.42
},
{
2.15
,
19.96
}}
a=ListPlot
[
data
,
PlotStyle->PointSize
[
0.02
]]
b=Fit
[
data
,{
1
,
x
},
x
]
c=Plot
[
b
,{
x
,
1
,
3
}]
Show
[
a
,
c
]
该程序第一行将测量数据输到一个名为
data
的列表中,第二行将数据所对应的点画在直角坐
标中,第三行利用函数
Fit
[]将数据进行直线拟
合并输出回归方程,第四行将拟合所得直线画在
直角坐标中,第五行程序将第二、四行所得到的图
形显现在同一个坐标中。该过程输出回归方程为
(输出的拟合直线见图
5
):
3.96736+7.4518x
图5
综上所述,利用
Mathematica
可以方便地描
绘平面或空间各种比较复杂的图形,借助图形可
以更好的观察函数的变化特征,分析其性质,从而
深刻理解相关概念性质的内涵,提高用数学方法
分析问题和解决问题的能力,为将来更好的从事
科研工作奠定良好的基础。
参考文献:
[
1
]章栋恩,许晓革
.高等数学实验
[
M
]
.北京
:高等教育
出版社。2004,7(12):4—11.
[
2
]华东师范大学
.数学分析
[
M
]
.北京
:高等教育出版
社,1991:276—277.
[
3
]陈丽安
.输出高品质MATLAB图形的方法与技巧
[
J
]
.计算机应用研究
,
2002
(
5
):
154
—
155.
67
第24卷第5期高等函授学报(自然科学版)Vol.24No.5
2011年10月JournalofHigherCorrespondenceEducation
(
NaturalSciences
)
2011
|
|