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

Вопросы с тегом
Численный анализ - это исследование алгоритмов, использующих численное приближение, в отличие от общих символьных манипуляций.
3 ответа

MATLAB рассчитать угол поворота

У меня 5 очков по кругу:.......... 1............................................ 2........ 3 4 5 Теперь у меня есть его повернутое изображение:............. ............................................................. 1........ 3 2 5 4 Примечание: …
21 июн '12 в 06:41
1 ответ

Собственные значения разреженной матрицы с использованием Eigen3

Существует ли четкий и эффективный способ нахождения собственных значений и собственных векторов реальной, симметричной, очень большой, скажем, 10000x10000, разреженной матрицы в Eigen3? Для плотных матриц существует решатель собственных значений, н…
14 май '15 в 19:52
6 ответов

Анализируя шумные данные

Недавно я запустил ракету с барометрическим высотомером, точность которого составляет примерно 10 футов (рассчитывается по данным, полученным во время полета). Записанные данные имеют временные приращения 0,05 с на выборку, и график зависимости высо…
24 дек '09 в 05:28
1 ответ

Волновое уравнение с FDM, Matlab

Я пытаюсь реализовать код Matlab для решения волнового уравнения, моя функция выглядит следующим образом: function [x,t,w] = wave_eqn(xl,xr,yb,yt,M,N,f,l,r,p) % input: space interval [xl,xr], time interval [yb,yt] % number of space steps M, number o…
26 май '17 в 10:01
2 ответа

Числовая точность в Фортране 95:

У меня есть следующий код Fortran: Program Strange Real(Kind=8)::Pi1=3.1415926535897932384626433832795028841971693993751058209; Real(Kind=8)::Pi2=3.1415926535897932384626433832795028841971693993751058209_8; Print*, "Pi1=", Pi1; Print*, "Pi2=", Pi2; …
04 май '13 в 03:11
5 ответов

Библиотека обработки сигналов в Java?

Я хотел бы вычислить спектральную плотность мощности временных рядов; выполнить полосовую, низкочастотную и высокочастотную фильтрацию; может быть, некоторые другие основные вещи. Есть ли хорошая библиотека Java с открытым исходным кодом, чтобы сдел…
11 мар '09 в 22:53
0 ответов

matlab fit exp2

Я безуспешно ищу документацию о функции подгонки, используя exp2 (сумма 2 показателей). Как работает функция, понятно: [curve, gof] = fit(x, y,'exp2'); Но поскольку существует несколько способов подбора суммы показателей, я пытаюсь выяснить, какой а…
02 сен '12 в 09:16
2 ответа

Реализация функции Ламберта W в Java

Я работаю над проектом и попал в ситуацию, когда мне нужна функция, чтобы иметь возможность получить хотя бы приближение значения W(x), функции W Ламберта, где x может быть любым действительным числом. Я работаю на Java. Я не мог найти какие-либо ре…
22 июл '11 в 16:19
1 ответ

Как найти числовую производную функции внутри системы уравнений Диффа?

Я хочу решить эту систему DE: function ydot= DRV(y) B=[1 0.5 4;7.2 0.6 5;3.3 8 0.2]; R1=[2 5 3;11 3 6;1.2 2 4]; Q1=eye(3,3); qv=[y(1);y(2);y(3)]; p=[y(4) y(5) y(6);y(7) y(8) y(9);y(10) y(11) y(12)]; qvdot=p*qv; pdot=p*B*inv(R1)*p+Q1; v=Q1*p*R1*qv; v…
23 июл '18 в 01:14
1 ответ

Численные подводные камни без FPU повышенной точности

В превосходной статье Моннио "Подводные камни проверки арифметики с плавающей точкой" приводятся примеры неожиданного численного поведения. Однако большинство примеров зависят от FPU с расширенной точностью или операций с повышенной точностью (напри…
26 июл '11 в 12:33
0 ответов

Кодирование кубического сплайна MATLAB

Вопрос: У меня практически нет опыта работы с MATLAB, но я попытался найти учебники в Интернете с примерами кубических сплайнов, но без особой помощи. Я попытался сделать следующее, но, честно говоря, я не знаю, правильно это или нет, или что делать…
17 ноя '15 в 21:06
1 ответ

C++ - обратная замена через Boost - ошибка неверного индекса

У меня есть следующий код для алгоритма обратной замены: #include <boost\numeric\ublas\matrix.hpp> typedef boost::numeric::ublas::matrix<double> dM; dM bSub(dM A, dM b) { unsigned int n = (int)b.size1(); assert(b.size2() == 1); assert(n …
27 апр '16 в 23:01
2 ответа

Построение интегральной таблицы Ромберга

Я пытаюсь написать программу на Фортране для генерации таблицы интеграции Ромберга. В книге " Численный анализ" Р.Л.Бурдена и JDFaires, 9-е изд., Есть алгоритм . в главе 4.5. Пока я написал это implicit none integer,parameter::n=4 real::a,b,f,r(n,n)…
30 авг '15 в 20:54
4 ответа

Рассчитать мощность с фиксированным (постоянным) показателем

У меня есть случай, когда мне нужно вычислить x^y a огромное количество раз, где y является константой, а x гарантированно будет действительным числом. Как это может быть сделано более эффективно, чем Pow(double x, double y), который будет выполнять…
23 окт '12 в 21:09
2 ответа

Ошибка чисел Python: большое число добавить небольшое число

Это упражнение из курса глубокого обучения Udacity. Может кто-нибудь объяснить, почему окончательный ответ не 1,0? v1 = 1e9 v2 = 1e-6 for i in range(int(1e6)): v1 = v1 + v2 print 'answer is', v1 - 1e9 # answer is 0.953674316406
17 мар '16 в 02:13
1 ответ

Оценка погрешности добавления вектора

Кахан в " Кресте" на странице 14 показывает, что угол между истинным перекрестным произведением и углом, вычисленным с использованием арифметики с плавающей запятой, ухудшается, когда угол между двумя скрещенными векторами приближается к нулю или π,…
09 дек '14 в 21:44
2 ответа

Алгоритм квадратного уравнения MATLAB

Какое условие я должен поставить в коде Matlab, чтобы получить точное решение квадратичного с этими формулами: x1=(-2*c)/(b+sqrt(b^2-4*a*c)) x2=(-2*c)/(b-sqrt(b^2-4*a*c)) Непосредственно применяя эти формулы, я не получаю правильного решения в некот…
14 мар '14 в 19:58
0 ответов

Программная реализация функции квантиля для t-распределения

Я хочу внедрить функцию квантиля в один из моих документов в LaTeX, В частности, квантильная функция для распределения t. Да, я говорю о фактическом численном приближении, например, qt() функция в Rвнутри моего документа. Как обычно осуществляется к…
07 июн '14 в 11:31
0 ответов

Java - численное интегрирование комплексной функции - дзета-функция, формула Абеля-Планы

(Примечание: я нашел частичное решение. В конце я вставил результаты теста) Я хочу численно интегрировать приближение дзета-функции. Это делается с помощью так называемой формулы Абеля-Планы. Формула Абеля-Планы может быть использована для численной…
22 авг '15 в 16:23
2 ответа

MATLAB - Адаптивный размер шага Рунге-Кутта

Я запрограммировал в MATLAB адаптивный размер шага RK4 для решения системы ODE. Код выполняется без ошибок, однако он не дает желаемой кривой, когда я пытаюсь построить x против y. Вместо того, чтобы иметь форму тороида, я просто получаю плоскую лин…
16 сен '16 в 05:49