3D packing 问题

使用Zoutendijk方法解决3D packing问题的matlab代码

function [ C t ] = combsin（ a b ）
% a*sin（s） + b*cos（s）

C = sqrt（a^2 + b^2）;

t1（1） = asin（b/C）;
if t1（1）>=0
    t1（2） = pi-t1（1）;
else
    t1（2） = -pi-t1（1）;
end

t2（1） = acos（a/C）;
t2（2） = -t2（1）;

if isempty（ find（t2 == t1（1） 1） ）
    t = t1（2）;
    return
else
    t = t1（1）;
    return
end

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5317  2013-11-18 14:29  3D with rotation\backup.txt

     文件        330  2013-11-15 13:45  3D with rotation\combsin.m

     文件        102  2013-11-11 21:42  3D with rotation\diff_val.m

     文件        233  2013-11-11 21:45  3D with rotation\fmin.m

     文件         77  2013-11-11 22:26  3D with rotation\fval.m

     文件       3387  2013-11-25 15:53  3D with rotation\genCollision.m

     文件         65  2013-11-15 14:12  3D with rotation\geneGpara.m

     文件        583  2013-11-25 14:09  3D with rotation\genFunction.m

     文件       1989  2013-11-25 14:04  3D with rotation\genGi.m

     文件       1139  2013-11-25 18:51  3D with rotation\main.m

     文件       1522  2013-11-20 23:19  3D with rotation\maxStep.m

     文件        587  2013-11-15 15:10  3D with rotation\maxStep5.m

     文件        740  2013-11-13 09:52  3D with rotation\numericSolve.m

     文件        681  2013-11-18 15:53  3D with rotation\plotRec.txt

     文件       2288  2013-11-25 14:59  3D with rotation\plotResult.m

     文件       1061  2013-11-25 15:53  3D with rotation\rt.m

     文件       1721  2013-11-25 14:52  3D with rotation\test.m

     文件       6506  2013-11-14 17:31  3D with rotation\zoutendijk.m

     文件       2949  2013-11-25 18:56  3D with rotation\zoutendijk1.m

     文件       5986  2013-11-14 20:10  3D with rotation\zoutendijk2.m

     文件       6189  2013-11-15 15:10  3D with rotation\zoutendijk5.m

     文件       3759  2013-11-19 20:00  3D with rotation\zoutendijk6.m

     文件       2794  2013-11-25 15:42  3D with rotation\zoutendijk7.m

     目录          0  2013-11-25 20:22  3D with rotation

----------- ---------  ---------- -----  ----

                50005                    24