1 / 14
文档名称:

东南大学Mathematica数学实验报告.doc

格式:doc   大小:570KB   页数:14页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

东南大学Mathematica数学实验报告.doc

上传人:pppccc8 2019/1/29 文件大小:570 KB

下载得到文件列表

东南大学Mathematica数学实验报告.doc

文档介绍

文档介绍::..Mathematica〔学实验报告姓名:于润滟学号:04213704成绩:实验七:空间曲线与曲面的绘制实验目的:学****利用Mathematica绘制三维图形来观测空间曲线和空间曲面图形的特点,并学****通过表达式判断不同的曲线类型。:观察二次曲面族z= +形。特别注意确定A的这样一些值,当々经过这些值时,曲面从一种类型变成了另一种类型。解:,y=rsinz,则二次曲曲族的方程口j变力z=r2+/r2cosfsiiV。输入以卜命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],r八2+k*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints^30]并赋予k不同的值:①k=-4时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],r八2+(-4)*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints一30]运行后符到图像:②k=-3时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],rA2+(-3)*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints->30]运行后得到图像:-10 - 00 05 ③k=-2时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],rA2+(-2)*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints->30]运行后得到阁像:-10 - 00 ④k=-1时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],r八2+(-1)*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints一30]运行后得到图像:⑤k=0时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],rA2+0*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints-30]运行£;得到图像:⑥k=1时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],rA2+1*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints-30]运行盾得到图像:10 05000510⑦k=2时:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],rA2+2*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints->30]运行后得到图像:-⑧k=3吋:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],rA2+3*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints->30]运行盾得到图像:⑨k=4吋:输入命令:ParametricPlot3D[{r*Cos[t],r*Sin[t],r八2+4*r*r*Sin[t]*Cos[t]},{r,0,1},{t,0,2Pi},PlotPoints->30]运行后得到图像:-- 00 o5 |o由上面不同k值对应的不同二次曲面图形可得:当0^|々|<2时,对应的二次曲面图形是椭圆抛物面;当|々|=2时,对应的二次曲面图形是抛物柱面;当|々|>2时,对应的二次曲面是双曲抛物面。实验八:无穷级数和函数逼近实验0的:学>」用Mathematica显示级数部分和的变化趋势;学会如何利用幂级数的部分和对函数进行逼近以及函数值的近似计算;展示Fourier级数对周期函数的逼近情况。题目:观察函数/W二-1,-;^1<0;1,0<1<兀展成的Fourier级数的部分和逼近/O)的情况。解:根据Fourier系数公式可得:“0=—£/{x)dx=1+—»71“ 2sinnxdx+f(-xsinz?%)^],an=—[£cosnxdx+ (-xcosnx)dx],bn=—[£故输入以下命令,从输出的图形屮观察Fourier级数的部分和逼近的情况::f[xJ:=Which[-2Pi<x<-Pi,1,-Pi<x<0,-x,0<x<Pi,1,Pi<x<2Pirx+2Pi];a[nJ:=(lntegrate[-x*Cos[n*x],{x,-Pi,0}]+lntegrate[Cos[n*x],{x,0,Pi}])/Pi;b[nJ:=(lntegrate[-x*Si