Python-сообщество

GJGUG · Фев. 27, 2021 00:42:23

Здравствуйте!

Подскажите пожалуйста, как реализовать в python интеграл, представленный на рисунке (или 18-19 строки в коде)?

Пытался через sympy.integrate, но потерпел крах - ValueError: Invalid limits given: (20.665,)

 import numpy as np
import sympy
from scipy import interpolate
 
b_2 = 6.67
l_k = 24.0
U_1 = [0.0, 4.3, 8.1, 11.4, 14.1, 16.3, 18.0, 19.3, 20.2, 20.7]
V_1 = [0.0, 0.6, 1.8, 3.5, 5.7, 8.3, 11.3, 14.8, 18.6, 22.8]
 
"Первая производная уравнения центроиды"
def Cn_1(x):
    x_points = U_1
    y_points = V_1
 
    tck = interpolate.splrep(x_points, y_points)
    return interpolate.splev(x, tck, der = 1)
 
"Расчет длины центроиды в зависимости от координаты X"
def l_cn(x):
        return sympy.integrate((1 + Cn_1(x)**2)**0.5, x)
 
"Полная длина центроиды"
L_max = l_cn(l_k-0.5*b_2)
print("L_max = ", L_max)

Прикреплённый файлы:
fig_3.PNG (43,4 KБ)

py.user.next · Фев. 27, 2021 03:29:41

На бумажке не пробовал его решить? Если под интегралом стоит число A, то ответ будет A * x + C, где C - константа.

Например:

integral(12345)dx = 12345 * integral(1)dx = 12345 * x + C

Отредактировано py.user.next (Фев. 27, 2021 03:30:08)

Python-сообщество

Уведомления

#1 Фев. 27, 2021 00:42:23

Неопределённый интеграл функции

#2 Фев. 27, 2021 03:29:41

Неопределённый интеграл функции

Board footer