MATLAB-三维柱状图

matlab三维柱状图bar3画出来总不尽如人意，尤其遇到横纵坐标间隔大小不一致的情况。前段时间看到Mickey Stahl的scatterbar3函数，受此启发，在原程序上进行了小修改。

选取的x和y分为[20:10:300]和[1:1:10]。两个坐标轴跨度较大。

处理后的图片如下：

图中柱子的大小可以通过参数调节

 代码如下（交流侵删）：

x=[20:10:300];%横轴
y=[1:1:10];%纵轴
%P1是要画的数据
T=[];
for i=1:length(x)
    for j=1:length(y)
        XY=[x(i) y(j) P(i,j)];
        T=[T;XY];
    end
end
 figure;scatterbar3(T(:,1),T(:,2),T(:,3),5,0.3)   %“5”是横轴柱宽，“0.3”是纵轴柱宽
 view(45,45)
 
 
 function scatterbar3(X,Y,Z,width,length)

% By Mickey Stahl - 2/25/02
% Engineering Development Group
% Aspiring Developer
% By Lyu-11/7/22

[r,c]=size(Z);
for j=1:r,
    for k=1:c,
        if ~isnan(Z(j,k))
            drawbar(X(j,k),Y(j,k),Z(j,k),width/2,length/2)
        end
    end
end

zlim=[min(Z(:)) max(Z(:))];
if zlim(1)>0,zlim(1)=0;end
if zlim(2)<0,zlim(2)=0;end
axis([min(X(:))-width max(X(:))+width min(Y(:))-length max(Y(:))+length zlim])
caxis([min(Z(:)) max(Z(:))])
view(3)
 end

function drawbar(x,y,z,width,length)

h(1)=patch([-width -width width width]+x,[-length length length -length]+y,[0 0 0 0],'b');
h(2)=patch(width.*[-1 -1 1 1]+x,length.*[-1 -1 -1 -1]+y,z.*[0 1 1 0],'b');
h(3)=patch(width.*[-1 -1 -1 -1]+x,length.*[-1 -1 1 1]+y,z.*[0 1 1 0],'b');
h(4)=patch([-width -width width width]+x,[-length length length -length]+y,[z z z z],'b');
h(5)=patch(width.*[-1 -1 1 1]+x,length.*[1 1 1 1]+y,z.*[0 1 1 0],'b');
h(6)=patch(width.*[1 1 1 1]+x,length.*[-1 -1 1 1]+y,z.*[0 1 1 0],'b');
set(h,'facecolor','flat','FaceVertexCData',z)
end