- •Введение
- •1. Введение в системное моделирование
- •1.1. Понятие системы
- •1.2. Структура, функция и эффективность системы. Управление системой
- •1.3. Системный подход к моделированию
- •1.4. Системный характер технологических объектов
- •1.5. Действующий элемент системы
- •1.6. Системы автоматизированного моделирования
- •1.7. Экспертные системы
- •Контрольные вопросы и задания
- •2. Общие вопросы математического моделирования
- •2.1. Понятие моделирования. Математическая модель
- •2.2. Оптимальное моделирование
- •2.3. Некоторые типовые оптимизационные модели
- •Контрольные вопросы и задания
- •3. Нелинейные модели оптимизации
- •3.1. Градиентные методы
- •3.2. Общая задача нелинейного программирования. Постановка задачи
- •3.3. Градиентные методы
- •3.4. Случайный поиск с локальной оптимизацией
- •Контрольные вопросы и задания
- •Заключение
- •Библиографический список
- •394006 Воронеж, ул. 20-летия Октября, 84
1.5. Действующий элемент системы
Как было ранее отмечено, объект, выполняющий определенные системные функции и не подлежащий дальнейшему расчленению в рамках морфологического описания, принимается в качестве «первичного» элемента системы (рис. 1.9). Связь элемента с внешней средой, к которой относятся и другие элементы системы, осуществляется через его входы ͞x=(x1,x2,….,xm) и выходы ͞y=(y1,y2,….yn).Через ͞х элемент подвергается внешним воздействиям, а через ͞y сам воздействует на внешнюю среду.
Рис. 1.9. Схематичное представление элемента системы
Такое представление предполагает наличие у элемента по крайней мере одного входа и одного выхода.
Поведение элемента будем рассматривать как последовательную смену его «внутреннего» состояния в некотором интервале времени. Такая динамика (динамический процесс), как правило, описывается системами дифференциальных или конечно-разностных уравнений.
В общем случае описание можно представить как
ut’=f[͞u(t), ͞x(t),t] , ͞u(0)=͞u0;
y(t)=h[͞u(t), ͞x(t),t],
где (t) – р-мерный вектор, компоненты которого описывают состояние элемента в момент времени t; ͞u0 – начальное состояние элемента; у(t) – n-мерный вектор выходов элемента; ͞x(t) – m-мерный вектор входов элемента.
В дискретные моменты времени динамика может быть представлена как система разностных уравнений вида
где k = 0,1,2,3,,.. – порядковый номер, характеризующий последовательное дискретное время.
Заметим, что состояние элемента в каждый момент времени можно интерпретировать как некоторую точку в р-мерном пространстве состояний (фазовом пространстве). Тогда поведение элемента во времени можно представить в виде некоторой траектории (фазовой траектории) в пространстве состояний, которая описывается вектором функцией ͞u(t).
Изучая элемент с позиций описания взаимодействия вход-выход, можно считать некоторым преобразователем входов в выходы в соответствии с правилами, определенными внутренним описанием:
͞y (t)=R(t) (͞x(t)), (1)
где R(t) – символическое обозначение совокупности преобразования входа в каждый выход (оператор преобразования).
В частном случае при t= const (t→∞):
͞y(t)= R(x). (2)
Системный элемент, обладающий вышерассмотренными свойствами, будем называть действующим элементом системы.
Например, предположим, что оператор R – линейный (линейное преобразование). Соотношения (1)-(2) можно записать как ͞у =A*k, где , ; – матрица линейного преобразования:
или
Связь действующих элементов. Пусть Е1и Е2 – два действующих элемента с операторами R1 и R2 соответственно, ͞x1=(x11,x12,….,x1m1) и ͞x2=(x21,x22,….,x2m2) – их входные векторы, а ͞y1=(y11,y12,….,y1m1) и ͞y1=(y21,y22,….,y2m1) – выходные векторы.
Элемент E1 воздействует на элемент Е2 только таким образом, что Е2 через свои входы «принимает» значения всех или нескольких выходов элемента Е1.
Тогда можно записать:
x2j= y1i (3)
для некоторых значений i и j , где ; . Очевидно, что равенство (3) должно быть выполнено хотя бы для одной пары значений ( i ,j ) иначе не было бы воздействия элемента E1 на элемент E2.
Введем квадратную матрицу S21, состоящую из m2 или n1 строк и столбцов в зависимости от того, m2≥ n1 или n1 ≥ m2. Элементы Sij матрицы определяются следующим образом:
Другими словами, элементы матрицы имеют значение I, если соответствующая составляющая вектора ͞y элемента Е1 принимается в качестве составляющей вектора входа ͞х2 элемента Е2 , и 0, если этого не происходит. Матрица S21 примерно имеет вид:
S21 – нуль-единичная матрица, если хотя бы один ее элемент не является нулем.
Матрицу S21 назовем матрицей связи элемента Е2 с элементом Е1 (рис. 1.10). С ее помощью система равенств (3) записывается как ͞x2=s21y1.
Рис. 1.10. Схематическое представление связи двух действующих элементов
Матрицу связи элементов можно обобщить для любых элементов системы Er и Ek независимо от того, связан Er с Ek или нет. Если элемент Er не связан с Ek , то соответственно Skr= 0 (0 – нулевая матрица, все ее элементы равны 0). С другой стороны, если Skr имеет все ненулевые элементы, то это означает, что каждая координата вектора ͞хk связана со всеми координатами вектора ͞yr. В общем случае
͞xr=Skr͞yk. (4)
Допустим, что система состоит из N элементов E1, E2 ,…, Ern, которые описаны входными векторами и выходными векторами, соответственно ͞х1 , ͞x2, … ͞xn и y1,͞y2,…,͞yn.Тогда для каждой пары элементов можно записать векторное соотношение связи типа (4). В результате получается N (N-1) векторных равенств:
͞xr=Skr͞yk,
, , (5)
Сформированные таким образом матрицы связи образуют в совокупности квадратную матрицу размером N на N вида
Матрицу S назовем матрицей связей системы в целом. Она характеризует сеть связей между всеми элементами системы и практически выражает структуру.
Матрица S является нуль-единичной, так как все нуль-единичные матрицы связей Skr выступают в виде ее элементов.
Исходя из вышеизложенного, взаимодействие любых двух элемента Er и Ek можно представить в виде векторных равенств:
xr=SkrRk(͞xk),
͞yr=Rr(͞xr), , , (6)
где Rk и Rг – операторы преобразований для элементов Ек и Er соответственно.
Обобщая способ действия (6) на всю систему, получим:
X*=R*S(X); Y*=S*R(Y), (7)
где X, Y – состояние входов и выходов системы до преобразования; X*, Y* – новое состояние входов и выходов системы после преобразования; R – матрица операторов преобразования отдельных элементов; S – матрица связей между элементами.
Таким образом, соотношения (6) и (7) позволяют формально описывать действие системы в целом через действия отдельных ее элементов, образующих структуру систем.