Описание тега matrix-decomposition

Описание тега Вопросы с тегом
В математической дисциплине линейной алгебры разложение матриц или факторизация матриц - это факторизация матрицы в произведение матриц.
1 ответ

Как сделать разложение ЛПНП, когда все, что у меня есть, это LU-решатель, который всегда применяет поворот?

Я надеюсь, что мне просто не хватает простого трюка с матричной арифметикой, но проблема, с которой я сталкиваюсь, заключается в том, что все, к чему у меня есть доступ, это LU-решатель (Matlab LU* или SuperLU), и мне нужно разложение LDL симметричн…
12 окт '17 в 12:03
0 ответов

4x4 преобразование матрицы преобразования

Я пытаюсь разложить однородную матрицу преобразования координат, которая является результатом этих матриц преобразований (примененных в порядке их записи): неравномерный Масштаб -> Повернуть -> Перевести У меня нет проблем с извлечением координат пе…
12 апр '14 в 15:18
1 ответ

Почему декомпозиция LU с использованием Parallel.For не работает?

Я пытаюсь решить разложение LU с помощью алгоритма Дулиттла - согласно этому документу. Без распараллеливания код работает нормально. Тем не менее, я хотел бы, чтобы этот код выполнялся параллельно - пытаясь создать параллельный внешний и два внутре…
01 дек '17 в 21:36
0 ответов

Написание функции R, которая выводит оценки максимального правдоподобия?

Как бы я написал функцию R, которая берет вектор отклика и ковариатную матрицу X и выводит вектор оценок максимального правдоподобия коэффициентов β, где μ = Xβ и E[Y] = µ, где Y - экспоненциально распределенный случайный вектор? Я могу использовать…
21 фев '16 в 23:40
2 ответа

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

Я ищу матричную операцию формы: B = M*A*N где A некоторая общая квадратная матрица и M а также N это матрицы, которые я хочу найти. Такие, что колонны B являются диагоналями A, В первом столбце главная диагональ, во втором диагональ смещена на 1 от …
04 янв '17 в 12:50
1 ответ

Решение для Lx=b и Px=b, когда A=LLt

Я разлагаю разреженную матрицу SPD A, используя Eigen. Это будет либо разложение LLt, либо LDLt (Холецкий), поэтому мы можем предположить, что матрица будет разложена как A = P-1 LDLt P где P - матрица перестановок, L - треугольная нижняя и D диагон…
21 мар '18 в 17:17
1 ответ

Ruby Matrix:: Ошибка разложения по собственным значениям

Я посмотрел на матрицу собственных векторов данной матрицы, но когда я пытаюсь инвертировать ее, у меня возникает ошибка в eigenvector_matrix_inv(), require 'matrix' m = Matrix[ [0.5703125, 1.8369140625, 0.0, 0.0], [-0.6875, -0.4609375, 0.0, 0.0], […
21 сен '12 в 00:22
0 ответов

Как обновить результаты функции разрешения в исходных матрицах L и U?

Я могу найти каждое уравнение для решения, но я не могу обновить результат каждого вычисления в исходной матрице. Я попытался задать набор элементов в матрице, заменив переменные u фактическим результатом из функции решения. Но тогда он выдает следу…
12 фев '19 в 08:36
1 ответ

Как работает приближение равенства Matlab?

В следующем коде я следовал процедуре для создания случайной положительно определенной матрицы P. Сначала я создал разложение по единственному значению [U,S,V] случайного массива A, и я пытаюсь проверить, что на самом деле U '* U == U * U' = I (где …
01 ноя '15 в 10:24
0 ответов

Как повторно использовать разложенную матрицу LU для решения Ax = b в JBLAS

Я начал работать с JBLAS, но столкнулся с проблемой, doubleMatrix x = Solve.solve(A,b); Это дает результаты просто отлично. Но если я захочу сделать перерасчет, чтобы найти Ax = b1, то он снова разложится, чтобы получить A = LU или, возможно, A = QR…
16 июл '15 в 09:20
0 ответов

Как наложить ограничения на разреженность в решателе компакт-дисков библиотеки Склеарн НМФ?

Недавно я обнаружил, что библиотека Sklearn NMF устарела в своем pg solver и теперь использует свой cd solver. С помощью cd solver я не думаю, что смогу применить ограничения разреженности. Или, может быть, уровень L1 неявно представляет разреженнос…
28 май '16 в 22:08
1 ответ

Функция chol() в R продолжает возвращать верхнюю треугольную (мне нужна нижняя треугольная)

Я пытаюсь получить разложение по нижнему треугольнику Холецского следующей матрицы в R, используя chol() функция. Тем не менее, он продолжает возвращать верхнюю треугольную декомпозицию, и я не могу найти способ получить нижнюю треугольную декомпози…
09 ноя '16 в 06:14
1 ответ

Ошибка атрибута X_transformed_fit_: AttributeError: у объекта "KernelPCA" нет атрибута "X_transformed_fit_"

Я пытаюсь выяснить, какие функции в моем наборе данных влияют на основные компоненты, и пытаюсь наблюдать, как мои данные вписываются в мой алгоритм PCA ядра. Я пытался использовать атрибут X_transformed_fit_, который существует в документальном фил…
24 май '18 в 11:47
0 ответов

Кальман Разложение в Python

В настоящее время я прохожу курс современного управления и хочу выполнить разложение Кальмана ( https://en.wikipedia.org/wiki/Kalman_decomposition) в системе, указанной в задании. Назначение предлагает использовать функцию minreal() MatLab ( https:/…
05 мар '18 в 22:28
0 ответов

Eigen3 Sparse Solver не копируемый

Я работаю над числовым кодом и хочу оценить, как различаются Sparse и Dense Matrix-LU (а также и другие, более поздние) для варианта использования кода. Eigens Dense Decomposition Objects могут быть копируемыми, и это используется для их кэширования…
23 ноя '16 в 12:37
1 ответ

Как вывести 2 или более массивов в функции фортрана?

Я пишу программу, которая вычисляет декомпозицию матрицы LU с частичным поворотом, и я хотел бы, чтобы функция выводила несколько (2 или 3) матриц без запуска программы несколько раз для вывода каждой из них по отдельности, что является пустой трато…
12 сен '16 в 08:38
3 ответа

LU разложение квадратной матрицы устранения Matlab Гаусса

Я пытаюсь создать программу, которая принимает квадратную (n-на-n) матрицу в качестве входных данных, и, если она обратима, LU будет разлагать матрицу с использованием гауссовского исключения. Вот моя проблема: в классе мы узнали, что лучше менять с…
21 июл '13 в 19:46
2 ответа

Лучший способ вычислить усеченное разложение по сингулярным значениям в Java

Я хочу сравнить лучшие 2 или 3 библиотеки для вычисления усеченного разложения по сингулярным значениям (SVD), то есть SVD, где хранятся только k самых больших значений в единственном числе. Кроме того, у меня есть эти ограничения: Это должна быть б…
13 ноя '13 в 15:01
2 ответа

Библиотеки для параллельного распределённого холесского разложения в c/ C++ в среде mpi?

Какие библиотеки доступны для параллельного распределенного холесского разложения плотных матриц в C/C++ в среде mpi? Я нашел библиотеку ScaLAPACK, и это может быть решением, которое я ищу. Кажется, что это немного неудобно для вызова, так как нужно…
08 июн '13 в 08:05
1 ответ

MATLAB LU Разложение Частичное вращение

Я пытаюсь работать с моей декомпозицией lu, в значительной степени основанной на декомпозиции LU с частичным поворотом Matlab function [L,U,P] = lup(A) n = length(A); L = eye(n); U = zeros(n); P = eye(n); for k=1:n-1 % find the entry in the left col…
30 мар '15 в 23:35