Навигация
ДРУГИЕ МЕТОДЫ ВЫЧИСЛЕНИЯ СОБСТВЕННЫХ ЗНАЧЕНИЙ
30010
знаков
20
таблиц
3
изображения
6. ДРУГИЕ МЕТОДЫ ВЫЧИСЛЕНИЯ СОБСТВЕННЫХ ЗНАЧЕНИЙ
В этом разделе мы рассмотрим два метода определения собственных значений, имеющие большое практическое значение. Оба разработаны в последние 20 лет и наиболее эффективны в тех случаях, когда требуется найти все собственные значения произвольной матрицы действительных или комплексных чисел. В обоих используются преобразования, позволяющие получить последовательность подобных матриц, сходящуюся к матрице блочной треугольной формы:
| X1 | * | … | … | * | * | * | |
| x2 | * | … | … | * | * | * | |
| x3 | … | … | * | * | * | ||
| … | … | * | * | * | |||
| … | * | * | * | ||||
| … | * | * | |||||
| 0 | … | * | |||||
| * |
где блоки Хm, представляют собой матрицы размерности 2 х 2, расположенные на главной диагонали. Собственные значения блоков Хm, являются в то же время собственными значениями исходной матрицы размерности п x п. Такая форма удобна, так как детерминант второго порядка блоков Хm позволяет определять комплексные собственные значения, не вводя комплексных элементов в окончательную матрицу. Если все собственные значения исходной матрицы действительные, то в окончательном виде она будет треугольной, причем собственные значения будут расположены на диагонали.
Метод LR
Этот метод первоначально был разработан Рутисхаузером в 1958 г. Метод основан на представлении матрицы A в виде произведения
А = LR,
где L — левая треугольная матрица с единичными диагональными элементами, а R — правая треугольная. Применяя преобразование подобия L-1 A R, видим, что,
A2 = L-1 A R = L-1 (RL)L = R L.
Следовательно,
Am-1 = L m-1 Rm-1,
Am = R m-1 Lm-1.
Этот процесс повторяется до тех пор, пока Ls не превратится в единичную матрицу Е, а Rs не приобретет квазидиагональную форму. Хотя этот метод очень удобен, он не всегда устойчив. Поэтому предпочтение часто отдают другому методу.
Метод QR
Метод QR. предложен Фрэнсисом в 1961 г. Соответствующий ему алгоритм определяется соотношением
Am = Q m Rm.
где Q m — ортогональная матрица, а Rm — верхняя треугольная матрица. При использовании метода последовательно получаем
Am+1 = QmTAm Qm = QmTQm Rm Qm = Rm Qm.
В пределе последовательность матриц А стремится к квазидиагональной форме. Этот метод сложнее предыдущего и требует больших затрат машинного времени. Однако его устойчивость,обусловленная использованием ортогональных преобразующих матриц, обеспечила ему прочную репутацию лучшего метода решения задач самой общей формы.
Пример 3
Пусть требуется найти все собственные значения произвольной матрицы размерности 6 x 6
| 2,3 | 4,3 | 5,6 | 3,2 | 1,4 | 2,2 |
| 1,4 | 2,4 | 5,7 | 8,4 | 3,4 | 5,2 |
| 2,5 | 6,5 | 4,2 | 7,1 | 4,7 | 9,3 |
| 3,8 | 5,7 | 2,9 | 1,6 | 2,5 | 7,9 |
| 2,4 | 5,4 | 3,7 | 6,2 | 3,9 | 1,8 |
| 1,8 | 1,7 | 3,9 | 4,6 | 5,7 | 5,9 |
Сделаем это в два приема, приведя сначала матрицу с помощью преобразования подобия к виду Гсссенберга, затем с помощью разновидности метода QR найдем собственные значения. В приведенной ниже программе использованы две подпрограммы из пакета программ для научных исследований фирмы IВМ. Подпрограмма НSВС преобразует матрицу размерности 6 x 6 к форме Гессенберга, а подпрограмма АТЕIG позволяет найти собственные значения.
{**********************************************************************}
Программа определение всех собственных значений произвольной матрицы размерности 6х5. Используются подпрограммы НSВС и АТЕIG из пакета программ для научных исследований фирмы IBM
{**********************************************************************}
DIMENSION A(6,6),RR(6),RI(6),IANA(6)
READ(5,100)((A(I,J),J=1,6),I=1,6)
WRITE(6,104)
104 FORMAT(///lX,’THE ORIGINAL MATRIX IS AS FOLLOWS’)
WRITE(6,103)
103 FORMAT(1X,65(-'--'))
WRITE(6,101)((A(I,J),J=1,6),I=1,6)
WRITE(6,103)
101FORMAT(6(1X,F10.5))
100 FORMAT(6F10.5)
CALL HSBG(6,A,6)
WRITE(6,105)
105 FORMAT(///1X,'THE MATRIX W HESSENBUR5 FORM IS') WRITE(6,103)
WRITE(6,101)((A(I,J),J=1,6),I=1,6)
WRITE(6,103)
CALL ATEIG(6,A,RR,RI,IANA,6)
WRITE(6,106)
106FORHAT(///1X,'THE EIGENVALUES ARE AS FOLLOUS')
WRITE(6,107)
Похожие работы
... 0135 0.7866 -0.5989 График значений собственного числа заданной матрицы за время итерационного процесса/ График значений собственного вектора, соответствующего собственному числу Рассмотрим другие примеры: Исходные данные: yn=[1,1,1]; L1= 0.01 edop=0.00001; a=[1 1 1; 2 3 4; 0 4 0]; Найдем собственные значения исходной матрицы, используя функцию eig. Получим L1= 6.2085 ...
... решения от численных методов расчёта. Для определения корней уравнения не требуется знания теорий групп Абеля, Галуа, Ли и пр. и применения специальной математической терминологии: колец, полей, идеалов, изоморфизмов и т.д. Для решения алгебраического уравнения n - ой степени нужно только умение решать квадратные уравнения и извлекать корни из комплексного числа. Корни могут быть определены с ...
... степень изношенности дисков, налипание к валу инородных предметов и так далее, то обратная задача состоит в диагностировании характеристик вала с дисками по собственным частотам колебаний вала. Известно, что изменения указанных значений характеристик вала проявляются в изменениях значений собственных частот его колебаний, что в свою очередь может привести к ненужным вибрациям, увеличению шума и т. ...
... собственных значений сводится к проблеме раскрытия определителя по степеням и последующему решению алгебраического уравнения m- й степени. Определитель называется характеристическим (или вековым ) определителем, а уравнение (2) называется характеристическим (или вековым ) уравнением. Различают полную проблему собственных значений, когда необходимо отыскать все собственные значения ...
0 комментариев