东南大学数学实验报告(1)

高等数学数学实验报告实验人员：院(系) 土木工程学院学号05A11210 姓名李贺__

实验地点：计算机中心机房

实验一空间曲线与曲面的绘制

一、实验题目：(实验习题1-2)

利用参数方程作图，做出由下列曲面所围成的立体图形：

2 2 2 2

⑴ Z 1 X y，x y X 及xOy平面；

⑵ z xy,x y 1 0 及z 0.

二、实验目的和意义

1、利用数学软件Mathematica绘制三维图形来观察空间曲线和空间曲面图形的特点，以加

强几何的直观性。

2、学会用Mathematica绘制空间立体图形。

三、程序设计

空间曲面的绘制

x x(u, V)

y y(u,v),u [u min , max ],V [V min , V max ]

作参数方程z z(u,v)所确定的曲面图形的Mathematica命令

为:

ParametricPlot3D[{x[u,v],y[u,v],z[u,v]},{u,umi n,umax}.

{v,vmi n,vmax}, 选项]

⑵

t2 = ParametricPlotJD [{u f 1 v}, [u^ ・0・§尸1}^ (v, 0F 1},

HxegLabel {"x" 11 y" J1 z" }. PlotPolnts t 5B,

Dlspla^unction -> Identity」：

t3 = ParametricPlotSD[{u f 0}* (u, -U.J5』1}^

{v z-0.5, 1} f AxesLabel {"x" 11y" 11 z"

PlotPoints 50, Display1 unction — Identity]: Slinw[tl z t2, t3 f DisplayFunction -> SDlsplajfunction]

四、程序运行结果

⑴

(2)

五、结果的讨论和分析

1、通过参数方程的方法做出的图形，可以比较完整的显示出空间中的曲面和立体图形。

2、可以通过mathematica软件作出多重积分的积分区域，使积分能够较直观的被观察。

3、从(1)中的实验结果可以看出，所围成的立体图形是球面和圆柱面所围成的立体空间

4、从(2) 中的实验结果可以看出围成的立体图形的上面曲面的方程是z xy，下底面的方程是z=o,右边的平面是x y 1 0。

实验一空间曲线与曲面的绘制

一、实验题目：(实验习题1-3 )

观察二次曲面族z x 2 y 2 kxy 的图形。特别注意确定k 的这样一些值,当过这些值时,曲面从一种类型变成了另一种类型。

二、实验目的和意义

1. 学会利用Mathematica 软件绘制三维图形来观察空间曲线和空间曲线图形的特点。

2. 学会通过表达式辨别不同类型的曲线。

三、程序设计

这里为了更好地分辨出曲线的类型,我们采用题目中曲线的参数方程来画图,即

z r 2 kr 2cos t sin t

输入代码：

ParametricPlot3D

[{r*Cos[t] , r*Sin[t]F2+ kW2*Cos[t]*Sin[t]},

{t, 0, 2*Pi}, {r, 0, 1} ，PlotPoints -> 30]

4、从(2) 中的实验结果可以看出围成的立体图形的上面曲面的方程是z xy，下底面的方式中k 选择不同的值：-4 到4 的整数带入

四、程序运行结果

k=4：

k=3：

k=2：

k=1:

k=0:

k=-1:

k=-2:

k=-3:

k=-4:

五、结果的讨论和分析

k 取不同值，得到不同的图形。我们发现，当|k|<2 时，曲面为椭圆抛物面；当|k|=2 时，曲面为抛物柱面；当|k|>2 时，曲面为双曲抛物面。

实验二无穷级数与函数逼近

一、实验题目：（实验习题2-2 ）

改变例2中m及x o的数值来求函数的幕级数及观察其幕级数逼近函数的情况。

、实验目的和意义

1. 利用Mathematica显示级数部分和的变化趋势。

2. 学会如何利用幕级数的部分和对函数进行逼近以及函数值的近似计算。

三、程序设计

若函数f(x) (1 x)m能展开成X-X。的幕级数(这里不验证)，则根据函数展开为幕级数的

展开公式，其展开式为f(x) 丄凹(x X o)n。因此首先定义f(x)的n阶导数的函数g(n, no n!

X。)，最后再构成和式即得f (x)的幕级数展开式。用Mathematica观察幕级数部分和逼近函数的情况。

m=- 2, x0=2 时

输入如下命令：

m=- 2;

f [ x_] : =( 1 +x) A m;

x0 = 2;

g[ n_, x0_]::= D[f[x],{x,n}]/. x x0;

s [ n_, x_] : =Su m[ g[k,x0] * ( x- x0) A k , { k, 0, n}]; _ 一k!

t=Tabl e[ s [n ,x] ,{n ,20}];

p1=PI ot [ Eval uat e[ t ] , { x, -1/2, 1/ 2}];

p2 = PI ot [ ( 1+x) A^ {x, - 1/ 2, 1/ 2} , PI ot St y I e RGBCol or [ 0, 0, 1]];

Show[ p 1 , p2]

四、程序运行结果