[笔记]MATLAB中的图形(3)

图4

例5 分别用ezmesh和ezsurf各绘制一个圆环面，尝试将两个圆环面放在一个图形界面内，观察他们有什么不同之处。

提示：圆环面的方程为

，而圆环面的参数方程

为：

程序如下： % fun1.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]

)) axis equal 图形如图5：

图5

例6 绘制黎曼函数图形。

提示：黎曼函数的定义为 y=

1/q (当p,q为正整数，p/q为既约分数，x=p/q属于(0,1)) 0 (当x=0,1及无理点，x属于[0,1]) 程序如下： % fun2.m n=100; x=[]; y=[]; k=1; for i=2:n for j=1:i-1

if gcd(i,j)==1 %求最大公约数 x(k)=j/i; y(k)=1/i; k=k+1; end end end

plot(x,y,'.b') axis([0,1,0,1]) 图形如图6：

图6

参考文献:

周品，等.MATLAB数学建模与仿真.北京：国防工业出版社，2009.