MATLAB自带有绘制三维柱状图的函数bar3，有版友反映不好控制柱子的颜色，特别是bar3不能绘制渐变色的柱子，并且bar3不能根据用户指定的散点数据绘制柱状图。为此，我编写了一个函数（scatterbar），用来绘制带渐变色的三维彩色柱状图。下面给出示例。

 

【例】根据散点数据绘制三维彩色柱状图

代码1：

 

[x,y] = meshgrid(-6:6,-3:0.5:3);

z = mvnpdf([x(:),y(:)],[0,0],[4,0;0,1]);

scatterbar(x,y,z)

效果图1：

http://attach.matlabsky.com/data/attachment/forum/201111/03/154251pyh5ux5csupxhhxu.jpg.thumb.jpg

代码2：

 

x = 0:0.1:2*pi;

y = sin(x);

scatterbar(x,y,y,100);

效果图2：

http://attach.matlabsky.com/data/attachment/forum/201111/03/154254umssno30csu4n0u0.jpg.thumb.jpg

代码3：

 

t = 0:0.1:2*pi;

x = cos(t);

y = sin(t);

z = abs(y);

scatterbar(x,y,z,100);

效果图3：

http://attach.matlabsky.com/data/attachment/forum/201111/03/1542566kp2b67om42sc48t.jpg.thumb.jpg

scatterbar函数源代码：

 

function scatterbar(x,y,z,scale)

%   根据散点数据绘制3维彩色柱状图

%   scatterbar(x,y,z,scale)  x,y,z是实值数组，用来指定柱子顶面中心点三维坐标。

%              scale是大于0的标量，用来指定柱子的粗细，scale越大，柱子越细，默

%              认情况下根据坐标点自动计算柱子的粗细。

%

%   CopyRight:xiezhh（谢中华）

%   2011.10.31

%   Example：

%       [x,y] = meshgrid(-6:6,-3:0.5:3);

%       z = mvnpdf([x(:),y(:)],[0,0],[4,0;0,1]);

%       scatterbar(x,y,z)

%       scatterbar(x,y,z,50);

 

% 输入参数类型判断

if nargin < 3

    error('至少需要三个输入参数');

end

if ~isreal(x) || ~isreal(y) || ~isreal(z)

    error('前三个输入应为实值数组');

end

 

% 提取x,y,z等长部分的元素

x = x(:);

y = y(:);

z = z(:);

n = min([numel(x) numel(y) numel(z)]);

x = x(1:n);

y = y(1:n);

z = z(1:n);

 

% 计算极差和差分值

rx = range(x);

ry = range(y);

dx = abs(diff(x));

dx = min(dx(dx>0));

dy = abs(diff(y));

dy = min(dy(dy>0));

 

% 自动计算柱子的粗细

if nargin == 3

    if ~isempty(dx)

        hx = dx/2;

    else

        hx = 0.5;

    end

    if ~isempty(dy)

        hy = dy/2;

    else

        hy = 0.5;

    end

end

 

% 根据用户输入参数scale计算柱子的粗细

if nargin == 4

    if ~isreal(scale) || scale < 0

        error('第四个输入应为正的标量');

    end

    if rx == 0 && ry == 0

        rx = 0.5*scale;

        ry = rx;

    elseif rx == 0 || ry == 0

        rx = max(rx,ry);

        ry = rx;

    end

    hx = rx/scale;

    hy = ry/scale;

end

 

% 通过循环绘制三维彩色柱状图

figure

hold on

Xp = [];

Yp = [];

Zp = [];

for i = 1:n

    [xp,yp,zp] = Vertices(x(i),y(i),z(i));

    Xp = [Xp;xp];

    Yp = [Yp;yp];

    Zp = [Zp;zp];

end

%通过surf函数生成彩色的立方体盒子

h = surf(Xp,Yp,Zp,Zp,'FaceColor','interp');

%set(h,'FaceAlpha',0.25);    %设置立方体盒子透明度

grid on

view(3)

hold off

 

%--------------------------------------------------

% 求柱子顶点的子函数

%--------------------------------------------------

function [xp,yp,zp] = Vertices(x,y,z)

    % 由长方体底面中心坐标求顶点坐标

    xp = [x-hx x-hx x+hx x+hx x-hx

          x-hx x-hx x+hx x+hx x-hx

          x-hx x-hx x+hx x+hx x-hx

          x-hx x-hx x+hx x+hx x-hx

          x    x    x    x    x

          NaN  NaN  NaN  NaN  NaN];

    yp = [y-hy y+hy y+hy y-hy y-hy

          y-hy y+hy y+hy y-hy y-hy

          y-hy y+hy y+hy y-hy y-hy

          y-hy y+hy y+hy y-hy y-hy

          y    y    y    y    y

          NaN  NaN  NaN  NaN  NaN];

    zp = [repmat(linspace(0,z,4)',[1,5]);z z z z z;NaN NaN NaN NaN NaN];

end

end

​