Вопросы по теме 'numerical-stability'
Стратегии устранения проблем с числовой стабильностью?
Я пытаюсь написать реализацию алгоритма факторизации спектральной плотности Уилсона [1] для Python. Алгоритм итеративно разлагает матричную функцию [QxQ] на ее квадратный корень (это своего рода расширение искателя квадратного корня Ньютона-Рафсона...
652 просмотров
schedule
05.05.2022
Вычисление обратной матрицы в Matlab
Я запускаю алгоритм оптимизации, который требует вычисления обратной матрицы. Цель алгоритма — исключить отрицательные значения из матрицы A и получить новую матрицу B. По сути, я начинаю с известных квадратных матриц B и C одинакового размера.
Я...
21552 просмотров
schedule
20.11.2022
Предотвращение числового переполнения при расчете значения И градиента функции логистических потерь
В настоящее время я пытаюсь реализовать алгоритм машинного обучения, который включает функцию логистических потерь . в МАТЛАБ. К сожалению, у меня возникли проблемы из-за числового переполнения.
В общем, для заданного входа s значение...
3700 просмотров
schedule
11.08.2022
ограничение ошибки в алгоритме аппроксимации функции
Предположим, у нас есть набор чисел с плавающей запятой с битовой мантиссой «m» и битами «e» для экспоненты. Предположим, что мы хотим аппроксимировать функцию "f".
Из теории мы знаем, что обычно используется «функция с уменьшенным диапазоном», а...
198 просмотров
schedule
25.06.2023
Взаимозаменяемость операций сложения и умножения с плавающей запятой IEEE 754
Взаимозаменяемо ли сложение x + x на умножение 2 * x в стандарте IEEE 754 (IEC 559) с плавающей запятой или, в более общем смысле, есть ли гарантия, что case_add и case_mul всегда дают точно тот же результат?
#include <limits>...
244 просмотров
schedule
17.05.2022
Один и тот же код интеграции в 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