function scatterbar(x,y,z,scale)
    % 根据散点数据绘制柱状图
    % x,y,z是等长的向量
        hx = range(x)/scale;
        hy = range(y)/scale;
        n = numel(x);

        figure
        hold on
        for i = 1:n
            vert = Vertices(x(i),y(i),z(i));
            %规定顶点顺序
            fac = [1 2 3 4;2 6 7 3;4 3 7 8;1 5 8 4;1 2 6 5;5 6 7 8];
            %通过patch对象生成蓝色的立方体盒子
            h = patch('faces',fac,'vertices',vert,'FaceColor','b');
            %set(h,'FaceAlpha',0.25);    %设置立方体盒子透明度
        end
        grid on
        view(3)
        hold off

        %--------------------------------------------------
        % 求柱子顶点的子函数
        %--------------------------------------------------
        function Ver = Vertices(x,y,z)
            % 由长方体底面中心坐标求顶点坐标
            Ver = [x-hx y-hy 0
                x-hx y+hy 0
                x+hx y+hy 0
                x+hx y-hy 0
                x-hx y-hy z
                x-hx y+hy z
                x+hx y+hy z
                x+hx y-hy z];
            end
    end