
- •1. Классификация моделей
- •4.Метод Эйлера
- •2.Дифференц. Ур-ие. Постановка задач.
- •5.Модифицированный метод Эйлера
- •6. Метод Рунге – Кутта
- •7. Метод Рунге-Кутта для систем дифференциальных уравнений.
- •11.Графы связей и структурные схемы динамических систем.
- •9. Выбор шага и погрешность решения
- •10. Жесткие задачи
- •11.Принципы Имитационное моделирование технических систем.
- •12. Сети связи и структруные схемы динамических систем.
- •13. Математическая модель гидравлического демпфера.
- •14. Математическая модель двигатель пост. Тока.
- •15. Преобразование передаточных функций звеньев в дифференциальные уравнения в форме Коши.
- •16. Синтез имитационной модели на основе структурной схемы.
- •17.Компонентные и топологические уравнения
- •18. Компонентные и топологические уравнения механической и электрической систем.
11.Графы связей и структурные схемы динамических систем.
В динамической системе, как правило, протекают процессы различной физической природы: механические, гидравлические, электрические, тепловые и т.д. Основой для построения структурной схемы такой системы служат физические законы. Они связывают между собой физические величины внутри системы и внешние силы. Основные уравнения динамической системы мы будем составлять с использованием понятия о динамических звеньях, изучаемых в ТАУ.
На первоначальном этапе составления структурной схемы связи между динамическими звеньями удобней представлять в виде графов. Подобно структурным схемам графы связей используются для наглядного изображения физических зависимостей в динамических системах.
Между структурной схемой и графом имеется прямое соотношение. Прямоугольник структурной схемы соответствует ребру, а линия передачи сигнала- вершине графа. Обычно при моделировании используют алгебраизированные графы связей. Сравним алгебраизированный граф связи и структурную схему одной и той же системы:
Правило составления графа связи:
Номер ребра соответствует номеру передаточной функции звена.
Ребро, соединяющее одноименные величины не имеет номера. Передаточная функция такого ребра равна 1.
Элементу сравнения на графе соответствует разность, а сумматору- сумма двух величин, равная третьей величине результату.
Результат вычитания или суммирования всегда записывают справа от знака равенства.
9. Выбор шага и погрешность решения
Одним из важнейших практических вопросов, которые встают
перед инженером, составляющим программы решения
дифференциальных уравнений, является выбор подходящей
величины шага. Если шаг слишком мал, то расчёт потребует
неоправданно много машинного времени. Если же наоборот - шаг выбран слишком большим, то
значительной будет погрешность, обусловленная усечением рядов, и
накопившаяся
суммарная ошибка будет также недопустимо
большой
Кроме
того, при больших значениях h
возникает
опасность
появления цифровой неустойчивости решения. Обычно, выбирая величину шага, стремятся, чтобы локальная
ошибка на шаге была меньше некоторой заданной допустимой
величины. Вообще говоря, если порядок точности метода Р, то
локальная
ошибка определяется выражением:
где
c некоторая
постоянная; h
- шаг
интегрирования.
Указанное
выражение с помощью О-символики
может быть записано более компактно
Рассмотрим
теперь глобальную ошибку дискретизации
в конечной точке t=tN.
По мере повышения требований к точности,
длина шага h будет убывать, а общее из числа N, необходимое для
достижения
tN,
будет возрастать:
Далее,
глобальная ошибка E(h)
может быть
представлена как
сумма
N локальных
ошибок и поэтому мы можем, округляя,
записать:
Для
метода Эйлера p
=1, так что
уменьшение длины шага в 2
раза уменьшает среднюю локальную ошибку примерно в 2p+1 = 4
раза. Но так как для достижения tN теперь потребуется вдвое больше
шагов, то глобальная ошибка уменьшится лишь в 2p = 2 раза.
Если используется один из методов прогноза и коррекции, то
ошибка на шаге определяется величиной последнего члена в формуле
коррекции.
При использовании же одношаговых методов Рунге - Кутта,
локальную ошибку не удается выразить в явной форме. Один из
методов оценки этой ошибки основан на двойном счете. Недостатком этого метода является то, что значение yn+1
приходится вычислять дважды, причем второй раз с помощью двух
шагов. С
другой стороны, главным достоинством методов прогноза и
коррекции является простота оценки ошибки на шаге.