CC BY
19
Серия «Математика»
Том 2 (2009), № 1, С. 221-232
Онлайн-доступ к журналу: http://isu.ru/izvestia
УДК 517.97
Минимаксное оценивание в задаче оптимизации динамики пучков заряженных частиц
Д. А. Овсянников
Санкт-Петербургский государственный университет
Аннотация. Рассматривается задача управления пучком взаимодействующих частиц. Исследуется минимаксный критерий оценки качества динамики пучка. Динамика частиц описывается интегро-дифференциальными уравнениями. Даются представления для вариации оптимизируемого функционала и необходимые условия оптимальности.
Ключевые слова: минимакс, пучок, взаимодействие частиц, оптимизация.
Математические проблемы управления в системах обыкновенных дифференциальных уравнений, включая управление пучками (ансамблями) траекторий, интенсивно исследуются уже многие годы. Особый интерес к задачам управления и оптимизации в последние годы вызван постоянно расширяющейся сферой приложений и, конечно, новыми существенно возросшими возможностями вычислительной техники. В данной статье рассматривается проблема управления пучками заряженных частиц с учетом их взаимодействия. Такие проблемы возникают при проектировании и расчете ускорителей заряженных частиц различного назначения. Специальному рассмотрению эти задачи обязаны и перспективности моделирования работы создаваемой современной электрофизической аппаратуры, в частности, моделирования динамики пучков заряженных частиц и плазмы в электромагнитных полях [6]-[10], [12]—[14]. Усложняющей моделирование особенностью последних является учет взаимодействия частиц. В математическом отношении подобные взаимодействия приводят к исследованию интегро-дифференциальных уравнений. В статье ставится задача оптимального управления ансамблем траекторий системы интегро-
Введение
дифференциальных уравнений по минимаксному критерию. Дается представление вариации функционала на траекториях таких систем и доказываются необходимые условия минимума максимума значений функционала по пучку траекторий. При этом используется техника введения сопряженной системы, позволяющая в дальнейшем строить приближенные методы оптимизации.
1. Постановка задачи
Пусть эволюция пучка частиц описывается уравнениями
^ж/^£ = / (£,ж,и), (1.1)
+ 1^/(£,ж,и) + р ё1уж / (£,ж,и), (1.2)
от дх
где
/ (£,ж,и) = /1 (£, ж, и) + / / (^, х, У*)р(^, У*)^У*, (1.3)
^ МЬ,и
с начальными условиями
ж(0) = жо € Мо С П, р(0,ж) = ро(ж). (1.4)
Здесь £ — время, £ € То = [0, Т], Т — фиксировано, ж — п-мерный вектор фазовых координат, р = р(£, ж) — плотность распределения частиц в фазовом пространстве в момент £, и = и(£) — г-мерная вектор-функция управления, Мо — компакт в Л” ненулевой меры, П — область в Л”, ро(ж)
— неотрицательная непрерывно-дифференцируемая функция, множество М*,щ есть сечение в момент £ пучка траекторий, исходящих из
множества Мо и соответствующих управлению и = и(£), т. е.
М*,„ = (ж* = ж(£, Жо, и) : Жо € Мо}, у* € М*,«. (1.5)
Считаем, что допустимые управления и = и(£) составляют класс О векторных кусочно-непрерывных на То функций со значениями в компактном множестве V € Ят. Предполагаем, что вектор-функции /1 (£, ж, и) и /2(£,ж,у) определены и непрерывны по совокупности своих аргументов на множествах То х П х V и То х П х П соответственно вместе с частными производными по ж и у до второго порядка включительно.
Решением системы (1.1), (1.2) при начальных условиях (1.4) и фиксированном управлении и будет совокупность вектор-функций Ж = ж(£,Жо,и), описывающих пучок траекторий, исходящих из множества Мо и функция р = р(£,ж), определенная на этом пучке траекторий.
Введенная математическая модель управления учитывает взаимодействие частиц в пучке. В выражении (1.3) вектор-функция /1 определяет воздействие внешних полей на частицу, функция /2 — взаимодействие частиц. Следует отметить, что введенную математическую модель динамики заряженных частиц можно рассматривать, как систему
уравнений Власова [4]. Однако, функция /2, выражающая взаимодействие частиц, в данном случае является достаточно гладкой и описывает некоторое сглаженное взаимодействие частицы со всем ансамблем частиц. Такие ’’сглаженные” модели получаются при различных численных решениях уравнения Власова, например, методом крупных частиц [3, 9, 10]. Существование и единственность решений уравнения Власова (классических и обобщенных) рассмотрены в работах [1, 2].
В данном случае можно говорить о классических решениях системы
(1.1), (1.2). Заметим, что решение уравнения (1.1) можно рассматривать независимо от уравнения (1.2). А именно, если в уравнении (1.1) с учетом (1.3) перейти к интегрированию по множеству Мо, то уравнение
(1.1) может быть представлено в виде
^Ж /*
— = / (£,ж,и) = /1 (£, ж, и) + / /2(£,ж(£,жо,и),ж(£,уо,и))ро(уо )^уо-
а£ ,]м0
(1.6)
Можно показать [9], что при сделанных предположениях относительно функции (1.3) решение уравнения (1.6), т. е. совокупность непрерывнодифференцируемых по £ вектор-функций ж = ж(£, Жо, и), Жо € Мо, существует и единственно на некотором отрезке [0, Л]. Далее предполагаем, что эти решения можно единственным образом продолжить на весь отрезок То при любом допустимом управлении и € О.
Существование и единственность классического решения уравнения
(1.2) следует из уравнения
^р(^ и)) = _р(£, ж(£, жо, и))ё1ух / (£, ж(£, жо, и), и(£)) (1.7)
с учетом единственности и непрерывной дифференцируемости ж(£, Жо,и) по начальным данным [11]. Заметим, что уравнения (1.1) и (1.7) являются уравнениями характеристик для уравнения (1.2). Отметим также, что для наших целей вместо уравнения (1.2) можно сразу рассматривать уравнение (1.7), определяющее плотности распределения частиц вдоль всего пучка траекторий. Однако, на наш взгляд, естественнее сначала записать уравнения (1.1), (1.2).
Рассмотрим функционал
I(и) = тах тах м(£, ж, р(£, ж)), (1.8)
характеризующий динамику управляемого процесса. Здесь ТЖ — дис-
кретное множество, состоящее из моментов £€То, =1, N, N фиксировано, ^>(£,ж,р) € С 1(То х П х Л1) — неотрицательная функция.
Используя определение сечения пучка траекторий (1.5), функционал (1.8) очевидно, можно представить также в виде
I(и) = тах тах м(£, ж(£, жо, и), р(£, ж(£, жо, и))). (1.9)
Ь&Их0£М0
Далее рассмотрим задачу минимизации функционала (1.9) по управлениям и € О.
2. Вспомогательные соотношения и леммы
Рассмотрим при каждом £ € ТЖ множество Л*(и), зависящее от управления и определяемое соотношениями
(2.1)
Rt(u) = (Xo : ж0 € Mo, ^>(t, ж0, u), p(t, x(t, ж0, u))) =
= maxxoeMo <^(t,x0,u),p(t,x(t,x0,u)))}.
Введем также множество
T(u) = (5: 5 € TN, ^>(t, ж(5, X0, u), p(5, ж(5, X0, u))) = /(u)}. (2.2)
Лемма 1. Пусть Rt(u) и Rt(u) - множества, определяемые соотношениями (2.1) при всех t € TN и соответствующие допустимым управлениям u = u(t) и u = u(t) + Au(t). Тогда имеет место соотношение
max min IIX) — ж0 II —► 0 (2.3)
xO'eRt(u) x'o€Rt(«)
при ||Au||l ^ 0 равномерно по t € TN.
Здесь и далее под ||Au||l ^ 0 понимаем, что u ^ u по норме в пространстве L[0,T].
Доказательство. Соотношение (2.3) означает, что по любому е > 0 найдется 5 > 0, такое, что при ||Au||l < 5 для произвольных t € TN и ж" € Rt(u) можно указать Х0 € Rt(u), при котором ||X) — x01 < е. Предположим противное. Тогда найдется последовательность допустимых управлений, такая, что ||Auk || < 5k, где 5k > 0 и 5k —► 0. При этом
k^-ж
найдутся такие Zk € Rt(u + Auk) при некотором t € TN, что
min ||zk — ж0|| > е. (2.4)
xo eRt(u)
По определению имеем ^>(t, x(t, zk, u), p(t, x(t, zk, u))) = max <^(t, x(t, ж0, u), p(t, x(t, ж0, u))),
xoeMo
(2.5)
где u = u(t) + Auk (t). Все точки Zk € M0. Поэтому, не ограничивая общности, будем считать последовательность (zk} сходящейся к точке z0 € M0. Пусть x(t) = x(t,x0,u); Ax(t) = x(t,x0,u) — x(t,x0,u);
Ap(t, x(t)) = p(t,X(t)) — p(t,x(t)). Можно показать [9], что
||Ax(t)||c = ||x(t, x0,u + Auk) — x(t,x0,u)||c ^ 0, (2.6)
||Ap(t,x(t))||c = Ip(t,x(t,x0,u + Auk)) — p(t,x(t,x0,u))||c ^ 0 (2.7)
равномерно по Ж0 € M0 при ||Auk ||l ^ 0. Учитывая это, перейдем в равенстве (2.5) к пределу при к ^ +то. Получим
^>(t,z0,u),p(t,x(t,z0,u))) = max <^(t,x(t,ж0,u),p(t,x(t,ж0,u))),
xoeMo
что означает принадлежность z0 множеству Rt(u) и противоречит тем самым неравенству (2.4). □
Лемма 2. Пусть T(u) и T(u + Au) - множества, определяемые соотношениями (2.2). Тогда существует 5 > 0, такое, что при ||Au||l < 5 множества T(u) и T(u + Au) совпадают, т. е. T(u) = T(u + Au).
Доказательство. Допустим противное. Тогда для любой последовательности положительных чисел (5j}, стремящейся к нулю при j ^ то, можно указать последовательности (tj} и (x0j}, такие, что
tj € T(u + Auj), x0j € (u + Auj), tj € T(u), (2.8)
где ||Auj ||l < 5j. Не ограничивая общности, считаем, что последовательность (tj} сходится к 5, а последовательность (x0j} — к Х0. В силу конечности множества TN и сходимости последовательности (tj}, начиная с некоторого номера к, имеем tj = tj+1 = t, j = к, к + 1,.... Отсюда и из (2.8) следует, что 5€T(u).
Кроме того, по определению
(tj,ж (tj, x0j, u), p(tj, x0j, u)) =
= max max <^(t, x(t, ж0, u), p(t, x(t, ж0, u))) = /(u) =
teTN xoeRt(U)
= max max w(t, x(t, ж0, u), p(t, x(t, ж0, u))),
teTN xoeMo V ' V V
где u = u+Auj. Переход в этом равенстве к пределу с учетом (2.6), (2.7) дает ^>(5, ж(5, Ж0, u)), р(5, X0,u) = /(u), а это означает по определению
(2.2), что 5 € T(u). Получили противоречие. □
Рассмотрим множества
Qt(u, Au) = Rt(u) U Rt(u + Au), S(u, Au) = T(u) U T(u + Au). Используя их, получим
/(u) = max max w(t, ж(1 ж0, u), p(t, ж(1 ж0, u))) =
V ' teTN xoeMo V V ' V V
= max max w(t, ж(1 ж0, u), p(t, ж(1 ж0, u))) = (2 9)
teT(U) xo€Rt(«) (2.9)
= max max w(t, ж(1 ж0, u), p(t, ж(t, ж0, u))).
teS(u,Au) xo€Qt(u,Au)
Представим функцию <^(£, ж(£, жо, и), р(£, ж(£, жо, и))) в виде ^>(£- ж(£) + Дж(£), р(£, ж(£)) + Др(£- ж(£))) =
= ж((), р((. ж(())) + «Ж(()+
+ д^-ж(*д р(М‘(*») ^р(£ ж(/)) + о((,.Жо, ||Дж(*)||с + ||Др(*,ж(())||с).
др
(2.10)
Здесь и далее вариация дж(£) = дж(ж*) = дж(£,ж(£)), ё1удж = ё1уХ(дж(ж*) и
Ди/1(£-ж(£),и(£)) = /1(£,ж(£),и(£)) _ Д(£- ж(£)и(£)), Д«й1ух/1(£, ж(£), и(£)) = ё1Ух Ди Л(£,ж(£), и(£)).
Вариация дж и дивергенция дж удовлетворяют уравнениям [9]
$дж д/ (£, ж(£), и(£))
дЖ "^Ж + Ди /1 (£-ж(£)-и(£))+
, [ д/2(£-ж(£),у*)
+ / -------«-------дж(у*)р(£,у*)$у*,
•/м4,„ ду
Йё1удж дё1ух/ (£, ж(£), и(£)) ^ . . . , „ , . . ...
—^ =------------хМ дд дж(£) + Ди ё1ух/1(£,ж(£),и(£)) +
Г дё1ух/2(£,ж(£),у*)
+ ---------о-----------дж(у*)р(£,у* )$у*-
./м*,м ду
(2.11)
(2.12)
с начальными условиями
дж(0) = 0, ё1удж(0) = 0. (2.13)
При этом др = др(ж(£)) удовлетворяет уравнению
^р ^тс-
= —драгух/ (£,ж(£)и(£)) _ р—агудж
или
$
— (др + рё1удж) = —(др + рё1удж)а1ух/.
Отсюда с учетом др(0) = 0 следует
др + р ё1ухдж = 0. (2.14)
Введем далее для краткости следующие обозначения: у*(£, жо,и) = (ж*(£, жо, и), р(£, жо, и)), ду*(£,уо) = (дж*(£, жо), др(ж(£, жо, и))).
МИНИМАКСНОЕ ОЦЕНИВАНИЕ В ЗАДАЧЕ ОПТИМИЗАЦИИ Из равенства (2.9) с учетом (2.10) имеем
/(u) < max max <^(t, y(t, ж0, u)) +
teS(u,Au) xoeQt(u,Au)
d^(t, y(t, ж0 ,u))
+ max max -------------------------- --------5y(t, ж0) + o1 =
teS(u,Au) xoeQt(u,Au) dy
= max max w(t, y(t, ж0, u)) +
teT(u) xoeRt(u)
d^(t, y(t, ж0, u)) _ , .
+ max max ------------------ -----5y(t, ж0)+ o1 + o2,
teT(u) xoeRt(u) dy
где
o1 = max max |o(t, ж0,11Aж|c + ||Ap||c)|,
te[0,T] xoeMo
d^(t, y(t, ж0, u))
02 = max max ------------------------- ------------5y(t,жo) —
teS(u,Au) xoeQt(u,Au) dy
d^(t, y(t, ж0, u))
— max max ---------------- ----------5y(t, ж0).
teT(u) xoeRt(u) dy
Используя (2.10), (2.16) получаем
/(u) > max max < w(t, y(t, ж0, u))+
V ' _ teTN xo eMo ' V V ''
+ dy(t, y(t, ж0, u))
dy — 01 >
> max max < w(t, y(t, ж0, u))+ teT(u) xo eRt(u) \
5у((,ж0)} — o, =
+
dy
Поскольку
= max max w(t, y(t, ж0, u)) +
teT(u) xoeRt(u)
d^(t, y(t, ж0, u))
+ max max -------------------- -----5y(t, ж0) — o1.
teT(u) xoeRt(u) dy
I(u) = max max <^(t, y(t, ж0, u)),
teT(u) xoeRt(u)
то из соотношений (2.15) и (2.18) следует
5/ — o1 < A/ < 5/ + o1 + o2, где A/ = /(u) — /(u) — полное приращение функционала (1.9);
5/ = max max 5у((,ж0 ).
teT(u) xoeRt(u) dy
(2.15)
(2.16)
(2.17)
(2.18)
(2.19)
(2.20)
Преобразуем правую часть равенства (2.20). В связи с этим при I € Т(и)Жо € Л*(и) рассмотрим
^(?,Жо) = У(дуХ0’и))ду(?,Жо). (2.21)
Введем вспомогательные вектор-функции ф(?ж) и ф(?), используя уравнения (2.11), (2.12) и равенства (2.13), (2.14), и преобразуем £^>. Обозначим Ж(?) = ж(?Жо’и), р (?) = р(?Жо(?)). Пусть
Л(М’Жо) = Л(?) = Л = д^(?’Ж®’р (^))р(^)’ ? € [0’4]’ (2.22)
а вектор-функции ф(?) и ф(? ж) удовлетворяют на пучке траекторий, соответствующих управлению и = и(?), уравнениям
#=_ ( д/(^’Ж(^)’и(^)) \*ф -/ д (^’Ж(^)’и(^)) \* (2 23)
V дЖ / V дЖ / .
^ф(?ж(?)) ( д/ (£’Ж(£)’и(£)) \
ф—
—
^ \ дж
д/2(£’Ж(£)’Ж(£))\* - (дdivx/2(? Ж(?)’ж(?)) 4 *
<9ж у \ дж
Г (д/2(^’У*’Ж*)
- (-----> Л- (2.24)
ф(?у)р(?у)%
•М.и V дж
при конечных условиях
ф -(у ’ (2.25)
дЖ
ф(?ж) = 0’ ж € (2.26)
Вариация (2.21) с учетом (2.22)—(2.26) преобразуется к виду
^>(?Жо) = - / [ф*(?)Дй/1(?’Ж(?)’и(?)) + ЛД^^ж/1 (^’Ж(^)’и(^))]^^-
о
Пф*(? ж*) Дщ/1 (^’ ж*’ и(?))р(? ж*)^ж^.
Дь,и
(2.27)
Заметим, что по построению функции ф(?) и Л(?) зависят от Жо и определены на интервале [0’?]. Подчеркивая это, будем записывать ф(?) = ф(?’?’Жо), Л(?) = Л(?’?Жо). Аналогично ф(?ж) определена на сечениях М*,щ при ? € [0’Л. Поэтому будем также использовать запись ф(?’?ж) для вектор-функции ф(?’Ж). Для упрощения записи в дальнейшем доопределим указанные выше функции на весь интервал [0’ Т]
следующим образом:
^(і,ЇЖ0) = 0 при і > і ;
Л(і, і ж0) = 0 при і > і;
ф(і, і ж) Е =0 при і > ?.
3. Необходимые условия оптимальности.
Используя представление вариации в виде (2.27) с учетом (2.19)— (2.21), докажем следующую теорему.
Теорема 1. Пусть и0 = и0 (і) — оптимальное управление, тогда при
і Є [0, Т] выполняется соотношение
тіп тах тах (ф0*(іД ж0)Ди°/і(і,ж0(і),и) +
«єи ієТ(и0) ж0ЄЯ((и°)
+Л0(і,і,Ж0)Ди0 div/l(і,x0(і),u)+ (3 1)
+ / ф0*(і,іу )Д«° /і(і,у,и)р0(і,у )гіу } = 0.
Здесь Жо(?) = ж(?’Жо’и); ро(?’Ж) соответствует управлению ио, а фо, Ло, фо удовлетворяют (2.22)-(2.26) на оптимальном процессе, т.е. на решении системы (1.1), (1.2) при и = ио(?).
Доказательство. Доказательство проводится стандартным образом с использованием игольчатой вариации
Д (?) = Iи - ио(?)’ ? € [0’ЛП [£’£ + 5)’
Й [0’ ? € [£’ £ + 5).
Очевидно, с учетом неравенств (2.19) и представления (2.27) достаточно показать, что величины 01 и 02, определенные равенствами (2.16), (2.17), имеют более высокий порядок малости, чем 5 из игольчатой вариации. Нетрудно показать, что о(?’Жо’ ||Дж||с + ||Др||с) из равенства
(2.10) более высокого порядка малости, чем
НДиЛЬ + НД^^ДН^
и, следовательно, более высокого, чем 5. Из непрерывной зависимости о(?’Жо’ ||Джс|| + ||Др|с) по ? и Жо, и компактности множеств То и Мо следует, что 01, определенное равенством (2.10) так же более высокого порядка малости, чем 5.
Мі.и°
Рассмотрим теперь величину 02. Отметим прежде всего, что вектор-функция д^>(?’ у(? Жо’ и))/ду непрерывна по совокупности переменных (?’ Жо) € То хМо и, следовательно, равномерно непрерывна на множестве То х Мо. Фиксируем е > 0. Тогда найдется 7 > 0, такое, что при ||жо — жон <7 выполняется
д^(?’У(?’Жо'’и)) д^(?’У(?’Жо’и))
ду
ду
<е
для всех ? € То и жо^о7 € Мо. По полученному 7 выберем согласно лемме 1 5 > 0 так, что при любом 5 < 5 для всякого жо; € Щ (ио +
Ди<5), ? € Т(ио), можно найти жо € Щ(и), при котором Ж — жо| <7. Будем считать 5 достаточно малым, таким, чтобы при 5 < 5 имело место утверждение леммы 2, а именно Т(ио) = Т(ио + Ди). Фиксируем некоторое 5 < 5 . Пусть 5 и Жо; таковы, что
тах
тах тах 5у((.жо)
д^(5 ’ у(5 ’Жо'’и))
ду
-5у(5 ’Жо).
В силу непрерывности 5у(?’ Жо) по Жо выбираем 7 так, что |5у(5’Жо') — 5у(5’Жо)| < е|5у(5’Жо')|.
Тогда, учитывая, что Т(ио) = 5(ио’ Ди), Щ(ио) С (и°’ Ди), имеем
д^(5 ’у(5 ^и))
д^(5 ’у(5 ’ жо ’ и))
5у(5’Жо) <
<
д^(5 ’у(5 ^и)) д^(5 ’у(5 ’ жо ’ и))
ду
д^(5 ’ у(5 ’Жо’и))
|5у(5 ’401+
ду
ду
|5у( 5’жо/)—5у( 5’жо)|.
В силу ограниченности
д^(? ’у(? ’Жо’и))
ду
< К’ (?’Жо) € То х Мо’
и предыдущих неравенств получим од < е(К + 1)|5у(5’Жо/)|, а так как |5у(?’ Жо)| имеет порядок малости, не меньший 5, то од = о(5), и, следовательно, 02 имеет такой же порядок малости.
Отметим также, что в точках разрыва управления ио = ио(? ) в соотношении (3.1) необходимо рассматривать соответствующие односторонние пределы. □
4. Заключение
На основе полученной вариации оптимизируемого функционала в виде (2.27), выбирая тот или иной способ варьирования управления (см. например [5]), можно получить необходимые условия оптимальности и в другой форме. Отметим также, что исследование динамики частиц с учетом их взаимодействия привело к необходимости введения интегро-дифференциального уравнения (2.24). Таким образом, в отличие от случая /2 = 0, не учитывающего взаимодействие частиц [8], необходимые условия оптимальности дополнились слагаемыми, содержащими распределенное влияние пучка через сопряженную переменную ф(£,ж). Однако, общая схема построения направленных методов оптимизации с использованием полученных условий оптимальности и вариации функционала, очевидно, сохраняются [8, 9].
Список литературы
1. Арсеньев А. А. Единственность и существование в малом классического решения системы уравнений Власова / А. А. Арсеньев // Докл. АН СССР. —
1974. — Т. 218, № 1. — С. 11—12.
2. Арсеньев А. А. Существование в целом слабого решения системы уравнений Власова / А. А. Арсеньев // Журн. вычисл. математики и мат. физики. —
1975. — Т. 16, № 1. — С. 136-147.
3. Березин Ю. А. Метод частиц в динамике разреженной плазмы / Ю. А. Березин, В. А. Вшивков. — Новосибирск, 1980. — 96 с.
4. Власов А. А. Теория многих частиц / А. А. Власов. — М.; Л., 1950. — 348 с.
5. Демьянов В. Ф. Негладкие задачи теории оптимизации и управления / В. Ф. Демьянов, Т. К. Виноградова, В. Н. Никулина и др. — Л.: Изд-во Ленингр. ун-та, 1982. — 324 с.
6. Дривотин О. И. О самосогласованных распределениях для пучка заряженных частиц в продольном магнитном поле / О. И. Дривотин, Д. А. Овсянников // Докл. РАН. — 1994. — Т. 33, № 3. — С. 284-287.
7. Овсянников А. Д. Управление пучком заряженных частиц с учетом их взаи-
модействия / А. Д. Овсянников // Вестн. С.-Петерб. ун-та. Сер. 10. — 2009. — Вып. 2. — С. 81-91.
8. Овсянников Д. А. Математические методы управления пучками / Д. А. Овсянников. — Л.: Изд-во Ленингр. ун-та, 1980. — 228 с.
9. Овсянников Д. А. Моделирование и оптимизация динамики пучков заряженных частиц / Д. А. Овсянников. — Л.: Изд-во Ленингр. ун-та, 1990. — 312 с.
10. Овсянников Д. А. Моделирование интенсивных пучков заряженных частиц / Д. А. Овсянников, О. И. Дривотин. — СПб.: Изд-во С.-Петерб. ун-та, 2003. — 174 с.
11. Степанов В. В. Курс дифференциальных уравнений / В. В. Степанов. — М.; Л.: Гостехиздат, 1945. — 406 с.
12. Ovsyannikov D. A. Modeling and Optimization Problems of Charged Particle Beams Dynamics / D. A. Ovsyannikov // Proc. of the 4th European Control Conference. Brussels, 1997. — P. 390-394.
13. Ovsyannikov D. Proceeding of the International Computational Accelerator Physics Conference (ICAP 2004). - St. Petersburg, Russia, 2004 / D. Ovsyannikov, M. Berz, K. Makino // Special Issue of Nuclear Instruments and Methods in Physics Reseach A. - 2006. - Vol. 558. - № 1.
14. Ovsyannikov D. Proceeding of the 15th International Workshop on Beam Dynamics and Optimization (BDO 2008). - St. Petersburg, Florida, USA, 2008 / D. Ovsyannikov, M. Berz, K. Makino // Special Issue of International Journal of Modern Physics A. — 2009. — Vol. 24. — № 5.
D. A. Ovsyannikov
Minimax in the optimization problem of charged particle beam dynamics
Abstract. The control problem of charged particle beam with taking into account their interaction is considered. Minimax criteria is investigated. Particle dynamics is described by a system of integro-differential equations. Functional variation and necessary optimality conditions are presented.
Keywords: minimax, beam, interaction of particles, optimization.
Овсянников Дмитрий Александрович, доктор физ.-мат. наук профессор, Санкт-Петербургский государственный университет, 198508 СПб, Старый Петергоф, Университетский пр.35, тел.: (812) 428-47-29, (dovs@mail.ru)
Ovsyannikov Dmitry Alexandrovich, Doctor of physical & mathematical sciences, Professor, Saint-Petersburg State University, 198504 SPb, St. Pe-tergoff, Universitetskij pr.35, Phone: (812) 428-47-29, (dovs@mail.ru)
