21Q!丛Q:!皇
TechnOIOgy
Conaultin暑Herald
Scienceand
学术论坛
应用MATLAB函数绘制=次曲面图①
袁玲
(株洲职业技术学院湖南株洲4'2001)
摘要:二次曲面方程是高等数学中重要的章节内容。其图像是主要研究对象,徒手画图有一定的难度,用matlab函数辅助绘图足最理想的工具,这篇文章用实例来介绍了6种matlab函数绘二次曲面图的方法,并进行了比较。关键词:MATLAB作图二次曲面中图分类号:TGl46文献标识码:A文章编号:1673—0534(2007)07(a)一0090—021引言
现代1=程上的许多问题都涉及到数学函数,数学函数在人们的日常生活中也是相当重要的。函数的具体表达形式是方程式,无论在《空间解析几何》中,还是《微积分》中,我们都曾经学过和用过,而二次曲面方程是数学中相当重要的方程,其图像是三维空间图像,是我们的主要研究对象,教学上老师们要经常画出它们的图像进行直观教学,画图是一个比较辣手的问题。针对画图问题,作者下面介绍了几种用MATLAB函数辅助绘曲面图的方法。2用MA'ILAB函数绘二次曲面图
绘制三维图形前,一般先用meshlzrid函数在平面(通常在XOY平面上)的一个矩形区间内打网格,标记网格点,再直接用相关函数完成作图【1。2),从下面的多种作图方法中町以看出meshgrid函数的重要性。2.1用surf或mesh函数绘图
SUrf函数绘制的是三维表面图,mesh函数绘制的是三维网格图¨J,当二次曲面方程是标准方程时,原方程式可化为z=rf’,)或V=t[x,z1或x=硝y,z)时,我们就用这两种函数完成绘图。
r2
v2
,2
>>subplot(233),surf(X,Y,z3)lholdon,
surf(x,Y,Z4)ltitle(‘surf函数作的双曲面’),
>>subplot(234)lmesh(x,Y,Z3)Iholdonlmesh(X,Y,24)Ititle(‘mesh函数作的双曲面’);
>>Z5=X. 謦}4. y.*Y/9I
曲面’),
>>Z5=X. X+4. Y. Y/9l
>>subplot(235)lplot3(X,Y,Z5)tgridOnI
title(‘plot3函数作的抛物面’)I
>>subplot(236),contour3(x,Y,Z5)Igridon}title(‘contour3函数作的抛物面’)I
结果见图2。
>>subplot(235)Isurf(x,Y,Z5) title(‘surf函数作的抛物面’)l
>>subplot(236)Imesh(x,Y,Z5)Ititle(‘mesh函数作的抛物面’)I
结果见图1。
图2
图l
例1分别绘曲面寻+寺+素21,
等+等一吾引和等+等=;在区域一2
《x<2,一3《1,≤3,一4<:<4内的图像。
2.2用plot3或contour3函数绘图
PlOt3函数绘制的是二维直角坐标曲线图,contour3函数绘制的是三维等高线图‘11,对于例l我们可以通过这两种函数绘曲线来组成曲面图。
>>x=一2:0.1:2lIt=--3:0.1:3l【Z,明
2.3用ezsurf或ezmssh函数绘图
ezsurf函数和ezmesh函数主要针对参数方程的三维作图函数,它们是专业作图函数‘34I,ezsurf函数绘制的是三维表面图,ezmesh函数绘制的是三维网格图,当二次曲面方程可化为参数方程时,就可以用这两种函数完成绘图。下面将例l中的方程化为参数方程,再用ezsurf函数和ezmesh函数绘曲面图。
z一斟愕+等以
峰掣z2±4J卜等一等、
广——:——彳
=meshgrid(x,’’),
>>z1=4.*sqrt(1-X. .L/4-Y.*y/9)I
Z2=-4.*sqrt(1一X. X./4-Y. y/9)l
椭球方程了X2+了y2+百Z2=l化为参数方
h=2c∞口cos夕
>>subplot(231),plot3(X,Y,ZI),hold
Onlplot3(X,Y,z2)
>>grid
onItitle(‘plot3函数作的椭球
】、z=J2+等,翮matlab命+iKq:
>>x=一2:0.1:2Iv=-3:0.1:3“X,明
=meshgrid(x,v),
>>zl=4.,sqrt(1一x. j,./4-Y.*r/9)lz2=-4.,sqrt(1一x. X./4一Y. Y/9)l
>>subplot(231)Isuff(X,Y,Z1)Ihold
onl
面’),
>>subplot(232)lcontour3(五r'ZI)Iholdonlcontour3(j,'y,z2)
>>gridOilIUUe(‘contour3函数作的椭球面’)I
>>23=4.*sqrt(x. Z./4+y.*y/9-1)l
-Iy=z3:si4n咖Ole口os,,其中(。《口<2 pi广
pi/2<#<pi/2)s双曲方程等+等一言.1
卜=2tcosat
化为参数方程1y=3ts广ina,其中(o《口
【z=|=44t‘一1
<2 Pi。,≥1或t‘一1)I抛物面方程
,
,
Z4=-4. sqrt(x. X./4+Y. ∥9一1)l
>>subplot(233)!plot3(兄Y,z3),hold
onl
surt(x,Y。z笛ltitle(‘surf函数作的椭球面’)。
>>subplot(232);mesh(x,Y,z1)lholdOrll
mesh(x,Y,22)。title(‘mesh函数作的椭球面’),
plot3(x,Y,z4)
>>gridOilItitle(‘plot3函数作的双曲面’)l
>>subplot(234) contour3(兄Y,Z3)Ihold
等+等2寺化为参数方程1yi嚣口 其
中(0<口<2*pi,一m<t<+*)。再按例I的要求绘曲面图:
f垆2f咖口
>>z3=4.*sqrt(x. x./4+Y. ∥9-1)I
Z4=m4.*sqrt(X. x./4+Y.}Y/9一1)l
onlcontour3(X,Y。24)
>>grid
On。title(‘contour3函数作的双
①2007年株洲职业技术学院院级立项课题:编号zZYKY0708,名称《应用MATLAB辅助理工专业课程教学的研究》I2006年株洲职业技术学院院级课题资助项目:编号zZYKY0603,名称《高职院校高等数学精品课程开发与建设研究》。
90科技咨询导报Science
and
TechnologyConsulting
Herald
万方数据