
- •Девиз:
- •Модели и программы
- •Литература:
- •Тема 1: «Экспериментальное обеспечение моделей»
- •Мы изучаем модели:
- •Физическое обоснование. Уравнение неразрывности
- •Вывод уравнения неразрывности
- •Расчетная схема
- •Основные уравнения
- •Основное уравнение переноса в конечно-разностной форме
- •Сеточная схема расчета (используется при численном моделировании методом конечных элементов, -метод Галеркина)
- •Таким образом, экспериментальным обеспечением моделей влагопереноса
- •Начальные условия
- •Модель движения воды в почве
- •Граничные условия
- •Продолжение темы 1
- •Нижнее граничное условие
- •Это было о потоках влаги, о динамике влажности почвы
- •Пункт 2 темы: Перенос тепла в почве qТ dTdz
- •Уравнение теплопереноса
- •Основное уравнение теплопереноса
- •Расчет переноса тепла внутри почвы
- •Уравнение теплопереноса в
- •Условия для применения расчетной схемы
- •Расчет переноса тепла внутри почвы
- •Уравнение теплопереноса в
- •Условия для применения расчетной схемы
- •Условия для применения расчетной схемы
- •стоп
- •Педотрансферные функции
- •«Педотрансферные функции превращают данные, которые у нас есть, в данные, которые нам нужны!»
- •Определяют педотрансферные функции
- •1-й способ (точечный)
- •2-й способ (параметрический)
- •Педотрасферные функции
- •Нейросетевые алгоритмы
- •Математическая модель нейрона
- •«Педотрансферные функции превращают данные, которые у нас есть, в данные, которые нам нужны!»
- •Тема 2. Аппроксимация эмпирических зависимостей
- •Виды функций
- •Убывающие функции
- •С одним максимумом
- •V.С несколькими экстремумами
- •Определение
- •Следующий вопрос темы 2
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •3.Определение параметров аппроксимации
- •Задача – найти минимум S, подобрав соответствующие значения b1 и b2
- •Определение
- •Расчет среднеквадратической
- •Метод сканирования
- •СТАТИСТИКИ ДЛЯ АНАЛИЗА ПАРАМЕТРОВ
- •Анализ ошибок (погрешностей)
- •Практическое задание 1.
- •Зависимость Y от X
- •«Выдача» после аппроксимации
- •Гистограмма погрешностей
- •График реальных значений от расчетных
- •График зависимости погрешностей от расчетной величины
- •Анализ на систематические погрешности
- •Сравнение параметров
- •Статистическое доказательство значимости различий параметров
- •Физический смысл сравнения параметров (пример с термофильными микроорганизмами)
- •Правила подбора функций
- •Практическое занятие №1
- •Модели и программы
- •Тема 2. Аппроксимация эмпирических зависимостей
- •Виды функций
- •Убывающие функции
- •С одним максимумом
- •V.С несколькими экстремумами
- •Определение
- •Следующий вопрос темы 2
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •3.Определение параметров аппроксимации
- •Задача – найти минимум S, подобрав соответствующие значения b1 и b2
- •Определение
- •Расчет среднеквадратической
- •Метод сканирования
- •СТАТИСТИКИ ДЛЯ АНАЛИЗА ПАРАМЕТРОВ
- •Анализ ошибок (погрешностей)
- •Практическое задание 1.
- •Зависимость Y от X
- •«Выдача» после аппроксимации
- •Гистограмма погрешностей
- •График реальных значений от расчетных
- •График зависимости погрешностей от расчетной величины
- •Анализ на систематические погрешности
- •Сравнение параметров
- •Статистическое доказательство значимости различий параметров
- •Физический смысл сравнения параметров (пример с термофильными микроорганизмами)
- •Правила подбора функций
- •Практическое занятие №1
- •Задание после аппроксимации
- •Заключительные вопросы (!)
- ••ЗАДАНИЕ ПО ПАРАМЕТРИЗАЦИИ ЭКСПЕРИМЕНТАЛЬНЫХ ЗАВИСИМОСТЕЙ (АППРОКСИМАЦИЯ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ)
- •Получить параметры аппроксимации для 2-й экспериментальной зависимости. Сравнить достоверность отличий полученных параметров по
- •Пример работы с программой STATISTICA
- •Аппроксимация в программе STATISTICA
- •Сравнение функций для описания экспериментальных данных и выбор лучшей
- •Сравнение моделей. Критерий Вильямса-Клюта
- •Из раздела: Аппроксимация
- •Пример расчета критерия Вильямса-Клюта
- •4. Функции, наиболее употребительные в почвоведении
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •А если мы не знаем вида? Или функция очень сложная? С несколькими экстремумами?
- •Пример Харитоновой
- •4. Функции, наиболее употребительные в почвоведении
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •В прошлый раз… Заканчиваем тему «Аппроксимация экспериментальных данных»
- •Впрограмме STATISTICA
- •Сравнение функций для описания экспериментальных данных и выбор лучшей
- •Сравнение моделей. Критерий Вильямса-Клюта
- •Для одних и тех же данных можно
- •• Критери й Вильямса-Клута (по Рыбалко
- •Пример расчета критерия Вильямса-Клюта
- •4. Функции, наиболее употребительные в почвоведении
- •Вхимии почв: уравнения сорбции
- •В биологии почв, экологии: физический смысл параметров гауссиады
- •Раздел почвоведения, использвание
- •А если мы не знаем вида? Или функция очень сложная? С несколькими экстремумами?
- •Общие правила подбора функций для экспериментальных данных
- •Новая тема
- •Работа с физически обоснованными моделями Модель водо- и солеобмена почв HYDRUS
- •Структура моделей
- •Типы наиболее употребительных моделей
- •2 типа моделей:
- •Составляющие водного баланса почв
- •Водный баланс почвы (одномерный)
- •Балансовые модели основаны на использовании гидрологических
- •Устройство балансовых моделей
- •Имитационные модели ……
- •HYDRUS
- •Что будем считать?
- •«Первая задача с прогнозной моделью HYDRUS»
- •Основные расчетные процессы модели HYDRUS
- •Geometry information
- •Заказ
- •Time information
- •Формирование препроцессора
- •Условия сеточной схемы
- •Параметры ван Генухтена и педотрансферные функции
- •NB! Запомнить, записать.
- •Граничные условия
- •Метеоусловия (условия на
- •Графический редактор Слои (материал 1 и 2)
- •Графический редактор (начальные условия: давление влаги по профилю)
- •Постпроцессор
- •Об умении читать выдачи
- •Пример работы с программой STATISTICA
- •Тема 3 «МОДЕЛИ РАВНОВЕСИЯ И ПЕРЕНОСА ВЕЩЕСТВ В ПОЧВАХ»
- •Ионные равновесия в растворах (гомогенный обмен)
- •Основные уравнения
- •Ионные равновесия при наличии твердой фазы (гетерогенный обмен)
- •3. Модели переноса веществ
- •Конвективный переносC веществ
- •Движение в почве с учетом
- •Явление гидродинамической дисперсии (по Holzbecher,1998)
- •Гидродинамическая дисперсия и
- •Конвективно-диффузионное уравнение
- •Определения
- •Конвективно-диффузионное уравнение
- •Это конвективно-дифф.уравнение, но без учета процессов, которые могут приводить к появлению дополнительных количеств
- •Конвективно-диффузионное уравнение
- •1.Ионный обмен
- •Итак, о слагаемом в конвективно-
- •Нерастворяющий объем (объем выталкивания, exclusive water)
- •1.Ионный обмен
- •Сорбция/десорбция кинетическая
- •Кинетика нулевого порядка
- •Кинетическая сорбция (разложение) веществ по кинетике 1-го порядка
- •Кинетики десорбции/разложения
- •Аналогично сорбции/десорбции описываются и процессы разложения, распада, роста (микроорганизмов) и пр.
- •Мы разобрались с описаниями всех процессов,
- •Основные параметры переноса солей
- •ВЫВОД:
- •Перенос веществ с преимущественными потоками растворов в почве
- •Полевой фильтрационный эксперимент
- •Основные уравнения переноса веществ по макропорам
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Переходим ко второму вопросу
- •Главное, - введение концепции
- •Исследования Оценка риска Управление
- •Вопрос 1 лекции 1
- •Статистический анализ кривой «доза-эффект»
- •Установление реперной концентрации (ВС5) для формальдегида
- •ПДК, ОБУВ и риски. Определения
- •Пример с агротехнологиями
- •ВЫВОД:
- •Полевой фильтрационный эксперимент
- •Основные уравнения переноса веществ по макропорам
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Конвективно-диффузионное уравнение
- •1.Ионный обмен
- •1.Ионный обмен
- •Сорбция/десорбция кинетическая
- •Кинетика нулевого порядка
- •Кинетическая сорбция (разложение) веществ по кинетике 1-го порядка
- •Кинетики десорбции/разложения
- •Аналогично сорбции/десорбции описываются и процессы разложения, распада, роста (микроорганизмов) и пр.
- •Основные параметры переноса солей
- •ВЫВОД:
- •Перенос веществ с преимущественными потоками растворов в почве
- •Полевой фильтрационный эксперимент
- •Основные уравнения переноса веществ по макропорам
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Переходим ко второму вопросу
- •Главное, - введение концепции
- •Исследования Оценка риска Управление
- •Вопрос 1 лекции 1
- •Статистический анализ кривой «доза-эффект»
- •Установление реперной концентрации (ВС5) для формальдегида
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Блок метеоусловий
- •Симулятор погоды
- •Сравнение моделей. Критерий Сайерта
- •Основные параметры модели
- •Схема передвижения веществ к корню растения
- •Основные уравнения модели
- •Определение чувствительности модели
- •Анализ чувствительности модели
- •Анализ чувствительности модели к параметрам ван Генухтена
- •Оценка
- •Тема IV
- •разговор

Практическое занятие №1
Аппроксимация экспериментальных данных

Задание после аппроксимации
•Открыть PARAMETERS&STAND. ERR.: списать значения параметров, их ошибок, t-критерия и уровня вероятности. Далее провести анализ погрешностей:
•DISTRIB.OF RESIDUALS. Посмотреть и проанализировать гистограмму распределения погрешностей. Выводы о погрешностях. Здесь просмотреть нормальность распределения погрешностей в услуге «Normal Probability Plot of Residuals».
•PREDICTEDvsOBSERVED. Проанализировать совпадение экспериментальных и расчетных.
•PREDICTEDvsRESIDUALS. Проанализировать изменение погрешностей с ростом функции. Вывод о наличии систематических ошибок.
•Если аппроксимация прошла неудачно, вернуться к п.3, и повторить те же операции с другим видом функции.
•Аппроксимировать значения 2-й функции: зависимости переменной У11 от Х1. Получить параметры аппроксимации для 2-й экспериментальной зависимости.
•Сравнить достоверность отличий полученных параметров по t-критерию (табличные t7=2.36; t8=2.31).

Заключительные вопросы (!)
Ответить на вопросы:
•(1) что указывает на хорошую аппроксимацию,
•(2) достоверны ли параметры у 2-х аппроксимаций экспериментальных зависимостей,
•(3) о чем говорят различия параметров, если эти функции отражают физический (биологический, и пр.) процесс,
•(4) о чем говорят распределения остатков.

•ЗАДАНИЕ ПО ПАРАМЕТРИЗАЦИИ ЭКСПЕРИМЕНТАЛЬНЫХ ЗАВИСИМОСТЕЙ (АППРОКСИМАЦИЯ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ)
Даны экспериментальные данные в виде значений независимой переменной (например, Х1) и двух зависимых (например, Y1 и Y11). Открыть пакет STATISTICA. Набить данные в столбцах 1, 2 и 3 (VAR1, VAR2 и VAR3).
Определить тип полученной функции. Уточнить вид зависимости (возрастающий, убывающий, с максимумом). Проделать аналогично для Y1 и Y11.
Статистика > Дополнительные линейные/нелинейные модели > Нелинейная оценка. Выбрать вид аппроксимации User Specified regression, lеast squares>OK> Function to be estimated. В поле ввода Estimated function
задать вид функции, в котором независимая и зависимая переменные обозначаются названиями соответствующих столбцов (V1, V2), а оцениваемые параметры — b1, b2 и т.д. Открыть таблицу Summary: Parameter estimation.
Переписать значения параметров, их ошибок, t-критерия и уровня вероятности. Далее провести анализ погрешностей, открыв вкладку Residuais:
Посмотреть и проанализировать гистограмму распределения погрешностей Histogram of residuals. Сделать выводы о погрешностях. Просмотреть нормальность распределения погрешностей на графике Normal Probability Plot of Residuals.
Проанализировать совпадение экспериментальных и расчетных данных
Observed vs. Predicted.
Проанализировать изменение погрешностей с ростом функции. Сделать вывод о наличии систематических ошибок Residuals vs. Predicted.

Получить параметры аппроксимации для 2-й экспериментальной зависимости. Сравнить достоверность отличий полученных параметров по t-критерию (табличные t7=2.36; t8=2.31).
Ответить на вопросы:
1)Что указывает на хорошую аппроксимацию?
2)Достоверны ли параметры у 2-х аппроксимаций экспериментальных зависимостей?
3)О чём говорят различия параметров, если эти функции отражают физический (биологический, и пр.) процесс?
4)О чём говорят распределения остатков?
9. Подберите 2 модели для Ваших данных. Оценить достоверность различий моделей, докажите, какая из 2-х лучше по критерию Вильямса-Клюта.
Ваш вариант 1 задачи :
Вариант 10 |
|
|
|
x10 |
y10 |
|
y1010 |
|
14,3 |
1,6 |
2,1 |
|
30 |
2 |
2,7 |
|
60 |
2,225 |
3 |
|
70 |
2,3 |
3,5 |
|
100 |
2,5 |
4 |
|
135 |
2,6 |
4,3 |
|
151 |
2,676 |
4,5 |
|
164 |
2,7 |
4,6 |
|
185 |
2,8 |
4,767 |
|
400 |
3,3 |
6 |

Пример работы с программой STATISTICA
Var1 Var2
11.000 0.600
22.500 1.010
34.000 2.980
45.000 4.960
57.000 6.000
6 |
9.500 |
4.100 |
7 |
10.500 |
2.500 |
8 |
11.500 1.500 |
|
9 |
13.000 |
0.500 |
10 |
15.000 |
0.210 |

Аппроксимация в программе STATISTICA

Сравнение функций для описания экспериментальных данных и выбор лучшей
Анализ погрешностей (ошибок) Yэксп – Yрасч.
yэ Yр
Эти ошибки находятся в файл выдачи в табличной форме. Эту таблицу надо сохранить и использовать для расета критерия Вильямса-Клюта.

Сравнение моделей. Критерий Вильямса-Клюта
Vi 12 ( 'i ''i )
Ui ( i'' i' )
Vi a Ui
строят регрессию. Чем больше а, тем значительнее
отличаются модели. Дальше достоверность а по критерию Стьюдент Если а достоверно, то модели достоверно различаются.
Если t-критерий >0, то 1-я модель лучше, если tа<0, то 2-я лучше.

Из раздела: Аппроксимация
(пример 1)
Реальные
данные
X |
Y |
x |
b2 |
|
|||
|
|
||||||
2 |
1.5 |
y ' |
|
|
|
|
|
|
|
||||||
|
|
b1 |
|
|
|||
5 |
2.4 |
y" b1 exp b2 |
x |
||||
10 |
3.0 |
||||||
|
|
|
|
|
|||
20 |
4.0 |
|
|
|
|
|

• |
Критери й Вильямса-Клута (по |
|
||
|
Рыбалко |
z Yr y1 |
y2 |
|
|
|
|||
«МАТЕМАТИЧЕСКИЕ МОДЕЛИ КОНТРОЛЯ НАДЁЖНОСТИ ОБЪЕКТОВ |
||||
|
ЭНЕРГЕТИКИ », 2010) |
|
|
|
|
|
|
2 |
|
• |
Рассчитывают параметр Z : t y1 y2 |
|||
• |
Рассчитывают параметр t: Z t |
|
||
• |
Находят регрессию Z от t: |
|
|
|
Если λ достоверна и положительна, то 1-я модель лучше
•При положительном и статистически значимом критерии λ первая регрессия лучше, при отрицательном и значимом – вторая. Если λ незначим, то модели не различабтся.