Описание тега numerical-methods

Алгоритмы, которые решают математические задачи с помощью численного приближения (в отличие от символьных вычислений).
0 ответов

Сопоставить числовые точки данных с заранее определенной сеткой

У меня есть данные из какого-то эксперимента, которые должны напоминать известную сеточную структуру. Тем не менее, данные выключен неизвестным углом искажается неизвестным образом может содержать некоторые смещенные точки Для примера, пожалуйста, п…
12 авг '14 в 14:06
2 ответа

Полное сингулярное разложение в R

В большинстве приложений (особенно статистических) достаточно тонкого СВД. Однако иногда требуется полный SVD, чтобы получить ортобазис нулевого пространства матрицы (и его сопряженного). Похоже, что svd() в R возвращает только тонкую версию. Можно …
05 июл '11 в 02:35
0 ответов

Каковы элементы вектора поворота в устранении Гаусса-Иордана

В числовых рецептах на стр. 39 (стр. 4 в pdf) был предложен следующий алгоритм для нахождения стержня: void gaussj(float **a, int n, float **b, int m) /* Linear equation solution by Gauss-Jordan elimination, equation (2.1.1) above. a[1..n][1..n] is …
21 фев '16 в 07:05
2 ответа

MATLAB: коэффициенты фильтра IIR

Я довольно новичок в обработке сигналов, и один из моих проектов - реализовать класс фильтра C++. Мне нужны коэффициенты более высокого порядка для типичных фильтров, таких как Чебышевские типы I и II, Баттерворт, Эллиптик, и, к сожалению, большинст…
1 ответ

Преобразовать числовую модель теплопроводности Matlab в Python

Я пытаюсь преобразовать свою модель Matlab для переходной теплопроводности в Python. К сожалению, результаты моего численного решения в Python не совпадают с результатами модели Matlab. Я использую Spyder IDE для написания своего кода. Основные разл…
2 ответа

Почему более эффективно реализовать log_sum, как это?

Согласно этой теме: Вычислить журнал (a+b) Иногда log_sum реализован так: log(a + b) = log(a * (1 + b/a)) = log a + log(1 + b/a) Я не понимаю, почему этот подход более эффективен. У кого-нибудь есть идеи по этому поводу?
20 ноя '16 в 07:17
1 ответ

Лучший способ добавить 3 числа (или 4, или N) в Java - Kahan Sums?

Я нашел совершенно другой ответ на этот вопрос, весь оригинальный вопрос больше не имеет смысла. Однако способ ответа будет полезен, поэтому я немного его модифицирую... Я хочу подвести итог три double числа, скажем a, b, а также c наиболее устойчив…
0 ответов

Решение наименьших квадратов для матрицы вместо вектора

Задача состоит в том, чтобы найти Z таким, чтобы эпсилон (E) (уравнение 21) был минимизирован. Z - это матрица MxN, которую мы должны найти. Zx и Zy также являются матрицами MxN, которые также уже известны. Dx и Dy - это матрицы NxN, которые проводя…
0 ответов

Определитель матрицы Якоби точек сетки

У меня есть функция, в которой ее якобиан является матрицей 2x2, я вычислил ее для любого вектора (teta1,teta5), а затем вычислил якобиан. from numpy import linspace, meshgrid from numpy import linalg as LA import numdifftools as nd tetas = [np.radi…
09 янв '18 в 16:28
3 ответа

Как сделать эффективную для памяти двумерную свертку на больших массивах

У меня есть проблема, когда мне нужно свернуть один очень большой 2D-массив (файл на диске) с меньшим массивом, который помещается в памяти. scipy.signal.fftconvolve замечательно, когда массивы умещаются в памяти, но не помогают, когда нет. Есть ли …
16 окт '14 в 16:47
1 ответ

Выявление ошибок в квадратных уравнениях

Предположим, вы пишете программу на Python для вычисления действительных корней из ax2 + bx + c = 0, где заданы коэффициенты a, b и c, действительные числа. Традиционные формулы для двух решений: x1 = (-b - sqrt(b*b-4*a*c)) / 2a, x2 = (-b + sqrt(b*b…
17 сен '13 в 19:55
2 ответа

Сюжет двух серий в матлабе

Я хотел бы наметить значения r1 а также r2 на одном участке, как я могу это сделать? for i=1:10 r1=rand(1) r2=max(rand(1,2)) end Спасибо за помощь!
24 апр '14 в 21:24
2 ответа

VHDL для научных вычислений

Мне было интересно, если люди используют VHDL/FPGA в научных вычислениях. Пример сценария, о котором я думал, был сказать: Построить сумматор с плавающей запятой произвольной точности Настройте плату ПЛИС, чтобы затем добавить такие номера Поэтому я…
1 ответ

Численный метод Рунге-Кутты Плохое приближение

Я пытаюсь использовать метод Рунге-Кутты, чтобы сравнить его с функцией lsode. Но он работает довольно плохо, все остальные методы, которые я использовал (Фовардс и Эйлер, Хеун) для сравнения с lsode, делают работу лучше, до такой степени, что они п…
12 ноя '18 в 00:18
1 ответ

Интегральная функция в R

Я определил следующую функцию в R (где a это заданная константа и это известная функция), но когда я пытаюсь вычислить ее значения, я получаю довольно странные числа. Я должен отметить, что для того, чтобы вычислить значения f(s)Я на самом деле реши…
1 ответ

Интеграция осциллирующей (интерполированной) функции с GSL и python

Я пытаюсь интегрировать сильно колебательные данные с функцией qawo из научной библиотеки GSL и pygsl в python. Поскольку я имею дело с данными, я думал, что интерполированная функция могла бы работать, но GSL дает мне неверный результат !!, Позволь…
1 ответ

Использование памяти odeint и время выполнения

Я заметил, что odeint использует очень мало памяти по сравнению с моей реализацией алгоритма RK4 или Mathematica. Для того же размера шага, odeint использует около 3,11 ГБ, в то время как моя программа использует 7 ГБ, а в Mathematica я должен вручн…
30 авг '12 в 18:04
1 ответ

Проблема с числовой интеграцией

В настоящее время я застрял на действительно простой ошибке, спрятанной где-то, и, надеюсь, кто-то сможет пролить свет на это. Я пытаюсь численно интегрировать набор данных с обычными методами, и получаю результаты, которые я не совсем ожидал. Затем…
1 ответ

Решение Ax =B с помощью UMFPACK

Я использую UMFPACK для решения Ax = B, где A и B - матрицы. В настоящее время в моей реализации я решаю x(:,i) = A\b(:,i) один за другим, чтобы получить x. Есть ли способ напрямую получить x без необходимости циклически перебирать каждый столбец B,…
11 фев '16 в 17:52
1 ответ

Заставить numy fsolve работать над кусочно-постоянными функциями

Я пытаюсь использовать fsolve, чтобы найти точки пересечения различных параметрических кусочно-определенных функций: max_price = 20.0 b = 0.5 consumption_func = lambda x,b : max_price if x <= b else max_price - (x-b) func = lambda x: consumption_…
09 дек '16 в 09:47