matlab四维图形,MATLAB二维、三维、四维绘图

二维图形绘制直角坐标极坐标统计图面域图填色图羽毛图、射线图Voronoi图彩带图二维半图三维图形绘制线状图数组图网格图面状图统计图瀑布图散点图函数图四维图形绘制切片图颜色权重图Reference

二维作图一般使用plot, polar函数，辅以axis, xlabel, ylabel, legend, legend, title等添加坐标轴标题、图表标题等信息。也可以打开数据表格，选中需要绘图的数据，使用Matlab提供的plot工具条，绘制图形后，在figure的工具栏中点选添加坐标轴标题、图表标题等信息。

1%二维绘图示例代码：直角坐标 2%准备数据 3x = 0:0.01:3*pi; %pi特殊变量，圆周率 4y1 = sin(x); 5y2 = sin(2*x); 6y3 = sin(3*x); 7%设置当前绘图区域 8figure; 9%绘图10plot(x,y1,x,y2,x,y3);11%设置坐标轴和网格属性12axis([0 8 -2 2]); %axis，坐标轴13grid on;14%标注图形15xlabel('x');16ylabel('y');17title('演示绘图基本步骤');18legend('sin(x)','sin(2x)','sin(3x)')极坐标1%二维绘图示例代码：极坐标2t = 0:pi/50:2*pi;3r = sin(t).*cos(t);4polar(t,r,'-*')统计图 1%二维绘图示例代码：统计图 2figure; 3subplot(221) 4x = -2.9:0.2:2.9; %条形图 5bar(x,exp(-x.^2)) 6subplot(222) 7x = 0:0.1:4; %针状图 8y = (x.^0.8).*exp(-x) 9stem(x,y)10subplot(223)11x = 0:0.25:10; %阶状图12stairs(x,sin(2*x)+sin(x))13subplot(224)14x = [43 78 88 43 21];%饼图15pie(x)面域图1%二维绘图示例代码：面域图2clf;x = -2:2; %注意：自变量要单调变化,clf为清除图形命令3Y = [3,5,2,4,1;3,4,5,2,1;5,4,3,2,5]; %各因素的相对贡献份额4%Cum_Sum = cumsum(Y); %各曲线在图上的绝对坐标5area(x',Y',0);6legend('因素A','因素B','因素C'),grid on,colormap(spring)填色图1%二维绘图示例代码：填色图2X = [0.5 0.5 0.5 0.5;0.5 0.5 0.5 0.5;0 1 1 0];3Y = [0.5 0.5 0.5 0.5;0.5 0.5 0.5 0.5;0 0 1 1];4Z = [1 1 1 1;0 0 0 0;0 0 0 0];5C = [1 0 0 1;0 1 0 1;0 0 1 0];6fill3(X,Y,Z,C),view([-10 55]),colormap cool7xlabel('x'),ylabel('y'),box on;grid on羽毛图、射线图1%二维绘图示例代码：羽毛图、射线图2t = -pi/2:pi/12:pi/2;3r = ones(size(t)); %单位半径4[x,y] = pol2cart(t,r); %极坐标转为直角坐标5subplot(1,2,1),compass(x,y),title('射线图')6subplot(1,2,2),feather(x,y),title('羽毛图')Voronoi图1%二维绘图示例代码：Voronoi图2clf;rand('state',111)3n = 30;A = rand(n,1)-0.5;B = rand(n,1)-0.5; %产生30个随机点4T = delaunay(A,B); %求相邻三点组5T = [T T(:,1)]; %为使三点剖分三角形封闭而采取的措施6voronoi(A,B) %画Voronoi图7hold on;axis square;8fill(A(T(10,:)),B(T(10,:)),'y'); %画一个剖分三角形9voronoi(A,B) %重画Voronoi图，避免线被覆盖彩带图 1%二维绘图示例代码：彩带图 2clear,clf 3zeta2 = [0.1 0.2 0.3 0.4 0.5 0.6 0.8 1.0]; 4n = length(zeta2); 5for k=1:n 6 Num{k,1}=1; 7 Den{k,1}=[1 2*zeta2(k) 1]; 8end 9S = tf(Num,Den); %产生单输入多输出系统10t = (0:0.4:30)'; %时间采样点'11[Y,x] = step(S,t); %单输入多输出系统的响应12tt = t*ones(size(zeta2)); %为画彩带图，生成为函数值Y维数相同的时间矩阵13ribbon(tt,Y,0.4) %画彩带图14%至此彩带图已经生成，以下命令都是为了使图形效果更好、标示更清楚而用。15view([150,50]),shading interp,colormap(jet) %设置视角、明暗和色图16light,lighting phong,box on %设置光源、照射模式和坐标框二维半图伪彩图pcolor；

等位线命令contour、contourf；

等位线标高命令clabel。

1%二维绘图示例代码：二维半图 2clf;clear;[X,Y,Z]=peaks(40); %获得peaks图形数据 3n = 4; %等高线分级数 4subplot(1,2,1),pcolor(X,Y,Z) %伪彩图 5colormap jet,shading interp 6hold on,C = contour(X,Y,Z,n,'k:'); %用黑虚线画等位线，并给出标识数据 7clabel(C) %随机标识法 8zmax = max(max(Z)); 9zmin = min(min(Z));10caxis([zmin,zmax]);11colorbar %画垂直色标尺12hold off,subplot(1,2,2)13[C,h,CF] = contourf(X,Y,Z,n,'k:');14clabel(C,h) %沿线标识法三维图形绘制线状图plot3：绘制三维曲线图形；

stem3：绘制三维枝干图；

axis：设置坐标属性；

grid on：打开坐标网络；

view(方位角，仰俯角)：设置视角。

1%三维绘图示例代码：线状图2t = 0:pi/50:10*pi;3plot3(sin(t),cos(t),t) %画三维线状图4axis square; %使各坐标轴的长度相等5grid on数组图根据三个互相对应的数组绘制曲面图。

1%三维绘图示例代码：数组图 2%Author: Jing-Xuan Yang 3 4e = syntheticopioidsdatatotal.YYYY; 5c = syntheticopioidsdatatotal.longitude; 6d = syntheticopioidsdatatotal.DrugReports; 7 8[md,nd] = size(d); 910for i = 1:md11 if d(i)>80012 d(i) = 100;13 end14end1516d = d + 1000;17[X,Y]=meshgrid(linspace(min(e),max(e),100),linspace(min(c),max(c),100)) ;18Z=griddata(e,c,d,X,Y);19figure(1);20surf(X,Y,Z,'FaceColor','interp','EdgeColor','none','FaceLighting','phong');21axis tight;22xlabel('Latitude');23ylabel('Longitude');24zlabel('Alarm Threshold Level');25colorbar;26title('Alarm Threshold Level of Synthetic Opioids');网格图mesh：绘制三维网格图。

meshc：绘制带有基本等高线的网格图。

meshz：绘制带有基准平面的网格图。

1%三维绘图示例代码：网格图 2x = -8:.25:8; 3y = x; 4[X,Y] = meshgrid(x,y); %生成网格点坐标 5R = sqrt(X.^2+Y.^2)+eps;%定义R 6Z = sin(R)./R; %生成函数Z 7subplot(1,3,1); 8mesh(Z) %画网格图 9title('mesh(Z)');10subplot(1,3,2);11meshc(Z) %画带有基本等高线的网格图12title('meshc(Z)');13subplot(1,3,3);14meshz(Z) %画带有基准平面的网格图15title('meshz(Z)');面状图surf：绘制面状图。

surfl：(是L不是1)绘制设定光源方向的面状图。

shading interp和shading flat：把曲面上的小格平滑掉，使光滑。

shading faceted：是默认状态，它使曲面上有小格。

colormap：设定图形的颜色。

1%三维绘图示例代码：面状图 2z = peaks(25); %定义一个25×25的高斯分布矩阵 3subplot(1,2,1) 4surf(z); %画面状图 5shading interp; 6colormap(hsv); 7title('surf(z)'); 8subplot(1,2,2) 9surfl(z)10shading interp;11colormap(hsv);12title('surfl(z)');统计图bar3：绘制垂直方向直方图。

bar3h：绘制水平方向直方图。

stem3：绘制三维枝干图。

contour3：绘制等高线图。

1%三维绘图示例代码：直方图 2y1 = [5 2 1;8 7 3;8 5 6;5 1 5;4 3 2]; 3subplot(2,2,1) 4bar3(y1) %创建三维垂直方向直方图 5subplot(2,2,2) 6bar3h(y1) %创建三维水平方向直方图 7 8%三维绘图示例代码：枝干图 9t = (0:127)/128*2*pi;10×2 = t.*sin(t);11y2 = t.*cos(t);12subplot(2,2,3)13stem3(x2,y2,t,'fill') %画枝干图，并指定点为实心14view([65 30]) %设置视角1516%三维绘图示例代码：等高线图17[x3,y3,z3] = peaks; %定义一高斯分布矩阵18subplot(2,2,4)19contour3(x3,y3,z3,30) %生成30条等高线瀑布图1%三维绘图示例代码：瀑布图2%绘制瀑布图3Z = peaks(65);4waterfall(Z) %画瀑布图5shading faceted,colormap(copper); %使曲面有小格，铜色散点图 1%三维绘图示例代码：散点图 2x = 3*pi*(-1:0.2:1); 3y = x; 4[X,Y] = meshgrid(x,y); 5R = sqrt(X.^2+Y.^2)+eps; 6Z = sin(R)./R; %引入eps避免0/0 7C = abs(del2(Z)); %求“五点格式”差分，反映曲面变化 8meshz(X,Y,Z,C) %由曲面变化决定用色 9hold on,scatter3(X(:),Y(:),Z(:),'filled')10hold off,colormap(hot);函数图 1%三维绘图示例代码：抛物面 2×1 = -4:.25:4; 3y1 = x1; 4[X1,Y1] = meshgrid(x1,y1); %生成网格点坐标 5Z1 = X1.^2+Y1.^2; 6subplot(1,2,1) 7mesh(X1,Y1,Z1) 8 9%三维绘图示例代码：球面10×2 = -1:.2:1;11y2 = x2;12[X2,Y2] = meshgrid(x2,y2);13Z2 = sqrt(4-X2.^2-Y2.^2)+eps;14subplot(1,2,2)15mesh(X2,Y2,Z2)四维图形绘制切片图 1%四维绘图示例代码：切片图 2clf;[X,Y,Z,V]=flow; %取4个方向的射流数据矩阵，V是射流速度 3×1 = min(min(min(X)));%取x坐标上下限 4×2 = max(max(max(X))); 5y1 = min(min(min(Y)));%取y坐标上下限 6y2 = max(max(max(Y))); 7z1 = min(min(min(Z)));%取z坐标上下限 8z2 = max(max(max(Z))); 9sx = linspace(x1+1.2,×2,5); %确定5个垂直于x轴的切面坐标10sy = 0,sz = 0; %在y=0,z=0处，取垂直切面11slice(X,Y,Z,V,sx,sy,sz); %画切面图12view([-12 30]),shading interp;colormap jet;axis off;colorbar;颜色权重图 1%四维绘图示例代码：颜色权重图 2%Author: Jing-Xuan Yang 3 4b = Part1Heroindatacombined.latitude; 5c = Part1Heroindatacombined.longitude; 6a = Part1Heroindatacombined.YYYY; 7d = Part1Heroindatacombined.DrugReports; 8 9[X,Y]=meshgrid(linspace(min(b),max(b),100),linspace(min(c),max(c),100)) ;10Z=griddata(b,c,a,X,Y);11V = griddata(b,c,d,X,Y);12figure(1);13surf(X,Y,Z,V,'FaceColor','interp','EdgeColor','none','FaceLighting','phong');14axis tight;15xlabel('Latitude');16ylabel('Longitude');17zlabel('DrugReports');18colorbar;19title('Distribution of Herion in Five States');ReferenceCSDN: