Автор работы: Пользователь скрыл имя, 28 Декабря 2011 в 16:44, курсовая работа
В данной работе целью моделирования является изучение изменения состояния объекта в фазовом пространстве в виде явных функции координат и времени.
Задачами курсовой работы является:
1. Разработка имитационной модели изменения пространственно-временного состояния объекта в трехмерном пространстве относительно неподвижной системы координат
2. Построение концептуальной модели изменения пространственно-временного состояния объекта в трехмерном пространстве
3. Построение модели изменения состояния объекта в фазовом и Гильбертовом пространствах
1. Введение
2. Теоретическая часть……………………………………………………………………….…...
3. Практическая часть
3.1. Разработка имитационной модели изменения пространственно-временного состояния объекта в трехмерном пространстве относительно неподвижной системы координат.
3.2. Построение концептуальной модели изменения пространственно-временного состояния объекта в трехмерном пространстве.
3.3 Построение модели изменения состояния объекта в фазовом и Гильбертовом пространствах.
3.4 Оценка математической модели пространственно-временного состояния объекта.
3.5 Оценка и анализ результатов моделирования.
3.6 Статистический метод оценки изменения пространственно-временного состояния объекта.
3.7 Прогнозирование функции отклика объекта на изменение его геометрических свойств.
4. Заключение..........................................................................................
5. Список литературы...................
Обычно предполагают, что ошибки исходных данных распределены по нормальному закону с нулевым математическим ожиданием и дисперсией . При этих предположениях для оценки точности можно применить метод Монте-Карло и, выполнив достаточное число опытов, оценить точность моделирования вертикальных движений сооружения. Более простой, но менее надёжный путь оценки точности состоит в применении известных методов теории ошибок. Основной недостаток в этом случае состоит в необходимости линеаризации оцениваемых функций. Чтобы избежать принятия гипотез о функции распределения ошибок и необходимости линеаризации оцениваемых функций, используем метод имитационного моделирования для оценки точности результатов моделирования, полагая, что нам известна предельная абсолютная погрешность определения исходных данных, связанная со средней квадратической погрешностью выражением .
Для этого положим, что
, (11)
где – вектор ошибок исходных данных;
– вектор ошибок результатов моделирования.
Следовательно,
, (12)
и для
выполнения оценки точности результатов
моделирования достаточно вычислить
вектор (12). Координаты вектора
зададим равным предельной погрешности
исходных данных. Задача оценки точности
результатов моделирования сводится к
определению предельных значений фазовых
координат
и
. Полагая предельную погрешность исходных
данных равной ±D, можно определить векторы
и
значений левой и правой границ интервалов
в пределах, в которых должны находиться
значения исходных данных
Имея
значения
исходных данных, можно вычислить
соответствующие значения фазовых координат
и
и графически или аналитически определить
неустойчивые состояния объекта, где реальные
значения фазовых координат превосходят
предельно допустимые.
Таблица 11. Результаты моделирования изменения пространственного состояния объекта
номер цикла | дата | μα(t) | μ(t) | μβ(t) | αα(t) | α(t) | αβ(t) |
0 | 0 | 147,620 | 147,690555 | 147,760976 | 0,000000 | 0,000000 | 0,000000 |
1 | 0,14 | 147,705 | 147,775886 | 147,846314 | 0,000484 | 0,000484 | 0,000483 |
2 | 1,02 | 147,698 | 147,768801 | 147,839226 | 0,000338 | 0,000338 | 0,000338 |
3 | 1,16 | 147,717 | 147,786944 | 147,857379 | 0,000545 | 0,000545 | 0,000545 |
4 | 2,09 | 147,720 | 147,790605 | 147,861030 | 0,000368 | 0,000368 | 0,000367 |
5 | 3,12 | 147,706 | 147,776202 | 147,846629 | 0,000410 | 0,000409 | 0,000409 |
6 | 4,17 | 147,717 | 147,787160 | 147,857591 | 0,000501 | 0,000501 | 0,000500 |
7 | 5 | 147,699 | 147,769743 | 147,840172 | 0,000510 | 0,000510 | 0,000509 |
8 | 6,1 | 147,712 | 147,782113 | 147,852538 | 0,000456 | 0,000456 | 0,000455 |
9 | 6,28 | 147,709 | 147,779571 | 147,849997 | 0,000404 | 0,000404 | 0,000403 |
10 | 7,18 | 147,716 | 147,786509 | 147,856939 | 0,000375 | 0,000374 | 0,000374 |
11 | 8,14 | 147,692 | 147,762811 | 147,833238 | 0,000381 | 0,000381 | 0,000381 |
12 | 10,21 | 147,710 | 147,780205 | 147,850637 | 0,000516 | 0,000516 | 0,000516 |
13 | 12 | 147,796 | 147,866621 | 147,937059 | 0,000546 | 0,000546 | 0,000545 |
14 | 12,60 | 147,7881 | 147,858549 | 147,928987 | 0,000540 | 0,000539 | 0,000539 |
Рисунок
12. График фазовых координат μ, α
для дискретных данных X
Таблица 12. Результаты моделирования изменения пространственного состояния объекта
номер цикла | дата | μα(t) | μ(t) | μβ(t) | αα(t) | α(t) | αβ(t) |
0 | 0 | 65,076 | 65,144071 | 65,212391 | 0,000000 | 0,000000 | 0,000000 |
1 | 0,14 | 65,171 | 65,239059 | 65,307404 | 0,001171 | 0,001169 | 0,001168 |
2 | 1,02 | 65,157 | 65,225732 | 65,294076 | 0,001031 | 0,001029 | 0,001028 |
3 | 1,16 | 65,135 | 65,203275 | 65,271615 | 0,000921 | 0,000920 | 0,000919 |
4 | 2,09 | 65,155 | 65,223240 | 65,291583 | 0,001060 | 0,001058 | 0,001057 |
5 | 3,12 | 65,151 | 65,219107 | 65,287459 | 0,001217 | 0,001215 | 0,001213 |
6 | 4,17 | 65,178 | 65,246758 | 65,315104 | 0,000961 | 0,000959 | 0,000958 |
7 | 5 | 65,168 | 65,235976 | 65,304314 | 0,000989 | 0,000987 | 0,000986 |
8 | 6,1 | 65,148 | 65,216344 | 65,284701 | 0,001224 | 0,001222 | 0,001220 |
9 | 6,28 | 65,179 | 65,247343 | 65,315686 | 0,001102 | 0,001100 | 0,001098 |
10 | 7,18 | 65,156 | 65,224428 | 65,292769 | 0,001038 | 0,001036 | 0,001035 |
11 | 8,14 | 65,152 | 65,220701 | 65,289043 | 0,000962 | 0,000960 | 0,000959 |
12 | 10,21 | 65,167 | 65,235162 | 65,303504 | 0,000842 | 0,000841 | 0,000839 |
13 | 12 | 65,247 | 65,314910 | 65,383282 | 0,001401 | 0,001398 | 0,001395 |
14 | 12,60 | 65,2391 | 65,307410 | 65,375779 | 0,001355 | 0,001352 | 0,001349 |
Рисунок
13. График фазовых координат μ, α для дискретных
данных Y
Таблица 13. Результаты
моделирования изменения пространственного
состояния объекта
|
Рисунок
14. График фазовых координат μ, α для дискретных
данных Z
Таблица 14. Предельные значения фазовых координат для дискретных данных X
номер цикла | дата | Xа(t) | X(t) | Xb(t) | βi-αi | μi-μ1 | S(t) |
0 | 0,00 | 5,199497 | 5,201978 | 5,204458 | 0,005 | 0,000000 | устойчиво |
1 | 0,14 | 5,203138 | 5,205619 | 5,208100 | 0,005 | 0,003641 | устойчиво |
2 | 1,02 | 5,202769 | 5,205249 | 5,207730 | 0,005 | 0,003271 | устойчиво |
3 | 1,16 | 5,203563 | 5,206043 | 5,208525 | 0,005 | 0,004066 | устойчиво |
4 | 2,15 | 5,203613 | 5,206093 | 5,208574 | 0,005 | 0,004116 | устойчиво |
5 | 3,12 | 5,203112 | 5,205592 | 5,208073 | 0,005 | 0,003615 | устойчиво |
6 | 3,30 | 5,203556 | 5,206036 | 5,208517 | 0,005 | 0,004059 | устойчиво |
7 | 4,17 | 5,202900 | 5,205381 | 5,207862 | 0,005 | 0,003403 | устойчиво |
8 | 5,23 | 5,203315 | 5,205796 | 5,208277 | 0,005 | 0,003818 | устойчиво |
9 | 6,10 | 5,203194 | 5,205675 | 5,208155 | 0,005 | 0,003697 | устойчиво |
10 | 6,28 | 5,203530 | 5,206011 | 5,208492 | 0,005 | 0,004033 | устойчиво |
11 | 7,08 | 5,202502 | 5,204983 | 5,207464 | 0,005 | 0,003005 | устойчиво |
12 | 7,18 | 5,203234 | 5,205715 | 5,208196 | 0,005 | 0,003737 | устойчиво |
13 | 8,14 | 5,206886 | 5,209367 | 5,211849 | 0,005 | 0,007390 | не устойчиво |
14 | 8,59 | 5,20400 | 5,20650 | 5,20890 | 0,005 | 0,004522 | устойчиво |