Несмищенная и смещённая оценки для решения задачи Коши для обобшенного уравнения неизотропной диффузии



Download 0,76 Mb.
bet1/3
Sana24.02.2022
Hajmi0,76 Mb.
#196351
  1   2   3
Bog'liq
1.2 параграф


1.2 Несмищенная и - смещённая оценки для решения задачи Коши для обобшенного уравнения неизотропной диффузии.
Рассмотрим задачу Коши в классической постоновке в мерном пространстве в слое
, (2.1)
, , (2.2)
где и постоянные коэффчциенты, .
Пусть - матрица, - матрица.
Предлогаем:

  • матрица симметрична и положительно определена.

  • матрица такая что, матрица Грамма положительно определена (невырождена).

Функции и будем считать непрерывными в . В педложении сушествования единственности решения задач (2.1)-(2.2) построим алгоритм для его численной реализации.
Для дальнейшего будут удобны следующие обозначения. Введем (в блочной записи) матрицы. Эти матрицы имееют следующий вид
, ,
, ,
блок размера , блок размера , блок размера , размера , матрицы аналогичных размеров. Здесь и далее - единичная матрица, - диагональная матрица .
Лемма 2.1. Матрицы и при положительно определены.
Доказательства. Легко проверить, что матрица в блочном разбиении имеет вид

Достаточно доказать положительную определенность матрицы . Соответствуюшая ей квадратичная форма от переменных ( - - мерное евклидово пространство строк) и имеет вид
Видим, что она равна сумме двух квадратичных форм с невырожденными матрицами Грамма и и поэтому . Равенство влечет равенства и , т.е. Лемма доказана.
Из доказанного следует, что матрица является симметричной и положительно определенной. Это означает что квадратичная форма для произвольного , поэтому матрица представима в виде произведения некоторой матрицы на ее транспонированную , а именно .
Тогда - фундаментальное решение уравнение (2.1) с сингулярностью в точке имеет вид
(2.3)
где , .
Спомошью фундаментального решения определим в области (шароиды) следующим образом. Пусть число (параметр) .
Область

будем называть шароидом радиуса с центром в точке . можно переписать в виде
(2.4)
Из (2.4) видно, что удовлетворяет следущим условиям , . Каждое сечение шароида горизонтальной плоскостью , , является - мерным эллипсоидом с центром в точке . Если . то . При и монотонно стягиваются к центру .
Поэтому сушествует такое что при , .
Пусть такое что тогда для решения задачи (2.1)-(2.2) справедливо следующие соотношение
, (2.5)
где
,
.
Здесь мерная единичная сфера с обычными ортогональными координатами , для , . единичный - мерный вектор.
, (2.6)
.
Из (2.6) следует, что если , при , то .
Приступим к построению марковский цепи , на которой будем строить несмещенную оценку решения задачи (2.1)-(2.2).
Для , преминев формулу (2.5) полочим

Отсюда следует, что ядро интегральрого ураинения (2.5) можно рассматривать как плотность распределения.
Рассмотрим интеграл
.
Сделав замену переменных , получим

где - гамма функция. Теперь можем представить следующим образом

где - плотность гамма-распределенной случайной величины с параметром , .
Алгоритм выбора матрицы , удовлетворяющей соотношению приводится в параграфе 1.1. Для наших целей будем моделировать случайной вектор с полностью распределения
, (2.7)
где - индикатор множества , - поверхность единичной сферы, - элементы матрицы размера - . Так как являются координатами единичного вектора и получим , где - наибольшее собственное значение матрицы . Тогда можно моделировать случайный вектор с плотностью распределения (2.7) методом Неймана. Приведем алгоритм моделирования.

Download 0,76 Mb.

Do'stlaringiz bilan baham:
  1   2   3




Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©hozir.org 2024
ma'muriyatiga murojaat qiling

kiriting | ro'yxatdan o'tish
    Bosh sahifa
юртда тантана
Боғда битган
Бугун юртда
Эшитганлар жилманглар
Эшитмадим деманглар
битган бодомлар
Yangiariq tumani
qitish marakazi
Raqamli texnologiyalar
ilishida muhokamadan
tasdiqqa tavsiya
tavsiya etilgan
iqtisodiyot kafedrasi
steiermarkischen landesregierung
asarlaringizni yuboring
o'zingizning asarlaringizni
Iltimos faqat
faqat o'zingizning
steierm rkischen
landesregierung fachabteilung
rkischen landesregierung
hamshira loyihasi
loyihasi mavsum
faolyatining oqibatlari
asosiy adabiyotlar
fakulteti ahborot
ahborot havfsizligi
havfsizligi kafedrasi
fanidan bo’yicha
fakulteti iqtisodiyot
boshqaruv fakulteti
chiqarishda boshqaruv
ishlab chiqarishda
iqtisodiyot fakultet
multiservis tarmoqlari
fanidan asosiy
Uzbek fanidan
mavzulari potok
asosidagi multiservis
'aliyyil a'ziym
billahil 'aliyyil
illaa billahil
quvvata illaa
falah' deganida
Kompyuter savodxonligi
bo’yicha mustaqil
'alal falah'
Hayya 'alal
'alas soloh
Hayya 'alas
mavsum boyicha


yuklab olish