Практическая работа №1

Иркутский государственный технический университет

НАДЕЖНОСТЬ ИНФОРМАЦИОННЫХ СИСТЕМ

Методические указания для студентов специальности

«Информационные системы в машиностроении»

ИРКУТСК

Г.

Надежность информационных систем. Методическое пособие по выполнению практических работ. Составил доцент Хрусталев Ю.П. - Иркутск, 2007

Рассмотрены методы анализа надежности информационных систем с учетом выхода из строя отдельных их элементов. Анализ основан на использовании результатов теории случайных процессов. Широко используется понятие «состояние системы». На основе графа состояний составляются системы дифференциальных уравнений – уравнения Колмогорова. Решения этих уравнений дают представления о динамике надежности систем и позволяют получить количественные оценки надежности. В пособии рассмотрены различные методы исследования надежности информационных систем:

- аналитические методы, основанные на применении преобразования Лапласа;

- метод численного решения систем ОДУ: метод Рунге – Кутты;

- методы имитационного моделирования.

На конкретных примерах показано, что все эти методы дают решения задач надежности, совпадающие с очень высокой точностью.

Предполагается, что студенты, изучающие данный предмет, прослушали курс «Моделирование информационных систем » и знакомы с методами имитационного моделирования таких систем.

Пособие состоит из 4 работ, содержит 28 страниц текста и содержит 10 рисунков. Список литературы содержит 4 наименования.

Практическая работа №1.

Решение задач на ЭВМ при наличии сбоев.

Потоки сбоев описываются потоками Пальма (потоками Эрланга порядка k). Если

точка Практическая работа №1 - student2.ru случайно попадает на интервал Практическая работа №1 - student2.ru между событиями потока, то закон распределения интервала Практическая работа №1 - student2.ru изменяется. При этом он по прежнему подчиняется закону Эрланга, но порядок становится равным k+1 [1].

Рассмотрим более подробно интервал Практическая работа №1 - student2.ru :

Практическая работа №1 - student2.ru

Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru

Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru

0 Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru

Рис.1 Схема попадания случайной точки τ на интервал T*

Практическая работа №1 - student2.ru – участок точки Практическая работа №1 - student2.ru от предыдущего события. Практическая работа №1 - student2.ru

Практическая работа №1 - student2.ru – интервал времени от точки до Практическая работа №1 - student2.ru последующего события (сбоя).

Если решение задачи начинается в случайный момент времени Практическая работа №1 - student2.ru , то задача будет решена, если время ее решения будет меньше, чем время, оставшееся до очередного сбоя

( т.е. меньше интервала R). Чтобы найти вероятность такого события - Р(А), необходимо знать плотность распределения случайной величины R.

В [1] показано, что

Практическая работа №1 - student2.ru (1)

где Практическая работа №1 - student2.ru Практическая работа №1 - student2.ru - плотность распределения случайной величины R,

F(t)= Практическая работа №1 - student2.ru - функция распределения случайной величины t;

k- порядок процесса Эрланга;

mt = Практическая работа №1 - student2.ru - математическое ожидание интервалов времени между событиями

потока;

Практическая работа №1 - student2.ru - интенсивность потока.

Таким образом, вероятность того, что задача, решение которой начато в момент t0, будет решена с первого раза, равна интегралу от плотности вероятности fR(t), вычисленному в пределах от t0 до бесконечности.

Практическая работа №1 - student2.ru (2)

Введём функцию: Практическая работа №1 - student2.ru (3)

Т.кПрактическая работа №1 - student2.ru ,[1]

то вероятность Р(А) можно достаточно просто вычислить с помощью функции R(m,a).

Эту задачу можно решить с помощью методов имитационного моделирования [2]. Программа, написанная в системе GPSS, состоит из двух сегментов. Первый сегмент имитирует процесс решения задачи. Транзакт вводится в систему в случайный момент времени (случайная величина выбирается из равномерного распределения: пусть для определенности в интервале, равном восьми часам ). Затем задача – транзакт – занимает компьютер, где и задерживается на время решения. После чего компьютер освобождается и задача покидает систему через оператор TERMINATE 1 с меткой GOOD (задача решена успешно).

Во втором сегменте имитируется поток сбоев. Генерация потоков Эрланга производится выбором из простейшего потока (интервалы между сбоями подчиняются экспоненциальному распределению) каждого k-го события. Данная процедура может быть реализована с помощью такой группы операторов:

num fvariable N$GO@3

GENERATE 120,FN$EXP

GO ADVANCE 1

TEST E V$num,0,TER

…. … …

TER TERMINATE

В данном фрагменте программы предполагается, что в исходном потоке сбои следуют в среднем через 120 интервалов времени. Интервалы между сбоями подчиняются экспоненциальному распределению (оператор GENERATE 120,FN$EXP). Оператор TEST пропускает далее только каждый третий транзакт (переменная num равна частному от деления по модулю на 3 числа транзактов, прошедших через метку GO). Транзакты , не кратные трем, покидают систему через оператор TERMINATE с меткой TER.

Таким образом формируется поток Эрланга 3-го порядка. Оператор ADVANCE необходим только для того, чтобы вставить метку GO (оператор GENERATE не может быть помеченным).

Транзакты, прошедшие блок TEST, захватывают компьютер, отсылая транзакт- задачу (если она в этот момент решалась на компьютере) на оператор с меткой AWAY:

AWAY TERMINATE 1

Многократно запустив программу с помощью оператора START (например, START 10000), необходимо подсчитать число транзактов, покинувших систему через метку GOOD, т. е. число успешных «прогонов» задачи, чтобы найти вероятность решения задачи с одного раза.

Пример решения задачи.

Интервал Практическая работа №1 - student2.ru между последовательными сбоями ЭВМ, устраняемыми практически мгновенно с помощью программных средств, имеет распределение Эрланга Практическая работа №1 - student2.ru порядка с параметром Практическая работа №1 - student2.ru (1/час). Для решения задачи требуется работа ЭВМ без сбоев в течение 2-х часов. Задачи начинают решать в произвольный момент Практическая работа №1 - student2.ru , никак не связанный с потоком сбоев.

Найти вероятность события:

A={задача будет решена с первого раза}.

Решение:

Событие Практическая работа №1 - student2.ru состоит в том, что с.в. Практическая работа №1 - student2.ru -время, оставшееся до очередного сбоя, принимает значение больше 2-х часов.

Практическая работа №1 - student2.ru

при Практическая работа №1 - student2.ru получаем:

Практическая работа №1 - student2.ru (4)

(т.к. Практическая работа №1 - student2.ru ) Практическая работа №1 - student2.ru

Практическая работа №1 - student2.ru

т.к. Практическая работа №1 - student2.ru

Практическая работа №1 - student2.ru

т.е. сделав замену переменных: Практическая работа №1 - student2.ru , получим:

Практическая работа №1 - student2.ru

(нижний предел в интеграле заменили: Практическая работа №1 - student2.ru )

Но Практическая работа №1 - student2.ru

Или можно Практическая работа №1 - student2.ru получить: Практическая работа №1 - student2.ru и т.д.

Практическая работа №1 - student2.ru

т.е.

Практическая работа №1 - student2.ru

Интеграл Практическая работа №1 - student2.ru можно вычислить непосредственно, т.е. без использования функций R(m,a) , используя возможности системы MathCAD.

Метод имитационного моделирования при достаточно большом числе прогонов(10000 и более) дает решение, совпадающее с аналитическим с точностью до одной тысячной.

Наши рекомендации