Впрочем, даже имея в своем распоряжении все необходимые инструменты, мне предстояло проделать немалую подготовительную работу. Частично это было обусловлено тем, что никто до меня не решал комплексные уравнения Монжа-Ампера для случая более чем одного измерения. Как альпинист, постоянно стремящийся к покорению новых высот, я стремился к покорению более высоких размерностей. Чтобы подготовить себя к схватке с многомерным уравнениям Монжа-Ампера, нелинейность которых сама собой подразумевалась, мы с моим другом Ш. Ченгом принялись за исследование различных многомерных случаев, начав с задач в вещественных числах с целью впоследствии перейти к более сложным комплексным уравнениям.
Для начала мы рассмотрели знаменитую задачу, выдвинутую на рубеже XX века Германом Минковским. Задача Минковского состояла в том, чтобы установить возможность или невозможность существования некоей структуры, удовлетворяющей определенному набору критериев. Рассмотрим простой многогранник. Его структуру можно охарактеризовать, подсчитав число граней и ребер и определив их размеры. Задача Минковского состояла в обратном: можно ли, зная форму, площадь, число и ориентацию граней, определить, существует ли в действительности многогранник, удовлетворяющий данным критериям, и если да, то будет ли он единственным?
Задача на самом деле была более общей, поскольку имела отношение не только к многогранникам, но и в принципе к любым выпуклым поверхностям. Вместо того чтобы говорить об ориентации граней, с равным успехом можно говорить о кривизне, указав для каждой точки поверхности направление перпендикулярных к ней — нормальных векторов, что соответствует ориентации поверхности в пространстве. После этого уже можно задаться вопросом, существует ли объект с указанной кривизной.
Удобно, что эта задача может быть представлена не только в геометрической форме. Она также может быть записана в виде дифференциального уравнения в частных производных. По словам Эрвина Лутвака из Политехнического института при Нью-Йоркском университете: «Если вы сможете решить геометрическую задачу, то автоматически получите дополнительный приз: решение сложнейшего дифференциального уравнения в частных производных. Такая взаимосвязь между геометрией и дифференциальными уравнениями в частных производных делает эту задачу столь важной».[49]
Мы с Ченгом нашли способ решения этой задачи, и наша статья, посвященная этому вопросу, вышла в 1976 году. Как выяснилось, другое решение было представлено несколькими годами раньше — в 1971 году российским математиком Алексеем Погореловым. Ни я, ни Ченг никогда не видели его статьи, поскольку она была опубликована на родном языке Погорелова. В конце концов, все свелось к решению сложнейшего дифференциального уравнения в частных производных из тех, которые никогда до этого не решались.
Несмотря на то что никому до нас не удавалось решить проблему данного типа, за исключением Погорелова, работа которого была нам неизвестна, процедура, позволяющая работать с нелинейными дифференциальными уравнениями в частных производных, на тот момент была уже хорошо разработана. Метод работы с подобными уравнениями, получивший название метода непрерывности, был основан на использовании последовательных приближений. И хотя этот общий подход никоим образом нельзя было назвать новым, особенность состояла в том, что каждая конкретная задача предусматривала разработку своей собственной стратегии, необходимой для ее решения. Основная идея заключалась в последовательной аппроксимации решения различными функциями так, чтобы каждое следующее приближение давало результаты лучше, чем предыдущее. Суть доказательства состояла в том, чтобы показать, что после достаточно большого числа итераций приближенная функция с большой точностью совпадет с решением искомого дифференциального уравнения. В случае удачи полученное путем аппроксимации приближение нужно рассматривать не как решение дифференциального уравнения, которое можно представить в виде определенной формулы, а как доказательство того факта, что решение существует. Для гипотезы Калаби и других задач того же типа существование решения дифференциального уравнения в частных производных эквивалентно доказательству существования определенной геометрии для заданных «топологических» условий. Это не означает, что вы ничего не знаете о решении, существование которого только что доказали. Схема, которая была использована для доказательства существования решения, зачастую может быть легко преобразована в численный метод для приближенного решения на компьютере. О численных методах речь пойдет в девятой главе.
Рис. 5.1. Математик Ш. Ю. Ченг (фотография Джорджа М. Бергмана)
Метод непрерывности был назван так потому, что он подразумевает непрерывное преобразование решения некоего известного уравнения вплоть до его полного совпадения с решением искомого. Процедуру преобразования, как правило, разбивают на две части, одна из которых работает только в непосредственной близости от известного решения.
Рис. 5.2. Наглядная иллюстрация метода Ньютона. Для того чтобы найти точку пересечения определенной кривой или функции с осью X, сначала нужно наугад подобрать некую точку x0 наиболее подходящую для этого. Затем необходимо провести касательную к кривой в точке x0 и отметить точку, в которой эта касательная пересечет ось X (это будет точка x1). В том случае, если наше изначальное предположение не было полностью ошибочным, продолжая этот процесс, мы будем получать точки все ближе и ближе к искомой
Одна из этих частей носит название метода Ньютона, так как она в определенной степени основана на методе, разработанном Исааком Ньютоном более трехсот лет назад. Для того чтобы продемонстрировать этот метод в действии, рассмотрим функцию y=x3-3x+1, которая описывает кривую, пересекающую ось X в трех различных точках, являющихся корнями этого полинома. Подход, предложенный Ньютоном, позволяет определить положение корней на оси X, что далеко не всегда можно сделать, просто взглянув на уравнение. Предположим, что напрямую решить уравнение нельзя, однако один из корней соответствующей функции можно найти вблизи точки x1. Касательная, проведенная к кривой в этой точке, пересечет ось X в другой точке — x2, находящейся ближе к искомому корню, чем точка x1. Если мы проведем касательную в точке x2, она пересечет ось X в точке x3, которая будет еще ближе к искомому корню. Таким образом, многократное повторение данной процедуры должно довольно быстро привести нас к искомому корню, если только начальная точка x1 была выбрана более-менее удачно.
В качестве еще одного примера рассмотрим набор уравнений Et только одно из которых, Е0 (для которого t = 0), мы способны решить. При этом в действительности нам нужно решить уравнение E1 (для которого t = 1). Мы могли бы использовать метод Ньютона, если мы находимся в непосредственной близости к точке t = 0, решение уравнения в которой хорошо известно, но этот подход не может привести нас к 1. В этом случае необходимо прибегнуть к другому методу оценки, обладающему большей применимостью.
Как же это сделать? Представим, что над Тихим океаном была запущена ракета, которая приземлилась в радиусе ста миль от атолла Бикини. Это дает нам некоторое представление о том, где ракета может быть, другими словами — ее общую позицию, но мы хотели бы знать больше, например ее скорость, или ее ускорение, или как это ускорение изменялось в течение полета. Это можно сделать при помощи дифференциального исчисления — путем взятия первой, второй и третьей производных от функции, описывающей зависимость положения ракеты от времени. С таким же успехом можно брать производные и более высоких порядков, но для эллиптических уравнений второго порядка того типа, которым я занимаюсь, третьей производной вполне хватает.
Одного лишь знания производных функции недостаточно, хотя задача по их нахождению сама по себе может быть чрезвычайно трудоемкой. Кроме того, производные нужно «контролировать». Иными словами, необходимо установить для них границы — удостовериться, что они не могут быть ни чрезвычайно велики, ни чрезвычайно малы. Только в этом случае полученные решения будут «стабильны» — то есть не будут бесконтрольно раздуваться, тем самым дисквалифицируя себя как решения и разрушая наши надежды на них. Итак, взяв для начала нулевую производную — то есть исходную функцию, описывающую изменение положения ракеты с течением времени, мы устанавливаем для нее наличие верхних и нижних границ — иными словами, делаем оценки, показывающие, что решение по крайней мере возможно. Та же самая операция проводится для всех производных более высоких порядков, что позволяет удостовериться в том, что они не являются ни бесконечно большими, ни бесконечно малыми, а функции, их описывающие, не флуктуируют совершенно беспорядочным образом. Это позволяет априори оценить скорость, ускорение, зависимость ускорения от времени и т. д. Если мы можем таким образом проверить все производные от нулевой до третьей, значит, у нас есть хороший способ оценить уравнение в целом и приличный шанс найти его решение. Подобный процесс оценки и доказательства того, что оценочные данные сами по себе находятся под контролем, как правило, представляют самую сложную часть всего процесса.