<think>我们使用matplotlib和numpy来绘制3D图形。为了创建动态效果，我们可以使用matplotlib的动画模块。 思路： 1. 使用参数方程表示玫瑰曲线。玫瑰曲线的参数方程通常为：r = a * sin(n*θ) 或 r = a * cos(n*θ)，这里我们使用正弦函数。 在3D中，我们可以将其扩展到球坐标，然后转换为笛卡尔坐标。 2. 一个常见的3D玫瑰曲线方程（球形玫瑰）可以表示为： r = sin(8 * θ) # 这里8控制花瓣数量 然后我们使用球坐标到笛卡尔坐标的转换： x = r * sin(φ) * cos(θ) y = r * sin(φ) * sin(θ) z = r * cos(φ) 但是，为了形成花瓣，我们可以将r与φ有关，例如：r = a * sin(n * φ) 然后与θ结合形成旋转。 3. 实际上，我们可以使用一个更常见的3D玫瑰公式（例如Maarten函数）： r = cos(4 * φ) * (1 + cos(θ)) [这里只是一个例子，我们可以调整] 4. 另一种方法是使用参数方程，例如： x = cos(θ) * sin(φ) y = sin(θ) * sin(φ) z = cos(φ) 然后通过乘以一个关于φ的函数来形成花瓣（例如乘以sin(4*φ)等）。 这里我们采用一个经典的3D玫瑰曲线（球形玫瑰）方程： r = sin(4 * φ) # 4表示花瓣数量，可以调整 然后： x = r * cos(θ) * sin(φ) y = r * sin(θ) * sin(φ) z = r * cos(φ) 但是，为了形成多层花瓣，我们可以使用更复杂的函数。这里我们参考一个常见的“玫瑰”形状： r = 1 - sin(φ) # 这只是一个基础，我们可以乘以一个花瓣函数 实际上，我们可以使用以下方程（来自网络）来生成玫瑰： r = sin(8 * φ) # 8个花瓣 θ = θ 5. 为了动态效果，我们可以让花瓣随时间旋转，即改变θ的值。 然而，我们也可以使用另一种方法：使用参数t（时间）来改变角度，从而形成动画。 具体步骤： - 创建参数φ和θ的网格。φ从0到π，θ从0到2π。 - 计算每个点的r（半径）值，这里我们使用r = sin(4*φ)（4个花瓣）或者8个花瓣。 - 将球坐标转换为笛卡尔坐标。 - 绘制曲面。 - 使用FuncAnimation来旋转图形（通过改变视角）。 但是，为了更逼真的玫瑰，可能需要更复杂的方程。这里我们采用一个较为简单的版本，主要展示动态旋转。 另一种更逼真的方法是使用笛卡尔坐标下的参数方程，例如： x = cos(θ) * cos(φ) * r y = sin(θ) * cos(φ) * r z = sin(φ) * r 其中，r = a + b * cos(c * θ) * sin(d * φ) 等。 这里我们参考一个网上常见的3D玫瑰代码（简化版）： u, v = np.mgrid[0:2*np.pi:20j, 0:np.pi:10j] x = np.cos(u)*np.sin(v) y = np.sin(u)*np.sin(v) z = np.cos(v) # 然后乘以一个形状函数，例如： # x *= (1 + 0.5*np.cos(5*u)) # y *= (1 + 0.5*np.cos(5*u)) # z *= (1 + 0.5*np.cos(5*u)) 但是，这样得到的形状并不像玫瑰。因此，我们采用以下方程（一种常见的3D玫瑰参数方程）： u = t * np.pi (t从0到1) v = t * np.pi * 2 x = cos(v) * sin(u) y = sin(v) * sin(u) z = cos(u) # 然后乘以一个形状函数：r = 1 - (1 - mod(v, 2*np.pi/8) * 2) ** 2 等等，比较复杂。 实际上，我们可以使用一个已经验证的玫瑰曲线方程（来源：http://www.mathematische-basteleien.de/rose.htm）： 在2D中，玫瑰曲线：r = cos(k*θ) 在3D中，我们可以考虑将多个玫瑰曲线叠加。 这里，我们采用一个在网络上找到的3D玫瑰的代码（经过修改）： u = t * 2 * np.pi (t在0到1之间) v = t * 2 * np.pi k = 5 # 花瓣数量 r = 1 - 0.5 * (1 - np.mod(k*v, 2*np.pi)/np.pi) ** 2 x = r * cos(v) * sin(u) y = r * sin(v) * sin(u) z = r * cos(u) 然而，这个方程可能不太正确。我们尝试另一个版本（参考：https://stackoverflow.com/questions/41495658/use-matplotlib-to-create-3d-rose）： u = np.linspace(0, 2*np.pi, 100) v = np.linspace(0, np.pi, 100) u, v = np.meshgrid(u, v) # 玫瑰方程 r = 0.5 * (1 + 0.25 * np.sin(12 * v)) # 12个花瓣 x = r * np.cos(u) * np.sin(v) y = r * np.sin(u) * np.sin(v) z = r * np.cos(v) 这个方程会生成一个球体上有12个花瓣的玫瑰。但是花瓣是凸起的。 我们选择一种更常见的：使用球谐函数。不过，为了简单，我们使用一个已知的玫瑰参数方程： x = cos(u) * cos(v) y = sin(u) * cos(v) z = sin(v) # 然后乘以一个形状函数：r = 1 + 0.2 * cos(5*u) # 5个花瓣 因此，完整的方程为： r = 1 + 0.2 * cos(5*u) # 5个花瓣，振幅0.2 x = r * cos(u) * cos(v) y = r * sin(u) * cos(v) z = r * sin(v) 这样，在u方向上（经度）会形成5个花瓣。 为了动态效果，我们可以让u和v随时间增加一个相位，或者旋转整个图形。 这里，我们使用旋转视角的方法来实现动态效果。 步骤： 1. 创建数据：u, v 网格（u: 0~2π, v: -π/2~π/2，但通常v从0到π，这里我们使用v从0到π） 2. 计算r = 1 + 0.2 * np.cos(5*u) 3. 计算x, y, z 4. 绘制曲面 5. 创建动画函数，每次更新视角（例如增加azim角度） 注意：为了让玫瑰更立体，我们可以使用多个层，但这里我们简化。 代码实现： 我们将使用mpl_toolkits.mplot3d中的Axes3D，以及matplotlib.animation.FuncAnimation。 注意：由于绘制3D动画可能比较慢，我们适当减少网格点。 另外，为了美观，我们可以给花瓣上色（使用colormap）。 我们使用'viridis'颜色映射，但也可以选择红色系的，如'Reds' 由于玫瑰通常是红色的，我们可以设置颜色为红色系。 具体代码： 注意：在动画中，我们不会改变形状，而是旋转整个图形（通过改变ax.view_init的elev和azim角度）。这里我们只改变azim（水平旋转）。 但是，FuncAnimation中更新函数需要更新图形，我们可以重新绘制整个曲面，但这样效率低。另一种方法是只改变图形的旋转角度（通过set_azim）。 然而，mplot3d不支持直接更新视角的动画（通过set_azim）并重绘，我们可以使用ax.view_init(elev, azim)并在每次更新时设置azim，然后返回绘图对象。 参考：https://stackoverflow.com/questions/37872986/animate-3d-plot-using-matplotlib 我们采用每次更新视角的方法。 代码结构： - 创建图形和3D轴 - 生成网格数据 - 计算x,y,z - 绘制初始曲面 - 定义动画更新函数：每次改变azim角度，并设置新的视角 - 使用FuncAnimation生成动画 注意：azim角度从0到360循环。 但是，我们也可以同时改变elev角度，形成上下旋转，但这里只做水平旋转。 为了效果，我们让azim每帧增加1度。 由于我们只需要旋转，所以不需要重新计算网格数据，只需要改变视角。 因此，更新函数中，我们只需要更新ax.view_init中的azim。 但是，在FuncAnimation中，每次回调需要返回所有需要更新的对象。我们这里返回整个图形对象（ax）？实际上，我们需要返回绘图对象（plot_surface返回的对象）？但是视角改变并不改变数据，所以不需要更新数据点。我们只需要重绘。 实际上，我们可以不返回任何对象，然后使用fig.canvas.draw()，但这样不符合FuncAnimation的规范。我们可以返回一个空列表，或者返回我们设置视角的轴对象。实际上，我们可以返回一个空元组，然后通过blit=False来重绘整个图形。 由于3D图形的重绘比较慢，我们只改变视角，所以使用set_azim方法，然后返回整个轴对象？但是，mplot3d没有提供直接设置azim后返回可重绘对象的方法。 另一种方法是：在更新函数中，使用ax.view_init(elev=10, azim=frame)，然后返回ax对象？但是FuncAnimation要求返回一个可迭代的艺术家对象。我们可以返回一个我们之前绘制的曲面对象，但是曲面对象并没有改变，所以实际上不需要更新。因此，我们可以返回一个空列表，然后设置blit=False。 代码： 我们将使用： def update(frame): ax.view_init(elev=10, azim=frame) return fig, 然后在FuncAnimation中，我们传入这个函数，并设置frames=np.arange(0,360,1)，interval=50（每50ms一帧） 注意：初始视角我们可以设置elev=10（仰角10度），azim=0。 现在，我们写完整代码。 </think> 以下是使用Python实现动态3D玫瑰花的完整代码，基于Matplotlib和NumPy库： ```python import numpy as np import matplotlib.pyplot as plt from matplotlib import cm from matplotlib.animation import FuncAnimation from mpl_toolkits.mplot3d import Axes3D # 设置图形参数 fig = plt.figure(figsize=(10, 8)) ax = fig.add_subplot(111, projection='3d') ax.set_axis_off() # 隐藏坐标轴 # 玫瑰参数方程 u = np.linspace(0, 2 * np.pi, 100) # 花瓣角度 v = np.linspace(0, np.pi, 50) # 花瓣高度 u, v = np.meshgrid(u, v) # 玫瑰形状参数 k = 5 # 花瓣数量 a = 1 # 花瓣高度 b = 0.2 # 花瓣宽度 # 创建动态玫瑰函数 def update(frame): ax.cla() # 清除当前帧 ax.set_axis_off() # 动态参数：随时间变化的旋转和大小 rotation = frame * 0.02 scale = 0.8 + 0.2 * np.sin(frame * 0.05) # 玫瑰参数方程 (球坐标转笛卡尔坐标) r = scale * (a * np.sin(k * u) + b) * (1 - 0.5 * np.cos(v)) x = r * np.cos(u) * np.sin(v) y = r * np.sin(u) * np.sin(v) z = r * np.cos(v) # 应用旋转 rot_x = x * np.cos(rotation) - y * np.sin(rotation) rot_y = x * np.sin(rotation) + y * np.cos(rotation) # 绘制3D曲面 (使用红色系渐变色) surf = ax.plot_surface(rot_x, rot_y, z, rstride=1, cstride=1, cmap=cm.RdPu, alpha=0.8, edgecolor='none') # 设置视角和范围 ax.view_init(elev=30, azim=frame*0.5) ax.set_xlim([-1.5, 1.5]) ax.set_ylim([-1.5, 1.5]) ax.set_zlim([-1.5, 1.5]) return surf, # 创建动画 ani = FuncAnimation(fig, update, frames=200, interval=50, blit=True) plt.tight_layout() plt.show() # 如需保存动画（取消注释） # ani.save('3d_rose.gif', writer='pillow', fps=20) ``` ### 代码说明： 1. **参数方程**： - 使用球坐标转换公式：$x = r \cos u \sin v$, $y = r \sin u \sin v$, $z = r \cos v$ - 花瓣形状由 $r = (a \sin(ku) + b)(1 - 0.5 \cos v)$ 控制，其中$k$控制花瓣数量 2. **动态效果**： - 通过`FuncAnimation`实现旋转动画 - 花瓣大小随时间周期性变化（`scale`参数） - 视角自动旋转（`azim`参数随时间变化） 3. **视觉优化**： - 使用紫红色渐变 (`cm.RdPu`) 模拟玫瑰颜色 - 隐藏坐标轴增强视觉效果 - 边缘柔化处理 (`edgecolor='none'`) 4. **参数调整**： - 修改`k`值改变花瓣数量（建议5-8） - 调整`a`控制花瓣高度，`b`控制花蕊大小 - 修改`frame`和`interval`调整动画速度 ### 运行要求： - 安装依赖库：`pip install numpy matplotlib pillow` - 建议在Jupyter Notebook中运行以获得最佳交互效果 - 取消最后注释可将动画保存为GIF 此代码生成的花朵具有真实花瓣的曲面特征，通过参数方程控制花瓣的波状边缘和立体形态，动态旋转展示花朵的3D结构[^1]。