小白求助各位大神

v = Tuples[{-1, 1}, 4];
e = Select[Subsets[Range[Length[v]], {2}], Count[Subtract @@ v[[#]], 0] == 3 &];
f = Select[Union[Flatten[#]] & /@ Subsets[e, {4}], Length@# == 4 &];
f = f /. {a_, b_, c_, d_} :> {b, a, c, d};
rotv[t_] = (RotationMatrix[ t, {{0, 0, 1, 0}, {0, 1, 0, 0}}].RotationMatrix[
2 t, {{1, 0, 0, 0}, {0, 0, 0, 1}}].#) & /@ v;
proj[t_] := Most[#]/(3 - Last[#]) & /@ rotv[t];
Animate[Graphics3D[GraphicsComplex[
proj[t], {Orange, Specularity[0.75, 10], Sphere[Range[16], 0.05],
Tube[e, 0.03], Opacity[0.3], Polygon@f}], Boxed -> False,
Background -> Black, ImageSize -> 390, PlotRange -> 1], {t, 0.0,
Pi/2.0, 0.075}]
有没有大神能简述一下这段代码是如何生成不停运动的超正方体

第一句是生成一个四元组，由-1和1轮换而成。
e = Select[Subsets[Range[Length[v]], {2}], Count[Subtract @@ v[[#]], 0] == 3 &];
第二句Subsets[Range[Length[v]], {2}]生成{1,2,3,4}的族集，但是元素个数不超过2
Count[Subtract @@ v[[#]], 0] == 3 &是一个函数，#是Subsets[Range[Length[v]], {2}]的元素之一v[[#]]是四元组里至多两个的元素Substract是让他们相减，不一定能减成。
Count是数里面有几个零。Count[Subtract @@ v[[#]], 0] == 3 &合起来就是若零的个数为3就是真。
Select就是选出这些含两个点邻接的集合。
e估计是边集的意思。

第三句估计就是面集合了
Union[Flatten[#]] & /@ Subsets[e, {4}]
是取最多四条边，然后把各边合起来的点集合并一下，也就是重复的删掉。
Length@# == 4 & 加上select就是选出有四个点的那种。如果是5个点以上就意味着不是个面了。
f是四个点构成的面集合
第四句是换一下位置，让第一个点和第二个点换一下。不太清楚为啥。

rotv[t_] = (RotationMatrix[ t, {{0, 0, 1, 0}, {0, 1, 0, 0}}].RotationMatrix[
2 t, {{1, 0, 0, 0}, {0, 0, 0, 1}}].#) & /@ v;
proj[t_] := Most[#]/(3 - Last[#]) & /@ rotv[t];
rotv应该是个旋转公式吧，proj应该是投影吧。这个书上自己查吧。
#和&看成整体。
/@
@@
的看成一个整体。
Animate就不说了

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频！

4回复贴，共1页

<<返回mathematica吧

分享到:

日	一	二	三	四	五	六