Уравнение в частных производных - Partial differential equation

Визуализация решения двумерной уравнение теплопроводности с температурой, представленной в вертикальном направлении

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

Функцию часто считают "неизвестной", которую нужно решить, аналогично тому, как Икс рассматривается как неизвестное число, которое необходимо решить в алгебраическом уравнении, таком как Икс2 − 3Икс + 2 = 0. Однако обычно невозможно выписать явные формулы для решений уравнений в частных производных. Соответственно, существует огромное количество современных математических и научных исследований методов численно приблизительный решения некоторых дифференциальных уравнений в частных производных с помощью компьютеров. Уравнения с частными производными также занимают большой сектор чисто математическое исследование, в котором обычные вопросы, вообще говоря, касаются выявления общих качественных особенностей решений различных уравнений в частных производных.[нужна цитата ]

Уравнения с частными производными повсеместно используются в математически ориентированных научных областях, таких как физика и инженерное дело. Например, они лежат в основе современного научного понимания звука, тепла, распространение, электростатика, электродинамика, динамика жидкостей, эластичность, общая теория относительности, и квантовая механика.[нужна цитата ] Они также возникают из многих чисто математических соображений, таких как дифференциальная геометрия и вариационное исчисление; среди других известных приложений, они являются фундаментальным инструментом в доказательстве Гипотеза Пуанкаре из геометрическая топология.

Частично из-за такого разнообразия источников существует широкий спектр различных типов дифференциальных уравнений в частных производных, и были разработаны методы для работы со многими возникающими отдельными уравнениями. Таким образом, обычно признается, что не существует «общей теории» уравнений в частных производных, а специальные знания в некоторой степени разделены между несколькими существенно различными подполями.[1]

Обыкновенные дифференциальные уравнения образуют подкласс дифференциальных уравнений в частных производных, соответствующих функциям одной переменной. Стохастические уравнения в частных производных и нелокальные уравнения являются, по состоянию на 2020 год, особенно широко изученными расширениями понятия «PDE». Более классические темы, по которым все еще ведется много активных исследований, включают: эллиптический и параболический уравнения в частных производных, механика жидкости, Уравнения Больцмана, и диспергирующий уравнения в частных производных.

Вступление

Говорят, что функция ты(Икс, у, z) трех переменных "гармонический "или" решение то Уравнение лапласа "если он удовлетворяет условию

Такие функции широко изучались в девятнадцатом веке из-за их актуальности для классическая механика. Если функция задана явно, то проверка гармоничности ее обычно сводится к простому вычислению. Например

оба гармоничны, в то время как

не является. Может показаться удивительным, что два приведенных примера гармонических функций имеют столь разительно отличающуюся друг от друга форму. Это отражение того факта, что они нет, в любом случае, оба частных случая «общей формулы решения» уравнения Лапласа. Это разительно отличается от случая обыкновенные дифференциальные уравнения (ОДУ) примерно похоже к уравнению Лапласа, с целью многих вводных учебников найти алгоритмы, приводящие к общим формулам решения. Для уравнения Лапласа, как и для большого числа дифференциальных уравнений в частных производных, такие формулы решения не существуют.

Более конкретно природу этого отказа можно увидеть в случае следующей PDE: для функции v(Икс, у) двух переменных рассмотрим уравнение

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

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

Чтобы обсудить такие теоремы существования и единственности, необходимо уточнить домен «неизвестной функции». В противном случае, говоря только в терминах типа «функция двух переменных», невозможно осмысленно сформулировать результаты. То есть область определения неизвестной функции должна рассматриваться как часть структуры самого PDE.

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

  • Позволять B обозначим круг единичного радиуса вокруг начала координат на плоскости. Для любой непрерывной функции U на единичном круге есть ровно одна функция ты на B такой, что
и чье ограничение на единичную окружность дается формулой U.
  • Для любых функций ж и грамм на реальной линии , есть ровно одна функция ты на ℝ × (−1, 1) такой, что
и с ты(Икс, 0) = ж(Икс) и ты/у(Икс, 0) = грамм(Икс) для всех значений Икс.

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

  • Если ты это функция на 2 с
тогда есть числа а, б, и c с ты(Икс, у) = топор + к + c.

В отличие от предыдущих примеров, это PDE нелинейный, благодаря квадратным корням и квадратам. А линейный PDE является таким, что сумма любых двух решений также является решением, и все постоянные кратные любого решения также являются решением.

Корректность

Корректность относится к общему схематическому пакету информации о PDE. Чтобы сказать, что PDE корректна, нужно иметь:

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

Это несколько расплывчато из-за необходимости применимости к нескольким различным PDE. В частности, требование «непрерывности» неоднозначно, поскольку обычно существует множество неэквивалентных средств, с помощью которых оно может быть строго определено. Однако довольно необычно изучать PDE без определения того, каким образом она правильно сформулирована.

Наличие локальных решений

В слегка слабой форме Теорема Коши – Ковалевского по существу заявляет, что если все члены в дифференциальном уравнении в частных производных состоят из аналитические функции, то в определенных областях обязательно существуют решения УЧП, которые также являются аналитическими функциями. Хотя это фундаментальный результат, во многих ситуациях он бесполезен, поскольку невозможно легко контролировать область получаемых решений. Кроме того, известны примеры линейных дифференциальных уравнений в частных производных, коэффициенты которых имеют производные всех порядков (которые, тем не менее, не являются аналитическими), но вообще не имеют решений: это удивительно пример был обнаружен Ганс Леви в 1957 г. Таким образом, теорема Коши-Ковалевского обязательно ограничивается в своем объеме аналитическими функциями. Этот контекст исключает многие явления, представляющие как физический, так и математический интерес.

Классификация

Обозначение

При написании PDE частные производные обычно обозначают индексами. Например:

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

Греческая буква Δ обозначает Оператор Лапласа; если ты является функцией п переменные, то

В физической литературе оператор Лапласа часто обозначается как 2; в математической литературе, 2ты может также обозначать гессианская матрица из ты.

Уравнения первого порядка

Линейные и нелинейные уравнения

PDE называется линейный если оно линейно относительно неизвестного и его производных. Например, для функции ты из Икс и у, линейное уравнение в частных производных второго порядка имеет вид

куда ая и ж являются функциями только независимых переменных. (Часто смешанные частные производные тыху и тыyx будет приравнено, но это не требуется для обсуждения линейности.) Если ая константы (не зависящие от Икс и у), то PDE называется линейный с постоянными коэффициентами. Если ж равен нулю всюду, то линейное уравнение в частных производных однородный, иначе это неоднородный. (Это отдельно от Асимптотическая гомогенизация, в которой изучается влияние высокочастотных колебаний коэффициентов на решения УЧП.)

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

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

Многие из фундаментальных УЧП в физике квазилинейны, например Уравнения Эйнштейна из общая теория относительности и Уравнения Навье-Стокса описание движения жидкости.

УЧП без каких-либо свойств линейности называется полностью нелинейный, и обладает нелинейностями на одной или нескольких производных высшего порядка. Примером может служить Уравнение Монжа – Ампера, возникающий в дифференциальная геометрия.[2]

Линейные уравнения второго порядка

Эллиптический, параболический, и гиперболический Уравнения в частных производных второго порядка широко изучаются с начала двадцатого века. Однако есть много других важных типов PDE, включая Уравнение Кортевега – де Фриза. Есть также гибриды, такие как Уравнение Эйлера – Трикоми, которые варьируются от эллиптических до гиперболических для разных областей области. Существуют также важные расширения этих базовых типов для PDE более высокого порядка, но такие знания являются более специализированными.

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

где коэффициенты А, B, C... может зависеть от Икс и у. Если А2 + B2 + C2 > 0 над регионом ху-плоскость, PDE второго порядка в этой области. Эта форма аналогична уравнению для конического сечения:

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

Так же, как классифицируют конические секции и квадратичные формы на параболические, гиперболические и эллиптические на основе дискриминант B2 − 4AC, то же самое можно сделать для PDE второго порядка в данной точке. Тем не менее дискриминант в PDE определяется как B2AC из-за соглашения ху срок 2B скорее, чем B; формально дискриминант (ассоциированной квадратичной формы) есть (2B)2 − 4AC = 4(B2AC), с пониженным коэффициентом 4 для простоты.

  1. B2AC < 0 (эллиптическое уравнение в частных производных ): Решения эллиптические уравнения в частных производных настолько гладкие, насколько позволяют коэффициенты, в пределах области, где определены уравнение и решения. Например, решения Уравнение Лапласа являются аналитическими в той области, в которой они определены, но решения могут принимать негладкие граничные значения. Движение жидкости с дозвуковыми скоростями можно аппроксимировать с помощью эллиптических УЧП, а уравнение Эйлера – Трикоми является эллиптическим, где Икс < 0.
  2. B2AC = 0 (параболическое уравнение в частных производных ): Уравнения, которые параболический в любой точке может быть преобразован в форму, аналогичную уравнение теплопроводности заменой независимых переменных. Решения сглаживаются по мере увеличения преобразованной временной переменной. Уравнение Эйлера – Трикоми имеет параболический тип на прямой, где Икс = 0.
  3. B2AC > 0 (гиперболическое уравнение в частных производных ): гиперболический уравнения сохраняют любые разрывы функций или производных в исходных данных. Примером может служить волновое уравнение. Движение жидкости со сверхзвуковой скоростью можно аппроксимировать с помощью гиперболических УЧП, а уравнение Эйлера – Трикоми является гиперболическим, где Икс > 0.

Если есть п независимые переменные Икс1, Икс2 ,… Иксп, общее линейное уравнение в частных производных второго порядка имеет вид

Классификация зависит от сигнатуры собственных значений матрицы коэффициентов ая,j.

  1. Эллиптический: все собственные значения либо положительные, либо отрицательные.
  2. Параболический: все собственные значения либо положительные, либо отрицательные, за исключением нулевого.
  3. Гиперболический: есть только одно отрицательное собственное значение, а все остальные положительные, или есть только одно положительное собственное значение, а все остальные отрицательны.
  4. Ультрагиперболический: существует более одного положительного собственного значения и более одного отрицательного собственного значения, и нет нулевых собственных значений. Существует лишь ограниченная теория ультрагиперболических уравнений (Курант и Гильберт, 1962).

Системы уравнений первого порядка и характеристические поверхности

Классификацию дифференциальных уравнений в частных производных можно распространить на системы уравнений первого порядка, в которых неизвестные ты сейчас вектор с м компоненты, а матрицы коэффициентов Аν находятся м к м матрицы для ν = 1, 2,… п. Уравнение в частных производных принимает вид

где матрицы коэффициентов Аν и вектор B может зависеть от Икс и ты. Если гиперповерхность S дается в неявной форме

куда φ имеет ненулевой градиент, то S это характерная поверхность для оператора L в данной точке, если характеристическая форма обращается в нуль:

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

  1. Система первого порядка Лу = 0 является эллиптический если для L: значения ты на S а дифференциальное уравнение всегда определяет нормальную производную от ты на S.
  2. Система первого порядка - это гиперболический в момент, если есть космический поверхность S с нормальным ξ в таком случае. Это означает, что для любого нетривиального вектора η ортогонален ξ, и скалярный множитель λ, уравнение Q(λξ + η) = 0 имеет м настоящие корни λ1, λ2,… λм. Система строго гиперболический если эти корни всегда различны. Геометрическая интерпретация этого условия такова: характерная форма Q(ζ) = 0 определяет конус (нормальный конус) с однородными координатами ζ. В гиперболическом случае этот конус имеет м листы, а ось ζ = λξ проходит внутри этих листов: он не пересекает ни один из них. Но при смещении от начала координат на η эта ось пересекает каждый лист. В эллиптическом случае нормальный конус не имеет реальных листов.


Аналитические решения

Разделение переменных

Линейные уравнения в частных производных могут быть сведены к системам обыкновенных дифференциальных уравнений с помощью важной техники разделения переменных. Этот метод основан на характеристике решений дифференциальных уравнений: если можно найти какое-либо решение, которое решает уравнение и удовлетворяет граничным условиям, то оно то решение (это также относится к ODE). Мы предполагаем в качестве анзац что зависимость решения от параметров пространства и времени может быть записана как произведение членов, каждый из которых зависит от одного параметра, а затем посмотреть, можно ли это сделать для решения проблемы.[3]

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

Это возможно для простых PDE, которые называются разделимые дифференциальные уравнения в частных производных, а область обычно представляет собой прямоугольник (произведение интервалов). Разделимые PDE соответствуют диагональные матрицы - думая о «стоимости фиксированного Икс«как координату, каждую координату можно понимать отдельно.

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

Метод характеристик

В особых случаях можно найти характеристические кривые, на которых уравнение сводится к ОДУ - изменение координат в области для выпрямления этих кривых позволяет разделить переменные и называется метод характеристик.

В более общем плане можно найти характерные поверхности.

Интегральное преобразование

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

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

Если область конечна или периодична, бесконечная сумма решений, таких как Ряд Фурье уместно, но интеграл таких решений, как Интеграл Фурье обычно требуется для бесконечных доменов. Приведенное выше решение для точечного источника уравнения теплопроводности является примером использования интеграла Фурье.

Замена переменных

Часто PDE можно преобразовать в более простую форму с помощью известного решения с помощью подходящего замена переменных. Например, Блэк – Скоулз PDE

сводится к уравнение теплопроводности

заменой переменных (подробнее см. Решение уравнения Блэка-Шоулза на Wayback Machine (архивировано 11 апреля 2008 г.))

Фундаментальное решение

Неоднородные уравнения[требуется разъяснение ] часто может быть решена (для PDE с постоянным коэффициентом всегда решается) путем нахождения фундаментальное решение (решение для точечного источника), затем взяв свертка с граничными условиями для получения решения.

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

Принцип суперпозиции

Принцип суперпозиции применим к любой линейной системе, включая линейные системы УЧП. Распространенной визуализацией этой концепции является взаимодействие двух синхронизированных по фазе волн, которые объединяются для получения большей амплитуды, например грех Икс + грех Икс = 2 греха Икс. Тот же принцип можно наблюдать в PDE, где решения могут быть действительными или сложными и аддитивными. суперпозиция Если ты1 и ты2 являются решениями линейных уравнений в частных производных в некотором функциональном пространстве р, тогда ты = c1ты1 + c2ты2 с любыми константами c1 и c2 также являются решением этого УЧП в том же функциональном пространстве.

Методы для нелинейных уравнений

Нет общеприменимых методов для решения нелинейных уравнений в частных производных. Тем не менее, результаты существования и уникальности (например, Теорема Коши – Ковалевского ) часто возможны, как и доказательства важных качественных и количественных свойств решений (получение этих результатов является важной частью анализ ). Вычислительное решение нелинейных УЧП, двухступенчатый метод, существуют для конкретных уравнений, таких как нелинейное уравнение Шредингера.

Тем не менее, некоторые методы можно использовать для нескольких типов уравнений. В час-принцип это самый мощный метод решения недоопределенный уравнения. В Теория Рикье-Жане эффективный метод получения информации о многих аналитических сверхопределенный системы.

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

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

Метод группы Ли

С 1870 г. Софус Ли Работа поставила теорию дифференциальных уравнений на более удовлетворительную основу. Он показал, что теории интеграции старых математиков могут, введя то, что сейчас называется Группы Ли, быть отнесенными к общему источнику; и те обыкновенные дифференциальные уравнения, которые допускают то же бесконечно малые преобразования представляют сопоставимые трудности интеграции. Он также подчеркнул тему трансформации контакта.

Общий подход к решению УЧП использует свойство симметрии дифференциальных уравнений, непрерывную бесконечно малые преобразования решений к решениям (Теория лжи ). Непрерывный теория групп, Алгебры Ли и дифференциальная геометрия используются, чтобы понять структуру линейных и нелинейных дифференциальных уравнений в частных производных для генерации интегрируемых уравнений, чтобы найти ее Слабые пары, операторы рекурсии, Преобразование Бэклунда и, наконец, нахождение точных аналитических решений для PDE.

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

Полуаналитические методы

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

Численные решения

Три наиболее широко используемых численные методы решения УЧП являются метод конечных элементов (FEM), методы конечных объемов (FVM) и методы конечных разностей (FDM), а также другие методы, называемые Meshfree методы, которые были созданы для решения проблем, в которых вышеупомянутые методы ограничены. FEM занимает видное место среди этих методов, особенно его исключительно эффективная версия более высокого порядка. hp-FEM. Другие гибридные версии методов FEM и Meshfree включают обобщенный метод конечных элементов (GFEM), расширенный метод конечных элементов (XFEM), спектральный метод конечных элементов (SFEM), бессеточный метод конечных элементов, разрывной метод конечных элементов Галеркина (DGFEM), Безэлементный метод Галеркина (EFGM), Интерполяция безэлементного метода Галеркина (IEFGM) и др.

Метод конечных элементов

Метод конечных элементов (FEM) (его практическое применение, часто известное как анализ конечных элементов (FEA)) представляет собой численный метод для поиска приближенных решений уравнений в частных производных (PDE), а также интегральных уравнений. Подход к решению основан либо на полном исключении дифференциального уравнения (задачи установившегося состояния), либо на преобразовании УЧП в аппроксимирующую систему обыкновенных дифференциальных уравнений, которые затем численно интегрируются с использованием стандартных методов, таких как метод Эйлера, Рунге – Кутта и т. Д.

Метод конечных разностей

Конечно-разностные методы - это численные методы аппроксимации решений дифференциальных уравнений с использованием конечная разница уравнения для аппроксимации производных.

Метод конечных объемов

Подобно методу конечных разностей или методу конечных элементов, значения вычисляются в дискретных точках сетчатой ​​геометрии. «Конечный объем» относится к небольшому объему, окружающему каждую узловую точку на сетке. В методе конечного объема поверхностные интегралы в дифференциальном уравнении с частными производными, которые содержат член дивергенции, преобразуются в объемные интегралы с использованием теорема расходимости. Затем эти члены оцениваются как потоки на поверхностях каждого конечного объема. Поскольку поток, входящий в данный объем, идентичен потоку, выходящему из соседнего объема, эти методы сохраняют массу конструктивно.

Энергетический метод

Энергетический метод - это математическая процедура, которая может использоваться для проверки корректности начально-краевых задач.[4] В следующем примере энергетический метод используется, чтобы решить, где и какие граничные условия должны быть наложены, чтобы результирующий IBVP был корректным. Рассмотрим одномерный гиперболический УЧП в виде

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

Используя это

где для второго отношения использовалось интегрирование по частям, получаем

Здесь обозначает стандартную L2-норму. для корректности потребуем, чтобы энергия решения не возрастала, т. е. чтобы , что достигается указанием в если и в если . Это соответствует только наложению граничных условий на притоке. Обратите внимание, что корректность допускает рост с точки зрения данных (начальных и граничных), и поэтому достаточно показать, что удерживается, когда все данные установлены в ноль.

Смотрите также

Основные примеры PDE

Типы граничных условий

Различные темы

Примечания

  1. ^ Клайнерман, Серджиу. PDE как единый предмет. GAFA 2000 (Тель-Авив, 1999). Геом. Функц. Анальный. 2000 г., специальный том, часть I, 279–315.
  2. ^ Клайнерман, Серджиу (2008), "Уравнения с частными производными", в Гауэрсе, Тимоти; Барроу-Грин, июнь; Лидер, Имре (ред.), Принстонский компаньон математики, Princeton University Press, стр. 455–483.
  3. ^ Гершенфельд, Нил (2000). Сущность математического моделирования (Перепечатано (с корр.) Под ред.). Кембридж: Cambridge Univ. Нажмите. п.27. ISBN  0521570956.
  4. ^ Густафссон, Бертил (2008). Методы разности высокого порядка для зависимых от времени PDE. Ряды Спрингера в вычислительной математике. 38. Springer. Дои:10.1007/978-3-540-74993-6. ISBN  978-3-540-74992-9.

Рекомендации

дальнейшее чтение

внешняя ссылка