我现在定义并多面体棱边并矢E=nfinei+nfjnej,其中nfi,nfj为一条棱边所属两个平面的外法向量，nei和nej为该棱边两对于两个平面的外法向量，用matlab读取boj格式文件，并计算每一条棱边的E

可以使用MATLAB中的函数readObj来读取BOJ格式文件，并使用向量叉乘来计算每一条棱边的E。具体的代码如下：

% 读取BOJ格式文件
[V,F] = readObj('model.boj');

% 计算每一条棱边的E
E = zeros(size(F,1)*3,3);
for i = 1:size(F,1)
    % 获取该面的三个顶点坐标
    v1 = V(F(i,1),:);
    v2 = V(F(i,2),:);
    v3 = V(F(i,3),:);
    
    % 计算该面的两个外法向量
    n1 = cross(v2-v1,v3-v1);
    n2 = cross(v3-v2,v1-v2);
    
    % 计算每一条棱边的E
    E((i-1)*3+1,:) = cross(n1,n2)/norm(cross(n1,n2));
    E((i-1)*3+2,:) = cross(n1,n2)/norm(cross(n1,n2));
    E((i-1)*3+3,:) = cross(n1,n2)/norm(cross(n1,n2));
end

解释一下上面的代码：

1. 使用函数readObj读取BOJ格式文件，返回顶点坐标矩阵V和面索引矩阵F。

2. 遍历每一个面，计算该面的两个外法向量n1和n2，使用向量叉乘计算每一条棱边的E。

3. 将所有计算出的E保存在矩阵E中，其中每一行对应一条棱边。注意，由于一个面有三条边，因此每个面会计算出三条边的E，所以矩阵E的行数应该是面数的三倍。