如何相加f(x,n)与f(x,m)并积分二者乘积？代码报错求助

• 函数相加：对于两个单变量函数 ( f(x,n) ) 和 ( f(x,m) )，它们的和函数通过逐点求和定义：( h(x) = f(x,n) + f(x,m) )，即对任意输入 ( x )，取两个函数在该点的函数值相加即可。
• 乘积的积分：
  1. 不定积分：( \int f(x,n) \cdot f(x,m) dx )，可根据分部积分、换元法等规则结合函数具体形式求解；
  2. 定积分：( \int_{a}^{b} f(x,n) \cdot f(x,m) dx )，在不定积分结果基础上代入上下限计算，或直接用数值积分方法（若函数无解析解）。

报错原因 #

混用了NumPy和SymPy的函数：NumPy的np.cos仅支持数值数组，无法处理SymPy的符号变量x，导致抛出AttributeError。即使替换为sp.cos，嵌套循环的逻辑错误仍会导致积分失败。

核心修正点 #

1. 全流程使用SymPy符号运算：常量（如π）和函数（如cos）都改用SymPy内置对象，彻底避免与NumPy混用。
2. 调整循环逻辑：嵌套循环会让同一个n对应的项被重复添加4次（对应m的4次循环），这不符合求和逻辑。应分开两个独立循环，分别生成n和m的所有项后再求和。

修正后的代码 #

import sympy as sp

x = sp.symbols('x')

dp_xi = []
dp_eta = []

# 生成n从1到4的所有xi项
for n in range(1, 5):
    coeff = (n - 2) * sp.pi / 180
    term = sp.cos(coeff * x) * coeff
    dp_xi.append(term)

# 生成m从1到4的所有eta项
for m in range(1, 5):
    coeff = (m - 2) * sp.pi / 180
    term = sp.cos(coeff * x) * coeff
    dp_eta.append(term)

# 对所有项求和
df_xi = sum(dp_xi)
df_eta = sum(dp_eta)

# 计算积分并化简结果
I = sp.integrate(df_xi * df_eta, (x, 0, 1))
print("I=", sp.simplify(I))

额外说明 #

• 若需将积分结果转为数值，可调用I.evalf()方法获取浮点值；
• SymPy的simplify()方法可简化复杂的解析积分结果，提升可读性。

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