### 李萨如图形动画绘制 以下是基于 `matplotlib` 的 Python3 示例代码，用于生成李萨如图形的动态动画。此代码通过调整参数 $p$ 和 $q$ 控制两个正弦波的比例关系，从而形成不同的李萨如图案。 #### 动态二维李萨如图 以下代码实现了动态更新的二维李萨如图： ```python import matplotlib.pyplot as plt import numpy as np from matplotlib.animation import FuncAnimation # 创建画布和子图 fig, ax = plt.subplots() ax.set_xlim([-1.2, 1.2]) ax.set_ylim([-1.2, 1.2]) ax.set_title('Dynamic Lissajous Figure') # 参数设置 theta = np.linspace(0, 2 * np.pi, 1000) p = 3 # 调整比例因子 p (x 方向频率倍数) q = 4 # 调整比例因子 q (y 方向频率倍数) # 初始化数据 x = np.sin(p * theta) y = np.sin(q * theta) line, = ax.plot(x, y) def update(frame): global p, q phase_shift = frame / 50 # 相位偏移量随帧变化 x = np.sin(p * theta + phase_shift) y = np.sin(q * theta) line.set_data(x, y) return line, ani = FuncAnimation(fig, update, frames=np.linspace(0, 2 * np.pi, 100), interval=50) plt.show() # 如果需要保存为 GIF 文件，可以取消下面这行注释 # ani.save('lissajous_animation.gif', writer='pillow') ``` 这段代码利用了 `FuncAnimation` 函数来创建动态效果[^1]。每一帧都会重新计算相位差并更新图像中的曲线位置。 --- #### 静态三维李萨如图 如果希望扩展到静态三维空间，则可以通过如下方式实现： ```python import matplotlib.pyplot as plt import numpy as np from mpl_toolkits.mplot3d import Axes3D # 设置绘图环境 fig = plt.figure(figsize=(8, 6)) ax = fig.add_subplot(111, projection='3d') # 数据准备 t = np.linspace(0, 2 * np.pi, 1000) A_x, A_y, A_z = 1, 1, 1 # 幅度系数 f_x, f_y, f_z = 3, 4, 5 # 频率比例因子 x = A_x * np.sin(f_x * t) y = A_y * np.sin(f_y * t) z = A_z * np.sin(f_z * t) # 绘制轨迹 ax.plot(x, y, z, label="Lissajous Curve", color="blue") # 图形属性配置 ax.set_xlabel("X Axis") ax.set_ylabel("Y Axis") ax.set_zlabel("Z Axis") ax.legend(loc="upper left") plt.title("Static 3D Lissajous Figure") plt.show() ``` 该部分展示了如何构建一个简单的三维李萨如图，并允许进一步自定义颜色、标签和其他样式选项[^2]。 --- #### 动态三维李萨如图 为了制作更加复杂的动态三维版本，可采用交互模式配合定时刷新机制完成渲染过程： ```python import matplotlib.pyplot as plt import numpy as np from mpl_toolkits.mplot3d import Axes3D # 开启交互模式 plt.ion() # 初始化图表对象 fig = plt.figure() ax = fig.add_subplot(111, projection='3d') # 时间变量范围 t = np.linspace(0, 2 * np.pi, 1000) # 初始条件设定 freq_ratio_x, freq_ratio_y, freq_ratio_z = 3, 4, 5 amp_x, amp_y, amp_z = 1, 1, 1 for i in range(100): # 循环模拟时间流逝 phi_t = i / 10 # 当前时刻对应的相角增量 # 更新坐标系下的三个维度分量 x = amp_x * np.sin(freq_ratio_x * t + phi_t) y = amp_y * np.sin(freq_ratio_y * t) z = amp_z * np.sin(freq_ratio_z * t) # 清除旧的数据点以便重绘新路径 ax.clear() ax.plot(x, y, z, c="green", lw=1.5) ax.set_xlim(-1.5, 1.5); ax.set_ylim(-1.5, 1.5); ax.set_zlim(-1.5, 1.5) ax.set_title(f"Animated 3D Lissajous at Time Step {i}") # 实现平滑过渡显示 plt.pause(0.05) # 关闭交互状态恢复默认行为 plt.ioff(); plt.close() ``` 这里运用到了多次调用清除命令 (`clf`) 或局部清理方法 (`cla`) 结合暂停指令 (`pause`) 达成连续播放的效果[^2]。 --- ### 总结说明 以上分别介绍了两种主要形式——平面与立体下李萨茹图形的表现手法及其对应程序逻辑结构设计思路。无论是基础版还是高级定制化需求都能找到合适的解决方案。