Python实现动态3D立方体：旋转的3D立方体动画

个人简介：某不知名博主，致力于全栈领域的优质博客分享 | 用最优质的内容带来最舒适的阅读体验！文末获取免费IT学习资料！

🍅 文末获取更多信息 🍅 👇🏻 精彩专栏推荐订阅收藏 👇🏻

专栏系列	直达链接	相关介绍
书籍分享	点我跳转	书籍作为获取知识的重要途径，对于IT从业者来说更是不可或缺的资源。不定期更新IT图书，并在评论区抽取随机粉丝，书籍免费包邮到家
AI前沿	点我跳转	探讨人工智能技术领域的最新发展和创新，涵盖机器学习、深度学习、自然语言处理、计算机视觉等领域的研究进展和趋势分析。通过深入解读前沿技术、案例研究和行业动向，为读者带来关于人工智能未来发展方向和应用前景的洞察和启发。
Elasticsearch	点我跳转	详解 Elasticsearch 搜索和数据分析引擎
科技前沿	点我跳转	本档是关于科技和互联网的专栏，旨在为读者提供有趣、有用、有深度的科技资讯和思考。从多个角度探讨科技与人类生活的关系，包括但不限于科技趋势、产品评测、技术解读、行业观察、创业故事等内容。希望通过本栏，与读者分享科技的魅力和思考，让科技成为我们生活的一部分，而不仅仅是一个陌生的词汇。
Java之光	点我跳转	本栏将带领读者深入探索Java编程世界的种种奥秘。无论你是初学者还是资深开发者，这里都将为你提供丰富的Java知识和实用的编程技巧。
Linux学习日志	点我跳转	本专栏致力于探索Linux操作系统的各个方面，包括基础知识、系统管理、网络配置、安全性等。通过深入浅出的文章和实践指南，帮助读者更好地理解和应用Linux，提高系统管理和开发技能。无论你是初学者还是有经验的Linux用户，都能在本专栏中找到有用的信息和解决方案。
MySQL之旅	点我跳转	专栏将带领读者进入MySQL数据库的世界，探索其强大的功能和应用。我们将深入探讨MySQL的基本概念、SQL语言的应用、数据库设计与优化、数据备份与恢复等方面的知识，并结合实际案例进行讲解和实践操作。
精通Python百日计划	点我跳转	我们将引领你踏上一段为期100天的编程之旅，逐步深入了解和掌握Python编程语言。无论你是编程新手还是有一定基础的开发者，这个专栏都会为你提供系统而全面的学习路径，帮助你在短短100天内成为Python高手。

引言

3D动画在数据可视化和图形学中具有重要意义，能够生动地展示复杂的三维结构和运动。在这篇博客中，我们将使用Python来实现一个动态旋转的3D立方体。通过利用Matplotlib库，我们能够轻松创建和动画化3D立方体。

准备工作

前置条件

在开始之前，你需要确保你的系统已经安装了Matplotlib库。如果你还没有安装它，可以使用以下命令进行安装：

pip install matplotlib

Matplotlib是一个强大的Python绘图库，支持生成各种静态、动态和交互式的图形。

代码实现与解析

导入必要的库

我们首先需要导入Matplotlib库和其他必要的模块：

import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
from matplotlib.animation import FuncAnimation

初始化3D立方体

我们需要定义3D立方体的顶点和边：

# 定义立方体的顶点
vertices = np.array([
    [1, 1, 1],
    [1, 1, -1],
    [1, -1, 1],
    [1, -1, -1],
    [-1, 1, 1],
    [-1, 1, -1],
    [-1, -1, 1],
    [-1, -1, -1]
])

# 定义立方体的边
edges = [
    [0, 1], [0, 2], [0, 4], [1, 3], [1, 5],
    [2, 3], [2, 6], [3, 7], [4, 5], [4, 6],
    [5, 7], [6, 7]
]

绘制立方体

我们定义一个函数来绘制立方体：

def draw_cube(ax, vertices, edges):
    for edge in edges:
        points = vertices[edge]
        ax.plot3D(*zip(*points), color="b")

旋转立方体

我们定义一个旋转矩阵来旋转立方体：

def rotate(vertices, angle_x, angle_y, angle_z):
    # 旋转矩阵
    rotation_x = np.array([
        [1, 0, 0],
        [0, np.cos(angle_x), -np.sin(angle_x)],
        [0, np.sin(angle_x), np.cos(angle_x)]
    ])
    rotation_y = np.array([
        [np.cos(angle_y), 0, np.sin(angle_y)],
        [0, 1, 0],
        [-np.sin(angle_y), 0, np.cos(angle_y)]
    ])
    rotation_z = np.array([
        [np.cos(angle_z), -np.sin(angle_z), 0],
        [np.sin(angle_z), np.cos(angle_z), 0],
        [0, 0, 1]
    ])
    # 应用旋转
    return np.dot(vertices, rotation_x).dot(rotation_y).dot(rotation_z)

创建动画

我们使用FuncAnimation创建动画效果：

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

def update(frame):
    ax.cla()
    angle = np.radians(frame)
    rotated_vertices = rotate(vertices, angle, angle, angle)
    draw_cube(ax, rotated_vertices, edges)
    ax.set_xlim([-2, 2])
    ax.set_ylim([-2, 2])
    ax.set_zlim([-2, 2])

ani = FuncAnimation(fig, update, frames=360, interval=30)
plt.show()

完整代码

import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
from matplotlib.animation import FuncAnimation

# 定义立方体的顶点
vertices = np.array([
    [1, 1, 1],
    [1, 1, -1],
    [1, -1, 1],
    [1, -1, -1],
    [-1, 1, 1],
    [-1, 1, -1],
    [-1, -1, 1],
    [-1, -1, -1]
])

# 定义立方体的边
edges = [
    [0, 1], [0, 2], [0, 4], [1, 3], [1, 5],
    [2, 3], [2, 6], [3, 7], [4, 5], [4, 6],
    [5, 7], [6, 7]
]

# 绘制立方体函数
def draw_cube(ax, vertices, edges):
    for edge in edges:
        points = vertices[edge]
        ax.plot3D(*zip(*points), color="b")

# 旋转立方体函数
def rotate(vertices, angle_x, angle_y, angle_z):
    # 旋转矩阵
    rotation_x = np.array([
        [1, 0, 0],
        [0, np.cos(angle_x), -np.sin(angle_x)],
        [0, np.sin(angle_x), np.cos(angle_x)]
    ])
    rotation_y = np.array([
        [np.cos(angle_y), 0, np.sin(angle_y)],
        [0, 1, 0],
        [-np.sin(angle_y), 0, np.cos(angle_y)]
    ])
    rotation_z = np.array([
        [np.cos(angle_z), -np.sin(angle_z), 0],
        [np.sin(angle_z), np.cos(angle_z), 0],
        [0, 0, 1]
    ])
    # 应用旋转
    return np.dot(vertices, rotation_x).dot(rotation_y).dot(rotation_z)

# 创建动画
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

def update(frame):
    ax.cla()
    angle = np.radians(frame)
    rotated_vertices = rotate(vertices, angle, angle, angle)
    draw_cube(ax, rotated_vertices, edges)
    ax.set_xlim([-2, 2])
    ax.set_ylim([-2, 2])
    ax.set_zlim([-2, 2])

ani = FuncAnimation(fig, update, frames=360, interval=30)
plt.show()