Поэтому я закодировал эту формулу введите описание изображения здесь

Который я получаю:

def sumAN(theta,CoefAN,n_cl):
    # this function give us the sumatory in the right side of the formula
    Sumatorio = np.array([])
    for count,i in enumerate(theta):
        sumatorio = 0
        for count2,j in enumerate(n_cl):
            sumatorio = sumatorio +CoefAN[count2]*sin(int(count2+1)*i)
        Sumatorio = np.append(Sumatorio,sumatorio)
    return Sumatorio

cl= 4*((np.radians(alpha)+A0)*tan(theta/2)+sumAN(theta,CoefAN,n_cl))

Чтобы немного объяснить это:
- Альфа: постоянная
- A0: постоянный
- AN: np.array ([]) (n значений)
- theta: независимая переменная
После этого мне нужно вычислить следующий интеграл: введите описание изображения здесь

Вот где у меня проблемы:

ch = integrate.quad(lambda theta:(4*((alpha_char+A0)*tan(theta/2)+sumAN(theta,CoefAN,n_charl)))*(cos(theta)-cos(xa))*sin(theta),0,xa)[0]

У меня есть все ограничения и все такое. Но получаю следующую ошибку:

объект с плавающей точкой не повторяется

Я не знаю, как продолжить. Итак, мой вопрос: как я могу интегрировать эту функцию с помощью метода integration.quad? Может, я поменяю способ составления суматориев? Как я могу написать функцию иначе, чтобы это работало? заранее спасибо

0
jaimebw 15 Апр 2020 в 19:31

1 ответ

Лучший ответ

Это должно работать

import numpy as np
from scipy.integrate import quad

def integrand(theta, theta_a, alpha, A):
    sum = 0
    # get sum
    for index, value in enumerate(A):
        if index == 0:
            sum += (alpha + A[index]) * np.tan(0.5 * theta)
        else:
            sum += A[index] * np.sin(index * theta)
    # note that multiplication with 4 and multiplication with 1/4
    # result in one as prefactor
    return -sum * (np.cos(theta) - np.cos(theta_a))

# calculate integral 
theta_a = 0
alpha = 0 
array_coefficients = np.array([1, 2, 3, 4])
integral = quad(integrand, 0, 1, args=(theta_a , alpha, array_coefficients))
1
MachineLearner 15 Апр 2020 в 16:43