Научная статья на тему 'Статистический анализ и обработка данных имитационного моделирования транспорта электронов в облученных материалах'

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

CC BY
49
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МОДЕЛЬ ТРАНСПОРТА ЭЛЕКТРОНОВ / МЕТОД МОНТЕ-КАРЛО / МЕТОД НАИМЕНЬШИХ КВАДРАТОВ / ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ / ПАРАМЕТРЫ АППРОКСИМАЦИИ / ОДНОФАКТОРНЫЙ ДИСПЕРСИОННЫЙ АНАЛИЗ

Аннотация научной статьи по математике, автор научной работы — Поправка С. Т., Павельчук А. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Статистический анализ и обработка данных имитационного моделирования транспорта электронов в облученных материалах»

2. Чучуева, И.А., Павлов, Ю.Н. Модель прогнозирования временных рядов по выборке максимального подобия: Дис. ...канд. техн. наук / МГТУ им. Н.Э. Баумана, 2012. - 153 с. - Режим доступа: http://www.mbureau.ru/ sites/default/flles/pdf/Chuchueva-Dissertation.pdf.

3. Временной ряд: Википедия. Свободная энциклопедия - Режим доступа: https://ru.wikipedia.org/wiki/ %00%92%01 %80%В0%В5%В0%ВС%В0%В5%В0%ВВ%В0%ВВ%В0%ВЕ%В0%В9_%Б 1 %80%В 1 %8Р%В0% В4

4. Метод экспоненциального сглаживания: 4апа1уйс8. Решения для бизнес-анализа. - Режим доступа: https://4analytics.ru/prognozirovanie/rn

5. Прогноз по методу экспоненциального сглаживания Хольта: 4апа1уйс8. Решения для бизнес-анализа. -Режим доступа: https://4analytics.ru/prognozirovanie/prognoz-po-metodu-eksponencialnogo-sglajivaniya-xolta.html

УДК 514.13

Т.А. Юрьева, А.П. Филимонова

АПРИОРНЫЕ ОЦЕНКИ РЕШЕНИЯ УРАВНЕНИЯ ТИПА МОНЖА - АМПЕРА

В МЕТРИКЕ Cm+2'a\S?)

В статье рассматривается доказательство существования априорных оценок решения дифференциального уравнения типа Монжа - Ампера в метрике Cm+2-a'(S?).

Ключевые слова: гиперболическое пространство, уравнение Монжа - Ампера, отрицательная эллиптичность, бельтрамиевы координаты.

В работе [1] мы ввели дифференциальное уравнение вида:

РпР22 ~ Pl2 ~ Pi 1 (20thP ' Pv + SHP ' СНР) + 2PnPuPvCtkP ~ Pll (20thP ' Pi + SHP ' CHP C0S2 V) "

2

~(pl cos2 v + )2 + 2pl+ 2pi cos2 v + sh2pcos2 v = cosv

- Kt{u,v,p) + Pi v + sh2p ■ cos2 v)2 _

COS V

здесь pn, pl2, p21 - вторые ковариантные производные функции p = p(u,v) относительно метрики единичной сферы Sf .

Напомним, что к данному уравнению приводит задача восстановления регулярной выпуклой гомеоморфной сфере (О - центр сферы, радиус равен 1), звездной относительно точки О поверхности F: р = р(и, v) в трехмерном пространстве Нъ постоянной отрицательной кривизны (пространство Лобачевского) по ее функции гауссовой кривизны. ¿Sf рассматривается как двумерное многообразие, атлас на ¿Sf выбран так, что в каждой его карте выполняется неравенство cosv>a>0; и, v — сферические координаты. Уравнение отрицательно эллиптично при условии Kint(u,v,p) = Ki(u,v,p) = Ki>-1 [1].

Исследование дифференциальных уравнений такого типа начинается с получения априорных оценок решения в соответствующей метрике. В работе [2] мы получили оценку решения р — p(u,v)

уравнения в метрике С°(5'12), то есть оценку самого решения.

Лемма 1. Пусть в Нъ зафиксированы две концентрические сферы с центрами в точке О и радиусами рх и р2 (рх< р2 )• Пусть функция Кы(и9у9р) = К1 определена в £2х7?+ и удовлетворяет

условиям: 1) Кг > -1; 2) Кг = —\— + 1а{и, V,/?), где к > О внутри сферы 52 и к < 0 - вне сферы 52 .

р

В этом случае всякое решение р = р(и, V) исследуемого уравнения задает поверхность Р:р = р(и,у)9 расположенную между сферами и 82 . Иными словами, при перечисленных выше условиях, наложенных на функцию Кш{и9у9р)-К^ есть априорная оценка решения р{и9у) уравнения рх < р(и,у) < р2.

Оценки решения р = р(и, V) исследуемого уравнения в метрике С1 (£2) нами получены в работе [3].

Лемма 2. Пусть К[п1(и9у9р) удовлетворяет условиям леммы 1, а /? = р(м,у)еС2(Х12) - решение исследуемого уравнения. Тогда первые производные функции р - р{и9у) ограничены некоторой

о о

постоянной, зависящей только от чисел рх и р2: |/?м(и,у)|<-|д,(г/,у)|<--.

В работе [4] показан процесс получения априорных оценок решения р - р{и, V) исследуемого уравнения в метрике С2(£2). Результат сформулируем в виде леммы 3.

Лемма 3. Пусть функция Кы(и9у9р) удовлетворяет условиям леммы 1 и принадлежит классу С2 (£2 х ) . Тогда любое решение р = р(и,у) исследуемого уравнения, принадлежащее классу С2(£2), ограничено в метрике С2(£2) постоянной, которая зависит от чисел рх, р2 и свойств функции Кы(и,г,р) = Кп а именно: М Кг{и9у9р)9 \Кг{и9у9р)\ 2 2 г

Теперь приведем доказательство существования априорных оценок решения р = р(и,у) исследуемого уравнения в метрике С/77+2'а (£2). Справедливо следующее утверждение.

Пусть функция Кы(и9у9р) = Кг принадлежит классу Ст,а($2 х7?+), т> 2, 0 < а < 1. Далее пусть функция Кы(и9у9р) = К1 удовлетворяет условиям леммы 1. Тогда решение р = р(и9у) исследуемого уравнения ограничено в метрике ст+2>а )? где а' <а, причем постоянная к= Ы Кг(г/,V,р) и норма функции Кг(щу9р): .

Из условий сформулированного утверждения, наложенных на р{и9у) и функцию Кг(и9у9р)9 и результата работы [4] следует, что решение г = г(х,у) уравнения:

2 1Ч[(1 + Р2 +Ч2)-(рх + ЧУ~г)2]2 ~ г

П-ъ =(л^+1)——--^——^-^—--, полученного в той же работе, ограничено в метрике

(1 х у % )

С2(К2). Напомним, что мы пользовались моделью Кэли - Клейна пространства Нъ, х, у, г - бель-трамиевые координаты, а круг К2 в плоскости хОу введен в указанной выше работе. Тогда из работы [5] решение 1 = 1(х9у) уравнения в бельтрамиевых координатах в предположении принадлежности функции К1(и9у9р) = К1 классу Ст,а($2 хЯ+) будет ограничено в круге К2 и в метрике ст+2,а ,

а'<а. Норма <С', где постоянная С зависит от чисел k= inf Кг{и,у,р) и

SlÁPuPi]

\\KXu9v9p)\\cM^2^+y

В работе [4] введены функции й = й(г/, v), v =v(u,v). Эти функции аналитические. Из равенств [4]:

- dv

COS V--= sin v0 cos v sin(w - Щ ) ,

du

- dv

cos v--= cos v0 cos v + sin v0 sin V cos(u - UQ ) ,

dv

- du . - . - dv

cosvcos и---sin v sin w--- cosvcos(w-w0),

du du

- du . - . - dv

cos v cos и---sin v sin и--= - sin v sm(w - щ),

dv dv

- d2v

Гdv^

cosv

du2 - d2v

smv-

\duJ

= sin v0 cos v sin(w - u0) ,

. - dv dv

cos v---sinv----= - sin v0 sin v sm(w - щ),

dudv dv du

- d2v . -

COSV--r-SinV'

dv2

Гdv^

= - cos v0 sin v + sin v0 cos v cos(w - u0) ,

- d и .- - dv du -.-

COS VCOS и--r-sinvcosw----COSVSinW'

du du du

fdu^

\duJ

- d2v -sinvsinw--

du2

-cosvsinw-

Гdv^

KdUJ

. - - dv du

-smvcosw---- -cosvsm(w-w0),

du du

- d и .- - dv du - . - du du . - . - d v

cos veos и---smvcosw----cosv sin г/----sin v sin w---

dudv dv du du dv dvdu

- . - dv du . - - du dv

-cosvsmw----smvcosw---= -smvcos(w-w0),

dv du dv du

- d и .- - dv du ._ _ cosvcosw--- - sin vcosw----sin и cosv

dvz

dv dv

/du\ . - . - <92v -sinvsmw--

dv¿

-COSVSinM

. - - dv du

- sin v cos и---- - cos v sm(w - щ),

dv dv

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

множителями функции cosv и eos veos^ , близкие к 1 и (-1) соответственно, в окрестности точки (л\0).

Тогда из соотношений [4]:

du dv _ d2u

Рии Рии

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

fdu^

\dUJ

+ 2 рш

ди ди ди

Гdv^

уди;

. a2v

du du du dv d2u

Puv Рш л л Pñv л л Рй ^ ^

ди dv ди dv dvdu

_ du dv _ dv dv _ <32v

" Pvñ Z ^ Pvv ~Z Г ^ Pv : dv du dv du dudv

Pvv = Рй

/ди"

vövy

_ _ dv du _ д2 и _ ov ov ov

í dv^ vÖvy

d2v

последовательным дифференцированием получаем равенства, из которых следует ограниченность функции р(и,у) в метрике ст+2,а в некоторой окрестности при условии ограниченности функции р(й9 V) в соответствующей при движении В окрестности.

Напомним, что р = р(й,у) и В введены в работах [3,4]. Движение В (модель Кэли - Клейна пространства Нъ) задано матрицей:

В =

COS V, 0

о

- sm vn sm щ - sm vn cos г/(

cos и.

о

о"

- sin и{

о

- sm V(

о

cosv0 sm и{

о

cosv0 cos г/,

о у

В работе [3] мы получили равенства:

р-(eos veos и -zx -sinv-z^ cos v sin и) = shpchpcosv(zx sin u + zy cos и),

/};(cosvcosw -zx -sinv-z^ cosvsinw) = shpchp(smvcosu + zx eosv-zy sinvsinw), а в работе [4]

путем дифференцирования этих равенств по й и v вторично получили равенства, приводить которые вследствие их объема в данной работе не будем. С помощью последовательного дифференцирования последних приходим к выводу, что из ограниченности функции z = z(x,y) в метрике ст+2,а в круге К2 следует ограниченность функции p(u,v) в метрике ст+2,а в окрестности со, так как при старших производных функции p(u,v) множителем всегда является выражение: eos veos и — zx -sinv-z^ cos v sin и . Это выражение в указанной выше окрестности близко к (-1).

Из приведенных нами рассуждений следует наличие «локальной» оценки решения р - р{и, v)

исследуемого уравнения в метрике ст+2,а . В силу того, что сфера S2 является компактным многообразием, из наличия «локальных» оценок решения уравнения следует ограниченность решения p = p(u,v) в метрике Cm+2'a\S?).

Постоянная ограничения зависит от чисел pl9 р2, k = inf Кг (г/, v, р) и нормы функции

Kt(u,v,p): IK¿u9v9p)\\c4s}xW^r

Утверждение, сформулированное нами в данной работе, доказано.

1. Филимонова, А.П., Юрьева, Т.А. Аналог теорем расположения замкнутых выпуклых поверхностей с заданной функцией внутренней кривизны в пространствах постоянной кривизны // Вестник АмГУ. - 2017. -Вып. 79.- С. 17-21.

2. Филимонова, А.П., Юрьева, Т.А. Априорные оценки решения в метрике С0^2) уравнения типа Монжа - Ампера на сфере как двумерном многообразии в пространстве постоянной кривизны // Международный научно-исследовательский журнал. - 2016. - № 9-2(51). - С. 132-136.

3. Филимонова, А.П., Юрьева, Т.А. Априорные оценки градиента решения уравнения некоторого класса Монжа - Ампера // Вестник Бурятского гос. ун-та. - Серия «Математика, информатика». - 2019. - № 1. -С. 49-55.

4. Филимонова, А.П., Юрьева, Т.А. Априорные оценки решения некоторого дифференциального уравнения типа Монжа - Ампера в метрике C2(S2) на сфере S2 как двумерном многообразии//Вестник Амурского гос. ун-та. - 2019. -№ 85. - С. 9-15.

5. Nirenberg, L. On nonlinear elliptic partial differential equations and H"older continuity // Comm. Pure Appl. Math. - 1953.-№ 1.

i Надоели баннеры? Вы всегда можете отключить рекламу.