Matlab数据与函数的作图a.pptVIP

第三章 数据与函数的作图 本章导读 作图功能是MATLAB的特长之一，MATLAB提供高端和低端作图命令，我们可以非常方便地将数据绘制成二维、三维的图形，而高端命令可以让我们随心所欲地画出各种函数图形。作图分为两大类，一类是数据作图，如给定数组x，y作数据的图形。一类是给一个函数直接作图，如做的图形。事实上图形正在演变成一种语言，它可以用更直观的形式展示问题的内涵。 在数学建模竞赛中，我们的论文广泛使用图文并茂的形式，可以更好地说明解决问题的思路。特别是在国际数学建模竞赛MCM、ICM中，由于中国学生的语言表达能力，特别是书面表达能力不如以英语为母语的美国、英国、澳大利亚的等国家，因此在论文中多插图形往往可以收到事半功倍的效果。 【例3.4.2】在人口预测的例子中，我们作出6子个图形。 3.6 特殊三维作图 有时在作立体图时我们需要挖去其中的某一部分或裁减某一部分，这时我们可以利用“非数”方法进行处理，从而获得挖去图形一块的效果。 3.6.1 图形的镂空 【例3.6.1】利用“非数”NaN，对图形的第三部分进行剪切处理。 clf; t=linspace(0,2*pi,100);r=1-exp(-t/2).*cos(4*t); %旋转母线 [X,Y,Z]=cylinder(r,60); % 产生旋转柱面数据 ii=find(X0Y0); % 确定x-y平面第三象限上的数据下标 Z(ii)=NaN; % 剪切 surf(X,Y,Z);colormap(spring),shading interp light(position,[-3,-1,3],style,local) %设置光源 material([0.5,0.4,0.3,10,0.3]) %设置表面反射 【例3.6.2】演示：如何利用“非数”NaN，对图形进行镂空处理。 P=peaks(30);P(18:20,9:15)=NaN; % 镂空 surfc(P);colormap(summer) light(position,[50,-10,5]),lighting flat material([0.9,0.9,0.6,15,0.4]) 3.6.2 裁切 【例3.6.3】将单位球x-y平面绝对值大于0.5的部分图形只置为0。 [X,Y,Z]=sphere(100) ii=find(abs(X).5|abs(Y).5); % 找x-y平面的带形部分 Z(ii)=zeros(size(ii)); %强制为0 surf(X,Y,Z) shading interp;colormap(hsv) % 对阴影进行平滑插值并置色 light(position,[0,-15,1]);lighting phong % 设置光线 material([0.8,0.8,0.5,10,0.5]) % 设置材质 3.7 统计图形的作图 在图形中另一类经常使用的图形是统计图形，统计模型在数学建模中所占比例不断加大。统计图形包括直方图、条形图、饼形图等等十分丰富。这里仅做简单介绍，我们将在统计工具箱中作详细介绍。 3.7.1 方差分析图 方差分析是对给定的样本推断样本是否来自一个总体，这仅介绍单因素方差分析的图表，它可以对方差分析表进行更直观的说明。结果是非常显著的不同。 【例3.7.1】Hogg 和 Ledolter (1987)对来自五个不同地方牛奶细菌的含量进行了研究，目的是研究不同地方牛奶的细菌含量是否有显著的差异。 load hogg % 调出hogg数据 hogg [p,tbl,stats] = anova1(hogg); % 产生方差分析表和图形 p % 显示显著性水平 hogg = 24 14 11 7 19 15 7 9 7 24 21 12 7 4 19 27 17 13 7 15 33 14 12 12 10 23 16 18 18 20 p = 1.1971e-004 3.7.2 正态分布检验图、经验分布图 检验一组随机样本是否为正态分布是统计分析中经常要使用的方法，用正态分布检验图可以直观地看出一组随机数是否为正态分布。 【例3.7.2】用正态分布随机数发生器产生20个随机数，再用卡方分布随机数发生器产生20个随机数，用正态分布检验图来判断那