WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

загрузка...
   Добро пожаловать!

Pages:     || 2 |

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

На правах рукописи

БОРЗЫХ Алексей Николаевич

О НОВЫХ МЕТОДАХ РЕШЕНИЯ ЧАСТИЧНОЙ

ПРОБЛЕМЫ СОБСТВЕННЫХ ЗНАЧЕНИЙ

01.01.07 — вычислительная математика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата физико-математических наук

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

2008

Работа выполнена на кафедре вычислительной математики математико-механического факультета Санкт-Петербургского государственного университета.

Научный руководитель: кандидат физико-математических наук,

доцент САМОКИШ Борис Андреевич

Официальные оппоненты: доктор физико-математических наук, профессор ДЕМЬЯНОВИЧ Юрий Казимирович (Санкт-Петербургский государственный университет)

доктор физико-математических наук,

профессор ХАЗАНОВ Владимир

Борисович (Санкт-Петербургский государственный морской технический университет)

Ведущая организация: Московский государственный

университет им. М. В. Ломоносова

Защита состоится 24 декабря 2008 г. в 13 часов на заседании совета Д 212.232.49 по защите докторских и кандидатских диссертаций при Санкт-Петербургском государственном университете по адресу: 199048, Санкт-Петербург, Васильевский остров, 14 линия, д. 29, математико-механический факультет, ауд. 13.

С диссертацией можно ознакомиться в Научной библиотеке им. М. Горького Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская набережная, д. 7/9.

Автореферат разослан "____"______2008 г.

Ученый секретарь диссертационного совета,

доктор физико-математических наук, профессор А. А. Архипова

Общая характеристика работы

Актуальность темы диссертационной работы

Проблема собственных значений в различных ее постановках (полная или частичная, симметричная или несимметричная) является задачей вычислительной математики, интерес к которой не угасает уже много десятилетий. Задачи подобного рода, с одной стороны, достаточно часто возникают в разнообразных инженерных расчетах, приводя к матрицам, как правило, больших размерностей; с другой стороны, имеют кубическую зависимость объема вычислений от размера задачи, что требует существенного времени счета даже на современных быстродействующих ЭВМ. Универсального алгоритма, обеспечивающего эффективное решение задачи в любой ее постановке, не существует. Многие из алгоритмов, существующие ныне, проходили длительный путь от их первоначального «зарождения» до современного вида, возникшего благодаря многократным совершенствованиям и теоретическим исследованиям (например, QR-алгоритм). На основании вышесказанного можно утверждать, что появление любого нового алгоритма, решающего какую-либо постановку проблемы собственных значений, должно вызывать научный интерес.

Цель диссертационной работы

Целью диссертационной работы является разработка и исследование нового алгоритма, вычисляющего наибольшее собственное значение симметричной вещественной матрицы, а также алгоритма, вычисляющего наибольшее сингулярное число несимметричной вещественной матрицы; получение теоретических доказательств линейной сходимости разработанных алгоритмов; формирование вычислительного процесса, оптимального для реализации на ЭВМ; проведение сравнительных численных экспериментов с предложенными в диссертации и ранее известными алгоритмами.

Методы исследования

Линейная сходимость предложенных алгоритмов доказывается теоретически. Построение экономичной вычислительной схемы основывается на минимизации количества арифметических операций, выполняемых на каждом шаге алгоритма. Эмпирическое сравнение предложенных в диссертации и ранее известных алгоритмов проводится на основе их программной реализации и тестовых испытаний. Эксперименты выполняются на матрицах различной размерности, для различных относительных точностей расчета, для матриц с различной близостью старших собственных (сингулярных) чисел, для вычислительных процессов с разным коэффициентом верхней релаксации. Характеристикой вычислительной сложности алгоритмов служит количество арифметических операций (флопов), затраченных на вычислительный процесс.

Достоверность и обоснованность

Достоверность результатов подтверждена строгими математическими доказательствами, результаты согласуются с проведенными численными экспериментами.

Результаты, выносимые на защиту

  1. В диссертации предложен новый алгоритм, вычисляющий наибольшее собственное число симметричной матрицы, а также новый алгоритм, вычисляющий наибольшее сингулярное число несимметричной матрицы.
  2. Линейная сходимость предложенных алгоритмов теоретически доказана, эмпирически проверена.
  3. Установлена связь между предложенными алгоритмами и методом релаксации отношения Релея.
  4. Построена оптимальная реализация предложенных алгоритмов, обеспечивающая в вычислительном процессе минимальное количество арифметических операций на каждом шаге.
  5. Установлены оптимальные критерии, позволяющие останавливать вычислительный процесс при достижении заданной точности.
  6. Предлагаются некоторые модификации, позволяющие ускорить сходимость (в частности, применение техники верхней релаксации).
  7. Построена вычислительная схема, позволяющая применять матрицы отражения вместо матриц вращения.
  8. Произведено сравнение предложенных алгоритмов с ранее известными алгоритмами на основе их программной реализации и тестовых испытаний. Эксперименты проводились для матриц различной размерности, для различных относительных точностей расчета, для матриц с различной близостью старших собственных (сингулярных) чисел, для вычислительных процессов с разным коэффициентом верхней релаксации. Эмпирически установлены различные зависимости, выявлены наиболее эффективные алгоритмы для каждой конкретной задачи.

Научная новизна

В диссертации предложены новые алгоритмы решения задачи на собственные и сингулярные значения. Новыми также являются доказательства сходимости этих алгоритмов. Проведены сравнительные эксперименты, отражающие поведение предложенных и ранее известных алгоритмов на разных задачах.

Практическая ценность

Предложенные в диссертации алгоритмы могут использоваться для решения конкретных практических задач, приводящих к проблеме собственных значений, а также для дальнейшего изучения и совершенствования рассматриваемых методов. Полученные в третьей главе результаты численных экспериментов носят рекомендательный характер по выбору того или иного алгоритма в зависимости от поставленной задачи.

Апробация результатов

Основные положения и результаты, включенные в диссертацию, докладывались на конференциях и семинарах Санкт-Петербургского государственного университета, Московского инженерного-физического института (государственного университета), Санкт-Петербургского политехнического университета, Санкт-Петербургского электротехнического университета, на международных конференциях «Математика. Компьютер. Образование».

Публикации

По материалам диссертации опубликовано 13 научных работ, среди них 4 полноценные статьи, 9 – тезисы и краткие статьи. Работа [2] опубликована в журнале, рекомендованном ВАК.

Объем и структура работы

Диссертационная работа состоит из трех глав, заключения, списка публикаций по теме диссертации, списка литературы. Диссертация изложена на 109 листах машинописного текста, включает 30 рисунков, список научных публикаций содержит 13 наименований, список основной литературы – 21 наименование.

Содержание диссертации

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

Вторая глава посвящается теоретическому изучению новых алгоритмов.

В §1 представляется идея новых алгоритмов, заключающаяся в следующем. Пусть для симметричной матрицы требуется решить частичную проблему собственных значений, например вычислить наибольшее собственное число. Предположим, что строковые суммы элементов матрицы взаимно равны. Тогда вектор, состоящий из единиц, будет собственным вектором, а нахождение одного из собственных чисел сводится к вычислению суммы элементов какой-либо из строк. Например:

6 – собственное число.

Предположим, что строковые суммы элементов матрицы являются не строго равными, а приближенно равными. Например:

.

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

.

Для оценки близости приближенного значения собственного числа к точному значению собственного числа можно использовать теорему об оценке невязки:

где.

Для приведенного выше примера эта оценка дает: 0.017.

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

В §2 предлагается новый алгоритм вычисления наибольшего собственного числа произвольной симметричной матрицы, основанный на последовательных преобразованиях подобия. В качестве преобразующих матриц используются ортогональные матрицы двумерного вращения, изменяющие на каждом шаге алгоритма две строки и два столбца текущей матрицы. Для преобразований выбираются строки с минимальной и максимальной строковыми суммами элементов. Угол поворота матрицы вращения выбирается так, чтобы уравнивать строковые суммы преобразуемых двух строк. Показывается, что в пределе выполнения данного процесса все строковых сумм становятся равными, а значит, локализуют одно из собственных чисел.

Доказательство сходимости данного процесса основывается на доказательстве следующих утверждений.

Утверждение 1. На каждом шаге алгоритма существует хотя бы одно собственное значение, принадлежащее интервалу , где

Доказательство утверждения непосредственно следует из теоремы об оценке невязки, если в качестве приближенного собственного числа брать, а в качестве приближенного собственного вектора брать.

Утверждение 2. На каждом шаге алгоритма выполняемое преобразование подобия обеспечивает равенство сумм элементов строк i и j:, где – преобразованная матрица, и – номера преобразуемых строк и столбцов.

Доказательство утверждения непосредственно следует из вычислительной схемы алгоритма, в которой угол поворота в матрицах вращения выбирается соответствующим образом.

Утверждение 3. На каждом шаге алгоритма выполняемое преобразование подобия влечет увеличение суммы всех элементов матрицы.

Доказательство утверждения также непосредственно следует из вычислительной схемы алгоритма. Отмечается, что задача уравнивания сумм элементов двух строк (т.е. задача ) и задача максимизации суммы всех элементов матрицы (т.е. задача, где – угол поворота в матрице вращения) имеют одинаковые решения (реализуются одной матрицей вращения).

Утверждение 4. Сумма всех элементов матрицы в ходе выполнения алгоритма достигает некоторого предельного значения.

Доказательство утверждения основано на ограниченности сверху суммы всех элементов матрицы, преобразуемой матрицами двумерного поворота, причем верхним пределом является ( – размер матрицы, – наибольшее собственное число).

Лемма 1. Пусть векторы Пусть – острый угол между векторами : Тогда существует ортогональная.

Доказательство леммы основывается на непосредственном построении и изучении матрицы.

Утверждение 5. Справедливо неравенство

для.

Доказательство тривиально.

Лемма 2. Пусть – симметричная матрица, – собственные числа, – соответствующие собственные векторы,. В утверждении 1 было показано, что на интервале существует хотя бы одно собственное число. Пусть мало настолько, что на интервале существует ровно одно собственное число, и пусть оно простое. Тогда существуют ортогональная и симметричная такие, что

Доказательство основано на рассмотрении матрицы где Данная матрица характерна тем, что имеет собственное число и собственный вектор. Матрица рассматривается в отношении матрицы как «возмущенная» матрица. Применение теоремы о возмущении собственного числа и собственного вектора позволяет утверждать, что матрица имеет собственное число и собственный вектор, удовлетворяющие следующим неравенствам:

1) ;

2)

где – угол между векторами и.

Согласно лемме 1 существует ортогональная матрица, переводящая в :

.

Строится. Показывается, что

Утверждение 6. Справедливо неравенство

Доказательство основывается на рассмотрении суммы элементов матрицы как функции от :

Выполнив замены по формулам Тейлора первого порядка, формируется неравенство, где.

Взяв и учитывая, что, формируется неравенство.

Утверждение 7. Пусть – ортогональная матрица, Тогда и представляются в виде

где - кососимметричная,

Определяется, где логарифм задается рядом по степеням

Отмечается, что, из чего следует.

Оцениваются.

Утверждение 8. Начиная с некоторого шага алгоритма, имеет место неравенство

Согласно лемме 2 существуют ортогональная и симметричная :

,

По утверждению 7 возможно представление:

где – кососимметричная.

Тогда представляется в виде:

где.

Оценивается норма.

Показывается, что.

Тогда.

Теорема 1. Алгоритм имеет линейную скорость сходимости.

Доказательство основывается на рассмотрении утверждений 6 и 8:

,

.

Из них следует:

,

что доказывает линейную сходимость.

Pages:     || 2 |






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»