Дифференциация с помощью sympy

Я пытаюсь дифференцировать функцию z. Но это дает мне ошибку, которая говорит: поднять TypeError («не удается преобразовать выражение в число с плавающей запятой»)

TypeError: невозможно преобразовать выражение в число с плавающей запятой

Я не могу понять свою ошибку.

from sympy import *
import numpy as np 
import math


T_wb =Symbol ('T_wb')
z =math.pow(10, 10.79574 * (1 - 273.16/T_wb)                    \          
    - 5.028 * math.log10(T_wb/273.16)                           \
            + 1.50475 * math.pow(10, -4)                        \
            * (1 - math.pow(10, -8.2969 * (T_wb / 273.16 - 1))) \
            + 0.42873 * math.pow(10, -3)                        \
            * (math.pow(10, 4.76955 * (1 - 273.16 / T_wb)) - 1) \
            + 2.78614)
zprime = z.diff(T_wb)
print (zprime)

python sympy differentiation
person Aditya Kodkani    schedule 04.09.2018    source источник

Ответы (2)


arrow_upward
2
arrow_downward

Математические функции модуля math отличаются от математических функций модуля sympy. Первые работают с числами (с плавающей запятой), а вторые работают с симпатичными выражениями и символами для выполнения аналитического вывода. Поэтому вместо math.log необходимо использовать sympy.log.

Функция pow отличается. Это встроенная функция Python (эквивалентная оператору **). Так же, как и другие операторы (+, -, *, /), нет необходимости вызывать специальную функцию. Например type(T_wb**2) дает хорошо sympy.core.power.Pow

import sympy as sp

T_wb = sp.Symbol('T_wb')

z =   pow(10, 10.79574 * (1 - 273.16/T_wb)) \
    - 5.028 * sp.log(T_wb/273.16, 10) \
    + 1.50475e-4 * (1 - pow(10, -8.2969 * (T_wb / 273.16 - 1))) \
    + 0.42873e-3 * (pow(10, 4.76955 * (1 - 273.16 / T_wb)) - 1) \
    + 2.78614

zprime = z.diff(T_wb)
print(zprime)

дает:

0.55857099968694*10**(4.76955 - 1302.850278/T_wb)*log(10)/T_wb**2 + 2948.9643384*10**(10.79574 - 2948.9643384/T_wb)*log(10)/T_wb**2 + 4.57049358434617e-6*10**(-0.0303737736125348*T_wb + 8.2969)*log(10) - 5.028/(T_wb*log(10))
person xdze2    schedule 04.09.2018

arrow_upward
0
arrow_downward

Его ошибка броска в math.log10(T_wb/273.16) T_wb является переменной, и когда он пытается вычислить журнал, sympy не может преобразовать его в число с плавающей запятой, потому что это переменная.

person Ankit Kumar Namdeo    schedule 04.09.2018