各位,我把这个问题解决了。
主要思路是利用fill3画出各个单元,并根据各个单元每一个节点的值赋予颜色。
其中利用到循环。
clc
clear all
clf
a=load('element.txt')
b=load('vertices.txt')
c=load('velocity.txt')
a=a(:,2:5)
b=b(:,2:4)
for i=1:900
for j=1:4
x(j)=b(a(i,j),1);
y(j)=b(a(i,j),2);
z(j)=b(a(i,j),3);
cc(j)=c(a(i,j));
end
hold on
fill3(x,y,z,cc)
end
view(3)
axis on
axis equal