Справочник функций

Ваш аккаунт

Войти через: 
Забыли пароль?
Регистрация
Информацию о новых материалах можно получать и без регистрации:

Почтовая рассылка

Подписчиков: -1
Последний выпуск: 19.06.2015

4. Квадратная матрица A называется ортогональной, если AT=A-1. Определить, явля-ется ли заданная матрица A(n,n) ортогональной.

78K
03 мая 2012 года
ninika
5 / / 03.05.2012
Квадратная матрица A называется ортогональной, если AT=A-1. Определить, является ли заданная матрица A(n,n) ортогональной.
82K
03 мая 2012 года
suroviy
3 / / 03.05.2012
Не понял вопроса: можно действовать по определению, проверив близость матрицы transpose(A) * A к единичной. Программируется это просто (нужно только оценивать близость к единичной матрице). Но если у вас матрицы индуцируются конкретным оператором, лучше убедиться в этом теоретически.
78K
05 мая 2012 года
ninika
5 / / 03.05.2012
не совсем поняла)))
82K
05 мая 2012 года
suroviy
3 / / 03.05.2012
В точной арифметике, если матрица ортогональная, то при умножении ее на транспонированную, получится единичная матрица. Ну например, исходная матрица:
A = {1/sqrt(2), 1/sqrt(2) ; -1/sqrt(2), 1/sqrt(2) } -- запись по строкам
Транспонированная от нее:
transpose(A) = {1/sqrt(2), -1/sqrt(2) ; 1/sqrt(2), 1/sqrt(2) }

Умножим transpose(A) на A, получим единичную матричку:
E={1, 0 ; 0, 1}

Проблема в том, что в вещественной арифметике (например, в double), получится не E, а что-то похожее на
E_={1.0000000000001, 1.E-12 ; 1.E-12, 1.0000000000001}


Т.е. в вещественной арифметике невозможно "точно" установить, что матрица ортогональнана, но в данном случае можно проверять, чтобы диагональные элементы отличались от 1 менее чем на 1.E-10, а внедиагональные были не более 1.E-10 (все сравнения по модулю).

Как правило, матрицы это не сферические кони в вакууме, а образы операторов в конкретном базисе, поэтому лучше теоретически доказывать интересные свойства данных матриц (ортогональность, СПО, ...).

Наверное, если будут еще вопросы, лучше в личку писать.

Знаете кого-то, кто может ответить? Поделитесь с ним ссылкой.

Ваш ответ

Реклама на сайте | Обмен ссылками | Ссылки | Экспорт (RSS) | Контакты
Добавить статью | Добавить исходник | Добавить хостинг-провайдера | Добавить сайт в каталог