搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  nodetemp.txt
资料下载链接地址: https://bbs.pinggu.org/a-1641229.html
附件大小:


matlab如何画四维图像。







〈本帖收录于〉——〖素质文库〗

(三大专题:〖Only one〗、〖On Spss〗、〖On Matlab〗)







法一


这个问题相当于有限元里面后处理的例子。一个数组一共四列,前三列为坐标xyz,第四列为温度。



  1. clear
  2. clc
  3. data=load('c:/nodetemp.txt')
  4. sx=0:0.2/5:0.2
  5. sy=0:0.3/5:0.3
  6. sz=0:0.2/5:0.2
  7. v=zeros(6,6,6)
  8. for index=1:216
  9. if (data(index,1)-0.00)<1e-5
  10. xbiao=1
  11. elseif (data(index,1)-0.04)<1e-5
  12. xbiao=2
  13. elseif (data(index,1)-0.08)<1e-5
  14. xbiao=3
  15. elseif (data(index,1)-0.12)<1e-5
  16. xbiao=4
  17. elseif (data(index,1)-0.16)<1e-5
  18. xbiao=5
  19. elseif (data(index,1)-0.2)<1e-5
  20. xbiao=6
  21. end

  22. if (data(index,2)-0.00)<1e-5
  23. ybiao=1
  24. elseif (data(index,2)-0.06)<1e-5
  25. ybiao=2
  26. elseif (data(index,2)-0.12)<1e-5
  27. ybiao=3
  28. elseif (data(index,2)-0.18)<1e-5
  29. ybiao=4
  30. elseif (data(index,2)-0.24)<1e-5
  31. ybiao=5
  32. elseif (data(index,2)-0.3)<1e-5
  33. ybiao=6
  34. end

  35. if (data(index,3)-0.00)<1e-5
  36. zbiao=1
  37. elseif (data(index,3)-0.04)<1e-5
  38. zbiao=2
  39. elseif (data(index,3)-0.08)<1e-5
  40. zbiao=3
  41. elseif (data(index,3)-0.12)<1e-5
  42. zbiao=4
  43. elseif (data(index,3)-0.16)<1e-5
  44. zbiao=5
  45. elseif (data(index,3)-0.2)<1e-5
  46. zbiao=6
  47. end
复制代码

读入的txt文件:


[hide][/hide]


结果图形:








法二:


这样的方法没代表性吧。 说白了,就是把这些node练成多边形,然后根据结点的温度填上结点的颜色。因为软件填色的本身设定,整个多边形就填上色。Matlab 可以做,Mathematica也可以做。因为你把那些结点关系省略了,我就用Interpolation来补充一下了。


data = Import["C/nodetemp.txt", "Table"]

data1 = Transpose[{data[[All, 1 ;; 3]], data[[All, 4]]}];

b = Max@data[[All, 4]];

a = Min@data[[All, 4]];

ff = Interpolation[data1, InterpolationOrder -> 1]

tt = RegionPlot3D[

0 <= x <= 0.2 && 0 <= y <= 0.3 && 0 <= z <= 0.2, {x, 0, 0.2}, {y,

0, 0.3}, {z, 0, 0.2}, Mesh -> {8, 5, 6}, PlotPoints -> {10, 7, 8},

MaxRecursion -> 0];

cls = Map[(ColorData["TemperatureMap"][(ff @@ # - a)/(b - a)]) &,

tt[[1, 1]]];

Row[{Graphics3D[{GraphicsComplex[

tt[[1, 1]], {Opacity[0.8], EdgeForm[],

Cases[tt[[1]], Polygon[__], {0, Infinity}]},

VertexColors -> cls]}, Axes -> True,

AxesLabel -> {"x", "y", "z"}, Boxed -> False, BoxRatios -> 1,

ViewPoint -> {-2, -2, 1}, ImageSize -> 350],

Graphics[DensityPlot[(y - a)/(b - a), {x, 0, 0.5}, {y, 22, 25},

ColorFunction -> (ColorData["TemperatureMap"][#1] &)][[1]],

Ticks -> {None, Range[22, 25, 0.5]}, Frame -> False, Axes -> True,

AspectRatio -> 10, ImageSize -> 45]}, Spacer[5]]








    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

GMT+8, 2025-12-26 18:29