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

Стратегии устранения проблем с числовой стабильностью?
Я пытаюсь написать реализацию алгоритма факторизации спектральной плотности Уилсона [1] для Python. Алгоритм итеративно разлагает матричную функцию [QxQ] на ее квадратный корень (это своего рода расширение искателя квадратного корня Ньютона-Рафсона...
652 просмотров

Вычисление обратной матрицы в Matlab
Я запускаю алгоритм оптимизации, который требует вычисления обратной матрицы. Цель алгоритма — исключить отрицательные значения из матрицы A и получить новую матрицу B. По сути, я начинаю с известных квадратных матриц B и C одинакового размера. Я...
21552 просмотров

Предотвращение числового переполнения при расчете значения И градиента функции логистических потерь
В настоящее время я пытаюсь реализовать алгоритм машинного обучения, который включает функцию логистических потерь . в МАТЛАБ. К сожалению, у меня возникли проблемы из-за числового переполнения. В общем, для заданного входа s значение...
3700 просмотров

ограничение ошибки в алгоритме аппроксимации функции
Предположим, у нас есть набор чисел с плавающей запятой с битовой мантиссой «m» и битами «e» для экспоненты. Предположим, что мы хотим аппроксимировать функцию "f". Из теории мы знаем, что обычно используется «функция с уменьшенным диапазоном», а...
198 просмотров

Взаимозаменяемость операций сложения и умножения с плавающей запятой IEEE 754
Взаимозаменяемо ли сложение x + x на умножение 2 * x в стандарте IEEE 754 (IEC 559) с плавающей запятой или, в более общем смысле, есть ли гарантия, что case_add и case_mul всегда дают точно тот же результат? #include <limits>...
244 просмотров

Один и тот же код интеграции в Matlab и python, Matlab стабилен, python взрывается
Вот алгоритм метода экспоненциальной разницы во времени с использованием исходного кода Matlab из Оксфорда. clc clear % viscosity nu = 0.5; % Spatial grid and initial condition: N = 128; x = 2*pi*(0:N-1)'/N; u = cos(x).*(1+sin(x)); v = fft(u);...
205 просмотров
schedule 08.03.2023