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

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

Вычисление R-матрицы в QR-разложении с поворотом столбца в R

У меня есть матрица A который можно получить по следующему коду: b = matrix(c(20, 33, 10, 12, 14, 22, 34, 55, 11, 40, 0, 0, 0, 0, 0, 33,40, 66, 78, 90, 11, 45, 32, 55, 65), nrow = 5, ncol= 5) A =crossprod(b) > A [,1] [,2] [,3] [,4] [,5] [1,] 1929…
24 янв '14 в 05:30
1 ответ

Множественный регрессионный анализ в R с использованием QR-разложения

Я пытаюсь написать функцию для решения множественной регрессии с помощью QR-разложения. Вход: вектор y и матрица X; выход: б, е, р ^2. До сих пор я получил это и ужасно застрял; Я думаю, я сделал все слишком сложным: QR.regression <- function(y, …
1 ответ

QR-разложение и разложение Холецкого в R

Недавно я читал о том, как R-матрица QR-разложения может быть вычислена с использованием разложения Choleski. Отношение это: R = разложение Холецкого (A^TA) Пример: > A=matrix(c(1,2,3,2,3,5,1,3,2), nrow=3) > A [,1] [,2] [,3] [1,] 1 2 1 [2,] 2 …
1 ответ

QR-разложение в MatLab

У меня есть следующая проблема. Моя задача - подогнать полином к данным. Я хочу имплантировать QR-алгоритм, используя процесс ортогонализации Грама-Шимдта. В эту функцию встроено: function [ Q,R ] = QRDec( A ) n = length(A(1,:)); for i=1:n Q(:,i) = …
14 окт '13 в 13:34
1 ответ

Как наиболее эффективно использовать QR-разложение в Юлии?

Избегание выделения массива хорошо для производительности. Тем не менее, мне еще предстоит понять, как наиболее эффективно можно выполнить QR-разложение матрицы A, (примечание: необходимы матрицы Q и R) Просто используя Q, R = qr(A) это, вероятно, н…
1 ответ

Имеет ли cuSolverDN или другая библиотека CUDA пакетную версию QR-разложения для плотных матриц для решения A*x = b?

Я пытаюсь решить A*x = b, где A имеет комплексные значения и является плотным. Я использовал метод cusolverDnCgeqrf() из библиотеки cuSolverDN, чтобы выполнить QR-разложение для одного линейного набора уравнений. Однако я хочу сделать это несколько …
14 июл '17 в 00:10
0 ответов

Что мы можем сделать с тензорной факторизацией?

У меня три вопроса о тензорной факторизации. как обстоят дела (или применение) для тензорной факторизации (разложения)? насколько вероятно, что это станет основной технологией в будущем? Как вы используете его?
1 ответ

Неожиданный конец ввода, R программирование

В настоящее время я создаю приложение на R для вычисления разложения QR-матрицы, разложения неотрицательной матрицы QR и вычисления ICA. На данный момент я работаю над первым заданием. Я получаю следующую ошибку: источник ("trial.R") Ошибка в источн…
06 июл '17 в 15:10
0 ответов

Изменить / уменьшить собственную перестановочную матрицу

У меня проблемы с решением того, что я считаю довольно простой проблемой. Основная проблема в том, что я хочу изменить Eigen PermutationMatrix, но я не знаю как. Я делаю QR-разложение какой-то матрицы X используя библиотеку C++ Eigen. Я делаю это на…
18 ноя '17 в 04:53
1 ответ

Использование QR-разложения для решения наименьших квадратов в Matlab

Я использую Matlab для оценки регрессионной модели с обычными наименьшими квадратами (OLS). Модель y = xB, где x очень разреженная матрица с размерностью 500000 x 2500, Я использую разложение QR: [C,R] = qr(x,y,0) а затем оценивая b с b = R\C Мой во…
22 апр '14 в 17:07
2 ответа

Различные результаты разложения QR с NumPy и CULA

Я выполняю QR-декомпозицию двумя различными способами: используя стандартный метод numpy и используя функцию GEQRF LAPACK, реализованную в библиотеке CULA. Вот простой пример в Python (PyCULA используется для доступа к CULA): from PyCULA.cula import…
24 апр '14 в 13:22
1 ответ

Невозможно получить коэффициенты линейной регрессии в R после успешного поиска вопросов и ответов через домохозяина

Я вручную пытаюсь вычислить коэффициенты регрессии, а не использую какие-либо значения по умолчанию для данных http://people.sc.fsu.edu/~jburkardt/datasets/regression/x31.txt Вот мой код, который правильно производит Q&R;, удовлетворяющий A=QR. Но я…
0 ответов

Обратная матрица Eigen SparseQR не является точной, как в Matlab

В моем алгоритме я использую разреженную матрицу обратной операции и решаю ее, используя метод A*x=b, используя метод декомпозиции QR. На Matlab операция QR работает нормально. Однако когда я попытался преобразовать код в C++ с помощью библиотеки Ei…
1 ответ

Функция qr в R и Matlab

У меня есть вопрос о преобразовании функции Matlab в R, и я надеялся, что кто-то может помочь. Стандартное QR-разложение, используемое как в matlab, так и в R, называется qr(). Насколько я понимаю, стандартный способ выполнения разложения qr на обои…
01 дек '17 в 15:31
1 ответ

QR-разложение в TensorFlow

Я видел, что существуют методы для разложения Холецкого и решения линейных систем с использованием метода QR в TensorFlow, однако я не могу найти способ реализовать разложение QR в TensorFlow. Как вы выполняете QR-разложение в TensorFlow?
28 сен '16 в 19:48
1 ответ

Функция ортонормирования Лапака для прямоугольной матрицы

Мне было интересно, есть ли в Лапаке функция для ортонормирования столбцов очень высокой и тощей матрицы. Аналогичный предыдущий вопрос задавал этот вопрос, предположительно в контексте квадратной матрицы. Моя настройка следующая: у меня есть матриц…
1 ответ

Поиск корней с помощью сопутствующей матрицы

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

Как получить Q из результатов QR-факторизации?

DGEQRF и SGEQRF из LAPACK возвращают часть Q факторизации QR в упакованном формате. Распаковка вроде требует O(k^3) шаги (k продуктов низкого ранга), и, кажется, не очень просто. Плюс числовая стабильность ведения k Последовательное умножение мне н…
2 ответа

Найти вращение и перекос преобразования матрицы

У меня есть следующая матрица преобразования в CSS // rotate the element 60deg element.style.transform = "matrix(0.5,0.866025,-0.866025,0.5,0,0)" И я могу найти вращение, используя это... // where a = [0.710138,0.502055,-0.57735,1,0,0] var rotation …
24 фев '11 в 16:02
0 ответов

Реализация псевдокода для обновления QR - какая часть моего кода неверна?

Я реализую некоторый псевдокод для обновления QR-разложения - я хочу удалить строку из матрицы A и обновить Q а также R соответственно. Мои результаты отличаются от результатов, полученных Scipy's qr_delete() хотя (см. строку 21 в прилагаемом фрагме…
31 дек '18 в 17:56