Содержание
ЛАБОРАТОРНАЯ РАБОТА № 1. 2
ЛАБОРАТОРНАЯ РАБОТА №3 5
ЛАБОРАТОРНАЯ РАБОТА №4 7
ЛАБОРАТОРНАЯ РАБОТА №5 9
ЛАБОРАТОРНАЯ РАБОТА № 1.
Тема: “Разработка математических моделей
и анализ динамики оптимизируемых систем”
Цель работы: Разработать математическую модель системы и проанализировать колебания оптимизируемой системы.
Задание:
-
Разработать программу и произвести численное интегрирование системы обыкновенных дифференциальных уравнений (1) на интервале времени при заданных начальных условиях и параметрах системы.
-
Исследовать свободные колебания системы, варьируя параметры () при заданных начальных условиях.
-
Исследовать вынужденные колебания системы, варьируя параметры () при различных частотах возмущающего воздействия .
Исходные данные
5000 |
Kg |
|
100000 |
N/m |
|
50000 |
N*s/m |
|
500 |
Kg |
|
1000000 |
N/m |
|
20000 |
N*s/m |
Исходные данные варьируются в соответствии с номером варианта (бригады) : .
Начальные значения переменных
0.1 |
М |
|
0 |
М/s |
|
0 |
M |
|
0 |
M/s |
Вычислить параметры, характеризующие колебания масс и , и вывести их на экран монитора в виде:
|
1/s |
|
|
Безр. |
|
|
1/s |
|
|
Безр. |
Параметры, характеризующие колебания масс и , рассчитываются по формулам:
; . . (3)
Параметры численного интегрирования
0 |
S |
|
50 |
S |
|
0.005 |
S |
Параметры гармонического возмущения
0.1 |
M |
|
20 |
1/s |
Вариант 1. Произвести вариацию параметра и оценить его влияние на частотные и диссипативные свойства динамической системы.
Выводы:
-
Параметр с1 характеризует жесткость первой пружины.
-
Чем больше значение параметра с1, тем медленнее система приходит в равновесное состояние.
ЛАБОРАТОРНАЯ РАБОТА № 2.
Тема: “ Моделирование внешней среды ”
Цель работы: Разработать программу для генерации случайного процесса и проанализировать корреляционную функцию и точность моделирования.
Задание:
1. Построить реализацию процесса заданной с корреляционной функцией на интервале времени с шагом и вывести ее график на экран монитора.
-
Вычислить для смоделированной реализации процесса оценку корреляционной функции и сравнить ее с исходной (заданной) корреляционной функцией . Оценить точность моделирования.
-
Произвести вариацию параметра и оценить его влияние на формируемый процесс.
Вариант 1. Произвести вариацию параметра и оценить его влияние на формируемый процесс.
Исходные данные
0.05 |
m |
|
0.15 |
1/m |
|
0.20 |
1/m |
|
0.3 |
1/m |
|
10.0 |
m/s |
|
0.1 |
Безр. |
|
50.0 |
s |
|
0.001 |
s |
|
1.5 |
s |
|
0.01 |
s |
Теоретические сведения:
При цифровом моделировании стохастических динамических систем возникает необходимость в формировании реализаций случайных возмущений по заданным статистическим характеристикам. На практике наиболее чаще всего встречаются нормальные случайные процессы. В этом случае корреляционная функция исчерпывающе описывает свойства случайного процесса.
Для формирования реализаций случайных процессов удобно использовать алгоритмы, основанные на линейном преобразовании последовательности [n] независимых чисел, распределенных по нормальному закону распределения, в последовательность q[n], коррелированную по определенному закону. Алгоритмы моделирования рассмотрены в работе [1, с.319-325].
Рассматривается нормально распределенный случайный процесс с корреляционной функцией:
,
где среднеквадратическое значение процесса;
,
где номер варианта;
;
; ; ;
, , и заданные константы.
Так, при алгоритм имеет вид:
q[n]=a0[n]+b1q[n-1],
где [n]-реализация независимых нормально распределенных чисел с параметрами m=0, =1.
Параметры алгоритма:
где h-заданный шаг дискретности аргумента; , , E – коэффициенты; среднеквадратичное значение формируемого процесса .
Необходимо вывести на экран монитора рассчитанные на интервале времени значения оценок математического ожидания и среднеквадратичного значения процесса . Сравнить их с заданными значениями и оценить точность.
Результат работы
Выводы: параметр с1 влияет на время кореляции, чем больше сила, тем больше время корелляции.
ЛАБОРАТОРНАЯ РАБОТА №3
Тема: «Статистический анализ систем методом Монте-Карло»
Цель работы: Разработать программу для статического анализа и получения зависимости критерия качества и ограничений.
Задание:
-
Произвести статистический анализ системы, используя разработанные ранее математическую модель системы (работа 1) и модель внешней среды (работа 2), варьируя указанный в задании параметр системы.
-
Построить зависимости критерия качества и ограничений от варьируемых параметров
Вариант 1. Произвести вариацию параметра .
Исходные данные
-0.1 |
M |
|
0.1 |
M |
|
|
M |
|
0.999 |
Безр. |
|
0.999 |
Безр. |
Формулы для вычисления оценки статистического критерия качества и ограничений второго рода
В качестве критерия принять среднеквадратичное значение ускорений массы в .
Ограничения вероятностей и невыбросов (невыходов) процесса , за уровни , а также процесса - за уровень являются ограничениями на функционалы от фазовых переменных системы (ограничения 2-го рода).
Ограничения 2-го рода задаются выражениями:
где заданные вероятности невыбросов;
;
;
заданные константы.
Вероятности невыбросов и вычисляются как математические ожидания соответствующих характеристических функций :
Результат работы
Выводы: параметр с1 влияет на критерии качества.
ЛАБОРАТОРНАЯ РАБОТА №4
Тема: «Параметрическая оптимизация»
Цель работы: Разработать программу для статического анализа и получения зависимости критерия качества и ограничений.
Задание
-
Для заданного метода минимизации функций многих переменных разработать алгоритм и программу минимизации заданного критерия качества с учетом ограничений первого и второго рода.
-
Протестировать алгоритм и программу на функции Розенброка. Оценить эффективность (скорость сходимости) алгоритма и программы.
-
Используя математическую модель системы (работа 1), модель внешней среды (работа 2), алгоритмы вычисления критерия и ограничений (работа 3) произвести оптимизацию указанных параметров системы (коэффициенты c1, c2, k1, k2) таким образом, чтобы критерий качества принимал наименьшее значение при удовлетворении ограничений первого и второго рода. На экран монитора вывести таблицу рассчитанных величин, отражающих динамику поиска.
-
Оценить влияние на точность оптимизации параметра, определяющего конец поиска (условия прекращения поиска).
-
Осуществить спуск из различных начальных точек. Сравнить результаты и сделать заключение о том, имеет ли задача локальные минимумы.
Замечания
-
Процедуру минимизации оформить в виде отдельного UNIT без использования формы.
-
Недостающие исходные данные принять такими же, как и в предыдущих работах.
В качестве критерия качества системы принять среднее квадратическое значение ускорения выходной переменной (среднеквадратическое значение ускорения массы , выраженное в долях ).
Вариант 1. Произвести оптимизацию по одному параметру , а затем по четырем параметрам:
Результат работы
Выводы: провели параметрическую оптимизацию.
ЛАБОРАТОРНАЯ РАБОТА №5
Тема: «Оптимальное проектирование бесконечномерных систем
(систем управления)»
Задание
Уравнения движения системы, линейной относительно управляющей функции, заданы в виде:
где заданные функции; константа:
Требуется найти управляющую функцию обеспечивающую для выходной переменной заданные желаемые свойства движения и кроме того, условия асимптотической устойчивости: при .
Вариант 1:
.
Результат работы
Вывод: провели оптимальное проектирование бесконечномерных систем.