积极 AI Logo 积极 AI

用单个均匀分布生成球面上两个不同均匀点：正四面体构型

用单个均匀分布生成球面上两个不同均匀点：正四面体构型
几何画板及绘图内容加载中，请稍等……
2021-01-22 02:43
向量立体几何概率统计几何作图
该范例展示如何在球面上仅借助“一次均匀分布”的随机机制，构造得到两个彼此不同且在球面上均匀分布的点，并用正四面体体积参数 V=4 的构型来组织采样差异。核心数学概念包括球面几何、向量与旋转变换、均匀分布采样、以及通过对极值/约束条件的比较来确保两点的对应关系满足“均匀且不同”的性质。作图要点：从构造摘要提炼可见它使用了 Sphere/Circle 来建立球面与辅助圆，并用 PointIn 判定点在球面相关区域的位置；通过 Vector、Angle、RotateText 与 PolyLine/Polygon 进行几何关系与角度刻画；同时结合 If、Max、Min、Length、Sum 等条件与量度逻辑，计算与筛选两点的参数与可行范围；用 TableText/Textfield 将中间量（如 t、α、fmin3D/fmax3D、Zmin/Zmax 等）可视化，并用 CurveCartesian 生成相关曲线用于检查对应约束。适合高中到大学初级阶段的几何与概率统计/计算数学学习者，也适合作为GeoGebra中“球面均匀采样与立体几何构型验证”的课堂演示或作业拓展。

GeoGebra 命令及参数

  1. Polygon({((-π), ((-π)) / 2), ((-π), π / 2), (π, π / 2), (π, ((-π)) / 2)}) → q_2
  2. PointIn(q_2) → Atest_2
  3. Function(f01, (-π), π, ((-π)) / 2, π / 2) → f0
  4. Zip((x(a), y(a), f0(a)), a, Zmin_2) → fmin3D
  5. Zip((x(a), y(a), f0(a)), a, Zmax_2) → fmax3D
  6. Vector((1.3 * Z3), (1.5 * Z3)) → AVZ3r
  7. Vector((1.3 * Z3), (1.3 * Z3) + UnitVector[Sum[Zip[UnitVector[a - Z3], a, lP]]]) → AVZ3
  8. Vector((1.3 * Z2), (1.5 * Z2)) → AVZ2r
  9. Vector((1.3 * Z2), (1.3 * Z2) + UnitVector[Sum[Zip[UnitVector[a - Z2], a, lP]]]) → AVZ2
  10. Vector((1.3 * Z1), (1.5 * Z1)) → AVZ1r
  11. Vector((1.3 * Z1), (1.3 * Z1) + UnitVector[Sum[Zip[UnitVector[a - Z1], a, lP]]]) → AVZ1
  12. Vector((1.3 * Atest_R), (1.3 * Atest_R) + UnitVector[Sum[Zip[UnitVector[a - Atest_R], a, lP]]]) → AVtest
  13. Vector((1.3 * Atest_R), (1.5 * Atest_R)) → AVtestr
  14. Sphere(O, R) → sphera
  15. Circle(O, R, xOyPlane) → c
  16. PointIn(q_2) → RZ1_2
  17. PointIn(q_2) → RZ2_2
  18. Angle(AVtest, AVtestr) → αA
  19. Angle(AVZ1, AVZ1r) → αZ1
  20. Angle(AVZ2, AVZ2r) → αZ2
  21. Zip(Segment[Atest_R, p], p, lP) → seg
  22. Sum(seg) → sum_seg
  23. PolyLine({((-π), ((-π)) / 2), ((-π), π / 2), (π, π / 2), (π, ((-π)) / 2), ((-π), ((-π)) / 2)}) → q1
  24. Flatten(Sequence[Sequence[((-π) + ((2 * π / na) * i), ((-π)) / 2 + (π / na * j)), j, 1, na], i, 1, na]) → ltest_2
  25. Zip((R; x(a); y(a)), a, ltest_2) → ltest_3
  26. Length(lZmax') → lälZmax'
  27. Textfield(n) → InputBox11
  28. alfa(R, lP, lZsad') → αsad
  29. If(αZ1 > 100°, 180° - αZ1, αZ1) → αZ1'
  30. Textfield(m_1) → InputBox12
  31. Textfield(t) → InputBox13
  32. Textfield(m_2) → InputBox14
  33. PointIn(q_2) → RZ3_2
  34. Sum(Zip[1 / Length[p - Z3], p, lP]) → ZahZ3
  35. Angle(AVZ3, AVZ3r) → αZ3
  36. Textfield(R) → InputBox1
  37. Textfield(na) → InputBox2
  38. Max(f(Zmax_2)) → max_f
  39. Min(f(Zmin_2)) → min_f
  40. Zip(Text["d_{" + (a) + "}", (Atest_R + Element[lP, a]) / 2], a, 1…n) → lPtxt
  41. Zip(Segment[F_o, p], p, lP) → seg0
  42. Sum(seg0) → sum_seg0
  43. RotateText(txt, 90°) → txt1
  44. Sum(lZ3sad) → sumlZ3sad
  45. Textfield(Points3D) → InputBox7
  46. Textfield(lZsad') → InputBox17
  47. Textfield(lZmin') → InputBox16
  48. Textfield(lZmax') → InputBox15
  49. Textfield(jmin) → InputBox18
  50. Textfield(jmax) → InputBox19
  51. Textfield(jsad) → InputBox20
  52. Zip(Text["" + (b), a - (0.1, 0.1)], a, Zmax_2, b, 1…Length[Zmax_2]) → l2maxtxt
  53. Zip(Text["" + (b), a - (0.1, 0.1)], a, Zmin_2, b, 1…Length[Zmin_2]) → l2mintxt
  54. Zip(Text["" + (b), a - (0.1, 0.1)], a, Zsad_2, b, 1…Length[Zsad_2]) → l2sadtxt
  55. Zip(Vector[r1, (1.153846154 * r1)], r1, (1.3 * Zip[(R * UnitVector[a]), a, lZmax'])) → lVmaxr
  56. Zip(Vector[r1, (1.153846154 * r1)], r1, (1.3 * Zip[(R * UnitVector[a]), a, lZmin'])) → lVminr
  57. Zip(Vector[r1, (1.153846154 * r1)], r1, (1.3 * Zip[(R * UnitVector[a]), a, lZsad'])) → lVsadr
  58. alfa(R, lP, lZmax') → αmax
  59. alfa(R, lP, lZmin') → αmin
  60. arrow(R, lP, lZmax') → lVmax
  61. arrow(R, lP, lZmin') → lVmin
  62. arrow(R, lP, lZsad') → lVsad
  63. CurveCartesian((-π) + ((2 * π) * a), ((-π)) / 2, f0((-π) + ((2 * π) * a), ((-π)) / 2), a, 0, 1) → par1
  64. CurveCartesian((-π) + ((2 * π) * a), π / 2, f0((-π) + ((2 * π) * a), π / 2), a, 0, 1) → par2
  65. CurveCartesian((-π), ((-π)) / 2 + (π * a), f0((-π), ((-π)) / 2 + (π * a)), a, 0, 1) → par3
  66. CurveCartesian(π, ((-π)) / 2 + (π * a), f0(π, ((-π)) / 2 + (π * a)), a, 0, 1) → par4
  67. Polygon({((-π), ((-π)) / 2), ((-π), π / 2), (π, π / 2), (π, ((-π)) / 2)}) → q_3
  68. Textfield(achs) → InputBox3
  69. Circle(Atest_2, 0.2) → ci
  70. TableText(Append["i", Sequence[Length[lZsad']]], Append["\frac{N°}{levels} ", msad], Append["P_{sad }", lZsad'], Append["(f'x,f'y)", lZf'xysad], Append[" Δα≟ π ", αsad], "v") → text2c
  71. TableText(Append["i", Sequence[Length[αmax]]], Append["\frac{N°}{levels} ", mmax], Append["P_{Max }", lZmax'], Append["(f'x,f'y)", lZf'xymax], Append[" Δα≟ π ", αmax], "v") → text2a
  72. TableText(Append["i", Sequence[Length[αmin]]], Append["\frac{N°}{levels} ", mmin], Append["P_{Min }", lZmin'], Append["(f'x,f'y)", lZf'xymin], Append[" Δα≟ π ", αmin], "v") → text2b
  73. Zip((x(a), y(a), f0(a)), a, Zsad_2) → fsad3D
  74. Zip(Atest_2 + (r_c; ((2 * π / n_c) * k)), k, 0…(n_c - 1)) → l0r
  75. Sum(Sequence[If[sum00_r(i) ≠ sum00_r(i + 1), 1, 0], i, 1, Length[sum00_r] - 1]) → sum01_r
  76. Sum(Unique[sgn(f(Atest_2) - f(l0r))]) → sum_r

AI 绘图提示语参考

在三维球面上以同一中心与半径画出球体，并给出以球面相关区域为依据的可视化平面轮廓：在球心与球面基础圆所在的平面内画出一条封闭四边形折线（对应经纬参数取值的边界），同时在该参数面上再画出另一条细分的四边形边界折线用于检验。选取一个位于封闭四边形区域内的参考方向点，并在球面上构造与参考方向形成不同夹角的第二方向点：通过从参考点到球心的向量关系建立两条方向向量，再用夹角测量得到约束角的大小。最后在球面与参数面上对这些角与边界对应的可行区进行标注。
以上提示语为根据原始绘图命令整理而来，仅供参考。