第三章 作图1ppt课件

1、3.1 作图命令的搜索作图命令的搜索在界面上点击在界面上点击help功能，然后进行功能，然后进行plot关键字的搜索见图关键字的搜索见图3.1.1在这里列出了几乎所有和在这里列出了几乎所有和plot有关的图形命令，单击任何你需要的有关的图形命令，单击任何你需要的命令，在右边的窗口中可以看到该命令的语法、使用方法和一些示命令，在右边的窗口中可以看到该命令的语法、使用方法和一些示例。例如打开例。例如打开plot，下面的例子初步展示，下面的例子初步展示MATLAB作离散图和连作离散图和连续图的例子。续图的例子。【例【例3.1.1】打开】打开plot，我们可以看到有关，我们可以看到有关plot的语法说

2、明，实例的语法说明，实例等。见图等。见图根据帮助功能我们就可以方便地编程，例如作函数的图形，这里用根据帮助功能我们就可以方便地编程，例如作函数的图形，这里用数据作图。数据作图。x = -pi:pi/10:pi; % 产生自变量数组产生自变量数组x（-pi，pi）y = tan(sin(x) - sin(tan(x); % 产生相应因变量点的数组产生相应因变量点的数组yplot(x,y,-rs,LineWidth,2, MarkerEdgeColor,k,. MarkerFaceColor,g,.MarkerSize,10) % 作数据作数据x,y的图形，并定义线的颜色、点的形状等。的图形，并定

3、义线的颜色、点的形状等。【例【例3.1.2】用图形表示连续调制波形。】用图形表示连续调制波形。t1=(0:11)/11*pi;% 产生产生12个自变量的点个自变量的点y1=sin(t1).*sin(9*t1); t2=(0:100)/100*pi; 产生产生101个自变量的点个自变量的点y2=sin(t2).*sin(9*t2);subplot(2,2,1),plot(t1,y1,r.),axis(0,pi,-1,1),title(子图子图 (1)subplot(2,2,2),plot(t2,y2,r.),axis(0,pi,-1,1),title(子图子图 (2)subplot(2,2,3)

4、,plot(t1,y1,t1,y1,r.)axis(0,pi,-1,1),title(子图子图 (3)subplot(2,2,4),plot(t2,y2),axis(0,pi,-1,1)title(子图子图 (4) 【例【例3.1.4】用曲面图表现函数。首先取】用曲面图表现函数。首先取x，y轴的点，利用轴的点，利用meshgrid命令产生命令产生x-y平面的网格点，在每个网格点上计算函数平面的网格点，在每个网格点上计算函数的值。的值。x=-4:4;y=x;X,Y=meshgrid(x,y); %生成生成 x-y 坐标坐标网格点矩阵网格点矩阵Z=X.2+Y.2;%计算格点上的函数值计算格点上的函

5、数值surf(X,Y,Z);hold on,colormap(hot)stem3(X,Y,Z,bo) %用来表现用来表现在格点上计算函数值在格点上计算函数值3.2 图形基本概念和图形设置命令图形基本概念和图形设置命令3.2.1 图形基本概念图形基本概念 图形事实上是我们使用计算机命令在计算机屏幕图形事实上是我们使用计算机命令在计算机屏幕root上绘制的各种对象上绘制的各种对象object的组合，这些对的组合，这些对象包括图形显示框及其子对象，如各种菜单、各种工具图象包括图形显示框及其子对象，如各种菜单、各种工具图标和正方形的空白绘图显示区，以及在这个区域中用标和正方形的空白绘图显示区，以及在这

6、个区域中用plot、surf等命令作出的图形对象、颜色对象、光线对等命令作出的图形对象、颜色对象、光线对象、注释对象等等。众多的对象组成的有机体形成了特色象、注释对象等等。众多的对象组成的有机体形成了特色鲜明的计算机图形。下面我们列出主要的对象及其注释，鲜明的计算机图形。下面我们列出主要的对象及其注释，见表见表3.2.1 对象是分层次的，在一个图形中有父对象、子对象、兄妹对象等。如对象是分层次的，在一个图形中有父对象、子对象、兄妹对象等。如Root是是Fiugre的父对象，的父对象，Image是是Figure的孙对象，而的孙对象，而Axes和和Uicontrol等是同辈对象。等是同辈对象。此图

7、形结构可以用构架数组来表示，如此图形结构可以用构架数组来表示，如MATLAB特有的句柄构架数组特有的句柄构架数组handle structure）。典型的层次见图）。典型的层次见图3.2.1 在一个图形中除了图形本身外还有各种修饰命令，例如图形的标题、坐标的名称、在一个图形中除了图形本身外还有各种修饰命令，例如图形的标题、坐标的名称、线条的颜色和粗细，在三维图形中还要考虑立体光线等等。合理地使用这些修饰命令将线条的颜色和粗细，在三维图形中还要考虑立体光线等等。合理地使用这些修饰命令将使我们的图形更好看，更具说服力。下例是由多个对象组成的图形。使我们的图形更好看，更具说服力。下例是由多个对象组成

8、的图形。3.2.2 图形各目标的设置命令图形各目标的设置命令set和信息传递命令和信息传递命令get当我们用高端作图命令如当我们用高端作图命令如plot、surf进行作图时，实际是进行作图时，实际是调用个图形的构架数组。包括图形框架、菜单、工具图标、调用个图形的构架数组。包括图形框架、菜单、工具图标、坐标、图形背景颜色以及图形等等各种层次对象坐标、图形背景颜色以及图形等等各种层次对象objects）。有关这些对象的设置都是系统内定的值，）。有关这些对象的设置都是系统内定的值，如果想改变某一个对象的属性的值，使用如果想改变某一个对象的属性的值，使用set命令来重新命令来重新设置。设置。 在一个图

9、形界面上对各对象之间进行信息传递，可以用在一个图形界面上对各对象之间进行信息传递，可以用get命令获取某窗口中的信息。命令获取某窗口中的信息。先介绍图形句柄的概念。我们来看几个术语：先介绍图形句柄的概念。我们来看几个术语：gcf：为当前图形句柄：为当前图形句柄current figure handle），即），即获得由作图命令产生的图获得由作图命令产生的图3.2.1的图形层次结构。的图形层次结构。 gca：为图形当前的坐标系统（：为图形当前的坐标系统（ the handle to the current axes）。）。设置图形参数命令设置图形参数命令set的语法和获取属性信息的命令的语法和获

10、取属性信息的命令get：set(H,PropertyName,PropertyValue,.) % 对句柄对句柄H各种属性重新设置各种属性重新设置set(H,a) % 和上面命令等价，而将参数全部放入构架数组和上面命令等价，而将参数全部放入构架数组aset(H,pn,pv.) % 和上面命令等价，而将参数全部放入元胞数组和上面命令等价，而将参数全部放入元胞数组pn和和pvset(H,pn,) % 对多个句柄进行同时设置对多个句柄进行同时设置a= set(h) % 对象对象h的属性和值放到构架数组的属性和值放到构架数组a中中a= get(0,Factory) % 将当前整个句柄所有属性的值放到构

11、架数组将当前整个句柄所有属性的值放到构架数组a中中a= get(0,FactoryObjectTypePropertyName) % 将当前句柄的某个属性和值放入将当前句柄的某个属性和值放入a中中a= set(h,Default) % 将对象将对象h的内定值放到构架数组中的内定值放到构架数组中a= set(h,DefaultObjectTypePropertyName) % 将某个属性放到构架数组中将某个属性放到构架数组中 = set(h,PropertyName) % 将某个属性值放到元胞数组中将某个属性值放到元胞数组中H：为句柄，如：为句柄，如gcf、gca等等。等等。PropertyNa

12、me：为该句柄下的属性名，如线型、点的颜色等。：为该句柄下的属性名，如线型、点的颜色等。PropertyValue：一个属性名的值，如线型实线或虚线等，点可：一个属性名的值，如线型实线或虚线等，点可以取红、蓝等等。以取红、蓝等等。a：为一个构架数组，用于存放参数：为一个构架数组，用于存放参数PropertyName,PropertyValue,.，当参数很多时用构架，当参数很多时用构架数组进行说明较为方便。数组进行说明较为方便。pn：为元胞数组，存放各种属性名：为元胞数组，存放各种属性名PropertyName1，PropertyName2，vp：为元胞数组，存放各种属性名相应的值：为元胞数组

13、，存放各种属性名相应的值PropertyValue1，PropertyValue2，【例【例3.2.1】作两幅图，并用】作两幅图，并用set命令重设一些属性，用命令重设一些属性，用get显示显示属性的改变。属性的改变。x=linspace(-pi,pi,10);y=sin(x)subplot(1,2,1),h=plot(x,y) % 用系统内定的参数作图用系统内定的参数作图a= get(h) % 显示图形的每种特性的参数显示图形的每种特性的参数subplot(1,2,2),h=plot(x,y) set(h,LineStyle,-,LineWidth,1,Marker,d)b=get(h) %

14、 显示修改过的特性的参数显示修改过的特性的参数属性参数为：属性参数为：b = BeingDeleted: off BusyAction: queue ButtonDownFcn: Children: 0 x1 double Clipping: onColor: 0 0 1 CreateFcn: DeleteFcn: EraseMode: normal HandleVisibility: on HitTest: on Interruptible: on LineStyle: - % 已改变已改变 LineWidth: 1 % 已改变已改变 Marker: diamond % 已改变已改变 Mar

15、kerEdgeColor: auto MarkerFaceColor: none MarkerSize: 6 Parent: 150.0005 Selected: off SelectionHighlight: on Tag: Type: line UIContextMenu: UserData: Visible: on XData: 1x10 double YData: 1x10 double ZData: 1x0 double【例【例3.2.2】对于例】对于例3.2.1可以用构架数组来改变特性的参数，以可以用构架数组来改变特性的参数，以下程序可以获得一样的结果。下程序可以获得一样的结果。x

16、=linspace(-pi,pi,10);y=sin(x)subplot(1,2,1),h=plot(x,y) % 用系统内定的参数作图用系统内定的参数作图a= get(h)subplot(1,2,2),h=plot(x,y) b=get(h)act.LineStyle=-; % 定义构架数组中线形的参数定义构架数组中线形的参数act.LineWidth=1; % 定义构架数组中线宽的参数定义构架数组中线宽的参数act.Marker=d; % 定义构架数组中点型的参数定义构架数组中点型的参数set(h,act)set和和get命令在作图，特别是以后介绍的界面设计中起着关键性命令在作图，特别是以

17、后介绍的界面设计中起着关键性的作用，在后面的部分还有介绍。进一步的学习可查阅的作用，在后面的部分还有介绍。进一步的学习可查阅help中的中的set和和get的关键字。的关键字。3.3 数据二维作图和图形编辑数据二维作图和图形编辑3.3.1 二维作图二维作图 表表3.2.1列出了作图的大概步骤列出了作图的大概步骤 【例【例3.3.1】以我国国民生产总值的数据为例，进行作图和拟合】以我国国民生产总值的数据为例，进行作图和拟合研究，首先作出数据的图形，然后我们进行拟合的各种实验，从研究，首先作出数据的图形，然后我们进行拟合的各种实验，从中找出满意的结果。中找出满意的结果。A=load(e:dataz

18、ggdp.txt) % 定义定义n行两列的时间序列数组行两列的时间序列数组x=A(:,1) % 将时间定义为数组将时间定义为数组xy=A(:,2) % 将将gdp国民生产总值定义为国民生产总值定义为yplot(x,y,r*) % 将数据点用红色，将数据点用红色，*形绘制图形形绘制图形结果我们可以看界面图结果我们可以看界面图3.3.1打开界面菜单打开界面菜单Tools，并点击，并点击Basic Fitting 我们就可以看到以下我们就可以看到以下界面。界面。见图见图3.3.2 在在Basic Fitting栏下我们可以选择各种拟合方法，选择线性栏下我们可以选择各种拟合方法，选择线性linear和

19、和7阶多项式进阶多项式进行拟合。对作残差图行拟合。对作残差图Plot residuals进行勾选，表示要做残差分析。我们可以明显地看进行勾选，表示要做残差分析。我们可以明显地看出出7阶多项式拟合明显好与线性拟合。在第二张图形上显示了两种拟合的残差，残差图更阶多项式拟合明显好与线性拟合。在第二张图形上显示了两种拟合的残差，残差图更证明了我们的判断。见图证明了我们的判断。见图3.3.3二维图形的编辑与完善二维图形的编辑与完善 另一方法是直接在图形窗口上对图形进行编辑。另一方法是直接在图形窗口上对图形进行编辑。一、用图形修饰命令对图形进行编辑一、用图形修饰命令对图形进行编辑用命令来完整地修饰图形，事

20、实上是编一套程序，其具体格用命令来完整地修饰图形，事实上是编一套程序，其具体格式为：式为：作图命令（作图命令（ 数据，子命令，子命令，数据，子命令，子命令，子命令，子命令 ）修饰命令修饰命令1修饰命令修饰命令2 修饰命令修饰命令n以下将二维作图的曲线线型选择和颜色选择命令罗列为表以下将二维作图的曲线线型选择和颜色选择命令罗列为表3.3.2 点的颜色、点边界的颜色、点的大小设置等见表点的颜色、点边界的颜色、点的大小设置等见表3.3.3 点型的选择见表点型的选择见表3.3.4 图形标注辅助修饰命令与图形子命令不同，它不在图形命令内，而是作为独立的命图形标注辅助修饰命令与图形子命令不同，它不在图形命

21、令内，而是作为独立的命令紧跟在图形命令后面。主要的图形辅助修饰命令见表令紧跟在图形命令后面。主要的图形辅助修饰命令见表3.3.5【例【例3.3.2】利用上面的修饰命令作一个完的整修饰图形，图形标】利用上面的修饰命令作一个完的整修饰图形，图形标题为题为“我的完善图我的完善图”。x = -pi:pi/10:pi;y = tan(sin(x) - sin(tan(x); %准备数据准备数据%作曲线图，曲线设为虚线，线宽为作曲线图，曲线设为虚线，线宽为2，红颜色。点的边界为黑色，红颜色。点的边界为黑色，点面色为绿色，点面色为绿色，%点的大小设为点的大小设为10。plot(x,y,-rs,LineWid

22、th,2,. MarkerEdgeColor,k,. MarkerFaceColor,g,. MarkerSize,10); legend(tan(sin(x)-sin(tan(x),1) % 在图内标注公式在图内标注公式xlabel(时间时间) % 坐标坐标x标为时间标为时间ylabel(函数函数) % 坐标坐标y标为函数标为函数 title(我的完善图我的完善图) % 图形标题注为图形标题注为“我的完善图我的完善图”clf;t=0:pi/50:2*pi;y=sin(t);plot(t,y);axis(0,2*pi,-1.2,1.2)text(pi/2,1,fontsize16leftarr

23、owitsin(t)fontname隶书隶书极大值极大值) 二、在窗口界面上对图形进行直接修饰二、在窗口界面上对图形进行直接修饰 利用命令进行修饰的大部分功能可以更方便、更直观地利用窗口利用命令进行修饰的大部分功能可以更方便、更直观地利用窗口编辑功能进行。免去了编程步骤。我们以下面的例子来说明。编辑功能进行。免去了编程步骤。我们以下面的例子来说明。【例【例3.3.4】用窗口编辑功能对图形进行修饰。首先仅将图形简单】用窗口编辑功能对图形进行修饰。首先仅将图形简单绘出。绘出。t=0:pi/50:2*pi;y=sin(t);plot(t,y)进入编辑状态进入编辑状态对图形加文字对图形加文字在图形中加

24、箭头在图形中加箭头 在编辑状态下我们选定曲线，即用鼠标单击曲线的任何部位，这时我们可以看见该在编辑状态下我们选定曲线，即用鼠标单击曲线的任何部位，这时我们可以看见该曲线对象变为当前状态。见图曲线对象变为当前状态。见图3.3.7双击该曲线，我们就可以看到设置曲线的界面，如图双击该曲线，我们就可以看到设置曲线的界面，如图3.3.8在这里我们可以设置曲线的类型如实线、各类虚线，线宽，颜色等，还可以设置点的类在这里我们可以设置曲线的类型如实线、各类虚线，线宽，颜色等，还可以设置点的类型、大小、边界颜色和点的颜色等等。现在我们来设置一下，同时加入一些注释，其效型、大小、边界颜色和点的颜色等等。现在我们来

25、设置一下，同时加入一些注释，其效果图见图果图见图3.3.93.4 数据三维图形绘制和各种修饰数据三维图形绘制和各种修饰 三维图形为立体图形，和二维图形有了本质的不同。三维图形为立体图形，和二维图形有了本质的不同。(1) 三维图形有了视角的概念，不同的视角对于同一物体会产生不同的视觉效果，因此三维图形有了视角的概念，不同的视角对于同一物体会产生不同的视觉效果，因此必须有刻画视角的必须有刻画视角的view修饰命令。修饰命令。(2) 光线的照射角度也会对物体的视觉效果产生影响，因此必须有刻画光线照射角度的光线的照射角度也会对物体的视觉效果产生影响，因此必须有刻画光线照射角度的命令。命令。(3) 颜色

26、的设置也会对视觉产生影响，因此要有刻画颜色的命令。颜色的设置也会对视觉产生影响，因此要有刻画颜色的命令。(4) 对于一个空间的物体，我们不但需要从外观上了解它的性状，同时需要对它的横截对于一个空间的物体，我们不但需要从外观上了解它的性状，同时需要对它的横截面进行了解。例如利用磁共振方法我们可以得到人的大脑的数据，我们可以利用解剖图面进行了解。例如利用磁共振方法我们可以得到人的大脑的数据，我们可以利用解剖图来检查脑部是否有病变，因此我们还需要有能够进行切片的命令来检查脑部是否有病变，因此我们还需要有能够进行切片的命令slice来描述这样的情况。来描述这样的情况。MATLAB为我们提供了丰富的立体

27、修饰命令，因此三维作图的内容多也较二维图形复杂为我们提供了丰富的立体修饰命令，因此三维作图的内容多也较二维图形复杂的多。的多。3.4.1 三维图形的绘制三维图形的绘制 MATLAB为我们提供了各种三维数据的曲面作图命令，它们可以作空间的曲线、空为我们提供了各种三维数据的曲面作图命令，它们可以作空间的曲线、空间的曲面、等高线图、各种方向上的投影图、空间的散点图等等。这些命令都是高端命间的曲面、等高线图、各种方向上的投影图、空间的散点图等等。这些命令都是高端命令，既可以根据数据矩阵作图，又可以根据函数进行作图。现将主要命令列表如下：令，既可以根据数据矩阵作图，又可以根据函数进行作图。现将主要命令列

28、表如下： 尽管三维图形的修饰命令比较多，然而其基本步骤和二维图形还是差不多的。表尽管三维图形的修饰命令比较多，然而其基本步骤和二维图形还是差不多的。表3.4.2列出了三维作图的基本步骤。列出了三维作图的基本步骤。 MATLAB中提供了很多实际例子的数据，这些例子都是广泛流行的教材中的实例，中提供了很多实际例子的数据，这些例子都是广泛流行的教材中的实例，为我们理解和学习提供了方便。如美国人口从为我们理解和学习提供了方便。如美国人口从1900年到年到1970的普查数据，的普查数据，peaks峰峰函数数据）、函数数据）、wind风向数据）、风向数据）、mri脑部磁共振数据等。下面我们来看一些实例。脑

29、部磁共振数据等。下面我们来看一些实例。 【例【例3.4.1】利用】利用plot3命令，对极坐标三维数据作空间曲线图。命令，对极坐标三维数据作空间曲线图。t = 0:pi/50:10*pi;plot3(t.*sin(t),t.*cos(t),t,LineWidth,2) % 作空间曲线图，线宽为作空间曲线图，线宽为2axis square; grid on % x,y,x坐标的比例是相等的坐标的比例是相等的peaks是是MATLAB提供的一个产生提供的一个产生peaks函数数据的命令，它可以根据用户的要求产函数数据的命令，它可以根据用户的要求产生生peaks数据点的稀密程度，从而产生精确度不同视

30、觉效果的图形。数据点的稀密程度，从而产生精确度不同视觉效果的图形。【例【例3.4.2】利用】利用MATLAB给定的数据矩阵给定的数据矩阵peaks作立体图。作立体图。X,Y,Z = peaks(50); % 用用peaks函数产生函数产生5050的的X-Y面上的点，并计算面上的点，并计算Zsurfc(X,Y,Z) % 用数据用数据X，Y，Z产生曲面图形和等高线图产生曲面图形和等高线图colormap hsv % 对图形着色对图形着色 axis(-3 3 -3 3 -10 5) % 给出作图的坐标范围给出作图的坐标范围注意，有关图形的其他修饰上面的程序没有写出，如图形的视角、光线的设置、曲面表注

31、意，有关图形的其他修饰上面的程序没有写出，如图形的视角、光线的设置、曲面表面的性质、曲面的材质等等，这时计算机对这些对象有一个系统内定值或者称为缺省。面的性质、曲面的材质等等，这时计算机对这些对象有一个系统内定值或者称为缺省。【例【例3.4.3】人体头部磁共振数据】人体头部磁共振数据mri为为MATLAB提供的案例数据，它是大脑内部结构提供的案例数据，它是大脑内部结构的影像数据。数据结构为四维的影像数据。数据结构为四维128128127数组，其中最后的维数是按数组，其中最后的维数是按z轴对轴对脑部切片的脑部切片的27幅横截面映像图。我们可以从各个角度、内部和外部来窥视这批数据。本幅横截面映像图

32、。我们可以从各个角度、内部和外部来窥视这批数据。本例只是观看该多维数据的外表，并显示对数据的平滑和光线的处理并作图。例只是观看该多维数据的外表，并显示对数据的平滑和光线的处理并作图。load mri % 读入脑磁共振数据并放如数组读入脑磁共振数据并放如数组D中中D = squeeze(D); % 抽去抽去D中第三维，数据变为立体数据中第三维，数据变为立体数据12812827）Ds = smooth3(D); % 对数据进行平滑处理对数据进行平滑处理hcap = patch(isocaps(D,5),. % 对脑上部用一片横截面盖住，对颜色进行插对脑上部用一片横截面盖住，对颜色进行插值平滑。值平

33、滑。 FaceColor,interp,. EdgeColor,none);hiso = patch(isosurface(Ds,5),. % 对脑周围数据建立外表面，并着色。对脑周围数据建立外表面，并着色。 FaceColor,1,.75,.65,. EdgeColor,none);view(45,30) % 建立视角建立视角axis tight % 设置坐标为紧凑参数值设置坐标为紧凑参数值daspect(1,1,.4) % 对坐标设置适当的比率对坐标设置适当的比率lightangle(45,30); % 定义光源的角度定义光源的角度set(gcf,Renderer,zbuffer); li

34、ghting phong % 对当前图形句柄进行设置对当前图形句柄进行设置isonormals(Ds,hiso) % 计算数据边缘的模计算数据边缘的模 set(hcap,AmbientStrength,.6) % 设置头顶部参数设置头顶部参数set(hiso,SpecularColorReflectance,0,SpecularExponent,50); % 设置设置表面的参数表面的参数3.4.2 图形颜色命令图形颜色命令colormap 现在我们开始学习图形修饰命令，三维图形修饰命令十分丰富，包括对图形上色、现在我们开始学习图形修饰命令，三维图形修饰命令十分丰富，包括对图形上色、设置图形的视

35、角、对图形进行平滑处理、设置图形目标的材质等等。设置图形的视角、对图形进行平滑处理、设置图形目标的材质等等。我们面对的是一个平面的屏幕，何以能显示一个立体的图形呢？例如我们在纸上画一个我们面对的是一个平面的屏幕，何以能显示一个立体的图形呢？例如我们在纸上画一个球。首先用圆规画一个圆，然后根据光学和透视原理在圆中画上深浅不同的色彩，这时球。首先用圆规画一个圆，然后根据光学和透视原理在圆中画上深浅不同的色彩，这时我们就可以在平面的纸上看出一个逼真的球形了。同样的道理，计算机显示一个立体物，我们就可以在平面的纸上看出一个逼真的球形了。同样的道理，计算机显示一个立体物，除了给定的数据点组成的图形框架外

36、，还对每一数据点映射一个颜色，或者说在作图时除了给定的数据点组成的图形框架外，还对每一数据点映射一个颜色，或者说在作图时每个数据矩阵对应一个颜色矩阵，我们就可以作出形象逼真的立体彩色图形了。颜色矩每个数据矩阵对应一个颜色矩阵，我们就可以作出形象逼真的立体彩色图形了。颜色矩阵是阵是3列列n行的，它对应行的，它对应n个数据点。个数据点。colormap就是对数据点加载某种颜色矩阵的函数，其示意图如下：就是对数据点加载某种颜色矩阵的函数，其示意图如下：nnnzyxzyxzyx222111321232221131211nnnccccccccc 数据矩阵数据矩阵DM 颜色矩阵颜色矩阵CMcolormap

37、(CM)这里，颜色矩阵的元素为界于这里，颜色矩阵的元素为界于0和和1之间的实数，因此颜色矩阵之间的实数，因此颜色矩阵CM的每一行都是由红、的每一行都是由红、绿、蓝调成的某种颜色。例如：表绿、蓝调成的某种颜色。例如：表3.4.3列出了部分颜色数据的调色结果。列出了部分颜色数据的调色结果。【例【例3.4.4】利用】利用colormap命令选择颜色趋势。读入液体流出的数据命令选择颜色趋势。读入液体流出的数据（1） 取灰色取灰色load flujetimage(X)colormap(gray)（2colormap选择参数为选择参数为jetload flujetimage(X)colormap(jet)

38、3.4.3 图形的视角选择图形的视角选择view 图形的视角可以很好地刻画从某个位置观察一个立体的图形，我们知道从不同的角图形的视角可以很好地刻画从某个位置观察一个立体的图形，我们知道从不同的角度看立体物体可以得到不同的图形感观，这个视角命令为度看立体物体可以得到不同的图形感观，这个视角命令为view。而视角还包括观察距离。而视角还包括观察距离的远近，这个命令为的远近，这个命令为zoom。从不同的视角观察物体，我们的视角事实上只要两个参数，。从不同的视角观察物体，我们的视角事实上只要两个参数，即水平旋转角和升降旋转角。因此视角的一般形式为：即水平旋转角和升降旋转角。因此视角的一般形式为：vie

39、w(az,el)view(az,el)view(x,y,z)这里这里az为为azimuth水平旋转角度，水平旋转角度，el为为elevation上下旋转角度，而上下旋转角度，而view(x,y,z)为在向量方向为在向量方向(x, y, z)上的视角。上的视角。我们来看以下示意图。我们来看以下示意图。【例【例3.4.5】当我们不使用】当我们不使用view命令时，系统有一个内定的视角，命令时，系统有一个内定的视角， 即即view(-37.5，30)，本例给出四个不同的角度的图形。，本例给出四个不同的角度的图形。X,Y = meshgrid(-2:.25:2); Z = X.*exp(-X.2 -Y

40、.2);subplot(2,2,1),surf(X,Y,Z),title(系统缺省视角系统缺省视角)subplot(2,2,2),surf(X,Y,Z),view(180,0),title(180,0)视角视角)subplot(2,2,3),surf(X,Y,Z),view(-37.5,-30),title(-37.5,-30)视角视角)subplot(2,2,4),surf(X,Y,Z),view(0,90),title(0,90)视角视角)3.4.4 光线对立体的影响光线对立体的影响 我们都知道光线对三维图形有影响，如一个球正对光线的部分发亮，而背我们都知道光线对三维图形有影响，如一个球正

41、对光线的部分发亮，而背光处则发暗。对光的反映还根据物体的材料不同而有不同，例如强光照在一光处则发暗。对光的反映还根据物体的材料不同而有不同，例如强光照在一个网球上没有非常强烈的反光，而如果强光照在一个抛光的金属球上则有非个网球上没有非常强烈的反光，而如果强光照在一个抛光的金属球上则有非常强烈的反光。所有这些现象在常强烈的反光。所有这些现象在MATLAB中都可以模拟，既有各种光线强中都可以模拟，既有各种光线强度、位置的命令，也有物体材质的刻画命令。度、位置的命令，也有物体材质的刻画命令。 创建光线对象命令创建光线对象命令light，light命令的语法为：命令的语法为：light(Propert

42、yName, PropertyValue, .)这里特性这里特性PropertyName 可以为：可以为：Color 光源的颜色光源的颜色 颜色的值可以用数组颜色的值可以用数组a,b,c，元素均为小于等于，元素均为小于等于1，大于等于，大于等于0的数，也可以用的数，也可以用y， r，分别表示黄，红等来表示。，分别表示黄，红等来表示。Style 光源的距离类型，即是无限远还是某距离的光源光源的距离类型，即是无限远还是某距离的光源 infinite | localPosition 光源的方位，某个角度光源的方位，某个角度 用用a,b,c表示光源的方位，系统内定为表示光源的方位，系统内定为1，0，1

43、创建材质对象的命令为创建材质对象的命令为material，其语法为：，其语法为： material shiny 抛光材质，即具有强烈的反光效应抛光材质，即具有强烈的反光效应material dull 具有吸光效应的材质，反光度差具有吸光效应的材质，反光度差 material metal 金属材质，具有反光效应金属材质，具有反光效应利用参数利用参数ka，kd，ks自定义反光效果的命令为：自定义反光效果的命令为：material(ka kd ks) material(ka kd ks n)material(ka kd ks n sc)【例【例3.4.6】作四个球，显示材质的不同，对光线的反映不同】

44、作四个球，显示材质的不同，对光线的反映不同cc=0,1,0; p=30subplot(2,2,1)； % 作第一幅子图，完全按系统内定的参数作第一幅子图，完全按系统内定的参数sphere(p); % 利用利用sphere函数产生圆球数据函数产生圆球数据title(material default )subplot(2,2,2) % 作第二幅子图作第二幅子图sphere(p);shading interp % 对阴影进行插值平滑处理对阴影进行插值平滑处理 set(gcf,Renderer,zbuffer) % 对图形句柄进行设置对图形句柄进行设置set(findobj(gca,type,surf

45、ace),. % 对图形对象进行表面处理，如光线对图形对象进行表面处理，如光线等等 FaceLighting,phong,. AmbientStrength,.3,DiffuseStrength,.8,. SpecularStrength,.9,SpecularExponent,25,. BackFaceLighting,unlit)light(Position,0 -2 1,color,cc) % 设置光源位置设置光源位置material dull % 设置图形对象的材质为吸收型材质设置图形对象的材质为吸收型材质title(material dull) subplot(2,2,3) sphe

46、re(p); shading interp set(gcf,Renderer,zbuffer)set(findobj(gca,type,surface),. FaceLighting,phong,. AmbientStrength,.3,DiffuseStrength,.8,. SpecularStrength,.9,SpecularExponent,25,. BackFaceLighting,unlit)light(Position,0 -2 1,color,cc)material shinytitle(material shiny )subplot(2,2,4) sphere(p); sh

47、ading interp set(gcf,Renderer,zbuffer)set(findobj(gca,type,surface),. FaceLighting,phong,. AmbientStrength,.3,DiffuseStrength,.8,. SpecularStrength,.9,SpecularExponent,25,. BackFaceLighting,unlit)light(Position,0 -2 1,color,cc)material metaltitle(material metal)colormap(jet)3.4.5 透视效果函数透视效果函数我们可以利用透

48、明函数对立体对象进行透明处理。透明处理的对象为物体表面、图像、我们可以利用透明函数对立体对象进行透明处理。透明处理的对象为物体表面、图像、边缘等。透明函数为边缘等。透明函数为alpha(p)，这里，这里p取取0，1之间的实数，当取之间的实数，当取0时图形对象完全时图形对象完全透明，而取透明，而取1时则完全不透明。时则完全不透明。【例【例3.4.7】对射流表面进行透明处理，射流数据为】对射流表面进行透明处理，射流数据为flow。计算结果见图。计算结果见图3.4.9x y z v = flow;subplot(1,2,1);p = patch(isosurface(x,y,z,v,-3);ison

49、ormals(x,y,z,v,p);set(p,facecolor,red,edgecolor,none);daspect(1 1 1);view(3); axis tight; grid on;camlight; lighting gouraud;title(系统内定值，不透明系统内定值，不透明)subplot(1,2,2);p = patch(isosurface(x,y,z,v,-3);isonormals(x,y,z,v,p);set(p,facecolor,red,edgecolor,none);daspect(1 1 1);view(3); axis tight; grid on;

50、camlight; lighting gouraud;alpha(.4);title(透明度为透明度为0.4)3.4.6 立体图形的切片处理立体图形的切片处理 在科学研究中我们常常将物体进行切片处理，观察物体内部的结构。典型的在科学研究中我们常常将物体进行切片处理，观察物体内部的结构。典型的有对人体的磁共振获得人体某部位的立体数据，然后进行截取某截面进行深有对人体的磁共振获得人体某部位的立体数据，然后进行截取某截面进行深入研究。或将人体组织进行冷冻，然后进行病理切片进行病理研究。入研究。或将人体组织进行冷冻，然后进行病理切片进行病理研究。2019年全国大学生数学建模竞赛就是有关切片问题的例子。

51、年全国大学生数学建模竞赛就是有关切片问题的例子。MATLAB提供了各种提供了各种命令可以对立体进行切片处理。切片命令的语法为：命令可以对立体进行切片处理。切片命令的语法为：slice(V,sx,sy,sz)slice(X,Y,Z,V,sx,sy,sz)slice(V,XI,YI,ZI)slice(X,Y,Z,V,XI,YI,ZI)slice(.,method)h = slice(.)slice命令可以分别按命令可以分别按x，y，z方向进行切片这里。方向进行切片这里。V：是立体数据：是立体数据sx，sy，sz：为数据点或向量，分别在三个轴上按向量给定的点作切面。：为数据点或向量，分别在三个轴上按

52、向量给定的点作切面。对对slice产生的切片作旋转处理，产生的切片作旋转处理，rotate命令可以对切片进行任意角度的旋转。命令可以对切片进行任意角度的旋转。rotate的语法为：的语法为：rotate(h,direction,alpha)rotate(.,origin)这里这里h：为由：为由slice建立的切片建立的切片direction：为二维或三维显示方向的数组：为二维或三维显示方向的数组alpha： 对给定的方向按右手准则旋转的角度对给定的方向按右手准则旋转的角度 【例【例5.4.8】 仍然是对射流数据进行处理，现在我们对它进行切片处理，并从不同的角仍然是对射流数据进行处理，现在我们对

53、它进行切片处理，并从不同的角度来看切片的图像。结果见图度来看切片的图像。结果见图5.4.10 x,y,z,v = flow;ax=-38.5,16;0,45 % 建立视角数据建立视角数据for m=1:2 % 建立切片大小建立切片大小xmin = min(x(:); ymin = min(y(:); zmin = min(z(:);xmax = max(x(:); ymax = max(y(:); zmax = max(z(:);% 同时作两幅不同方向的切片图同时作两幅不同方向的切片图subplot(1,2,m) hslice = surf(linspace(xmin,xmax,100),.

54、linspace(ymin,ymax,100),. zeros(100); % 建立切片建立切片rotate(hslice,-2,0,0,-45) %旋转切片旋转切片45度，形成新的切片。度，形成新的切片。% 获取新切片的数据获取新切片的数据xd = get(hslice,XData);yd = get(hslice,YData);zd = get(hslice,ZData);delete(hslice) % 将原切片删除将原切片删除% 对切片上色、插值加密等附加修饰对切片上色、插值加密等附加修饰h = slice(x,y,z,v,xd,yd,zd);set(h,FaceColor,inter

55、p, EdgeColor,none, DiffuseStrength,.8)daspect(1,1,1)axis tightbox offview(ax(m,:)camzoom(1.4)camproj perspectivelightangle(-45,45)colormap (jet(24)set(gcf,Renderer,zbuffer)end 图图 3.4.10【例【例3.4.9】本例为】本例为2019全国大学生数学建模竞赛全国大学生数学建模竞赛A题，题例见第十章。再现血管的原题，题例见第十章。再现血管的原状并切片，进行比较来验证恢复的正确性。图状并切片，进行比较来验证恢复的正确性。图3

56、.4.11为血管的再现图，图为血管的再现图，图3.4.12为切为切片比较图。片比较图。flops(0);tic load (f:cmcmtuResult) %读取数据，中轴线上读取数据，中轴线上100点的坐标，半径点的坐标，半径imgnumber=100R=sum(Result(:,1)/imgnumber;%Result(:,2:3)=Result(:,2:3)z=1:imgnumber;Poly_zx=polyfit(z,Result(z,2),6);Poly_zy=polyfit(z,Result(z,3),6);x=polyval(Poly_zx,z);y=polyval(Poly_z

57、y,z);% 用多项式拟合中轴方程并在上面取用多项式拟合中轴方程并在上面取100个点个点% 作三维图作三维图for k=1:imgnumber X,Y,Z=ellipsoid(x(k),y(k),z(k),R,R,R,10); h=surf(X,Y,Z); set(h,FaceColor,interp,. EdgeColor,none,. DiffuseStrength,.8) hold onendhold onalpha(0.2)hslice = surf(linspace(0,512,100),. linspace(0,512,100),. zeros(100)+30) % 建立切片建立切

58、片colormap (jet(24)delete(h)图图3.4.11是由给定的是由给定的100个等距离横截面图像，计算出血管的中轴线方程。然后在计算个等距离横截面图像，计算出血管的中轴线方程。然后在计算机上进行血管的原样恢复。为验证恢复的正确性，我们可以对血管进行切片并跟原截面机上进行血管的原样恢复。为验证恢复的正确性，我们可以对血管进行切片并跟原截面进行比较。图进行比较。图3.4.12为切片示意图，为说明问题对血管进行了透明处理。为切片示意图，为说明问题对血管进行了透明处理。事实上我们对事实上我们对100片截面进行了比较，这里仅取片截面进行了比较，这里仅取5、40、70、90四幅图像的比较

59、图，四幅图像的比较图，在图在图3.4.13中，左上图为第中，左上图为第5个截面与复原血管切片的图形，我们可以看出拟合的很好。个截面与复原血管切片的图形，我们可以看出拟合的很好。右上的小图为第右上的小图为第40幅截面与复原血管切片的图形，左下图为第幅截面与复原血管切片的图形，左下图为第70幅截面与复原血管切片幅截面与复原血管切片的图形，而右下图为第的图形，而右下图为第90幅截面与复原血管切片的图形。幅截面与复原血管切片的图形。我们可以看到第我们可以看到第90幅图像拟合的不好，这是因为复原的血管末端没有延伸出去。如果再幅图像拟合的不好，这是因为复原的血管末端没有延伸出去。如果再近一步进行处理，如对

60、血管进行外延，则结果会大大地改善。近一步进行处理，如对血管进行外延，则结果会大大地改善。3.5 多子图绘图多子图绘图MATLAB为我们提供了多图共显的功能，这为我们进行比较研为我们提供了多图共显的功能，这为我们进行比较研究提供了方便。前面我们已经大量使用了多图共显的图形。究提供了方便。前面我们已经大量使用了多图共显的图形。例如图例如图3.4.13，同时显示四幅截面拟合图形，比较好地说明，同时显示四幅截面拟合图形，比较好地说明了血管恢复的整体情况。了血管恢复的整体情况。多子图绘图的基本命令为：多子图绘图的基本命令为：subplotm, n, k） 表示在表示在mn个图形中的第个图形中的第k个子个子图。图