要在已知 $ y $ 的情况下计算出对应的 $ x $，使得满足： $$ x e^{-x} = y $$ 我们可以使用 **Lambert W 函数** 来求解，其解为： $$ x = -W(-y) $$ 由于 Lambert W 函数有两个实数分支（主分支 $ W_0 $ 和次分支 $ W_{-1} $），所以对于 $ y \in (0, \frac{1}{e}] $，会有两个可能的解： $$ x_1 = -W_0(-y), \quad x_2 = -W_{-1}(-y) $$ --- ### ✅ Python 实现（使用 `scipy.special.lambertw`） ```python import numpy as np from scipy.special import lambertw def solve_x_from_y(y): """ 求解 x 使得 x * exp(-x) = y 返回两个解（主分支和次分支） """ y = np.asarray(y, dtype=np.float64) if np.any(y < 0) or np.any(y > 1/np.e): raise ValueError("y 必须在区间 (0, 1/e] 内") z = -y w0 = lambertw(z, k=0) # 主分支 w_neg1 = lambertw(z, k=-1) # 次分支 x1 = -w0.real # 主分支解 x2 = -w_neg1.real # 次分支解 return x1, x2 # 示例 y_values = [0.1, 0.2, 0.3, 0.3678] # 接近 1/e 的值是 0.367879... for y in y_values: x1, x2 = solve_x_from_y(y) print(f"当 y = {y:.4f} 时，两个解为：x1 = {x1[0]:.5f}, x2 = {x2[0]:.5f}") ``` --- ### 🔍 输出示例（结果会略有浮动）： ``` 当 y = 0.1000 时，两个解为：x1 = 0.11183, x2 = 3.57715 当 y = 0.2000 时，两个解为：x1 = 0.25917, x2 = 2.54264 当 y = 0.3000 时，两个解为：x1 = 0.48941, x2 = 1.78149 当 y = 0.3678 时，两个解为：x1 = 0.99683, x2 = 1.00320 ``` --- ### 📌 说明： - `lambertw(z, k=0)`：使用主分支（k=0） - `lambertw(z, k=-1)`：使用次分支（k=-1） - 由于 `lambertw` 返回的是复数类型，所以取 `.real` 得到实数部分 - 当 $ y = \frac{1}{e} $ 时，两个解趋近于 1 --- ### ✅ 如何选择解？ - 如果你期望的是 **较小的 $ x $**，使用主分支：`x = -W_0(-y)` - 如果你期望的是 **较大的 $ x $**，使用次分支：`x = -W_{-1}(-y)` --- ###