jetzt

(2).球面图

球面图绘制由函数sphere来实现

[X,Y,Z]=sphere(N)             此函数生成3个(N+1)*(N+1)的矩阵,利用函数        surf(X,Y,Z) 可产生单位球面.

[X,Y,Z]=sphere         此形式使用了默认值N=20.

Sphere(N)             只是绘制了球面图而不返回任何值.

例.绘制地球表面的气温分布示意图.

[a,b,c]=sphere(40);

t=abs(c);

surf(a,b,c,t);

axis('equal')   %此两句控制坐标轴的大小相同.

axis('square')

colormap('hot')

jetzt

例.旋转柱面图.

r=abs(exp(-0.25*t).*sin(t));

t=0:pi/12:3*pi;

r=abs(exp(-0.25*t).*sin(t));

[X,Y,Z]=cylinder(r,30);

mesh(X,Y,Z)

colormap([1 0 0])

jetzt

3.       三维旋转体的绘制

为了一些专业用户可以更方便地绘制出三维旋转体,MATLAB专门提供了2个函数:柱面函数cylinder和球面函数sphere

(1)   柱面图

柱面图绘制由函数cylinder实现.

[X,Y,Z]=cylinder(R,N)  此函数以母线向量R生成单位柱面.母线向量R是在单位高度里等分刻度上定义的半径向量.N为旋转圆周上的分格线的条数.可以用surf(X,Y,Z)来表示此柱面.

[X,Y,Z]=cylinder(R)或[X,Y,Z]=cylinder此形式为默认N=20且R=[1 1]


例.柱面函数演示举例

x=0:pi/20:pi*3;

r=5+cos(x);

[a,b,c]=cylinder(r,30);

mesh(a,b,c)

jetzt

2.       裁剪处理

利用不定数NaN的特点,可以对网图进行裁剪处理

例.图形裁剪处理

P=peaks(30);

subplot(2,1,1);

mesh(P);

title('裁剪前的网图')

subplot(2,1,2);

P(20:23,9:15)=NaN*ones(4,7);        %剪孔

meshz(P)                         %垂帘网线图

title('裁剪后的网图')

colormap([0 0 1])                   %蓝色网线

注意裁剪时矩阵的对应关系,即大小一定要相同.

jetzt

1. 消隐处理

例.比较网图消隐前后的图形

z=peaks(50);

subplot(2,1,1);

mesh(z);

title('消隐前的网图')

hidden off

subplot(2,1,2)

mesh(z);

title('消隐后的网图')

hidden on

colormap([0 0 1])

jetzt

我们可用subplot来同时画出数个小图形於同一个视窗之中：   

subplot(2,2,1); plot(x, sin(x));   

subplot(2,2,2); plot(x, cos(x));   

subplot(2,2,3); plot(x, sinh(x));   

subplot(2,2,4); plot(x, cosh(x));

jetzt

此外，MATLAB也可对图形加上各种注解与处理：   

xlabel('Input Value'); % x轴注解   

ylabel('Function Value'); % y轴注解   

title('Two Trigonometric Functions'); % 图形标题   

legend('y = sin(x)','y = cos(x)'); % 图形注解   

grid on; % 显示格线

jetzt

小整理：plot绘图函数的叁数 字元 颜色字元 图线型态y 黄色. 点k 黑色o 圆w 白色x  xb 蓝色+ +g 绿色* *r 红色- 实线c 亮青色: 点线m 锰紫色-. 点虚线-- 虚线  

图形完成後，我们可用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围：   

axis([0, 6, -1.2, 1.2]);

jetzt

若要同时改变颜色及图线型态（Line style），也是在座标对後面加上相关字串即可：  

plot(x, sin(x), 'co', x, cos(x), 'g*');

jetzt

若要改变颜色，在座标对後面加上相关字串即可：   

plot(x, sin(x), 'c', x, cos(x), 'g');