Вопросы по теме 'numerical-methods'
Как сделать топографическую карту из разреженных выборочных данных?
Мне нужно сделать топографическую карту местности, для которой у меня есть только довольно редкие выборки данных (x, y, высота) . Очевидно, что я не могу сделать полностью точную карту, но хотелось бы, чтобы она была в каком-то смысле «гладкой»....
3679 просмотров
schedule
08.10.2022
Библиотека автоматического дифференцирования в Scheme / Common Lisp / Clojure
Я слышал, что одной из первоначальных мотиваций Маккарти для изобретения Лиспа было создание системы для автоматического распознавания. Несмотря на это, мой поиск в Google не дал никаких библиотек / макросов для этого. Существуют ли какие-либо...
3264 просмотров
schedule
16.10.2022
Численное интегрирование функции со значениями, известными в заданном наборе точек (конечном и дискретном) по области, ограниченной дискретными точками?
Пусть D будет областью, ограниченной серией точек {x_i,y_i} (1‹=i‹=N). (Площадь не обязательно должна быть выпуклой, и предполагается, что точки идут вдоль граничной кривой.)
Пусть f — функция, определенная на D, но мы знаем только ее значения в...
790 просмотров
schedule
04.01.2024
МАТЛАБ; Многочлен Лагранжа, интересная ошибка
В школьном домашнем задании мне дали задачу на вычисление многочлена Лагранжа. Набор данных (x,y) включал x = 0,1,2....20, а y = x - 0,3*rand() + 1; Затем нас попросили наблюдать за феноменом Рунге . После этого нас попросили вычислить узлы...
1887 просмотров
schedule
29.02.2024
Численное интегрирование сложной функции с помощью пакета Cubature C
Я хочу использовать пакет C Cubature для выполнения многомерного интеграла сложная функция. Я попытался сделать это следующим образом для очень простой функции f (x, y) = x + y * i над квадратом [0,1] x [0,1]. Точный результат - 0,5 + 0,5i....
1175 просмотров
schedule
14.11.2022
Подогнать дистрибутив к счетчику в scipy
У меня есть объект collections.Counter со счетчиком вхождений разных значений, например:
1:193260
2:51794
3:19112
4:9250
5:6486
Как я могу подогнать распределение вероятностей к этим данным в scipy? scipy.stats.expon.fit(), кажется, хочет...
302 просмотров
schedule
25.03.2024
Производительность JavaScript Math.sqrt
Профилируя код, я нашел файл Math.sqrt. функция конкретно является основным узким местом в большом двойном вложенном цикле, который выполняется на каждом временном шаге в моей программе. Есть ли способ улучшить его производительность? Должен ли...
1473 просмотров
schedule
31.03.2024
Как я могу получить все решения этого уравнения в MATLAB?
Я хотел бы решить следующее уравнение: tan(x) = 1/x
Что я сделал:
syms x
eq = tan(x) == 1/x;
sol = solve(eq,x)
Но это дает мне только одно численное приближение решения. После этого я прочитал о следующем:
[sol, params, conds] =...
757 просмотров
schedule
13.12.2022
Одномерный поиск корня с минимальным количеством вычислений функции
Я ищу алгоритмы поиска корней, которые используют очень мало вычислений функций (цель - минимум). Задача поиска корней имеет следующие характеристики:
f(x) = 0, R -> R
оценка функции ( f(.) ) чрезвычайно затратна*;
для начала доступен...
167 просмотров
schedule
26.12.2023
Процедура минимизации быстро уменьшается, а затем увеличивается
Я запустил алгоритм "внутренней точки" FMINCON от MATLAB. Алгоритм использует BFGS и, таким образом, учитывает информацию в числовом гессиане.
Я оптимизирую пространство параметров, которое заведомо невыпукло по отношению к уникальному решению. Я...
12 просмотров
schedule
16.09.2022
Рунге-кутта для спаренных ОДУ
Я создаю функцию в Octave, которая может решать N связанное обыкновенное дифференциальное уравнение типа:
dx/dt = F(x,y,…,z,t)
dy/dt = G(x,y,…,z,t)
dz/dt = H(x,y,…,z,t)
Любым из этих трех методов (Эйлера, Хойна и Рунге-Кутта-4)....
781 просмотров
schedule
24.01.2024
Модифицированный алгоритм метода секущей Javascript?
Мне нужна ваша помощь, чтобы закодировать этот новый метод секущих в Javascript.
В методе секущих используются два начальных параметра, тогда как для этого нового метода требуется только один начальный параметр.
Заранее спасибо.
Спасибо...
366 просмотров
schedule
04.01.2024
Наименьшие квадраты с использованием общего матричного векторного умножения, а не разреженных матриц
Есть ли способ вычислить
\argmin_{x}\|Ax-b\|_2
на основе функции, которая вычисляет матричные векторные произведения Ax без явного сохранения, разреженного или неразреженного, A в памяти?
В Python я бы использовал для этого...
70 просмотров
schedule
24.10.2022
Python resolve_bvp дифференциальное уравнение четвертого порядка
Я пытаюсь применить solve_bvp scipy к следующей проблеме
T''''(z) = -k^4 * T(z)
С граничными условиями для домена размером l и некоторой константой A :
T(0) = T''(0) = T'''(l) = 0
T'(l) = A
До сих пор я свел уравнение...
175 просмотров
schedule
07.02.2024
Нормы в Python для чисел с плавающей запятой и десятичных чисел (с фиксированной точкой)
Рекомендуется ли использовать собственную реализацию Python с плавающей запятой или ее decimal для случаи использования, где важна точность?
Я думал, что на этот вопрос будет легко ответить: если накопленная ошибка имеет значительные последствия,...
415 просмотров
schedule
09.12.2023
Тело этой функции Matlab работает, но не сама функция (ошибка interp1)
Я написал следующий фрагмент подкода (с комментариями параметров) для алгоритма итерации политики Эйлера. Когда я пытаюсь запустить тело функции (все, что ниже глобального), например, a1 = 1, оно работает и возвращает скаляр. Однако, когда я вызываю...
27 просмотров
schedule
09.12.2023
Задача с корнями нелинейного уравнения
У меня есть гиперболическая функция, и мне нужно найти ее 0. Пробовал различные классические методы (деление пополам, Ньютон и т.д.).
Вторые производные непрерывны, но недоступны аналитически, поэтому я должен исключить методы, использующие их....
61 просмотров
schedule
06.06.2024
Лучший способ решить систему нелинейных уравнений, где неизвестными являются векторы?
У меня есть следующие уравнения, которые описывают геометрическую обратную задачу, которую я хотел бы решить:
В общем, задача трехмерная, поэтому для каждого из перечисленных векторных уравнений я получаю три уравнения (по одному на каждое...
55 просмотров
schedule
09.03.2024