Приближенное решение нелинейного параболического и обыкновенного дифференциального уравнений и приближенный расчет функционала качества при известных управляющих воздействиях Текст научной статьи по специальности «Математика»

атематические проблемы управления

УДК 519.3:62-50

ПРИБЛИЖЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНОГО ПАРАБОЛИЧЕСКОГО И ОБЫКНОВЕННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЙ И ПРИБЛИЖЕННЫЙ РАСЧЕТ ФУНКЦИОНАЛА КАЧЕСТВА ПРИ ИЗВЕСТНЫХ УПРАВЛЯЮЩИХ ВОЗДЕЙСТВИЯХ

Т.К. Юлдашев

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

Ключевые слова: параболическое уравнение, начальные и граничные условия, оптимальное управление, обобщенная разрешимость, приближенное решение, минимизация функционала.

ВВЕДЕНИЕ

Развитие теории оптимального управления связано с ростом требований к быстродействию и точности систем регулирования. На основе математической теории оптимального управления разработаны способы построения оптимальных по быстродействию систем и процедуры аналитического конструирования оптимальных регуляторов. Современные методы решения задач управления в значительной степени основаны на концепции оптимальности, что определяет широкое применение методов и алгоритмов теории оптимизации при проектировании и совершенствовании систем управления. Многие задачи управления формулируются как конечномерные оптимизационные задачи. К таким задачам, в частности, относятся и задачи адаптивных систем управления [1—5].

Сложность задач теории оптимального управления потребовала более широкой математической базы для ее построения. Теория оптимального управления для систем с распределенными пара-

метрами, описываемых уравнениями с частными производными, стала разрабатываться сравнительно недавно. За короткое время она получила бурное развитие, все шире проникая в различные области техники и технологических процессов. К системам с распределенными параметрами относятся задачи аэрогазодинамики, химических реакций, диффузии, фильтрации, процессов горения, нагрева и др. [6—11].

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

применяется широким спектр разных методов [12-14].

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

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

1. ПОСТАНОВКА ЗАДАЧИ

Рассмотрим нелинейную задачу управления процессом распространения тепла по стержню конечной длины

^^^ - ^^ = Д', X) + Ж X, И(*, X), т(1)) (1)

dt

Зх

при смешанных условиях

u(t, x)|t=0 = ф(х), u(t, x)|x=0 = u(t, x)|x =, = 0,

(2) (3)

где f — функция внешнего источника, f(t, x, u, т) е е C(D х R х DT), u(t, x) — функция состояния управляемого процесса, ф(х) — функция распределения тепла по стрежню в начальный момент времени, <p(x)|x=0 = <p(x)lx = i = 0, <p(x) е C3(Di), D - DTх Dp DT = [0, T], Dl = [0, /], 0 < I< ю, 0 < T< ю, P(t, x) — управляющая функция, а функция т(?) определяется как решение нелинейного обыкновенного дифференциального уравнения

т'(1) = h(t, T(t), 9(0) (4)

при начальном условии

т(0) = т0 = const, (5)

9(t) е C(DT) — вторая управляющая функция.

В настоящей работе при фиксированном управлении P(t, x) применяется метод разделения переменных, основанный на поиске решения смешанной задачи (1)—(3) в виде ряда Фурье

да

u(t, x) = X a(i)b(x),

i = 1

где b.(x) = Z2- sinX.x, X. = ^ , i = 1, 2, ...

1 l ii i

Рассмотрим пространства С. Л. Соболева 1 0 1

W2 (D), W2 (D) (см., например, работы [15—17]). Следуя О.А. Ладыженской, рассмотрим обобщенное решение задачи (1)—(3) из класса W1 (D) (см. например, работу [15]).

Определение. Если функция u(t, x) е (D) удовлетворяет интегральному тождеству

T I

00

J J м?, У)

дФ( t, y) _ д2Ф( t, y)

dt

Зу

[P (?, У) +

+ /(t, у, u(t, у), т(?))]Ф(?, у) 1 dydt = J ф[Ф(?, y)]t=ody

О 1

для любого Ф(?, x) е W2 (D), то она называется обобщенным решением смешанной задачи (1)—(3). ♦ Задача. Найти состояния: u *(t, x) — обобщенное решение смешанной задачи (1)—(3) и T*(t) — решение начальной задачи (4), (5) при заданных управляющих функциях P*(t, x) е {P*:|P*(t, x)| < Ql, (t, x) е D}, 9*(t) е {9*:|9*(t)| < Q2, t е DT}, что доставляют минимум функционалу

T

J[P, 9] = jg(t, x, u(t, x), T(t), P(t, x), 9(t))dt,

0

где 0 < Q* = const, i = 1, 2.

2. УКОРОЧЕННОЕ РЕШЕНИЕ СМЕШАННОЙ ЗАДАЧИ (1)—(3)

Во множестве {a(t) = (a/t))k(t) е C(Dt), i = 1, 2, 3, ...}

с помощью нормы

IIa (t )|| b2(T)

£ niax |аг< t)f i = 1 T

1/2

вводится банахово пространство В2(Т). Наряду с ним рассматривается и конечномерное банахово

пространство (T) с нормой

N.

( ? ^ <( T) =

' N

:

i = 1

max I a

t e

N(?)|2

1/2

0

30

Для произвольной функции £(х), х е в пространстве Х2(^) вводится норма

Г 1 1 Х/2 Ш х)11 х2( д) = | £( У )2 ^ | < ю.

Для числовой последовательности фг в пространстве /2 используется норма

U = f £ КГ

г = 1

1/2

< да.

Также рассмотрим и норму

f N ^1/2

N _ f V N2 I ^ Ф /N - | £ Фг I ^ <

Будем пользоваться обозначениями: ) —

класс непрерывных функций, ограниченных по норме с положительным числом М; у } —

класс непрерывных функций, удовлетворяющих условию Липшица по переменным и, V, ... и с положительным коэффициентом X. Для функций одной переменной индекс опускается.

Из работы [19] в частном случае следует, что задача Коши (4), (5) при фиксированных значениях управления 9(?) имеет единственное решение на отрезке Бт, если выполняются условия:

Щ т(?), 9(0) е е С(£тхX хХ2) п ВиА(М0) п Х/р{#0(?)|т},

т

где 0 < | |#0(0|Л < ю; Хх, Х2 — отрезки прямой,

Предполагается также, что Р(?, х) допускает разложение в ряд Фурье по собственным функциям

0 < M0 — const. Это решение может быть найдено методом последовательных приближений

b(x): x) — £ p,.(t)b,.(x), где p,.(t) — JP(t, y)b,.(y)dy.

г = 1 0

Функцию g(t, y, u(t, y), x(t), P(t, y), e(t)) разложим в ряд Фурье:

да T l / да

J[p, е] — U JgI t, y, £ a/Ob/y), T(t),

г = 1 0 0

j = 1

£ P/(t)bj(y), (t) Ib,(y)dydf b(x). (6)

1 = 1

Из работы [18] в частном случае следует, что обобщенное решение смешанной задачи (1)—(3) представимо в виде

г = 1

t I

u(t, x) — £ Wt) + J J G(t, s)

00

P,;(s) +

ds ¡> b.(x), (7)

+ /5, У, £ а/*)Ъ/у), т(5) IЪ,.(у)Ау

V 1 = 1 V .

где аг(?) определяется как решение счетной системы нелинейных интегральных уравнений (ССНИУ):

T0(t) — Т0, T(k + 1)(t) — Т0 + Jh(s, Tk(s), e(s))ds,

0

k — 0, 1, 2

? A ?

(9)

Вместо уравнений (7) и (8) рассмотрим укороченную систему интегральных уравнений

Nt

uN(t, x) — £ foN (t) + J

i = 1

N

Pi (s) +

JV

+ Jf I s, y, £1 j(s) bjN(y), T(s) I bN(y)dy

s GN(t, s)dsl bN(x),

(10)

N / ч

где аг (?) определяется как решение конечной

системы нелинейных интегральных уравнений (КСНИУ):

Л (t) — 4 (t) + J

N

PN (s) +

t l

a,.(t) — o,.(t) + J J G,.(t, s) s

00

Pi(s) + f s, y, £ a/s)b/y), t(s) | bi(y)dy

v j = 1

ds, (8)

22 ю/t) — Ф,ехр{—t}, G(t, s) — exp{—(t - s)},

Фг — J Ф(У)Ь,(УМУ.

+ Jf fs, У, £ j(s) bf (y), T(s) 1 bN(y)dy

j = 1

s gN (t, s)ds,

(11)

где щ (?) = фf ехр{— ?}, (?, 5) = ехр{-(? - 5)}

А начальные данные фN подбираются из условия

N

(2) так, что при N ^ ю сумма ф^х) = £ фг- Ъг- (х)

г = 1

0

0

X

да

да

0

X

0

да

X

0

аппроксимирует функцию p(x) е L2(Dl), а сумма

N

II Nk + 1,л Nk,J| „

II а (t) - а (ОН T) <

Р^, л) = X л (?)Ьт (*) - функцию Р(?, л).

г = 1

Аналогично вместо выражения (6) рассмотрим укороченный функционал

k--+----- 1

< мк +1 M22k +1 l 2 А

illH(s, x)|| L2(D) ds

J^,(15)

IbN(l) , Mik +1 Mf + 1

N T l

N

J[pN, 9] = X jjgIt, y, X aj(t)bj(y), T(t), i = 10 0 1 j = 1

N

j = 1

где М = Це^/, я)||Т), М2 =

= = (М/ + !(М2)2к +

Существование решения КСНИУ (11) следует из справедливости оценок (14) и (15), так как при к ^ да

N

X pj (t)bj(y), 9(t)bj (y) I dydt- bj (x). (12) последовательность функций {а (t)}k = 1 сходится рав-

3. ОДНОЗНАЧНАЯ РАЗРЕШИМОСТЬ КСНИУ (11)

Теорема 1. Пусть выполняются условия:

T

1) А = j||f(t, x, uN, т)|1 п dt < ю, 0 < А = const;

0 2

2) f(t, x, uN, т) е Lip{H(t, x)| N}, где 0 <

номерно по / к функции а^(г) е В^ (Т). Для доказательства единственности решения в пространстве В^ (Т) предположим, что КСНИУ (11) имеет два решения: е В^ (Т) и 3^(0 е В^ (Т). Тогда для их разности справедлива оценка

ИЛ/) - о!вж(Г) *

< j||H(s, x)||i2(щds < ю;

1 t

< MM22 l2 il|H(s,x)||i7(D;) IIaN(s)-0N(s)||BN(t)ds. (16) 0

u2(D[) " Wll^'(t)'

Применение к оценке (16) неравенства Гронуолла — 3) ||ф^|| м < да. Беллмана дает ||аЖ(/) - о!^ Т) = 0 для всех / е Бт. От-

12 2

Тогда при фиксированных значениях управления сюда следует единственность ре0ения КСНИу (11) в

р^ (?) и функции т(?) КСНИУ (11) имеет единственное решение на отрезке Бт.

Доказательство. Рассмотрим итерационный процесс:

пространстве B2 (T). Теорема доказана. ♦

Подставляя КСНИУ (11) в формулу х) =

N

= X а (?)Аг- (х), получаем систему уравнений (10).

г = 1

Если е ^ (Т), то справедлива оценка

а2 (t) = ш2 (t) + Jp; (s) G;( t, s) ds,

N

ti

00

N

i k +'(t) = aN0(t) + JJf I s, y, X ау k(s)bjN(y), т(s) I x

j = 1

(13)

x bN(y) GN(t, s)dyds, k = 0, 1, 2, 3, ..., t e DT,

где ш^ (/) = ф^ ехр{— /}, (/, я) = ехр{-(/ - я)},

мк

а; (/) является к-м приближением КСНИУ (11).

В силу условий теоремы, из формул (13) по индукции получаем оценки

|uN(t, x)| < X |aN(t)| • |bN(x)| < MjlaN(t)||bn(T) < ю. i= 1 2

4. СХОДИМОСТЬ РЯДА ФУРЬЕ (7) ПРИ ИЗВЕСТНЫХ УПРАВЛЯЮЩИХ ФУНКЦИЯХ, В ТОМ ЧИСЛЕ ОПТИМАЛЬНЫХ

Пусть Р*(?, х) и 9*(?) — известные оптимальные управляющие воздействия. Тогда из формул (9)—(12) приходим к итерационным процессам

< J

0

N

X

i = 1

II N1,., N0,J| „

IIа (t) - а (t)||BN(T) <

Jf ^s, У, X ауУ0( s) bjN(y )| bN(y) dy

0 ^ j = 1

1/2

uN(k + 1)*(t, x) = X К (t) + j pjk* (s) + i = 1 I 0 .

+ jf is, y, X afk* (s) bj(y), Tk*(s)I X

0 1 j = 1 )

S j X I Gf(t, s)|2 [ ds < MjM2T/ А, (14)

1/2

x bN (y)dy

(t, s)dsl bN(x),

k

0

а

где

Nk* / ч |

P (s) +

N(k + 1)* ,л N /л _i_ f

a (t) = юг- (t) + J

0

+ Of (s, y, I ajNk* (s) (y), Tk*(s) j s

«1t« = Mo

k +1

j#01( s) ds

k 1 t ■1 + 2 |H02(s)qfc+i(s)ds. (25) Ю. o

В силу оценки (22), из формулы (25) следует, что

s bN(y)dy

(t, s)ds,

lim S1t(0 = 0.

(18)

k ^ ю

(26)

Применяя неравенства Гронуолла — Беллмана к (24), получаем

т(к + 1)*(t) = to + Jh(s, Tk*(s), ek*(s))ds, (19) - T(k + 1)*(t)|| < S1k(t)exp{jH01 (s)dsl. (27)

N T l ( N

J[pNk*, ek*j = I J Jg(t, y, I afk* (t) bf(y), Tk*(t), i= 10 0 ^ j= 1

I j* (t) j(y), ek*(t) 1 bN(y)dydt- bN(x), j = 1 j

k = 1, 2, 3, ... . (20)

Предположим, что для оптимальных управлений справедливы оценки

|P*(t, x) - PNk*(t, x)| < Qn, (t), lim Qn, (t) = 0, (21)

k N ^» k k ^ да

|e*(t) - ek*(t)| < qk(t), lim qk(t) = 0. (22)

k ^ да

Для оценки допускаемой погрешности по состоянию T*(t) докажем, что справедлива

Теорема 2. Пусть выполняются условие (22) и

h(t, T(t), e(t)) e Bnd(M0) n L/p{H01(t)|T; H02(t)|T},

T

где 0 < M0 = const, 0 < J |#0i(t)|dt < да, i = 1, 2.

0

Тогда справедливо соотношение

lim ||T*(t) - T(k + 1)*(t)|| = 0. (23)

k ^ да

Доказательство. Рассмотрим итерационный процесс (19). В силу условия теоремы, имеем

И0 - т(к + 1)4(t)|| < ||T*(t) - + 2)*(t)|| +

+ ||T(k + 2)*(t) - + 1)4(t)|| < j Ho1(s)||T*(s) - + 1)*(s)||ds +

0

+ j Ho2(s)||9*(s) - e(k + 1)*(s)||ds + j Ho1(s)|| T(k + »*(s) -

00 t

- Tk*(s)||ds + jHo2(s)||e(k + »*(s) - ek*(s)||ds < 0 t

< S1k(t) + j Ho1(s)||T*(s) - T(k + 1)4(s)||ds, (24) 0

В силу выражения (26) отсюда следует соотношение (23). Теорема доказана. ♦

Теперь оценим допускаемую погрешность по состоянию и*(?, х).

Теорема 3. Пусть выполняются условия (21), (22) и:

1) /(?, х, и, т) е ВпА(А) п Р/р{#п(?, х)м; Н12(?)||т},

где

0 < J||Hn(s,x)||i2(öi)ds< да, 0 < J||H12(s)||С(0()ds < да;

2) M?)llb2(T) < да Mf + 1 Mf + 11(k + 1)/2 < 1. Тогда справедливо соотношение

lim |и*(?, x) - иж*(?, x)| = 0. ♦ (28)

N ^ да k ^ да

Доказательство см. в Приложении.

5. СХОДИМОСТЬ ФУНКЦИОНАЛА КАЧЕСТВА

Теперь покажем сходимость функционала (20). Здесь справедлива

Теорема 4. Пусть выполняются условия теорем 2 и 3. Если

g(t, x, и, т, P, 8) е е Z/p{H2i(t)|H; H22(t)|x; H23(t)|P; H24(t)|0},

T

где 0 < {H2i.(i)d? < да, i = 1, 4 , то справедливо соот-

0

ношение

lim |J[p*, 8*] - J[pNk*, 8k*]| = 0. (29)

N ^ да k ^ да

Доказательство. В силу условий теоремы из выражений (6) и (20) имеем

т

\J[P*] - J[Pж*]| < \ [Я21(/)|и*(t, x) - uNk*(t, x)| +

+ H22(t)|T*(t) - Tk*(t)| + H23(t)|P*(t, X) - PNk*(t, X)| +

0

0

0

+ H24(t)|e*(t) - ek*(t)|]dt < j H21(t)VNk(t)dt +

+ | я22(г)81к(/)ехр {|Н01( я )&} + 0 10 ] ТТ

+ I ^(^ж^ + \ ЩАМ^Л. (30)

00

С учетом соотношений (21)—(23) и (28) переход к пределу в выражении (30) при N ^ да, к ^ да дает соотношение (29). Теорема доказана. ♦

ЗАКЛЮЧЕНИЕ

Аналитическое решение нелинейных задач оптимального управления очень сложно. На практике широко применяются различные приближенные методы построения программного и синтезирующего оптимального управления. В данной работе решение смешанной задачи (1)—(3) ищется в виде ряда Фурье по собственным функциям

Ь,(х) = Я бш!,.х, X. = Щ , I = 1, 2, ... . Предполага-

1 А/ I II I

ется, что управляющая функция Р(?, х) разлагается в ряд Фурье по собственным функциям Ьг(х). В результате смешанная задача (1)—(3) сводится к изучению счетной системы нелинейных интегральных уравнений (8) при фиксированных значениях управляющих функций. Вместо счетной системы (8) исследуется однозначная разрешимость конечной системы нелинейных интегральных уравнений (11), для чего применяется метод последовательных приближений. Для конечной системы (11) рассматривается итерационный процесс Пикара (12).

В предположении выполнения оценок (21) и (22) доказывается сходимость решения итерационного процесса (12) к решению счетной системы (8) при N ^ да, к ^ да для известных, в том числе оптимальных, управляющих функций Р *(?, х) и 9*(?). Далее доказывается сходимость функционала качества (6), для чего используются последовательности функций (17)—(19) и последовательность функционала (20).

ПРИЛОЖЕНИЕ

Доказательство теоремы 3. Сначала докажем,

что

lim ||aN*(t) - +*(t)||BN_ = 0. (31)

k ^ю B2 ( T)

Аналогично оценкам (14) и (15) из выражения (18) получаем оценку

l|aN*( t) - aN(k +1) *( t)|| ^ <

< ||aN(k + 2)*(t) - aN(k +1)*(t)||Bf(^ + + ||aN*(t) - aN(k+2)*(t)||Bf(^ <

t

< Mi Ж2 /1/2 J[| Hii( s, x„ |£2( щ

+ l|Hl2(s)| C(Dt) l|T(k + 1}*(s) - Tk*(s)| C(D,)] ds + t

+ Mi M22 /1/2 {Г| |Hn (s,

N(k + 1)*, , Nk*, J ,

а^ (s) - а (s)|| BN(t) +

N*, , N(k + 1)>

'11 (X)||£2(D;) Иа (s) - а

(s )l<(t) +

+ llH12(s)| C(Dt) h*(s) - T(k + 1)*

II C( л,)] ds < s0fc(t) +

+ Mi M22 /1/2 {Г| |HU (s,

N*, , N(k + 1)

'11 (s X)||£2(D;) Иа (s) - а

'(*)! <(t) +

+

llH12(s)|| C(Dt) Иs) - T

(k + 1) *

|C( Dt)

]ds, (32)

где

k +1

k + 1 ,.,-2k + 1 , 2

s0k(t) = Mj M2 /

jllH11(s, XЯ£2(D,)ds

k!

+ Mn

jllH12(sC(D)ds

0

_1_ k-!

Теперь, подставляя неравенства (27) в оценку (32), имеем

II aN* (t) - aN(k +1)*( t )|| Bf( ^ < Sok(t) +

t

+ MiM22/1/2 j||H12(s)||c(d)Sik(t)expJ jHo1(s)dsLds +

0

+ MiM22 /1/2 j||HU(s, X)||L2(D1) S

iV*/ ч N

II а (s) - а

ii<(ods.

(33)

Применяя к оценке (33) неравенства Гронуолла — Беллмана, получаем

"N*o "N(k + 1)*(t)llBf(ч <

II а (t) - а

v« +

+ Mi M22 /2 j||H12(s)|| C(Dt) Sik(t)exp j jHo1(s)ds (■ ds

x exp j Mi M22 /1/2 J ||H(s, x)|| £2(d;) ds L. (34)

Так как в силу теоремы lim S0k(t) = 0, то в силу со-

k ^ ю

отношения (26), из выражения (34) получаем соотношение (31).

o

0

0

k

1

А

+

0

0

Теперь оценим величину

VNk = \u*(t, x) - мж*(г, x)\.

(35)

м*

С этой целью оценим разности: и*(г, х) — и (г, х) и ид*(г, х) — иж*(г, х). Сначала рассмотрим следующее соотношение

д».

\u*(t, X) - и (г, x)\ < nN,

где

nn

= J

N

T l

+ и

0 0

ф(г, y)

ф(у) - E ФГbf(y)J [ф(г, y)]t=,

P*(t, y) - E |pn*(t, z)bN(z)dz

i = 10

dy +

bN(У)

dydt +

т l

+ и

00

ф(г, y)

/I t, y, I «* (г) by(y)| -

N l ( N

Ei/It, z, E

i = 10 V * = 1

- E/|t,z, E «N*(t)bN(y)|bN(z)dz

bN (У)

dydt. (36)

n* N

Если a (t) e B2 (T) является решением КСНИУ

(18), то покажем, что

lim nN = 0.

N -

(37)

Действительно, так как aN*(t) e BN (T), то из ра-

венства

N

lim и *(t, x) = lim E а,- * (t)- b; (x) = u*(t, x)

N

ч = 1

в силу условий теоремы следует, что

lim /(t, x, uN*(t, x)) = /(t, x, u*(t, x)) (38)

N — ю

в смысле метрики L2(D).

Тогда первый интеграл в выражении (36) стремится к нулю при N ^ да. Сходимость второй разности в выражении (36) при N ^ да, в силу оценки (21), следует из соотношения

N

lim PN*(t, x) = lim E pN* (t) bN(x) = P*(t, x).

N — ю N — ю ^ _ 1

Сходимость последней разности в выражении (36) при N ^ да следует из соотношения (38). Следовательно, имеет место формула (37). Отсюда следует, что

lim \и*(t, x) - uN*(t, x)\ = 0.

N— ю

(39)

В силу того, что

N

\uN*(t, x) - uNk*(t, x)\ < E \aN* (t) - «N (t)\ • \ bN (x)\ < i = 1

< ИЛ aN* (t) - aNk*( t )||

(T)'

из соотношения (31) получаем

lim \uN*(t, x) - uNk*(t, x)\ = 0.

k — ю

(40)

В силу условий теоремы, соотношений (39) и (40), из выражений (7) и (17) следует, что

lim VNk = lim \u*(t, x) - uNk*(t, x)\ < lim \u*(t, x) -

N— ю nk N— ю N— ю

k — ю k — ю

n* n*

- uN*(t, x)\ + lim \uN*(t, x) - uNk*(t, x)\ = 0,

k — ю

т. е. получаем соотношение (28). Теорема доказана. ♦

ЛИТЕРАТУРА

1. Александров А.Г. Оптимальные и адаптивные системы. — М.: Высшая школа, 1989. — 263 с.

2. Андреев Ю.Н. Управление конечномерными линейными объектами. — М.: Наука, 1976. — 424 с.

3. Вязгин В.А., Федоров В.В. Математические методы автоматизированного проектирования. — М.: Высшая школа, 1989. — 184 с.

4. Кротов В.Ф., Гурман В.И. Методы и задачи оптимального управления. — М.: Наука, 1973. — 448 с.

5. Куропаткин П. В. Оптимальные и адаптивные управления. — М.: Наука, 1980. — 288 с.

6. Бутковский А.Г. Теория оптимального управления системами с распределенными параметрами. — М.: Наука, 1965. — 474 с.

7. Евтушенко Ю.Г. Методы решения экстремальных задач и их применение в системах оптимизации. — М.: Наука, 1982. — 432 с.

8. Егоров А.И. Оптимальное управление тепловыми и диффузионными процессами. — М.: Наука, 1978. — 464 с.

9. Лионс Ж.Л. Оптимальное управление системами, описываемыми уравнениями с частными производными. — М.: Мир, 1972. — 412 с.

10. Лурье К.А. Оптимальное управление в задачах математической физики. — М.: Наука, 1975. — 480 с.

11. Рапопорт Э.Я. Оптимальное управление системами с распределенными параметрами. — М.: Высшая школа, 2009. — 680 с.

12. Срочко В.А. Итерационные методы решения задач оптимального управления. — М.: Физматлит, 2000. — 160 с.

13. Тятюшкин А.И. Численные методы и программные средства оптимизации управляемых систем. — Новосибирск: Наука, 1992. — 193 с.

14. Федоренко Р.П. Приближенное решение задач оптимального управления. — М.: Наука, 1978. — 488 с.

15. Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического типа. — М.: Наука, 1967. — 736 с.

16. Мазья В.Г. Пространства С. Л. Соболева. — Л.: Изд-во ЛГУ, 1985. — 416 с.

17. Соболев С.Л. Некоторые применения функционального анализа в математической физике. — М.: Наука, 1988. — 336 с.

18. Юлдашев Т.К. Смешанная задача для нелинейного интег-ро-дифференциального уравнения с параболическим оператором высокой степени // Журн. вычисл. матем. и ма-тем. физики. — 2012. — Т. 52, № 1. — С. 112—123.

19. Юлдашев Т.К. Развитие теории нелинейных дифференциальных уравнений с максимумами: дис. ... канд. физ.-мат. наук. — Бишкек: Институт математики НАН Кыргызской Республики, 1993. — 121 с.

Статья представлена к публикации членом редколлегии

Е.Я. Рабиновичем.

Юлдашев Турсун Камалдинович — канд. физ.-мат. наук,

Сибирский государственный аэрокосмический университет

им. академика М.Ф. Решетнева, г. Красноярск,

® (391) 291-91-19, И tursunbay@rambler.ru.

0

—ю