2011111MATLAB绘图及其矩阵运算课件.pptxVIP

MATLAB绘图及其矩阵运算主讲教师：叶春生华中科技大学材料学院csye@mail.hust.edu.cnTel:027ATLAB讲座MATLAB绘图及其矩阵运算一、数据和函数的可视化二、数值计算三、牛顿迭代法及其几何意义四、工程应用一、数据和函数的可视化 数据可视化的目的在于：通过图形，从一堆杂乱的离散数据中观察数据间的内在关系，感受由图形所传递的内在本质。 主要阐述：离散数据表示成图形的基本机理；曲线、曲面绘制的基本技法和指令；特殊图形的生成和使用示例；如何使用线型、色彩、数据点标记凸现不同数据的特征；如何利用着色、灯光照明、反射效果、材质体现和透明度处理渲染、烘托表现高维函数的性状；如何生成和运用标识，画龙点睛般地注释图形；如何显示和转换unit8、unit16、double三种不同数据类型所体现的变址、灰度和真彩图象，如何读写各种标准图象格式文件；如何通过图形窗的交互操作对图形进行修饰、调整；如何打印和输出图形文件。1 离散数据和离散函数的可视化【例1】用图形表示离散函数n=0:12;y=1./abs(n-6);plot(n,y,r*,MarkerSize,20)grid on Warning: Divide by zero.2 连续函数的可视化【例2】用图形表示连续调制波形t1=(0:11)/11*pi; y1=sin(t1).*sin(9*t1);t2=(0:100)/100*pi;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),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 plot的基本调用格式【例3】用图形表示连续调制波形及其包络线。 t=(0:pi/100:pi); y1=sin(t)*[1,-1]; y2=sin(t).*sin(9*t); t3=pi*(0:9)/9; y3=sin(t3).*sin(9*t3);plot(t,y1,r:,t,y2,b,t3,y3,bo)axis([0,pi,-1,1]) 【例4】采用模型 画一组椭圆。th = [0:pi/50:2*pi];a = [0.5:.5:4.5];X = cos(th)*a;Y = sin(th)*sqrt(25-a.^2);plot(X,Y),axis(equal), xlabel(x), ylabel(y)title(A set of Ellipses) 4 三维绘图的基本操作【例5】简单例题。t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,b-,x,y,z,bd),view([-82,58]),box on,legend(链,宝石) 5 网线图、曲面图基本指令格式【例6】用曲面图表现函数clf,x=-4:4;y=x;[X,Y]=meshgrid(x,y);Z=X.^2+Y.^2;surf(X,Y,Z);hold on,colormap(hot)stem3(X,Y,Z,bo) 6 图形的透视【例7】透视演示[X0,Y0,Z0]=sphere(30);X=2*X0;Y=2*Y0;Z=2*Z0;surf(X0,Y0,Z0);shading interphold on,mesh(X,Y,Z),colormap(hot),hold offhidden offaxis equal,axis off 7 图形的镂空【例8】演示：如何利用“非数”NaN，对图形进行剪切处理。t=linspace(0,2*pi,100); r=1-exp(-t/2).*cos(4*t);[X,Y,Z]=cylinder(r,60);ii=find(X0Y0);Z(ii)=NaN;surf(X,Y,Z);colormap(spring),shading interplight(position,[-3,-1,3],style,local)material([0.5,0.4,0.3,10,0.3]) 8 色彩控制【例9】演示：（A）用MATLAB预定义的两个色图矩阵，构成一个更大的色图阵。（B）把色轴的范围设置得比着色阵 C的数据范围小，使色图CM的两端色彩展宽，中间色彩