Экспериментальное исследование численных методов решения систем линейных алгебраических уравнений

Санкт-Петербургский Государственный

Политехнический Университет

Отчет о лабораторной работе № 1.

«ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ  ЧИСЛЕННЫХ  МЕТОДОВ  РЕШЕНИЯ  СИСТЕМ ЛИНЕЙНЫХ  АЛГЕБРАИЧЕСКИХ  УРАВНЕНИЙ»

                                                                                                                                   Выполнил:

                                                                                                                студент группы 2091/1

                                                                                                                         Журавлев Иван

                                                                                                                                       Преподаватель

Анисимов Андрей

Санкт-Петербург

2005 г.

1 Сравнение естественного и стандартного числа обусловленности матрицы

После проведения ряда численных экспериментов была составлена следующая таблица:

Вывод: Проанализировав ее,  можно сделать некоторые предварительные выводы:

 Относительное различие стандартного и естественного числа обусловленности  невелико—всегда в пределах одного порядка. Чем хуже обусловлена матрица(чем выше µ1 иµ2) тем меньше они отличаются. Для матриц с хорошей обусловленностью точность оценки стандартного числа с помощью процедуры DECOMP  зависит от того, каким алгоритмом генерировалась матрица. Например, µ2 матрицы, сгенерированных по алгоритму 2 оценивались достаточно точно. А для матриц, сгенерированных по алгоритму1 оценка cond отличалась от µ2 иногда даже на один порядок.  Естественно, что  на практике стараются сводить поставленные задачи к системам с хорошо обусловленными  матрицами (). Так как в этом случае чувствительность решения к погрешности исходных данных мала. Если же невозможно обойти использование плохо обусловленной матрицы, то необходимо обеспечить высокую точность исходных данных, в противном случае мы рискуем получить  результаты, не  имеющие к действительности никакого отношения.

2 Исследование точности решения как функции от порядка матрицы при постоянном числе обусловленности(хорошо обусловленные матрицы)

После проведения многочисленных испытаний была получена масса материала. Однако  в ходе анализа были оставлены как наиболее показательные следующие данные:

Из таблицы видно, что при практически постоянном числе  обусловленности ошибка  возрастает при возрастании порядка.

Вывод: Такой  характер зависимости можно объяснить тем, что в методе Гаусса при возрастании порядка системы возрастает число операций сложения и деления и, как следствие,  число округлений. 

3 Исследование точности решения как функции от порядка матрицы при постоянном числе обусловленности (плохо обусловленные матрицы).

Были проведены аналогичные Пункту№2 исследования и получены следующие результаты:

Вывод: Нетрудно заметить, что возрастание относительной погрешности при росте порядка системы имеет более быстрый характер. Что, собственно, и предсказывается теорией.

Необходимо отметить, что оценка погрешности с помощью процедуры DECOMP дает сильно завышенное значение. Конечно же, в случае плохо обусловленных матриц применение данных алгоритмов решения нецелесообразно, так как погрешность достигает положительных степеней. То есть необходимо использовать какие-то иные алгоритмы. Какие?

4  Исследование зависимости точности решения от числа обусловленности при постоянном порядке системы.

Были проведены эксперименты с матрицей порядка 7.

Результаты измерений могут быть представлены следующими таблицей и графиками. График зависимости ошибки от обусловленности?

Вывод: Как и следовало ожидать, относительная ошибка возрастает с ростом числа обусловленности. При достаточно высоких числах обусловленности(порядка 10^15) использование алгоритма Гаусса становиться нецелесообразным, так как погрешность имеет уже положительные степени. Какую либо зависимость вектора невязки от числа обусловленности установить не удалось.Почему?

5.Исследование возможности улучшения обусловленности задачи путем внесения в матрицу малого случайного возмущения.