Question

我已经定义了一个函数integrate_boole，如下所示：

def integrate_boole(f,l,r,N):
N = 4 * int(N // 4)
h=((r-l)/N)
xN = np.linspace(l,r,N+1)
fN = f(xN)
return ((2*h)/45)*(7*fN[0]+32*(np.sum(fN[1:-1:2]))+12*. (np.sum(fN[2:-2:4]))+14*(np.sum(fN[4:-4]))+7*fN[-1])

我现在正在尝试为f（x）在0和pi之间进行积分，其中f（x）定义为：

def f(x):
return x*np.sin(x)/(1+np.cos(x)**2)

当我用f（x）代替函数时：

integrate_boole(f(x),0,np.pi,8)

我收到以下错误消息：

-----------------------------------------------------------------------  ----
AttributeError                            Traceback (most recent call    last)
<ipython-input-23-91f1a95793b5> in <module>()
----> 1 integrate_boole(f(x),0,np.pi,8)

 <ipython-input-20-165d275ae26c> in f(x)
  1 def f(x):
----> 2     return x*np.sin(x)/(1+np.cos(x)**2)
  3 myx = np.linspace(0,np.pi,1000)
  4 plt.plot(myx,f(myx),label=r"$f(x)=\frac{x\sin{x}}{1+\cos^2{x}}$")
  5 plt.ylabel("$f(x)$")

AttributeError: 'Symbol' object has no attribute 'sin'

我最初在同一笔记本的开头导入了以下库：

import numpy as np
import matplotlib.pyplot as plt
import sympy as sym
from ipywidgets.widgets import interact
sym.init_printing(use_latex="mathjax")
x, y, z, t = sym.symbols('x y z t')

我已经检查了其他有关大致相同主题的文章，但是我不相信我的numpy和sympy库会像它们为其他人所做的那样发生冲突。但是它们可能是，我不知道为什么。任何帮助将不胜感激。

Answer 1

integrate_boole的第一个参数应该是一个函数。

f是一个函数
f(x)是函数f为输入x返回的数字。

您需要致电

integrate_boole(f,0,np.pi,8)

不是

integrate_boole(f(x),0,np.pi,8)

这是完整的代码，带有正确的缩进：

import numpy as np
import matplotlib.pyplot as plt
import sympy as sym

x, y, z, t = sym.symbols('x y z t')

def integrate_boole(f,l,r,N):
    N = 4 * int(N // 4)
    h=((r-l)/N)
    xN = np.linspace(l,r,N+1)
    fN = f(xN)
    return ((2*h)/45)*(7*fN[0]+32*(np.sum(fN[1:-1:2]))+12*(np.sum(fN[2:-2:4]))+14*(np.sum(fN[4:-4]))+7*fN[-1])

def f(x):
    return x*np.sin(x)/(1+np.cos(x)**2)

integrate_boole(f,0,np.pi,8)
# 2.470207745145361

'Symbol'对象没有属性'sin'

1 个答案: