Введение в вычислительную физику - Федоренко Р.П.
ISBN 5-7417-0002-0
Скачать (прямая ссылка):
Вторая часть книги точнее соответствует содержанию, термина «вычислительная физика». В ней собраны описания методов приближенного решения частных задач, имеющих, однако, важные области приложения в современной науке и технике. Каждый параграф посвящен одной из таких задач. Принят следующий способ изложения. Вначале дается замкнутая математическая формулировка задачи, указывается ее «прикладное происхождение». Физическая терминология используется для «оживления» изложения, но никакого физического обоснования постановки задачи не проводится — это дело физика, а не вычислителя. При этом указываются те особые • обстоятельства, которые делают задачу нестандартной, требующей разработки специальных вычислительных методов.
Затем описывается метод приближенного решения, оказавшийся достаточно эффективным. Основное внимание уделяется именно тем деталям метода, которые учитывают специфическую нестандартность данной задачи и которым метод обязан своей эффективностью. Попутно обсуждаются те трудности, с которыми сталкиваются при стандартном подходе к задаче (формально не только возможном, но иногда даже строго обоснованном). В некоторых случаях приводятся и обсуждаются характерные численные результаты. Стандартные детали вычислительной методики описываются бегло, а иногда и совсем опускаются.
Материал второй части книги несет двойную нагрузку. Во-первых, Описываемые задачи достаточно интересны в приложениях и опыт их успешного решения представляет прямой интерес в связи с задачами именно этого типа. Во-вторых, разработка эффективного
8
ПРЕДИСЛОВИЕ
алгоритма частной задачи обычно связана с использованием приемов, имеющих более широкое, выходящее за рамки данной задачи значение. Автор предпочитает знакомить читателя с такими приемами на примерах конкретных задач, в которых они были использованы с большим эффектом. Есть и другой путь — выделить эти приемы как отдельные самостоятельные сюжеты, дать абстрактное описание ситуаций, в которых их применение целесообразно. Подобный способ изложения представляется нам чуждым духу вычислительной физики.
Отметим некоторые технические детали изложения. Текст книги разбит на параграфы, каждьій из которых имеет свою нумерацию формул. При ссылке на формулы другого параграфа используется двойной номер (параграфа и формулы). Впрочем, автор стремился свести к минимуму подобные ссылки. В тексте опускаются и библиографические ссылки. Этот недостаток компенсируется библиографическим комментарием, тем более необходимым, что во многих местах излагаются результаты, еще не вошедшие прочно в учебную литературу и часто освещенные лишь в журнальных, а то и ротап-ринтных публикациях. Курсивом в тексте выделены общеупотребительные термины вычислительной математики.
Оба курса, на основе которых написана эта книга, читались по предложению академика О. М. Белоцерковского, много сделавшего для внедрения компьютерных наук в «систему физтеха». Пользуюсь случаем высказать Олегу Михайловичу свою искреннюю благодарность.
Автор должен отметить и неоценимое влияние, оказанное на него коллегами по Институту прикладной математики им. М. В. Келдыша. Нет возможности упоминать их здесь, автор постарался должным образом отметить их вклад в развитие предмета книги в библиографическом комментарии к списку литературы. Это будет полезно и для будущих историков науки, которым рано или поздно предстоит изучать историю становления отечественной вычислительной математики.
ЧАСТЬ ПЕРВАЯ
ОСНОВЫ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ
§ 1. Решение систем нелинейных уравнений
*
В самых различных задачах вычислительной физики часто возникает необходимость решать системы нелинейных уравнений. Такую систему будем записывать кратко в виде
/(х)=0, (1)
имея в виду, что х есть точка п-мерного пространства, т.е. х =
= {xj, х2, хп}, а / — «-мерная вектор-функция, т.е. / = = {/i« /г> •••» /«}• Таким образом, (I) есть система п уравнений с п неизвестными:
/,(х,, Xv ..., хп) = 0,
(2)
fn(xр х2, ..., х„) =0.
Конечно, когда речь идет о нелинейных уравнениях в общем случае, мы не имеем ни теорем о существовании, ни теорем о единственности решения. Тем не менее, имея дело с системой (1), мы
предполагаем, что искомое решение существует. Оно, быть может, не единственно, и методы, которые будут рассмотрены ниже, не имеют целью найти все решения; обычно достаточно будет какого-то одного. Более того, мы предполагаем, что из каких-то содержательных соображений нам известно примерное расположение этого решения, некоторая не очень большая область, в которой оно находится. Таким образом, лучше говорить не о решении системы нелинейных уравнений, а об уточнении имеющегося весьма грубого приближения к некоторому решению. Ниже это будет должным образом конкретизировано.
Метод Ньютона. Основная вычислительная конструкция, применяемая для решения системы (1), по традиции приписывается Ньютону, хотя теоретические исследования этого алгоритма были выполнены лет на сто позже (Фурье, Коши). Основу метода состав-
