Question

2dshepp-logan模型matlab代码 不要用phantom函数

Answer 1

sp = Sphere[{0, 0, 0}, 1];
E0 = Scale[sp, {0.69, 0.92, 0.9}];
E1 = Scale[sp, {0.6624, 0.874, 0.880}];
E2 = Translate[
Rotate[Scale[sp, {0.4100, 0.16, 0.21}], -72 Degree, {0, 0,
1}], {-0.22, 0.0, -0.25}];
E3 = Translate[
Rotate[Scale[sp, {0.3100, 0.11, 0.22}],
72 Degree, {0, 0, 1}], {0.22, 0.0, -0.25}];
E4 = Translate[Scale[sp, {0.21, 0.25, 0.5}], {0.0, 0.35, -0.25}];
E5 = Translate[Scale[sp, {0.046, 0.046, 0.046}], {0.0, 0.1, -0.25}];
E6 = Translate[
Scale[sp, {0.046, 0.023, 0.020}], {-0.08, -0.605, -0.25}];
E7 = Translate[
Rotate[Scale[sp, {0.046, 0.023, 0.020}],
90 Degree, {0, 0, 1}], {0.06, -0.605, -0.25}];
E8 = Translate[
Rotate[Scale[sp, {0.056, 0.04, 0.1}],
90 Degree, {0, 0, 1}], {0.06, -0.105, 0.625}];
E9 = Translate[Scale[sp, {0.056, 0.056, 0.1}], {0.0, 0.1, 0.625}];
E10 = Translate[Scale[sp, {0.046, 0.046, 0.046}], {0.0, -0.1, -0.25}];
E11 = Translate[
Scale[sp, {0.023, 0.023, 0.023}], {0.0, -0.605, -0.25}];
Show[Graphics3D[{Opacity[.3], E0, E1, E2, E3, E4, E5, E6, E7, E8, E9,
E10, E11}]]