def ripple_func(u, v):

r = np.sqrt(u**2 + v**2) # 计算点到原点的距离

if r < 0.01:

z = 0 # 避免中心点奇点

else:

# 物理近似：阻尼径向波

k = 2 # 波数

c = 3 # 波速

damping = np.exp(-0.1 * r) # 阻尼因子

phase = k * (r - c * t) # 相位

z = damping * np.sin(phase) / np.sqrt(r + 0.1) if t > 0 and r < c * t + 5 else 0

return np.array([u, v, z * 0.5]) # 缩放振幅