试验定积分的近似计算48693.pdf

《试验定积分的近似计算48693.pdf》由会员分享，可在线阅读，更多相关《试验定积分的近似计算48693.pdf（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、.word 格式,专业.专注 .实验一 特殊函数与图形 一、问题背景与实验目的 著名的 Riemann 函数大家都很熟悉了，但是关于它的图像你是否清楚呢？除了最上面那几点，其他都很难画吧？你想不想看看下面那些“挤在一起”的点是怎样分布的呢？还有几何中的马鞍面、单叶双曲面等是怎样由直线生成的，是不是也想目睹一下呢？这些，都离不开绘图 实际上绘图一直是数学中的一种重要手段，借助图形，往往可以化繁为简，使抽象的对象得到明白直观的体现比如函数的基本性质，一个图形常可以使之一目了然，非常有效它虽不能代替严格的分析与证明，但在问题的研究过程中，可以帮助研究人员节约相当一部分精力此外，它还可以使计算、证明、

2、建模等的结果得到更明白易懂的表现，有时，这比科学论证更有说服力 同时，数学的教学与学习过程也离不开绘图借助直观的图形，常可以使初学者更容易接受新知识如数学分析中有不少函数，其解析式着实让人望而生畏，即使对其性质作了详尽的分析，还是感到难明就里；但如果能看到它的图形，再配合理论分析，则问题可以迎刃而解又如在几何的学习中，会遇到大量的曲线与曲面，也离不开图形的配合 传统的手工作图，往往费力耗时，效果也不尽理想计算机恰恰弥补了这个不足，使你可以方便地指定各种视角、比例、明暗，从各个角度进行观察 本实验通过对函数的图形表示和几个曲面（线）图形的介绍，一方面展示它们的特点，另一方面，也将就Matlab软

3、件的作图功能作一个简单介绍 大家.word 格式,专业.专注 .将会看到，Matlab 的作图功能非常强大 二、相关函数（命令）及简介 1平面作图函数：plot，其基本调用形式：plot(x,y,s)以 x 作为横坐标，y 作为纵坐标s 是图形显示属性的设置选项例如：x=-pi:pi/10:pi;y=sin(x);plot(x,y,-rh,linewidth,2,markeredgecolor,b,markerfacecolor,g)图 1 在使用函数plot时，应当注意到当两个输入量同为向量时，向量x与y必须维数相同，而且必须同是行向量或者同是列向量 绘图时，可以制定标记的颜色和大小，也可以

4、用图形属性制定其他线条特征，这些属性包括：linewidth 指定线条的粗细 markeredgecolor 指定标记的边缘色 markerfacecolor 指定标记表面的颜色 markersize 指定标记的大小.word 格式,专业.专注 .若在一个坐标系中画几个函数，则 plot 的调用格式如下：plot(x1,y1,s1,x2,y2,s2,)2 空间曲线作图函数：plot3，它与 plot 相比，只是多了一个维数而已其调用格式如下：plot3(x,y,z,s)例如：x=0:pi/30:20*pi;y=sin(x);z=cos(x);plot3(x,y,z)得到三维螺旋线：图 2 3空

5、间曲面作图函数：（1）mesh 函数绘制彩色网格面图形调用格式：mesh(z)，mesh(x,y,z)和 mesh(x,y,z,c)其中，mesh(x,y,z,c)画出颜色由 c 指定的三维网格图若 x、y 均为向量，则length(x)=n，length(y)=m，m,n=size(z)（2）surf 在矩形区域内显示三维带阴影曲面图调用格式与 mesh 类似.word 格式,专业.专注 .（3）ezmesh 用符号函数作三维曲面网格图 调用格式：ezmesh(x,y,z)其中 x=x(s,t),y=y(s,t),z=z(s,t)画图区域默认为：-2*pi s 2*pi 且-2*pi t 7

6、)输出：b=3 4 6 7 三、实验内容 数学分析中，特别是积分部分，我们接触了不少有趣的函数，由于其中有的不是一一对应的，用上面的方法无法画出它们的图像，这时就只能用参数了 此外还有些图形只能用参数来画，比如空间曲线，在计算机上不接受“两个曲面的交线”这种表示，所以也只能用参数来实现 用参数方式作图的关键在于找出合适的参数表示，尤其是不能有奇点，最好也不要用到开方所以要找的参数最好是有几何意义的当然这也不可一概而论，需要多积累经验 1利用函数 plot 在一个坐标系中画以下几个函数图像，要求采用不同颜色、不同线形、不同的符号标记函数为：.word 格式,专业.专注 .sin(),cos(),

7、sin(2),(0,2)xtytztt 程序如下：t=0:pi/20:2*pi;x=sin(t);y=cos(t);z=sin(2*t);plot(t,x,-k*,t,y,-rs,t,z,:bo)图像如下：图 3 2绘制类似田螺线的一条三维螺线(方程自己设计)程序如下：t=0:.1:30;x=2*(cos(t)+t.*sin(t);y=2*(sin(t)-t.*cos(t);z=1.5*t;.word 格式,专业.专注 .plot3(x,y,-z)%取 z 主要是为了画图看起来更清楚 axis equal 图像如下：图 4 3利用函数2222sinxyzxy，绘制一个墨西哥帽子的图形 程序如下

8、：a,b=meshgrid(-8:.5:8);%先生成一个网格 c=sqrt(a.2+b.2)+eps;z=sin(c)./c;mesh(a,b,z)axis square 图像如下：.word 格式,专业.专注 .图 5 思考：这里的 eps 是什么？其作用是什么？4利用 surf 绘制马鞍面图形（函数为：2294xyz）程序如下：x,y=meshgrid(-25:1:25,-25:1:25);z=x.2/9-y.2/4;surf(x,y,z)title(马鞍面)grid off 图像如下：.word 格式,专业.专注 .图 6 5分别用 ezmesh 和 ezsurf 各绘制一个圆环面，尝

9、试将两个圆环面放在一个图形界面内，观察它们有什么不同之处 提示：圆环面的方程为：2 ,6 ,)(22222rRrzRyx，而圆环面的参数方程为：2,0 ,2,0 ,sinsin)cos(cos)cos(vuurzvurRyvurRx 程序参见附录 1 图像如下：图 7.word 格式,专业.专注 .6绘制黎曼函数图形，加深对黎曼函数的理解 说明：黎曼函数的定义为 1 (0,1)0101pppqxqqqyxx,当、为正整数，为既约分数,0，当，及无理点，程序参见附录 2 图像如下：图 8 四、自己动手 1 作出下图所示的三维图形：.word 格式,专业.专注 .图 9 提示：图形为圆环面和球面的

10、组合.2 作出下图所示的墨西哥帽子及其剪裁图形：图 10 3画出球面、椭球面、双叶双曲面、单叶双曲面.word 格式,专业.专注 .4若要求田螺线的一条轴截面的曲边是一条抛物线：0y时25xz试重新设计田螺线的参数方程，并画出该田螺线 5 作出下图所示的马鞍面（颜色为灰色，并有一个标题：“马鞍面”）：图 11 6 绘制图 8 所示的黎曼函数图形，要求分母的最大值n的数值由键盘输入（提示：使用 input 语句）五、附录 附录 1：（fulu1.m）程序如下：subplot(1,2,1)ezmesh(6+2*cos(u)*cos(v),(6+2*cos(u)*sin(v),2*sin(u),0,2*pi,0,2*pi)axis equal subplot(1,2,2)ezsurf(6+2*cos(u)*cos(v),(6+2*cos(u)*sin(v),2*sin(u),0,2*pi,0,2*pi).word 格式,专业.专注 .axis equal 附录 2：（fulu2.m）程序如下：n=100;x=;y=;k=1;for i=2:n for j=1:i-1 if gcd(i,j)=1%用函数 gcd(m,n)可求 m 和 n 的最大公约数 x(k)=j/i;y(k)=1/i;k=k+1;end end end plot(x,y,.b);axis(0,1,0,1)