Уточнение метода кригинга для исследования геоповерхностей
Одним из основных геостатистических методов исследования окружающей среды является метод кригинга, согласно которому оцениваемое значение регионализованной переменной величины ψр в точке p определяется как взвешенное среднее известных наблюдений в соседних точках по формуле:
где Wi — вес i-го значения регионализованной переменной величины ψi по отношению к оцениваемой точке p из k соседних точек.
Метод кригинга предусматривает решение системы уравнений [1]:
где γ(ξij) — значение полувариограммы для расстояния ξij между точками i и j; γ(ξip) — значение полувариограммы для расстояния ξip между известной точкой i и оцениваемой точкой p, λ— множитель Лагранжа.
Однако применение системы уравнений (2) нередко приводит к тому, что некоторые веса Wi оказываются либо меньше нуля, либо больше 1. В частности, такие результаты приведены Дэвисом [1]. Очевидно, это существенно искажает физический и математический смысл весов Wi и приводит к неправильной оценке величины ψp.
Подобные результаты могут быть вызваны либо неточным подбором модели полувариограммы, либо наличием среди точек наблюдений величины ψi точек, не оказывающих пpaктического влияния на оцениваемое значение ψр. Такие точки должны иметь веса, равные нулю.
Дисперсия оценки методом кригинга может быть выражена формулой Матерона [2]:
где - дисперсия ПЭС относительно ее среднего значения; - ковариация среднего значения и значения ; - ковариация значений и .
Система уравнений (2) выводится из условия минимума дисперсии в формуле (3), условия для весов Wi :
а также заменой ковариаций и соответственно на разности и .
Однако, если подставить в формулы (3) и (4) вместо Wi () и выводить систему уравнений кригинга относительно , а затем сделать обратную подстановку (Wi вместо ), то получается система уравнений кригинга
позволяющая учесть ограничения.
Так как решение системы уравнений (5) может быть затруднено из-за ее нелинейности, то для ее компьютерного решения предлагается следующий алгоритм:
Решить систему линейных уравнений (2).
Если полученное решение не отвечает ограничениям , то исключить из рассмотрения точку наблюдения с наименьшим отрицательным весом (или присвоить ему значение ноль), иначе - закончить вычисления.
Перейти после указанных изменений к пункту 1.
Полученное при реализации данного алгоритма (или использования формулы (5)) решение позволяет учесть ограничения и точнее оценить значение регионализованной переменной величины ψр.
Литература
- Дэвис Дж. С. Статистический анализ данных в геологии. Пер. с англ. В 2 кн. Пер. В.А. Гoлyбевой; Под ред. Д.А. Родионова. Кн. 1. — М.: Недра, 1990. — 319 с; Кн. 2. — М.: Недра, 1990. — 427 с.
- Матерон Ж. Основы прикладной геостатисти- ки. — М.: Мир, 1968. — 408с.
Статья в формате PDF 120 KB...
25 04 2024 20:20:36
Статья в формате PDF 576 KB...
24 04 2024 7:47:36
Проводился анализ изменений биоэлектрической активности головного мозга и сверхмедленной активности в нервной, дыхательной и сердечно-сосудистой системах в процессе адаптивного биоуправления с биологической обратной связью по параметрам церебральной гемодинамики и медитации. Осуществлялась регистрация сверхмедленной активности нервной и сердечно-сосудистой систем и локализация биоэлектрической активности нервной системы. Выявлено вовлечение различных мозговых структур в реализацию поведенческих стратегий в группах обучившихся различным видам самоуправления, что говорит о различии механизмов достижения конечного результата. Полученные результаты свидетельствуют о вовлечении кардиореспираторной синхронизации в изменение биоэлектрической активности только при релаксации с помощью адаптивного биоуправления. Осуществлена проверка резонансной гипотезы релаксации, согласно которой при совпадении частот изменения дыхания, биоэлектрической активности мозга, сердечного ритма и сосудистого тонуса происходит усиление активности в вовлекаемых в резонансный ответ структурах. ...
23 04 2024 16:51:40
Статья в формате PDF 119 KB...
22 04 2024 4:48:52
Статья в формате PDF 142 KB...
21 04 2024 1:58:13
Статья в формате PDF 112 KB...
19 04 2024 17:20:18
Статья в формате PDF 139 KB...
18 04 2024 2:53:45
Статья в формате PDF 263 KB...
17 04 2024 14:15:39
16 04 2024 13:21:36
Статья в формате PDF 112 KB...
15 04 2024 10:12:35
Статья в формате PDF 121 KB...
14 04 2024 8:19:28
Статья в формате PDF 269 KB...
12 04 2024 4:18:40
11 04 2024 23:49:58
Статья в формате PDF 202 KB...
10 04 2024 0:24:52
Статья в формате PDF 286 KB...
09 04 2024 4:44:52
Статья в формате PDF 150 KB...
08 04 2024 8:14:29
Статья в формате PDF 106 KB...
07 04 2024 11:50:23
Статья в формате PDF 102 KB...
06 04 2024 20:46:30
Статья в формате PDF 105 KB...
04 04 2024 14:47:36
Статья в формате PDF 112 KB...
03 04 2024 8:49:21
Статья в формате PDF 101 KB...
01 04 2024 22:21:30
Статья в формате PDF 108 KB...
31 03 2024 3:36:19
Статья в формате PDF 125 KB...
30 03 2024 18:19:59
Статья в формате PDF 266 KB...
29 03 2024 1:36:59
Боль является одним из самых распространенных симптомов, встречающихся в медицинской пpaктике. Было изучено влияние фототерапии на интенсивность боли при невропатиях тройничного нерва травматического происхождения. Для лечения использовались фотонные матрицы Коробова «Барва -флекс/КИК» в сочетании с магнитной матрицей «Барва-флекс/МАГ». Ежедневно интенсивность боли оценивалась по визуальной аналоговой шкале боли. Фототерапия оказывает положительное влияние в виде сокращения интенсивности и длительности болевого синдрома. ...
27 03 2024 14:40:23
Статья в формате PDF 287 KB...
26 03 2024 23:17:30
Статья в формате PDF 113 KB...
25 03 2024 10:45:12
Статья в формате PDF 221 KB...
24 03 2024 4:13:38
Статья в формате PDF 206 KB...
22 03 2024 5:50:26
Статья в формате PDF 109 KB...
21 03 2024 21:14:32
Статья в формате PDF 216 KB...
20 03 2024 2:29:12
Статья в формате PDF 108 KB...
19 03 2024 21:31:33
18 03 2024 5:53:34
Статья в формате PDF 126 KB...
17 03 2024 9:44:50
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::