|
|||||||
АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция |
Метод вращения (Якоби)Основан на вращении системы координат в n -мерном пространстве. Матрица вращения имеет вид: , при этом выполняется С2+S2=1,C=cos α, S=sin α. Цель метода – найти такие вращения, которые обращают определенные элементы исходной матрицы в ноль. Замечание. Любую симметричную матрицу можно привести к трехдиагональному виду путем цепочки преобразования подобия с матрицей вида Pij, k=1,…,n. Указанный процесс можно применить и к несимметричной матрице, при этом получится почти треугольная матрица. Процесс преобразования исходной матрицы путем элементарного вращения на любом k –шаге можно представить в виде рекуррентного соотношения , k=1,2,… 1 шаг: Вычислим . При этом отличными являются элементы, стоящие в i и j -столбцах , k=1,2,…,n. 2 шаг: отличными являются элементы, стоящие в i и j - строчках , k=1,2,…,n.
Неопределенными параметрами являются C и S, но они связанны соотношением С2+S2=1. Недостающее одно уравнение получаем из условия обращения в ноль некоторого элемента матрицы А, в общем случае алгоритм представляет последовательное обращение всех ненулевых элементов, лежащих вне трех диагоналей исходной матрицы.
3 шаг: После соответствующих алгебраических преобразований, получим , для определенности считают . Можно показать, что в построенном процессе все внедиагональные элементы становятся сколь угодно малыми при достаточном числе подобных итераций. Собственными векторами являются столбцы матрицы .
Пример. Симметричную трехдиагональную матрицу свести путем подобных преобразований к диагональной. При этом на диагонале получим собственные значения. Обнулим a12. Исходная матрица А0 = .
Матрица подобия .
При этом элементы a12 =a21 = 2SC-C2+S2-2CS=0, C2=S2. Учитывая основное тригонометрическое соотношение, получаем , угол равен α = π/4.
Поиск по сайту: |
Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.004 сек.) |