发布时间 : 星期五 文章[笔记]MATLAB中的图形(3)更新完毕开始阅读8f127ab485868762caaedd3383c4bb4cf7ecb7f9
图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.