|
Текст подпрограммы и версий ( Фортран ) aeg2c.zip , aeg2p.zip |
Тексты тестовых примеров ( Фортран ) taeg2c.zip , taeg2p.zip |
|
Текст подпрограммы и версий ( Си ) aeg2c_c.zip |
Тексты тестовых примеров ( Си ) taeg2c_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) aeg2c_p.zip , aeg2z_p.zip |
Тексты тестовых примеров ( Паскаль ) taeg2c_p.zip , taeg2z_p.zip |
Вычисление всех собственных значений комплексной матрицы QR - алгоpитмом со сдвигом.
Исходная матрица приводится унитарными преобразованиями к верхней форме Хессенберга, которая затем используется для вычисления собственных значений.
Дж.Х.Уилкинсон, Алгебраическая проблема собственных значений, "Наука", М., 1970.
SUBROUTINE AEG2C ( N, AR, AI, EVR, EVI, IERR)
Параметры
| N - | порядок исходной матрицы (тип: целый); |
| AR, AI - | вещественные двумерные массивы размерности N на N, содержащие соответственно действительную и мнимую части исходной матрицы; |
|
EVR - EVI | вещественные одномерные массивы размерности N, содержащие соответственно действительную и мнимую части вычисленных собственных значений; |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; значение IЕRR полагается равным номеру собственного значения, для вычисления которого потребовалось более 30 итераций, при этом собственные значения с индексами IЕRR + 1, IЕRR + 2, ..., N вычислены правильно. |
Версии
| AEG2P - | вычисление всех собственных значений комплексной матрицы с повышенной точностью. Массивы АR, АI, ЕVR, ЕVI имеют тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы
| UTAE10 - | подпрограмма выдачи диагностических сообщений при работе подпрограмм АЕG2С и АЕG2Р. |
Замечания по использованию
| Подпрограммы АЕG2С(P) не сохраняют исходную матрицу. | |
| Подпрограмма АЕG2C использует служебную подпрограмму AA02C. Подпрограмма АЕG2P использует служебные подпрограммы AA01P, AA02P, AA03P. |
DIMENSION AR(4, 4), AI(4, 4), EVR(4), EVI(4)
DATA AR /1., 0.42, 0.54, 0.66, 0.42, 1., 0.32, 0.44, 0.54, 0.32,
* 1., 0.22, 0.66, 0.44, 0.22, 1./
DATA AI /16*0./
CALL AEG2C (4, AR, AI, EVR, EVI, IERR)
Результаты:
EVR(1) = 2.3227488 ,
EVR(2) = 0.2422607 ,
EVR(3) = 0.6382838 ,
EVR(4) = 0.79670669 ,
EVI = (0., 0., 0., 0.) ,
IERR = 0