%?此函數用來輸出位于邊界面上的點構成的面
function?[C]?=?faceinbound（A）
C=cell（1100）;?%存放所有滿足要求的面的元胞數組
bound1=0;
bound2=1;?%邊界為0到1的正方體
[nm]=size（A）;
B0=cell（1m）;%存放位于0邊界上的點
for?i=1:m
????B0{i}=[];
end
B1=cell（1m）;%存放位于1邊界上的點
for?i=1:m
????B1{i}=[];
end
for?i=1:n
????for?j=1:m
????????if?abs（A（ij）-bound1）<0.00001
????????????B0{j}=[B0{j};A（i:）];
????????end
????????if?abs（A（ij）-bound2）<0.00001
????????????B1{j}=[B1{j};A（i:）];
????????end
????end
end
k=1;
for?i=1:m
????if?size（B0{i}1）>=3?%3點以上構成一個面
????????C{k}=B0{i};
????????k=k+1;
????end
????if?size（B1{i}1）>=3
????????C{k}=B1{i};
????????k=k+1;
????end
end
[nm]=size（C）;
for?i=1:m
????if?isempty（C（i））?%獲取C中非空元素數量
????????break
????end
end
C=C（1:i）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????421??2015-01-25?16:48??intersectfordouble.m

?????文件?????333301??2015-01-25?16:28??test.xls

?????文件???????2095??2015-01-25?17:00??Voronoin_3D.m

?????文件????????821??2015-01-25?16:42??faceinbound.m

-----------?---------??----------?-----??----

???????????????336638????????????????????4