Курсовая работа: Теория распределения информации

Курсовая работа

Министерство науки и высшего образования Республики Казахстан

Алматинский институт энергетики и связи

Кафедра Автоматической электросвязи

г. Алматы, 1999 г.

Задание 1.

Построить огибающую распределения вероятности занятия линии в пучке из V, на каждую из которых поступает интенсивность нагрузкиа при условии, что:

а) N >> V; б) N V; в) N, V

Для каждого используемого распределения рассчитать среднее число занятых линий и их дисперсию.

Для расчета число линий в пучке определить из следующего выражения:

V= ;

целая часть полученного числа, где NN– номер варианта.

Средняя интенсивность нагрузки, поступающей на одну линию:

а = 0,2+0,01 * NN

Примечания:

Для огибающей распределения привести таблицу в виде:

В распределении Пуассона привести шесть – восемь составляющих, включая значение вероятности для i = (целая часть А)

А = а * V

Решение:

Случайной называют такую величину, которая в результате эксперимента принимает какое то определенное значение, заранее не известное и зависящее от случайных причин, которые наперед предугадать невозможно. Различают дискретные и непрерывные случайные величины. Дискретная случайная величина определяется распределением вероятностей, непрерывная случайная величина – функцией распределения основными характеристиками случайной величины являются математическое ожидание и дисперсия.

Определим исходные данные для расчета:

V=

a = 0.2 + 0.01 * 11 = 0.31 Эрл (средняя интенсивность нагрузки)

А = а * V = 0,31 * 11 = 3,41 » 4 Эрл (нагрузка)

а) Определим вероятности занятия линий в пучке из V = 11, при условии N >> V (N– число источников нагрузки).

Для этого используем распределение Эрланга, представляющее собой усеченное распределение Пуассона, в котором взяты первые V+1 значения и пронумерованы так, чтобы сумма вероятностей была равна единице.

Распределение Эрланга имеет вид:

P_i (V) = , ,

где P_i (V) – вероятность занятия любых i линий в пучке из V.

Для определения составляющих распределения Эрланга можно применить следующее реккурентное соотношение:

Математическое ожидание и дисперсия числа занятых линий соответственно равны:

где P_v –вероятность занятости всех линий в пучке из V.

Произведем расчет:

Р₀ =

Р₁ = Р₀ * = 0,072 Р₂ = Р₁ * = 0,144

Р₃ = Р₂ * = 0,192 Р₄ = Р₃ * = 0,192

Р₅ = Р₄ * = 0,153 Р₆ = Р₅ * = 0,102

Р₇ = Р₆ * = 0,058 Р₈ = Р₇ * = 0,029

Р₉ = Р₈ * = 0,012 Р₁₀ = Р₉ * = 4,8 * 10^-3

Р₁₁ = Р₁₀ * = 1,7 * 10^-3

M( i ) = 4 * (1 - 1,7 * 10^-3 ) = 3,99

D( i ) = 3,99 – 4 * 1,7 * 10^-3 * (11 – 3,99) = 3,94

Данные результаты вычислений сведем в таблицу 1:

Таблица1

б) Определим вероятность занятия линий в пучке из V=11, при условии N@V. Применим распределение Бернулли (биноминальное распределение), которое имеет вид:

где: P_i (V) – вероятность занятия любых i линий в пучке из V;

- число сочетаний из V по i (i = 0, V)

,

а – средняя интенсивность поступающей нагрузки на одну линию

V-линейного пучка от N источников.

Для вычисления вероятностей можно воспользоваться следующей рекурентной формулой:

Математическое ожидание и дисперсия числа занятых линий соответственно равны:

M( i ) = V*a; D( i ) = V * a * (1-a)

Произведем расчет:

;

Р₁ = 16,8*10^-3 *

Р₂ = 16,8*10^-3 *

Р₃ = 16,8*10^-3 *

Р₄ = 16,8*10^-3 *

Р₅ = 16,8*10^-3 *

Р₆ = 16,8*10^-3 *

Р₇ = 16,8*10^-3 *

Р₈ = 16,8*10^-3 *

Р₉ = 16,8*10^-3 *

Р₁₀ = 16,8*10^-3 *

Р₁₁ = 16,8*10^-3 *

M( i ) = 11 * 0,31 = 3,41; D( i ) = 11 * 0,31 * (1 – 0,31) = 2,35

Результаты вычислений сведем в таблицу 2:

Таблица2

в) Определим вероятность занятия линий в пучке из V=11 , при условии N,V®¥.

Используем распределение Пуассона, как вероятность занятия i линий в бесконечном пучке линий за промежуток времени t:

, ,

где: l - параметр потока, выз/час

lt – средняя интенсивность нагрузки поступающей на пучок линий (А=lt).

Легко показать, что:

,

Произведем расчет:

Р₀ = * е^-4 = 0,018 Р₁ = 0,018 * = 0,036

Р₄ = * 0,018 = 0,192 Р₆ = 0,018 * = 0,102

Р₈ = 0,018 * = 0,029 Р₁₀ = 0,018 * = 0,0052

Р₁₂ = 0,018 * = 0,0006

M( i ) = D( i ) = 4

Результаты вычислений сведем в таблицу 3:

Таблица3

По данным таблиц 1, 2, 3 построим графики огибающей вероятности для трех случаев: а) N>>V, б) N@V, в) N, V ®¥ ; рис. 1.

Задание 2.

На коммутационную систему поступает простейший поток вызовов с интенсивностью А.

Рассчитать вероятность поступления не менее к вызовов за промежуток времени [ 0, t^* ]:

Р_к (t^* ), где t^* = 0,5; 1,0; 1,5; 2,0

Построить функцию распределения промежутков времени между двумя последовательными моментами поступления вызовов:

F(t^* ), t^* = 0; 0,1; 0,2; …

Рассчитать вероятность поступления не менее к вызовов за интервал времени [ 0, t^* ]:

P_{i ³ k} (t^* ), где t^* = 1

Примечание: 1. Для расчета значений A и V взять из задания 1.

2.Число вызовов к определить из выражения: к = [V/2] - целая часть числа.

Для построения графика взять не менее пяти значений F(t^* ). Результаты привести в виде таблицы:

Расчет P_{i ³ k} (t^* ) провести не менее чем для восьми членов суммы.

Решение:

Потоком вызовов называют последовательность однородных событий, поступающих через случайные интервалы времени. Поток вызовов может быть задан тремя эквивалентными способами:

Вероятностью поступления к вызовов за интервал времени [0,t).

Функцией распределения промежутков времени между двумя последовательными моментами поступления вызовов.

Вероятность поступления не менее к вызовов за интервал времени [0,t).

Свойства потоков: станционарность, ординарность и полное или частичное отсутствие последействия. Потоки классифицируются с точки зрения наличия или отсутствия этих свойств.

Основными характеристиками потоков вызовов являются: интенсивность m и параметр l.

Простейшим потоком называется ординарный стационарный поток без последействия.

Рассчитаем вероятность поступления не менее к вызовов за интервал времени [0,t).

,

где: к = 0, 1, …;

t^* = t /`t ; где `t – средняя длительность обслуживания вызова.

Определим данные для расчетов:

К = 11/2 = 6; А = 4; V = 11;

Производим расчеты для t^* = 0,5 с.

P₂ (0,5) = 0,13 P₃ (0,5) = 0,18 P₄ (0,5) = 0,09

P₅ (0,5) = 0,03 P₆ (0,5) = 0,012

Производим расчеты для t^* = 1,0 с.

P₂ (1) = 0,14 P₃ (1) = 0,19 P₄ (1) = 0,19

P₅ (1) = 0,15 P₆ (1) = 0,1

Производим расчеты для t^* = 1,5 с.

P₂ (1,5) = 0,044 P₃ (1,5) = 0,089 P₄ (1,5) = 0,13

P₅ (1,5) = 0,16 P₆ (1,5) = 0,16

Производим расчеты для t^* = 2 с.

P₂ (2) = 0,01 P₃ (2) = 0,028 P₄ (2) = 0,057

P₅ (2) = 0,91 P₆ (2) = 0,122

Рассчитаем функцию распределения промежутков времени между двумя последовательными моментами поступления вызовов:

где Z_k – промежуток времени между ( к-1 )-м и к-м вызовами.

F(0) = 1 – e^{-4 * 0} = 0 F(0,1) = 1 – e^{-4 * 0,1} = 0,32 F(0,2) = 1 – e^{-4 * 0,2} = 0,55

F(0,3) = 0,69 F(0,4) = 0,79 F(0,5) = 0,86

F(0,6) = 0,9 F(0,7) = 0,93

Результаты вычислений занесем в таблицу 4:

Таблица 4

Рассчитаем вероятность поступления не менее к вызовов за промежуток времени [0, t^* ):

, при t^* =1.

P_{6 ³ 6} (1) = 1 – 0,84 = 0,16 P_{10 ³ 6} (1) = 1 – 0,005 = 0,995

P_{7 ³ 6} (1) = 1 – 0,05 = 0,95 P_{11 ³ 6} (1) = 1 – 0,001 = 0,999

P_{8 ³ 6} (1) = 1 – 0,02 = 0,98 P_{12 ³ 6} (1) = 1 – 0,0006 = 0,9994

P_{9 ³ 6} (1) = 1 – 0,013 = 0,987 P_{13 ³ 6} (1) = 1 – 0,0001 = 0,9999

Интенсивность простейшего потока вызовов m численно равна параметру l, а при t = `t =1: m = l = А = 4.

Задание 3.

Рассчитать интенсивность поступающей нагрузки на входы I ГИ для АТСКУ – А _{вх. I ГИ} .

Рассчитать средние интенсивности удельных абонентских нагрузок для абонентских лини народно-хозяйственного и квартирного секторов : А_НХ и А_КВ , а так же среднюю удельную интенсивность нагрузки на абонентскую линию АТС - А_{ИСХ .}

Пересчитать интенсивность нагрузки на выход ступени I ГИ.

Исходные данные, таблица 5:

Таблица 5

Решение:

1. Основными параметрами интенсивности нагрузки являются:

N_i – число источников нагрузки i-й категории.

C_i – среднее число вызовов, поступающих от одного источника i-й категории в ЧНН (час наибольшей нагрузки).

t_i – средняя длительность одного занятия для вызова от источника i-й категории.

Различают следующие категории источников нагрузки: абонентские линии народнохозяйственного сектора (НХ), абонентские линии квартирного сектора индивидуального пользования (кв.и.), абонентские линии квартирного сектора коллективного сектора (кв.к.), таксофоны (т). Для расчета используем две категории: абонентские линии народнохозяйственного сектора (НХ) и абонентские линии квартирного сектора (кв).

Интенсивность поступающей нагрузки:

,

Средняя длительность одного занятия зависит от типа системы коммутации и определяется выражением:

где: Р_р – доля вызовов из общего числа, для которых соединения закончились разговором; Р_з – доля вызовов из общего числа, для которых соединения не закончились разговором из-за занятости линии вызываемого абонента; Р_но –то же из за неответа вызываемого абонента; Р_ош – то же из-за ошибок в наборе номера; Р_техн - то же из-за технических неисправностей в узлах коммутации (при расчетах Р_техн =0); t_{р i} , t_з , t_но , t_ош , t_техн – средние длительности занятий соответствующие этим случаям. Их можно определить из следующих выражений:

t_Pi = t_y + t_пв + T_i + t₀

t_з = t_y + t_сз + t₀

t_но = t_y + t_пвн + t₀

t_ош = 18 с.

где: t_у – средняя длительность установления соединения; t_пв и t_пвн средняя длительность слушания сигнала «КПВ» (t_пв =7 с. в случае разговора между абонентами; t_пвн =30 с. в случае неответа вызываемого абонента);

T_i – продолжительность разговора для вызова i-й категории;

t_о – продолжительность отбоя;

t_сз – продолжительность слушания сигнала “Занято”

t_у = 0,5*t_МАВИ + t_МРИ + t_МРИ + t_СО + n*t_{Н +} t_{I ГИ} + t_{М I ГИ} + t_{МС D} + t_{МС D}

где t_j – время ожидания обслуживания маркером j-й ступени; t_j = 0,1 с.

t_МАВИ – время установления соединения маркером АВ на ступени АИ при исходящей связи; t_МАВИ = 0,3 с.

t_МРИ - время установления соединения маркером ступени РИ; t_МРИ = 0,2 с.

t_{М I ГИ} - время установления соединения маркером ступени IГИ; t_{М I ГИ} = 0,65 с.

t_{МС D} - время установления соединения маркером CD; t_{МС D} = 1 С.

t_СО – средняя длительность слушания сигнала «Ответ станции»; t_СО = 3 с.

t_Н – средняя длительность набора одного знака номера; t_Н = 1,5 с.

n – значность номера.

Значения t_о и t_сз для АТСКУ следующие: t_сз = 0,6 с., t_о = 0.

Р_Р = 0,6; Р_з = 0,2; Р_но = 0,15; Р_ош = 0,05;

t_у = 0,5 * 0,3 + 0,1 + 0,2 + 3 + 5 * 1,5 + 0,1 + 0,65 + 0,1 + 1 = 12,8 с.

t_рнх = 12.8 + 7 + 100 + 0.6 = 120,4 с.

t_ркв = 12,8 + 7 + 130 + 0,6 = 150,4 с.

Р_Р * t_рнх = 0,6 * 120,4 = 72,24

Р_Р * t_ркв = 0,6 * 150,4 = 90,24

t_з = t_у + t_сз + t_о = 12,8+0+0,6 = 13,4 с.

Р_з * t_з = 0,2*13,4 = 2,68

t_но = t_у + t_пвн + t_о = 12,8+30+0,6 = 43,4 с.

Р_но * t_но =0,15*43,4 = 6,51

Р_ош * t_ош = 0,05*18 = 0,9

t_нх = 72,24+2,68+6,51+0,9+0 = 82,33 с.

t_кв = 90,24+2,68+6,51+0,9+0 = 100,33 с.

А_{ВХ I ГИНХ} = = 434,5 Эрл

А_{ВХ I ГИКВ} = = 167,2 Эрл

А_{ВХ I ГИ} = 434,5 + 167,2 = 601,7 Эрл

2. Рассчитаем средние интенсивности удельных абонентских нагрузок для абонентских линий народнохозяйственного и квартирного секторов:

, Эрл

, Эрл

Средняя удельная интенсивность нагрузки на абонентскую линию АТС:

, Эрл

А_НХ = = 0,087 Эрл А_КВ = = 0,042 Эрл

А_ИСХ = = 0,07 Эрл

3. Пересчитаем нагрузку со входа ступени I ГИ на ее выход:

,

где t_{вх I ГИ} и t_{вых I ГИ} – соответственно среднее время занятия входа ступени I ГИ и среднее время занятия выхода ступени I ГИ:

t_{вых I ГИ} = t_{вх I ГИ} - Dt,

где Dt – разница между временами занятия на входе и выходе ступени I ГИ. Для АТСКУ:

Dt = 0,5*t_МАВИ + t_МРИ + t_МРИ + t_СО + n*t_{Н +} t_{М I ГИ} + t_{М I ГИ}

t_{ВХ I ГИ} = А_{ВХ I ГИ} / N_нх * С_нх + N_кв * С_кв

Dt = 0,5 * 0,3 + 0,1 + 0,2 + 3 + 5 * 1,5 + 0,1 + 0,65 = 11,7 с.

t_{ВХ I ГИ} = = 86,6 с.

t_{ВЫХ I ГИ} = t_{ВХ I ГИ} - Dt = 86,6 – 11,7 = 74,9 с.

А_{ВЫХ I ГИ} = 74,9/86,6 * 601,7 = 520,4 Эрл

Задание 4.

Рассчитать и построить зависимость числа линий V и коэффициента использования h (пропускная способность) от величины интенсивности нагрузки при величине потерь Р = 0,0N_В , где N_В – номер варианта.

Результаты расчета представить в виде таблицы при Р = const (постоянная).

Решение:

Вероятность занятия любых i линий в полнодоступном пучке из V при обслуживании простейшего потока вызовов определяется распределением Эрланга:

Различают следующие виды потерь: потери от времени P_t , потери по вызовам P_в , потери по нагрузке P_н . Потери по времени P_t - доля времени, в течение которого заняты все V линии пучка. Потери по вызовам определяются отношением числа потерянных вызовов С_пот к числу поступивших С_пост :

P_в = С_пот / С_пост

Потери по нагпрузке определяются отношением интенсивности потерянной нагрузки Y_пот к интенсивности поступившей А :

P_н = Y_пот / А

При обслуживании простейшего потока вызовов перечисленные выше три вида потерь совпадают P_t = P_в = P_н и равны вероятности занятия V линий в пучке:

Р_V = P_t = P_в = P_н = E_V,V (A) =

Обслуженной нагрузкой называют нагрузку на выходе коммутационной схемы, ее интенсивность определяют из выражения:

Y = F - Y_ПОТ = A * (1 - E_V (A))

Среднее использование одной линии в пучке равно:

h = Y / V

При Р = 0,011 (11 вариант), по известным А, используя таблицы вероятности потерь определим соответствующие V и рассчитаем для каждого значения А интенсивность Y и среднее использование h.

А = 1, Эрл V₁ =5 Y₁ =1(1-0,011) = 0,989 h = 0,197

А = 3, Эрл V₃ =8 Y₃ =3(1-0,011) = 2,96 h = 0,986

А = 5, Эрл V₅ =11 Y₅ =5(1-0,011) = 4,94 h = 0,449

А = 10, Эрл V₁₀ =18 Y₁₀ =10(1-0,011) = 9,89 h = 0,549

А = 15, Эрл V₁₅ =24 Y₁₅ =15(1-0,011) = 14,83 h = 0,617

А = 20, Эрл V₂₀ =30 Y₂₀ =20(1-0,011) = 19,78 h = 0,659

А = 25, Эрл V₂₅ =36 Y₂₅ =25(1-0,011) = 24,73 h = 0,686

А = 30, Эрл V₃₀ =42 Y₃₀ =30(1-0,011) = 29,67 h = 0,706

А = 40, Эрл V₄₀ =53 Y₄₀ =40(1-0,011) = 39,56 h = 0,746

А = 50, Эрл V₅₀ =64 Y₅₀ =50(1-0,011) = 49,45 h = 0,772

Результаты расчетов занесем в таблицу 6:

Таблица 6

Построим график зависимости числа линий V и коэффициента использования h от величины интенсивности нагрузки Y при величине Р=0,011.

Задание 5.

1. Построить оптимальную равномерную неполнодоступную (НПД) схему, имеющую следующие параметры: V – емкость пучка, g – число нагрузочных групп, d – доступность. Привести матрицу связности.

Исходные данные:

V = 25*N_гр + N_В

D = 10*N_гр

где N_гр – номер группы , N_В – номер варианта.

8, если N₈ =1-10;

g = 10, если N₈ =11-21

12, если N₈ =21-…

2. Рассчитать и построить зависимость числа линий V от величины потерь Р неполнодоступного пучка при значении A и D=10 по формуле Эрланга, О Делла, Пальма-Якобеуса. Результаты привести в виде таблицы и графика:

*- Модифицированная формула Пальма-Якобеуса.

Исходные данные: А – поступающая нагрузка взять в задании 1.

Решение:

Неполнодоступное включение это когда входу доступны не все, а часть выходов (d-определяет количество доступных выходов, d<V). Главная особенность НПД схем в том, что при одних и тех же параметрах можно построить множество различных схем, отличающихся пропускной способностью. Основными параметрами схемы являются: g – число нагрузочных групп, d – доступность, V – количество подключаемых к выходам соединительных устройств. Нагрузочной группой называется совокупность источников вызовов, обслуживаемых одними и теми же d-соединительными устройствами в НПД схеме. НПД схемы бывают трех видов ступенчатая, равномерная и идеально-симметричная. По типу соединений: прямое, перехваченное и со сдвигом. При прямом включении объединяются одноименные выходы соседних нагрузочных групп. При перехваченном включении выходы каждой нагрузочной группы соединяются по возможности равномерно с одноименными выходами остальных нагрузочных групп. При включении со сдвигом выходы одной нагрузочной группы соединяются с разноименными выходами других нагрузочных групп.

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

Цилиндр – это циклосхема, у которой обязательно равенство V=g (число выходов совпадает с числом нагрузочных групп). Размер цилиндра d представляет собой число охватываемых выходов каждой нагрузочной группы. Цилиндр размера d называется d-шаговым. Кроме размера цилиндр характеризуется наклоном.

Для построения оптимальной схемы нужно построить матрицу связности. Матрица связности – квадратная (g,g), симметричная относительно главной диагонали (по диагонали стоит d доступность), элементы матрицы связности показывают число связей между нагрузочными группами. Для оптимальности схемы необходимо чтобы матрицы связности были однородными и не отличались не более чем на единицу.

1.

V = 25*1+11 = 36

D = 10*1 = 10

G = 10

1) Определим размер цилиндров:

r = [(g*d)/V] (целая часть)

r = [(10*10)/36] = 2

2) Наша схема будет состоять из r и r+1 шаговых цилиндров

r+1 = 2 + 1 = 3

3) Определяем общее количество цилиндров:

k » V / g k » 36 / 10 » 4

4) Определим количество двух шаговых цилиндров:

5) Определим количество трех шаговых цилиндров:

k_r+1 = k – k_r

k_r+1 = 4 – 1 = 3

6) Определим наклон цилиндров. Для этого строим матрицу связности (табл. 7):

Таблица 7

7) Построим схемы цилиндров:

1 2 3 4 5 6 7 8 9 10

I

II

11 12 13 14 15 16 17 18 19 20

III

IV

V

21 22 23 24 25 26 27 28 29 30

VI

VII

VIII

31 32 33 34 35 36

XIX

X

XI

2. Для практических расчетов пропускной способности однозвенных НПД коммутационных схем используют приближенные методы.

Упрощенная формула Эрланга:

где У₀ – интенсивность обслуженной нагрузки пучком линий;

Р – вероятность потерь;

средняя пропускная способность одной линии пучка.

Формула О¢ Делла:

где У_D – нагрузка, обслуженная полнодоступным пучком из d линий при потерях и приблизительно определяемая с помощью 1-й формулы Эрланга.

Формула Пальма-Якобеуса:

где А – интенсивность поступающей нагрузки на пучок линий.

В модифицированной формуле Пальма-Якобеуса вместо поступающей нагрузки А в формулу Пальма-Якобеуса подставляется значение фиктивной нагрузки А_ф определяемой из выражения:

А_ф = Y / (1 - E_V (А_ф ))

P = E_V (А_ф ) / (E_{V - d} (А_ф ))

где Y = А(1-Р)

Рассчитаем по формуле Эрланга:

Р = 0,001

У_О = А(1-Р) = 4(1-0,001) = 3,996

V=3,996 / = 7,99 » 8

Р = 0,002

У_О = 3,992 V = 7,43 » 8

Р = 0,003

У_О = 3,988 V = 7,12 » 8

Рассчитаем по формуле О¢ Делла:

Р = 0,001

У_О = 3,996 У₁₀ = 3,089

V = 10 + = 15,79 » 16

Р = 0,002

У_О = 3,992 У₁₀ = 3,420 V = 14,78 » 15

Р = 0,003

У_О = 3,988 У₁₀ = 3,637 V = 14,1 » 15

Список литературы

Корнышев Ю.Н., Фань Г.Л. «Теория распределения информации». М., Радио и связь, 1985 г.

Башарин Г.Л. Таблицы вероятностей и средних квадратичных отклонений потерь на полнодоступном пучке линий. М., 1962 г.

Ионин Г.Л., Седол Я.Я. Таблицы вероятностных характеристик полнодоступного пучка при повторных вызовах. М., Наука, 1970 г.

Айтуова Р.Ч., Туманбаева К.Х. Методические указания к выполнению курсовой работы. Алматы, АИЭС, 1998 г.