Python求解水文EI与DI双曲函数方程遇阻，求解决方案

首先，你的方程存在括号匹配错误，原方程右侧的平方根号部分括号不完整，正确的表达式应该是：
EI = sqrt(DI * tanh(1/DI) * (1 + sinh(DI) - cosh(DI)))

利用双曲函数恒等式cosh(x) - sinh(x) = e^{-x}，可以简化表达式，减少数值计算误差：
EI = sqrt(DI * tanh(1/DI) * (1 - e^{-DI}))

问题分析 #

1. 报错原因：sympy.nsolve需要指定初始猜测值，你没有提供，且该函数对非线性方程的数值稳定性不如专门的数值求解库；
2. 定义域限制：DI必须大于0，且EI的取值范围是(0,1)（当DI→0+时EI→0，DI→+∞时EI→1）。如果你的EI不在这个区间，方程无解。

解决方案 #

使用scipy.optimize.root_scalar（推荐brentq方法，无需初始猜测但需要根所在区间）来求解，步骤如下：

1. 先确认你的EI值在(0,1)范围内；
2. 定义目标函数（令f(DI)=0）；
3. 调用数值求解器在合理区间内找根。

修正后的代码 #

import numpy as np
from scipy.optimize import root_scalar

# 假设mean_ET和mean_perc是已定义的数组
mean_ET = [...]  # 你的实际蒸发数据
mean_perc = [...]  # 你的实际降水数据

EIlist = []
DIlist = []

# 定义目标函数：f(DI) = 计算得到的EI - 给定的EI
def target_func(DI, given_EI):
    term1 = DI * np.tanh(1/DI)
    term2 = 1 - np.exp(-DI)
    calc_EI = np.sqrt(term1 * term2)
    return calc_EI - given_EI

for ET, perc in zip(mean_ET, mean_perc):
    EI = ET / perc
    EIlist.append(EI)
    print(f"EI: {EI}")
    
    # 检查EI是否在有效范围
    if not (0 < EI < 1):
        print(f"EI={EI}超出有效范围(0,1)，无解")
        DIlist.append(None)
        continue
    
    # 使用brentq方法求解，区间设为(1e-6, 1000)覆盖所有可能的DI值
    sol = root_scalar(target_func, args=(EI,), method='brentq', bracket=[1e-6, 1000])
    
    if sol.converged:
        DIlist.append(sol.root)
        print(f"求解得到DI: {sol.root}")
    else:
        print(f"求解失败: {sol.flag}")
        DIlist.append(None)

关键说明 #

• brentq方法是基于区间的数值求解器，只要函数在区间两端符号相反（这里因为EI∈(0,1)，f(1e-6)≈0-EI<0，f(1000)≈1-EI>0），就能保证找到根；
• 简化后的表达式减少了双曲函数的计算量，同时降低数值误差；
• 增加了EI范围检查，避免无效求解。

内容的提问来源于stack exchange，提问作者Oskar kistrup