%%
%1.四个对数图
x = logspace(-1,1,100); 
y = x.^2;
subplot(2,2,1); 
plot(x,y); 
title('Plot');
subplot(2,2,2); 
semilogx(x,y); 
title('Semilogx');
subplot(2,2,3); 
semilogy(x,y); 
title('Semilogy');
subplot(2,2,4); 
loglog(x, y); 
title('Loglog');
%绘制X轴垂直网格线
set(gca,'XGrid','on');

%%
%2.plotyy绘制双纵坐标图
x = 0:0.01:20;
y1 = 200*exp(-0.05*x).*sin(x);
y2 = 0.8*exp(-0.5*x).*sin(10*x);
[AX,H1,H2] = plotyy(x,y1,x,y2);
set(get(AX(1),'Ylabel'),'String','Left Y-axis')
set(get(AX(2),'Ylabel'),'String','Right Y-axis')
title('Labeling plotyy');
set(H1,'LineStyle','--'); 
set(H2,'LineStyle',':');

%%
%3.柱状图（hist表示柱形数量）
y = randn(1,1000);
subplot(2,1,1);
hist(y,10);
title('Bins = 10');
subplot(2,1,2);
hist(y,50);
title('Bins = 50');

%%
%4.三种条形图（普通条形图、带颜色的对比条形图、对比3D条形图）
x = [1 2 5 4 8]; y = [x;1:5];
subplot(1,3,1); bar(x); title('A bargraph of vector x');
subplot(1,3,2); bar(y); title('A bargraph of vector y');
subplot(1,3,3); bar3(y); title('A 3D bargraph');

%%
%5.堆积条形图和水平条形图
x = [1 2 5 4 8]; 
y = [x;1:5];
subplot(1,2,1); 
bar(y,'stacked'); 
title('Stacked');
subplot(1,2,2); 
barh(y); 
title('Horizontal');

%%
%6.饼状图（普通饼状图，切割饼状图、3D切割饼状图）
a = [10 5 20 30];
subplot(1,3,1); pie(a);
subplot(1,3,2); pie(a, [0,0,0,1]);
subplot(1,3,3); pie3(a, [0,0,0,1]);

%%
%7.极坐标图
x = 1:100; theta = x/10; r = log10(x);
subplot(1,4,1); polar(theta,r);
theta = linspace(0, 2*pi); r = cos(4*theta); 
subplot(1,4,2); polar(theta, r);
theta = linspace(0, 2*pi, 6); r = ones(1,length(theta));
subplot(1,4,3); polar(theta,r);
theta = linspace(0, 2*pi); r = 1-sin(theta);
subplot(1,4,4); polar(theta , r);

%%
%8.梯形图
x = linspace(0, 4*pi, 40); y = sin(x);
subplot(1,2,1); stairs(y);
subplot(1,2,2); stem(y);

%%
%9.盒图
load carsmall
boxplot(MPG, Origin);

%%
%10.误差条形图
x=0:pi/10:pi; y=sin(x); 
e=std(y)*ones(size(x)); 
errorbar(x,y,e)

%%
%11.STOP标志（fill填充颜色）
t =(1:2:15)'*pi/8; x = sin(t); y = cos(t);
fill(x,y,'r'); axis square off;
text(0,0,'STOP','Color', 'w', 'FontSize', 80, ...
'FontWeight','bold', 'HorizontalAlignment', 'center');

%%
%12.三个国家2012年奥林匹克奖牌数量对比图
G = [46 38 29 24 13]; S = [29 27 17 26 8];
B = [29 23 19 32 7]; h = bar(1:5, [G' S' B']);
title('Medal count for top 5 countries in 2012 Olympics');
ylabel('Number of medals'); xlabel('Country');
legend('Gold', 'Silver', 'Bronze')

%%
%13.三维图像
[x, y] = meshgrid(-3:.2:3,-3:.2:3); 
z = x.^2 + x.*y + y.^2; surf( x, y, z); box on; 
set(gca,'FontSize', 16); zlabel('z'); 
xlim([-4 4]); xlabel('x'); ylim([-4 4]); ylabel('y'); 
%上图颜色分布图（将数据可视化）
imagesc(z); axis square; xlabel('x'); ylabel('y');

%%
%14.颜色条和方案
colorbar;%颜色普通分布图
colormap(hot);%热成像图
colormap(cool);%将图像凉爽化
colormap(gray);%灰阶图像

%%
%15.将图像绿色渐变
x = [1:10; 3:12; 5:14];
imagesc(x);
colorbar;
map = zeros(256,3);
map(:,2) = (0:255)/255;
colormap(map);

%%
%16.2D图像拉成3D图像
x=0:0.1:2*pi;
plot(x,sin(x));

%%
%17.plot3绘制3D图像
x=0:0.1:3*pi; z1=sin(x); z2=sin(2.*x); z3=sin(3.*x);
y1=zeros(size(x)); y3=ones(size(x)); y2=y3./2;
plot3(x,y1,z1,'r',x,y2,z2,'b',x,y3,z3,'g'); grid on; 
xlabel('x-axis'); ylabel('y-axis'); zlabel('z-axis');

%%
%18.更多3D线段图
%平行螺旋上升图
t = 0:pi/50:10*pi;
plot3(sin(t),cos(t),t)
grid on; axis square;
%螺旋上升至顶尖图
turns = 40*pi;
t = linspace(0,turns,4000);
x = cos(t).*(turns-t)./turns;
y = sin(t).*(turns-t)./turns;
z = t./turns;
plot3(x,y,z); grid on;

%%
%19.曲面绘图
x = -3.5:0.2:3.5; y = -3.5:0.2:3.5;
[X,Y] = meshgrid(x,y);
Z = X.*exp(-X.^2-Y.^2);
subplot(1,2,1); mesh(X,Y,Z);%未填充颜色
subplot(1,2,2); surf(X,Y,Z);%填充颜色

%%
%20.contour()画出3D转2D的轮廓
x = -3.5:0.2:3.5; 
y = -3.5:0.2:3.5;
[X,Y] = meshgrid(x,y); 
Z = X.*exp(-X.^2-Y.^2);
subplot(2,1,1); 
mesh(X,Y,Z); 
axis square;
subplot(2,1,2); 
contour(X,Y,Z); 
axis square;

%%
%21.等高线图
x = -3.5:0.2:3.5; y = -3.5:0.2:3.5;
[X,Y] = meshgrid(x,y); Z = X.*exp(-X.^2-Y.^2);
subplot(1,3,1); contour(Z,[-.45:.05:.45]); 
axis square;
subplot(1,3,2); [C,h] = contour(Z); 
clabel(C,h); 
axis square;
subplot(1,3,3); contourf(Z); 
axis square;

%%
%22.表面/网格和轮廓的组合
x = -3.5:0.2:3.5; y = -3.5:0.2:3.5;
[X,Y] = meshgrid(x,y); Z = X.*exp(-X.^2-Y.^2);
subplot(1,2,1); meshc(X,Y,Z);
subplot(1,2,2); surfc(X,Y,Z);

%%
%23.球形图
sphere(50); shading flat;
light('Position',[1 3 2]);
light('Position',[-3 -1 3]);
material shiny;
axis vis3d off;
set(gcf,'Color',[1 1 1]);
view(-45,20);

%%
%24.高亮球形图
[X, Y, Z] = sphere(64); h = surf(X, Y, Z); 
axis square vis3d off;
reds = zeros(256, 3); reds(:, 1) = (0:256.-1)/255;
colormap(reds); shading interp; lighting phong;
set(h, 'AmbientStrength', 0.75, 'DiffuseStrength', 0.5);
L1 = light('Position', [-1, -1, -1]);
set(L1, 'Position', [-1, -1, 1]);
set(L1, 'Color', 'g');

%%
%25.包含多边形的图形对象
v = [0 0 0; 1 0 0 ; 1 1 0; 0 1 0; 0.25 0.25 1; ...
0.75 0.25 1; 0.75 0.75 1; 0.25 0.75 1];
f = [1 2 3 4; 5 6 7 8; 1 2 6 5; 2 3 7 6; 3 4 8 7; 4 1 5 8];
subplot(1,2,1); patch('Vertices', v, 'Faces', f, ...
'FaceVertexCData', hsv(6), 'FaceColor', 'flat');
view(3); axis square tight; grid on;
subplot(1,2,2); 
patch('Vertices', v, 'Faces', f, ...
'FaceVertexCData', hsv(8), 'FaceColor', 'interp');
view(3); axis square tight; grid on;

%%
%26.3D山景图
load cape
X=conv2(ones(9,9)/81,cumsum(cumsum(randn(100,100)),2));
surf(X,'EdgeColor','none','EdgeLighting','Phong',...
'FaceColor','interp');
colormap(map); caxis([-10,300]);
grid off; axis off;