1626435584-7c6402f545ecf856225d6cf8d21519c9 (844233), страница 2
Текст из файла (страница 2)
В настоящее время ЭВМ стали обычным оборудованием многих институтов и конструкторских бюро. Это позволило от простейших расчетов и оценок различных конструкций или процессов перейти к новой стадии работы — детальному математическому моделированию (вычислительному эксперименту), которое существенно сокращает потребность в натурных экспериментах, а в ряде случаев может их заменить. В основе вычислительного эксперимента лежит решение уравнений математической модели численными методами. Изложению численных методов посвящено немало книг.
Однако большинство этих книг ориентировано на студентов и научных работников математического профиля. Поэтому в настоящее время ощущается потребность в книге, рассчитанной на широкий круг читателей различных специальностей и сочетающей достаточную полноту изложения с разумной степенью строгости при умеренном объеме. Предлагаемая книга отвечает этим требованиям. Она достаточно полно освещает тот круг вопросов, знание которого наиболее часто требуется в практике вычислений, и содержит ряд разделов, которые редко включают в учебные пособия.
Умеренный объем достигнут за счет тщательного отбора материала и включения в книгу только наиболее эффективных и часто используемых на практике методов. Материал изложен четко и сжато, при этом большое внимание уделено рекомендациям по практическому применению алгоритмов; изложение пояснено рядом примеров. Для обоснования алгоритмов использован несложный математический аппарат, знакомый студентам физических и инженерных специальностей. Книга рассчитана на читателя, который занимается не столько разработкой численных методов, сколько их применением к прикладным проблемам. Однако в процессе работы над книгой читатель знакомится с основными идеями построения вычислительных алгоритмов и с их обоснованием и приобретает знания, достаточные для разработки новых алгоритмов.
Эта книга является по существу введением в численные методы. Овладев ею, читатель затем может углубить свои знания, обратившись к руководствам по теории разностных схем и по методам численного решения отдельных классов задач. Книга написана специалистом по теоретической и математической физике. Она возникла в результате работы автора над рядом актуальных проблем физики в Институте прикладной математики АН СССР и преподавания на физическом факультете МГУ. Несомненно, книга окажется полезной широкому кругу читателей — студентам, аспирантам, научным сотрудникам и инженерам математических, физических и технических специальностей. А. А.
Самарский ПРЕДИСЛОВИЕ Сложные вычислительные задачи, возникающие при исследовании физических и технических проблем, можно разбить на ряд элементарных — таких как вычисление интеграла, решение дифференциального уравнения и т. п. Многие элементарные задачи являются несложными и хорошо изучены. Для этих задач уже разработаны методы численного решения, и нередко имеются стандартные программы решения их на ЭВМ. Есть и достаточно сложные элементарные задачи; методы решения таких задач сейчас интенсивно разрабатываются (например, решение уравнений бесстолкновительной плазмы). Поэтому полная программа обучения численным методам должна состоять из ряда этапов. Вопервых, это освоение логарифмической линейки, клавишных вычислительных машин и программирования на ЭВМ.
Во-вторых, основы численных методов, содержащие изложение классических элементарных задач (включая основные сведения о разностных схемах). В-третьих, курс теории разностных схем. И в-четвертых — ряд специальных курсов, которые сейчас нередко называют методами вычислительной физики: численное решение задач газодинамики, аэродинамики, переноса излучения, квантовой физики, квантовой химии и т. д. Эта книга является введением в численные методы. Она начинается с простейших задач интерполирования функций и кончается недавно возникшим разделом вычислительной математики— методами решения некорректно поставленных задач. Книга написана на основе годового курса лекций, читавшихся автором сначала инженерам-конструкторам, а после переработки †студент физического факультета МГУ.
Для каждой задачи существует множество методов решения. Например, хорошо обусловленную систему линейных уравнений можно решать методами Гаусса, Жордана, оптимального исключения, окаймления, отражений, ортогонализации и рядом других. Интерполяционный многочлен записывают в формах Лагранжа, Ньютона, Грегори — Ньютона, Бесселя, Стирлинга, Гаусса и Лапласа— Эверетта. Подобные методы обычно являются вариациями одного-двух основньгх методов, и если даже в каких-то частных случаях имеют преимущества, то незначительные.
Кроме того, многие методы создавались до появления ЭВМ, и ряд из них в качестве существенного элемента включает интуицию вычислителя. Появление ЭВМ потребовало переоценки старых методов, что до конца еще не сделано, и до сих пор по традиции большое количество неэффективных методов кочует из учебника в учебник. Отчасти это объясняется тем, что эффективность многих методов сильно зависит от мелких деталей алгоритма, почти не поддающихся теоретическому анализу; поэтому окончательный отбор лучших методов можно сделать только на основании большого опыта практических расчетов.
В этой книге сделана попытка такого отбора, опирающаяся на многолетний опыт решения большого числа разнообразных задач математической физики. Для большинства рассмотренных в книге задач изложены только наиболее эффективные методы с широкой областью применимости. Несколько методов для одной и той же задачи даны в том случае, если они имеют существенно разные области применимости, или если для данной задачи еще не разработано достаточно удовлетворительных методов. Часто приходится слышать, что наступила эпоха ЭВМ, а «ручные» расчеты являются архаизмом.
На самом деле это далеко не так. Прежде чем поручать ЭВМ большую задачу, надо сделать много оценочных расчетов и на их основе понять, какие методы окажутся эффективными для данной задачи. Конечно, даже в мелких расчетах ЭВМ с хорошим математическим обеспечением и набором периферийных устройств (телетайп, дисплей, графико-построитель) оказывает большую пользу. Однако логарифмическая линейка и клавишные машины еще долго будут необходимы. Поэтому большинство методов, изложенных здесь, в равной мере пригодны для ЭВМ и «ручных» расчетов.
Основное внимание в книге уделено выработке практических навыков у читателя. Поэтому в первую очередь изложены алгоритмы, даны рекомендации по их применению и отмечены «маленькие хитрости» вЂ” те незначительные на первый взгляд практические приемы, которые сильно повышают эффективность алгоритма. Теоретическое обоснование методов приведено лишь в той мере, в какой оно необходимо для лучшего усвоения и практического применения. В книгу включен ряд сведений, не относящихся к необходимому минимуму, но полезных читателю для лучшего понимания тонких деталей вычислительных процессов. Чтобы не увеличивать объем книги и избежать сложных выкладок, эти сведения приведены, как правило, без доказательств, но со ссылками на дополнительную литературу.
Некоторые сведения даны в форме задач в конце каждой главы. Предполагается, что читатели знакомы с основами высшей математики, включая краткие сведения об уравнениях в частных производных. Необходимые дополнительные сведения, которые не содержатся в обязательных курсах университетов и втузов, сообщаются здесь в соответствующих разделах.
Книга разделена на главы; параграфы и пункты. В начале каждой главы кратко изложено ее содержание. Нумерация таблиц и рисунков — единая по всей книге, а нумерация формул— самостоятельная в каждой главе. Если ссылка не выходит за пределы данной главы, то указывается только номер формулы; если выходит — то номер главы и номер формулы. В конце книги дан список литературы.
Приведенные в нем учебники и монографии рекомендуются для углубленного изучения отдельных разделов. Журнальные статьи даны для указания на оригинальные работы, их список не претендует на полноту; более полная библиография имеется в рекомендованных учебниках. Общий подход к теории и практике вычислений, определивший стиль этой книги, сложился у меня под влиянием А. А. Самарского и В.
Я. Гольдина за много лет совместной работы. Ряд актуальных тем был включен по инициативе, А. Г. Свешникова и В. Б. Гласко. Много ценных замечаний сделали А. В. Гулин, Б. Л. Рождественский, И. М. Соболь, И. В. Фрязинов, Е. В. Шикин и сотрудники кафедры прикладной математической физики МИФИ. В оформлении рукописи мне помогли Л. В. Кузьмина и В. А. Кра-сноярова. Я пользуюсь случаем искренне поблагодарить всех названных лиц, и в особенности Александра Андреевича Самарского.
Н. Н. Калиткик ГЛАВА! ЧТО ТАКОЕ ЧИСЛЕННЫЕ МЕТОДЫ? Глзва 1 является вводной. В 1 1 рассмотрены роль математики при решении физико-технических задач и место чисзенных методов среди других математических методов и кратко наложена история численных методов. В 1 2 разобраны основные понятия приближенного анализа: корректность постановки задач, определение близости точного и приближенного решений, структура погрешности. В 1. Математические модели и численные методы 1. Решение задачи, Физиков математика интересует не сама по себе, а как средство решения физических задач.
Рассмотрим поэтому, как решается любая реальная задача — например, нахождение светового потока конструируемой лампы, производительности проектируемой химической установки или себестоимости продукции строящегося завода. Одним из способов решения является эксперимент. Построим эту лампу, установку или завод и измерим интересующую нас характеристку.
Если характеристика оказалась неудачной, то изменим проект и построим новый завод и т. д. Ясно, что мы получим достоверный ответ на вопрос, но слишком медленным и дорогим способом. Другой способ в математический анализ конструкции или явления. Но такой анализ применяется не к реальным явлениям, а к некоторым математическим моделям этих явлений.
Поэтому первая стадия работы †э формулировка математической модели (постановка задачи). Для физического процесса модель обычно состоит из уравнений, описывающих процесс; в эти уравнения в виде коэффициентов входят характеристики тел или веществ, участвующих в процессе. Например, скорость ракеты при вертикальном полете в вакууме определяЕтся уравнением где М вЂ” начальная масса ракеты, т(г) — заданный расход горю- 14 [ГЛ. ! ЧТО ТАКОЕ ЧИСЛЕННЫЕ МЕТОДЫ? чего, д — ускорение поля тяготения, а с — скорость истечения газов, зависящая от калорийности топлива и среднего молекулярного веса продуктов сгорания, Любое изучаемое явление бесконечно сложно.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
