〖绘 • Matlab〗画四维图的例子

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

      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.                
    
23.              if (data(index,2)-0.00)<1e-5
    
24.                 ybiao=1
    
25.             elseif (data(index,2)-0.06)<1e-5
    
26.                ybiao=2
    
27.             elseif (data(index,2)-0.12)<1e-5
    
28.                ybiao=3
    
29.              elseif (data(index,2)-0.18)<1e-5
    
30.                ybiao=4  
    
31.             elseif (data(index,2)-0.24)<1e-5
    
32.                ybiao=5
    
33.             elseif (data(index,2)-0.3)<1e-5
    
34.                ybiao=6
    
35.              end
    
36.             
    
37.             if (data(index,3)-0.00)<1e-5
    
38.                 zbiao=1
    
39.             elseif (data(index,3)-0.04)<1e-5
    
40.                zbiao=2
    
41.             elseif (data(index,3)-0.08)<1e-5
    
42.               zbiao=3
    
43.              elseif (data(index,3)-0.12)<1e-5
    
44.                zbiao=4  
    
45.             elseif (data(index,3)-0.16)<1e-5
    
46.                zbiao=5
    
47.             elseif (data(index,3)-0.2)<1e-5
    
48.                zbiao=6
    
49.             end

复制代码

读入的txt文件：

本帖隐藏的内容

nodetemp.txt (8.86 KB)

2014-9-30 20:04:51 上传

结果图形：

---

     法二：

     这样的方法没代表性吧。 说白了，就是把这些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]]

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB matla atlab Atl Lab matlab 有限元 如何

读入的txt文件：

相当好啊                                                                        
                                       
                                                     
                                                     
                                             
                                                                 
                                                                                 

   不错的分享，学习一下。

谢谢分享，非常感谢。

vergy good............................

厉害
学无止境啊

真好！！！！！

代码少了一半啊？？？？还有后面一个是什么软件的代码？

这个敢情非不错，对维度作图有了一个全新的认识