Matlab系列之二维图形（下）

x1=0:0.1:4;
x2=0:0.2:4;
y1=x1.*sin(x1);
y2=x2.*cos(x2);
plotyy(x1,y1,x2,y2,@stem,'plot');
text(2,-1,'\ity=x*sinx');
text(3,1,'\ity=x*cosx');

subplot(221);fplot(@sin,[0,4*pi],'-*');
subplot(222);fplot(@(x)sin(1/x),[0.01,0.1],'-*');
subplot(223);fplot(@(x)abs((exp(-j*x*(0:9)))*ones(10,1)),[0,2*pi],'-*');
subplot(224);fplot(@(x)[sin(x),cos(x),tan(x)],[-2*pi,2*pi],'-*');

ezplot(f,[min,max])%设置x的自变量范围[min,max]，y方向默认
ezplot(f,[xmin,xmax,ymin,ymax])%x方向的变量范围[xmin,xmax]，y方向的变量范围[ymin,ymax]
ezplot(x,y)%(-2*pi到2*pi)默认范围绘制x和y
ezplot(x,y,[min,max])%[min,max]范围内绘制x,y

ezplot('x+y=3',[-5 5 -1 3])
axis square

semilogx(...)%在x轴上采用常用对数进行标定，使用格式和plot相同，以下两个也一样
semilogy(...)%在y轴上采用常用对数进行标定
loglog(...)%x和y都采用常用对数进行标定

x=0.001:0.001:20;
y=5+log(x)+x;
semilogx(x,y,'b');
hold on
semilogx(x,x+3,'g');

%y=x^3函数图
x=[1:100];
subplot(311);
plot(x,x.^3);
title('plot y=x^3')

%y=x^3对数坐标图
subplot(312)
loglog(x,x.^3);
title('loglog logy=logx^3')

%y=x^3半对数坐标图
subplot(313)
semilogy(x,x.^3);
title('semilogy logy=x^3')

polar(theta,rho)%theta：夹角 rho：极径
polar(theta,rho,s)%s指定极坐标图形的样式和颜色

theta=0:0.01:2*pi;
rho=sin(2*theta).*cos(2*theta);
polar(theta,rho,'.r');

调用：bar(x,y,width,param);
说明：x,y分别为横纵坐标向量，x在默认状态下取值为y的向量长度，假设其值分别为m和n，则构成m x n的矩阵，所得的图形即：m组柱状图，每组有n条；width为竖条的宽度，默认是0.8，若值超过1，则会竖条之间会发生重叠；param为参数，有分组式（grouped）和累加式（stacked）,默认是分组式；
备注：bar是二维的直柱状图，还有bar3（三维垂直柱状图，参数多了给分离式【detached】），barh（水平柱状图），bar3h（三位水平柱状图）

Y=round(rand(5,3)*10);%产生5*3的随机矩阵数据
subplot(221);
bar(Y,'grouped');%分组式
title 'Group'
subplot(222);
bar(Y,'stacked');%累加式
title 'Stack'
subplot(223);
barh(Y,'stacked');%水平柱状图，累加式
title 'Stack'
subplot(224);
bar(Y,1.5);%1.5的宽度，发生图形叠加
title 'Width=1.5'

调用：area(x,y)
说明：绘制 Y 对 X 的图，并填充 0 和 Y 之间的区域。X 的值可以是数值、日期时间、持续时间或分类值。
          如果 Y 是向量，则将 X 指定为由递增值组成的向量，其长度等于 Y。如果 X 的值不增加，则 area 将在绘制之前对值进行排序。
          如果 Y 是矩阵，则将 X 指定为由递增值组成的向量，其长度等于 Y 的行数。area 将 Y 的列绘制为填充区域。对于每个 X，最终结果是 Y 行的相应值的和。还可以将 X 指定为大小等于 Y 的矩阵。为了避免 X 为矩阵时出现意外输出，一般将 X 的列指定为重复列。

Y = [1, 5, 3;
     3, 2, 7;
     1, 5, 3;
     2, 6, 1];%构建矩阵
figure
area(Y)
set(gca,'layer','top');
title 'Area'

调用：fill(X,Y,C)
说明：根据 X 和 Y 中的数据创建填充的多边形（顶点颜色由 C 指定）。C 是一个用作颜色图索引的向量或矩阵。如果 C 为行向量，length(C) 必须等于 size(X,2) 和 size(Y,2)；如果 C 为列向量，length(C) 必须等于 size(X,1) 和 size(Y,1)。必要时，fill 可将最后一个顶点与第一个顶点相连以闭合多边形。X 和 Y 的值可以是数字、日期时间、持续时间或分类值。

t = (1/16:1/8:1)'*2*pi;
x = cos(t);
y = sin(t);
fill(x,y,'r')%创建一个红色八边形
axis square

调用：compass(z)或compass(u,v);feather(z)或feather(u,v)
说明：z代表复数向量，u和v分别代表复数向量的实部和虚部

%创建一个由随机矩阵的特征值构成的罗盘图。
rng(0,'twister') %初始化随机数生成
M = randn(20,20);
Z = eig(M);
figure
compass(Z)

theta = -pi/2:pi/16:pi/2;%将 theta 定义为介于 −2π 和 2π 之间的值。
r = 2*ones(size(theta));%定义 r 为与 theta 同样大小的向量。
[u,v] = pol2cart(theta,r);%创建显示 theta 的方向的羽毛图。由于 feather 使用笛卡尔坐标，使用 pol2cart 将 theta 和 r 转换为笛卡尔坐标。
feather(u,v)

调用：pie(X,explode,'label') 
说明：X 是向量。explode和x同长度的向量，决定是否从饼图中分离出对应的块，非零就表示该部分饼块要分离；label标注饼图的字符串数据 
备注：三维饼图用pie3，使用格式和pie的相同

x=[1 3 0.5 2.5 3];%为方便理解，直接设置了和为10的五个数，然后按百分比就分别是：10% 30% 5% 25% 30% 
explode=[0 1 0 0 1];%将第2、5块饼图进行偏移，即将两个3分别代表的饼块从原本的饼状图分离出来 
subplot(211); 
pie(x,explode);%直接绘制 
labels={'1','2','3','4','5'};%设置五块饼对应的标注 
subplot(212); 
pie(x,explode,labels);%直接绘制+标注 
colormap prism

调用：hist(x)、hist(x,nbins) 
说明：x是一个待划入统计的元素组，nbins是指定直方图分布间隔

x = randn(1000,3); %产生1000*3的随机数组合 
subplot(221) 
hist(x)%直接绘图，自动等距间隔 
subplot(222) 
hist(x,(-4:4)); 
subplot(223) 
hist(x,(-2:2)); 
subplot(224) 
hist(x,[-4 -2 -1.5 0 0.5 3]);