用matlab读取obj格式文件并求每一个三角形法向量

以下是一个读取obj格式文件并求每个三角形法向量的例子：

filename = 'example.obj';

fid = fopen(filename, 'r');

vertices = [];
normals = [];
faces = [];

while ~feof(fid)
    line = fgetl(fid);
    
    if line(1) == 'v'
        vertex = sscanf(line, 'v %f %f %f');
        vertices = [vertices; vertex'];
    elseif line(1) == 'f'
        face = sscanf(line, 'f %d//%d %d//%d %d//%d');
        faces = [faces; face'];
    end
end

fclose(fid);

num_faces = size(faces, 1);
tri_normals = zeros(num_faces, 3);

for i = 1:num_faces
    v1 = vertices(faces(i,1), :);
    v2 = vertices(faces(i,3), :);
    v3 = vertices(faces(i,5), :);
    
    edge1 = v2 - v1;
    edge2 = v3 - v1;
    
    normal = cross(edge1, edge2);
    normal = normal / norm(normal);
    
    tri_normals(i, :) = normal;
end

这个例子中，我们首先打开文件并读取其中的顶点、法向量和面数据。然后我们逐个迭代每个三角形面，并计算它的法向量。计算法向量的方法是通过叉积计算两个边的向量，然后标准化得到单位向量。最后，我们将每个三角形面的法向量存储在一个数组中。