TyperError: неподдерживаемые типы операндов для -: 'float' и 'builtin_function_or_method' - PullRequest
Новая
       84

TyperError: неподдерживаемые типы операндов для -: 'float' и 'builtin_function_or_method'

1 голос
/ 25 апреля 2020

Я пытаюсь реализовать алгоритм Булирша-Стоера: я создал функцию, которая содержит его упрощенный алгоритм. Но сначала я создал функцию midpoint: в основном, первое применение этого метода средней точки использует два этапа интеграции. Число шагов увеличивается на два при последовательной интеграции, за каждой интеграцией следует экстраполяция Ричардсона: процедура останавливается, когда два последовательных решения отличаются (в смысле root -значного квадрата) менее чем на заданный допуск. 

def integrate(F,x,y,xStop,tol):
def midpoint(F,x,y,xStop,nSteps):
#Midpoint formulas    
    h = (xStop -x)/nSteps
    y0 = y
    y1 = y0 + h*F(x,y0)

    for i in range(nSteps-1):
        x = x + h
        y2 = y0 + 2.0*h*F(x,y1)
        y0 = y1
        y1 = y2
    return 0.5*(y1 + y0 + h*F(x,y2))
def richardson(r,k):
#Richardson's extrapolation    
    for j in range(k-1,0,-1):
        const = (k/(k - 1.0))**(2.0*(k-j))
        r[j] = (const*r[j+1] - r[j])/(const - 1.0)
    return

kMax = 51
n = len(y)
r = np.zeros((kMax,n))
#Start with two integration steps
nSteps = 2
r[1] = midpoint(F,x,y,xStop,nSteps)
r_old = r[1].copy
#Increase the number of integration by 2
#and refine result by Richardson extrapolation
for k in range(2,kMax):
    nSteps = 2*k
    r[k] = midpoint(F,x,y,xStop,nSteps)
    richardson(r,k)
    #Compute RMS change in solution
    e = sqrt(np.sum((r[1] - r_old)**2)/n)
    #Check for convergence
    if e < tol: return r[1]
    r_old = r[1].copy()
print("Midpoint method did not converge")

И вот алгоритм: 

def bulStoer(F,x,y,xStop,H,tol = 1.0e-6):
X =[]
Y =[]
X.append(x)
Y.append(y)
while x < xStop:
    H = min(H,xStop -x)
    y = integrate(F,x,y, x + H, tol) #Midpoint method
    x = x + H
    X.append(x)
    Y.append(y)
return np.array(X),np.array(Y)

Я пытаюсь проверить это на примере: 

def F(x,y):
F =np.zeros(2)
F[0] = y[1]
F[1] =(-y[1] - y[0]/0.45 + 9.0)/2.0
return F

H = 0.25
xStop = 10.0
x = 0.0
y = np.array([0.0,0.0])
X,Y = bulStoer(F,x,y,xStop,H)
plt.plot(X,Y[:,1],'o-')
plt.xlabel('Time(s)')
pkt.ylabel('Current(A)')
plt.grid(True)
plt.show
input("\nPress return to exit")

, но у меня есть такая ошибка: 

unsupported operand type(s) for -: 'float' and 'builtin_function_or_method'

В этой строке: e = sqrt(np.sum((r[1] - r_old)**2)/n)

Может ли кто-нибудь помочь мне с этой ошибкой: Я вообще не понимаю ее происхождение !!

Спасибо много!

1 Ответ

0 голосов
/ 25 апреля 2020

изменить строку r_old = r[1].copy на r_old = r[1].copy()

...