Вопросы по теме 'numerical-integration'
Как выполнить численное интегрирование вектора в MATLAB?
У меня есть вектор из 358 чисел. Я хотел бы произвести численное интегрирование этого вектора, но я не знаю функции этого.
Я обнаружил, что мы можем использовать trapz или quad, но я не совсем понимаю, как интегрироваться без этой функции.
60594 просмотров
schedule
23.04.2023
Можно ли выполнить этот интеграл численно в Matlab или Mathematica?
Я хочу иметь возможность выполнить приведенный ниже интеграл полностью численно.
где , и , и являются константами, которые для простоты можно установить до 1 .
Интеграл по x можно выполнить аналитически вручную или с помощью...
1072 просмотров
schedule
23.04.2024
Численное интегрирование главного значения Коши
Я использую функцию интегрирования R, но подынтегральное выражение доставляет мне некоторые проблемы (у него есть особенность). Учти это,
distrib <- function(x){
dnorm(x, mean=700,sd=50)
}
phi <- function(nu, epsilon=20, maxi=Inf){...
1718 просмотров
schedule
03.02.2023
Пользовательская периодическая функция без счетчика
Я использую ode45 для решения простой ОДУ:
function dCdt=u_vent(t,C)
if t> 600 && t<= 720
Q=Q2;
elseif t> 1320 && t<= 1440
Q=Q2;
elseif t> 2040 && t<=...
575 просмотров
schedule
20.08.2023
Проблемы с интеграцией(), возвращающей интеграл 0
У меня возникли проблемы с интегралом, который, насколько я могу судить, должен вести себя хорошо, но неправильно возвращает 0 , когда я пытаюсь интегрировать от -Inf до Inf . Рассмотрим следующий код R:
t = 466
N = 50
df = 198
rscale =...
51 просмотров
schedule
05.09.2023
Ошибки в составном правиле Симпсона в C ++
Я пишу небольшую программу для аппроксимации интегралов на C ++ с использованием составного правила Симпсона. Кажется, это работает, но когда я вычисляю интеграл для разного количества интервалов (отсюда и цикл for ), появляются некоторые ошибки,...
1463 просмотров
schedule
21.02.2022
Изменить константу в вычислениях ОДУ при определенных условиях с помощью флага
У меня есть ODE для расчета изменения кислотности. Все работает нормально, только хотелось бы менять константу всякий раз, когда кислотность достигает критической точки. Предполагается, что это какой-то необратимый эффект, который я хочу...
706 просмотров
schedule
04.06.2022
Python: аппроксимация ln(x) с использованием ряда Тейлора
Я пытаюсь построить приближение для ln (1.9) с точностью до десяти цифр (так что 0,641853861).
Я использую простую функцию, которую я построил из ln[(1 + x)/(1 - x)]
Вот мой код:
# function for ln[(1 + x)/(1 - x)]
def taylor_two(r, n):...
3651 просмотров
schedule
16.08.2022
Вычисление трапециевидной интеграции от начала до каждой временной точки в R
У меня есть кадр данных x :
head(x)
# time Qfr
#1 1 0.004751271
#2 2 0.005405618
#3 3 0.005785781
#4 4 0.006028213
#5 5 0.006179973
#6 6 0.006263814
Я пытаюсь вычислить численное интегрирование от time = 0 до...
1544 просмотров
schedule
03.03.2023
2-кратный (повторный) интеграл в R
Я хочу вычислить двукратный повторяющийся интеграл (не двойной интеграл) в R, например,
куда
На практике и f(x), и g(x) довольно сложны, но для проведения эксперимента упростим g(x)=1 и f(x)=cos(x), в R я использую integrate для...
143 просмотров
schedule
15.11.2022
Решение двойного интеграла, связанного с многомерным нормальным числом
Я пытаюсь решить двойной интеграл, связанный с многомерной нормальной плотностью с известным средним вектором и ковариационной матрицей:
library(cubature)
mu1 <- matrix(c(3,3), nrow=2)
sigma1 <- rbind(c(4,-1), c(-1,6))
quadratic <-...
585 просмотров
schedule
19.11.2023
Нет значения в функции числового интегрирования Python
Я пытаюсь написать код, который вычисляет интегралы с использованием прямоугольного правила, а также позволяет пользователю вводить интегральные пределы и количество делений (прямоугольников). Я написал функцию, но для определенных значений она...
121 просмотров
schedule
24.08.2022
Распараллелить интегрировать в R
Мне было интересно, возможно ли и как распараллелить функцию интегрирования в R, чтобы ускорить численное интегрирование.
Я пишу здесь небольшой пример на случай, если кто-то может мне помочь.
b0=function(time){0.05*sin(0.1*time)+0.1}
d0=0.14...
312 просмотров
schedule
24.09.2022
Двойное интегрирование дискретных данных по времени
Допустим, у меня есть набор точек данных с именем signal , и я хочу дважды интегрировать его по времени (т. е. если signal было ускорением, я хотел бы интегрировать его дважды по времени, чтобы получить положение). Я могу интегрировать его один...
1087 просмотров
schedule
07.04.2024
уменьшение количества аргументов в функциях в C
Я хочу выполнить численное интегрирование функции f с помощью функции qtrap, определенной в «Численных рецептах в C».
double qtrap(double (*func)(double), double a, double b);
Как показано, это одномерное интегрирование переменной типа...
125 просмотров
schedule
25.04.2024
Численное интегрирование с использованием Fortran 90
Я пытаюсь использовать правило Симпсона для вычисления интеграла sqrt(1-x^2) в интервале от -1 до 1. Однако сумма, представленная переменной "s", в разработанном мной коде, не вообще не сходятся к пи более 2. Я абсолютный новичок в фортране и...
1032 просмотров
schedule
09.11.2023
Как интегрировать с помощью scipy массив как часть функции
Я хотел бы использовать интегральную команду в scipy и иметь функцию, которая один раз умножается на каждый элемент массива.
import math
import matplotlib.pyplot as plt
import numpy as np
import scipy.integrate as integrate
from scipy.integrate...
246 просмотров
schedule
24.10.2023
Каким будет вычислительно более быстрый способ реализовать это двумерное численное интегрирование?
Меня интересует двумерное численное интегрирование. Сейчас я использую scipy.integrate.dblquad , но он очень медленный. Пожалуйста, смотрите код ниже. Мне нужно оценить этот интеграл 100 раз с совершенно разными параметрами. Поэтому я хочу...
562 просмотров
schedule
16.02.2023
Реализация одномерной интеграции gsl в C++ с параметрами в качестве границы интеграции
Я хотел бы выполнить следующую одномерную интеграцию в С++ с использованием gsl,
I(x) = int_{x/4}^1 dy y^{3/2+a} (1-y)^{1/2} exp(1.13*sqrt(log(4y/x))), , где a и x рассматриваются как константы интегрирования, и я хотел бы иметь возможность...
134 просмотров
schedule
10.09.2023
Интеграция логистической дистрибуции в R
Функция плотности вероятности логистического распределения $f(x) = e^{-x} (1+e^{-x})^{-2}$
LL_cdf <- function(x) ((1+exp(-x))^(-1) # cumulative distribution function of logistic distribution
LL_pdf <- function(x) ( (exp(-x))*(1+exp(-x))^-2...
61 просмотров
schedule
20.07.2022