Об одной модели оптимального управления температурным режимом теплицы Текст научной статьи по специальности «Математика»

УДК 517.977.56

И.В. Асташова, Д.А. Лашин, А.В. Филиновский1

ОБ ОДНОЙ МОДЕЛИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ ТЕМПЕРАТУРНЫМ РЕЖИМОМ ТЕПЛИЦЫ

При выращивании растений в промышленных теплицах требуется поддерживать температуру в точке роста растений, находящейся на фиксированной высоте, в соответствии с заданным суточным графиком температур, допуская малые отклонения. При этом можно увеличивать температуру, увеличивая подогрев пола теплицы и уменьшать температуру, открывая форточки на ее потолке. Далее поставим задачу поддержания на некоторой заданной высоте c температуры z(t) в течение промежутка времени 0 < t < T. Для решения задачи предлагается и анализируется математическая модель, использующая уравнение теплопроводности. Физический смысл данной задачи заключается в том, что на одном конце бесконечно тонкого стержня длины l (высота теплицы) в течение времени T поддерживают температуру ф({) (управляющая функция), а на другом конце задан тепловой поток ф(t). Требуется найти такую управляющую функцию ф0(Ь), при которой температура в определенной точке c была бы максимально близка к заданной температуре z(t). Оценка качества управления осуществляется с помощью квадратичного интегрального функционала.

Ключевые слова: оптимальное управление, температурный режим, теплица, уравнение теплопроводности, квадратичный интегральный функционал.

Введение

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

!© Асташова И.В., Лашин Д.А., Филиновский А.В., 2016

Асташова Ирина Викторовна (ast@diffiety.ac.ru), кафедра дифференциальных уравнений, механико-математический факультет, Московский государственный университет имени М.В.Ломоносова, 119991, Российская Федерация, г. Москва, Ленинские горы, 1; кафедра высшей математики, факультет МЭСИ, Российский экономический университет им.Г.В.Плеханова, 117997, Россия, Москва, Стремянный переулок, 36.

Лашин Дмитрий Александрович (dalasin@gmail.com), ООО НПФ ФИТО, 142784, Российская Федерация, г. Москва, Московский, 35-12.

Филиновский Алексей Владиславович (flnv@yandex.ru), кафедра высшей математики, Московский государственный технический университет им. Н.Э.Баумана, 105005, Российская Федерация, г. Москва, 2-я Бауманская ул., 5; кафедра дифференциальных уравнений, механико-математический факультет, Московский государственный университет имени М.В.Ломоносова, 119991, Российская Федерация, г. Москва, Ленинские горы, 1.

получены в работах [1] - [3]. Некоторые методы доказательства основных результатов содержатся в [4] и [5]. Подобные экстремальные задачи с интегральными функционалами рассматривались различными авторами ([6] - [9]). Обзор ранее полученных результатов можно найти в работе [10], библиография последних работ содержится в [11]. Задача минимизации функционала с финальным наблюдением и задача минимизации времени управления рассматривались в работах [7] - [11]. См. также [12], [13].

1. Основные обозначения и постановка задачи

1. Рассмотрим смешанную задачу для уравнения теплопроводности

щ = (а(х)их)х, 0 < х < I, г > 0, (1.1) с достаточно гладким коэффициентом а(х), удовлетворяющим условию

0 < ао ^ а(х), х € [0,1], включающую краевые условия

и{0и1=Ф(% *> 0> (1.2)

и начальные условия

и(х, 0) =0, 0 <х <1, (1.3)

где Ф € (0,Т), ф € W1(0,T) для любого Т > 0.

В этой работе мы будем предполагать, что ф(Ф) — заданная функция, а ф(Ф) — искомая управляющая функция.

Пусть Qт = (0,1) х (0,Т), а У^^т) (см. [14], ч. 1, §1(2)) — банахово пространство, состоящее из элементов W1'0(Qт) с конечной нормой

(I \ 1/2 1/2

/ и(х,г)2йх\ + ( / их(х,г)2 ¿хЛ\ , (1.4)

¿0 у/ )

имеющих непрерывно меняющиеся по г € [0,Т] следы из Ь2(0,1) на сечениях (0,1).

Обозначим через W2l(Qт) пространство, состоящее из таких элементов п(х,г) € W21 ^т), что п(',Т) = 0, п(0, ■) = 0.

Будем рассматривать обобщенное (слабое) решение задачи (1.1)-(1.3) из энергетического класса, то есть такую функцию и € У2'0^т), что и(0,г) = ф(г), удовлетворяющую интегральному тождеству

/ (а(х)их(х,г)пх(х,г) — и(х,г)пь(х,г)) ¿хА = а(1) 'ф(1)п(1,г)А (1.5) JQT ./о

для любой функции п(х,г) € W2(QT).

Имеет место следующая теорема, которую приведем без доказательства. Теорема 2.1. Обобщенное решение задачи (1.1)-(1.3) из класса У-2'0^т) существует, единственно и удовлетворяет неравенству

\икт < С(\ф\ш1(о,т) + \\ф\\шЦо,т)), (1.6)

где постоянная С не зависит от ф и ф.

Отметим, что доказательство существования решения проводится методом Га-леркина, а доказательство единственности и оценка нормы решения получены с использованием энергетического неравенства. Обозначим это единственное решение иф.

2. Далее поставим задачу поддержания на некоторой заданной высоте с температуры г(Ь) в течение всего промежутка времени 0 ^ £ ^ Т.

Рассмотрим задачу (1.1)-(1.3). Пусть г € Ь2(0,Т). Обозначим через Фм, где М > 0, множество функций

Фм = [ф € W1(0,T), \\ф\\Ш1(о,т) < М}.

Для произвольного с € (0,1] определим функционал

3[ф]= [ (пф(с,1) - г(1))23,1 (1.7)

■)о

Рассмотрим задачу минимизации данного функционала. Обозначим

Ы 3[ф] = т. (1.8)

ФеФы

Физический смысл данной задачи заключается в том, что на одном конце бесконечно тонкого стержня длины I в течение времени Т поддерживают температуру ф(1) (управляющая функция), а на другом конце задан тепловой поток ф(1). Требуется найти такую управляющую функцию фо(£), при которой температура в фиксированной точке х = с была бы максимально близка к заданной температуре Оценка качества управления осуществляется с помощью функционала 3 [ф].

2. Существование и единственность оптимальной управляющей функции

Теорема 3.1. Существует единственная функция фо(£) € Фм, для которой т = 3 [фо].

Доказательство теоремы.

В доказательстве существования и единственности минимизирующей функции функционала 3 используется следующая лемма.

Лемма 3.1 Пусть А — выпуклое замкнутое подмножество гильбертова пространства Н. Тогда для любого х € Н существует единственный элемент у € А, для которого

\\х - У\\ = Ш \\х - г\\. (2.1)

геА

Обозначим

Вм = [иф(с, ■), ф € Фм} С Ь2(0,Т)

и докажем, что Вм — выпуклое замкнутое подмножество в (0,Т). Для любых у^ = и,ф. (с, ■) € Вм, ] = 1, 2, имеем \\ф^ ^¿(от) ^ М, ] = 1, 2. Следовательно, для всех а € (0,1) выполняется неравенство

\\аф1 + (1 - а)ф2\ш1(0,т) < ^Лшцот) + (1 - а)\\ф2\\ш^{0,т) < M,

так что ау1 + (1 - а)у2 € Вм и множество Вм выпукло. Теперь докажем, что Вм замкнуто в Ь2(0,Т).

Пусть функции yk = Пфк (с, ■) £ BM образуют фундаментальную последовательность в Ьэ(0,Т), стремящуюся к y £ Ьэ(0,Т). Соответствующие функции [фк} £ ФМ образуют слабо предкомпактное множество в W2;(0,T). Таким образом, для некоторой подпоследовательности фи6 существует слабый предел ф £ W2(0,T) при j ^ ж. Из свойств слабо сходящихся последовательностей в гильбертовом пространстве (см. [15], Гл. 1, §1, T. 1.1) следует

UWwi(0,T) < limsup \\фкз \ \w1(0,t) < M, (2.2)

j^ж

так что ф £ ФМ. Теперь, по теореме Банаха-Сакса (см. [16], Ch. 2, Sec. 3) существует такая подпоследовательность kjn, что

lim \\фп - ф\^}(0,т) =0, (2.3)

п

где фп = пЕ фиj,. Следовательно,

i=i

1 п

\\фп\^(0,Т) < \\фкп\w1(0,T) < M, (2.4)

n

l=l

и, используя (2.2), мы получим

1 п

yn = nYsyk>l £ BM.

n i=i

Таким образом, для соответствующей последовательности решений

1 п

Щп = иФк:)1 1 = 1

из оценки (1.6) мы получим неравенства

\ит — Un\Qт < с\\фт — фп\\ш^(0,т) ^ 0, т,п ^ ж. (2.5)

Функция ип является слабым решением задачи

(Щп)ь = (а(х)(ип)х)х, 0 <х<1, Т >г> 0, йп(0,г) = фп(г), Т >г> 0,

(йп)х(1,г) = Ф(г), Т >г> 0,

йп(х, 0) = 0, 0 < х < I. Это означает, что йп(0,г) = фп(г) и интегральное тождество

/ (а(х)(йп(хЛ))хцх(х,г) — йп(х.г)%(х.г)) ¿хА = а(1) ф(1)п(1,г)А (2.6)

JQT Jо

выполняется для каждой функции п € ^т). Принимая во внимание соотношения (2.3), (2.5) и (2.6), мы видим, что существует предельная функция = = иф € У21,0^т), которая является слабым решением задачи (1.1)—(1.3) с граничной функцией ф и

\и — йп^т < С\\ф — фп\ш1(0,т). (2.7)

Таким образом, по неравенству, следующему из теорем вложения (см. [15], гл. 1, § 6, (6.15)), мы получим

\\у — й'п^ ■)\\ь2(0,т) = \\и(С ■) — у-п^ ■)\\ь2(0,т)

< \и — йп^т < С\\ф — (рп^ш1 (0,т). (2.8)

Из (2.8) следует, что у = и(с, ■) € Вм и Вм является замкнутым подмножеством в 12(0,Т).

Таким образом, по лемме 3.1 существует единственная функция у(Ь) = и(с,Ь), где и € У21,0(ЯТ) — решение задачи (1.1)—(1.3) с такой функцией ф = ф0 € Фм, что

3[ф] = 3ф]. (2.9)

ФеФы

Докажем, что такая функция ф0 € Фм единственна. Если это не так, то рассмотрим две такие функции ф1, ф2 и соответствующую пару решений и1, и2. Функция и = и1 - и2 является решением задачи

щ = (а(х)их)х, 0 <х<1, 0 <КТ, (2.10)

й(0,г) = ф1(г) - ф2(ь) = Ф(г), 0 <г<Т, (2.11)

и(с,г) = 0, 0 <г<Т, (2.12)

их(1,г) = 0, 0 <г<Т, (2.13)

и(х, 0) = 0 0 < х < I. (2.14)

Теперь мы докажем, что функция и в прямоугольнике QT) = (с,1) х (0,Т) совпадает с решением задачи

щ = (а(х)их)х, с < х < I, 0 <г <Т, (2.15)

и(с,г) = 0, 0 <г<Т, (2.16)

их(1,г) = 0, 0 <г<Т, (2.17)

и(х, 0)=0, с < х < I. (2.18)

Чтобы доказать это, достаточно положить в интегральном тождестве (1.5) функцию п равной 0 на [0, с] х [0, Т]. Но решение задачи (2.15)—(2.18) равно нулю на [с,1] х [0, Т], следовательно, мы имеем

и(х,г) = 0, с<х <1, 0 <г<Т. (2.19)

Теперь мы докажем, что

и(х,г) = 0, 0 < х < I, 0 <г<Т. (2.20)

Заметим, что по теореме 2 из [17], §11, слабое решение и является классическим решением уравнения (2.10) в Qт. Теперь мы используем теорему 5 из [18], §3, где доказано следующее утверждение. Рассмотрим такую функцию и € С2,1 (О), О С Д2, что

щ = (а(х)их)х, (х,1) € О С Д2.

Пусть О0 — связная компонента множества О П [£ = ¿0}, а О — связное открытое подмножество множества О0. Тогда если = 0, то и\а0 =0.

Применяя эту теорему к решению и задачи (2.10)—(2.14) для любого ¿0 € (0,Т) и О0 = (0,1) х {¿0}, (О = (с,1) х {¿0}, мы получим, что из (2.19) следует (2.20). Таким образом, имеем и(х,Ь) = 0 для всех х € (0,1) и £ € (0,Т). Это означает, что ф(Ь) = и(0,Ь) = 0. Теорема доказана.

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

Обозначим через Фм класс управляющих функций

Ф0м = [ф € Ш1(0,Т) : ЦфЦ^т) < М, ф(0) = 0}.

Теорема 3.2. Существует единственная функция ф0(г) € Фм, для которой

Ш 3[ф] = 3ф].

Феф0ы

Обозначим через Фм класс управляющих функций

Ф0м = [ф € W1(0,T) : \\ф\^№,т) < М,ф1 < ф(г) < ф2}

с некоторыми постоянными ф1 и ф2.

Теорема 3.3. Существует единственная функция ф0(г) € Фм, для которой

Ш 3[ф] = 3[ф0].

феф Ы

3. О точной управляемости

Наряду с вопросом о существовании и единственности решения экстремальной задачи, важным вопросом является точная управляемость. Под точной управляемостью будем понимать возможность получения в точке х = с следа и(с,г) почти всюду совпадающего на [0, Т] с заданной функцией г (г). Соответственно, точным управлением будем называть такую функцию ф0 € Фм, при которой функционал 3[ф0] обращается в нуль:

С т

3[ф0]= (иф0 (с, г) - г(г))2 ¿г = 0.

0

Следующая теорема показывает, что множество функций г € Ь2(0,Т), допускающих точную управляемость, достаточно "мало".

Теорема 4.1. Множество всех функций г € Ь2(0,Т), для которых возможна точная управляемость, т.е. 3[ф] =0 для некоторого ф € Фм, является множеством первой категории в Ь2 (0,г).

Доказательство. Рассмотрим уравнение (1.1) для функции и1 € У2'0^т) с краевыми условиями

и1 (0,г) = ф1(г), 0 <г<Т, (и\)х(1,г) = Ф(г), 0 <г<Т, и такое же уравнение для функции и2(х,г) € У2'0^т) с краевыми условиями

и2(0,г) = ф2(г), 0 <г<Т, (и2)х(1,г) = Ф(г), 0 <г<Т.

Функция и = и1 - и2 является решением уравнения (1.1) с краевыми условиями

й(0,г) = ф(г) = ф1(г) - ф2(г), 0 <г<Т, (3.1)

йх(1,г) = 0, 0 <г<Т, (3.2) и начальным условием

й(х, 0) = 0, 0 < х < I. (3.3) Теперь рассмотрим в прямоугольнике QT2l) = (0, 21) х (0, Т) задачу

щ = (а(х)их)х, 0 <х < 21, 0 <г<Т, (3.4)

а(0,г) = ф>(г), 0 <г<Т, (3.5)

а(21,г) = ф>(г), 0 <г<Т, (3.6)

а(х, 0)=0, 0 <х< 21. (3.7)

где а(х) = а(х), х € (0,1), а(х) = а(21 — х), х € (1,21). Слабое решение задачи (3.4)-(3.7) — это функция йа(х,€) € У~2'°(Ят^), удовлетворяющая краевым условиям и(0,1) = и(21,Ь) = ф(Ь) и интегральному тождеству

IqT21

(a(x)ux(x, t)nx(x, t) — u(x, t)nt(x, t))dxdt = 0 (3-8)

для любой такой функции п & W^Qt), что ц(^,Т) = 0, п(0, •) = 0, n(2l, •) = 0. Из равенства(3-8) следует

u(x, t) = U(x, t), 0 < x < l, 0 <t <T. (3-9)

По принципу максимума для слабых решений (см. [14], ч. 3, §7, теорема 7-2), решение U(x,t) удовлетворяет неравенствам

ess inf 4>(t) < u(x,t) < ess sup фСЬ). (3.10)

te[0,T] te[0,T]

Из (3.10) получим

esssup \U\ ^ ess sup ^i(t) — ф2^)\, (3.11)

Qt 0<t<T

и, следовательно,

ess sup \U(c,t)\ ^ ess sup \ф\(t) — ф2^)\. (3.12)

0<t<T 0<t<T

Интегрируя неравенство (3.12), мы получим

f u(c,t)dt < t( sup \ф1 (t) — ф2(Ь)\) . (3.13)

0 0<t<T

Пусть функции ф1(Ь) и ф2(t) являются точными управляющими функциями для данных zi(t) и z2(t). Это означает, что

сТ

Аф1]= ! (ифг (c,t) — zi(t))2dt = 0,

0

J[ф2]= f (иф2 (c,t) — Z2(t))2dt = 0.

0

В этом случае неравенство (3.13) влечет неравенство

f (zi(t) — Z2(t))2 dt < Т (ess sup ^i(t) — ф2(t) A (3.14)

0 0<t<T

для произвольных функций zi (t) и Z2(t), допуская точную управляемость.

Пусть Z С L2(0,T) — множество функций z(t), допускающих точное управление. Имеем Z = U^=iZm, где Zm С L2(0,T) — множество функций, точно управляемых с помощью ф & Фм. Рассмотрим произвольную последовательность управляющих функций {фи} С Фм, M = 1, 2,..., и соответствующую последовательность {zk} С Zm. Множество Фм — ограниченное множество W2i(0,T). По теореме вложения для W)L (0, Т) имеем

ess sup \фи1 — фи, \ ^ 0, l,j ^ ж, (3.15)

0<t<T

для некоторой подпоследовательности фи-. Из (3.14) и (3.15) мы получим для последовательности {zkj (t)} С ZM соотношение

f (zkl (t) — zkj (t))2dt < T (ess sup \фк1 (t) — фи- (t)A ^ 0, j,l ^ ж. (3.16)

0 0<t<T

Из (3.16) следует, что Zм — предкомпактное множество в Ь2(0,Т). Таким образом, Zм нигде не плотно в Ь2(0,Т). А так как Z = , мы заключаем, что Z — множество первой категории в Ь2(0,Т). Теорема 4.1 доказана.

Далее будут сформулированы утверждения, показывающие, что точная управляемость может не иметь места не только для функций г € Ь2(0,Т), но и для г € С([0,Т]).

Рассмотрим вопрос о точной управляемости в задаче (1.1)—(1.3) в случае, когда Ф(г) = 0 (отсутствует поток через правый конец).

аь = (а(х)ах)х, 0 < х < I, 0 <г<Т, (3.17)

а(0,г) = ф(г),

ах(1,г) = 0, а(х, 0) = 0,

Теорема 4.2. Для любого М > 0 существует такая функция г € С([0,Т]), что для любой функции ф € Фм для задачи (3.17)-(3.19) справедливо неравенство 3ф > 0.

Рассмотрим теперь более общий случай, когда Ф(г) = 0, т. е. будем рассматривать задачу (1.1)-(1.3). Тогда имеет место следующее утверждение.

Теорема 4.3 Для любых М > 0 и М1 > 0 существует такая функция г € С([0,Т]), что для любой функции ф(г) € Фм и любой такой функции Ф(г) € W2(0,T), что ЦФ(г)\\Ш1(0,т) ^ М1, для задачи (1.1)-(1.3) выполняется неравенство

С т

3[ф]= (иф(с,г) - г(г))2 ¿г> 0.

0

Литература

[1] Лашин Д.А. Стратегия управления микрокламатом в теплицах // Гавриш. Москва, 2005. № 1. C. 33-35.

[2] Лашин Д.А. Об оптимальном управлении температурным режимом // Дифференц. уравнения. 2008. T. 44. Вып. 6. C. 853.

[3] Lashin D. A. On the existence of optimal control of temperature regimes // J. of Math. Sci. 2009. V. 158. № 2. P. 219-227.

[4] Some Problems in the Qualitative Theory of Differential Equations / I.V. Astashova [et al.] // J. of Natural Geometry. Jnan Bhawan. London. 2003. V. 23. № 1-2. P. 1-126.

[5] Качественные свойства решений дифференциальных уравнений и смежные вопросы спектрального анализа / под ред. И.В. Асташовой, 2012. М.: ЮНИТИ-ДАНА. 647 с.

[6] Лионс Ж.-Л. Оптимальное управление системами, описываемыми уравнениями с частными производными. М.: Мир. 1972.

[7] Бутковский А.Г. Оптимальное управление в системах с распределенными параметрам // Автоматика и телемеханика. 1961. T. 22. № 1. C. 17-26.

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

[9] Егоров Ю.В. Некоторые задачи теории оптимального управления // Ж. Выч. мат. и мат. физ. 1963. T. 3. № 5. C. 887-904.

0 <t<T, 0 <t<T,

0 < x < l.

(3.18)

(3.19)

Butkovsky A.G., Egorov A.I., Lurie K.A. Optimal control of distributed systems // SIAM J. Control. 1968. Vol. 6. № 3. P. 437-476.

Фурсиков А.В. Оптимальное управление распределенными системами. Теория и приложения., Новосибирск: Научная книга, 1999.

Гудвин Г.К., Гребе С.Ф., Сальгадо М.Э. Проектирование систем управления. Москва.: БИНОМ. 2004.

Farag M.H., Talaat T.A., Kamal E.M. Existence and uniqueness solution of a class of quasilinear parabolic boundary control problems // Cubo. 2013. V. 15. № 2. P. 111-119.

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

Ладыженская O.A. Краевые задачи математической физики, Москва: Физматлит, 1973.

16] Riesz F., Szokefalvi-Nagy B., Functional Analysis, New-York: Dover, 1990.

17] Ильин А.А., Калашников А.С., Олейник О.А. Линейные уравнения второго порядка параболического типа // Успехи матем. наук. 1962. 17. Вып. 3. C. 3-146.

[18] Ландис E. M., Олейник O.A. Общенная аналитичность и некоторые связанные с ней свойства решений эллиптических и параболических уравнений // Успехи матем. наук. 1974. Т. 29. Вып. 2. P. 190-206.

References

[1] D. A. Lashin. Strategy of management to microclimate in hothouses Gavrish. Moscow, 2005. no. 1, p.33-35. (in Russian)

[2] D. A. Lashin. On the optimal control of a temperature regime. Differ. Equ., 2008, V. 44, no. 6, P. 853.

[3] D. A. Lashin. On the existence of optimal control of temperature regimes. J. of Math. Sci., 2009, V. 158, no. 2, P. 219-227.

[4] I. V. Astashova, A. V. Filinovskiy, V. A. Kondratiev, L. A. Muravei, Some Problems in the Qualitative Theory of Differential Equations // J. of Natural Geometry. Jnan Bhawan. London. 2003, V. 23, no. 1-2, P. 1-126.

[5] Qualitative Properties of Solutions to Differential Equations and Related Topics of Spectral Analysis: scientific edition, edited by I. V. Astashova, Moscow: UNITY-DANA, 2012, 647 p. (in Russian)

[6] J. L. Lions, Optimal control of systems governed by partial differential equations, Berlin: Springer, 1971.

[7] A. G. Butkovsky. Optimal Control in the Systems with Distributed Parameters. Avtomatika i Telemechanika, 1961, 22, no. 1, P. 17-26.

[8] A. I. Egorov, Optimal Control by Heat and Diffusion Processes, Moscow: Nauka, 1978. (in Russian)

[9] Yu. V. Egorov, Some Problems of Theory of Optimal Control Zhurnal. Vych. Mat. i Mat. Fiziki. 1963, V. 3, no. 5, P. 887 - 904. (in Russian)

[10] A. G. Butkovsky, A. I. Egorov, K. A. Lurie. Optimal control of distributed systems. SIAM J. Control. 1968, Vol. 6, no. 3, P. 437 - 476.

[11] A. V. Fursikov, Optimal Control of Distributed Systems. Theory and applications, Novosibirsk: Nauchnaya Kniga, 1999. (in Russian)

[12] G. C. Goodwin, S. F. Graebe, M. E. Salgado, Control system design, London -New-York: Pearson, 2000.

[13] M. H. Farag, T. A. Talaat, E. M. Kamal. Existence and uniqueness solution of a class of quasilinear parabolic boundary control problems. Cubo, 2013, V. 15, no. 2, P. 111-119.

[14] O. A. Ladyzhenskaya, V. A. Solonnikov, N. N. Ural'seva, Linear and quasi-linear equations of parabolic type, Translations of Mathematical Monographs, 23, Providence, RI: American Mathematical Society, 1968.

[15] O. A. Ladyzhenskaya, Boundary value problems of mathematical physics, Moscow: Fizmatlit, 1973.

[16] F. Riesz, B. Szokefalvi-Nagy, Functional Analysis, New-York: Dover, 1990.

[17] A. M. Ilin, A. S. Kalashnikov, O. A. Oleinik. Linear equations of second order of parabolic type. Russ. Math. Surways, 1962, V. 17, issue 3, P. 3-146.

[18] E. M. Landis, O. A. Oleinik. Generalized analiticity and some connected properties of solutions of elliptic and parabolic equations. Russ. Math. Surways, 1974, V. 29, issue 2, P. 190-206.

I.V. Astashova, D.A. Lashin, A.V. Filinovskiy2

ON A MODEL OF OPTIMAL TEMPERATURE CONTROL

IN HOTHOUSES

While growing plants in industrial hothouses it needs to keep the temperature according to round-the-clock graph at the point of growth of plant located at the fixed height. Only small deviations are admitted. To obtain this it is possible to increase the temperature by heating the floor and to decrease the temperature by opening the ventilator windows at the ceiling. We propose and analyse the model based on the heat equation. Physical sense of this problem is that at one end of the infinitely thin rod of length l (the height of the hothouse) we keep during the time T the temperature (control function), while at the other

end we have the given heat flow ^(t). It requires to find the control function (t) such that the temperature at the fixed point c be maximally closed to the given temperature z(t). For the estimation of the control quality we use a quadratic integral functional.

Key words: optimal control, temperature control, hothouse, heat equation, quadratic integral functional.

Статья поступила в редакцию 18/VT/2016. The article received 20/VI/2016.

2Astashova Irina Viktorovna (ast@diffiety.ac.ru), the Dept. of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow, 119991, Russian Federation; Dept.of MESI, Plekhanov Russian University of Economics, 117997, Stremyanny lane, 36, Moscow, Russian Federation.

Lashin Dmitriy Alexandrovich (dalasin@gmail.com), FITO research and production company, 142784, Russian Federation, Moscow, Moscovskiy, 35-12.

Filinovskiy Alexey Vladislavovich (flnv@yandex.ru), Bauman Moscow State Technical University, 105005, Baumanskaya 2nd st., 5; Lomonosov Moscow State University, Faculty of Mechanics and Mathematics, 119991, GSP-1, Leninskiye Gory, 1, Moscow, Russian Federation.