_ ПРОБЛЕМЫ ЯДЕРНОЙ, РАДИАЦИОННОЙ _
И ЭКОЛОГИЧЕСКОЙ БЕЗОПАСНОСТИ -
УДК 51-73, 517.965
ОДНОЗНАЧНОСТЬ РЕШЕНИЙ ПРИ ИСПОЛЬЗОВАНИИ ЛИНЕЙНОГО ФУНКЦИОНАЛЬНОГО УРАВНЕНИЯ В МОДЕЛИ
РАДИАЦИОННОЙ ЗАЩИТЫ
© 2019 В.П. Чернявский
Саровский физико-технический институт - филиал Национального исследовательского ядерного университета МИФИ, Саров, Нижегородская обл., Россия
В статье рассматривается функциональное линейное уравнение со сдвигом в модели радиационной защиты для процессов переноса заряженных частиц и ионизирующего излучения. Целью работы является изучение вопросов существования и однозначности решений для различных случаев, которые возникают при изменении начальных параметров модели. Анализ системы, сопутствующей функциональному уравнению, проводится методами линейной алгебры. Для случая неравенства нулю главного определителя сопутствующей системы показана корректность полученных решений; функциональное уравнение имеет единственное решение. В случае равенства нулю определителя задача полностью решена для циклов длины 2. Функциональное уравнение не имеет решений, если определитель равен нулю, а ранг расширенной матрицы равен 2. Для случая совместной системы с вырожденной матрицей получены аналитические формулы общего решения однородного и неоднородного функциональных уравнений. Эти решения зависят от коэффициентов исходного уравнения, начальной функции, порождающей цикл, и содержат произвольно выбранную функцию. Для устранения возникающей неоднозначности можно перейти к модели с невырожденной матрицей, изменив систему весовых коэффициентов модельного уравнения, или привлечь дополнительные начальные условия.
Ключевые слова: линейное функциональное уравнение, уравнение со сдвигом, итерация, цикл, сопутствующая система, ранг, однородное и неоднородное функциональное уравнение.
Поступила в редакцию 16.09.2019 После доработки 05.11.2019 Принята к публикации 06.11.2019
Введение
Методы математического моделирования применяются при изучении различных вопросов радиационной безопасности [1-3]. В одной из моделей используется линейное функциональное уравнение со сдвигом, разработанное Карлеманом [4]. Такая модель описания впоследствии развивалась [5, 6] и применялась во многих задачах, в том числе, при решении уравнений гиперболического типа [7], к которым относятся уравнения переноса заряженных частиц и ионизирующего излучения в различных материальных средах. В настоящей работе рассматриваются вопросы однозначности решений линейного функционального уравнения с циклическими итерациями аргумента (уравнения со сдвигом) [8]:
а0(х)и(х) + а!(»(/О)) + а2(х>(/(/(х))) + ... + ап_1 (х)м(/(/(.../(х)...)))= Ъ(х), (1),
я-1
где ай{х),а1{х),...,ап_1{х),Ь{х\/{х) - заданные функции, определённые на множестве Ха Я, причём У а X, где У =/(X) - множество значений Дх); и(х) - искомая
© Национальный исследовательский ядерный университет «МИФИ», 2019
функция. Предполагается, что функция fx) имеет обратную, а её последовательные итерации f (х), f (f (х)), f (f (f (х))),... образуют цикл длины n (условие Карлемана), т.е.
/(/(.../(х)...)) ф X при \<к<п и /(/(.../(х)...)) = X.
(2)
Введём обозначения: /о(X) = х, / (х) = /(х), Л(х) = ./(./(.,../(х)...)). к = 1,2,..., п -1,
следовательно,
/и (х) = х = /0 (х) - условие цикличности. Обозначим также
аЫ (х) = а, (Л (*))> мк (х) = "(Л (*))> Ьк (х) = К/к (*)), /',Л: = 1,2,..и -1.
Тогда уравнение (1) принимает следующий вид:
а0.0 (Ж (х) + а0.1 (х)^ (х) + а 0.2 (х)и2 (х) +... + а0. п _ 1 (х)ип _1(х) = Ь0 (х). (3)
Подставляя последовательно в уравнение (3) вместо х итерации /к(х) (к = 1,2,...,и -1), вместе с уравнением (3) получим сопутствующую ему систему равенств:
а 0.0 (х)и0 (х) + а0.1 (хХ (х) + а0.2 (х)и2 (х) +... + а 0. п _ х (х)ип _х(х) = Ь0 (х)
а 1; п - 1 (х)и0 (Х) + а 1; 0 (ХХ (Х) + а 1; 1 (ХК (Х) + ■■■+ а 1; п - 2 (ХК - 1 (Х) = Ь1 (Х)
............................................... (4)
а п- 2; 2 (х)и 0 (х) + а ^ 2; 3 (х)М1 (х) + а п - 2; 4 (х)и 2 (х) + ■■■ + а „- 2; 1 (ХК - 1 (х) = 2 (х) ап- 1; 1 (ХК( ХН 1; 2 ( Х)М1( Х) + 1; 3 ( Х)М 2( х) + ■■■ + а „- 1;0( ХК - 1( Х) = Ъп- 1( Х)
Данная система уравнений является линейной относительно величин и0 (х), их (х),..., ип_х (х). Пусть А(х) - главный определитель системы:
А =
a
0; 0
a
0; 1
a
1;и-1
a 1; 0
ап-2;2 ап-2;3
an-hl an-V2
a
0;п-2
a
0;и-1
a
1;и-3
а1\п-2
a
и-2; 0
a
п - 1; п -1
a
a
и-2:1
и-1:0
(5)
Частного вида уравнения для случая А Ф 0 периодически встречались в печати (см. [9, с.82], [10, с.12], [11], [12, с.413] [13, с.470]). Однако условие А ф 0 не является необходимым для разрешимости исходного уравнения и, как отмечено [8, с. 58], требуются дополнительные исследования. В настоящей работе для А Ф 0 показана непротиворечивость формул, составляющих решение системы (4). Исследование вопросов однозначности решений при условии А = 0 полностью выполнено для циклов с п = 2.
Линейное уравнение для случая А Ф 0
При условии А Ф 0 решение системы, и притом единственное, находится по формулам Крамера [14, 15]:
k
n
k
мк(х) =
Д(х)
к = ОД,...,« —1,
(6)
где Дк (х) - определитель, полученный из (5) заменой к + 1-го столбца на столбец свободных членов. Решение уравнения (1) следует из (6) при к = 0: и(х) = и0(х) = До(х)/Д(х), но т.к. функции и к (х) получаются последовательными итерациями: ик + Лх) = ик(/(х)\ т0 формулы (6) в указанном смысле не должны противоречить друг другу.
Лемма 1
Пусть Д Ф 0, тогда решения (6) системы (4) корректны, т.е. удовлетворяют условиям
Д*+1(*) А* (/(*))
Д(х) Д(/(х))
Лг = 0,1.....и-1.
Доказательство. Подставим в (5) /(х) вместо х:
А(/(*)) =
а
1; 0
а
а
а
2; и-1
а
1;1
2; 0
и-1; 2 а 0; 1
а
и-1:3
а 0; 2
а
1; и - 2
а
1; и-1
а2\п-Ъ а2:п-2
ап-\\0 Ди-1;1 Д0;и-1 а0\0
Если последнюю строку сделать первой и последний столбец также сделать первым, то полученный в результате перестановок определитель совпадёт с определителем (5). Таким образом, Д(/(х)) = (-1)" • (-1)" Д(х) = Д(х). Точно так же доказывается, что если в определителе Ак(/(х)) последнюю строку сделать первой, а последний столбец также сделать первым, то полученный определитель совпадёт с определителем Д^+1(х), значит ДД/(х)) = (-1)" •(-1)иДА+1(х) = к = \,2,...,п-2
(при к = п - 1 различие лишь в том, что Дп_1(/(х)) = (-1)" -(-1)" Д0(х) = Д0(дс)). Лемма доказана.
Таким образом, если Д Ф 0, то уравнение (1) имеет единственное решение, определяемое формулой (7):
и{х) =
Д0(х)
Д(х)
(7)
Пример 1
Решить уравнение:
и{х) + 2и(/(х)) + МКт)) + ■■■ + пи{/{/{.. ./(х)...))) = ад, (8)
ч_ _;
и-1
где /х) - произвольная заданная функция, порождающая цикл (2) длины п. В обозначениях (1) ак (х) = к + 1, к = 0,1, ..., п - 1; главный определитель Д системы (4) является циркулянтом [16] и отличен от нуля [17]:
А =
2 1
4 3
п -1 п п-2 п-1
1 и
2 1
(-1)
и-1й„-1
(и + 1)
Для нахождения решения (7) определитель
А0 =
Ьо 2 3 . .. п-1 п
Ь1 1 2 . .. п-2 п-1
Ь2 п 1. п- 3 п-2
2 4 5 . . . 1 2
3 4. . . п 1
разложим по первому столбцу: где А
к.0 алгебраическое дополнение элемента Ьк, к — 0,1,...,п — 1. Непосредственный подсчёт даёт:
_(-!)" п"5 (п1 +п-2) _(-\)"п"-3(п2 +п + 2)
Аоо „ > Аю — .
2
2
Вычисления остальных алгебраических дополнений проводятся по единой схеме и приводят к одному и тому же результату:
Ак0=(-\у-1п"-\ к = 2,Ъ,...,п-\.
В итоге получаем единственное решение уравнения (8):
и{х) =
2
( п2+п-2
п2(п +1)
60(х) +
п +п + 2 2
Ъ1(х) + Ъ2(х) + ... + Ъп_1(х)
где Ъ0(х) = Ъ(х\ Ък(х) = ¿(/(.../(х)...)), к = 1, 2, ...,п — \. Отметим также, что уравнение
' к '
(8) с системой весовых коэффициентов {1, 2, ..., п} равносильно уравнению с системой коэффициентов {1, 1/2, ..., 1/п} (для установления равносильности достаточно обе части поделить на п и заменить последнее слагаемое и($^(...Ах).))) = Чх)).
Линейное уравнение для случая А = 0, п = 2
Проанализируем особенности нахождения решения задачи с А = 0 при п = 2. Уравнение (1) при п = 2 имеет следующий вид:
а0 (х)и(х) + ах (х)г/(/(х)) = Ь(х),
2
2
где функция Дх) инволютивна, т.е. порождает цикл длины 2: /(/(х)) = х, а каждая из функций а0 (х), а (х) тождественно отлична от нуля (иначе уравнение (9) становится тривиальным). Для сопутствующей уравнению системы
Г а0(х)и(х) + а1(х)и(/(х)) = Ь(х)
\ а1(/(х))и(х) + а0(/(х))и(х) = Ь(/(х)) (10)
главный определитель равен 0, значит, его строки пропорциональны с некоторым коэффициентом к(х): а0(х) - к(х)а] (/(л;)), ал (л;) = к(х)а0 (/(л;)). Если Ь(х) Ф к(х)Ь(/(х)),
тогда ранг расширенной матрицы
^ а (х) а (х) ь(х) ^
равен двум и система
Л(/(х)) а0(/(х)) Ь(/(х))/
(10) несовместна, следовательно, уравнение (9) не имеет решений. Пусть теперь Ь(х) = к(х)Ь(/(х)), т.е.:
а0(х) _ а,{х) _ Ь{х) =
«,(/(*)) Яо(Дх)) *(/(*)) ' л (П>
Ранг равен 1 и от системы (10) остаётся только уравнение (9). Метод решения, используемый в п.1, не применим. Введём функционал:
Ц_у(х)\ = а0 (х) у(х) + (х)у(/(х)\ (12)
который, очевидно, является линейным:
Ь[аи(х) + /?у(х)] = аЦи(х)~\ + /?Цу(х)] (V«,р е Л). Запишем с помощью (12) неоднородное уравнение (9):
Дк(*)] = 6(х), хеХ, (13)
и соответствующее ему однородное уравнение:
Ци(х)] = 0, хеХ (14)
Определение
Множество функций £2 о называется общим решением уравнения (14) наХ, если:
1) любая функция у(х), хеХ., принадлежащая £2 о, является решением уравнения
(14);
2) любое решение у(х), хеХ уравнения (14) принадлежит £2 о-Аналогично определяется общее решение неоднородного уравнения (13).
Лемма 2
Общее решение однородного уравнения (14) определяется формулой (15):
О) = «о (Дх)Мх) - аг ОМ Д *)), (15)
или равносильной ей:
Щ = (16)
2 2 а0 (л;)
где ^(х), у(х) - произвольные на X функции. Доказательство.
Будем искать решение уравнения (14) в виде:
и0 (х) = А(х)м>(х) + В(х)м>(/(х)),
где ^(х) - произвольная функция, а функции А(х) и В(х) подлежат определению. Подставим и0(х) в (14) и перегруппируем слагаемые:
к (х)А(х) + а1 (х)В(Ях))Мх) + (а0 (х)В(х) + а, (х) Д/(х))М/(х)) = О
Т.к. функция ^(х) произвольна, то отсюда следуют условия:
| а0 (х)А(х) + аг (х)В(/(х)) = О [ао( х)В( х) + аД х) А(/ ( х)) = 0,
которые можно записать в следующем виде:
А(х) _ В(х) _ аДх) В(/(х))= А(/(х))=~^)'
Величины А(х) и В(х) определяются неоднозначно. Если взять А(/(х)) = а0(х), тогда А(х) - а0 (/(х)), В(х) - -а] (х) и получаем формулу (15). Если же взять А(/(х)) = 1, тогда А(х) = 1, В(х) = -а,(х)/а0(х), что соответствует равенству (16), и т.п. Но это будут равносильные представления одной и той же формулы. Так, равносильность (15) и (16) становится очевидной, если связать у(х) и м>(х) соотношением: \'(х) = 2а0 (/(х))м>(х). Тем самым, для любых у(х) и ^(х) (15) и (16) являются решениями уравнения (14).
Осталось показать, что любое решение уравнения (14) содержится, например, в (16). Пусть g(x) - произвольная функция, удовлетворяющая (14), тогда - ах ОО^С/ОО) / а0 (х) - §(х). Теперь, чтобы убедиться, что g(x) содержится в (16), достаточно в качестве у(х) взять g(x). Лемма доказана.
Лемма 3
Пусть а0(х) Ф 0, тогда общее решение неоднородного уравнения (13) находится по формуле
«(*)=+Ы*) (17)
2а0(х) 2 2а0(х)
где у(х) - произвольная на X функция. Доказательство.
С помощью (11) легко убедиться, что функция ин (х) = Ь(х)/(2а0(х)) является частным решением неоднородного уравнения (13). Добавляя к ин (х) общее решение (16) однородного уравнения, получаем: и(х) = ин (х) + и0(х) - решение уравнения (13) для любой функции у(х), входящей в и0(х), т.к. в силу линейности функционала (12):
Ци„ (х) + Щ (*)] = Цин (х)] + Ци0 (х)] = ад.
Осталось доказать, что (17) содержит любое решение уравнения (13). Пусть g(x) произвольная функция, удовлетворяющая (13):
Ъ{х) = а0 (х^(х) + а1 (х^(Дх)).
(18)
Рассматривая теперь уравнение (13) с правой частью Ь(х), определяемой равенством (18), и подставляя (18) в (17), получим:
ф) = т + _а^л){т)
2
2а0 (х)
2
2 а0 (х)
Очевидно, что решение g(x) содержится в и(х), если в качестве у(х) взять g(x). Лемма доказана.
Пример 2
Найти все функции, удовлетворяющие уравнению:
х
х-1
и(х) + XII
х
х-1
х 2 - 2х + 2 х-\ ''
Х = Я\{ 1}.
(19)
х х~ — 2х + 2 х
Здесь /(х) =- и /(/(х)) = х - цикл длины 2; Ь(х) =-,ап(х) =-,а\
х-1 х-1 х-1
(х) = х. Для коэффициентов уравнения выполняются условия (11) (к(х) = 1), но А = 0 и
формула (7) не применима. Множество всех решений уравнения (19) находим по
формуле (17):
х2 - 2х + 2 1 ... х -1 ( х 4
2х 2 2
х-1
, (х ф 0, х ф 1),
V
где у(х) - произвольная функция.
Таким образом, в рамках примененной модели рассмотрен вопрос существования и единственности получаемых решений. Показано, что решения совместной сопутствуюшей системы в случае вырожденной матрицы уже при п = 2 содержат неоднозначность - произвольно выбранную функцию. Для устранения возникающей неоднозначности можно перейти к задаче с невырожденной матрицей, изменив систему весовых коэффициентов модельного уравнения, или привлечь дополнительные начальные условия.
СПИСОК ЛИТЕРАТУРЫ
1. Кучин, Н. Л. Математическое моделирование радиационного воздействия атомных объектов морской техники на окружающую среду и человека : диссертация доктора физико-математических наук / Н. Л. Кучин. - Санкт-Петербург, 2002. - 297 с.
2. Чирская, Н. П. Математическое моделирование свойств неоднородных структур для систем радиационной защиты / Н. П. Чирская, Е. Н. Воронина, В. Н. Милеев, Л. С. Новиков, В. В. Синолиц // Труды XXI Международной конференции «Радиационная физика твердого тела», т. 2. - Москва : ГНУ НИИ ПМТ, 2011. - С. 436-443.
3. Крюк, Ю. Е. Математические методы моделирования в оптимизации радиационной защиты / Ю. Е. Крюк, И. Е. Кунец // Вестник НТУ ХПИ. Серия: Информатика и моделирование. -Харьков : НТУ ХПИ. - 2011. - № 36. - С. 95-100.
4. Carleman, Т. Sur la theorie des equations integrates et ses applications, Verhandl. Internat. Math. Kongr. Zurich, 1 (1932), P.138 - 151.
5. Литвинчук, Г. С. Краевые задачи и сингулярные уравнения со сдвигом / Г. С. Литвинчук. -Москва : Наука, 1977. - 448 с.
6. Карапетянц, Н. К. Уравнения с инволютивными операторами и их приложения / Н. К. Карапетянц, С. Г. Самко. - Ростов-на-Дону : Издательство Ростовского университета, 1988. - 187 с.
7. Василевский, Н. Л. Об одном классе сингулярных интегральных уравнений с инволюцией и его применениях в теории краевых задач для дифференциальных уравнений в частных производных. II / Н. Л. Василевский, А. А. Карелин, П. В. Керекеша, Г. С. Литвинчук // Дифференциальные уравнения. - 1977. 13:11. - С. 2051-2062.
8. Антоневич, А. Б. Линейные функциональные уравнения: операторный подход / А. Б. Антоневич. - Минск : Издательство «Университетское», 1988. - 232 с.
9. Московский государственный университет. Справочник 2000. - Москва : Издательство Московского университета, 2000. - 240 с.
10. Агаханов, Н. X. Всероссийские олимпиады школьников по математике 1993-2006 / Н. Х. Агаханов, И. И. Богданов, П. А. Кожевников, О. К. Подлипский, Д. А. Терешин. -Москва : МЦНМО, 2007. - 472 с.
11. Бродский, Я. С. Функциональные уравнения / Я. С. Бродский, А. К. Слипенко. - Киев : Вища школа, 1983. - 96 с.
12. Полянин, А. Д. Справочник по интегральным уравнениям: Точные решения / А. Д. Полянин, А. В. Манжиров. - Москва : Издательство «Факториал», 1998. - 432 с.
13. Прасолов, В. В. Задачи по алгебре, арифметике и анализу / Прасолов В. В. - Москва : МЦНМО, 2007. - 608 с.
14. Мальцев, А. И. Основы линейной алгебры /А. И. Мальцев. - Москва : Наука, 2005. - 470 с.
15. Курош, А. Г. Курс высшей алгебры / А. Г. Курош. - Санкт-Петербург : Лань, 2019. - 432 с.
16. Фаддеев, Д. К. Лекции по алгебре / Д. К. Фаддеев. - Санкт-Петербург : Лань, 2002. - 416 с.
17. Фаддеев, Д. К. Сборник задач по высшей алгебре / Д. К. Фаддеев, И. С. Соминский. -Москва : Наука, 1977. - 288 с.
REFERENCES
[1] Kuchin N.L. Matematicheskoe modelirovanie radiatzionnogo vozdeistvia atomnikh objectov morskoi tecniki na okrudjaiushuiu sredu [Mathematical Modeling of the Radiation Effect of Atomic Objects of Marine Technology on the Environment and Humans]. Dissertaciya doktora fiziko-matematicheskix nauk [Dissertation of Doctor of Physical and Mathematical Sciences]. Sankt-Peterburg [Saint-Petersburg]. 2002. 297 p. (in Russian).
[2] Chirskaya N.P., Voronina E.N., Mileev V.N., Novikov L.S., Sinolits V.V. Matematicheskoe modelirovanie svoistv neodnorodnikh struktur dlia system radiatzionnoi zashiti [Mathematical Modeling of the Properties of Heterogeneous Structures for Radiation Protection Systems]. Moskva. Trudi XXI Mezdunarodnoi konferentsii «Radiatsionnaia phyzika tverdogo tela» [Moscow. Proceedings of the XXI International Conference «Radiation Solid State Physics»]. 2011. v. 2. P. 436-443 (in Russian).
[3] Kriuk J.E., Kunets I.E. Matematicheskie metodi modelirovania v optimizatsii radiatzionnoi zashiti [Mathematical Modeling Methods in Optimization of Radiation Protection]. Kharkov. Vestnik NTU KhPI. Seria: Informatika i modelirovanie [Bulletin of the National Technical University Kharkov Polytechnic Institute. Series: Computer Science and Modeling]. 2011. № 36. P. 95-100 (in Russian).
[4] Сarleman Т. Sur la theorie des equations integrates et ses applications, Verhandl. Internat. Math. Kongr. Zurich. 1 (1932). P.138-151 (in French).
[5] Litvinchuk G.S. Kraevie zadachi i singuliarnie uravnenia so sdvigom [Boundary Value Problems And Singular Equations with Shift]. Moscow. Pub. Nauka [Publishing House Science]. 1977. 448 p. (in Russian).
[6] Karapetiants N.K., Samko S.G. Uravnenia s involiutivnimi operatorami i ikh prilodjenia [Equations with Involutive Operators and their Applications]. Postov-na-Donu. Izdatelstvo Rostovskogo universiteta [Rostov on Don. Rostov University Publishing House]. 1988. 187 p. (in Russian).
[7] Vasilevsky N.L., Karelin A.A., Kereksha P.V., Litvinchuk G.S. Ob odnom klasse singuliarnikh uravneni' i ego primeneniakh v teorii kraevikh zadach dlia differentsialnikh uravneni' v chastnikh proizvodnikh. II [A Class of Singular Integral Equations with Involution and its Applications in the
Theory of Boundary Value Problems for Partial Differential Equations. II]. Differencial ny e uravneniya [Differential equation]. 1977, 13:11. P. 2051-2062 (in Russian).
[8] Antonevich A.B. Lineinie functional'nie uravnenia: operatorni podhod [Linear Functional Equations: Operator Approach]. Minsk. Izdatelstvo Universitetskoe [Publishing House University]. 1988. 232 p. (in Russian).
[9] Moskovsky Gosudarstvenny Universitet. Spravochnik 2000 [Moscow State University. Handbook 2000]. Moskva. Izdatelstvo Moskovskogo universiteta [Moscow. Moscow University Press]. 2000. 240 p. (in Russian).
[10] Agakhanov N.H., Bogdanov I.I., Kodjevnicov P.A., Podlipsky O.K., Tereshin D.A. Vserossi'skie olympiadi shkolnikov po matematike 1993-2006 [All-Russian Mathematics Olympiads 1993-2006]. Moscow. MCCME. 2007. 472 p. (in Russian).
[11] Brodsky J.S., Slipenko A.K. Funktsional'nie uravnenia [Functional Equations]. Kiev. Pub. Visha shkola [High school]. 1983. 96 p. (in Russian).
[12] Poliayin A.D., Manzhirov A. V. Spravochnik po integral'nim uravneniam: Tochnie reshenia [Handbook of Integral Equations: Exact Solutions]. Moscow. Pub. Factorial. 1998. 432 p. (in Russian).
[13] Prasolov V.V. Zadachi po algebre, arifmetike i analizu [Problems in Algebra, Arithmetic and Analysis]. Moscow. MCCME. 2007. 608 p. (in Russian).
[14] Mal'tsev A.I. Osnovi lineinoi algebra [The Basics of Linear Algebra]. Moscow. Pub. Nauka [Publishing House Science]. 2005. 470 p. (in Russian).
[15] Kurosh A. G. Kurs vishei algebri [Course in Higher Algebra]. Sankt-Peterburg [Saint-Petersburg]. Pub. Lan' [Publishing House «Doe»]. 2019. 432 p. (in Russian).
[16] Faddeev D.K. Lektsii po algebre [Lectures on Algebra]. Sankt-Peterburg [Saint-Petersburg]. Pub. Lan' [Publishing House "Doe"]. 2002. 416 p. (in Russian).
[17] Faddeev D.K., Sominsky I.S. Sbornik zadach po vishei algebre [Collection of Problems in Higher Algebra]. Moscow. Pub. Nauka [Publishing House Science]. 1977. 288 p. (in Russian).
Unambiguity of Decisions When Using Linear Functional Equation in the
Radiation Protection Model
V.P. Cherniavsky
Sarov Physical Technical Institute the branch of National Research Nuclear University "MEPhI", Dukhov St., 6, Sarov, Nizhny Novgorod region, Russia 607186 ORCID iD: 0000-0002-1628-2271 e-mail: [email protected]
Abstract - The paper considers a functional linear equation with a shift in the radiation protection model for the transport of charged particles and ionizing radiation. The aim of the work is to study the questions of the existence and uniqueness of solutions for various cases that arise when the initial parameters of the model change. The analysis of the system accompanying the functional equation is carried out by linear algebra methods. For the case of the inequality to zero of the main determinant of the accompanying system, the correctness of the obtained solutions is shown; the functional equation has a unique solution. If the determinant vanishes, the problem is completely solved for cycles of length 2. The functional equation has no solutions if the determinant is zero and the rank of the extended matrix is 2. For the case of a joint system with a degenerate matrix, analytical formulas for the general solution of a homogeneous and inhomogeneous functional equation are obtained . These solutions depend on the coefficients of the initial equation, the initial function generating the cycle, and contain an arbitrarily chosen function. To eliminate the ambiguity that arises, one can go to a model with a non-degenerate matrix by changing the system of weight coefficients of the model equation, or use additional initial conditions.
Keywords: linear functional equation, shift equation, iteration, cycle, associated system, rank, homogeneous and inhomogeneous functional equation.