Построение математической модели и анализ характеристических свойств системы методом планирования эксперимента
Формулировка задачи управления системой. Выявление входных воздействий, выходных переменных, возмущающих и управляющих воздействий и управляемых переменных. Получение математической модели, описывающей динамику изучаемой системы. Анализ свойств системы.
Рубрика | Математика |
Вид | задача |
Язык | русский |
Дата добавления | 27.03.2016 |
Размер файла | 17,8 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Новочеркасский инженерно-мелиоративный институт имени А.К. Кортунова
ФГБОУ ВПО «ДОНСКОЙ ГОСУДАРСТВЕННЫЙ АГРАРНЫЙ УНИВЕРСИТЕТ»
Кафедра менеджмента и информатики
Задание
для выполнения курсовой работы
по дисциплине "ИССЛЕДОВАНИЕ СИСТЕМ УПРАВЛЕНИЯ"
Тема: "Построение математической модели и анализ характеристических свойств системы методом планирования эксперимента"
ИСХОДНЫЕ ДАННЫЕ ДЛЯ РАСЧЕТА
№ 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Yоп
SЭ 2 = ; = 0,05 fЭ = 2; t = 4,303
На основе предложенной информации:
- определить тип изучаемого процесса, а также факторы, влияющие на уровень изменения этого процесса;
- рассчитать коэффициенты уравнения линейной модели;
- проверить полученную модель на адекватность;
- достроить полученную модель до квадратичной модели;
- оценить значимость коэффициентов регрессии;
- сделать перерасчет коэффициентов модели, если есть незначимые коэффициенты при квадратах;
- проверить квадратичную модель на адекватность;
- определить координаты функции отклика в зоне min и max;
- построить двумерные сечения поверхности отклика для трех вариантов;
- провести анализ проведения функции отклика по полученным графикам с использованием дополнительных сечений;
- произвести анализ проведения функции отклика в зоне минимума, максимума и в центре факторного пространства;
- произвести анализ степени риска рассчитанных коэффициентов модели;
- сделать основные выводы о произведенных расчетах, возможности использования модели, принятых ограничениях и т.д.
Введение
Существует несколько направлений для изучения различных систем, и в том числе сложных систем. Эти системы характеризуются значительным числом взаимосвязанных параметров. Задача исследования таких систем заключается в установлении зависимости между входными параметрами-факторами и выходными параметрами-показателями качества функционирования системы и определении уровней факторов, оптимизирующих выходные параметры системы.
В условиях неполного знания механизма явлений задачи идентификации и оптимизации, т.е. отыскания оптимальных условий протекания процессов, решаются с помощью экспериментально-статистических методов -- математическое планирование эксперимента. Математическое планирование эксперимента -- это процедура выбора числа и условий проведения опытов, необходимых и достаточных для решения данной задачи с требуемой точностью, методов математической обработки их результатов и принятия решений.
Достоинством этого математического аппарата является его универсальность, пригодность в большинстве областей исследования (промышленности, медицине, биологии, сельском хозяйстве и др.).
Основная задача математического планирования эксперимента -- переработка исходной информации для целей управления процессами и системами. При исследовании систем управления необходимо:
- сформулировать задачу управления системой;
- выявить входные воздействия (факторы), выходные переменные, возмущающие и управляющие воздействия и управляемые переменные;
- получить математическую модель, описывающую динамику изучаемой системы;
- провести анализ характеристических свойств системы как объекта управления (управляемость, устойчивость, наблюдаемость и т.п.).
1. Постановка задачи
переменный динамика математический возмущающий
Исследовать урожайность фруктового сада, при влиянии следующих факторов: расстояние между деревьями, объем вносимых удобрений, норма полива.
Значения факторов представлены в таблице 1.
Таблица 1 - Факторное пространство
Факторы |
Значения |
Единицы измерения |
||||||
максимальное |
минимальное |
среднее |
||||||
н.ед.из |
код |
н.ед.из |
код |
н.ед.из |
код |
|||
X1 - расстояние между деревьями |
5 |
+1 |
1 |
-1 |
3 |
0 |
м |
|
X1 - норма полива |
3000 |
+1 |
1000 |
-1 |
2000 |
0 |
м3\га |
|
X1 - объем удобрений |
2,5 |
+1 |
0,5 |
-1 |
1,5 |
0 |
т\га |
2. Построение линейной регрессионной модели
Таблица 2 - Исходная матрица
N |
X0 |
X1 |
X2 |
X3 |
X1 X2 |
X1 X3 |
X2 X3 |
X1 X2 X3 |
YОП |
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
1 |
+1 |
-1 |
+1 |
+1 |
-1 |
-1 |
+1 |
-1 |
24,3 |
|
2 |
+1 |
-1 |
-1 |
+1 |
+1 |
-1 |
-1 |
+1 |
41,5 |
|
3 |
+1 |
+1 |
-1 |
+1 |
-1 |
+1 |
-1 |
-1 |
14,2 |
|
4 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
34,6 |
|
5 |
+1 |
+1 |
+1 |
-1 |
+1 |
-1 |
-1 |
-1 |
27,9 |
|
6 |
+1 |
-1 |
+1 |
-1 |
-1 |
+1 |
-1 |
+1 |
32,5 |
|
7 |
+1 |
-1 |
-1 |
-1 |
+1 |
+1 |
+1 |
-1 |
25,6 |
|
8 |
+1 |
+1 |
-1 |
-1 |
-1 |
-1 |
+1 |
+1 |
20,3 |
|
? |
220,9 |
-26,9 |
17,7 |
8,3 |
38,3 |
-7,1 |
-11,3 |
36,9 |
Таблица 3 - Определение значимости коэффициентов bi
Параметры |
Коэффициенты |
||||||||
b0 |
b1 |
b2 |
b3 |
b12 |
b13 |
b23 |
b123 |
||
bi |
27,61 |
-3,36 |
2,21 |
1,04 |
4,79 |
-0,89 |
-1,41 |
4,61 |
|
bi кр |
1,92 |
||||||||
Значимые bi |
27,61 |
-3,36 |
2,21 |
4,79 |
4,61 |
В результате получаем окончательный вид уравнения регрессии, включающего значимые коэффициенты:
Yтеор= b0-b1x1+b2x2+ b12x1x2+ b123x1x2 x3
Yтеор=27,61-3,36х1+2,21x2+4,79 x1x2+4,61 x1x2 x3
Таблица 4 - Расчет дисперсии неадекватности
№ |
YОП |
Yтеор |
(YОП - Yтеор)2 |
% ошибки |
|
1 |
24,3 |
23,78 |
0,27 |
2,14 |
|
2 |
41,5 |
38,16 |
11,16 |
8,05 |
|
3 |
14,2 |
12,64 |
2,43 |
10,99 |
|
4 |
34,6 |
35,86 |
1,59 |
-3,64 |
|
5 |
27,9 |
26,64 |
1,59 |
4,52 |
|
6 |
32,5 |
33,00 |
0,25 |
-1,54 |
|
7 |
25,6 |
28,94 |
11,16 |
-13,05 |
|
8 |
20,3 |
21,86 |
2,43 |
-7,68 |
|
Сумма |
30,874 |
Дисперсия неадекватности S2 на = 23,81
Опытное значение критерия Фишера Fоп=
Теоретическое значение критерия Фишера Fтеор =
3. Построение квадратичной регрессионной модели
Таблица 5 - Расширенная матрица планирования
№ |
X0 |
X1 |
X2 |
X3 |
X1 2 |
X22 |
X32 |
X1 X2 |
X1 X3 |
X2 X3 |
X1 X2 X3 |
Уоп |
|
1 |
+1 |
-1 |
+1 |
+1 |
+1 |
+1 |
+1 |
-1 |
-1 |
+1 |
-1 |
24,3 |
|
2 |
+1 |
-1 |
-1 |
+1 |
+1 |
+1 |
+1 |
+1 |
-1 |
-1 |
+1 |
41,5 |
|
3 |
+1 |
+1 |
-1 |
+1 |
+1 |
+1 |
+1 |
-1 |
+1 |
-1 |
-1 |
14,2 |
|
4 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
34,6 |
|
5 |
+1 |
+1 |
+1 |
-1 |
+1 |
+1 |
+1 |
+1 |
-1 |
-1 |
-1 |
27,9 |
|
6 |
+1 |
-1 |
+1 |
-1 |
+1 |
+1 |
+1 |
-1 |
+1 |
-1 |
+1 |
32,5 |
|
7 |
+1 |
-1 |
-1 |
-1 |
+1 |
+1 |
+1 |
+1 |
+1 |
+1 |
-1 |
25,6 |
|
8 |
+1 |
+1 |
-1 |
-1 |
+1 |
+1 |
+1 |
-1 |
-1 |
+1 |
+1 |
20,3 |
|
9 |
+1 |
+1 |
0 |
0 |
+1 |
0 |
0 |
0 |
0 |
0 |
0 |
33,8 |
|
10 |
+1 |
-1 |
0 |
0 |
+1 |
0 |
0 |
0 |
0 |
0 |
0 |
27,4 |
|
11 |
+1 |
0 |
+1 |
0 |
0 |
+1 |
0 |
0 |
0 |
0 |
0 |
26,0 |
|
12 |
+1 |
0 |
-1 |
0 |
0 |
+1 |
0 |
0 |
0 |
0 |
0 |
17,1 |
|
13 |
+1 |
0 |
0 |
+1 |
0 |
0 |
+1 |
0 |
0 |
0 |
0 |
14,7 |
|
14 |
+1 |
0 |
0 |
-1 |
0 |
0 |
+1 |
0 |
0 |
0 |
0 |
26,3 |
|
? |
366,2 |
-20,5 |
26,6 |
-3,3 |
282,1 |
264 |
261,9 |
38,3 |
-7,1 |
-11,3 |
36,9 |
Таблица 6 - Определение значимости коэффициентов квадратичной модели
Параметры |
Коэффициенты |
|||||||||||
b0 |
b1 |
b2 |
b3 |
b11 |
b22 |
b33 |
b12 |
b13 |
b23 |
b123 |
||
Расчетные значения коэффициентов bi |
22,48 |
-2,05 |
2,66 |
-0,33 |
8,103 |
-0,947 |
-1,997 |
4,79 |
-0,89 |
-1,41 |
4,61 |
|
Критические значения коэффициентов bi |
3,47 |
1,72 |
3,47 |
1,92 |
||||||||
Значимые коэффициенты модели |
22,48 |
-2,05 |
2,66 |
8,103 |
4,79 |
4,61 |
Повторный расчет коэффициентов b0 и bii:
b0=21,025
b11=7,185
bii kp=3,22
Уравнение квадратичной регрессионной модели:
Yтеор=21,025-2,05х1+2,66x2+7,185 x12+4,79 x1x2 +4,6179 x1x2 x3
Таблица 7 - Расчетные данные для определения адекватности квадратичной модели
№ опыта |
YОП |
Yтеор |
(YОП - Yтеор)2 |
% ошибки |
|
1 |
2 |
3 |
4 |
5 |
|
1 |
24,3 |
23,52 |
0,61 |
3,21 |
|
2 |
41,5 |
37 |
20,25 |
10,84 |
|
3 |
14,2 |
14,1 |
0,01 |
0,70 |
|
4 |
34,6 |
38,22 |
13,10 |
-10,46 |
|
5 |
27,9 |
29 |
1,21 |
-3,94 |
|
6 |
32,5 |
32,74 |
0,06 |
-0,74 |
|
7 |
25,6 |
27,78 |
4,75 |
-8,52 |
|
8 |
20,3 |
23,32 |
9,12 |
-14,88 |
|
9 |
33,8 |
26,16 |
58,37 |
22,60 |
|
10 |
27,4 |
30,26 |
8,18 |
-10,44 |
|
11 |
26,0 |
23,685 |
5,36 |
8,90 |
|
12 |
17,1 |
18,365 |
1,60 |
-7,40 |
|
13 |
14,7 |
21,025 |
40,01 |
-43,03 |
|
14 |
26,3 |
21,025 |
27,83 |
20,06 |
|
? |
190,45 |
S2 на = 23,81
Fоп= 14,88
Fтеор =19,365
Вывод
Для уmax (x1x2) фактор x1 принимает максимальное значения от 0,7 до 1, минимальные - от -1 до -0,6. Фактор x2 принимает максимальные значения от 0,1 до 1, минимальные - от -1 до -0,7 (рисунок 1).
Для уmax (x1x3) максимальное значение находится для фактора x1 в диапазоне от 0,7 до 1, минимальное - от 0 до 1. Для фактора x3 максимальное значение в диапазоне от -0,9 до 1, минимальное - от -0,5 до 1 (рисунок 2).
Для уmax (x2x3) максимальное значение фактора x2 находится в диапазоне от 0,1 до 1, минимальное также принимает такие же значения. Для фактора x3 максимальное значение в диапазоне от 0,7 до 1, минимальное - от -1 до -0,7 (рисунок 3).
Размещено на Allbest.ru
Подобные документы
Решение дифференциальных уравнений математической модели системы с гасителем и без гасителя. Статический расчет виброизоляции. Определение собственных частот системы, построение амплитудно-частотных характеристик и зависимости перемещений от времени.
контрольная работа [1,6 M], добавлен 22.12.2014Знакомство с особенностями построения математических моделей задач линейного программирования. Характеристика проблем составления математической модели двойственной задачи, обзор дополнительных переменных. Рассмотрение основанных функций новых переменных.
задача [656,1 K], добавлен 01.06.2016Создание математической модели движения шарика, подброшенного вертикально вверх, от начала падения до удара о землю. Компьютерная реализация математической модели в среде электронных таблиц. Определение влияния изменения скорости на дальность падения.
контрольная работа [1,7 M], добавлен 09.03.2016Составление гамильтониан Н с учетом необходимых условий оптимальности для задачи Майера. Определение оптимального управления из условия максимизации. Получение конической системы уравнений и ее разрешение. Анализ необходимых условий оптимальности.
курсовая работа [113,1 K], добавлен 13.09.2010Построение сигнального графа и структурной схемы системы управления. Расчет передаточной функции системы по формуле Мейсона. Анализ устойчивости по критерию Ляпунова. Синтез формирующего фильтра. Оценка качества эквивалентной схемы по переходной функции.
курсовая работа [462,5 K], добавлен 20.10.2013Составление математической модели задачи. Приведение ее к стандартной транспортной задаче с балансом запасов и потребностей. Построение начального опорного плана задачи методом минимального элемента, решение методом потенциалов. Анализ результатов.
задача [58,6 K], добавлен 16.02.2016Изучение актуальной задачи математического моделирования в биологии. Исследование модифицированной модели Лотки-Вольтерра типа конкуренция хищника за жертву. Проведение линеаризации исходной системы. Решение системы нелинейных дифференциальных уравнений.
контрольная работа [239,6 K], добавлен 20.04.2016Способы решения системы линейных алгебраических уравнений: по правилу Крамера, методом матричным и Жордана-Гаусса. Анализ решения задачи методом искусственного базиса. Характеристика основной матрицы, составленной из коэффициентов системы при переменных.
контрольная работа [951,8 K], добавлен 16.02.2012Байесовские алгоритмы оценивания (фильтр Калмана). Постановка задачи оценивания для линейных моделей динамической системы и измерений. Запись модели эволюции и модели измерения в матричном виде. Составление системы уравнений, описывающей эволюцию системы.
курсовая работа [3,0 M], добавлен 14.06.2011Поиск собственных чисел и построение фундаментальной системы решений. Исследование зависимости жордановой формы матрицы А от свойств матрицы системы. Построение фундаментальной матрицы решений методом Эйлера, решение задачи Коши и построение графиков.
курсовая работа [354,7 K], добавлен 14.10.2010