Публикации по теме 'numerical-integration'


Численное интегрирование: интегрирование Ромберга
Интегрирование с использованием численных методов Я ценю эту концепцию интеграции прежде всего потому, что среди методов, которые мы обсуждали в классе, это единственный метод, который рассматривает возможность устранения своих ошибок. Хоть и не все, но устранение очень хорошо помогает. Прежде всего, давайте подойдем к этому, используя также подход сверху вниз. Сначала мы обсуждаем численное интегрирование в целом, а затем переходим к интегрированию Ромберга. Численное..

Вопросы по теме 'numerical-integration'

как использовать scipy.integrate, чтобы получить объем усеченной сферы?
Я борюсь с использованием scipy.integrate, я использовал tplquad, но как я могу использовать integrate , чтобы получить объем (усеченной) сферы? Спасибо import scipy from scipy.integrate import quad, dblquad, tplquad from math import* from numpy...
3118 просмотров
schedule 12.06.2023

Точность численного интегрирования в Matlab
Я пытаюсь интегрировать аналитическую функцию (составную функцию sqrt и trig) в области прямоугольника. У него нет точки сингулярности в этом районе, и он кажется идеальным кандидатом на использование dblquad . Мой вопрос заключается в том, как...
880 просмотров

Несколько выходов решателя Matlab ODE
У меня есть следующий код Matlab ODE: [t,y,~,~,ie] = ode23tb(@(t,y) RHSODE(t,y),[0,t_end], [i0;v0],options); Я хочу, чтобы решатель ODE также мог дать мне результат z, который является функцией y и dy / dt, такой, что z = f (y, dy / dt)....
1354 просмотров
schedule 15.06.2023

Matlab: вывод ode45 неверен для демпфера принудительной пружины
Вот мой код Matlab для решения ODE второго порядка для системы масса-пружина-приборная панель: function Spring clear all; close all; options=odeset('RelTol',1e-6); p0 = [1 0]; %initial position and velocity [t,p] = ode45(@SpringFunction, [0 20],...
1775 просмотров
schedule 24.04.2022

Метод квадратурного численного интегрирования Tanh-sinh сходится к неправильному значению
Я пытаюсь написать программу Python, чтобы использовать квадратуру Tanh-sinh для вычисления значения: но хотя программа сходится к разумному значению без ошибок в каждом случае, она не сходится к правильному значению (которое равно пи для...
1725 просмотров

Почему нецелые показатели степени вызывают появление nan в Python?
Я пытаюсь численно решить уравнение Лейна-Эмдена в Python, используя scipy.integrate.ode . По какой-то причине мой код работает с целыми значениями n (индекс политропы), такими как 3, но не работает с нецелыми значениями, такими как 2,9 и 3,1....
823 просмотров

Путаница на Фортране, модули, объявления, фиктивные аргументы
моя проблема состоит в том, чтобы интегрировать функцию с помощью интеграции romberg. Я пересдаю этот курс и успешно закончил последний семестр, но на этот раз мой профессор требует модули, которые я никак не мог понять. Другой проблемой всегда были...
234 просмотров

Интеграл внутри экспоненты внутри интеграла
Я пытаюсь численно интегрировать несколько разных выражений, но все они имеют интеграл внутри экспоненты, окруженный другим интегралом, т.е. Кажется, я не могу решить эту проблему, используя приведенный ниже код: syms s p t W =...
188 просмотров
schedule 10.10.2022

Построение функции путем определения X и Y, а затем интегрирования в R
Мне нужно построить функцию со значениями x, исходящими из первого столбца этой матрицы ниже, и значениями y, исходящими из второго столбца из той же матрицы, с целью последующего вычисления интеграла в желаемом диапазоне.:...
290 просмотров

Численный расчет тройного интеграла с использованием Python
Мне очень нужно вычислить этот интеграл. Я пытался сделать это в течение нескольких месяцев, используя пакет numpy в Python и, в частности, функциюintegr.tplquad. from __future__ import division from math import * import numpy as np import...
1332 просмотров

Неправильный результат при выполнении простой интеграции Монте-Карло в R
Я делаю часть презентации по численному интегрированию. В то время как сам разговор пойдет о более совершенных формах численного интегрирования (в основном выборка по важности и стратифицированная выборка), я упоминаю в части своего раздела выборку...
457 просмотров

Числовая интеграция в numpy
Я хочу сделать что-то совсем простое, но не могу найти в глубинах numpy . Я хочу численно и непрерывно интегрировать функцию, заданную ее значениями ( не ее формулой!). Это означает, что мне просто нужен массив, содержащий суммы начала входного...
14159 просмотров

Используйте интеграл2 или интеграл3, чтобы вычислить CDF в MATLAB
Я только что столкнулся с проблемой, используя integral2 или integral3 для расчета CDF в MATLAB. Предположим, у меня есть две независимые нормальные случайные величины X и Y, вектор среднего значения равен mu = [5;50] , а ковариационная матрица...
77 просмотров

Уравнение Шредингера не эволюционирует должным образом со временем?
Я пишу код на питоне для развития зависящего от времени уравнения Шредингера с использованием схемы Крэнка-Николсона. Я не знал, как справиться с этим потенциалом, поэтому огляделся и нашел путь от -crank-nicolson-method">этот вопрос , который я...
333 просмотров

Рекурсия в CUDA возвращает неправильный доступ к памяти
Я пишу программу численного интегрирования, которая реализует правило трапеций с адаптивным размером шага. Не вдаваясь в подробности, алгоритм использует рекурсию для вычисления интеграла закодированной математической функции в заданном интервале с...
391 просмотров
schedule 27.05.2023

Интегрирование по постоянной функции
Я пытаюсь интегрировать постоянную функцию в MATLAB 2017a, но я застрял. Прежде всего, когда я интегрирую с помощью следующего скрипта, я получаю правильный результат. Итак, скрипт работает для x0 , который зависит от t . function E=sol(n,k)...
194 просмотров

Вложенный интеграл в Python
У меня есть M -мерный интеграл, где внешние пределы над x_M равны [0, y] , следующие пределы над x_{M-1} равны [0, max(x_M, y-x_M)] , .... а внутренний интеграл превышает x_1 с пределами [0, max(x_2, y-x_M-...-x_2)] . Функция /...
259 просмотров

Поиск областей полигонов: интегрирование при вертикальных и горизонтальных геометрических ограничениях
Я пытаюсь рассчитать области выше нуля и ниже кривой. Моя кривая имеет дискретные значения x и y , которые выглядят как в примере ниже. y <- c(-1, 5, 2, 3.5, 1, 5.5, -2, 3, -1) plot(1:length(y), y, type = "l") abline(h = 0) Я пытаюсь...
228 просмотров

Как интегрировать (AUC) модель nls и доверительный интервал Монте-Карло в R
Я пытаюсь интегрировать (разрешить область под) нелинейную функцию (от nls() ) от x = 0 до бесконечности в R. Однако функция интеграции R вызывает функцию ( f ). Короче, хотелось бы сделать что-то приближенное: integrate(my.nls, lower = 0L,...
235 просмотров

Как решить уравнение пучка Эйлера – Бернулли в FiPy?
Чтобы понять, как работает FiPy, я хочу решить уравнение балки Эйлера-Бернулли. с фиксированными конечными точками: w''''(x) = q(x,t), w(0) = w(1) = 0, w'(0) = w'(1) = 0. Для простоты пусть q(x,t) = sin(x) . Как я могу определить и...
161 просмотров