第三章 作图2.pptVIP

x=3*pi*(-1:0.2:1);y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R; %引入eps避免0/0， C=abs(del2(Z)); %求“五点格式”差分，反映曲面变化。 meshz(X,Y,Z,C) %由曲面变化决定用色 hold on,scatter3(X(:),Y(:),Z(:),filled) hold off, colormap(hot) 例：矩阵作图指令plotmatrix有两种基本调用方式：（1）对于数据矩阵维的X和维的Y，调用格式plotmatrix(X,Y)将画出一个分割成个子散点图。其中第个子散点图是根据Y第列和X第列数据画出的。（2）对于数据矩阵维的X，调用格式plotmatrix(X)将画出分割成个子块的图。该图的对角块，画出的是X每列的数据的频数直方图；而其他子块是相应列构成的散点图。 该指令可用来观察数据矩阵（或同一矩阵列向量）间的统计关系。 randn(seed,1111),X=randn(100,2); Y=randn(100,2); subplot(1,3,1),plotmatrix(X) subplot(1,3,2),plotmatrix(X,X) subplot(1,3,3),plotmatrix(X,Y) 6. 14 ? 切片图和切片等位线图 例：利用slice和contourslice 表现MATLAB提供的无限大水体中水下射流速度数据flow 。flow是一组定义在三维空间上的函数数据。 在本例中，从图中的色标尺可知，深红色表示“正速度”（向图的左方），深蓝表示“负速度”（向图的右方）。 %以下指令用切面上的色彩表现射流速度 clf;[X,Y,Z,V]=flow; %取4个的射流数据矩阵，V是射流速度。 x1=min(min(min(X)));x2=max(max(max(X))); %取x坐标上下限 y1=min(min(min(Y)));y2=max(max(max(Y))); %取y坐标上下限 z1=min(min(min(Z)));z2=max(max(max(Z))); %取z坐标上下限 sx=linspace(x1+1.2,x2,5); %确定5个垂直x轴的切面坐标 sy=0; %在y=0处，取垂直y轴的切面 sz=0; %在z=0处，取垂直z轴的切面 slice(X,Y,Z,V,sx,sy,sz); %画切片图 view([-12,30]); shading interp;colormap jet; axis off;colorbar %以下指令用等位线表现射流速度 clf;v1=min(min(min(V)));v2=max(max(max(V))); %射流速度上下限 cv=linspace(v1,v2,15); %在射流上下限之间取15条等位线 contourslice(X,Y,Z,V,sx,sy,sz,cv); view([-12,30]) colormap jet;colorbar;box on 6.15 动态图形 例：卫星返回地球的运动轨线示意。 shg;R0=1; %以地球半径为一个单位 a=12*R0;b=9*R0;T0=2*pi; %T0是轨道周期 T=5*T0;dt=pi/100;t=[0:dt:T]; f=sqrt(a^2-b^2); %地球与另一焦点的距离 th=12.5*pi/180; %卫星轨道与x-y平面的倾角 E=exp(-t/20); %轨道收缩率 x=E.*(a*cos(t)-f);y=E.*(b*cos(th)*sin(t)); z=E.*(b*sin(th)*sin(t)); plot3(x,y,z,g) %画全程轨线 [X,Y,Z]=sphere(30);X=R0*X;Y=R0*Y;Z=R0