Mathematica中文论坛-非官方

标题: matlab surf 维数不对  [打印本页]

作者: meatball1982    时间: 2016-3-2 14:18
标题: matlab surf 维数不对 
本帖最后由 meatball1982 于 2016-3-2 14:24 编辑
  1. clear all
  2. clc
  3. close all

  4. [x,y,z]=peaks(60);

  5. h=bar3(z,1)
  6. for i = 1:length(h)
  7.     zdata = get(h(i),'ZData');
  8.     set(h(i),'CData',zdata)
  9.     % Add back edge color removed by interpolating shading
  10.     set(h,'EdgeColor','w')
  11.     alpha(h,0.9)
  12. end

  13. axis equal
  14. grid off
  15. view(0,-90)
  16. hidden off
  17. colorbar
  18. % axis off

  19. surf(z,'edgecolor','w')
  20. axis equal
  21. view(0,90)
  22. hidden off
  23. colorbar
复制代码


之前,画3D的平面都是用surf,然后再view(0,90)
一个小问题,当数据维数不大的时候,无太多差异。



这是两组60*60的数据,上面是bar3,下面是surf。

当数据比较少的时候,如果只有6*6





差异就明显了,上图bar3,每行有6个格,下图surf,每行有5个格,是两个值的差异来画的。


所以,当数据多的时候,或是只关注曲面的变化形势时,用surf没问题,esle,如果要说明一个图有多少个点,比如神经网络的节点个数,则最好使用bar3。

更新:
  1. col=z;

  2. h=imagesc(col);
  3. axis equal
  4. axis tight
  5. % k=[0 max(z(:))-min(z(:))]
  6. k=[min(z(:)) max(z(:))]
  7. colormap(jet)
  8. caxis(k);
  9. view(0,90)
  10. colorbar
复制代码




又更新: 
02-Mar-2016 14:13:41
每次用surf画图都少一行,少一列,也不是那个事儿.
于是,写了一个小函数.生成新的x,y,z,可以直接surf,然后,维数还是正确的.
  1. function [ X,Y,z_new ] = fun_mm_shift( z )
  2. %UNTITLED Summary of this function goes here
  3. %   Detailed explanation goes here

  4. [m,n]=size(z);
  5. [X,Y]=meshgrid([1:(n+1)]-0.5,[1:(m+1)]-0.5);
  6. tm=z;
  7. tm=[z;z(end,:)];
  8. z_new=[tm';tm(:,end)']';

  9. end
复制代码


把z扔里面.就出来对应的新的x,y,z了.老好用了.






欢迎光临 Mathematica中文论坛-非官方 (http://ilovemathematica.com/) Powered by Discuz! X3.2