- •Конспект лекций
- •Экспериментальные статистические методы моделирования
- •230100 Информатика и вычислительная техника
- •Тема 7. Динамический регрессионный анализ
- •Тема 8. Построение моделей динамики типа "вход-выход" нелинейных объектов
- •Тема 9. Модели пространства состояний
- •Тема 16. Управление экспериментом при построении нелинейно-параметризованной d-оптимальной линейно-комбинационной модели
- •Тема 17. Планирование эксперимента в условиях неоднородностей
- •Тема 1. Введение
- •1.1.Типы математических моделей
- •1.2. Представление объекта управления в реальных условиях эксплуатации типы математических моделей
- •Представление объекта управления в реальных условиях эксплуатации
- •Тема 18. Планирование эксперимента при наличии количественных и качественных факторов
- •Постановка задачи оценивания неизвестных парметров моделей статики
- •Элементы теории случайных величин
- •Свойства оценок
- •Векторные случайные величины и их характеристики
- •Тема 3. Методы оценивания неизвестных параметров регрессионных моделей
- •Метод наименьших квадратов
- •- Квадратичная формa. Корреляционная функция случайного процесса и ее свойства
- •Свойства автокорреляционных функций.
- •Спектральная плотность случайного процесса
- •Обобщенный метод наименьших квадратов
- •Омнк при малых выборках
- •Oмнк при большом объеме выборки
- •Метод максимального правдоподобия
- •Цель и этапы статистического анализа уравнений регрессии
- •Проверка значимости.
- •Проверка адекватности.
- •Расчет доверительных интервалов.
- •Тема 5. Построение моделей статики в случае нелинейной параметризации
- •Нелинейные методы оценивания
- •Метод Гаусса-Ньютона
- •Метод Хартли
- •Метод Марквардта
- •Тема 6. Построение моделей динамики объектов управления по экспериментальным данным
- •Понятие моделей динамики и их типы
- •Способы построения моделей импульсной переходной функции
- •2.Разложение импульсной переходной функции по заданной системе базисных функций.
- •Корреляционный метод идентификации
- •Методы решения уравнения винера-хопфа в дискретной форме
- •Вычисление оценки вектора неивестных параметров на основе первых уравнений
- •Вычисление оценки вектора неизвестных параметров на основе уравнений.
- •Использование корреляционного метода при подаче на вход псевдослучайного двоичного сигнала (псдс).
- •Задачи динамического регрессионного анализа.
- •Представление импульсной переходной функции набором ординат, вычисление оценок ординат.
- •Разложение импульсной переходной функции по заданной системе базисных функций.
- •Построение модели при известном значении параметра .
- •Построение модели при неизвестном значении параметра .
- •Тема 8. Построение моделей динамики типа "вход-выход" нелинейных объектов
- •Тема 14. D - оптимальное планирование
- •14.2. Процедура непрерывного планирования (пнп)
- •Тема 15. Последовательное планирование эксперимента. Планирование эксперимента в динамике
- •Построение d-оптимальных моделей, соответствующих разложению импульсной переходной функции по системе базисных функций
- •Планирование эксперимента в условиях дискретного дрейфа
- •Библиографический список
Проверка адекватности.
Проверка
адекватности сводится к проверке
гипотезы о равенстве дисперсии
адекватности и дисперсии внешнего
шума:
.
Дисперсия адекватности характеризует точность предсказания выхода по модели, то есть это дисперсия случайной величины, которая является разностью между экспериментальным значением выхода и значением выхода, рассчитанным по модели.
1.
-оценка
дисперсии внешнего шума.
2.
Отношение
- F-распределение (по закону Фишера)
Если
выполняется неравенство
,
то модель признается адекватной
экспериментальным данным. В противном
случае, модель не адекватна, то есть
структура модели выбрана неверно. Ее
необходимо изменить и для новой структуры
провести все исследование.
Квантиль
(
)
должен удовлетворять трем условиям
,
где
-уровень значимости;
-
число степеней свободы числителя
дисперсионного отношения;
-
число степеней свободы знаменателя
дисперсионного отношения.
Для определения используется справочник, но предварительно выбирается конкретное значение .
Расчет доверительных интервалов.
Задаем
конкретное значение входа
:
Доверительный интервал- это интервал возможных значений нашей величины.
В
рассчитанный интервал с вероятностью
попадает истинное значение
.
Ширина доверительного интервала определяет точность модели регрессии. Чем уже доверительный интервал, тем более точная модель регрессии.
Тема 5. Построение моделей статики в случае нелинейной параметризации
5.1. Нелинейные методы оценивания
5.2. Метод Гаусса-Ньютона
5.3. Метод Хартли
5.4. Метод Марквардта
Нелинейные методы оценивания
Нелинейно-параметризованная
модель может быть представлена в виде:
Оценивание неизвестных параметров нелинейно-параметризованных моделей возможно нелинейными методами наименьших квадратов. Наиболее распространенными являются:
метод Гаусса-Ньютона; метод Хартли; метод Марквардта.
Метод Гаусса-Ньютона
Имеется модель
Задается исходное начальное
приближение
.
Функция регрессии раскладывается
в ряд Тейлора в окрестности начального
приближения, при этом в разложении
используют только свободные член и
линейный член, то есть осуществляют
линеаризацию модели в окрестности
.
По
МНК находят оценку
,
то есть оценку разности между истинным
значением
и
:
.
Находится приближение:
.
Далее функция раскладывается в ряд
Тейлора в окрестности
.
Находится
и так далее.
Процедура заканчивается
на той операции, на которой норма
вектора
становится меньше некоторой, достаточно
малой, величины
.
Получена
линеаризованная по параметрам модель,
в которой неизвестны коэффициенты
.
Оценка
возможна МНК.
Каждый элемент суммы представляет собой матрицу размером q*q, являющуюся информационной матрицей Фишера для одного измерения входа-выхода.
В общем виде алгоритм Гаусса-Ньютона можно записать следующим образом:
Метод Хартли
Процедура Хартли
отличается от процедуры Гаусса-Ньютона
тем, что используются не все значения
,
а только часть его.
Метод Марквардта
Метод
Марквардта отличается тем, что вместо
матрицы М
используется матрица
,
где I-
единичная матрица. Говорят, что матрица
лучше обусловлена (имеет определитель
больше, чем у матрицы М).
