border_vertex = [];
for i = 1:size(edges,1)
FI = find(sum(ismember(faces, edges(i,:)),2) == 2);
if( dot(normalf(FI(1),:),normalf(FI(2),:)) == -1)
border_vertex = [border_vertex edges(i,:)];
end
end
对于每一边, 然后对每个面执行ismember