Вопросы с тегом «eigensystem»

РЕДАКТИРОВАТЬ: я проверяю, если какие-либо собственные значения имеют величину один или больше. Мне нужно найти наибольшее абсолютное собственное значение большой разреженной несимметричной матрицы. Я использовал eigen()функцию R , которая использует алгоритм QR из EISPACK или LAPACK, чтобы найти все собственные значения, а затем я использую, abs()чтобы получить абсолютные значения. Однако …

27 linear-algebra  algorithms  eigensystem  sparse-matrix 

Я хочу вычислить спектр ( все собственные значения) большой разреженной матрицы (сотни тысяч строк). Это трудно. Я готов согласиться на приближение. Существуют ли методы приближения для этого? Хотя я надеюсь получить общий ответ на этот вопрос, я также был бы удовлетворен ответом в следующем конкретном случае. Моя матрица - нормализованный …

14 linear-algebra  eigensystem  graph-theory  approximation 

Давайте начнем с проблемы формы (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 с набором заданных граничных условий ( Дирихле , Неймана , Робина , Периодического , Блох-Периодического ). Это соответствует нахождению собственных значений и собственных векторов для некоторого оператора LL\mathcal{L} при некоторой геометрии и граничных условиях. Подобную проблему можно получить, например, в акустике, …

13 pde  finite-element  eigensystem  verification 

Я должен решить обобщенные задачи на собственные значения где A и B являются трехдиагональными, B симметрично положительно определенным и вещественным, но A является только комплексно симметричным (не определенным или эрмитовым). Кроме того, мне нужно полное собственное разложение. В настоящее время я просто называю обобщенный eigensolver Лапака, но мне интересно, есть …

13 linear-algebra  eigensystem 

У меня есть программа, которая вычисляет наибольшее собственное значение из многих вещественных симметричных матриц 50x50, выполняя разложения по сингулярным числам для всех из них. SVD является узким местом в программе. Существуют ли алгоритмы, которые намного быстрее находят наибольшее собственное значение, или оптимизация этой части не даст большой отдачи от инвестиций?

13 linear-algebra  eigensystem 

Я пытаюсь выяснить, существует ли более быстрый способ вычисления всех собственных значений и собственных векторов очень большой и разреженной матрицы смежности, чем использование scipy.sparse.linalg.eigsh. Насколько я знаю, этот метод использует только разреженность и атрибуты симметрии матрицы. Матрица смежности также является двоичной, что заставляет меня думать, что есть более быстрый способ …

12 linear-algebra  python  performance  eigensystem  scipy 

Предположим, что A ∈ Rn × nA∈Rn×nA\in\mathbb{R}^{n\times n} - симметричная положительно определенная матрица. AAA достаточно большой, чтобы решить дорого A x = bAx=bAx=b. Существует ли итерационный алгоритм для наименьшего собственного значения AAA , который не включает инвертирование AAA в каждой итерации? То есть мне нужно было бы использовать итерационный алгоритм, …

11 linear-algebra  eigensystem  eigenvalues  iterative-method 

Я делаю диагонализацию Ланцоша большой разреженной матрицы (~ 2 миллиона элементов). Почти все шаги в алгоритме Ланцкоса выполняются параллельно на графическом процессоре, за исключением диагонализации матрицы Ланцоша для проверки сходимости. Для этого я использовал алгоритм TQLI из Numeric Recipes. Существуют ли методы нахождения собственной системы трехдиагональной матрицы, которые параллельны или …

11 linear-algebra  algorithms  eigensystem 

Предположим следующую матрицу дается [ 0,500 - 0,333 - 0,167 - 0,500 0,667 - 0,167 - 0,500 - 0,333 0,833 ] с транспонированной T . Продукт A T A = G дает [ 0,750 - 0,334 - 0,417 - 0,334 0,667 - 0,333 - 0,417 - 0,333 0,750 ] ,AAA⎡⎣⎢0.500−0.500−0.500−0.3330.667−0.333−0.167−0.1670.833⎤⎦⎥[0.500−0.333−0.167−0.5000.667−0.167−0.500−0.3330.833] …

11 linear-algebra  matrix  eigensystem 

Каждая действительная матрица можно свести к вещественной форме Шура Т = U T A U с помощью ортогонального преобразования similiary U . Здесь матрица T имеет квазитреугольную форму с 1 на 1 или 2 на 2 блока на главной диагонали. Каждый 1 по 1 блок соответствует реальному собственному значению А …

10 eigenvalues  eigensystem  lapack  numerics  scalapack 

Я пытаюсь диагонализировать некоторые плотные, плохо обусловленные матрицы. В машинной точности результаты являются неточными (возвращая отрицательные собственные значения, собственные векторы не имеют ожидаемой симметрии). Я переключился на функцию Eigensystem [] Mathematica, чтобы использовать преимущества произвольной точности, но вычисления очень медленные. Я открыт для любого количества решений. Существуют ли пакеты / …

10 linear-algebra  parallel-computing  eigensystem  precision  dense-matrix 

У меня плотная вещественная симметричная квадратная матрица. Размер составляет около 1000x1000. Мне нужно вычислить первый главный компонент и подумать, каким может быть лучший алгоритм для этого. Похоже, что MATLAB использует алгоритмы Арнольди / Ланцоша (для eigs). Но читая о них, я не уверен, есть ли у них какие-либо преимущества перед …

10 linear-algebra  algorithms  eigensystem 

У меня есть большая проблема с кубическим собственным значением: (A0+ λA1+λ2A2+λ3A3) х =0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Я мог бы решить это путем преобразования в линейную задачу на собственные значения, но это привело бы к системе 32323^2 как большой: ⎡⎣⎢-A0000я000я⎤⎦⎥⎡⎣⎢ИксYZ⎤⎦⎥= λ⎡⎣⎢A1я0A20яA300⎤⎦⎥⎡⎣⎢ИксYZ⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 …

9 c++  eigensystem  eigenvalues 

У меня есть проблема, когда мне нужно найти все положительные (так как собственное значение положительно) собственные пары небольшой (обычно меньше, чем 60x60) несимметричной матрицы. Я могу прекратить вычислять, когда собственное значение меньше определенного порога. Я знаю, что собственные значения реальны. Какие-нибудь предложения по алгоритмам, которые я мог бы использовать, чтобы …

9 performance  eigensystem  gpu 

Я знаю, что для решения симметричной задачи на собственные значения A x = λ xAx=λxAx = \lambda xмы можем использовать закон инерции Сильвестра, то есть число собственных значений AAA меньше, чем aaaравно количеству отрицательных элементов в где диагональная матрица получается из факторизации ЛПНП . Затем с помощью метода деления пополам …

9 linear-algebra  eigensystem