烟花秀是一个令人惊叹的视觉效果，它能够为观众带来无尽的欢乐和惊喜。在Python中，我们可以使用matplotlib和numpy库来模拟和绘制出类似的效果。下面是一个简单的Python代码示例，演示如何实现一个基本的烟花模拟。
首先，确保你已经安装了matplotlib和numpy库。你可以使用以下命令来安装它们：

pip install matplotlib numpy

接下来，创建一个Python脚本（例如fireworks.py），并将以下代码粘贴到脚本中：
```python
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

设置模拟的时间步数和总粒子数

num_steps = 100
num_particles = 50

创建模拟的空间范围

x_min, x_max = -10, 10
y_min, y_max = -10, 10
z_min, z_max = -10, 10
x = np.linspace(x_min, x_max, num_particles)
y = np.linspace(y_min, y_max, num_particles)
z = np.linspace(z_min, z_max, num_particles)
X, Y, Z = np.meshgrid(x, y, z)

初始化粒子位置和速度数组

positions = np.zeros((num_particles, 3))
velocities = np.zeros((num_particles, 3))

随机初始化粒子位置和速度

for i in range(num_particles):
positions[i, 0] = np.random.uniform(x_min, x_max)
positions[i, 1] = np.random.uniform(y_min, y_max)
positions[i, 2] = np.random.uniform(z_min, z_max)
velocities[i, 0] = np.random.uniform(-1, 1)
velocities[i, 1] = np.random.uniform(-1, 1)
velocities[i, 2] = np.random.uniform(-1, 1)

设置模拟的时间步长和重力加速度

dt = 0.1
g = 9.8

进行模拟的循环

for step in range(num_steps):

计算每个粒子的加速度和速度变化量

accelerations = np.zeros((num_particles, 3))
for i in range(num_particles):
for j in range(num_particles):
distance = np.sqrt((positions[i, 0] - positions[j, 0])2 + (positions[i, 1] - positions[j, 1])2 + (positions[i, 2] - positions[j, 2])2)
accelerations[i, 0] += (positions[j, 0] - positions[i, 0]) / distance3 velocities[j, 0] dt
accelerations[i, 1] += (positions[j, 1] - positions[i, 1]) / distance3 velocities[j, 1] dt
accelerations[i, 2] += (positions[j, 2] - positions[i, 2]) / distance3 velocities[j, 2] dt
accelerations[i] -= g np.sign(velocities[i]) dt # 加速度减去重力加速度的影响
velocities[i] += accelerations[i] dt # 根据加速度更新速度
positions[i] += velocities[i] dt # 根据速度更新位置
if positions[i][2] < z_min: # 如果粒子触及地面，将其速度设为0并向上移动一定距离以重新开始动画效果
velocities[i][2] = 0 # 重置速度的z分量
positions[i][2] = z_min + np.random.uniform(0, 0.5)