Описание тега numerical-integration

Описание тега Вопросы с тегом
Алгоритмы, которые объединяют функции по одному или нескольким измерениям с использованием методов аппроксимации, вместо точных решений в замкнутой форме с использованием символической алгебры и исчисления. Включает такие концепции, как адаптивная квадратура, методы Монте-Карло, анализ конечных элементов, цепи Маркова.
1 ответ

Числовой интеграл от 0 до бесконечности

Моя цель состоит в том, чтобы вычислить числовой интеграл функции распределения вероятности (PDF) расстояния электрона от ядра атома водорода на языке программирования Си. Я написал пример кода, однако он не может правильно найти числовое значение и…
01 ноя '13 в 09:58
1 ответ

Интеграция Монте-Карло с использованием выборки по важности с учетом функции предложения

Учитывая предложение о распространении Laplace: g(x) = 1/2*e^(-|x|) и размер выборки n = 1000, Я хочу провести интеграцию Монте-Карло (MC) для оценки θ: через выборку важности. В конце концов я хочу вычислить среднее и стандартное отклонение этой оц…
05 дек '16 в 22:49
1 ответ

Интеграция числовых данных

Завтра я пишу тест и собираюсь сделать все для Matlab, чтобы сэкономить время. Некоторые вопросы требуют численного интегрирования точек данных (точек, не обязательно функций). Например C=[0 1 5 8 10 8 6 4 3 2.2 1.5 0.6 0]; Я использовал trapz(C) оп…
22 июн '15 в 17:58
1 ответ

Ошибка Matlab при вычислении двойного интеграла

Я пытаюсь вычислить двойной интеграл, но получаю ошибку. Мой код: clear Gamma = rand([22,1]); data_i = rand([1,10]); Y_1 = 1; fun = @(D_star, Y_0 ) fun1(D_star , Y_1 , Y_0 , Gamma , data_i); p = integral2 ( fun , 0 , 1 , 0 , 1); и функция fun1 опред…
28 ноя '17 в 02:32
1 ответ

Интегральная функция в R

Я определил следующую функцию в R (где a это заданная константа и это известная функция), но когда я пытаюсь вычислить ее значения, я получаю довольно странные числа. Я должен отметить, что для того, чтобы вычислить значения f(s)Я на самом деле реши…
12 окт '14 в 01:02
0 ответов

Использование числовой интеграции для поиска области, заключенной в петлю координат

Я пытаюсь найти или написать некоторый код, который может дать мне область между двумя строками, и если нижняя пересекает верхнюю, учитывается отрицательная область. Для этого (не могу опубликовать изображение:() . Я написал код, но не уверен на 100…
03 дек '18 в 15:48
1 ответ

Проблема с числовой интеграцией

В настоящее время я застрял на действительно простой ошибке, спрятанной где-то, и, надеюсь, кто-то сможет пролить свет на это. Я пытаюсь численно интегрировать набор данных с обычными методами, и получаю результаты, которые я не совсем ожидал. Затем…
27 июн '17 в 23:01
1 ответ

Решение двойного интеграла, связанного с многовариантным нормальным числом

Я пытаюсь решить двойной интеграл, связанный с многомерной нормальной плотностью с известным средним вектором и ковариационной матрицей: library(cubature) mu1 <- matrix(c(3,3), nrow=2) sigma1 <- rbind(c(4,-1), c(-1,6)) quadratic <- function…
08 ноя '16 в 00:02
4 ответа

Метод квадратурного численного интегрирования Тань-синь, сходящийся к неправильному значению

Я пытаюсь написать программу на Python, чтобы использовать квадратурную Tanh-sinh для вычисления значения: но хотя программа сходится к разумному значению без ошибок в каждом случае, она не сходится к правильному значению (которое является пи для да…
27 июл '14 в 23:14
1 ответ

Ищем конкретный алгоритм для численного интегрирования

Рассмотрим следующее дифференциальное уравнениеf(x) = g'(x) У меня есть код сборки, который выплевывает значения функции f(x) для переменной x, где x идет от 0 в very large, Теперь я ищу схему, которая будет анализировать эти значения f(x) чтобы опр…
15 июн '15 в 23:01
1 ответ

Остановить интеграцию через определенное время в Matlab

Я хочу прекратить решение дифференциального уравнения в Matlab, если это займет больше времени. Я попробовал следующее, но это не работает... options = odeset('AbsTol',1e-8,'RelTol',1e-5); RUNTIME=5; timerID = tic; while (toc(timerID) < RUNTIME) …
30 авг '14 в 11:33
1 ответ

Построение функции путем определения X и Y, а затем интеграция в R

Мне нужно построить функцию со значениями x, поступающими из первого столбца этой матрицы ниже, и значениями y, поступающими из второго столбца из той же матрицы, с целью последующего вычисления интеграла в требуемом диапазоне. matrix=structure(c(0.…
26 окт '16 в 23:51
1 ответ

Численное интегрирование свертки интерполирующей функции

У меня есть x и y массивы точек данных, которые я использовал для создания интерполяционной функции func_spline как показано ниже. import numpy from numpy import loadtxt from scipy.interpolate import * x_given = numpy.arange(1,21400,21400/25000) y_g…
25 янв '18 в 13:22
1 ответ

Вызов функции в качестве переменной в другой функции

Я немного новичок в c-программировании, я пытаюсь интегрировать функцию с правилом трапеции. Моя функция заключается в следующем: extern double integrate_trapezoid(double function,double min, double max, double subdivisions) { double h = (max-min)/s…
07 ноя '13 в 09:36
1 ответ

Флаг итерации в Matlab ODE Solvers

Я решаю ODE, как y''(t) + a y(t) + b = 0 с матлабом ode45 решатель. Он повторяется до тех пор, пока не будет найдено решение, которое удовлетворяет критериям ошибки. Я хочу прочитать предыдущее значение y на каждом этапе для сравнения. Каков был бы …
06 авг '13 в 12:43
1 ответ

Ошибка численного интегрирования

Когда я запускаю следующий скрипт ниже, я получаю сообщение об ошибке Too many input arguments (Я понял, что это конкретно происходит при расчете intNH(5,2) во-вторых for петля). В основном происходит то, что fun2 за i=5 равняется 0 (потому что соот…
29 июн '16 в 13:46
3 ответа

Объем 3d-формы с использованием численного интегрирования со Scipy

Я написал функцию для вычисления объема пересечения куба и полупространства, и сейчас я пишу тесты для него. Я попытался вычислить объем численно следующим образом: integral = scipy.integrate.tplquad(lambda z, y, x: int(Vector(x, y, z).dot(normal) &…
17 мар '18 в 15:12
1 ответ

Может ли этот интеграл быть сделан численно в Matlab или Mathematica?

Я хочу быть в состоянии сделать интеграл ниже численно. где , а также , а также константы, которые для простоты могут быть установлены в 1, Интеграл по x можно сделать аналитически вручную или с помощью Mathematica, а затем y можно сделать численно …
17 авг '11 в 23:21
2 ответа

Как векторизовать функцию для интеграла2?

Я хочу оценить двойной интеграл вида $$\int_{-\infty}^a \int_{-\infty}^b \sum_{i,j}^K a_ia_jx^iy^j\exp(-x^2 - y^2 + xy)dx dy $$ где $ a_i $ и $ a_j $ являются константами. Поскольку интеграл линейный, я могу поменять местами суммирование и интегриро…
21 мар '13 в 19:08
1 ответ

Численное интегрирование в R: ошибка поведения плохого подынтегрального выражения

Это моя матрица func=structure(c(-14.7690673280818, -14.5543581356252, -12.1406211639974, -10.7200919648493, -9.55507848352318, -9.20790894914246, -8.74647670464071, -8.26548763467919, -7.3962484443768, -6.94590909664862, -6.63903257406218, -6.40598…
02 ноя '16 в 20:46