Результаты поиска по запросу "linear-algebra"

2 ответа

Решение неквадратной линейной системы с помощью R

Как решить неквадратную линейную систему с R:A X = B ? (в случае, если система не имеет решения или бесконечно много решений) Пример : A=matrix(c(0,1,-2,3,5,-3,1,-2,5,-2,-1,1),3,4,T) B=matrix(c(-17,28,11),3,1,T) A [,1] [,2] [,3] [,4] [1,] 0 1 ...

1 ответ

Numpy: преобразовать массив в треугольную матрицу

Я искал встроенный метод для преобразования линейного массива в треугольную матрицу. Поскольку мне не удалось найти один, я прошу помощи в его реализации. Вообразите массив как: In [203]: dm Out[203]: array([ 0.80487805, 0.90243902, 0.85365854, ...

3 ответа

NumPy векторизация двойного питона для цикла

V это (n, p) numpy массив, как правило, размеры n ~ 10, р ~ 20000 Код у меня сейчас выглядит A = np.zeros(p) for i in xrange(n): for j in xrange(i+1): A += F[i,j] * V[i,:] * V[j,:]Как бы я переписал это, чтобы избежать двойного цикла python for?

ТОП публикаций

2 ответа

Почему Einsum Numpy медленнее, чем встроенные функции Numpy?

Я обычно получаю хорошую производительность от функции einsum numpy (и мне нравится ее синтаксис). @ Ответ Офионаэтот вопрос [https://stackoverflow.com/questions/18365073/why-is-numpys-einsum-faster-than-numpys-built-in-functions] показывает, что ...

2 ответа

Вычислить жорданову нормальную форму матрицы в Python / NumPy

В MATLAB вы можете вычислить нормальную жорданову форму матрицы с помощью функцииjordan. Есть ли в NumPy и SciPy эквивалентная функция?

4 ответа

Построить матрицу полного ранга, добавив векторы из стандартной базы

У меня есть nxn особая матрица. Я хочу добавить к этой матрице k строк (которые должны быть из стандартного базиса e1, e2, ..., en) так, чтобы новая (n + k) xn матрица имела полный ранг столбца. Количество добавляемых строк k должно ...

2 ответа

PCA Уменьшение размерности

Я пытаюсь выполнитьPCA [http://en.wikipedia.org/wiki/Principal_component_analysis]уменьшение 900 габаритов до 10. Пока что у меня есть: covariancex = cov(labels); [V, d] = eigs(covariancex, 40); pcatrain = (trainingData - ...