### 使用 `matplotlib.animation` 绘制 3D 动态图 要使用 Python 的 `matplotlib.animation` 库绘制基于 \( \text{sin}(\sqrt{x^2 + y^2} + t) \) 的 3D 动态图，可以按照以下方式实现： #### 实现思路 1. **导入必要的库**：需要引入 `numpy`, `matplotlib.pyplot` 和 `mpl_toolkits.mplot3d` 来支持三维绘图。 2. **定义数据范围**：创建二维网格用于表示 \(x\) 和 \(y\) 轴的数据点。 3. **计算函数值**：对于每一帧的时间步长 \(t\)，计算对应的 \(z = \sin(\sqrt{x^2 + y^2} + t)\) 值。 4. **设置动画更新函数**：利用 `FuncAnimation` 更新每一帧中的 \(z\) 数据并重新渲染表面。 以下是完整的代码示例： ```python import numpy as np from mpl_toolkits.mplot3d import Axes3D import matplotlib.pyplot as plt from matplotlib.animation import FuncAnimation # 定义xy平面的范围和分辨率 x = np.linspace(-5, 5, 50) y = np.linspace(-5, 5, 50) X, Y = np.meshgrid(x, y) # 初始化时间变量 T = 0 # 创建图形对象 fig = plt.figure(figsize=(8, 6)) ax = fig.add_subplot(111, projection='3d') def update(frame): global T T += 0.1 # 时间增量 Z = np.sin(np.sqrt(X**2 + Y**2) + T) # 计算新的Z值 ax.clear() # 清除当前轴的内容 surf = ax.plot_surface(X, Y, Z, cmap='viridis', edgecolor='none') # 设置坐标轴标签和其他属性 ax.set_title(f"Frame {frame}: $\\sin(\\sqrt{{x^2 + y^2}} + t)$", fontsize=14) ax.set_xlabel('X axis') ax.set_ylabel('Y axis') ax.set_zlabel('Z axis') ax.set_xlim([-5, 5]) ax.set_ylim([-5, 5]) ax.set_zlim([-1, 1]) return surf, # 创建动画 ani = FuncAnimation(fig, update, frames=range(100), interval=50, blit=False) plt.show() ``` --- #### 关键点解析 1. **网格化数据** 利用 `np.meshgrid` 将一维数组扩展为二维网格，便于后续计算 \(z\) 值[^1]。 2. **动态更新机制** 每次调用 `update()` 函数时，都会增加全局变量 \(T\) 并重新计算 \(z\) 值。这使得整个过程具有动态效果[^3]。 3. **清除与重绘** 在每次更新前清空之前的绘图内容 (`ax.clear()`)，从而避免旧数据残留影响新画面的效果[^2]。 4. **色彩映射** 使用 `'viridis'` 颜色映射来增强视觉表现力，使波形变化更加直观[^1]。 --- #### 注意事项 - 如果运行速度较慢，可以通过减少网格密度（如将 `linspace` 中的点数降低至 20 或更少）优化性能。 - 对于保存 GIF 文件的情况，需安装额外工具包（如 ImageMagick），并通过 `ani.save("output.gif", writer="imagemagick")` 方法完成导出操作[^1]。 --- ###