Mathematica中文论坛-非官方

标题: 豆粑粑，matab 我画你个饼 ，3D等值面 [打印本页]
作者: meatball1982    时间: 2018-1-16 14:02
标题: 豆粑粑，matab 我画你个饼 ，3D等值面
本帖最后由 meatball1982 于 2018-1-16 14:04 编辑

之前，陌生女人的提问，
解答过程中觉得思路还不错（当时是在车上，临时的一个小程序。
今天，突然想到可以用patch的方法画各个3D的等值面。
效果，还不错。



  1. clear all
  2. clc
  3. clf

  5. %% outline
  6. % plot BING
  7. % contact meatball1982@163.com
  8. % if u really, really got problem with this m file.


  11. %% main
  12. % generate data
  13. [x,y,z]=peaks(100);
  14. % get the contour line
  15. [C,h]=contour(x,y,z,18);

  17. % get the structure, close the contour
  18. S=contourdata(C);
  19. clf


  22. %% plot BING
  23. for i = 1: length(S)
  24.     s_lev(i,1)=S(i).level;
  25. end

  27. tm_s = unique(s_lev);
  28. n_col = length(tm_s);
  29. [loca,locb]=ismember(s_lev,tm_s); % get the lin color position in col_mm
  30. col_mm=jet(n_col);

  32. %% outplut
  33. % the normal one
  34. subplot(1,2,1)
  35. surf(x,y,z,'edgecolor','none')
  36. colormap(col_mm)
  37. view(-20,20)
  38. axis tight
  39. colorbar
  40. % grace BING plot
  41. subplot(1,2,2)
  42. hold on
  43. for i = 1:length(S)
  44.     x_lin = S(i).xdata;
  45.     y_lin = S(i).ydata;
  46.     z_va  = S(i).level;
  47.     z_lin = z_va.*ones(size(x_lin));
  48.     alp_va = abs(z_va./max(abs(z(:))))/2+0.5;
  49. %     h=patch(x_lin,y_lin,z_lin,col_mm(locb(i),:),...
  50. %                                          'edgecolor',col_mm(locb(i),:)); % no transpancy
  51.     h=patch(x_lin,y_lin,z_lin,col_mm(locb(i),:),...
  52.                                           'edgecolor',col_mm(locb(i),:),...
  53.                                           'facealpha',alp_va); % transpancy
  54. end
  55. axis tight
  56. view(-20,20)
  57. grid on
复制代码



work_plot_BING.tar.gz

1.31 KB, 下载次数: 4



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