Вестник КРАУНЦ. Физ.-мат. науки. 2018. № 2(22). C. 33-44. ISSN 2079-6641
DOI: 10.18454/2079-6641-2018-22-2-33-44
ИНФОРМАЦИОННЫЕ И ВЫЧИСЛИТЕЛЬНЫЕ ТЕХНОЛОГИИ
УДК 519.622.2
ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧИ ПАРАМЕТРИЧЕСКОЙ ИДЕНТИФИКАЦИИ ДЛЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ
В.М. Абдуллаев1,2
1 Азербайджанский государственный университет нефти и промышленности, AZ1010, г. Баку, пр-т Азадлыг, 20
2 Институт систем управления НАН Республики Азербайджан, AZ114, г. Баку, ул. Вахабзаде, 9
E-mail: vaqif_ab@rambler.ru
В статье предлагается численный метод решения задачи, основанный на использовании метода прямых для сведения задачи к системе обыкновенных дифференциальных уравнений с неизвестными параметрами. Далее используется специальное представление решения полученной краевой задачи относительно линейной системы дифференциальных уравнений с нелокальными условиями, с помощью которого задача параметрической идентификации сводится к решению вспомогательных краевых задач и одной системы алгебраических уравнений.Важно отметить, что, в отличие от оптимизационных подходов, в данной работе не используется построение каких-либо итерационных процедур или минимизирующих последовательностей.
Ключевые слова: обратная задача, метод прямых, параболическое уравнение.
© Абдуллаев В.М., 2018
Введение
В данной работе предлагается подход к численному решению задач параметрической идентификации относительно дифференциальных уравнений с частными производными. Для идентификации постоянных во времени или по пространственной переменной имеются результаты дополнительных экспериментов, в процессе которых проводились замеры состояния объекта.
Отметим, что с подобными обратными задачами приходится сталкиваться на этапе параметрической идентификации математических моделей практически для всех динамических процессов, для которых предполагается строить системы автоматического или автоматизированного управления. В связи с этим различным аспектам исследования коэффициентно - обратных задач посвящено большое число публикаций [1-14].
Наиболее часто используемым методом параметрической идентификации является приведение обратной задачи к задаче параметрического оптимального управления и применение для ее решения численных методов первого порядка. Проблема применения такого подхода связана с необходимостью построения итерационных процедур, например, градиентного спуска, на каждой итерации которой необходимо решать саму исходную и сопряженную систему дифференциальных уравнений [4,12].
В данной работе предлагается численный метод решения задачи, основанный на использовании метода прямых для сведения задачи к системе обыкновенных дифференциальных уравнений с неизвестными параметрами. Далее используется специальное представление решения полученной краевой задачи относительно линейной системы дифференциальных уравнений с нелокальными условиями, с помощью которого задача параметрической идентификации сводится к решению вспомогательных краевых задач и одной системы алгебраических уравнений.Важно отметить, что, в отличие от оптимизационных подходов, в данной работе не используется построение каких-либо итерационных процедур или минимизирующих последовательностей.
Приводятся результаты численных экспериментов и их анализ.
Постановка задачи
Рассмотрим следующую задачу восстановления коэффициентов для параболического уравнения:
= < (х, г) + < (х,г)^Ц^Хт - <2(х, г)и(х, г) + (х, г; С) + /(х, г),
(х, г) е О = {(х, г) :0 < х < а, 0 < г < Т}. (1)
Здесь функция ^(х,г; С) имеет один из видов:
I
^ (х, г; С) = £ Бг (х, г) С(), (2)
г'=1 I
^ (х, г; С) = £ Бг(х, г) С,-(х), (3)
г'=1
£(х,г) > 0, ^(х,г), £2(х,г), f (х,г) е С2,1(П)- заданные непрерывные функции, £г-(х,г), г = 1,2, ...1, заданные непрерывные линейно-независимые функции. Выполнены также условия ограниченности в П значений функций: 0 < а1(х,г) < с1, |а2(х,г)|, |а3(х,г)| ,
дaз(x,г)/дх < с2.
С целью идентификации С(г) или С(х) имеются результаты N экспериментов, проведенных при различных начальных и краевых условиях:
и* (х, 0) = ф* (х), 0 < х < а, (4)
и* (0, г) = у* (г), и*(а, г) = у{(г), 0 < х < Т, (5)
здесь функция ф*(х), у0(г), у1 (г) непрерывны по своим аргументам, * = 1,2,...,N.
При каждом эксперименте проводились наблюдения за состоянием процесса. Обозначим через и*(х,г) состояние процесса в точке х в момент времени г при том эксперименте, * = 1,2,...,N.
Дополнительные условия, необходимые для определения вектор - функции С(г), полученные по результатам наблюдений за процессом при каждом *-том эксперименте, могут иметь различный вид, что связано с характером проводимых экспериментов и результатов наблюдения [3,11,14]. В частности, пусть в точках хг е (0;а), г = 1,2, ..X, ведутся наблюдения:
и*(хг,г)= ф*(г), г = 1,2,..., Ц, 0 < г < Т, (6)
где Ф*(г), г = 1,2,..., /- заданные непрерывно дифференцируемые функции.
А для определения вектор - функции С(х) = (С1(х), С2(х),..., С/ (х))* дополнительные условия являются результатом наблюдения заданные моменты времени, г* е (0,Т], г = 1,2,...,Ц:
и*(х, гг)= Ф* (х), г = 1,2,..., Ц, 0 < х < а. (7)
Здесь: Ф*(х), г = 1,2,..., /- заданные дважды непрерывно дифференцируемые функции.
Кроме того функции ф0(х), ф{(х) у/(г), у2(г) удовлетворяют условиям согласования:
ф0 (0) = у1 (0), ф0 (а) = (0), ф0 (0) = Ф* (0).
Задача состоит в нахождении непрерывной вектор-функции С(г), удовлетворяющей условиям (1),(2),(4)-(6) (задача А) или вектор-функцию С(х) в случае условий (1),(3),(4)-(5), (7) (задача В).
Рассматриваемая задача идентификации исследовалась многими авторами [1-14]. Были получены необходимые условия существования и единственности решения [5-7]. Предложенные в работах [3,4] методы решения этой задачи использовали сведение ее к задаче оптимального управления, для решения которой применялись итерационные методы.
В работах [8,9] использовалось сведение задачи к интегральному уравнению, решение которого представляет вычислительную сложность. Ниже предлагается использовать приведение задачи (17)-(24) к задаче параметрической идентификации относительно системы обыкновенных дифференциальных уравнений.
Метод решения
В данной работе предлагается численный метод решения задачи А и задачи В основанный на использовании метода прямых для сведения задачи к системе обыкновенных дифференциальных уравнений с неизвестными параметрами.
Решение задачи А. С целью приведения задачи А к обыкновенной система дифференциальных уравнений используем метод прямых. В области О проведем прямые г = гк = кт, к = 0,1,...,М, т = Т/М.
Производную ^х,г^ |г=гк в (1) аппроксимируем разностным отношением
Эф,?), ф,^)-ФЛ-ъ) , п^ ^ 1 о л*
дг |г=гк =-Т-+ °(т) к = 1,2,...,M,
и далее используем обозначения:
и(к)(х) = И(х,гк), С(к) = С(гк), ^(к) = ^ = У20к),
'<2(х, гк) 1
B (X)- § (x, tk)' * (x)- § (x, tk) ' (x)- l
х,г^' ^ <(х,гк)' ^ (х,гк) т■ <(х,гк)
/(к)(х)_ (/М + и(к-1)(х)^
/ (х) ^ < (х, гк) + т ■ < (х, г^ .
В результате получим уравнения второго порядка с обыкновенными производными:
*2Ц(к2(х) = 1« (х) ^^ +1« (х) и (к) (х)+
их
+ ^Бг(к)(х)с(к) + /(к)(х), к = 1,2,...М, и(0)(х) = ф(х). (8)
г'=1
Из (5) имеем следующие условия:
и(к)(0) = ^(к), и(к)(а) = ^2к), к = 1,2,...М, (9)
а дополнительные условия имеют вид
и(хг)= Ф(к), г = 1,2,..., I, (10)
Уравнения (8) при каждом к, к = 1,2,...,М, приведем к системе двух дифференциальных уравнений первого порядка:
^ = и2к)(х), ^ = <2к) (х) и1к)(х) +< « (х) и2к)(х)+
+ ¿Б(к)(х)С(к) + /(к)(х), к = 1,2,...М, (11)
г'=1
и{к)(0) = ^(к), и(к)(а) = ^2к), к = 1,2,...М. (12)
Для простоты задачи (11)-(12),( 10), при каждом к, к = 1,2,...,М, запишем в следующем виде:
^^ = А(х)и(х)+ В(х)С + ^(х),х е (0,а], (13)
ах
с краевыми условиями
а и (0) = у, в и (а) = у2. (14)
а дополнительные условия имеют вид
уи(х,-) = Ф,-, г = 1,2,...,Ц, (15)
Здесь и(х) е Я2-фазовое состояние системы; С = (С1,С2,...,С/)* е Я1- искомые параметры; А(х), В(х),^(х)-заданные матричные непрерывные функции по х, х е (0,а] соответственно размерности (2 х 2), (2 х /), (2 х 1), а, /3, '/-заданные матрицы размерности (1х 2). *-знак транспонирования.
Решение задачи (13), (14) -и*(х) для каждого того эксперимента будем искать в виде:
/
и*(х) = и0* (х) + £ иг*(х)Сг, х е (0, а* , * = 1, 2,...,N, (16)
г=1
где пока произвольные вектор-функция и0*(х) и иг*(г), должны удовлетворять условиям
а и0* (0) = у1, /3 и0* (а) = у2, а и* (0) = 0, /3 и''* (а) = 0, (17)
Несложно проверить, что в этом случае и*(х) из (16) удовлетворяет условиям (14) для всех * = 1,2,...,N и произвольных Сг, г = 1,2....,/.
Теорема 1. £сли функции и**(х), г = 0,1,...,/ являются при хе (0,а* решениями следующих задач:
аи°*(х) = А(х)и0*(х) + ^(х), аи0*(0) = у1, /3и0*(а) = у2, (18)
dUij (x) dx
= A(x)Uij(x) + B'(x), а(0) = 0, /3(a) = 0,
(19)
тогда функция и* (х), определенная формулой (16), удовлетворяет условиям (13), (14) для произвольных значений вектора С, * = 1,2,...,N.
Доказательство. Действительно, продифференцировав (16) и подставив в (13), после несложных преобразований получим:
- A(x)U0 j (x) - F (x)
dx
+
+ E
i=1
dUij (x) dx
- A(x)Uij (x) - B'(x)
Ci = 0 ,x e (0,aj] , j = 1,2,...,N.
Учитывая произвольность функций и0*(х) и иг*(х), г = 1,2....,/, требуя от них равенства нулю выражений, стоящих в квадратных скобках, получим, что эти функции должны быть решениями задач (18), (19). □
Для всех видов измерения будем предполагать, что проводимые эксперименты независимы, причем выполнено условие L > I . Дополнительные условия (15) называют условиями переопределения [3].
Используя значения полученных решений задач (18), (19) в точках, которые участвуют в проводимых измерениях вида (15), получим систему L линейных алгебраических уравнений относительно вектора C € Rl, которую в общем случае запишем так
QC = а (20)
Здесь матрица Q = (^)), г = 1,...,^ ] = 1,...,l, вектор G = g2,...,gL)* определяются видом и результатами проводимых наблюдений. В случае, если L = I, то решение определяется непосредственно из (20): С = Q-1G. В случае, если L > I под решением (20) будем понимать вектор:
С Q*G, (21)
называемый нормальным решением переопределенной системы ([3]).
Решение задачи В. В области П проведем прямые х = хк = кк, к = 0,1,...,М, к = а/М.
д /чч д2и(х,г) | ди(х,г) |
Аппроксимируя в (1) производные дХ2 |Х=Хк и дх |Х=Хк разностными отношениями
^ |х=хк = Фк+Ьг) - и(хк-1,') + 0(к2), к = 1,2,...,М- 1, (УХ 2к
д^фО, и(хк+ьг) - 2и(хк,г) + и(хк-ьг) . , 0 м л дх2 |х=хк = -к2-+ 0(к ), к = 1, 2,...,М - 1,
и далее используем обозначения:
и(к)(г) = и(хк,г),В(к)(г) = В(хк,г), % = Ф(хк), С(к) = С(хк), §(к)(г) = §(хк,г), §1(к)(г) = §1 (хк,г), §2к)(г) = §2(хк,г), /(к)(г) = /(хк,г), к = 1,2,...,М- 1.
В результате получим систему дифференциальных уравнений с обыкновенными производными:
= ^ (и(к+1)(г) - 2и(к)(г) + и(к-1)(г)) + (и(к+1}(г) - и(к-1)(г)) +
+§2(к) (г)и(к) (г) + £ в(к) (г)с(к) + /(к) (г), к = 1, 2,..., м - 1,
г=1
и(0)(г) = ^(г), и(М)(0 = ^(0, г € (0;Г], (22)
с условиями
и (к)(0) = ф(к), к = 1,2,..., М - 1. (23)
После некоторых обозначений получим задачу:
) = A(t)U(k)(t) + B(k)(t)C + F(t), k = 1,2,...,M- 1, U(k)(0) = ф
dt
Дополнительные условия (7) примут вид
и^(г,-) = Фг(к)*, г = 1,2,..., /. (25)
Задача (24) при наблюдениях (25) является частным случаем рассмотренной выше задачи параметрической идентификации относительно (13), (14) при каком-либо одном из условий (13)-( 14), т.к. вместо краевых условий имеем начальные условия Коши (23). Поэтому для ее численного решения можно использовать предложенной подход.
Результаты численных экспериментов
Были проведены многочисленные численные эксперименты на тестовых задачах с применением предложенных в данной работе формул и схем численного решения. Результаты экспериментов показали достаточную практическую эффективность описанного подхода.
Задача 1. Рассмотрим следующую задачу параметрической идентификации для параболического уравнения (1):
^ = ^^ (2х+г) ^ЦМ - 2и(х,,) + (х+, )С, (г) + (2х + 3, )С2(г)+ д г дх2 дх
+пе-2г (п8ш(пх) - (х + 2г)ео8(пх)) + (х + г)(1.5е-2г + 0.5) - 0.25(2х + 3)е3г,
(x, t) е П = {(x, t) : 0 < x < 1, 0 < t < 1}. (26)
Здесь "(x,t)- состояние процесса, C(t) е R2 - идентифицируемый вектор параметров, причем функция u(x,t) = (e-2t sin(nx)) и вектор-функция C(t) = (|e-2t + 4e3t) , как несложно проверить, удовлетворяют уравнениям (26).
Для идентификации вектора C(t) проведен один эксперимент (N = 1) при начально-краевых условиях:
u(x,0) = sin(nx), 0 < x < 1, u(0,t) = 0, u(1,t) = 0, 0 < t < 1, (27)
при котором проводились наблюдения в двух точках x = 0.25 и x = 0.5 (т.е. L = l = 2) и в результате были получены дополнительные разделенные многоточечные условия вида (6):
л/2 e_2t
u(0.25,—, "(0.5,t) = e-2t, 0 < t < 1. (28)
Для сведения задачи (26), (27) к задаче (13)—( 14) методом прямых при проведении численных экспериментов использовались различные значения шага т (т.е. число прямых). Решения для т е [0.02;0.04] практически не различались и были достаточно близки к точному решению. В таблице 1 приведены результаты, полученные при т = 0.04, т.е. M = 25.
В табл.1 приведены точные и полученные значения параметров C(t) при х = 0; 0.01; 0.03, что соответствует замерам соответственно без наличия помех, при помехах 1% и 3% от измеряемой величины в условиях (28).
Задача 2. Пусть относительно процесса (26), (27) дополнительно к наблюдениям (28) проводились наблюдения и в точках х = 0.125 и х = 0.75:
ц(0.125,г) = е-2г 8т(0.125п) , и(0.75,г) = е-2г 8т(0.75п), 0 < г < 1, (29)
т.е. в данной задаче N = 1,Ц = 4, / = 2, Ц > /. Параметры численных методов были аналогичны параметрам, используемым в задаче 1. Ясно, что точные решения задач 1 и 2 совпадают. В табл.2 приведены полученные результаты решения задачи 2 при различных уровнях помех на измеряемые величины (28), (29) состояния процесса.
Были проведены многочисленные другие численные эксперименты, приведение результатов которых заняло бы много места. Отметим лишь, что эти эксперименты показали возможность получения решения задач предлагаемым методом с требуемой высокой степенью точности и достаточно высокую его устойчивость к помехам в исходных данных задачи.
Из результатов, приведенных в таблицах, следует, что предлагаемый подход позволяет достаточно точно определять значения идентифицируемых параметров, если дополнительная информация о состоянии процесса известна точно.
В случае наличия помех при проведении замеров, как и следовало ожидать. полученные значения параметров точно соответствуют полученной искаженной информации, а следовательно не соответствуют искомым значениям. С целью более точного определения искомых значений идентифицируемых параметров необходимо увеличение точности замеров или количества информации, т.е. точек или моментов времени замера параметров состояния процесса.
Таблица 1.
Точные и полученные решения Задачи 1 при разных уровнях помех
Точные X = 0 X = 0.01 X = 0.03
t Ci(t) C2(t) Ci(t) C2(t) Ci(t) C2(t) Ci(t) C2(t)
0.02 1.9412 0.5309 1.9408 0.5308 1.9266 0.5269 1.8835 0.5151
0.10 1.7281 0.6749 1.7278 0.6748 1.7151 0.6698 1.6768 0.6549
0.20 1.5055 0.9110 1.5052 0.9108 1.4942 0.9042 1.4588 0.8828
0.30 1.3232 1.2298 1.3229 1.2296 1.3133 1.2206 1.2862 1.1954
0.40 1.1740 1.6600 1.1738 1.6597 1.1652 1.6476 1.1376 1.6085
0.50 1.0518 2.2408 1.0516 2.2404 1.0439 2.224 1.0181 2.1691
0.60 0.9518 3.0248 0.9516 3.0242 0.9447 3.0021 0.9232 2.9341
0.70 0.8699 4.0830 0.8697 4.0822 0.8634 4.0524 0.8464 3.9728
0.80 0.8028 5.5116 0.8026 5.5105 0.7968 5.4703 0.7795 5.3518
0.90 0.7479 7.4387 0.7478 7.4372 0.7423 7.3829 0.7255 7.2155
1.00 0.7030 10.0427 0.7029 10.0424 0.6977 9.9674 0.6812 9.7314
Таблица 2.
Полученное решение Задачи 2 при разных уровнях помех
X = 0 X = 0.01 X = 0.03
t Ci(t ) C2(t ) Ci(t ) C2(t ) Ci(t ) C2(t )
0.02 1.9408 0.5308 1.9292 0.5276 1.8907 0.5176
0.10 1.7278 0.6748 1.7174 0.6704 1.6849 0.6574
0.20 1.5052 0.9108 1.4962 0.9044 1.4618 0.8864
0.30 1.3229 1.2296 1.315 1.2224 1.2822 1.1929
0.40 1.1738 1.6597 1.1667 1.6492 1.1423 1.6168
0.50 1.0516 2.2404 1.0453 2.2271 1.0245 2.187
0.60 0.9516 3.0242 0.9459 3.0057 0.9261 2.9401
0.70 0.8697 4.0822 0.8645 4.0565 0.8499 3.985
0.80 0.8026 5.5105 0.7978 5.4752 0.7835 5.3463
0.90 0.7478 7.4372 0.7433 7.3933 0.7285 7.2453
1.00 0.7029 10.0424 0.6986 9.9824 0.6854 9.7816
Заключение
В статье исследуется численное решение коэффициентно-обратных задач относительно параболического уравнения. Предложен подход к решению рассматриваемых задач, основанный на использовании метода прямых для сведения задачи к системе обыкновенных дифференциальных уравнений с неизвестными параметрами. Далее используется специальное представление решения полученной краевой задачи относительно линейной системы дифференциальных уравнений с краевыми условиями, с помощью которого задача параметрической идентификации сводится к решению вспомогательных краевых задач и одной системы алгебраических уравнений. Приводятся результаты численных экспериментов и их анализ.
Список литературы
[1] Камынин В. Л., "Об обратной задаче определения правой части в параболическом уравнении с условием интегрального переопределения", Матем. Заметки, 77:4 (2005), 522-534. [Kamynin V. L., "Ob obratnoj zadache opredelenija pravoj chasti v parabolich-eskom uravnenii s usloviem integral'nogo pereopredelenija", Matem. Zametki, 2005 vol 77, №4, 522-534].
[2] Belov Yu. Ya., "Inverse problems for parabolic equations", J. Inv. Ill Posed Problems, . 1:4 (1993), 283-305.
[3] Сергиенко И. B., Дейнека В. С., "Решение некоторых обратных задач теплопроводности для составной пластины с использованием псевдообратных матриц", Доклады НАН Украины, 2011, № 12, 28-34. [Sergienko I. B., Dejneka V. S., "Reshenie nekotoryh obratnyh zadach teploprovodnosti dlja sostavnoj plastiny s ispol'zovaniem psevdoobratnyh matric", Doklady NAN Ukrainy, 2011, № 12, 28-34].
[4] Hasanov A., Otelbaev M., Akpayev B., "Inverse heat conduction problems with boundary and final time measured output data", Inverse Probl. Sci. Eng., . 19 (2011), 895-1006.
[5] Ivanchov M.I., "Inverse problems for the heat-conduction equation with nonlocal boundary conditions", Ukrainian Mathematical Journal, 45:8 (1993), 1186-1192.
[6] Саватеев Е.Г., "О задаче идентификации коэффициента параболического уравнения", Сиб. матем. журн., 36:1 (1995), 177-185. [Savateev E.G., "O zadache identifikacii kojef-ficienta parabolicheskogo uravnenija", Sib. matem. zhurn. 1995, 36:1, 177-185].
[7] Прилепко А.И., Костин А.Б., "О некоторых обратных задачах для параболических уравнений с финальным и интегральным наблюдением", Матем. сб., 183:4 (1992), 4968. [Prilepko A. I., Kostin A.B., "O nekotoryh obratnyh zadachah dlja parabolicheskih uravnenij s final'nym i integral'nym nabljudeniem", Matem. sb., 183:4 (1992), 49-68].
[8] Yan L, Fu C.L, Yang F.L., "The method of fundamental solutions for the inverse heat source problem", Eng. Anal. Boundary Elements, 32 (2008), 216-222.
[9] Ismailov M.I, Kanca F, Lesnic D., "Determination of a time-dependent heat source under nonlocal boundary and integral overdetermination conditions", Appl. Math. Comput., 218 (2011), 4138-4146.
[10] Cannon J.R, Duchateau P., "Structural identification of an unknown source term in a heat equation", Inverse Probl., 14 (1998), 535-551.
[11] Aida-zade K.R., Abdullayev V.M., "Solution to a class of inverse problems for system of loaded ordinary differential equations with integral conditions", J. Inv. Ill Posed Problems, 24:5 (2016), 543-558.
[12] Abdullayev V. M. , Aida-zade K. R., "Optimization of loading places and load response functions for stationary systems", Comput. Math. Math. Phys., 57:4 (2017), 634-644.
[13] Abdullaev V. M. , Aida-zade K. R., "On the numerical solution of loaded systems of ordinary differential equations", Comput. Math. Math. Phys, 44:9 (2004), 1505-1515.
[14] Aida-zade K.R., "A Numerical method of restroring the parameters of a dynamic system", Cybern. Syst. Analysis., 40:3 (2004), 392-399.
Список литературы (ГОСТ)
[1] Камынин В.Л. Об обратной задаче определения правой части в параболическом уравнении с условием интегрального переопределения // Матем. Заметки. 2005. Т.77. №4. С.522-534.
[2] Belov Yu.Ya., Inverse problems for parabolic equations // J. Inv. Ill Posed Problems. 1993. vol. 1. no. 4. pp.283 - 305.
[3] Сергиенко И.В., Дейнека В.С. Решение некоторых обратных задач теплопроводности для составной пластины с использованием псевдообратных матриц // Доклады НАН Украины. 2011. no. 12. С.28-34.
[4] Hasanov A., Otelbaev M., Akpayev В. Inverse heat conduction problems with boundary and final time measured output data // Inverse Probl. Sci. Eng., 2011. vol. 19. pp. 895-1006.
[5] Ivanchov M.I. Inverse problems for the heat-conduction equation with nonlocal boundary conditions // Ukrainian Mathematical Journal. 1993. vol. 45. № 8. С.1186-1192.
[6] Саватеев Е.Г. О задаче идентификации коэффициента параболического уравнения // Сиб. матем. журн., 1995. Т.36. № 1. C.177-185.
[7] Прилепко А.И., Костин А.Б. О некоторых обратных задачах для параболических уравнений с финальным и интегральным наблюдением // Матем. сб., 1992. Т. 183. № 4. C. 49-68.
[8] Yan L, Fu C.L, Yang F.L. The method of fundamental solutions for the inverse heat source problem // Eng. Anal. Boundary Elements. 2008. Vol.32. Pp. 216-222.
[9] Ismailov M.I, Kanca F, Lesnic D. Determination of a time-dependent heat source under nonlocal boundary and integral overdetermination conditions // Appl. Math. Comput., 2011. Vol. 218. Pp. 4138-4146.
[10] Cannon J.R, Duchateau P. Structural identification of an unknown source term in a heat equation // Inverse Probl., 1998. vol.14. pp. 535-551.
[11] Aida-zade K.R., Abdullayev V.M. Solution to a class of inverse problems for system of loaded ordinary differential equations with integral conditions // J. Inv. Ill Posed Problems. 2016. vol. 24. no. 5. pp.543-558.
[12] Abdullayev V. M. , Aida-zade K. R. Optimization of loading places and load response functions for stationary systems // Comput. Math. Math. Phys., 2017. vol. 57. no. 4. pp. 634-644.
[13] Abdullaev V. M. , Aida-zade K. R. On the numerical solution of loaded systems of ordinary differential equations // Comput. Math. Math. Phys., 2004. vol. 44. no. 9. pp.1505-1515.
[14] Aida-zade K.R. A Numerical method of restroring the parameters of a dynamic system // Cybern. Syst. Analysis., 2004. Vol. 40. no. 3. pp. 392-399.
Для цитирования: Абдуллаев В. М. Численное решение задачи параметрической индентификации для дифференциальных уравнений с частными производными // Вестник КРАУНЦ. Физ.-мат. науки. 2018. № 2(22). C. 33-44. DOI: 10.18454/2079-66412018-22-2-33-44
For citation: Abdullayev V. M. Numerical solution to parametric identification problems partial differential equations, Vestnik KRAUNC. Fiz.-mat. nauki. 2018, 22: 2, 33-44. DOI: 10.18454/2079-6641-2018-22-2-33-44
Поступила в редакцию / Original article submitted: 03.04.2018 В окончательном варианте / Revision submitted: 08.04.2018
Vestnik KRAUNC. Fiz.-Mat. Nauki. 2018. no.2(22). pp. 33-44. ISSN 2079-6641
DOI: 10.18454/2079-6641-2018-22-2-33-44 INFORMATION AND COMPUTER TECHNOLOGIES
MSC 65M32
NUMERICAL SOLUTION TO PARAMETRIC IDENTIFICATION PROBLEMS FOR PARTIAL DIFFERENTIAL EQUATIONS
V. M. Abdullayev12
1 Azerbaijan State Oil and Industry University, Azadlig avenue 20, AZ1010, Baku, Azerbaijan
2 Institute of Control Systems, Azerbaijan National Academy of Sciences, 9 B. Vahabzade str. Baku, AZ1141, Azerbaijan
E-mail: vaqif_ab@rambler.ru
The paper proposes a numerical method for solving the problem, based on the use of the method of lines to reduce the problem to a system of ordinary differential equations with unknown parameters. Next, we use a special representation of the solution of the obtained boundary value problem for a linear system of differential equations with nonlocal conditions, by means of which the problem of parametric identification reduces to solving auxiliary boundary value problems and a system of algebraic equations. It is important to note that, unlike optimization approaches, The construction of any iterative procedures or minimizing sequences is used.
Key words: inverse problem, method of lines, parabolic type equation, parametric identification.
© Abdullayev V.M., 2018