Вопросы по теме 'quad'

Как определить последовательность точек разрыва для метода integte.quad
У меня есть подынтегральное выражение, подобное дельта-функции Дирака. Но он конечно и с очень узкой базой. Мне нужно интегрировать это, но метод quad не распознал эту особенность. Посмотрев в документы python: точки : (последовательность...
859 просмотров
schedule 10.10.2021

Численное интегрирование сильно колеблющегося 1-D подынтегрального выражения (содержащего функции Бесселя) в Python
Я пытаюсь численно оценить подынтегральное выражение с действительным знаком, состоящее из нескольких функций Бесселя (первого и второго рода). Подынтегральное выражение колеблется и затухает, и его необходимо оценивать в диапазоне от 0 до + ∞. До...
1536 просмотров

Распараллелить вычисление матричных элементов с помощью scipy.integrate.quad
Мне нужно оценить каждый элемент матрицы, используя функцию с числовым интегралом (scipy.integrate.quad). Элементами матрицы являются пиксели серого изображения размером 5202 x 3465. У меня есть доступ к графическому процессору, и я хотел бы...
181 просмотров

Использование scipy.quad с трюком iε: плохие результаты
Чтобы обойти значение принципа Коши, я попытался интегрировать интеграл, используя небольшой сдвиг iε в комплексную плоскость, чтобы избежать полюса. Однако, как видно из рисунка ниже, результат довольно плохой. Код этого результата показан ниже. У...
186 просмотров

Определенный интеграл по одной переменной в функции с двумя переменными в Scipy
Я пытаюсь вычислить определенный интеграл функции с несколькими переменными только по одной переменной в scipy. Это похоже на то, как выглядит мой код- from scipy.integrate import quad import numpy as np def integrand(x,y): return...
8012 просмотров
schedule 06.08.2022

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

Проблема с интеграцией scipy quad в python
Я новичок в python, поэтому надеюсь, что моя формулировка имеет смысл. В настоящее время я пытаюсь смоделировать набор уравнений, которые требуют, чтобы произведение интегрирования умножалось на число с плавающей запятой. Я получаю вывод Nan только...
1582 просмотров
schedule 08.06.2023

Python: оценка интеграла для массива
import numpy import matplotlib.pyplot as plt from scipy import integrate def f(x,y): return x*y + x**2 def integral(x,y): I = integrate.quad(f, 0, x, args=(y,))[0] return I def gau(x,y): return (1+x)*integral(x,y) xlist =...
520 просмотров
schedule 16.04.2023

Есть ли способ использовать scipy.optimize.fsolve с jit_integrand_function и scipy.integrate.quad?
На основе приведенного здесь объяснения 1 , я пытаюсь использовать ту же идею для ускорения следующего интеграла: import scipy.integrate as si from scipy.optimize import root, fsolve import numba from numba import cfunc from numba.types import...
252 просмотров
schedule 10.07.2023

Есть ли способ численно интегрировать произведение двух функций в Python?
У меня есть две функции, которые принимают несколько аргументов: import numpy as np from scipy.integrate import quad gamma_s=0.1 #eV gamma_d=0.1 #eV T=298 #K homo=-5.5 #eV Ef=-5 #eV mu=0 #eV just displaces the function #Fermi-Dirac distribution...
60 просмотров
schedule 29.03.2023