如何用matlab计算体积

tm038

matlab里面有专门的函数吗？

ChaChing

Ref: http://www.mathworks.com/matlabc ... of-a-3d-convex-hull

1. function [volume,area]=area3d(v)
   
2. % compute volume and area of a convex hull of points v
   
3. % Malcolm A. MacIver, 2003
   
4. 
   
5. [K volume]=convhulln(v);
   
6. 
   
7. %
   
8. % Basic formula for computing triangle area
   
9. % || = 2-norm, VN = vertex of triangle
   
10. % ||V1 X V2 + V2 X V3 + V3 X V1||/2
    
11. 
    
12. area= ...
    
13. sum(sqrt(sum(( ...
    
14. [v(K(:,1),2).*v(K(:,2),3) - v(K(:,1),3).*v(K(:,2),2) ...
    
15.   v(K(:,1),3).*v(K(:,2),1) - v(K(:,1),1).*v(K(:,2),3)  ...
    
16.   v(K(:,1),1).*v(K(:,2),2) - v(K(:,1),2).*v(K(:,2),1)] + ...
    
17. [v(K(:,2),2).*v(K(:,3),3) - v(K(:,2),3).*v(K(:,3),2) ...
    
18.   v(K(:,2),3).*v(K(:,3),1) - v(K(:,2),1).*v(K(:,3),3)  ...
    
19.   v(K(:,2),1).*v(K(:,3),2) - v(K(:,2),2).*v(K(:,3),1)] + ...
    
20. [v(K(:,3),2).*v(K(:,1),3) - v(K(:,3),3).*v(K(:,1),2) ...
    
21.   v(K(:,3),3).*v(K(:,1),1) - v(K(:,3),1).*v(K(:,1),3)  ...
    
22.   v(K(:,3),1).*v(K(:,1),2) - v(K(:,3),2).*v(K(:,1),1)]).^2,2))) ...
    
23.   /2;

复制代码