Публикации по теме 'numerical-methods'
Числовая интерполяция: естественный кубический сплайн
Изучение возможностей кусочной интерполяции
Я был очарован этим методом, когда впервые услышал о нем, особенно о его происхождении. Этот пост покажет вам, как я решил это и сделал проще. Мы будем использовать нисходящий подход и стараемся, чтобы вы визуализировали во время чтения, чтобы лучше понять это.
Числовая интерполяция
Вкратце он спрашивает:
Как мы можем уместить все эти точки данных с помощью функции?
Функция, о которой мы говорим, также называется..
Численное интегрирование: интегрирование Ромберга
Интегрирование с использованием численных методов
Я ценю эту концепцию интеграции прежде всего потому, что среди методов, которые мы обсуждали в классе, это единственный метод, который рассматривает возможность устранения своих ошибок. Хоть и не все, но устранение очень хорошо помогает.
Прежде всего, давайте подойдем к этому, используя также подход сверху вниз. Сначала мы обсуждаем численное интегрирование в целом, а затем переходим к интегрированию Ромберга.
Численное..
Вопросы по теме 'numerical-methods'
Реализация производной на C/C++
Как производная f(x) обычно рассчитывается программно, чтобы обеспечить максимальную точность?
Я реализую метод Ньютона-Рафсона , и он требует получения производной функции .
60103 просмотров
schedule
07.08.2023
проблема с арифметикой с использованием логарифмов, чтобы избежать числового недополнения
У меня есть два списка фракций;
скажи A = [ 1/212, 5/212, 3/212, ... ]
и B = [ 4/143, 7/143, 2/143, ... ] .
Если мы определим A' = a[0] * a[1] * a[2] * ... и B' = b[0] * b[1] * b[2] * ...
Я хочу рассчитать значения A' / B' ,
Моя...
448 просмотров
schedule
08.08.2022
Ограничение этой программы для определения суммы обратных целых чисел, не содержащих нуля
Пусть A обозначает набор положительных целых чисел, десятичное представление которых не содержит цифры 0. Сумма обратных величин элементов в A , как известно, 23.10345.
Ex....
3592 просмотров
schedule
27.10.2022
Scipy optimise.curve_fit иногда не сходится
Я пытаюсь использовать numpy.optimize.curve_fit для оценки частоты и фазы последовательности включения/выключения. Это код, который я использую:
from numpy import *
from scipy import optimize
row = array([0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0, 1.0,...
2610 просмотров
schedule
05.10.2023
получение журнала очень маленьких значений с использованием numpy/scipy в Python
У меня есть массив Nx1, который соответствует распределению вероятностей, т.е. сумма элементов равна 1. Это представлено как обычный массив numpy. Поскольку N может быть относительно большим, например. 10 или 20, многие из отдельных элементов...
9490 просмотров
schedule
01.08.2022
как я могу оценить производную сплайн-функции в R?
R может сгенерировать сплайн-функцию, используя splinefun() в библиотеке сплайнов. Однако мне нужно оценить эту функцию по ее первой и второй производной. Есть ли способ сделать это?
Например
library(splines)
x <- 1:10
y <- sin(pi/x)...
6489 просмотров
schedule
29.08.2022
Как решить задачу линейного программирования с помощью DotNumerics?
Меня очень интересует численный анализ. Я использую приложение с открытым исходным кодом DotNumerics. Моя линейная система выглядит следующим образом:
1 * x + 3 * y <= 150
2 * x + 1 * y <= 100
где x >= 0, y >= 0
z = 10 * x...
4690 просмотров
schedule
05.09.2022
Как этот линейный решатель может быть связан с Mathematica?
Вот хороший линейный решатель с именем GotoBLAS . Он доступен для загрузки и работает на большинстве вычислительных платформ. Мой вопрос: есть ли простой способ связать этот решатель с ядром Mathematica, чтобы мы могли называть его как LinearSolve?...
531 просмотров
schedule
06.07.2023
Решение переопределенной системы ограничений
У меня есть n числовых переменных (не знаю, пофиг), назовем их X[n] . У меня также есть отношения m >> n между ними, назовем их R[m] , вида:
X[i] = alpha*X[j] , alpha — ненулевое положительное действительное число, i и j...
6060 просмотров
schedule
14.05.2022
Размытие по Гауссу с переменным радиусом, приближающееся к ядру
Я пишу размытие по Гауссу с переменным радиусом (стандартное отклонение), то есть каждый пиксель изображения свернут с использованием другого ядра. Стандартные методы вычисления размытия по Гауссу здесь не работают: БПФ, разделение осей,...
1856 просмотров
schedule
18.08.2022
Численная библиотека для минимизации многомерной функции на окнах
Нам нужно численно минимизировать функцию, которая принимает трехмерный вектор в качестве входных данных. Функция гладкая, поэтому хорошим выбором будет градиентный алгоритм.
Однако я привык к GSL, который, к сожалению, требует gcc. Однако мы...
256 просмотров
schedule
08.05.2023
Численные рецепты / Многомерный поиск корня (с использованием тритона): Как минимизировать максимальную ошибку
Этот вопрос связан с книгой «Численные рецепты на C++», поэтому он будет предназначен для людей, которые немного разбираются в этом, а также в многомерной оптимизации.
Я пишу программу, которая должна искать многомерный корень, и для ее решения я...
1361 просмотров
schedule
13.08.2022
Числовая библиотека для Scala
Я ищу библиотеку для численных вычислений в Scala (или Java, хотя что-то, что может использовать функции Scala, было бы намного лучше!) По крайней мере, со следующими возможностями:
L-БФГС
Минимизаторы (Пауэлл, КвазиНьютон, ...)
Численное...
2318 просмотров
schedule
15.05.2022
F# Microsoft Solver Foundation — класс NelderMeadSolver
Может ли кто-нибудь показать мне пример кода для использования класса NelderMeadSolver в F #?
Например, я хочу минимизировать следующую функцию: F(X, Y)
F = (X-1)^2 + (y-1)^2, где 0‹ X ‹ 2 , 0‹ Y ‹ 2 Ответ очевиден: X = 1, Y = 1
Я нашел...
900 просмотров
schedule
06.11.2022
Как построить функцию, используя код, который я написал с отправной точкой?
Я написал код алгоритма secant , и теперь у меня есть функция:
f(x) = 2x^3 - 4x^2 + 3x , с двумя начальными точками : x0 = -1 , x1 = 2 .
Мой вопрос в том, как я могу построить функцию, которую я написал, т.е. secant , с функцией выше ,...
776 просмотров
schedule
22.07.2023
Преобразование векторного ввода в алгоритм из матричного вывода
Я написал код для реализации алгоритма, который принимает на вход вектор q действительных чисел и возвращает на выходе сложную матрицу R . Приведенный ниже код Matlab создает график, показывающий входной вектор q и выходную матрицу R ....
313 просмотров
schedule
04.08.2022
Интеграция функции с использованием неравномерной меры (python/scipy)
Я хотел бы интегрировать функцию в python и указать плотность вероятности (меру), используемую для выборки значений. Если это не очевидно, интегрирование f(x)dx в [a,b] неявно использует равномерную плотность вероятности по [a,b] , и я хотел бы...
440 просмотров
schedule
17.06.2022
Линейная свертка и взаимная корреляция двух векторов в Matlab с использованием БПФ и вывода одинаковой длины
Я реализую алгоритм, в котором мне нужно вычислить линейную свертку и взаимную корреляцию между двумя одномерными векторами a и b с помощью БПФ. Если длина a равна m , а длина b равна n , то общая длина результирующей свертки будет m + n -...
3467 просмотров
schedule
28.02.2023
Численное решение двух связанных нелинейных дифференциальных уравнений второго порядка
Я столкнулся со следующей системой дифференциальных уравнений лагранжевой механики. Можете ли вы предложить численный метод с соответствующими ссылками и ссылками на то, как я могу его решить. Кроме того, есть ли более короткая реализация на Matlab...
2091 просмотров
schedule
12.02.2023
Преобразовать двоичную строку длиной более 52 бит в одну в Matlab?
Я пытаюсь преобразовать очень длинные двоичные строки, часто превышающие 52 бита, в числа. У меня не может быть фиксированного окна просмотра вперед, потому что я делаю это, чтобы вычислить версию сложности Лемпеля-Зива для нейронных данных.
Когда...
2199 просмотров
schedule
22.10.2022