ВВЕДЕНИЕ 4
ГЛАВА 1. Обзор литературы 8
§1. Метод последовательных приближений 8
§2. Метод ускорения сходимости монотонных приближений к решению
уравнения вида х = Ах + f 15
§3. Метод однопараметрического итеративного агрегирования решения линейных операторных уравнений вида х = Ax + f, где оператор A - матрица n - го порядка 21
§4. Метод однопараметрического итеративного агрегирования решения нелинейных операторных уравнений вида х = F(х) + f, где F(х) - нелинейный оператор 27
ГЛАВА 2. Построение приближений, сходящихся к спектральному радиусу и собственному вектору линейного оператора 32
§5. Построение приближений, сходящихся к спектральному радиусу линейного оператора 32
§6. Построение приближений, сходящихся к собственному вектору линейного оператора 41
ГЛАВА 3. Развитие методов построения приближений, сходящихся к
точному решению операторного уравнения вида х = Ах + f 56
§7. Об одном итерационном методе решения системы линейных алгебраических уравнений вида х = Ах + f с квадратной матрицей А, в случае,
когда спектральный радиус матрицы А, больше чем единица 56
§8. Получение двусторонних оценок точного решения х* операторного уравнения вида х = Ах + f, в случае, когда спектральный радиус опера¬тора А не обязательно меньше единицы 65
§9. О некоторых подходах к уточнению границ решения операторных уравнений вида х = Ах + f в случае, когда спектральный радиус операто- 73 ра А не обязательно меньше единицы.
§10. “Гибрид” методов ускорения сходимости монотонных приближений к решению х* уравнения вида х = Ах + f и однопараметрического итеративного агрегирования 86
§11. Об одном варианте метода ускорения сходимости монотонных приближений к решению уравнения вида х = Ax + f 93
§12. Об одном варианте метода Зейделя 100
ЗАКЛЮЧЕНИЕ 112
ЛИТЕРАТУРА 114
ПРИЛОЖЕНИЕ 121
При решении широкого класса задач математического анализа, алгебры, экономики требуется находить решение операторных уравнений. В тех случаях, когда процесс отыскания точного решения является затруднительным, на помощь приходят итерационные методы, позволяющие найти приближенное решение с определенной степенью точности. Соответствующий класс задач можно представить с помощью операторного уравнения вида
х = Ax+f (1)
с линейным или нелинейным оператором А, действующим в банаховом пространстве Е, и свободным членом f из этого пространства.
Большое практическое значение приобретает возможность строить приближения un и, соответственно, vn к решению х* операторного уравнения вида (1), такие что
*
un < х ^ vn
При этом, оказывается, параллельно решаются две важные задачи теории приближенных методов решения операторных уравнений - задача об оценке погрешности приближенного решения, а также задача об априорной оценке относительной погрешности приближенного решения.
Использование современных ЭВМ открывает широкие возможности для решения таких задач. Математическое моделирование стало активно внедряться в практику научных и прикладных разработок при исследовании сложных явлений и процессов, происходящих в экономике. Лишь с помощью современных ЭВМ удается проводить численное моделирование достаточно сложных экономических процессов.
Существенные прикладные преимущества, проведенные исследования практической скорости сходимости итерационных процессов к точному решению х операторного уравнения вида (1), возможность контроля результатов в процессе решения и наличие математических обоснований, дают повод обратить серьезное внимание на итерационные методы не только как на объект интересных теоретических исследований, но и как на новые подходы к организации и проведению реальных расчетов плановых задач большой размерности и сложной структуры.
Актуальность проблемы. Балансовая модель производства является одной из наиболее простых математических моделей. Она записывается в виде системы уравнений, каждое из которых выражает требование равенства (баланса) между количеством продукции, производимой отдельным экономическим объектом, и совокупной потребности в этом продукте. Отсюда происходит название модели.
Впервые балансовые модели начали использоваться в СССР в 20-х годах. В более или менее законченном виде теория балансовых моделей была разработана американским ученым В.В. Леонтьевым в середине 30-х годов. Однако в те годы ни уровень развития математической науки, ни качество вычислительной техники не позволили широко распространить балансовый метод.
За разработку и внедрение в практику метода межотраслевого баланса группа советских экономистов под руководством академика А.Н. Ефимова в 1968 году была удостоена Г осударственной премии СССР. В настоящее время большое количество работ посвящается этой модели и ее применению для решения различных задач. Такой интерес к балансовой модели определяется тем, что, как оказалось, эта модель хорошо отображает многие существенные особенности современного производства и в то же время легко поддается расчету. Во многих странах мира балансовый метод используется для экономического анализа, планирования и прогнозирования.
В связи с внедрением ЭВМ в научные разработки, значительно повысился интерес к различным численным методам и алгоритмам, реализация которых граничит с проведением вычислительного эксперимента. Потребность в таком под¬ходе к решению задач математической экономики диктуется все усложняющимися запросами практики, а также связана с попыткой создания более рациональных и более общих теоретических моделей для изучения сложных экономических явлений. Активное использование методов численного моделирования позволяет резко сократить сроки научных и конструкторских разработок.
Проведенные в диссертационной работе исследования направленные на разработку новых методов решения операторных уравнений, описывающих экономические модели (модель межотраслевого баланса). Получены следующие научные и практические результаты.
1. Разработан и апробирован на большом количестве примеров итерационный метод решения системы линейных алгебраических уравнений вида х = Ах + f с квадратной матрицей A, в случае, когда наибольшее по модулю собственное значение матрицы А, больше чем единица.
2. Предложен метод получения двусторонних оценок точного решения х* операторного уравнения вида х = Ах + f, в случае, когда спектральный радиус не обязательно меньше единицы, а также подходы к уточнению полученных оценок. Метод проиллюстрирован соответствующими примерами.
3. Получен синтез методов ускорения сходимости монотонных приближений к решению х* уравнения вида х = Ах + f и однопараметрического итеративного агрегирования.
4. Предложен вариант метода ускорения сходимости монотонных приближений к решению уравнения вида х = Ах + f, в котором упрощена задача поиска начальных приближений.
5. Разработан и апробирован на большом количестве примеров вариант метода Зейделя, позволяющий строить двусторонние приближения к точному решению уравнения вида х = Ах + f.
6. Составлена библиотека программ на языке программирования TURBO PASCAL, которая позволяет реализовывать полученные в данной работе методы и алгоритмы.
Таким образом:
- Разработаны новые методы решения операторных уравнений, описывающих экономические модели (модель межотраслевого баланса), обладающих высокой скоростью сходимости последовательностей к точному решению данных уравнений, а также способностью сходится к точному решению даже в тех случаях, когда спектральный радиус оператора больше единицы.
- Разработан комплекс программ на языке программирования TURBO PASCAL, реализующих эти алгоритмы.
1. Архангельский, Ю.С. Численные исследования методов итеративного агрегирования для решения задачи межпродуктового баланса /Ю.С. Архангельский, И.А. Вахутинский, Л.М. Дудкин и др. //Автоматика и телемеханика. - 1975. - №7. - С.75-82.
2. Ашманов, С. А. Введение в математическую экономику /С. А. Атттманов. - М.:Наука, 1984. - 296с.
3. Бабаджанян, А. А. О скорости сходимости метода однопараметрического итеративного агрегирования/А. А. Бабаджанян //Автоматика и телемеханика. - 1982. - №11. - С.171-173.
4. Бахвалов Н. Численные методы /Н. Бахвалов, Н. Жидков, Г. Кобельков - М.: Лаборатория Базовых Знаний, 2000. - 624 с.
5. Бахтин, И. А. Исследование уравнений с положительными операторами: Дис. ... д-ра физ.-мат. наук /И.А. Бахтин. - Ленинград, 1967. - 320с.
6. Бахтин, И.А. Метод последовательных приближений в теории уравнений с вогнутыми операторами /И.А. Бахтин, М.А. Красносельский // Сибирский математический журнал . - 1961.- Т.2, № 3.- С.313-330.
7. Бахтин, И.А. О непрерывности положительных операторов /И.А. Бахтин, М.А. Красносельский, В.Я. Стеценко // Сибирский математический журнал . - 1962.- Т.3, № 1.- С.8-17.
8. Беллман Р. Квазилинеаризация и нелинейные краевые задачи /Р. Беллман, Р. Калаба. - М.: Мир, 1968. - 270с.
9. Вен, В.Л. Некоторые вопросы агрегирования линейных моделей /В.Л. Вен, А.И. Эрлих //Известия АН СССР. Сер. техническая ки¬бернетика. - 1970.- №5. - С.3-8.
10. Вержбицкий, В.М. Численные методы (линейная алгебра и нели¬нейные уравнения): Учеб. Пособие для вузов /В.М. Вержбицкий. - М.: Высш. шк., 2000. - 266 с.
11. Вулих, Б.З. Введение в теорию полуупорядоченных пространств /Б.З. Вулих. - М.: Наука, 1961. - 407 с.
12. Вулих, Б.З. Введение в функциональный анализ /Б.З. Вулих. - М.: Физматгиз, 1967.- 415с.
13. Вулих, Б.З. Специальные вопросы геометрии конусов в нормиро¬ванных пространствах: Учебное пособие /Б.З. Вулих. - Калинин: Издательство калининского университета, 1978. - 84 с.
14. Гантмахер, Ф.Р. Теория матриц /Ф.Р. Гантмахер. - М: Наука, 1966. - 576с.
15. Гробова, Т. А. О методе однопараметрического итеративного агре¬гирования для решения систем линейных и нелинейных алгебраиче¬ских уравнений, интегральных уравнений /Т. А. Гробова
115
//Ставропольский государственный университет, Ставрополь, 2001.
- 24с. - Деп. в ВИНИТИ 19.11.01 №2392 - В2001.
16. Гробова, Т.А. Об одном новом варианте метода Зейделя /Т.А. Гро- бова //Научно-методическая конф. преподавателей и студентов «XXI век - век образования». Материалы 46 научно-метод. конф. преподавателей и студентов «XXI век - век образования». - Ставро¬поль, 2001. - С.4-9.
17. Гробова, Т.А., Об одном аналоге метода однопараметрического итеративного агрегирования /Т.А. Гробова, В.Я. Стеценко //Вестник СГУ. -2001. - Выпуск 28. - С.12-16.
18. Гробова, Т.А. Об одной новой схеме реализации вариантов метода Зейделя /Т.А. Гробова, В.Я. Стеценко // Вестник молодых ученых. - Санкт - Петербург, 2001. - С.34-39.
19. Данфорд Н., Шварц Дж. Линейные операторы. Общая теория /Н. Данфорд, Дж. Шварц. - М.: Иностранная литература, 1962. - 895c.
20. Демиденко, Н.А. Применение метода итеративного агрегирования к расширенной модели межотраслевого баланса /Н.А. Демиденко //Экономика и математические методы. - 1977. - Т.13, №3. - С.594- 598.
21. Дудкин, Л.М. Межотраслевой баланс и материальные балансы от¬дельных продуктов /Л.М. Дудкин, Э.Б. Ершов // Плановое хозяйст¬во. - 1965. - №5. - С.59-63.
22. Есаян, А.Р. Локализация спектра линейного оператора /А.Р. Есаян,
В.Я. Стеценко //Междунар. Конгресс математиков (1966; Москва). Тезисы кр. науч. сообщений Междунар. Конгресса математиков. Секция 5. - М., 1966. - С.45-74.
23. Есаян, А.Р. О разрешимости уравнений второго рода /А.Р. Есаян,
В.Я. Стеценко // Труды семинара по функциональному анализу. Во¬ронежский гос. ун-т. - Воронеж, 1963. - Вып. 7. - С.36-41.
24. Есаян, А.Р. Оценки спектра интегральных операторов и бесконеч¬ных матриц /А.Р. Есаян, В.Я. Стеценко // Докл. АН СССР. - 1964. - Т. 157, №2. - С.12-19.
25. Итеративное агрегирование и его применение в планировании /Под ред. Л.М. Дудкина. - М.: Экономика, 1979. - 328 с.
26. Канторович, Л.В. Функциональные анализ в нормированных про¬странствах /Л.В. Канторович, Г.П. Акилов. - М.: Наука, 1977. - 496 с.
27. Канторович, Л.В. Функциональный анализ в полуупорядоченных пространствах /Л.В. Канторович, Б.З. Вулих, А.Г. Пинскер. - М.: Физматгиз, 1959. - 684 с.
28. Канторович, Л.В. Приближенные методы высшего анализа /Л.В. Канторович, В.И. Крылов. - М.-Л.: Физматгиз, 1962. - 708с.
116
29. Коллатц Л. Функциональный анализ и вычислительная математика /Л. Коллатц. - М.: Мир, 1969. - 421с.
30. Колмогоров, А.Н. Элементы теории функций и функционального анализа /А.Н. Колмогоров, С.В. Фомин. - М.: Наука, 1981. - 543с.
31. Коршунова Н. Математика в экономике /Н. Коршунова, В. Плясу¬нов. - М.: Издательство «Вита-Пресс», 1996. - 368с.
32. Костенко, Т.А. К вопросу о существовании и единственности реше¬ния операторного уравнения, нелинейного относительно параметра Л. /Т.А. Костенко // «Международная школа-семинар по геометрии и анализу памяти Н.В. Ефимова» (1998; Ростов-на-Дону). Тезисы док¬ладов «Междунар. школы-семинара по геометрии и анализу памяти
Н.В. Ефимова» (5-11 сент., 1998г.). - Ростов-на-Дону, 1998. - С.104- 105.
33. Костенко, Т.А. О разрешимости операторных уравнений второго рода с линейными и нелинейными операторами /Т.А. Костенко // «Проблемы физико-математических наук», науч.-метод. конф. пре¬подавателей и студентов «Университетская наука - региону» (43; 1998; Ставрополь). Материалы XLIII науч.-метод. конф. преподава¬телей и студентов «Университетская наука - региону». - Ставро¬поль, 1998. - С.111-122.
34. Красносельский, М.А. Положительные решения операторных урав¬нений /М. А. Красносельский. - М.: Физматгиз, 1962. - 394с.
35. Красносельский, М.А. Правильные и вполне правильные конусы /М.А. Красносельский // Докл. АН СССР.-1960. - Т.135. - № 2. - С.241-255.
36. Красносельский, М.А. Приближенное решение операторных урав¬нений /М.А. Красносельский, Г.М. Вайникко, П.П. Забрейко и др. - М: Наука, 1969. - 456с.
37. Красносельский, М.А. Геометрические методы нелинейного анализа /М.А. Красносельский, П.П. Забрейко. - М.: Наука, 1965. - 624с.
38. Красносельский, М.А. Позитивные линейные системы: метод поло¬жительных операторов /М.А. Красносельский, Е.А. Лифшиц, В.И. Соболев. - М: Наука, 1985. - 256с.
39. Красносельский, М.А. Положительно обратимые линейные опера¬торы и разрешимость линейных уравнений /М.А. Красносельский, Е.А. Лифшиц, В.В. Покорный, В.Я. Стеценко // Докл. АН Таджик¬ской ССР. -1974. - Т.ХУ11, № 1. - С.12-15.
40. Красносельский, М.А. О сходимости метода однопараметрического агрегирования /М.А. Красносельский, А.Ю. Островский, А.В. Собо¬лев // Автоматика и телемеханика. - 1978. - №9. - С.102-109.
117
41. Красносельский, М.А. Замечания о методе Зейделя /М.А. Красно¬сельский, В. Я. Стеценко // Журнал вычислительной математики и математической физики. - 1969. - Т.9, №1. - С.177-182.
42. Крейн, М.Г. Линейные операторы, оставляющие инвариантным ко¬нус в пространстве Банаха /М.Г. Крейн, М.А. Рутман // Успехи ма¬тематических наук. - 1948. - Т.1, №3. - С.3-95.
43. Крукиер, Л.А. Численные методы решения задач конвекции- диффузии со смешанными производными /Л.А. Крукиер, Т.С. Мар¬тынова. - г. Ростов-на -Дону: Изд-во РГУ, 2003. -156с.
44. Кубекова, Б.С. Об уточнении оценок решения операторного уравне¬ния в полуупорядоченном пространстве с и0-ограниченным снизу оператором /Б.С. Кубекова // «Понтрягинские чтения -XI», науч. конф. (2000; Воронеж). Тезисы докладов науч. конф. «Понтрягин¬ские чтения -XI» (3-9 мая, 2000г.). - Воронеж, 2000. - С.95-98.
45. Кубекова, Б.С. Отыскание приближений по недостатку и по избытку к решению операторного уравнения с монотонно разложимым опе¬ратором /Б.С. Кубекова //«Математическое моделирование в науч¬ных исследованиях», Всероссийская науч. конф. (2000). Тезисы докладов Всероссийской науч. конф. «Математическое моделирова¬ние в научных исследованиях» (21-30 сент., 2000г.). - Ставрополь,
2000. - С.47-49.
46. Кубекова, Б.С. О методе однопараметрического итеративного агре¬гирования /Б.С. Кубекова, В.Я. Стеценко, Т. А. Гробова //«Математика. Компьютер. Образование», междунар. конф. (8; 2001; Пущино). Тезисы докладов восьмой междунар. конф. «Мате¬матика. Компьютер. Образование» (31 янв. - 5 февр., 2001г.).- Пу¬щино, 2001. - С.230-232.
47. Кубекова, Б.С. Об одном методе построения двусторонних прибли¬жений к решению операторного уравнения с монотонно разложи¬мым оператором /Б.С. Кубекова, В.Я. Стеценко, М.Н. Павлова //Журнал вычислительной математики и математической физики.-
2001. - Т.41, № 6. - С.846-854.
48. Кузнецов, Ю.А. К теории итерационных процессов /Ю.А. Кузнецов //Докл. АН СССР. - 1969. - Т.184, №4, -С.863-866.
49. Леонтьев, В.В. Экономика и математические методы /В.В. Леонтьев, Д. Форд. - М: Наука, 1972. - 242с.