МатВектор

Command Palette

Search for a command to run...

🔢 Численные методы

Степенной метод

англ. Power method

Многократное умножение на матрицу вытягивает вектор вдоль главного направления: при |λ₁| > |λ₂| итерации сходятся к доминирующему собственному значению, а отношение Рэлея его вычисляет.

Деньги, разложенные по вкладам с разными ставками, со временем стягиваются к самому доходному: остальное обгоняется и тонет в шуме. Степенной метод делает это с вектором и матрицей. Берём ненулевой и считаем : умножение на матрицу каждый раз растягивает компоненту вдоль доминирующего направления сильнее остальных, и при направление вектора ползёт к собственному вектору с максимальным по модулю значением . Само значение достают отношением Рэлея ; скорость сходимости — с коэффициентом за итерацию.

Числовой прогон: матрица со строками и имеет собственные значения и . Старт : отношения Рэлея , затем , , , — примесь второго направления делится на три за итерацию, потому что . Ограничения видны отсюда же: при близких по модулю собственных значениях сходимость вязнет, а при комплексной паре замирает вовсе; лечится сдвигом спектра . Место метода среди прочих — в уроке собственные значения: численные методы, точная теория — в собственных значениях, объекты — в терминах собственный вектор и диагонализация матрицы.

отношение Рэлея для матрицы со строками (2, 1) и (1, 2), где , : старт (1, 0) даёт — гашение втрое за итерацию

Частые вопросы

Что делать, если доминирующих значений несколько?

Если два наибольших по модулю значения равны — кратное или комплексная пара , — классическая итерация не выберет направление: предел зависит от старта или осциллирует. Стандартное лекарство — сдвиг : переход к меняет отношения модулей и делает доминирующим нужное значение. Для симметричных матриц работают блочные версии, вылавливающие несколько направлений сразу, — из них вырастают методы подпространств Крылова.

Зачем степенной метод, если есть точные разложения?

Потому что для огромных разреженных матриц — графов, сеток уравнений математической физики — точное разложение и не нужно, и невозможно: оно забивает память плотной структурой. Степенному методу доступна одна операция — умножение на вектор, а для разреженных матриц она стоит числа ненулевых элементов. Именно на этой операции построены методы Крылова, включая Ланцоша и GMRES; PageRank поисковиков — по сути степенной метод в масштабе миллиардов страниц.