Исходя из условия стационарности функционала запишем вариацию
функционала в следующем виде:
![](pics/va/45_.gif) |
(4) |
решение модели
варьируется только по начальным данным.
Интегрируя
по частям с учетом краевых условий и
выбирая пока ещё произвольные множители Лагранжа как решение следующей
сопряженной задачи:
![](pics/va/45_1.gif) |
(5) |
Сопряженная задача является формально сопряженной к задаче (1) в смысле
тождества Лагранжа, т.е. оператор сопряженной задачи является
сопряженным к оператору основной задачи.
Исходя из определения вариации функционала: |
![](pics/va/6_1.gif) |
Получаем
достаточно простое выражение для градиента функционала по начальным
данным:
![](pics/va/7.gif) |
(7) |
Этот градиент
оказывается решением построенной нами сопряженной задачи. Следующее
приближение начальных данных получим путем итерационного спуска в
направлении этого градиента функционала по формуле:
![](pics/va/8.gif) |
(8) |
В процессе
итераций мы инициализируем (восстанавливаем) такое начальное поле при
котором решение нашей модели хорошо соответствует данным измерений.
Алгоритм
идентификации в
трехмерном
случае является аналогичным одномерному.
|