- •1 Выделение подсистем на основе некоторой меры
- •1.1. Построить граф системы
- •1.3. Определитель системы
- •1.4. Определить передаточную функцию из узла X в узел y по формуле Мезона.
- •1.5 Расчет изменения энтропии системы и вероятности нахождения в каждом из узлов.
- •2. Выбор типов шкал
- •3. Определение порядка проведения работ
- •4. Построение моделей систем
- •5 Анализ иерархий
- •5.1. Составление матрицы критериев
- •Список используемой литературы.
- •Приложение 1 – Функциональная модель системы idef0
- •Приложение 2 - Информационная модель системы idef1x .
- •Приложение 3 – Динамическая модель системы в виде сети Петри
1.3. Определитель системы
Δ=
∆=1-(K1+K2+K3+K4+ K5)+(K2*K3+K2*K4+ K3*K4+ K3*K5)-(K2*K3*K4) =
=1-((-0.004-0.135-0.72-0.72-0.017)+(-0.135*(-0.72)+(-0.135)*(-0.72)+(-0.72)*
*(-0.72)-(-0.135)*(-0.72)*(-0.72)) = 1.8142
∆= 1.8142
1.4. Определить передаточную функцию из узла X в узел y по формуле Мезона.
Рассчитаем веса контуров:
K1=R12*R24*R45*R56*R67*R71=0.8*0.1*0.8*0.2*(-0.9)*0.4=-0.004 K2=R67*R73*R36=(-0.9)*0.5*0.3=-0.135
K3= R45*R54=0.8*(-0.9)=-0.72
K4=R12*R21=0.8*(-0.9)=-0.72
K5=R12*R23*R36*R67*R71=0.8*0.2*0.3*(-0.9)*0.4=-0.017
Определим пути (x=1, y=6):
P1=R12*R23*R36 =0.8*0.2*0.3=0.048;
P2= R12*R24*R45*R56=0.8*0.1*0.8*0.2=0.0128
Таблица 5 - Матрица путей
|
V1 |
V2 |
V3 |
V4 |
V5 |
V6 |
V7 |
P1 |
1 |
1 |
1 |
|
|
1 |
|
P2 |
1 |
1 |
|
1 |
1 |
1 |
|
Для того чтобы установить эквивалентный оператор учитывающий связь от i-го узла к j-му узлу используем формулу Мезона.
минор
i-го
пути
Pi = вес i-го путь от входа к выходу
=Р1=1-К3
=1+0,72=1,72
=Р2=1
=
1.5 Расчет изменения энтропии системы и вероятности нахождения в каждом из узлов.
Используя граф системы (рисунок 1) и данные таблицы 1, рассчитаем изменение энтропии системы.
Начальная энтропия:
Р7=Р6=Р5=Р4=Р3=Р2=Р1=1/7;
Hнач.= -7(1/7log21/7)=2,807
Формула Колмогорова для скорости изменения вероятности нахождения системы в том или ином состоянии:
dPi/dt=Σλj*Pj–Σλi*Pi,
где λ – интенсивность перехода их одного состояния в другое (равна массе соответствующих дуг).
P1:
=0,4P7-0,8*P1=0
P2: =0,8*P1- 0,2*P2-0,1*P2=0
P3: =0,5*P7+0,2*P2-0,3*P3=0
P4: =0,1*P2-0,9*P5-0,8* P4=0
P5: =0,8*P4-0,2*P5=0
P6:
=0,2*
P5+0,3*Р3+0,9*Р6,
так как неизвестных 7 а уравнений 8
вычеркиваем одно из уравнений.
P7: =-0,9*P6-0,5*P7-0,4* P7=0
∑ Рi=1
Решаем с помощью пакета прикладных программ для решения задач технических в SCILAB-5.4.0 3 способами.
1 способ: при помощи символа \ :
2 способ: методом обратной матрицы:
3 способ - метод Гаусса:
Отсюда получаем:
Р1=0.0812903,
Р2=0.2167742,
Р3=0.4154839,
Р4=-0.0077419,
Р5=-0.0309677,
Р6=0.1625806,
Р7=0.1625806.
Hуст(s)= - ∑(Pi*log2Pi)= -(0.0812903*log20.0812903+0.2167742*
*log20.2167742+0.4154839*log20.4154839+ (-0.0077419)*log20.0077419+
+ (-0.0309677)*log20.0309677+0.1625806*log20.1625806+0.1625806*
*log20.1625806=1.9416
H(s)=2,807-1.9416=0.8654>0.
Суммарная энтропия элементов установившегося состояния равна 1.9416, соответственно, изменение энтропии составило0.8654, и так как энтропия системы больше нуля, то система стремится к порядку.
Вывод:
С помощью матриц смежности, инцидентности и контуров была определена структура графа.
Были определены все элементарные пути из узла x в узел y, также определили вес пути от начала x до конца y.
С помощью формулы Мезона была установлена связь от i-го узла к j-му.
Используя веса контуров, был рассчитан определитель графа.
С помощью расчета энтропии методом Колмогорова была определена энтропия системы.
