Описание тега linear-algebra
Линейная алгебра - это основная часть математики, изучающая векторные пространства и линейные отображения между этими пространствами. Ключевые темы включают: матрицы, векторы, одновременные системы линейных уравнений, диагональную форму, собственные значения, собственные векторы, проекции, точки, линии, плоскости, гиперплоскости, вращения и отражения.
3
ответа
Повторная настройка уравнения с плавающей точкой
Я хотел бы знать, есть ли способ повысить точность расчета наклона. (Это появилось несколько месяцев назад здесь). Кажется, изменив: float get_slope(float dXa, float dXb, float dYa, float dYb) { return (dXa - dXb)/(dYa - dYb); } в float get_slope(fl…
27 май '09 в 18:51
2
ответа
Полное сингулярное разложение в R
В большинстве приложений (особенно статистических) достаточно тонкого СВД. Однако иногда требуется полный SVD, чтобы получить ортобазис нулевого пространства матрицы (и его сопряженного). Похоже, что svd() в R возвращает только тонкую версию. Можно …
05 июл '11 в 02:35
0
ответов
Как использовать ускорение BLAS для решения линейного уравнения разреженной матрицы?
В приложении для iOS мне нужно решить линейное уравнение Ax = B, где A - разреженная матрица с 40K строк и 10K столбцов. У Accelerate есть пакет Sparse Solver, но он все еще находится в бета-версии: https://developer.apple.com/documentation/accelera…
07 сен '17 в 06:44
2
ответа
Решение систем уравнений над конечными полями
Мне нужно инвертировать довольно большую (n > 1000) систему уравнений в конечном поле (по модулю 2^32-5, если быть точным). Есть ли библиотеки, которые могут сделать это эффективно? Желательно в Java, но с C все в порядке (у нас уже есть некоторые J…
27 окт '11 в 13:24
0
ответов
Как получить все собственные значения из разреженной матрицы с помощью eigs_gen
Я использую броненосец. eigs_gen функция (для SpMat разреженных матриц) имеет параметр k для числа собственных значений для вычисления. У меня матрица 3х3 my_matrix когда я бегу arma::cx_fvec values; arma::cx_fmat vectors; arma::eigs_gen (values, ve…
25 окт '17 в 20:31
1
ответ
Инверсия матриц, умножение и собственные вычисления в Hadoop
Мне нужно выполнить эти три шага: - Найти инверсию квадратной симметричной матрицы - Умножить результат на другую квадратную симметричную матрицу - Наконец, найти собственные векторы и собственные значения результирующей матрицы Код Python: S_i = np…
26 сен '14 в 05:38
1
ответ
Есть ли у Numpy встроенная поэлементно-матричная реализация модульного возведения в степень
Есть ли у numpy встроенная реализация для модульного экспонирования матриц? (Как указывает user2357112, я на самом деле ищу поэлементное модульное сокращение) Один из способов модульного возведения в степень для обычных чисел выполняется с помощью в…
22 июл '17 в 18:05
2
ответа
Как распределить записи матрицы в большую с использованием продуктов матрицы
Учитывая матрицу A, где A есть x_11 x_12 x_13 x_21 x_22 x_23 x_31 x_32 x_33 x_41 x_42 x_43 Как я могу эффективно создать эту вторую матрицу, используя продукты матрицы? x_11 0 0 x_21 0 0 x_31 0 0 x_41 0 0 0 x_12 0 0 x_22 0 0 x_32 0 0 x_42 0 0 0 x_13…
06 июн '13 в 11:23
2
ответа
Как получить реальные собственные значения и собственные векторы для трехдиагональной матрицы Тёплица?
Я строю матрицу 100*100 k и хочу использовать numpy.linalg.eig диагонализировать это. k=np.zeros((100,100)) np.fill_diagonal(k,-2) np.fill_diagonal(k[1:,:-1],1.5) np.fill_diagonal(k[:-1,1:],0.5) когда я пытаюсь уменьшить матрицу, такую как w,v=np.…
27 фев '18 в 02:01
2
ответа
Обратная матрица и умножение
Я новичок в мире матриц, извините за этот основной вопрос, который я не мог понять: У меня четыре матрицы (одна неизвестная). Матрица Х x <- c(44.412, 0.238, -0.027, 93.128, 0.238, 0.427, -0.193, 0.673, 0.027, -0.193, 0.094, -0.428, 93.128, 0.673…
09 ноя '11 в 12:35
1
ответ
Определение 3D-сцены из фотографии круга
Учитывая фото, содержащее круг, например, это фото фонтана: Можно ли определить 3D положение и поворот фонтана относительно камеры? Я понимаю, что мы должны определить масштаб, поэтому предположим, что фонтан имеет ширину 2 метра (диаметр круга, сос…
18 янв '16 в 15:07
1
ответ
Поддерживает ли Intel MKL Armadilo, Eigen или Biltz++?
Я имею в виду: может ли он определить некоторый массив, как это делает Blitz ++ (я ищу похожий простой синтаксис, нет необходимости быть точно таким же.): Array<double, 3> A(NX,NY,NZ) Благодарю. Я слышал, что Armadilo, Eigen, Boost и Biltz++ с…
07 май '13 в 04:10
1
ответ
Как добавить регуляризацию в линейном программировании Scipy (неотрицательные наименьшие квадраты)
Вот мой код LP с использованием NNLS Scipy: import numpy as np from numpy import array from scipy.optimize import nnls def by_nnls(A=None, B=None): """ Linear programming by NNLS """ #print "NOF row = ", A.shape[0] A = np.nan_to_num(A) B = np.nan_to…
16 фев '16 в 02:24
1
ответ
Как выразить c[i,j,k] = a[i,j] * b[i,k] в Numpy/Theano?
Определение c[i,j,k] = a[i,j] * b[i,k] for any i, j, k является поэлементным продуктом по отношению к iи внешний продукт по отношению к j а также k, Есть ли способ выразить это в NumPy/Theano без петель?
02 июн '15 в 12:15
2
ответа
Решение наименьших квадратов для матрицы вместо вектора
Задача состоит в том, чтобы найти Z таким, чтобы эпсилон (E) (уравнение 21) был минимизирован. Z - это матрица MxN, которую мы должны найти. Zx и Zy также являются матрицами MxN, которые также уже известны. Dx и Dy - это матрицы NxN, которые проводя…
02 июн '15 в 19:43
1
ответ
Как вставить столбец единиц и нулей в матрицу, используя Octave?
Предположим, у меня есть матрица с набором целых чисел. Я хочу использовать чек rand > 0.5 добавить к моей матрице случайный вектор 1 и 0. Как я мог это сделать?
20 сен '18 в 01:52
1
ответ
Как я могу узнать, какие размеры являются основным компонентом?
Я использую функцию princomp в Matlab для PCA. Насколько я понимаю, я мог проверить латентность, чтобы решить, сколько измерений мне нужно. [coeff, score, latent, t2] = princomp(fdata); cumsum(latent)./sum(latent); И с помощью trainMatrix = coeff(:,…
24 апр '13 в 22:00
4
ответа
Решение системы линейных уравнений в неквадратной матрице
У меня есть система линейных уравнений, которые составляют NxM матрица (т.е. не квадратная), которую мне нужно решить - или, по крайней мере, попытаться решить, чтобы показать, что нет решения для системы. (скорее всего, не будет никакого решения) Н…
25 окт '11 в 17:02
2
ответа
Как решить A*X - X*A' = 0
У меня есть уравнение вида A*X = X*A'где A и X - действительные, квадратные матрицы (3x3 в этом случае), и A известен, а A'представляют транспонирование A. Как решить для X использование MATLAB? (с точностью до масштабного коэффициента)
30 авг '16 в 06:20
1
ответ
Как извлечь только 3 собственных вектора изображения nxn в opencv?
Я пытаюсь преобразовать RGB-изображение в оттенки серого, используя следующую статью. Основной алгоритм, используемый в статье: алгоритм на основе нового PCA для преобразования изображений в оттенки серого Однако, когда я пытаюсь извлечь собственные…
07 июн '16 в 10:21