Matlab中怎样画这样复杂的柱状图？？

zxb1219

17 分钟前 上传 下载附件 (97.5 KB) bone cartilage 0.5 97.9916 0 1 98.8735 0 2.5 86.232 13.6252 5 39.6524 60.3476 bone cartilage 0.5 69.3858 12.9519 1 75.6398 12.9519 2.5 71.6691 22.3528 5 58.1054 37.7719 bone cartilage 0.5 99.9521 0 1 99.9581 0 2.5 78.7006 21.2575 5 59.3743 40.5838 bone cartilage 0.5 99.2139 0 1 99.2268 0 2.5 81.8138 17.413 5 60.6121 38.6147 如图所示：我想解决的问题是：1.把这四种工况的数据分成4组显示：即0.5这一组有4个柱子，1这组也是4个柱子......每一组的柱子之间没有空隙，而组与组的柱子之间    是有空隙的                                                  2.每个柱子都是堆叠图，堆叠的上下两部分用不同的图案来表示。并添加标签：如图上的cartilage和bone                                                  3.X轴只显示0.5，1，2.5，5这样4个数据                                                  4.每一种工况都用图例标签显示出来

liuyalong008

1. function mybargroup(d)
   
2. %%
   
3. %{
   
4. d =
   
5. 
   
6.     0.5000   97.9916         0
   
7.     1.0000   98.8735         0
   
8.     2.5000   86.2320   13.6252
   
9.     5.0000   39.6524   60.3476
   
10.     0.5000   69.3858   12.9519
    
11.     1.0000   75.6398   12.9519
    
12.     2.5000   71.6691   22.3528
    
13.     5.0000   58.1054   37.7719
    
14.     0.5000   99.9521         0
    
15.     1.0000   99.9581         0
    
16.     2.5000   78.7006   21.2575
    
17.     5.0000   59.3743   40.5838
    
18.     0.5000   99.2139         0
    
19.     1.0000   99.2268         0
    
20.     2.5000   81.8138   17.4130
    
21.     5.0000   60.6121   38.6147
    
22. %}
    
23. [~,i]=sort(d)
    
24. data=d(i(:,1),:)
    
25. x=sum(data(:,1)==data(1))
    
26. xdata=1:x*(x+1);
    
27. xtick=mean(reshape(xdata,x+1,[]));
    
28. xdata(x+1:x+1:end)=[];
    
29. xdata=reshape(repmat([xdata;xdata+1;xdata+1;xdata]',1,2)',4,[])'
    
30. ydata=reshape([zeros(length(data(:,1)),2),repmat(data(:,2),1,2),...
    
31.     repmat(data(:,2),1,2) repmat(data(:,2),1,2)+repmat(data(:,3),1,2)]',4,[])'
    
32. cdata=repmat(colormap(jet(2*x)),x,1)
    
33. arrayfun(@(i)patch(xdata(i,:),ydata(i,:),cdata(i,:)),1:size(xdata,1))
    
34. set(gca,'xtick',xtick,'xticklabel',unique(data(:,1)))
    
35. axis tight
    
36. box on
    
37. grid on
    
38. h=colorbar
    
39. set(h,'ytick',linspace(1/(4*x),1-1/(4*x),2*x),'yticklabel',1:2*x)
    
40. 
    

复制代码