scipy.integrate.quad内部被积函数返回NoneType异常排查与解决

需求背景 #

需要计算以下积分：
$$A_0 = \alpha - \frac{1}{\pi} \int_0 ^ c \left[ \frac{2 z'}{c} \cdot \frac{1}{\sqrt{1 - \left( \frac{2x}{c} - 1 \right) ^ 2}} \right] dx$$
其中$z'$是分段函数：
$$ z' = \begin{cases}
\frac{2m}{p^2} (p - x), & 0 \leq x \leq p\
\frac{2m}{(1 - p)^2} (p - x), & p < x \leq 1
\end{cases}$$

初始代码与第一个错误 #

初始实现代码 #

import numpy as np
import scipy.integrate as integrate

alpha = np.linspace(-10, 10, 100) # deg
NACA = '2412'
c = 1.2 # m
   
def z_prime(x, naca, c):
   m = float(naca[0]) / 100
   p = float(naca[1]) * c / 10

   if 0 <= x <= p:
      return (2 * m) * (p - x) / (p ** 2)
   if p < x <= 1:
      return (2 * m) * (p - x) / ((1 - p) ** 2)

def integrand_0(x, naca, c):
   return 2 * z_prime(x, naca, c)/ c / np.sqrt(1 - (2 * x / c - 1) ** 2)

for a in alpha:
   A0 = np.deg2rad(a) - integrate.quad(integrand_0, 0, c, args = (NACA, c)) / np.pi

错误信息 #

return 2 * z_prime(x, naca, c)/ c / np.sqrt(1 - (2 * x / c - 1) ** 2)
TypeError: unsupported operand type(s) for *: 'int' and 'NoneType'

错误原因 #

z_prime函数的判断条件存在范围漏洞：积分上限是c=1.2，当x>1时，两个if条件都不满足，函数无返回值，默认返回None。后续计算2 * z_prime(...)时，就会出现int与NoneType相乘的类型错误。

修改代码后出现的第二个错误 #

修改后的z_prime判断逻辑 #

if 0 <= x / c <= p:
   return (2 * m) * (p - x) / (p ** 2)
if p < x / c <= 1:
   return (2 * m) * (p - x) / ((1 - p) ** 2)

新错误信息 #

A0 = np.deg2rad(a) - (integrate.quad(integrand_0, 0, c, args = (NACA, c))) / np.pi
TypeError: unsupported operand type(s) for /: 'tuple' and 'float'

错误原因 #

scipy.integrate.quad函数返回的是元组，包含两个元素：(积分结果, 误差估计值)，不能直接对元组执行除法运算。

最终解决方法 #

提取quad返回的第一个元素（即积分结果）进行计算：

for a in alpha:
   integral_result = integrate.quad(integrand_0, 0, c, args=(NACA, c))[0]
   A0 = np.deg2rad(a) - integral_result / np.pi

另外建议给z_prime函数添加默认返回值，避免后续出现类似None的问题，比如在函数末尾添加：

else:
   return 0.0  # 或根据实际需求返回合理值

内容的提问来源于stack exchange，提问作者Jacob Ivanov