CC BY
33ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2020. Т. 28. № 4
Информатика, вычислительная техника и управление
УДК 681.5.015
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ В ОБРАТНЫХ ЗАДАЧАХ ТЕПЛОПРОВОДНОСТИ В УСЛОВИЯХ ИНТЕРВАЛЬНОЙ НЕОПРЕДЕЛЕННОСТИ НА ОСНОВЕ НЕЙРОННЫХ СЕТЕЙ*
А.Н. Дилигенская, А.В. Самокиш
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Аннотация. Работа посвящена развитию методов параметрической идентификации, применяемых для решения обратных задач технологической теплофизики в условиях интервальной неопределенности параметров. При решении таких задач в условиях действия возмущающих факторов использование искусственных нейронных сетей, воспроизводящих структуру оператора прямой задачи при небольшой размерности искомого вектора параметров, позволяет получить удовлетворительные результаты. На примере граничной обратной задачи теплопроводности предложен вариант построения нейронной сети для решения задачи с известной структурой математического оператора, характеризующегося небольшим числом параметров. Метод позволяет использовать априорную информацию о допустимом диапазоне принадлежности идентифицируемых характеристик или их параметров. Входными данными для нейросетевой модели являются все возможные реализации температурных состояний, полученные как реакции на входное воздействие, в качестве которого выступает идентифицируемая характеристика, удовлетворяющая всем допустимым вариантам сочетания параметров.
Решение задачи представляет итерационную процедуру, на каждом шаге которой с помощью нейронной сети уточняется значение одного из параметров, задающих параметрическое представление искомой величины. Сочетание нейросетевых технологий с качественной информацией о закономерностях протекания процессов нестационарной теплопроводности позволяет повысить точность решения задачи параметрической идентификации в условиях действия возмущений. Продемонстрированы варианты решения граничной обратной задачи теплопроводности при одно- и двухпараметрическом представлении искомой характеристики в условиях действия возмущений. Исследована точность решения задачи в зависимости от интенсивности возмущающего воздействия.
*Работа выполнена при финансовой поддержке РФФИ в рамках научных проектов № 18-08-00048 и № 18-08-00565.
Дилигенская Анна Николаевна, д.т.н., профессор кафедры «Автоматика и управление в технических системах».
Самокиш Александр Валерьевич, аспирант кафедры «Информационные и развивающие системы и технологии». 6
Ключевые слова: обратная задача теплопроводности, параметрическая
идентификация, интервальная неопределенность, нейросетевая модель.
Введение
Теория идентификации широкого ряда технологических процессов или технических объектов и систем, их свойств и отношений находит свое применение в различных областях науки и техники и является базовым компонентом для эффективного решения задач управления производствами или для проектирования промышленного оборудования. В каждой производственной сфере разрабатываются наиболее пригодные для использования подходы и методы структурной, параметрической или непараметрической идентификации, определенные спецификой предметной области исследований.
В области технологической теплофизики задачи идентификации процессов нестационарной теплопроводности сводятся к определению по измеряемым температурным полям неизвестных характеристик, среди которых наиболее типичными являются задачи восстановления внешних воздействий (тепловых потоков на поверхности, коэффициентов теплообмена между взаимодействующими средами), функций мощности внутренних тепловых источников, теплофизических характеристик материалов или начального состояния системы [1]. Такие задачи, предусматривающие нахождение одной или нескольких характеристик процесса по результатам измерения температурных полей, относятся к обратным задачам теплопроводности (ОЗТ), которые в зависимости от идентифицируемой характеристики подразделяются на граничные, внутренние, коэффициентные и ретроспективные ОЗТ.
Разработка методов решения обратных задач математической физики, в том числе обратных задач теплопроводности, является актуальным, быстро развивающимся направлением, и на сегодняшний день существует большое число методов [2, 3], основанных на различных приемах и подходах теории решения некорректных и обратных задач математической физики, идентификации, вычислительной математики, теории оптимального управления, искусственного интеллекта.
В большинстве подходов решение задачи структурной идентификации, как правило, сводится к формализации и структуризации исследуемых явлений в рамках совокупности конечного числа математических моделей, разработке алгоритмов генерации структурных зависимостей (альтернатив) и их перебору с целью выбора альтернативы, наиболее удовлетворяющей поставленному многофакторному критерию. Наличие априорной информации об операторе искомой характеристики и о математической модели объекта позволяет сразу перейти к этапу параметрической идентификации.
Если математическая модель, содержащая искомую характеристику, в остальном известна с высокой степенью достоверности и уровень возмущающих воздействий не является высоким, то значительными преимуществами обладают методы параметрической идентификации, основанные на использовании аналитических условий оптимальности [4, 5]. Эти методы показали эффективность их использования и высокую достоверность результатов в сочетании с относительно невысокой трудоемкостью реализации.
В ситуации, когда анализ объекта исследования может быть отнесен к слабоструктурированным проблемам вследствие сочетания неопределенных факторов различной природы, возникает потребность в разработке альтернативных
подходов, позволяющих находить регулярные решения в условиях неполноты и неточности исходной информации [6-9].
В данной статье предлагается подход к решению задачи параметрической идентификации на основе применения искусственных нейронных сетей (ИНС) [10-12], использующих информацию о структуре оператора прямой задачи при небольшой размерности искомого вектора параметров, в условиях неполноты и неточности исходной информации, обусловленных действием возмущающих факторов.
Постановка обратной задачи технологической физики
в условиях неполноты информации
Рассматривается базовая математическая модель относительно температурного поля Т(х,г) в процессе нестационарной теплопроводности одномерного
объекта канонической формы, заданная одномерным линейным уравнением Фурье в относительных координатах
дТ (х,г ) д2Т ( х,г )
—= —Ц!^, о < х < 1,0 < г < г , (1)
дг дх2
дополненная краевыми условиями
Т(х,г) = 0, х ф1]; дТ(0) = 0; = ч(г). (2)
На основе представленного математического описания формулируется модельная граничная ОЗТ, предусматривающая нахождение неизвестной функции плотности теплового потока д({) , действующего на границе х = 1 заданной пространственной области, при известных прочих параметрах и характеристиках по полученной экспериментальной информации Т = Т (х ,г) о температуре в некоторой точке х* е[0,1]. Результаты эксперимента неизбежно содержат возмущения, о которых известно, что они удовлетворяют нормальному закону распределения плотности вероятностей.
Выбор структуры математической модели основывается на физических особенностях процессов технологической теплофизики. Вследствие инерционности тепловых объектов идентифицируемые характеристики и результирующая температура, как правило, не содержат резких локальных скачков, что позволяет в зависимости от типа исследуемой ОЗТ ограничить выбор семейства структур искомых величин классами достаточно гладких функций.
Для основополагающих граничных ОЗТ изменение плотности граничного теплового потока может быть описано согласно [1] с помощью следующих типовых функций: постоянных, монотонных, линейных (возрастающих или убывающих), параболических, кубических, экспоненциальных, а также гармонических, или их комбинаций. Каждый альтернативный вариант структуры математической модели характеризуется определенным вектором параметров
=(/\ ,7 ,'"7п), значения которых могут быть найдены в результате параметрической оптимизации с использованием сформулированного определенным образом функционала качества. Для исследования базовых основ и возможностей использования нейросетевых технологий для решения задачи параметрической
идентификации процессов нестационарной теплопроводности в работе в качестве математических операторов искомой функции - плотности теплового потока - рассматриваются простые базовые функциональные зависимости от времени, характеризующиеся малым числом параметров: постоянные и экспоненциальные функции.
Построение нейросетевой модели процесса нестационарной
теплопроводности
Для решения обратной задачи теплопроводности в условиях действия возмущающих факторов использованы искусственные нейронные сети [13-15], обеспечивающие непосредственное восстановление искомых параметров на основе входной информации, что соответствует переходу от процедуры минимизации невязки к аппроксимации идентифицируемой характеристики путем обучения построенной ИНС. В основу такого подхода положена способность нейронных сетей аппроксимировать искомые зависимости с требуемой точностью.
Предполагается, что решение ОЗТ осуществляется при наличии некоторого объема априорной информации о значениях параметров восстанавливаемой характеристики, что, как правило, соответствует минимальным сведениям об условиях реализации конкретного процесса. Эта информация используется для формирования обучающей выборки нейросетевой модели. Априорная информация может принимать форму интервальной неопределенности, когда доступные сведения заключаются в информации о границах диапазона возможного изменения идентифицируемых параметров у е /п): уШп < у < у;тах. Заданные диапазоны
' Утах) изменения значений для каждого из параметров у, г =1,п с помощью равномерного шага дискретизации Ду разбиваются на одномерные массивы
(уо,уо,~уо,); уг) =у уг) =У1Шх;уг+)1 -уг) =ДУг Л = 1>к-1, кот°рые в дальнейшем используются в качестве входных данных ИНС. Выбор шага дискретизации может выбираться исследователем и в общем случае зависит от достоверности априорной информации. Так, если интервал возможного изменения параметров слишком широк, то возможно поэтапное решение задачи, предусматривающее постепенное сужение допустимого диапазона и, соответственно, уменьшение шага дискретизации. В таких условиях при использовании линейной математической модели, заданной в относительных координатах, приведение исходных данных к относительным единицам измерения позволит достичь унификации процедуры обучения ИНС.
Обучающая выборка, содержащая результирующую температуру Т (х',¡), полученную в некоторой точке х* е[0,1] на интервале наблюдения и соответствующую искомой функции д(/), формируется для всех возможных сочетаний значений параметров (у^у^у ) )х(у<2)'у22)'—у(22) )х •••х
х(у(п),у(п),...у(п)) на базе основного интегрального соотношения общего решения тепловой задачи
г
Т(х,/) = |С(х,У-т)д(т)ёт, (3)
где G(x,^,t-г) - функция Грина, соответствующая краевой задаче (1), (2) .
При постоянной плотности теплового потока q(t) = q0 = const выражение (3) принимает вид
T(x,t) = q0t + 2q0 £ (l-e^1)cos(nmx), (4)
m=1 П m \ '
а при экспоненциальной зависимости q (t) = q0 (l -e~pt) температура T (x,t) описывается соотношением
T (x,t ) = qot - | (l - e~p ) +
^l_e-^2m2t g-n2m2i_ -яЛ (5)
„ ~ / 1\m 1-e e -e / \
+2qo Z(-1) ,2 + ? 2 . cos(nmx)-n m n m - p
v ' /
Для построения адекватной модели процесса нестационарной теплопроводности использовалась двухслойная нейронная сеть с двадцатью пятью нейронами в скрытом слое. Процесс обучения ИНС, обеспечивающий нахождение значений синаптических весов нейронной сети, реализован на основе алгоритма обучения «байесовская регуляризация» с использованием сигмоидной функции активации. Оценка эффективности полученной нейросетевой модели проводилась на основе среднеквадратичного отклонения между зашумленной температурной зависимостью, построенной на основе (4) или (5), и модельным распределением температуры, в котором в качестве q(t)
использовалась постоянная или экспоненциальная функция, параметры которой получены с помощью ИНС.
Параметрическая идентификация граничного воздействия
в условиях неполноты экспериментальной информации
Решение ОЗТ, позволяющее провести оценивание параметров неизвестной характеристики на основе экспериментальной информации при априорно заданной структуре искомого воздействия, является задачей параметрической идентификации.
В качестве экспериментальной информации T' (x" ,t) использовались данные T * (x* ,t ) = T (x* ,t ) + £• randn, полученные на основе выражения (4) или (5) в точке x" на интервале идентификации t е 0,t в выборочные моменты времени
t е (tj,t2), ^ = 0;t^ = t* с учетом случайных возмущений, распределенных
по нормальному закону с нулевым математическим ожиданием и среднеквадратичным отклонением (СКО) 8.
Решение задачи в классе постоянных функций q(t) = q0 = const можно рассматривать как начальный этап для определения начального приближения для последующей итерационной процедуры поиска решений на множестве с большим числом параметров. В этом случае само значение q и является единственным искомым параметром у = q0.
Эффективность решения поставленной задачи существенно зависит от априорной информации об идентифицируемом воздействии. Задача решена в условиях интервальной неопределенности величины д , когда полагается, что имеются
сведения о границах интервала возможного изменения д е[^п, дтх ], что, как
правило, доступно на основе очевидных соображений или существующих ограничений. Обучающая выборка заполняется из допустимых значений
до е[д1,д2,-дк];д! =дтп,дк = дтах и значений температуры Т0(х",2,до), полученных согласно (4) в те же моменты времени 2 ,22 ,•••,2^) для всех возможных значений д . Таким образом, входные данные представлены в виде массива размерностью [# х к].
Идентификация тепловых процессов в случае, когда параметрическое представление искомой характеристики содержит единственный параметр у , является частным случаем общей многопараметрической задачи [16]. Возможности применения ИНС для решения граничной ОЗТ, когда параметрическая форма д(2) описывается двумя параметрами, были исследованы при поиске идентифицируемой функции в классе экспоненциальных зависимостей
д(г) = до (1-) = у (1-в722 ),у = (у1?у2), (6)
учитывающих инерционность тепловых процессов.
Одним из распространенных и универсальных подходов к решению многопараметрической задачи является автономная идентификация каждого из п параметров на основе построения п нейронных сетей с одним выходом каждая. Для учета взаимного влияния параметров друг на друга может быть реализовано рекурсивное поэтапное выполнение алгоритма автономной идентификации параметров, при котором на каждом очередном этапе на вход НС подаются значения параметров, полученные на предшествующем этапе, и соответствующие им модельные температурные распределения. Последовательный расчет этапов завершается, когда ошибка аппроксимации экспериментальных данных модельной зависимостью становится меньше некоторого априорно заданного значения функционала идентификации, формируемого, как правило, на основе абсолютного или среднеквадратичного значения температурной невязки.
Таким образом, решение сформулированной ОЗТ в классе функций (6) сводилось к последовательному расчету двух параметров д0 и р . Первый этап расчета д выполняется вышеописанным способом при зафиксированном на основе априорной информации значении р'0'. Найденное в результате первого этапа значение д() фиксируется и используется на следующем этапе для формирования обучающих данных ИНС, восстанавливающей значение Р. В этом случае формируется выборка допустимых значений искомого параметра Ре[Р,Р2 ,-Рк];Р =Рт1п,Рк =Рщах, аналогично первому этапу, и соответствующих
им модельных температурных кривых Т° (х ,2, д2гР<г)) . Завершение итерационной
процедуры происходит на основе абсолютного или относительного среднеквадратичного значения температурной невязки в соответствии с выполнением усло-
ВИй £(0 <8^ или )
•(') _ е(1 1 Л < ^ОН. При этом найденные на очередном I -ом этапе расчета значения %г) и могут быть использованы для сужения области допустимых решений на следующем I + 1-ом этапе.
Вычислительный эксперимент и полученные результаты
Решение задачи в классе постоянных функций
Изложенная методика была опробована на следующих модельных данных. Искомая постоянная величина плотности теплового потока, используемая для
расчета «экспериментальной» температуры, была принята % = 0.97. Для оценки временного интервала идентификации использовались относительные единицы t е[0,0.01,0.02, ...1], где шаг дискретизации составлял & = 0.01, а размер выборки
по временной переменной N=101. Предполагалось, что априорная информация о границах изменения искомой величины была известна лишь приближенно, в связи с чем допустимый интервал был принят как % е[^тт,9ти]5 9тп = 0,Ят,* = 5, что соответствует более чем пятикратному превышению истинного значения. В зависимости от интенсивности помехи шаг дискретизации был принят от = 0.05 (при 5 = 0%) до = 0.01 (при £=15%), что соответствовало от к=101 до к = 501 значению обучающей выборки.
Результаты одного из возможных решений задачи с применением нейросете-вых технологий, демонстрирующие их свойства, представлены в табл. 1. Здесь полученное на основе нейронной сети значение искомой величины плотности теплового потока обозначено как %ИНС) . Среднеквадратичная погрешность температурного распределения рассчитывалась по формуле
£ =
N 1(г (XV )-т0 (XV. %0ИНС> ))2,
(7)
где Т0 (XV, %(0ИНС)) соответствовало найденному значению %(ИНС).
Результаты показывают влияние уровня случайного возмущения в каждом из опытов, но при этом отражают сглаживающие свойства нейросетовых моделей в условиях действия неопределенных факторов. Восстановленные значения
%0
(ИНС )
в серии экспериментов зависят от случайной составляющей, при этом зна-
чения среднеквадратичной погрешности аппроксимирующей температурной кривой отражают зависимость от уровня СКО возмущающего воздействия.
Таблица 1
Точность параметрической идентификации в классе постоянных функций в зависимости от величины СКО погрешности измерения
5, % 0 1 5 10 15
(ИНС ) %0 0.9700 0.9725 0.9305 0.9993 0.9642
е 2.6988е-11 1.0909е-04 0.0027 0.0094 0.0215
Решение задачи в классе экспоненциальных функций
В данном модельном примере искомые значения параметров были приняты следующими: д* = 0.97; у3=47, временной интервал соответствовал предыдущему
случаю. На основе априорной информации об экспоненциальной зависимости искомой плотности теплового потока выбраны границы интервала изменения параметра фе[4,800] как предельные значения, определяющие требуемый характер зависимости на заданном временном интервале (рис. 1). В случае ф< 4 зависимость (6) на интервале ^ е[0,1] не выходит на установившееся значение, а при 3 > 800 участок изменения д^) является столь малым, что им можно пренебречь и искать д^) в классе постоянных функций. В зависимости от интенсивности помехи полученный диапазон был поделен с использованием шага дискретизации от Д3=4 до Д3=1, на основе чего был сформирован входной массив значений у и соответствующих им температурных распределений Г0 (XV,3) для обучения второй ИНС. Начальное приближение ф(0) выбрано равным 10. Характер кривых для всех трех случаев показан на рис. 1.
Рис. 1. Экспоненциальный характер д^) для предельных значений 3 -(1), 3тх "(2) и среднего етачения ф - (3)
На каждом из расчетных этапов точность определения каждого из параметров / = 1,2 существенно зависит от априорно задаваемого значения другого параметра, что особенно сказывается при возрастании уровня возмущающего воздействия. При отсутствии априорной информации для обеспечения сходимости решения к искомому значению рекомендуется сочетание нейросетевых технологий с численными алгоритмами, устанавливающими дополнительные ограничения на принадлежность искомых параметров определенным диапазонам на основе анализа поведения температурного поля исходя из его качественных свойств. Такой подход позволяет получить адекватные результаты при решении
двухпараметрической задачи (табл. 2, рис. 2-4). Итерационное выполнение указанных двух этапов по идентификации параметров д0 и ф проводилось до выполнения условия £(,) -£(,-1) <£доп .
Рис. 2. Зависимость точности решения £,% от уровня СКО случайного возмущения
1.4
1.2
1
0.8 ^ 0.6 0.4 0.2 0 -0.2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
í
Рис. 3. Зашумленная температурная кривая Т (х ,/) при £ = 10%, используемая в качестве экспериментальной информации - (1),
и полученное модельное распределение
Т0 (х>, д(0ИНС),ф( ИНС)) - (2)
0.25 0.2 0.15 0.1 0.05
0
-0.05
0 0.2 0.4 0.6 0.8 1
í
Рис. 4. Погрешность д^^^) -при 5=0%-1,5=3%-2,5 = 5%-3,5=10%-4
Таблица 2
Точность параметрической идентификации в классе экспоненциальных функций в зависимости от величины СКО погрешности измерения
5, % 0 1 3 5 10
(ИНС ) %0 0.9696 0.9730 0.9632 0.9765 0.9882
р(ИНС) 48.4582 63.4085 59.0270 64.4264 82.9217
е 2.4076е-07 1.4406е-04 9.5059е-04 0.0023 0.0112
Выводы
В работе предложен подход к нейросетевому решению задачи параметрической идентификации при заданной структуре математического оператора искомой характеристики. Проведенные исследования подтверждают возможность использования интеллектуальных методов для решения обратных задач математической физики в условиях неполноты информации в ряде случаев, когда искомая зависимость может быть аппроксимирована функцией с небольшим числом параметров.
Качество работы алгоритма существенно зависит от погрешности исходных данных. С увеличением числа параметров при действии возмущающих факторов для получения удовлетворительной точности решения нейросетевые технологии рекомендуется дополнять алгоритмами, контролирующими поведение искомых и измеряемых характеристик на основе их базовых свойств. Использование априорной информации позволяет повысить точность решения и уменьшить вычислительную сложность.
А М
Р/
I 1
; 1
¡Л I
Л*
Предложенный подход может быть в дальнейшем распространен на структурные операторы более сложной формы, характеризующиеся вектором параметров большей размерности.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Алифанов О.М. Обратные задачи теплообмена. М.: Машиностроение, 1988. 280 с.
2. Бек Дж., Блакуэлл Б., Сент-Клэр Ч. мл. Некорректные обратные задачи теплопроводности. М.: Мир, 1989. 312 с.
3. Алифанов О.М., Артюхин Е.А., Румянцев С.В. Экстремальные методы решения некорректных задач. М.: Наука, Гл. ред. физ.-мат. лит, 1988. 288 с.
4. Дилигенская А.Н., Рапопорт Э.Я. Аналитические методы параметрической оптимизации в обратных задачах теплопроводности с внутренним тепловыделением // Инженерно-физический журнал. 2014. Т. 87. № 5. С. 1082-1089.
5. Diligenskaya A.N. Methods of Sequential Parametric Optimization in Inverse Problems of Technological Thermophysics // 2019 XXI International Conference Complex Systems: Control and Modeling Problems (CSCMP). Samara, Russia, 2019. pp. 267-270.
6. Прангишвили И.В., Потоцкий В.А., Гинсберг К.С., Смолянинов В.В. Идентификация систем и задачи управления: на пути к современным системным методологиям // Проблемы управления. 2004. № 4. С. 2-15.
7. Воронова Н.П., Ковалев С.М., Шабельников А.Н. Идентификация и оценивание состояний нечетких динамических систем // Известия ЮФУ. Технические науки. 2016. № 6 (179). С. 128-130.
8. Бердник В.В., Мухамедяров Р.Д. Применение метода нейронных сетей для решения обратной задачи теплопереноса // ТВТ. 2003. Т. 41. № 6. С. 942-947.
9. Круглов В.В., Дли М.И., Голунов Р.Ю. Нечеткая логика и искусственные нейронные сети. М.: Физматлит, 2001. 221 с.
10. Хайкин С. Нейронные сети: полный курс. 2-е изд. М.: И.Д. Вильямс, 2006. 1104 с.
11. Аксёнов С.В., Новосельцев В.Б. Организация и использование нейронных сетей (методы и технологии). Томск, 2006. 128 c.
12. Лысенко Д.С., Данилушкин И.А. Динамическая модель котла-утилизатора на базе рекуррентной нейронной сети // Вестник Самарского государственного технического университета. Сер. Технические науки. 2020. Т. 28. № 2. С. 59-72.
13. Замятин Н.В., Медянцев Д.В. Методика нейросетевого моделирования сложных систем // Известия Томского политехнического университета. 2006. Т. 309. № 8. С. 100-106.
14. Gorbachenko V.I., Lazovskaya T.V., Tarkhov D.A., Vasilyev A.N., Zhukov M.V. Neural Network Technique in Some Inverse Problems of Mathematical Physics. In: Cheng L., Liu Q., Ronzhin A. (eds) Advances in Neural Networks - Lecture Notes in Computer Science, 2016. V. 9719.
15. Pavlov D.Yu. Neural networks in solution of inverse coefficient heat conduction problem // Вестник МГУ. Серия 15: ВМК. 1994. № 4. С. 44-51.
16. Eremenko Y.I., Poleshchenko D.A., Tsygankov Y.A. Development of neural network model of the multiparametric technological object. J. Fundam. Appl. Sci., 2017. 9 (7S). Рр. 706-721.
Статья поступила в редакцию 3 октября 2020 г.
PARAMETRIC IDENTIFICATION IN INVERSE HEAT CONDUCTION PROBLEMS UNDER INTERVAL UNCERTAINTY BASED ON NEURAL NETWORKS
A.N. Diligenskaya, A.V. Samokish
Samara State Technical University
244, Molodogvardeyskaya str., Samara, Russian Federation
Abstract. The application of parametric identification methods used to solve inverse problems of technological thermophysics under conditions of interval uncertainty of parameters is considered. To solve such problems under the action of disturbing factors, artificial neural networks can be used that reproduce the structure of the operator of the direct problem with a small dimension of the desired vector of parameters. A variant of constructing a neural network for solving the boundary inverse heat conduction problem with a known structure of a mathematical operator, characterized by a small number ofparam-eters, is proposed. The method allows the use of a priori information about the admissible range of belonging of the identified characteristics or their parameters. The input data for the neural network model are all possible realizations of temperature states, obtained as a reaction to the input action, which is an identifiable characteristic that satisfies all permissible options for the combination ofparameters.
The solution of the problem is an iterative procedure, at each step of which, using a neural network, the value of one of the parameters specifying the parametric representation of the desired value is refined. The combination of neural network technologies with high-quality information about the regularities of non-stationary heat conduction processes makes it possible to increase the accuracy of solving the problem ofparametric identification under conditions of disturbances. Variants of solving the boundary inverse problem of heat conduction with one- and two-parameters representation of the desired characteristic under the conditions of disturbances are demonstrated. The accuracy of solving the problem is investigated depending on the intensity of the disturbing effect.
Keywords: inverse heat conduction problem, parametric identification, interval uncertainty, neural network model
REFERENCES
1. Alifanov O.M. Obratnye zadachi teploobmena [Inverse Heat Transfer Problems]. Moscow: Mashi-nostroenie, 1988, 280 p. (In Russian).
2. Beck J.V., Blackwell B., C.R. Jr St. Clair. Inverse Heat Conduction: Ill-Posed Problems, New York: Wiley, 1985.
3. Alifanov O.M., Artyukhin E.A., Rumyantsev S.V. Jekstremal'nye metody reshenija nekorrektnyh zadach [Extremum Methods of Solution of Ill-Posed Problems and Their Application to Inverse Heat-Transfer Problems]. Moscow: Nauka, 1988, 288 p. (In Russian).
4. Diligenskaya A.N., Rapoport E.Ya. Analytical methods of parametric optimization in inverse heat-conduction problems with internal heat release // J. Eng. Phys. Thermophys., 2014. Vol. 87, no. 5, p. 1126.
5. Diligenskaya A.N. Methods of Sequential Parametric Optimization in Inverse Problems of Technological Thermophysics // 2019 XXI International Conference Complex Systems: Control and Modeling Problems (CSCMP). - Samara, Russia, 2019. pp. 267-270.
6. Pranghishvili I.V., Lototsky V.A., GhinsbergK.S., Smolyaninov V.V. System identification and control problems: on the way to modern system methodologies // Control Sciences. 2004. № 4. Pp. 2-15 (In Russian).
Anna N. Diligenskaya (Dr. Sci. (Techn.)), Professor. Alexandr V. Samokish, Postgraduate Student.
7. Voronova N.P., KovalevS.M., Shabelnikov A.N. Identification and estimation of states for fuzzy dynamical systems // Izvestiya SFedU. Engineering Sciences. 2016. № 6 (179). Pp. 128-130 (In Russian).
8. Berdnik V.V., Mukhamedyarov R.D. Application of the Method of Neural Networks to Solution of the Inverse Problem of Heat Transfer // High Temperature, 2003, V. 41, pp. 839-843.
9. Kruglov V.V., Dli M.I., Golunov R.Y. Nechetkaja logika i iskusstvennye nejronnye seti [Fuzzy logic and artificial neural networks]. Moscow: Fizmatlit, 2001. 221 p. (In Russian).
10. Haykin S. Nejronnye seti: polnyj kurs. 2-e izd [Neural networks: full course. 2nd ed.], Rev. M.: I.D. Williams, 2006. 1104 p. (In Russian).
11. Aksenov S.V., Novoseltsev V.B. Organizacija i ispol'zovanie nejronnyh setej (metody i tehnologii) [Organization and use of neural networks (methods and technologies)]. Tomsk, 2006 (In Russian).
12. Lusenko D.S., Danilushkin I.A. Dynamic model of waste heat boiler based on recurrent neural network // Vestn. Samar. Gos. Tehn. Un-ta. Tehn. Nauki, 2020. V. 28, No. 2. Pp. 59-72 (In Russian).
13. Zamyatin N. V., Medyantsev D.V. Methods of neural network modeling of complex systems // Izves-tia Tomsk. Pol-teh Un-ta, 2006. No. 8. Pp. 100-106. (In Russian).
14. Gorbachenko V.I., Lazovskaya T.V., Tarkhov D.A., Vasilyev A.N., Zhukov M.V. Neural Network Technique in Some Inverse Problems of Mathematical Physics. In: Cheng L., Liu Q., Ronzhin A. (eds) Advances in Neural Networks - Lecture Notes in Computer Science, 2016. V. 9719.
15. Pavlov D.Yu. Neural network solution of the coefficient inverse heat conduction problem // Vestn. Mosk. Univ., Ser. 15: Vych. Mat. Kibern., 1994, no. 4, pp. 51-56.
16. Eremenko Y.I., Poleshchenko D.A., Tsygankov Y.A. Development of neural network model of the multiparametric technological object // J. Fundam. Appl. Sci., 2017. 9 (7S). Pp. 706-721.
