В последнее десятилетие основным направлением в области развития телекоммуникационных сетей является создание интегрированной универсальной мультисервисной сети (МС), объединяющей различные виды связи на основе единых организационных и технологических принципов. Решением этой проблемы занимаются во многих информационно развитых странах, а также в нашей стране. Такая сеть предоставляет пользователям возможность мультисервисного обслуживания, т. е. возможность передавать, принимать и обрабатывать в цифровом виде различную по характеру и объему информацию.
Традиционные информационные потоки, порождаемые речевыми сообщениями, дополняются потоками факсимильных сообщений, звуковой почты, мобильной' связи, верификации кредитных" карточек, соединений с сетью Интернет и т. д. '.
В соответствии с принятой классификацией сетевые приложения мультисервисных сетей связи можно разбить на три основные группы: передача данных, пакетная телефония и потоковое видео. Нагрузка коммуникационных приложений, относящихся к первой группе, как правило, передается по принципу Best Effort. Она не чувствительна, если ее величина лежит в разумных пределах.
Нагрузка коммуникационных приложений второй и третьей групп принадлежит к категории мультимедийных нагрузок (Stream Traffic). Для ее передачи необходимо выделить гарантированную полосу пропускания и обеспечить характеристики качества обслуживания (QoS).
Если в начальный период использования сети Интернет главным достоинством пакетной передачи информации была возможность создавать надежные сети, способные передавать нагрузку на большие расстояния, то сейчас на первый план выходит способность современных пакетных технологий обеспечить заданное QoS.
К технологиям, обладающим данными характеристиками, необходимо, в первую очередь, отнести Frame Relay, ATM, а также MPLS.
Применение технологий Frame Relay, ATM, MPLS дает возможность внести в архитектуру протокола IP механизм образования виртуальных путей, позволяющих рассматривать процесс их предоставления для поступающих потоков, анализируемых на уровне соединения, аналогично процессу, имеющему место при занятии маршрута в сетях с коммутацией каналов.
Это означает, что модели и методы оценки показателей передачи нагрузки, развитые в теории телетрафика при анализе классических систем связи, можно переносить на модели, появляющиеся при описании процесса передачи мультимедийной нагрузки современных пакетных сетей.
В настоящее время известен ряд работ, касающихся анализа и синтеза мультисервисных сетей. В частности, можно отметить работы В. А. Ершова, Э. Б. Ершовой, В. С. Лагутина, С. Н. Степанова, Н. А. Кузнецова, A.A.Fredericks, K.W.Ross и других.
Вместе с тем, вопросы расчета пропускной' способности участков мультисервисных сетей и сейчас остаются открытыми, что объясняется новизной и сложностью проблемы, а также многообразием возможных вариантов источников мультимедийной нагрузки и количества объединяемых каналов для предоставления услуги пользователю цифровой сети.
Таким образом, актуальность проблемы создания методов анализа и синтеза мультисервисной сети обусловлена необходимостью:
• обеспечения оптимальных режимов работы и определения пропускной способности и приемлемой загрузки участков сети;
• эффективного использования оборудования сети связи при обеспечении требуемого качества обслуживания. Требуется как можно более точное определение вероятности потерь, так как ошибки, допущенные в расчетах, приводят к завышенным требованиям к оборудованию, необходимому для обслуживания трафика.
Целью диссертационной работы является:
1. Разработка и исследование приближенного метода, пригодного для расчетов пропускной способности элементов мультисервисной сети при полнодоступной и неполнодоступной схемах доступа к среде передачи данных по технологии ATM.
2. Разработка алгоритма определения допустимой величины нагрузки для отдельно взятого участка мультисервисной сети при заданных требованиях к вероятности потери вызова.
3. Развитие математических моделей для участка МСС с учетом ограниченного числа источников нагрузки.
Методы исследования. Для решения поставленных задач использовались аналитические методы теории вероятностей, теории массового обслуживания и теории телетрафика. Достоверность полученных результатов диссертационной работы подтверждается строгостью применяемого математического аппарата, сравнением точных и приближенных значений характеристик качества обслуживания вызовов, а также положительными результатами апробации и внедрения. Численные расчеты проведены с использованием современных средств вычислительной техники.
Научная новизна работы заключается в следующем:
1. Осуществлен сравнительный анализ уровня развития основных средств телекоммуникаций в России и Республике Бурятия. Для этой цели ' произведен расчет индексов развития по отдельным направлениям и средневзвешенного индекса, что позволило получить количественные параметры, характеризующие динамику развития средств связи в республике Бурятия и возможность внедрения фрагментов мультисервисной сети.
2. Исследована скученность мультисервисной нагрузки в зависимости от параметров, характеризующих её состав. В результате установлено, что коэффициент скученности наиболее существенно зависит от профиля нагрузки, с помощью которого задается относительный вес каждой группы источников.
3. Исследована точность приближенных моделей, основанных на Z-аппроксимации, при исходных данных, которые характерны для современных потоков мультисервисной нагрузки, передаваемых с помощью технологии ATM.
4. Решена задача определения допустимой величины мультисервисной нагрузки для участка сети при заданных требованиях к вероятности потерь вызова.
5. Разработана новая математическая модель для мультисервисной системы при ограниченном числе источников нагрузки. По сравнению с аналогичными результатами, полученными ранее Ершовым В. А. и его учениками, предлагаемая модель носит более общий характер и доведена до конечных математических соотношений, выражающих вероятности состояний системы.
Практическая ценность результатов диссертационной работы заключается в том, что исследование перечисленных выше вопросов создает основу для более достоверного анализа и прогнозирования качества обслуживания пользователей современных мультисервисных сетей. В работе излагаются принципы построения моделей, которые могут быть использованы для вычисления характеристик обслуживания и оценки величины ресурса, необходимого для передачи с заданным качеством нагрузки мультисервисных сетей. Основные алгоритмы и методы доведены до реализации в виде программ, что способствует более широкому применению полученных результатов в практической деятельности по разработке, проектированию и эксплуатации сетей связи. С использованием простых вычислительных средств они позволяют:
1) определять пропускную способность участков’мультисервисной сети в случае полнодоступного и неполнодоступного пучка каналов- 2) осуществить расчет допустимого числа потоков нагрузки разного типа при условии, что известно количество каналов на данном участке сети и соблюдаются установленные требования к качеству обслуживания вызовов- 3) рассчитывать потери по вызовам на участках сети при ограниченном количестве источников мультимедийной нагрузки или при комбинированном характере источников, когда объединяются простейший и примитивный потоки.
Реализация результатов исследований. Основные результаты, изложенные в диссертации, включены в отчеты по НИР «Разработка моделей и алгоритмов оптимизации стратегии модернизации телекоммуникационных сетей» (шифр «Аспект-Сибирь»), а также применяются в учебном процессе кафедры АЭС СибГУТИ и БФ СибГУТИ в лекционных курсах «Теория телетрафика» и «Сети связи». Программы, с помощью которых реализуются теоретические результаты исследований, использовались для решения задач анализа технического состояния и планирования развития телефонной сети г. Улан-Удэ. Это подтверждается актами внедрения, которые приведены в приложении диссертации.
Апробация работы. Основные результаты, изложенные в I диссертационной работе, докладывались и обсуждались на: Международной научно-технической конференции «Информатика — и проблемы телекоммуникаций» (Новосибирск, 2002 г.), Российских научно-технических конференциях (Новосибирск, 1999, 2000, 2004 г. г.), научно-практических конференциях ВСГТУ (Улан-Удэ, 1999, 2001 г. г.), Региональной научно-технической школе-семинаре «Современные проблемы радиотехники СПР-2003» (Новосибирск, 2003 г.).
Публикации. Основные результаты диссертации опубликованы в 9 печатных работах, их них 1 статья, тезисы 7 докладов в материалах Международных, Российских и Региональных научно-технических конференций, 1 учебное пособие.
Структура и объем работы. По структуре диссертация состоит из введения, четырех глав, заключения. Объём работы составляет 153 страницы,.
4.3 Выводы.
1 В известных математических моделях для решения задач теоретической оценки качества предоставляемых услуг и проектирования необходимых ресурсов мультисервисной сети предполагается, что нагрузка поступает от бесконечного числа источников, то есть является пуассоновской. Это можно считать оправданным только для магистральных участков сети.
2 В работах Ершова В. А. с соавторами предлагаются способы расчета пропускной способности звена мультисервисной сети при ограниченном количестве источников нагрузки, что вполне оправдано. Такая постановка вопроса весьма актуальна, так как число абонентов некоторых классов может быть небольшим. В отличие от этих работ, где рассматривается не более трех классов пользователей, предлагаемая модель с ограничением числа источников носит самый общий характер и содержит системы уравнений статистического равновесия.
3 По сравнению с опубликованными раннее работами, которые посвящены исследованиям полнодоступного пучка каналов при обслуживании мультисервисной нагрузки указанного типа, в диссертации рассматривается также наиболее общий случай с точки зрения состава источников, когда объединяются простейший и примитивный потоки поступающих вызовов, формируемых соответствующими группами источников нагрузки. При этом для разных групп источников нагрузки характерны следующие особенности:
— поток поступающих вызовов является простейшим с заданной суммарной интенсивностью, что теоретически соответствует бесконечному числу источников нагрузки:
— количество источников нагрузки в группе ограничено и для вызовов определенного типа все источники характеризуются некоторой постоянной в группе интенсивностью посылки вызова в свободном состоянии.
Для вероятностей состояний рассматриваемого пучка составлена система уравнений статистического равновесия.
4 В результате решения системы уравнений статистического равновесия для каждого из рассмотренных случаев получены явные формулы, которые позволяют вычислить:
1) индивидуальную вероятность потери вызова каждого типа;
2) средние потери по вызовам и по нагрузке для системы в целом.
При разных профилях поступающей нагрузки представлены примеры численных расчетов, с помощью которых анализируются зависимости характеристик системы от ее параметров.
5 Разработанные алгоритмы и программы позволяют рассчитывать потери по вызовам в пучке заданной емкости от источников разных видов с помощью выведенных формул. Полученные графические иллюстрации показывают зависимости вероятности потерь вызовов на отдельном участке мультисервисной сети от параметров системы и характера нагрузки.
ЗАКЛЮЧЕНИЕ
.
Проблема эффективного использования пропускной способности телекоммуникационных систем и средств связи в условиях интегрированного разнородного трафика, а также задача определения сетевых ресурсов, необходимых для предоставления пользователю мультисервисного обслуживания с гарантированным уровнем качества обслуживания являются актуальными и важными, имеющими большое теоретическое и практическое значения.
В процессе решения указанных задач в диссертационной работе получены следующие результаты:
1 На основе собранных данных произведен подробный анализ уровня развития новых технологий по республике Бурятия и внедрения услуг мультисервисной сети в г. Улан-Удэ, и, как результат, рассчитан средневзвешенный индекс развития, являющийся критерием создания фрагментов мультисервисной сети.
Сделан вывод о создании хороших условий для построения и развития мультисервисной сети связи, которая является одной из основных целей,, поставленных Министерством связи России.
2 Исследован один из важных параметров мультисервисной нагрузкиэто коэффициент скученности, выражающийся через математическое ожидание и дисперсию нагрузки. Скученность обусловлена неординарностью потока занятий каналов, поэтому при поступлении вызовов от источников разных категорий коэффициент скученности мультисервисной нагрузки равен средневзвешенному числу каналов, которые требуются для обслуживания вызовам отдельных категорий, с весами, равными интенсивности нагрузки на каналы, создаваемой вызовами этих категорий. В результате исследования зависимости скученности от различных параметров, выявлена следующая закономерность: — коэффициент скученности не зависит от конкретной численности источников нагрузки, а зависит только от профиля нагрузки, с помощью которого задается относительный вес каждой группы источников (причем скученность является возрастающей функцией от относительного веса группы источников).
3 Произведена оценка точности приближенного метода расчета потерь вызовов при различных значениях параметров. Показано, что точность приближенного метода является удовлетворительной. Отмечено, что относительная погрешность приближенного метода будет минимальна для тех потоков вызовов, для которых кратность вызова близка к скученности объединенной нагрузки.
4 Разработан алгоритм проектирования нагрузки на участке мультисервисной сети с помощью метода эквивалентных замен и метода условного разделения пропускной способности пучка. Его применение предполагает вычисление коэффициента скученности поступающей нагрузки, расчета нормы на средние потери вызовов, определение пропускной способности пучка каналов относительно источников нагрузки /-он категории, расчет коэффициента пропорциональности, через который устанавливается взаимосвязь между профилем мультисервисной нагрузки и численностью соответствующих источников.
5 Разработана математическая модель для мультисервисной системы при ограниченном числе источников нагрузки. Данная модель носит самый общий характер и содержит системы уравнений статистического равновесия.
Результаты расчета вероятностных характеристик и определения пропускной способности звена мультисервисной сети в соответствии с предложенным методом наглядно показывают, что для проведения практических расчетов интенсивностей потоков вызовов, нагрузок и потерь могут быть использованы достаточно простые расчетные формулы при высокой вычислительной эффективности.
Для более точного анализа пропускной способности участка мультисервисной сети разработана математическая модель при комбинированном характере обслуживаемой нагрузки, когда часть нагрузки является пуассоновской (от бесконечного числа источников вызовов), а другая часть создается ограниченным числом источников. Проведенные численные расчеты показали, что значения вероятностей потерь, вычисленных с учетом ограниченного числа источников нагрузки существенно ниже. Разработанные методы оценки показателей совместного обслуживания отличаются хорошей точностью и приемлемы для инженерных расчетов. program BROAD22- uses graphlabel LI, L2- const vmax=120- N=3- Dcl=0.05;
Q:array[l.N] of integer =(1,2,3) — Pr: array[l.N] of real=(0.8,0.15,0.05) — var gr, V, g, c, i, il J, k, iv, M, W, Nv: integerp2, ac, cl, p0, p3, pl, pmy, y, z, D, yl, zn, SA: realx:array[l.N] of integerAl, P: array[l.N]of realPx, Py: array[l.N*20] of integery2: array [l.N, 0. vmax] of real, flhtextgrdriver: integergrmod: integergrppath: stringsize: word;
Pk: Pointerprocedure graf (cl: realp2: real) — label mlvar jl, j, u: integerxl, yl, x2, y2:integerx, y: reals: stringbegin setcolor (15)-u:=0- line (50,50,50,450) — line (50,450,600,450) — x:=cl* 1200- y:=ln (p2)*0.4343*80- xl :=round (x)-yl :=round (y) — Px[gr]: =xl-Py[gr]-=ylStr (cl:0:2,s) — line (xl-190,453,xl-190,447)-outtextxy (xl-200,455,s) — setlinestyle (userbitln,$ 1010, l)-setcolor (6)-line (xl-190,50,xl-190,447);
И1″ setcolor (6)-setlinestyle (userbitln,$ 1010,l)-line (530,50+j*80,55,50+j*80) — setlinestyle (solidln, 0,2)-setcolor (15) — line (45,50+j*80,55,50+j*80);
Str (-j:0,s)-outtextxy (35,45+j*80,s)-Str (10:0,s)-outtextxy (20,55+j*80,s) — for jl:=l to 10 do line (48,50+jl*8+j*80,52,50+jl*8+j*80)-endx2:=Px[gr-6]-y2:=Py[gr -6]- u:=round (gr/6)-lu:=gr-(u*6)-2- if gr<7 then goto mlif u=l then begin setcolor (l)-setlinestyle (solidln, 0,2)-endif u=2 then begin setcolor (2)-setlinestyle (solidln, 0,2)-endif u=3 then begin setcolor (3)-setlinestyle (solidln, 0,2)-endif u=4 then begin setcolor (9)-setlinestyle (Dashedln, 0,2)-endif u=5 then begin setcolor (10)-setlinestyle (Dashedln, 0,2)-endif u=6 then begin setcolor (l l)-setlinestyle (Dashedln, 0,2)-endIine (x2−190,50-y2,xl-190,50-yl) — ml: endfunction ERL (v, y-real):realvar i, iv: integerp, del: realbegin iv:=trunc (v)-del:=v-ivp:=y*(2-del+y)/(del+y*(2+y)) — for i:=l to iv do p:=p*y/(i+del+p*y) — ERL:=p endbegin grdriver:=VGAgrmod:-VGAHigrppath: -d: tpbgi';
InitGraph (grdriver, grmod, gфpath) — assign (fH,'resbr4.txt')-rewrite (fll) — v:=30-while v<=vmax do begin gr:=0- write (fll,'V=', V:3) — cl:==0.2-while cl<=0.65 do begin writeln (fll) — write (fll,' cl=', cl:4:3)-writeln (fll) — ас:=сРУfor i:=l to N do Al[i]: =ac*Pr[i]- zn:=0- for i:=l to N do begin x[i]: =0-p[i]:=0-y2[i, 0]:=1-M:=V Div Q[i]- for k:=l to M do y2[i, k]: =y2[i, k-l]*Al[i]/k endL2: W:=0- for i:=l to N do W:=W+X[i]*Q[i]- If W>V then begin {1} c:=nLI: ifX[c]=0 then begin c:=c-lgoto LIendif c>l then begin {3} X[c]: =0- X[c-l]: =X[c-l]+lGoto L2-end end else{l} begin{4} pl:=lfor i:=l to n do pl*y2[i, X[i]]-zn:zn+plfor i:=l to N do if W>v-Q[i] thenp[i]: =p[i]+plX[n]: =X[n]+l-Goto L2 end{4};
SA:=0-pmy:=0- write (fU,' pt=') — for i=l to N do begin p[i]: =p[i]/znwrite (fll, p[i]: 9:6)-p2:=p[i]-gr:=gr+lgraf (cl, p2) — pmy:=pmy+p[i]*Al [i]*Q[i]- SA:=SA+Al[i] *Q[i] end-writeln (fll) — pmy:=pmy/SAy:=0-D:=0- for i:=l to N do begin y:=y+Q[i]*Al[i]-D:=D+sqr (Q[i])*Al[i] endz:=D/y-yl :=y/zpO:=Erl (V/z, yl) — write (fll,' Plind-) — for i:=l to N do begin p3:=(l+(l-pO)/(l+pO)*(Q[i]/z-l))*pOwrite (fll, p3:9:6)-gr:=gr+l-p2:=p3- graf (cl, p2) — end-writeln (fll) — cl:=cl+Dcl endreadln-InitGraph (grdriver, grmod, gфpath) — v:=v+30- end-gr-=0- close (fll)-end. program BROADdпотери в ЦСИО по точной и приближенной формуле} label L1, L2- const vmax=50- {максимальное число каналов} N=3- (число потоков} Dcl=0.05;
Q:array [1.N] of integer = (1,2,3) — Pr: array [1.N] of real = (0.8,0.15,0.05) — var c, i, il, j, k, v, iv, M, W, EV: integerd, pl, p2, p0,p3,pmy, y, z, DD, vl, v2, yl, zn, SA, g, vz, RV, cl, ac: realx: array [I. N] of integer;
A1,P:array [1.N] of real;
B, gm: array [0.vmax] of realy2: array [1.N, 0. vmax] of realtl: textprocedure geifa (vl, v2: integerd, yl, y2: realvar pb, pf: real) — var Ck, Cl, gk, gl, s, sb, sf, x: realv, k, l, kl: integerbegin v:= vl+v2- Cl:=ls:=l-gl:=lsb:=0- sf:=v2- forl:=0 to v2 do begin if 1 >0 then begin Cl:=((v2-l+l)/y2)*Clx:=v-l+lif x >d then begin gl:=((x-d)/x)*glCl:= Cl/(l-gl) — sb:= sb + gl *C1 ends:=s+Cl-if v2>l then sf:=sf + (v2-l)*Cl endfor k:=l to vl do begin if k=l then begin Ck:= CIgk:= gl endkl:=k+l;
Ck:= ((vl-k+1) *kl/k/y 1) * Ckx:= v-kl+1- if x >d then begin gk:=((x-d)/x)*gkCk:= Ck/(l-gk) — sb:= sb+gk*Ck ends:=s+Ckif v2>l then sf: = sf + (v2-l)*Ck end endpb: =(sb + l)/spf: = sf/s/y2 endbegin assign (fl.'rbrdl.pas') — rewrite (fl) — v:= 30- while v <- Vmax do begin d:= 0.9*vwrite (fl, 'v-, v:3) — cl:=0.3- while cl<0.75 do begin writeln (fl) — write (fl,'cl=', cl :4:3) — writeln (fl) — ac:= cl*vfor i:=l to N do Al[i]: = ac * Pr[i]- gm[v]: =lfor i: = v downto 1 dogm[i-l]: = gm[i]*(i-d)/i;
B[0]: =1- for i:=l to v do B[i]: =B[i-l]*(l-gm[i-l]);
Zn:=0- for i:=l to N do begin x[i]: = 0- p[i]: =0- y2[i, 0]: =lM:= v Div Q[i]- for k:=l to M do y2[i, k]: = y2[i, k-l]*Al[i]/k endL2: W:=0- for i:= 1 to N do W:=W+X[i]*Q[i]- if W>v then begin{l} c:=nLI: if X[c]=0 then begin.
LI: ifX[c]=0 then begin c:= c-1- goto LIendif c>l then begin {3} X[c]: =0- X[c-1]: = X[c-1]+1- Goto L2- end end else {1} begin{4} pl:= B[W]- for i:=l to n do pl:=pl*y2[i, X[i]]-zn:=zn+plfor i:=l to N do begin g:=lif W d then g:=l — B[W+Q[i]]/B[W] else g-=0 endp[i]: =p[i]+pl*g end;
X[n]: = X[n]+1- Goto L2 end {4};
SA:=0- pmy:=0- for i:=l to N do begin p[i]: =p[i]/znpmy: = pmy + p[i]*Al[i]*Q[i]- SA:=SA+Al[i]*Q[i]- writeln (fl,' Рточн-, p[i]: 9:6) — endpmy:=pmy/SA- {приближенный метод} y:=0- DD:=0- for i:=l to N do begin y:=y+Q[i]*Al[i]- DD:=DD+sqr (Q[i])*Al[i] endz:= DD/yyl:=y/zvz:=V/zEV:=trunc (vz) — RV:=Vz-EVgeifa (EV, 0, d/z, yl, yl, pl, pmc) — geifa (EV+ l, 0, d/z, yl, yl, p2, pmc) — pO:=pl*(lRV)+p2*RVfor i:=l to N do begin p3:= (l+(l-pO)/(l+pO)*(Q[i]/z-l))*pOwriteln (fl, 'Plind-, p3:9:6) — endcl:=cl+Dclend;
V:=V+10-cndclose (fl) — end. program Mult2- const na=5- {number of intensities} n=3- {number of call types} vmax=18- {max.number of channels} Mmin=3- {min.number of call sources} type arn=array [l.n] of realain=array[l.n] of integerconst R: ain=(1,2,4) — {kratnost}.
T:arn=(l.0,1.0,1.0) — {mean holding times} var i, Sa, il, ivar, v: integerY, YO, PC, PY, del: realM0:ain- {numbers of call sources} A, Al, La, Pi: arnfl :text;
Procedure MultFin (k:integer-var Y, PC, PY: real-var La, Pi: arn) — label LI, L2- var B, i, j, ind: integer;
GL, WP, W, WY, SW. GP:realml, x: ain;
M, L, SL, SP: arn;
W0:array[l.n, 0. vmax] of real;
Procedure PLUS (il:integervar ind: integer) — var i, j: integerbegin ind:=0- for i:=il downto 1 do if x[i].
SW: =0 — GL: =0 — GP: =0 — WP: =0 — WY: =0 — for i:=l to n do begin x[i]: =0- SL[i]: =0-SP[i]:=0- M[i] :=M0[i] * (k+1E6 * (к-1)) — A[i] :=A1 [i]/(k+lE6 * (к-1)) — ml[i]: =v div R[i] - if ml[i]>M[i] then ml[i]: =trunc (M[i]) — W0 [ i, 0 ]: =1 — for j:=l to ml[i] do.
W0[ij]: -W0[i, j-1]*(M[i]-j +l)/j*A[i]*T[i] end {i};
LI: B:=0- for i:=l to n do B:=B+x[i]*R[i] - if B>v then begin i:=nwhile (x[i]=0) do i:=i-lif i=l then goto L2 else begin.
PLUS (i-l, ind) — if ind=0 then goto L2- goto LI end end else begin {3} W:=lfor i:=l to n do W:=W*W0[i, x[:i] ] - SW:=SW+Wfor i:=l to n do begin.
L[i] :=A[i] *(M[i]-x[i]) *WSL[i]: =SL[i]+L[i] ;
GL:=GL+L[i] - if B>v-R[i] then begin SP[i]: =SP[i]+L[i]- GP:=GP+L[i] end endPLUS (n, ind) — if ind=l then goto LI end {3}- L2: for i:=l to n do begin La[i]: =SL[i]/SWPi[i]: =SP[i]/SL[i];
WY:= WY+SL[i]*T[i]*R[i]- Wp:=Wp+SP[i]*T[i]*R[i] end;
Y:=WY/SW-PC:=GP/GL-PY:=WP/WY end {MiltFin}- begin assign (fl, 'resmul2.pas') — rewrite (fl) — for i:=I to n do MO[i]: =Mminfor il:=0 to n do begin if il>0 then M0[il]: =M0[il]*2- Sa:=0- for i:=l to n do Sa:=Sa+R[i]*MO[i]- writeln (fl,' il—, il:2,' Sa=', Sa:3) — v:=6- while (v<=vmax) do begin write (fl,' v—, v:2) — for ivar:=l to n do begin for ii=l to n do Al[i]: =v/Sa*ivar/naMultFin (l, Y, PC, PY, La, Pi) — write (fl,' M0-) — for i:=l to n do write (fl, M0[i]: 3) — write (fl/A=-, A[l]: 9:6) — write (flY-, Y:8:4,! PC=', PC:8:6,' PY=', PY:8:6) — for i:=l to n do write (fl,' Pi=', Pi[i]: 9:6) — writeln (fl) ;
MultFin (2,Y0,PC, PY, La, Pi) — del:=(YO-Y)/YOwrite (flYO-, Y0:8:4,' deKdel:7:4,' PC=', PC:8:6,' PY=1, PY: 8: 6) — for i:=l to n do write (fl,' Pi=', Pi[i]: 9:6) — writein (fl);
Y^J*itdn (f 1 j endv:=v+6- writeln (fl) end endclose (fl) end. proqram MultTestconst na=5- {number of intensities} n=3- { number of call types} vmax=12- {max. number of channels} type arn=array[l.n]of realarn^array [l.n] of integerconst M0: ain= (10,5,2) — {numbers of call sources} R: ain= (1,2,4) — {kratnost}.
A0:arn= (0.2,0.1,0.02) — {call intensities for one free source} T: arn= (1.0,1.0,1.0) — {mean holding times} var i, ia, il, ivar, v: integerY, Y0, PC, PY, del: realA, Al, La, Pi: arnfl: text-Procedure MultFin (k, ia: integervar Y, PC, PY: realvar La. Pi: arn) — label LI. L2- varB, i, j, ind: integer;
GL, WP, W, WY, SW, GP: realml, x: ainM, L, SL, SP: arnW0: airay [l.n, 0. vmax] of realProcedure PLUS (il: integervar ind: integer) — var ij: integerbeqin ind:=0- for i:=il downto 1 do if x[i]M[i] then ml[i]: =trunc (M[i]) — W0[i.0]: =lforj:=l to ml[i] do.
W0[i, j]-=WO[ij — l]*(M[i] -j +1) /j *A[i] *T[i] end{i}- L1: B:=0- for i:=l to n do B:=B+x[i] *R[i]- if B>v then beqin i:=nwhile (x[i]=0) do i:=i- 1- if i=l then doto L2 else beqin PLUS (i-l, ind);
If ind=0 then doto L2- doto LI end end else beqin {3} W:=lfor i:=l to n do W:=W*W0 [i, x[i]]- SW:=SW+Wfor i:=l to n do beqin.
L[i]: =A[i] * (M[i] -x[i]) *WSL[i]: =SL[i]+L[i]- GL:=GL+L[i]- If B>v — R[i] then beqin SP[i]: =SP[i]+L[i]- GP:=GP+L[i] end endPLUS (n, ind);
Ifind=l then goto LI end{3}- L2: for i i=l to n do beqin La[i]: =SL[i] / SWPi[i]: =SP[i] / SL[i]- WY:=WY+SL[i] *T[i] *R[i]- WP:=WP+SP[i] *T[i] *R[i] endY:=WY / SWPC:=GP / GLPY:=WP / WY end {MultFin}- beqin assign (fl, 1 resmult. pas ') — rewrite (fl) — write (fl,' число источников M0= ') — for i:=l to n do write (fl, M0[i]: 3) — writeln (fl) — write (fl, 1 кратность R= ') — for i:=l to n do write (fl, R[i]: 3) — writeln (fl) — for i:=l to n do Al[i]: =A0[i]- for il:==0 to n do beqin {il}- if il >0 then A1 [il]: =Al[il]*2- v:=4- while (v<=vmax) do beqin {v} writein (fl,' v-, v:2) — for ivar:==l to na do beqin {ivar} ia:=ivar*v div 4;
MultFin (1, ia, Y, PC, PY, La, Pi) — write (fl,' A=*) — for i:=l to n do write (fl, A[i]: 9:6) — writein (fl,' Y= Y:10:6,' PC= PC:8:6,' PY= ', PY:8:6) — for i:=1 to n do write (fl,' La=', La[i]: 9:6,' Pi= Pi[i]: 9:6) — writein (fl);
MultFin (2, ia, Y0, PC, PY, La, Pi) — del:=(Y0-Y)/Y0- writein (fl,' Y=Y0:10:6,' del= del:7:4,' PC= PC:8:6,' PY=PY:8:6) — for 1=1 to n do write (fl, 1 La= La[i]: 9:6,' Pi= Pi[i]: 9:6) — writeln (fl);
•. 1 / pi (kLa ь|> kb «1* ф Ф ф чЬ ф ф ф ф ф ф ф ф ф ф ф Ф Ф Ф Ф Ф Ф Ф ф Ф ф ф ф Ф Ф Ф ф ф I.
X/t*1TPl Ты TI i end {ivar}- v:=v+4- writeln (fl) end{v} end{il}- close (fl) end.