(03214746秦泽天)14-15-3高等数学数学实验报告

高等数学数学实验报告

实验人员：院（系） 能源与环境学院 学号 03214746 姓名 秦泽天 实验地点：计算机中心机房

实验七

一、实验题目： 2.利用参数方程作图，作出有下列曲面所围成的立体：

(1) 面；

及,12222xOy x y x y x z =+--= (2) .0及01,==-+=z y x xy z

3.观察二次曲面族.22kxy y x z ++=的图形。特别注意确定k 的这样一些值，当k 经过这个值时，曲面从一种类型变成了另一种类型。

二、实验目的和意义

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

合的直观性。

三、程序设计

计算：

2.（1）解得相应曲面部分的参数方程为

⎪⎩⎪⎨⎧-===21sin *cos *r z t

r y t r x ，t r cos =，0=z

2.（2）解得相应曲面部分的参数方程为

⎪⎩⎪⎨⎧===t t r z t

r y t r x sin *cos *sin *cos *2，01sin *cos *=-+t r t r ，0=z

3.函数表达式

;22kxy y x z ++=

程序：

2.(1)

s1=ParametricPlot3D[{r*Cos[t],r*Sin[t],2

1r

},{r,0,1} ,{t,0,2*Pi},PlotRange→{-1,1},AxesLabel→{"x","y","z"}]; s2=ParametricPlot3D[{Cos[t]*Cos[t],Cos[t]*Sin[t],Z},{Z ,0,1},{t,0,2*Pi},PlotRange→{-1,1},AxesLabel→{"x","y"," z"}];

s3=ParametricPlot3D[{x,y,0},{x,-1,1},{y,-1,1},PlotRang e→{-1,1},AxesLabel→{"x","y","z"}];

Show[s1,s2,s3]

2.(2)

s1=ParametricPlot3D[{r*Cos[t],r*Sin[t],r2*Cos[t]*Sin[t ]},{r,0,1},{t,0,2*Pi},PlotRange→Automatic,AxesLabel→{" x","y","z"}];

s2=ParametricPlot3D[{Cos[t]*1/(Sin[t]+Cos[t]),1/(Sin[t ]+Cos[t])*Sin[t],Z},{Z,0,1},{t,0,2*Pi},PlotRange→{0,1} ,AxesLabel→{"x","y","z"}];

s3=ParametricPlot3D[{x,y,0},{x,-1,1},{y,-1,1},PlotRang e→{0,1/2},AxesLabel→{"x","y","z"}];

Show[s1,s2,s3]

3.

For[i=-5,i≤5,i++,Plot3D[z=x^2+y^2+5∗x∗y,{x,−10,10},{y,−10,10}]]

四、程序运行结果

五、结果的讨论和分析

2.在实验过程中，通过对参数方程的求解，绘制出各个曲面的图形以及叠加曲面的图形，在（1）中使用参数方程，简化了运算，使绘制的图形更确切，而在（2）中使用参数方程效果不太理想，对于参数的取值范围掌握不准，绘制显得模糊，考虑到原方程很简单，或许直接绘制会得到更好的效果。

3.观察到图像随k值变化而变化。双击其中一幅图片，可以看到所有图片连续播放的动画。动画中可明显看到一种曲面的渐变方式，以及变为另外一种曲面的过程。

实验八

一、 实验题目：

1.观察函数⎪⎩⎪⎨⎧<≤<≤--=ππx x x x f 0,10,)(展成的傅里叶级数的部分和逼近)(x f 的情

况。

2.观察级数 !1n

n n n ∞=∑ ,的部分和序列的变化趋势，并求和。

二、实验目的和意义

用Mathematica 显示级数部分和的变化趋势；学会如何利用幂级数的部分

和对函数的逼近以及进行函数值的近似计算；展示傅里叶级数对周期函数的逼近情况。

三、程序设计

计算：

1.根据傅里叶系数公式可得：

21)(10ππππ+==

⎰-dx x f a ]cos )cos ([1

00⎰⎰+-=-ππ

πnxdx dx nx x a n ，

]sin )sin ([100⎰⎰+-=-πππnxdx dx nx x b n 程序：

1.

f[x_]:=Which[-2Pi ≤x ≤-Pi,1,-Pi ≤x<0,-1,0≤x

a[n_]:=(Integrate[-x*Cos[nx],{x,-Pi,0}]+Integrate[Cos[nx],{x,0,Pi}])/Pi;

b[n_]:=(Integrate[-x*Sin[nx],{x,-Pi,0}]+Integrate[Sin[nx],{x,0,Pi}])/Pi;

s[x_,n_]:=a[0]/2+Sum[a[k]*Cos[kx]+b[k]*Sin[kx],{k,1,n}];

g1=Plot[f[x],{x,-2Pi,2Pi},PlotStyle →RGBColor[0,0,1],D

isplayFunction→Identity];

m=16;

For[i=1,i≤m,i+=2,g2=Plot[Evaluate[s[x,i]],{x,-2Pi,2Pi} ,––DisplayFunction→Identity];

Show[g1,g2,DisplayFunction→$DisplayFunction]]

2.

s[n_]:=Sum[k!/kk,{k,1,n}];

data=Table[s[n],{n,1,100}];

NSum[n!nn,{n,Infinity}]

四、程序运行结果