odeint - (TypeError: невозможно преобразовать выражение в float) - импорт выражения в функцию для выполнения odeint - PullRequest
Новая
       114

odeint - (TypeError: невозможно преобразовать выражение в float) - импорт выражения в функцию для выполнения odeint

0 голосов
/ 01 мая 2020

Это ошибка, которую я получаю 

Traceback (most recent call last):

  File "C:\Users\user\.spyder-py3\Numerical Methods Problems\FreeFall.py", line 40, in <module>
    ans=odeint(vel,0,t)

  File "C:\ProgramData\Anaconda3\lib\site-packages\scipy\integrate\odepack.py", line 245, in odeint
    int(bool(tfirst)))

  File "C:\ProgramData\Anaconda3\lib\site-packages\sympy\core\expr.py", line 325, in __float__
    raise TypeError("can't convert expression to float")

TypeError: can't convert expression to float

А вот мой код - я довольно новичок в кодировании и учусь использовать его для численных расчетов: 

from scipy.integrate import odeint
from sympy import *
import numpy as np
import math

def diff_cd(re_no):
    Re=Symbol('Re')
    expr=(24/Re)+(6/(1+(Re**0.5)))+0.4
    ans=diff(expr,Re).subs(Re,re_no)    
    return ans

def diff_re(k,u_no):
    u=Symbol('u')
    expr=k*u
    ans=diff(expr,u).subs(u,u_no)
    return ans    

ans = [diff_cd(20),diff_re(11,15)]

rhog=1.2
mug=1.872*(10**(-5))
a=0.3

u=Symbol('u')
pi=math.pi

k=(2*rhog*a/mug)
Re=k*u

p1=(rhog*pi*(a**2)*u*Re*((24/Re)+(6/(1+(Re**0.5)))+0.4))+(0.5*rhog*pi*(a**2)*(u**2)*diff_cd(Re)*diff_re(k,u))

ansfu=p1*(-1/24)


def vel(y,t):
    dudt = ansfu
    return dudt

t=np.linspace(1,100,100)
ans=odeint(vel,0,t)


print(ans)

Мне просто нужно получить ответ без этой ошибки. Также есть ли способ сделать все это в одной функции

1 Ответ

1 голос
/ 02 мая 2020

Если я добавлю print(vel(0,0)) перед вызовом ode, я получу 

-271.868595022194*u**2*(-3.97723522060237e-7*u**(-0.5)/(u**0.5 +
0.00509901951359279)**2 - 1.6224e-8/u**2) - 
543.737190044387*u**2*(0.4 + 6/(196.116135138184*u**0.5 + 1) + 0.000624/u)

То есть выражение sympy, а не число. odeint не может работать с этим!

u определяется как Symbol, и поэтому любое выражение python, использующее его, также будет выражением sympy.

Особенно, если вы новичок в кодировании, вы должны придерживаться одного пакета или другого. Если вы определяете функцию символически, используйте sympy и ее собственные решатели. Но если важно числовое решение c, которое производит scipy, определите функцию с помощью python / numpy. Не пытайтесь смешивать sympy и numpy (без особого опыта).

...