Python - Как построить аргумент в интеграле, который не является интегрируемым значением

Question

Python - Как построить аргумент в интеграле, который не является интегрируемым значением

1

Я хочу интегрировать функцию, которая не имеет решения закрытой формы с неизвестной переменной, а затем построит зависимость от неизвестной переменной. Чтобы попробовать более простой тест, я попытался использовать интеграл от f(x,c) = (x^2+c), интегрированный по x и график с разными значениями c. Однако приведенный ниже код получает ошибку

только массивы size-1 могут быть преобразованы в сканы Python

хотя интеграл от числа, например, integral(5), кажется, возвращает правильное скалярное значение.

import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate


def f(x,c):
    return x**2+c

def integral(c):
    return integrate.quad(f,0,10, args = (c,))[0]
y = np.linspace(0,20,200)


plt.plot(y, integral(y))

Cor 15 окт. 2018, в 10:22

Источник

1

Вы действительно не должны писать функцию с именем int . Это создает ненужную путаницу со встроенной функцией .
khelwood 15 окт. 2018, в 07:56
0

Какая у вас переменная test ?
Sharku 15 окт. 2018, в 07:58
1

Я думаю, что test должен быть f так как это функция, определенная выше
onno 15 окт. 2018, в 08:04
0

отредактировано для изменения функции int на целое и исправлено опечатка теста
Cor 15 окт. 2018, в 15:54

Показать ещё 2 комментария

Теги:

python

scipy

plot

integral

2 ответа

Ещё вопросы

Вы действительно не должны писать функцию с именем int . Это создает ненужную путаницу со встроенной функцией .
Я думаю, что test должен быть f так как это функция, определенная выше
отредактировано для изменения функции int на целое и исправлено опечатка теста

onno · Answer 1 · 2018-10-15T06-42-00.000Z

Вы передаете массив numpy как аргумент c то время как вы хотите интегрировать по x для всех элементов c. Поэтому вы можете использовать это:

def f(x,c):
    return x**2+c

def integrate_f(c):
    result = np.zeros(len(c))
    counter = 0
    for item in c:
        result[counter] = integrate.quad(f,0,10, args = (item))[0]
        counter +=1
    return result

c_array = np.linspace(0,1,200)
plt.plot(c_array, integrate_f(c_array))

Sharku · Answer 2 · 2018-10-15T05-56-00.000Z

onno было немного быстрее. Но вот мое подобное решение. Вам нужно перебрать все разные c:

import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate


def f(x,c):
    return x**2+c

def getIntegral(c_list):
    result = []
    for c in c_list:
        integral = integrate.quad(f,0,10,args = c)[0]
        result.append(integral)
    return result

if __name__ == "__main__":


    c_list = np.linspace(0,20,200)

    plt.plot(c_list, getIntegral(c_list))
    plt.show()