MOGUĆNOSTI PRIMENE FRAKCIONOG RAČUNA U MODELOVANJU TELEKOMUNIKACIONOG SAOBRAĆAJA

Size: px
Start display at page:

Download "MOGUĆNOSTI PRIMENE FRAKCIONOG RAČUNA U MODELOVANJU TELEKOMUNIKACIONOG SAOBRAĆAJA"

Transcription

1 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 ПРЕГЛЕДНИ ЧЛАНЦИ REVIEW PAPERS ОБОЗОРНЫЕ СТАТ ЬИ MOGUĆNOSTI PRIMENE FRAKCIONOG RAČUNA U MODELOVANJU TELEKOMUNIKACIONOG SAOBRAĆAJA Braka D. Mikavica, Aleksadra M. Kostić-Ljubisavljević, Vesa M. Radojić-Đogatović Uiverzitet u Beogradu, Saobraćaji fakultet b.mikavica@sf.bg.ac.rs, a.kostic@sf.bg.ac.rs, v.radojic@sf.bg.ac.rs OBLAST: telekomuikacije VRSTA ČLANKA: pregledi člaak JEZIK ČLANKA: srpski DOI: /vojtehg Sažetak: Frakcioi raču je oblast matematičke aalize koja se bavi izučavajem i primeom izvoda i itegrala proizvoljog reda. Ovom teorijom bavili su se mogi pozati matematičari među kojima su Ojler, Rima, Liuvil, Abel i Furije. Predložeo je više defiicija za izračuavaje izvoda i itegrala ecelog reda. U ovom radu daje se pregled ekih predložeih defiicija, kao i osove postavke frakcioog račua sa posebim aglaskom a mogućostima jegove primee u domeu modelovaja telekomuikacioog saobraćaja. Čijeica je da frakcioi raču posledjih deceija alazi sve veću primeu u razim aučim oblastima. Modeli zasovai a frakcioom račuu pokazali su se korisim u fizici, mehaici, elektrotehici, biohemiji, medicii, ekoomiji, teoriji verovatoće. U ovom radu aalizira se mogućost primee frakcioog račua u modelovaju telekomuikacioog saobraćaja. Istraživaja pokazuju da se karakteristike telekomuikacioog saobraćaja a lokalom i globalom ivou u mreži, kao što su samosličost i zavisost u dugom opsegu, efikasije mogu opisati pomoću frakcioog račua umesto kovecioalih stohastičkih procesa. U radu su prikazai predložei modeli zasovai a frakcioom račuu koji modeluju feomee prisute u savremeim telekomuikacioim mrežama. Ključe reči: modeli, saobraćaj, telekomuikacije, frakcioi raču. ZAHVALNICA: Ovaj rad je deo rezultata istraživaja a projektu TR32025 koji fiasira Miistarstvo prosvete, auke I tehološkog razvoja Republike Srbije 64

2 Uvod Telekomuikacioe mreže beleže eorma razvoj. Potrebe korisika u pogledu ovih servisa i aplikacija koje zahtevaju veće propuse opsege i veći kvalitet servisa svakim daom rastu. Rast tražje utiče a eophodost izgradje adekvate arhitekture mreže, plairaje, održavaje, kotiuirai razvoj i uapređivaje. U tom smislu, od velikog začaja je evaluacija efikasosti telekomuikacioih mreža korišćejem metoda mereja, aalize i simulacije ovih mreža. Merejem i statističkom aalizom telekomuikacioog saobraćaja otkriveo je da saobraćaj u mreži pokazuje velike eregularosti (burstiess), kako u pogledu varijabilosti iteziteta saobraćaja, tako i u pogledu oblika fukcije autokorelacije. Uočeo je da saobraćaj ima fraktale karakteristike samosličost i zavisost u dugom opsegu (Devetskiotis, da Foseca, 2005). Usled toga, potreba je velika širia propusog opsega, a vrlo često je to jeda od uzroka eefikasosti mreže. Za razliku od ustaljeih modela, zasovaih a Poasoovoj raspodeli, koji se sreću u mrežama sa komutacijom kola, kod modela zasovaih a samosličosti telekomuikacioog saobraćaja javljaju se problemi koje je teško progozirati, izmeriti i kotrolisati saobraćaj u mreži. Stoga, mereje, aaliza i modelovaje samosličog telekomuikacioog saobraćaja još uvek predstavlja izazov. Različita istraživaja su pokazala da se saobraćaj u savremeim telekomuikacioim mrežama može adekvato opisati statističkim modelima zasovaim a frakcioom račuu. Frakcioi raču predstavlja oblast matematike koja se bavi izučavajem i primeom izvoda i itegrala ecelog reda. Posledjih deceija alazi sve veću primeu u razim sferama auke, od medicie, fizike, biohemije, ekoomije, teorije verovatoće, mehaike, elektrotehike. Težište ovog rada je a mogućosti primee frakcioog račua u savremeim telekomuikacioim mrežama u modelovaju telekomuikacioog saobraćaja. Rad je kocipira a sledeći ači. Nako uvodog razmatraja, daju se osove postavke frakcioog račua. Broji pozati matematičari koji su se bavili izučavajem frakcioog račua dali su svoje defiicije i u ovom radu se daje pregled ekih od predložeih defiicija. U aredom poglavlju se opisuju feomei koji se sreću u savremeim telekomuikacioim mrežama samosličost (eg. self-similarity) i zavisost u dugom opsegu (eg. Log ru depedece - LRD). Zatim se daje pregled ekih predložeih modela saobraćaja zasovaih a frakcioom račuu. Na kraju se daju zaključa razmatraja. Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

3 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 Osove postavke frakcioog račua Izučavaje frakcioog račua počije u 7. veku sa Lopitalom i Lajbicom. Lakroa je prvi objavio eceli red izvoda. Tako, za y =, a R +, a pokazao je da /2 d y Γ ( a + ) a /2 =. /2 d Γ ( a + /2) Zapravo, dobio je da /2 d / d = 2 / π (isti rezultat je dobio i Rima-Liuvil, koji važi i daas). Iako je aziv frakcioi raču zapravo pogreša, a ozaka itegracija i diferecijacija proizvoljog reda adekvatija, uobičaje je aziv frakcioi raču, koji se koristi još od doba Lopitala. Furije, koji je 822. izveo vredost itegrala fukcije f = d cos π dp, formalo je izrazio verziju sa izvodima 2 R R ν d ν νπ f = f d p cos p( ) dp ν d 2π + 2 R R, gde se vredost ν odosi a bilo koju vredost, pozitivu ili egativu. Abel je primeio frakcioi raču za rešavaje itegrale jedačie koja se javlja u formulaciji tautohroog problema (odrediti oblik krive tako da je vreme spuštaja tela zaemarljive mase iz datu krivu, bez treja, pod uticajem gravitacije, ezaviso od početog položaja tela). Abelova itegrala jedačia ima oblik: /2 k = ( t) f ( t) dt 0. () Abel je izučavao opštije forme itegralih jedačia sa jezgrom oblika ( t). Itegral (), sa izuzetkom faktora / Γ ( /2), poseba je slučaj koačog itegrala koji defiiše frakcioo itegraljeje reda ½. U itegralim jedačiama oblika () fukcija f je epozata. Abel je desu strau jedačie () apisao kao π /2 d / d /2 f. Zatim je primeio operaciju d / d a obe strae jedačie i dobio /2 /2 je: /2 d k π f /2 d =, (2) 66

4 jer ovi frakcioi operatori imaju osobiu da D 2D 2 f = D f = f. Stoga, kada se izračua frakcioi izvod reda ½ kostate k u (2), f ( ) je određea. To je izuzeto dostiguće u oblasti frakcioog račua. Važo je uočiti da frakcioi izvod kostate ije uvek jedak uli (Miller, Ross, 993). Prvi ozbilja pokušaj da se dâ logiča defiicija frakcioih izvoda a a pripada Liuvilu. Liuvil je pošao od pozatog rezultata D e = a e gde je d D =, N, i proširio ga, ajpre a poseba slučaj kada je, 2 d ν = 2 a = ν a ν a a zatim ga proširio a proizvolja red, ν R+ sa D e = a e. Prikazao je a f c k ke k = 0 ν ν a D f = c k kake. k = 0 fukciju f ( ) kao = i defiisao izvod proizvoljog reda ν sa Njegov drugi metod odosio se a eksplicitu fukciju I a u e u du 0 a a t a 0 ν a ν Γ a + D = Γ a a. Posmatrao je itegral I = t e dt = Γ a 0 =, a zatim zameom u = t dobio rezultat a (za Rea > 0). Možeći obe strae = I / Γ ( a) sa D ν dobio je ( ν ) a ν S obzirom a to da proizvolja diferecijala jedačia oblika d y = 0 -tog reda ima opšte rešeje y = c 0 + c c, Liuvil je d d y zaključio da frakcioa jedačia oblika = 0, R + takođe treba da d ima odgovarajuće opšte rešeje. Neka je f lokalo itegrabila u itervalu ( a, ), tada je -ta iteracija itegraljeja data kao: u u a : 2... (! ) a a a a. Ι f = du du f u du = u f u du Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

5 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 za skoro svako uz a < < i N. Ako se uzme da je (! ) =Γ, dobija se eposreda geeralizacija itegrala fukcije f ecelog reda > 0, Ι f = ( u) f ( u) du Γ, (desi itegral) (3) a ( ) i sličo za < < b, ( ) b a Ι f = ( u ) f ( u) du Γ (levi itegral), b oboje defiisao za odgovarajuće f. Ideksi uz I ozačavaju graice itegracije u datom redosledu. Može se uočiti da je za = jedačia (3) jedistveo rešeje početog problema ( ) ( ) y = f, y a = y' a =... = y a = 0 Kada je a = jedačia (3) je ekvivaleta Liuvilovoj defiiciji, a za a = 0 dobija se Rimaova defiicija (bez komplemetare fukcije). Najčešće se kaže da je a Ι f Rima-Liuvilov frakcioi itegral reda fuk- cije f. Sa druge stae, jedačie W f = aι f = ( u ) f ( u) du Γ ( ) W f = Ι f = ( u) f ( u) du Γ ( ) azivaju se Vejlovi frakcioi itegrali reda, defiisai za odgovarajuću fukciju f. Levi i desi itegral frakcioih itegrala ( Ι ) i Ι ( ) a f b f povezai su preko Parsevalove jedakosti (frakcioa itegracija po delovima) koja se dobija za a = 0 i b = : f ( 0Ι g) d = ( W f ) g d 0 0. Sledeća svojstva važe za desi itegral frakcioih itegrala (uvođejem adekvatih zamea važi slučaj levog itegrala). Neka je f L a, Ι je koača svuda u itervalu loc ( ). Tada, ako je a >, a f ( a, ) i pripada L (, ) loc a. Ako je a = pretpostavlja se da se f poaša 68

6 tako da za takav itegral kovergira. Pod ovim pretpostavkama frakcioi itegral zadovoljava svojstvo aditivosti: ( β, 0) β + β aι aι = aι f >. Uzimajući u obzir redosled itegraljeja, dobija se: u β β aι aι f = ( u) du ( u t) f ( t) dt Γ Γ( β) a a. β = f () t dt ( u) ( u t) du Γ Γ( β) a t u t Drugi itegral je, uvodeći zameu y =, jedak t + β + β β ( t) ( y) y dy =Β( β, )( t) 0 Γ = Γ + gde je (, β ) Γ( β) ( t ) β + β Β Beta fukcija. Kada se to uvrsti u gorje rezultate, dobija se: (, 0) + aι f = aι aι f N >, što implicira -tu diferecijaciju: d a d + Ι f = Ι f, ( N, 0) a > za svako. (4) Prethodi rezultati važe i za kompleksi parametar, ako se uslov > 0 zamei sa Re > 0. Tada se operacija aι može smatrati holomorfom fukcijom od za Re > 0, što se može proširiti a celu kompleksu rava. Kako bi se ovo pojasilo, polazi se od pretpostavke da je f beskoačo diferecijabila u skupu realih brojeva. Ako je a > tada je ( f ) ( a ) = 0 za = 0,,2... Tada za svako fikso > a itegral u (9) je holomorfa fukcija od za Re > 0. Neka je R i r > 0. Skup Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp (, ) { : (, ) } : 2 Κ r = y R d y < r = y R i yi < r i= 69

7 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 aziva se otvorea kugla oko radijusa r. Skup A R je otvore ako važi A, r > 0, Κ(, r) A. Fukcija f je diferecijabila ili holomorfa a otvoreom skupu Ω ako poseduje izvod u svakoj tački skupa Ω. Važi sledeće: ako je fukcija diferecijabila, oda postoji izvod svakog reda a celom skupu Ω. Takva diferecijabila fukcija aziva se holomorfa. Parcijala itegracija -puta daje: + ( Ι f = Ι f ), ( Re 0, N) a a > (5) Primejujući aditivo svojstvo (4) a izraz sa dese strae (5) i diferecirajući rezultat -puta po, dobija se: d d ( Ι f = Ι f ), ( Re 0, N) a a > (6) čime se pokazuje da su operacije itegracije frakcioog reda i difereciraja reda ad fukcijom f komutative. Ako se poovo pogleda formula (5) može se uočiti da je desa straa C;Re > i izraza holomorfa fukcija od sa većim domeom { } d + jedak je Ι f po formuli (6). Ovako se može proširiti Ι ( ) a d C;Re 0, aalitički defiišući za C i Re 0 a dome { } + d + a a a ( ) a f Ι f = Ι f = Ι f (7) d Za svaki ceo broj > Re, dobija se 0 ( Ι f = f, Ι f = f ), ( N) a a 70 (8) Neka je kompleksi broj takav da je Re > 0 i = [ Re ], gde [ Re ] ozačava ceo deo od Re. Desi itegral frakcioog izvoda reda defiisa je kao d a a = Ι, [ ] D f f d za svako f L ( a, ) loc ( Re ) = + (9) za koje izraz a desoj strai važi.

8 Defiicije itegrala i izvoda proizvoljog reda, Re 0, N : (ekvivaleto) mogu se objediiti a sledeći ači, za ad f ( u) f ( u) du, ( Re < 0) Γ ( ) =, d aι f,( Re > 0; Re < ) d što se često aziva diferitegral fukcije f reda. Ovaj proces često se aziva i frakcioo itegro-difereciraje. Može se uočiti da je levi izraz frakcioog izvoda reda defiisa d b d b ( Re ) kao: D f = ( ) Ι f, [ ] = +. Frakcioi izvod potpuo imagiarog reda = iθ, θ 0 defiisa je kao f ( u) iθ d ad f = = du Γ( iθ ) d, iθ a u a odgovarajući itegral reda = iθ defiiše se kao: iθ d + iθ d iθ aι f = aι f = ( u) f ( u) du d Γ ( + iθ ) d a. Defiicija frakcioog itegro-difereciraja za svako a C je kompleta 0 0 uvođejem idetičkog operatora ad f : = aι f = f za svako = 0. Takođe, može se uočiti da su frakcioi operatori lieari: + = + ad cf c2f2 cad f c2 ad f2, gde su c i c 2 kostate. Pri defiisaju dovoljih uslova za postojaje frakcioih izvoda u relaciji (9) razmatra se slučaj 0< Re <, ako je > (Hifler, 2000). Pretpostavlja, f AC a, b, što se da je f eprekida u koačom itervalu [ ab ], formalo [ ] zači da je f diferecijabila svuda u itervalu ( ab, ) i f ' L ( a, b) a. a [ ab, ] jedaka je f = f '( u) du+ f ( a) = Ι f ' + f ( a) i u itervalu Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

9 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 Zameom ovoga u formuli Ι Ι i Ι komutativi, dobija se: ' a f f ( a) ( 2 ) aι f = aι aι f + Γ Diferecirajem po, dobija se:, uz čijeicu da su operatori. f ( a) d ad f = aι f = aι f ' + a d Γ (0) što pokazuje da, u opštem slučaju, operatori aι i d isu komutativi. d Formula (0) može se proširiti uvođejem kompleksog broja, takvog da je Re. Rezultati su obuhvaćei pretpostavkom koja sledi. Najpre se uvode sledeće ozake: za N, AC [ a, b] ozačava skup od ( ) diferecijabilih fukcija f a [ ab, ], takvih da su f, f ',..., f 0 eprekide a [ ab, ]. Može se uočiti da je AC [ a, b ] jedako AC [ a, b ]. Pretpostavka: ako je f AC [ a, b] f je defiisaa samo u koačom itervalu [ ab, ], tada postoji ad f, Db f za 0< Re <. Štaviše, r f ( a) f '( u) ad f L ( a, b) za r < i ad f = + du Re ( ) Γ ( a) a ( u). Aalogo važi za D b f. Ako je f AC [ a, b], = [ Re ] + tada postoji D a f za Re 0 i važi ( k ) ( ) ( ) f a f u ad f = + du Γ + k Γ + k = 0 a u () Drugi ači da se defiiše frakcioi izvod reda, takođe prema Liuvilu, jeste: ( : = Ι ), [ ] ad f a f ( Re ) = +. (2) Očigledo, f mora biti diferecijabila kako bi izraz sa dese strae (2) postojao. Relacija između dva frakcioa izvoda (9) i (2) data je formulom (), aime ( ) f ( a) k ad f = ad f + a. Γ + k k = 0 72

10 Defiicija (9) ajčešće je korišćea u matematičkim krugovima, dok je defiicija (2) češće korišćea u slučajevima kada su početi uslovi defiisai preko celih izvoda (Hifler, 2000). Neki predložei modeli primee frakcioog račua za modelovaje telekomuikacioog saobraćaja Izvodi i itegrali celog reda imaju jasu fizičku i geometrijsku iterpretaciju. Međutim, prihvatljiva geometrijska i fizička iterpretacija itegraljeja i difereciraja frakcioog reda dugo ije bila pozata. To je oblast koja se daas sve više razvija, kako u teoriji, tako i u praksi. Posledjih godia frakcioi raču alazi primeu u studijama o viskoelastičim materijalima, tokovima fluida, elektromagetoj teoriji i verovatoći, saobraćaju i medicii (Dalir, Bashour, 200). U ovom delu rada težište je a mogućosti primee u telekomuikacioim mrežama Modeli Markova i jegova teorija redova često su izučavai i itezivo primejivai kao alat za evaluaciju performasi. Međutim, aalize zasovae a modelima Markova imaju određee edostatke. To su: u mrežama sa velikim protocima verovatoća odbacivaja paketa u baferima sa ograičeim kapacitetom teško se može odrediti i često se aproksimira verovatoćom odbacivaja u sistemu sa beskoačim baferima (Kim, Shroff, 200), kada se multipleksira veliki broj izvora Markova, veliki kapacitet agregatih dolazih procesa rezultira emogućošću izračuavaja, ed-to-ed aaliza performasi zahteva precizo modelovaje odlazih procesa iz redova za aalizu sledećeg hopa (modelovaje odlazih procesa za Markove izvore takođe je veliki problem), sveobuhvata mereja i aalize saobraćaja u mreži ukazuju a to da paketski saobraćaj ima svojstva samosličosti ili zavisosti u dugom opsegu (LRD), koja se e mogu opisati modelima Markova koje karakteriše zavisost u kratkom opsegu (Cheg et al, 2007). Samosličost podrazumeva da je određeo svojstvo objekta, poddome ekog diamičkog sistema ili vremeska serija, očuvao u odosu a vremesko i/ili prostoro skaliraje. Ako je objekat samosliča ili fraktala, jegovi delovi, kada se uvećaju, u određeom smislu imaju oblik celie. Determiistička samosličost podrazumeva samo jaku formu rekurzive regularosti. Međutim, za potrebe modelovaja saobraćaja, gde je stohastička varijabilost osova karakteristika, determiistička samosličost mora se dalje uopštiti. Za razliku od determiističkih frak- Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

11 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 tala, stohastički fraktali e poseduju egzaktu sličost jihovih delova u poređeju sa celiom. Mera sličosti je oblik grafa sa adekvato izabraim stepeom sličosti. Ako se usvoji da su vremeske serije u saobraćaju uzorci stohastičkog procesa, pri čemu se zaemaruje mera sličosti fokusirajem a određee statistike ili, recimo, reskalirajem vremeskih serija, tada je moguće postići egzaktu sličost objekata. Statistike drugog reda su statistička svojstva koja prikazuju burstiess ili varijabilost. Oblik autokorelacioe fukcije ima začaju ulogu. Autokorelacioa fukcija je kriterijum u odosu a koji se defiiše ivarijatost. Korelacija, kao fukcija kašjeja, opada zato sporije u odosu a ekspoecijalu. Postojaje ovakve zavisosti aziva se zavisost u dugom opsegu (Park, Williger, 200). Zavisost u dugom opsegu je svojstvo vremeskih serija da pokazuju jaku zavisost tokom velikih vremeskih itervala. Formalo se može izraziti a više ekvivaletih ačia. Neka je Xt, t Z diskreta promeljiva stacioare serije. Serija Xt, t Z se aziva log rage depedece (LRD), ako za jeu fukciju kovarijase važi γ () 2 2H t = E X0 EX0 Xt EXt cγ t, t (3) i važi H, 2 gde je c γ kostata veća od 0. Uslov (3) može se ekvivaleto apisati pomoću Furijeove trasformacije kao: 2H f λ c f λ, λ 0, fukcija spektrale gustie X t iλ p gde je f ( λ) ( 2 π) = e γ ( p), λ [ π, π] i f p c je kostata veća od 0. Parametar H aziva se Hurstov parametar. Što je H veće, vremeska zavisost je veća, jer fukcija kovarijase sporije opada ka beskoačosti (Park et al, 20). Postoje različite teorije o uzroku pojave samosličosti. Po prvoj, uzrok samosličosti leži u aplikativom sloju. Druga teorija tvrdi da je poreklo samosličosti a trasportom sloju. U ovoj teoriji posmatraju se efekti primee TCP algoritama kotrole zagušeja a saobraćaj u mreži, kao i situacija koja astaje kada više izvora saobraćaja rutira saobraćaj preko istih putaja. Treća teorija razmatra uzroke a mrežom sloju. Ova teorija, pre svega, povezuje samosličost sa kritičim feomeima koji astaju u tačkama u kojima sloboda tok saobraćaja prerasta u zagušeje (Smith, 20). 74

12 Vremesko-prostoro modelovaje preosa paketa u TCP/IP mrežama U astavku će biti prikaza predložei model koji predstavlja telekomuikacioi saobraćaj kao stohastički proces sa mogućošću beskoačog sredjeg kašjeja. Ovim modelom može se objasiti pojava zavisosti u dugom opsegu i fraktalih osobia tokova podataka. Takođe, prikazaa je formala veza između heavy-tailed raspodela kašjeja, hiperboličkog opadaja kašjeja paketa, fukcije autokovarijase i frakcioih izvoda. Nedava mereja telekomuikacioog saobraćaja pokazala su da se jegove karakteristike efikasije mogu opisati pomoću ecelih izvoda umesto kovecioalih stohastičkih procesa (Zaborovsky et al, Nova iterpretacija frakcioog račua otvara ova područja upotrebe ovih razvijeih matematičkih alata radi razumevaja lokalih i globalih karakteristika saobraćaja u mrežama sa komutacijom paketa. Fraktala dimezija i zavisost u dugom opsegu mogu se uočiti u privatim, lokalim i WAN mrežama. Defiisaje kvaliteta servisa (QoS) za Iteret servise zahteva adekvate modele virtuelih koekcija, u zavisosti od tokova podataka. U opštem slučaju, takva koekcija se sastoji od ekoliko trazitih čvorova i likova. Bez smajeja opštosti pretpostavlja se da se a koekciji koja se razmatra primejuje TCP protokol i da postoji ed-to-ed kotrola. Pored toga, pri defiisaju modela tokova podataka razmatraju se procesi u trazitim čvorovima i komuikacioim likovima, kao što je prikazao a slici. TCP baferovaje kašjeje pri preosu odbacivaje paketa TCP čvor 0 čvor čvor + čvor M Slika Kašjeje/odbacivaje paketa u TCP koekciji Figure Packet delay/dropp processes i a TCP coectio Схема Задержка/потеря пакетов в ТСР соединении Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp Ovi procesi javljaju se u različitim oblicima kašjeja paketa i akcijama jihovog odbacivaja. Posmatra se vreme propagacije paketa. Detaljom elaboracijom sve do ivoa trazitih čvorova (slika 2b) vreme propagacije T može se predstaviti kao suma kašjeja paketa u svakom 75

13 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 M hopu, T = ( t L B i + ti ), gde je t L vreme propagacije paketa između čvo- t kašjeje procesiraja/ baferovaja paketa u čvoru, M je ukupa broj čvorova u koekciji. i rova i + i (kašjeja duž lika), i B a) b) c) predajik prijemik predajik prijemik predajik + prijemik t B t L RTT odbacivaje paketa beskoačo kašjeje t Slika 2 Preos paketa: a) ed-to-ed model; b) čvor-čvor model; c) skok model Figure 2 Packet trasmissio: a) ed-to-ed model; b) ode-ode model; c) jump model Схема 2 Передача пакета: а) модель ed-to-ed б) модель узел-узел с) модель прыжка U savremeim TCP/IP mrežama kašjeje duž lika t L u datoj virtueloj koekciji može se smatrati kostatim (Hifler, 2000). Stoga, varijacija vremea propagacije, što ima izuzeta uticaj a parametre kvaliteta servisa, uzrokovaa je varijacijama kašjeja baferovaja t B. Paket se može odbaciti usled popujeosti bafera ili greške pri rutiraju u trazitim čvorovima. Sa aspekta prijemika traziti čvor se poaša kao zamka za pakete i zadržava ih zauvek. U ovom slučaju, t B uzima vredost beskoačo. Može se uočiti da će uprkos gubitka paketa ili beskoačog kašjeja protokol a trasportom sloju osigurati preos određeih paketa koristeći retrasmisiju. Rezultat toga je porast ukupog broja prosleđeih paketa i, stoga, utiče a efektivu produktivost virtuele koekcije. Virtueli kaal određe je koačim brojem trazitih čvorova i i likova. Svaki traziti čvor karakteriše se prosečim vremeom preosa paketa do sledećeg čvora. Promeljiva t uključuje kašjeje paketa u baferu i vreme propagacije paketa do sledećeg čvora. U ovom slučaju se preos 76

14 paketa između trazitih čvorova može predstaviti kao sekveca diskretih prostorih skokova duž likova (slika 2c). Vremeski iterval između skokova je slučaja promeljiva koja je jedaka vremeu preosa paketa do sledećeg čvora. Takvi skokovi formiraju statistički diamički proces komutacije paketa i direkto utiču a karakteristike saobraćaja u mreži. Za defiisaje matematičkog modela takvog procesa uvodi se fukcija gustie raspodele verovatoće f(t) preosa paketa od čvora do čvora +. Pretpostavlja se da svi traziti čvorovi u virtuelom kaalu imaju ista prava u jihovom dopriosu vredosti ed-to-ed vremea preosa. Mogućost odbacivaja paketa odgovara staju u kojem paket e apušta traziti čvor i ostaje u ovom čvoru zauvek. Tada sredje kašjeje paketa u čvoru zadovoljava uslov < t >= tf () t dt =, 0 f () t > 0, f () t dt =. 0 Odgovarajući izraz za f(t) može se apisati kao: γ f () t =,0 t γ < γ <. + + Zaključuje se da upravo ovakva stuktura fukcije gustie verovatoće karakteriše dugoroču statističku zavisost, koja je detaljo aaliziraa u račuarskim mrežama (Babic et al, 998). Za dalju kostrukciju modela uvodi se fukcija F ( τ ), F τ = f t dt =, γ 0 ( + τ ) ( τ ) gde je τ vreme postojaja paketa u trazitom čvoru virtuelog kaala. Odredimo sada ajverovatiji broj paketa u čvoru u treutku t, koristeći formulu: t ( ; t) = ( ; t τ) f ( τ) dτ + 0 F( t), gde je 0 0 broj paketa u čvoru pre pristizaja paketa iz čvora -. Vredost () ije ograičea veličiom bafera i može uključiti sve pakete u virtuelom kaalu koji se prosleđuju i odbacuju u čvoru. Koristeći ove ozake, jedačia prosleđivaja paketa može se prikazati kao: ( ; ) γ t 0 Γ( γ ) Dt ( ; t) = +, γ t Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

15 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 gde levi deo jedačie predstavlja frakcioi izvod fukcije (;t) sa ekspoetom γ, t γ ( ; τ ) Dt ( ; t) = dτ Γ( γ ) γ. t τ 0 Uzimajući u obzir da promeljiva uzima diskrete vredosti, prethoda jedačia rešava se uz počete uslove: 0( 0) = 0 i 0 ( k) = 0, k =,2,... Tako se dobija 2 Γ ( γ ) ( k; t) = 0 k kγ( γ ) γ 2γ γ+, t Γ( 2γ) t Γ( γ) t ili 2 ( γ) ( γ) Γ Γ ( k; t) = 0 k + γ ( 2 ) 2γ γ t γ t γ +. Γ Γ t Uzimajući u obzir asimptotska svojstva ovog rešeja, moguće je apisati ( 0; t) = 0 t γ. Za k = važi sledeći izraz 2 ( γ) ( γ) Γ Γ ( ; t) = 0 ( 2 ) 2 t γ γ t γ + γ t γ +. Γ Γ Stoga se izračuavaje može astaviti za svako k. Ovaj pristup omogućava izračuavaje korelacioe fukcije. Za počete uslove ( 0; t) = 0δ ( t) moguće je apisati: Γ( γ ) c( m; t) =Γ( γ ) m ( 2γ ) 2γ ( 3 ) 3 t γ γ Γ + Γ + t 2 2γ Γ( γ ) = 0 Γ( γ ) t m ( 2 2γ) ( 2 3γ) γ Γ Γ t Ovaj izraz pokazuje da korelacioa fukcija hiperboličo opada sa promeom t. Zato, pod uslovom da je γ <, takvi procesi pokazuju svojstva koja se opisuju frakcioim račuom. Sa m = 0 izraz za varijasu glasi: 0Γ( γ ) 2γ D() t = c( 0; t) = t, Γ ( 2 2γ ) što je karakterističo za procese sa zavisošću u dugom opsegu i asimptotskom osobiom samosličosti drugog reda. 78

16 Stoga, karakeristike telekomuikacioog saobraćaja, koje su dosta izučavae u moderoj literaturi, mogu se direkto dovesti u vezu sa metodama komutacije paketa, koje se koriste u račuarskim mrežama. Proces preosa paketa može se opisati jedačiama čija se rešeja dobijaju alatima frakcioog račua. Koristeći vremeska asimptotska rešeja, diferecijale jedačie sa frakcioim izvodima mogu se formulisati tako da se mogu dobiti rešeja drugog reda za statističke momete procesa. Specifičost ove karakteristike saobraćaja jeste da jea korelacioa fukcija hiperboličo opada, što se mora uzeti u obzir pri defiisaju metoda za kotrolu saobraćaja (Zaborovsky et al, Modelovaje telekomuikacioog saobraćaja frakcioom aalizom mogućosti i edostaci Empirijska istraživaja telekomuikacioog saobraćaja pokazala su ivarijata fraktala svojstva saobraćaja koja imaju začaju ulogu u kotekstu aalize kvaliteta servisa. U IP telefoiji i video streamig-u preos paketa mogo zavisi od algoritama zasovaih a modelima saobraćaja. Osovi problem u dizaju takvih modela jeste da se stvore uslovi pod kojima je moguća kotrola garatovaih parametara kvaliteta servisa. U ovom delu rada razmatraju se predosti i ograičeja metoda frakcioe aalize za opisivaje fraktalih svojstava izmereog saobraćaja. Fraktala svojstva ukazuju a postojaje perioda kocetrovaja visokih ili iskih aktivosti izvora saobraćaja u velikim vremeskim itervalima. Rezultat toga jeste da je korelacioa struktura stream-ova saobraćaja u savremeim mrežama u kotradiktorosti sa tradicioalim modelima, pa se Poasoovim modelima e mogu opisati empirijski utvrđea fraktala svojstva agregatog saobraćaja. Svi stream-ovi se statistički multipleksiraju sa lokalim eregularostima u čvorovima. Stoga, svaki paket ima stohastičko kašjeje i dolazi do etrivijalih stohastičkih pojava. Statistike prvog reda kašjeja paketa pružaju iformacije u pogledu dužie burst-a, a statistike drugog reda pružaju iformacije o eregularostima u okružeju i jihovoj spektraloj gustii. U takvim tokovima paketa egativa razlika između vremea međudolazaka i međuodlazaka paketa odgovara formiraju klastera, što povećava verovatoću odbacivaja paketa i utiče a parametre kvaliteta servisa. Kako bi se obezbedio potpu opis takvih procesa, eophodo je u aalizu uključiti svojstva koja se ispoljavaju u kraćim i dužim vremeskim itervalima. Frakcioa aaliza tokova paketa u različitim itervalima omogućava opis frakcioe prirode saobraćaja u mreži i kvalitativo i kvatitativo (Zaborovsky, Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

17 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 Izmerei saobraćaj u mreži je kozisteta u velikim vremeskim itervalima ili asimptotski samosliča i može se opisati jedostavim modelom sa jedim parametrom. Ovo svojstvo važi globalo i u vremeskom domeu i u obimu i ezato izmejeo u kraćim vremeskim itervalima. Dakle, za potrebe modelovaja saobraćaja eophodo je u model ugraditi mogućost lokalizacije u velikim vremeskim itervalima. Mera u kojoj se proces u mreži može smatrati fraktalim može se zapisati a sledeći ači (Zaborovsky, D L( δ) = aδ (4) gde je D fraktala dimezija, δ parametar skaliraja izabraog elemeta, a faktor koji opisuje meru. Neka je g () t fukcija kojom se opisuju svojstva procesa. Ako važi sledeća jedakost ( ϕ ) ϕ g t = g t (5) proces ima svojstvo ivarijatosti obima. Jedakost (5) omogućava defiisaje klase modela koji opisuju svojstva samosličosti. Na primer, za model gubitala diamičkog sistema važi: t ( m) ( m) m u() t = g( τ ) ( t ) q τ < tq + ξ t dτ p, ( 2ε ) t 0 gde je =,4,...2 m, a jezgro itegraljee trasformacije može se račuati kao: ( m) m, ( m) ( m) m τ ( t + ξ t) t τ < t + ξ t =. 0, ( m), ( m) m τ < t τ t + ξ t Kada p ovaj izraz postaje kovolucioi itegral iz Katovog skupa: t β u() t = Bξ ( τ) g( τ) dτ Γ( β ) β, t 0 gde je l 2 / l 3 Γ β Gama fukcija, B ξ kostata koja zavisi od karakteristika Katovog skupa. Uzimajući u obzir opšta svojstva procesa u mreži, β = fraktala dimezija, { } 2 {( ) } mb = M B t B t0 = 0, 2 DB = δb = M B t B t0 = N0 t t0 t t0, moguće je apisati: { ( 2 ) + ( 2 ) ( 2 ) + } 2 M B () t M BH t BH t BH t BH t BH t BH t 2H rh = = 2 { H } 80

18 { } = i () 2 2H 2H 2H Ako je M BH () t t k2h t = 2 t, izraz za takav proces može se zapisati kao sledeća kovolucioa jedačia: gde je: BH () t = h( t t' ) db( t' ), Γ H + 2 H /2 ',0 ' H /2 t t t t ht ( t' ) =. t t', t' < 0 t Može se uočiti da je jezgro ovog itegrala fukcija koja zadovoljava uslov (4). Tako, defiicija fraktalog svojstva ili ivarijatog od obima za kotiuali proces Z = Z() t, t T zadovoljava jedakost: d { } H (),, 0,0 Z t = a Z at t T a > < H <. Ako za svako celo p važi X Z( k ) Z( k) k = +, tada je d H ( m X = m X ) (6) Ako X ima pozitivu vredost i sredja vredost ije jedaka uli, i X i X M{ X} e može biti precizo samosliča proces u skladu sa defiicijom (6). Međutim, X M{ X} može biti asimptotski samosliča proces. Za ispitivaje fraktalih svojstava i karakteristika saobraćaja u mreži, moguće je koristiti apsolute statističke momete fukcija defiisae kao: q q m m ( m) μ ( q) = M X = M Xk. m k= Ako proces X ima svojstvo samosličosti (6), tada je vredost ( m ) μ ( q ) proporcioala ( q m β ). Stoga, za fiksu vredost q važi sledeća jedakost: log μ ( m) ( q ) β ( q ) log m C ( q ) = +. Prema (6) moguće je zapisati ( q) q( H ) β =. U ovom slučaju, svojstvo samosličosti može se formulisati kao lieara zavisost od promee ( m log ) μ ( q ). Kada proces X () t eliearo zavisi od q, koristi se Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

19 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 multifraktali kocept. Proces X () t aziva se multifraktali ako se logaritam jegovog apsolutog mometa ( m ) μ ( q ) liearo meja sa promeom logaritma agregatog ivoa m, ali e zavisi liearo od q. Za aalizu ovih svojstava ije dovoljo koristiti samo prvi ili drugi momet. U tom slučaju koriste se eparametarske statistike za agregate serije: ˆ μ ( m ) ( q) m ( m) = Xk. N / m k = ˆ ( m ) q Ako se log μ q meja liearo, proces X može se smatrati multifraktalim (Zaborovsky, Zaključak Frakcioi raču ima broje primee u mogim aučim discipliama, a u radu je posebo istakuta primea frakcioog račua za modelovaje telekomuikacioog saobraćaja. Feomei koji se sreću u savremeim telekomuikacioim mrežama, kao što su samosličost i zavisost u dugom opsegu, mogu se opisati primeom frakcioog račua. Pokazalo se da pretpostavka o mogućosti beskoačog kašjeja paketa u trazitom čvoru u mrežama sa komutacijom paketa adekvato modeluje odbacivaje paketa u čvoru. Heavy-tailed fukcija gustie verovatoće baferovaja paketa zadovoljava ovu pretpostavku. Pokazalo se da se diamika prosleđivaja paketa u virtuelim koekcijama u TCP/IP mrežama opisuje pomoću frakcioih izvoda i odgovara procesima sa zavisošću u dugom roku. Takođe, frakcioa aaliza može biti vrlo korista alat pri izučavaju lokalih i globalih feomea koji se sreću u telekomuikacioim mrežama. Literatura Babic, G., Vadalore, B., Jai, R., 998, Aalysis ad Modelig of Traffic i Moder Data Commuicatio Networks, Ohio State Uiversity Techical Report, OSU-CISRC-/98-TR02 Cheg, Y., Zhuag, W., Lig, X., 2007, Towards a FBM Model Based Network Calculus Framework with Service Differetiatio, Mobile Networks ad Applicatio, 2(5-6), pp Dalir, M., Bashour, M., 200, Applicatios of Fractioal Calculus, Applied Mathematical Scieces, 4(2), pp Devetskiotis, M., L.S. da Foseca, N., 2005, Modelig etwork traffic with log rage depedece: characterizatio, visualizatio ad tools, Computer Networks, 48(3), pp Hifler, R., 2000, Applicatio of fractioal calculus i physics, Sigapore, New Jersey, Lodo, Hog Kog, World Scietific. 82

20 Kim, H., Shroff, N., 200, Loss Probability Calculatios ad Asymptotic Aalysis for Fiite Buffer Multipleers, IEEE/ACM Trasactios o Networkig, 9(6), pp Miller, K., Ross, B., 993, A itroductio to the fractioal calculus ad fractioal differetial equatios, New York, Joh Willey ad sos. Park, C., Heradez-Campos, F., Le, L., Marrod, L.S., Park, J., Papiras, V., Smith, F.D., Smith, R.L., Trovero, M., Zhu, Z., 20, Log Rage Depedece Aalysis of IteretTraffic, Joural of Applied Statistics, 38(7), pp Park, K., Williger, W., 200, Self-Similar Network Traffic ad Performace Evaluatio, New York, Joh Willey ad sos. Smith, R.D., 20, The Dyamics of Iteret Traffic: Self-Similarity, Self-Orgaizatio, ad Comple Pheomea, Advaces i Comple Systems, 4(6), pp Zaborovsky, V., Podgurski, Y., Yegorov, S., New traffic model o the base of fractioal calculus, [Iteret], Dostupo a: < Preuzeto: Zaborovsky, V., Fractal Aalysis ad Wavelet trasform: Potetial ad Limitatio for Traffic Models Desig, [Iteret], Dostupo a: < Preuzeto: ВОЗМОЖНОСТЬ ПРИМЕНЕНИЯ МЕТОДА ДРОБНОЙ ПРОИЗВОДНОЙ ПРИ МОДЕЛИРОВАНИИ ТЕЛЕКОММУНИКАЦИОННОГО ТРАФИКА ОБЛАСТЬ: телекоммуникация ВИД СТАТЬИ: обзорная статья ЯЗЫК СТАТЬИ: сербский Резюме: Дробная производная (или производная дробного порядка) является обобщением математического понятия производной. Существует несколько разных способов обобщить это понятие, но все они совпадают с понятием обычной производной в случае натурального порядка. Данной проблемой занимались многие известные математики, среди которых Эйлер, Риман, Лиувилль, Абель и Фурье. В данной статье приведен обзор некоторых определений, наряду с основными положениями дробного исчисления, которые представлены с акцентом на возможность применения данного метода при моделировании телекоммуникационного трафика. Дело в том, что метод производной дробного порядка, нашел широкое применение в различных областях науки. Модели, основанные на данном методе, оказались востребованными в физике, механике, электротехнике, биохимии, медицине, экономике, а также в теории вероятности. В данной статье мы анализируем возможность применения метода производной дробного порядка для моделирования телекоммуникационного трафика. Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

21 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 Исследования показывают, что характеристики телекоммуникационного трафика в сети на локальном и глобальном уровнях, а также самоподобие и зависимость в широком диапазоне, могут быть эффективно описаны с помощью метода дробной производной, вместо обычной модели мультипликативного стохастического каскада. В данной статье представлены предлагаемые модели, основанные на методе дробной производной, воспроизводящие явления, присутствующие в современных телекоммуникационных сетях. Ключевые слова: моделирование, трафик, телекоммуникации, дробная производная. POSSIBILITY OF FRACTIONAL CALCULUS APPLICATION FOR TELECOMMUNICATION TRAFFIC MODELLING FIELD: Telecommuicatios ARTICLE TYPE: Review Paper ARTICLE LANGUAGE: Serbia Summary: Fractioal calculus is a field of mathematical aalysis cocered with research ad applicatio of derivatives ad itegrals of a arbitrary order. May famous mathematicias studied the theory of fractioal calculus such as Euler, Riema, Liouville, Abel, Fourier ad others. There are may proposed defiitios for calculatig derivatives ad itegrals of o-iteger order. I this paper, several proposed defiitios alog with basic statemets of fractioal calculus are preseted with a emphasis o a possibility of fractioal calculus applicatio i telecommuicatio traffic modellig. The fact is that fractioal calculus is widely used i various scietific disciplies i recet decades. Models based o fractioal calculus have proved to be very useful i physics, mechaics, electrical egieerig, biochemistry, medicie, ecoomy, ad probability theory. This paper aalyses a possibility of applicatio of fractioal calculus for modellig telecommuicatio traffic. May research studies have show that traffic characteristics at a local ad global level, such as self-similarity ad log rage depedece, ca efficietly be described by fractioal calculus istead of usig covetioal stochastic processes. Some proposed models based o fractioal calculus that describe pheomea preset i moder telecommuicatio etworks are preseted i this paper. Itroductio Measuremets ad statistic aalyses of telecommuicatio traffic have discovered that traffic i packet switched etworks shows sigificat irregularities burstiess, both i terms of traffic itesity variability ad shape of autocorrelatio fuctio. It is oticed that traffic has frac- 84

22 tioal characteristics self-similarity ad log rage depedece. As a result, a large badwidth is required, ad very ofte, this is oe of the causes of etwork iefficiecy. I compariso with covetioal models based o Poisso distributio widely used i circuit switched etworks, i models based o selfsimilarity property there are problems that are difficult to predict, measure ad cotrol i telecommuicatio traffic. Thus, measuremet, aalysis ad modellig of self-similar etwork traffic are a sort of challege. Differet research studies have show that traffic i moder telecommuicatio etworks ca efficietly be described with statistical models based o fractioal calculus. Basic statemets of fractioal calculus Although the term fractioal calculus, is actually a misomer, ad a term itegratio ad differetiatio of a arbitrary order is more suitable, fractioal calculus has bee a commo term sice l'hopital s era. Fourier also studied derivatives ad itegrals of a arbitrary order. Abel applied fractioal calculus for a itegral equatio which appears i a tautochroe problem formulatio (fid a shape of a curve for which the time take by a object slidig without frictio i uiform gravity to its lowest poit is idepedet of its startig poit). The first serious attempt to obtai a logical defiitio of a fractioal derivative belogs to Liouville. Let f be a locally itegrable o (a, ). Usually, a -fold iterated itegratio is marked as a Ι f, ad refers to as Riema-Liouville fractioal i- tegral of order of f. Defiitios of itegrals ad derivatives of a arbitrary order ca be cosolidated ito a differitegral. The process of obtaiig a differitegral is referred to as fractioal itegro-differetiatio. Some proposed models of fractioal calculus applicatio for telecommuicatio traffic modellig Recet research of telecommuicatio traffic has show that it ca efficietly be described by derivatives of a arbitrary order istead of covetioal stochastic processes. A ew iterpretatio of fractioal calculus creates ew fields of applicatio of these mathematical tools i order to achieve uderstadig of local ad global characteristics of etwork traffic. A fractioal dimesio ad log rage depedece ca be see i private, LAN ad WAN etworks. Temporal-spatial model for packet trasmissio i TCP/IP etworks Defiig Quality of Service parameters for Iteret services requires appropriate virtual coectios, depedig o traffic flows. Geerally, such a coectio cosists of several trasit odes ad liks. I moder TCP/IP etworks, delay i a give virtual coectio ca be cosidered as a costat value. A bufferig delay variatio causes a propagatio delay. The trasmissio process ca be described by equatios that ca be solved by fractioal calculus. Mikavica, B. i dr., Mogućosti primee frakcioog račua u modelovaju telekomuikacioog saobraćaja, pp

23 VOJNOTEHNIČKI GLASNIK / MILITARY TECHNICAL COURIER, 205., Vol. LXIII, No. 2 Fractal aalysis modellig of telecommuicatio traffic-potetials ad limitatios The fractal properties of telecommuicatio traffic i packet switched etworks idicate the eistece of periods of low or high activity cocetratios i log time itervals. As a result, the correlatio structure of traffic streams is i cotradictio with commo models. Therefore, the empirically cofirmed fractal properties of aggregate traffic caot be described by Poisso models. All streams are statistically multipleed with local irregularities i odes. Thus, each packet is stochastically delayed ad o-trivial processes occur. I order to achieve a complete descriptio of these processes, it is ecessary to ivolve the properties eistig i short ad log itervals i such a aalysis. The fractioal aalysis of packet flows i differet time itervals esures a qualitative ad quatitative descriptio of the telecommuicatio traffic fractioal ature. Coclusio Fractioal calculus has may applicatios. This paper is primarily cocered with fractioal calculus applicatio for telecommuicatio traffic modellig. It was show that a idefiite packet delay i the trasit ode i packet switchig etworks adequately models the droppig of the packet i the ode. The dyamics of packet trasmissio i virtual coectios i TCP/IP etworks ca be described by the equatios i fractioal derivatives ad it correspods to processes with log-ragedepedece. Also, fractioal calculus ca be a very useful tool i the research of local ad global pheomea i telecommuicatio etworks. Key words: models, traffic, telecommuicatios, fractioal calculus. Datum prijema člaka / Paper received o / Дата получения работы: Datum dostavljaja ispravki rukopisa / Mauscript correctios submitted o / Дата получения исправленной версии работы: Datum koačog prihvataja člaka za objavljivaje / Paper accepted for publishig o / Дата окончательного согласования работы:

Red veze za benzen. Slika 1.

Red veze za benzen. Slika 1. Red veze za benzen Benzen C 6 H 6 je aromatično ciklično jedinjenje. Njegove dve rezonantne forme (ili Kekuléove structure), prema teoriji valentne veze (VB) prikazuju se uobičajeno kao na slici 1 a),

More information

Projektovanje paralelnih algoritama II

Projektovanje paralelnih algoritama II Projektovanje paralelnih algoritama II Primeri paralelnih algoritama, I deo Paralelni algoritmi za množenje matrica 1 Algoritmi za množenje matrica Ovde su data tri paralelna algoritma: Direktan algoritam

More information

Geometrijsko mesto korena

Geometrijsko mesto korena Geometrijsko mesto korea U dosadašjem delu kursa su, između ostalog, bile razmatrae karakteristike SAU i povezaost tih karakteristika sa položajem polova sistema u kompleksoj ravi. Uočea je direkta zavisost

More information

On the convergence rates of Gladyshev s Hurst index estimator

On the convergence rates of Gladyshev s Hurst index estimator Noliear Aalysis: Modellig ad Cotrol, 2010, Vol 15, No 4, 445 450 O the covergece rates of Gladyshev s Hurst idex estimator K Kubilius 1, D Melichov 2 1 Istitute of Mathematics ad Iformatics, Vilius Uiversity

More information

Trougaone norme i primena u fazi skupovima

Trougaone norme i primena u fazi skupovima Sadržaj Predgovor... 3 1. Trougaoe orme i koorme... 5 1.1 Trougaoe orme... 5 1.2 Trougaoe koorme... 10 1.3 Neprekidost... 13 1.4 Algebarski aspekt... 15 1.5 Polugrupe i t-orme... 21 2. Fazi aritmetika...

More information

ZANIMLJIV NAČIN IZRAČUNAVANJA NEKIH GRANIČNIH VRIJEDNOSTI FUNKCIJA. Šefket Arslanagić, Sarajevo, BiH

ZANIMLJIV NAČIN IZRAČUNAVANJA NEKIH GRANIČNIH VRIJEDNOSTI FUNKCIJA. Šefket Arslanagić, Sarajevo, BiH MAT-KOL (Banja Luka) XXIII ()(7), -7 http://wwwimviblorg/dmbl/dmblhtm DOI: 75/МК7A ISSN 5-6969 (o) ISSN 986-588 (o) ZANIMLJIV NAČIN IZRAČUNAVANJA NEKIH GRANIČNIH VRIJEDNOSTI FUNKCIJA Šefket Arslanagić,

More information

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL

More information

Numerical Methods for Ordinary Differential Equations

Numerical Methods for Ordinary Differential Equations Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html

More information

Reliability and Queueing

Reliability and Queueing Copyright 999 Uiversity of Califoria Reliability ad Queueig by David G. Messerschmitt Supplemetary sectio for Uderstadig Networked Applicatios: A First Course, Morga Kaufma, 999. Copyright otice: Permissio

More information

TEORIJA SKUPOVA Zadaci

TEORIJA SKUPOVA Zadaci TEORIJA SKUPOVA Zadai LOGIKA 1 I. godina 1. Zapišite simbolima: ( x nije element skupa S (b) d je član skupa S () F je podskup slupa S (d) Skup S sadrži skup R 2. Neka je S { x;2x 6} = = i neka je b =

More information

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet

More information

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed

More information

VIEWPOINTS. Slavica Jovetic* s comment on Correlation analysis of indicators of regional competitiveness: The case of the Republic of Serbia (2013)

VIEWPOINTS. Slavica Jovetic* s comment on Correlation analysis of indicators of regional competitiveness: The case of the Republic of Serbia (2013) Ecoomic Horizos May - August 2014 Volume 16 Number 2 161-163 Faculty of Ecoomics Uiversity of Kragujevac UDC: 33 eissn 2217-9232 www. ekfak.kg.ac.rs VIEWPOINTS Slavica Jovetic* s commet o Correlatio aalysis

More information

1 6 = 1 6 = + Factorials and Euler s Gamma function

1 6 = 1 6 = + Factorials and Euler s Gamma function Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio

More information

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

A 2nTH ORDER LINEAR DIFFERENCE EQUATION A 2TH ORDER LINEAR DIFFERENCE EQUATION Doug Aderso Departmet of Mathematics ad Computer Sciece, Cocordia College Moorhead, MN 56562, USA ABSTRACT: We give a formulatio of geeralized zeros ad (, )-discojugacy

More information

FUZZY ALTERNATING DIRECTION IMPLICIT METHOD FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS IN THREE DIMENSIONS

FUZZY ALTERNATING DIRECTION IMPLICIT METHOD FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS IN THREE DIMENSIONS FUZZY ALTERNATING DIRECTION IMPLICIT METHOD FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS IN THREE DIMENSIONS N.Mugutha *1, B.Jessaili Jeba #2 *1 Assistat Professor, Departmet of Mathematics, M.V.Muthiah

More information

Chapter 4 : Laplace Transform

Chapter 4 : Laplace Transform 4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic

More information

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations Global Joural of Sciece Frotier Research Mathematics ad Decisio Scieces Volume 3 Issue Versio 0 Year 03 Type : Double Blid Peer Reviewed Iteratioal Research Joural Publisher: Global Jourals Ic (USA Olie

More information

Research Article A New Second-Order Iteration Method for Solving Nonlinear Equations

Research Article A New Second-Order Iteration Method for Solving Nonlinear Equations Abstract ad Applied Aalysis Volume 2013, Article ID 487062, 4 pages http://dx.doi.org/10.1155/2013/487062 Research Article A New Secod-Order Iteratio Method for Solvig Noliear Equatios Shi Mi Kag, 1 Arif

More information

KLASIFIKACIJA NAIVNI BAJES. NIKOLA MILIKIĆ URL:

KLASIFIKACIJA NAIVNI BAJES. NIKOLA MILIKIĆ   URL: KLASIFIKACIJA NAIVNI BAJES NIKOLA MILIKIĆ EMAIL: nikola.milikic@fon.bg.ac.rs URL: http://nikola.milikic.info ŠTA JE KLASIFIKACIJA? Zadatak određivanja klase kojoj neka instanca pripada instanca je opisana

More information

The Sumudu transform and its application to fractional differential equations

The Sumudu transform and its application to fractional differential equations ISSN : 30-97 (Olie) Iteratioal e-joural for Educatio ad Mathematics www.iejem.org vol. 0, No. 05, (Oct. 03), 9-40 The Sumudu trasform ad its alicatio to fractioal differetial equatios I.A. Salehbhai, M.G.

More information

VERTICAL MOVEMENTS FROM LEVELLING, GRAVITY AND GPS MEASUREMENTS

VERTICAL MOVEMENTS FROM LEVELLING, GRAVITY AND GPS MEASUREMENTS rd IAG / 2th FIG Symposium, Bade, May 22-24, 26 VERTICAL MOVEMENTS FROM LEVELLING, GRAVITY AND GPS MEASUREMENTS N. Hatjidakis, D. Rossikopoulos Departmet of Geodesy ad Surveyig, Faculty of Egieerig Aristotle

More information

Introduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT

Introduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory

More information

Modified Logistic Maps for Cryptographic Application

Modified Logistic Maps for Cryptographic Application Applied Mathematics, 25, 6, 773-782 Published Olie May 25 i SciRes. http://www.scirp.org/joural/am http://dx.doi.org/.4236/am.25.6573 Modified Logistic Maps for Cryptographic Applicatio Shahram Etemadi

More information

Average Number of Real Zeros of Random Fractional Polynomial-II

Average Number of Real Zeros of Random Fractional Polynomial-II Average Number of Real Zeros of Radom Fractioal Polyomial-II K Kadambavaam, PG ad Research Departmet of Mathematics, Sri Vasavi College, Erode, Tamiladu, Idia M Sudharai, Departmet of Mathematics, Velalar

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio

More information

Science & Technologies COMMUTATIONAL PROPERTIES OF OPERATORS OF MIXED TYPE PRESERVING THE POWERS - I

Science & Technologies COMMUTATIONAL PROPERTIES OF OPERATORS OF MIXED TYPE PRESERVING THE POWERS - I COMMUTATIONAL PROPERTIES OF OPERATORS OF MIXED TYPE PRESERVING TE POWERS - I Miryaa S. ristova Uiversity of Natioal ad World Ecoomy, Deartmet of Mathematics Studetsi Grad "risto Botev", 17 Sofia, BULGARIA

More information

Algoritam za množenje ulančanih matrica. Alen Kosanović Prirodoslovno-matematički fakultet Matematički odsjek

Algoritam za množenje ulančanih matrica. Alen Kosanović Prirodoslovno-matematički fakultet Matematički odsjek Algoritam za množenje ulančanih matrica Alen Kosanović Prirodoslovno-matematički fakultet Matematički odsjek O problemu (1) Neka je A 1, A 2,, A n niz ulančanih matrica duljine n N, gdje su dimenzije matrice

More information

NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS THERMAL SCIENCE, Year 07, Vol., No. 4, pp. 595-599 595 NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS by Yula WANG *, Da TIAN, ad Zhiyua LI Departmet of Mathematics,

More information

Gamma Distribution and Gamma Approximation

Gamma Distribution and Gamma Approximation Gamma Distributio ad Gamma Approimatio Xiaomig Zeg a Fuhua (Frak Cheg b a Xiame Uiversity, Xiame 365, Chia mzeg@jigia.mu.edu.c b Uiversity of Ketucky, Leigto, Ketucky 456-46, USA cheg@cs.uky.edu Abstract

More information

Structural Functionality as a Fundamental Property of Boolean Algebra and Base for Its Real-Valued Realizations

Structural Functionality as a Fundamental Property of Boolean Algebra and Base for Its Real-Valued Realizations Structural Fuctioality as a Fudametal Property of Boolea Algebra ad Base for Its Real-Valued Realizatios Draga G. Radojević Uiversity of Belgrade, Istitute Mihajlo Pupi, Belgrade draga.radojevic@pupi.rs

More information

μ are complex parameters. Other

μ are complex parameters. Other A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.

More information

Reconstruction of the Volterra-type integro-differential operator from nodal points

Reconstruction of the Volterra-type integro-differential operator from nodal points Keski Boudary Value Problems 18 18:47 https://doi.org/1.1186/s13661-18-968- R E S E A R C H Ope Access Recostructio of the Volterra-type itegro-differetial operator from odal poits Baki Keski * * Correspodece:

More information

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty

More information

Lecture 27: Optimal Estimators and Functional Delta Method

Lecture 27: Optimal Estimators and Functional Delta Method Stat210B: Theoretical Statistics Lecture Date: April 19, 2007 Lecture 27: Optimal Estimators ad Fuctioal Delta Method Lecturer: Michael I. Jorda Scribe: Guilherme V. Rocha 1 Achievig Optimal Estimators

More information

Advanced Stochastic Processes.

Advanced Stochastic Processes. Advaced Stochastic Processes. David Gamarik LECTURE 2 Radom variables ad measurable fuctios. Strog Law of Large Numbers (SLLN). Scary stuff cotiued... Outlie of Lecture Radom variables ad measurable fuctios.

More information

Research Article Some E-J Generalized Hausdorff Matrices Not of Type M

Research Article Some E-J Generalized Hausdorff Matrices Not of Type M Abstract ad Applied Aalysis Volume 2011, Article ID 527360, 5 pages doi:10.1155/2011/527360 Research Article Some E-J Geeralized Hausdorff Matrices Not of Type M T. Selmaogullari, 1 E. Savaş, 2 ad B. E.

More information

OPTIMAL STOPPING AND EXIT TIMES FOR SOME CLASSES OF RANDOM PROCESSES. Vladyslav Tomashyk

OPTIMAL STOPPING AND EXIT TIMES FOR SOME CLASSES OF RANDOM PROCESSES. Vladyslav Tomashyk NATIONAL TARAS SHEVCHENKO UNIVERSITY OF KYIV UKRAINE OPTIMAL STOPPING AND EXIT TIMES FOR SOME CLASSES OF RANDOM PROCESSES Vladyslav Tomashyk Mechaics ad Mathematics Faculty Departmet of Probability Theory,

More information

On the Variations of Some Well Known Fixed Point Theorem in Metric Spaces

On the Variations of Some Well Known Fixed Point Theorem in Metric Spaces Turkish Joural of Aalysis ad Number Theory, 205, Vol 3, No 2, 70-74 Available olie at http://pubssciepubcom/tjat/3/2/7 Sciece ad Educatio Publishig DOI:0269/tjat-3-2-7 O the Variatios of Some Well Kow

More information

Generalized Semi- Markov Processes (GSMP)

Generalized Semi- Markov Processes (GSMP) Geeralized Semi- Markov Processes (GSMP) Summary Some Defiitios Markov ad Semi-Markov Processes The Poisso Process Properties of the Poisso Process Iterarrival times Memoryless property ad the residual

More information

Chapter 4. Fourier Series

Chapter 4. Fourier Series Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,

More information

wavelet collocation method for solving integro-differential equation.

wavelet collocation method for solving integro-differential equation. IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma

More information

DISTRIBUTION LAW Okunev I.V.

DISTRIBUTION LAW Okunev I.V. 1 DISTRIBUTION LAW Okuev I.V. Distributio law belogs to a umber of the most complicated theoretical laws of mathematics. But it is also a very importat practical law. Nothig ca help uderstad complicated

More information

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS http://www.paper.edu.c Iteratioal Joural of Bifurcatio ad Chaos, Vol. 1, No. 5 () 119 15 c World Scietific Publishig Compay AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC

More information

Chapter 9: Numerical Differentiation

Chapter 9: Numerical Differentiation 178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical

More information

Zadatci sa ciklusima. Zadatak1: Sastaviti progra koji određuje z ir prvih prirod ih rojeva.

Zadatci sa ciklusima. Zadatak1: Sastaviti progra koji određuje z ir prvih prirod ih rojeva. Zadatci sa ciklusima Zadatak1: Sastaviti progra koji određuje z ir prvih prirod ih rojeva. StrToIntDef(tekst,broj) - funkcija kojom se tekst pretvara u ceo broj s tim da je uvedena automatska kontrola

More information

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu

More information

Uvod u relacione baze podataka

Uvod u relacione baze podataka Uvod u relacione baze podataka Ana Spasić 2. čas 1 Mala studentska baza dosije (indeks, ime, prezime, datum rodjenja, mesto rodjenja, datum upisa) predmet (id predmeta, sifra, naziv, bodovi) ispitni rok

More information

RAINFALL PREDICTION BY WAVELET DECOMPOSITION

RAINFALL PREDICTION BY WAVELET DECOMPOSITION RAIFALL PREDICTIO BY WAVELET DECOMPOSITIO A. W. JAYAWARDEA Departmet of Civil Egieerig, The Uiversit of Hog Kog, Hog Kog, Chia P. C. XU Academ of Mathematics ad Sstem Scieces, Chiese Academ of Scieces,

More information

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

Generating Functions for Laguerre Type Polynomials. Group Theoretic method It. Joural of Math. Aalysis, Vol. 4, 2010, o. 48, 257-266 Geeratig Fuctios for Laguerre Type Polyomials α of Two Variables L ( xy, ) by Usig Group Theoretic method Ajay K. Shula* ad Sriata K. Meher** *Departmet

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Analysis of composites with multiple rigid-line reinforcements by the BEM

Analysis of composites with multiple rigid-line reinforcements by the BEM Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100

More information

Fall 2013 MTH431/531 Real analysis Section Notes

Fall 2013 MTH431/531 Real analysis Section Notes Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters

More information

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations . Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a

More information

THE REPRESENTATION OF THE REMAINDER IN CLASSICAL BERNSTEIN APPROXIMATION FORMULA

THE REPRESENTATION OF THE REMAINDER IN CLASSICAL BERNSTEIN APPROXIMATION FORMULA Global Joural of Advaced Research o Classical ad Moder Geometries ISSN: 2284-5569, Vol.6, 2017, Issue 2, pp.119-125 THE REPRESENTATION OF THE REMAINDER IN CLASSICAL BERNSTEIN APPROXIMATION FORMULA DAN

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

SOME ASPECTS OF THE STIC SYSTEM STABILITY CALCULATION 1 UDC : (045)

SOME ASPECTS OF THE STIC SYSTEM STABILITY CALCULATION 1 UDC : (045) FACTA UNIVERSITATIS Series: Architecture ad Civil Egieerig Vol. 7, N o 1, 9, pp. 35-41 DOI: 1.98/FUACE9135B SOME ASPECTS OF THE STIC SYSTEM STABIITY CACUATION 1 UDC 64.46:64.73.5(45) Emra Bujar 1, Dragoslav

More information

Control chart for number of customers in the system of M [X] / M / 1 Queueing system

Control chart for number of customers in the system of M [X] / M / 1 Queueing system Iteratioal Joural of Iovative Research i Sciece, Egieerig ad Techology (A ISO 3297: 07 Certified Orgaiatio) Cotrol chart for umber of customers i the system of M [X] / M / Queueig system T.Poogodi, Dr.

More information

Vector Quantization: a Limiting Case of EM

Vector Quantization: a Limiting Case of EM . Itroductio & defiitios Assume that you are give a data set X = { x j }, j { 2,,, }, of d -dimesioal vectors. The vector quatizatio (VQ) problem requires that we fid a set of prototype vectors Z = { z

More information

Research Article A Note on Ergodicity of Systems with the Asymptotic Average Shadowing Property

Research Article A Note on Ergodicity of Systems with the Asymptotic Average Shadowing Property Discrete Dyamics i Nature ad Society Volume 2011, Article ID 360583, 6 pages doi:10.1155/2011/360583 Research Article A Note o Ergodicity of Systems with the Asymptotic Average Shadowig Property Risog

More information

Orthogonal Gaussian Filters for Signal Processing

Orthogonal Gaussian Filters for Signal Processing Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios

More information

arxiv: v1 [math.nt] 5 Jan 2017 IBRAHIM M. ALABDULMOHSIN

arxiv: v1 [math.nt] 5 Jan 2017 IBRAHIM M. ALABDULMOHSIN FRACTIONAL PARTS AND THEIR RELATIONS TO THE VALUES OF THE RIEMANN ZETA FUNCTION arxiv:70.04883v [math.nt 5 Ja 07 IBRAHIM M. ALABDULMOHSIN Kig Abdullah Uiversity of Sciece ad Techology (KAUST, Computer,

More information

Introductions to HarmonicNumber2

Introductions to HarmonicNumber2 Itroductios to HarmoicNumber2 Itroductio to the differetiated gamma fuctios Geeral Almost simultaeously with the developmet of the mathematical theory of factorials, biomials, ad gamma fuctios i the 8th

More information

Lecture 19: Convergence

Lecture 19: Convergence Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may

More information

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad

More information

ON SOME TRIGONOMETRIC POWER SUMS

ON SOME TRIGONOMETRIC POWER SUMS IJMMS 0: 2002 185 191 PII. S016117120200771 http://ijmms.hidawi.com Hidawi Publishig Corp. ON SOME TRIGONOMETRIC POWER SUMS HONGWEI CHEN Received 17 Jue 2001 Usig the geeratig fuctio method, the closed

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

577. Estimation of surface roughness using high frequency vibrations

577. Estimation of surface roughness using high frequency vibrations 577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia

More information

Basis for simulation techniques

Basis for simulation techniques Basis for simulatio techiques M. Veeraraghava, March 7, 004 Estimatio is based o a collectio of experimetal outcomes, x, x,, x, where each experimetal outcome is a value of a radom variable. x i. Defiitios

More information

Research Article Nonautonomous Discrete Neuron Model with Multiple Periodic and Eventually Periodic Solutions

Research Article Nonautonomous Discrete Neuron Model with Multiple Periodic and Eventually Periodic Solutions Discrete Dyamics i Nature ad Society Volume 21, Article ID 147282, 6 pages http://dx.doi.org/1.11/21/147282 Research Article Noautoomous Discrete Neuro Model with Multiple Periodic ad Evetually Periodic

More information

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,

More information

Numerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials

Numerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials IOSR Joural of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 11, Issue 6 Ver. IV (Nov. - Dec. 15), PP 1-11 www.iosrjourals.org Numerical Solutios of Secod Order Boudary Value Problems

More information

Fourier Series and their Applications

Fourier Series and their Applications Fourier Series ad their Applicatios The fuctios, cos x, si x, cos x, si x, are orthogoal over (, ). m cos mx cos xdx = m = m = = cos mx si xdx = for all m, { m si mx si xdx = m = I fact the fuctios satisfy

More information

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System The Mathematical Model ad the Simulatio Modellig Algoritm of the Multitiered Mechaical System Demi Aatoliy, Kovalev Iva Dept. of Optical Digital Systems ad Techologies, The St. Petersburg Natioal Research

More information

Fourier-Laplace Transforms of Some Special Functions

Fourier-Laplace Transforms of Some Special Functions Fourier-Lalace Trasforms of Some Secial Fuctios V. D. Sharma ad A. N. Ragari Mathematics Deartmet, Arts, Commerce ad Sciece College, Amravati- 444606(M.S), Idia. vdsharma@hotmail.co.i Mathematics Deartmet,

More information

Subject: Differential Equations & Mathematical Modeling-III

Subject: Differential Equations & Mathematical Modeling-III Power Series Solutios of Differetial Equatios about Sigular poits Subject: Differetial Equatios & Mathematical Modelig-III Lesso: Power series solutios of differetial equatios about Sigular poits Lesso

More information

MDIV. Multiple divisor functions

MDIV. Multiple divisor functions MDIV. Multiple divisor fuctios The fuctios τ k For k, defie τ k ( to be the umber of (ordered factorisatios of ito k factors, i other words, the umber of ordered k-tuples (j, j 2,..., j k with j j 2...

More information

Teaching Mathematics Concepts via Computer Algebra Systems

Teaching Mathematics Concepts via Computer Algebra Systems Iteratioal Joural of Mathematics ad Statistics Ivetio (IJMSI) E-ISSN: 4767 P-ISSN: - 4759 Volume 4 Issue 7 September. 6 PP-- Teachig Mathematics Cocepts via Computer Algebra Systems Osama Ajami Rashaw,

More information

Inverse Nodal Problems for Differential Equation on the Half-line

Inverse Nodal Problems for Differential Equation on the Half-line Australia Joural of Basic ad Applied Scieces, 3(4): 4498-4502, 2009 ISSN 1991-8178 Iverse Nodal Problems for Differetial Equatio o the Half-lie 1 2 3 A. Dabbaghia, A. Nematy ad Sh. Akbarpoor 1 Islamic

More information

Advanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology

Advanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4

More information

Application of Homotopy Perturbation Method for the Large Angle period of Nonlinear Oscillator

Application of Homotopy Perturbation Method for the Large Angle period of Nonlinear Oscillator Applicatio of Homotopy Perturbatio Method for the Large Agle period of Noliear Oscillator Olayiwola, M. O. Gbolagade A.W., Adesaya A.O. & Akipelu F.O. Departmet of Mathematical Scieces, Faculty of Sciece,

More information

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity

More information

THE POTENTIALS METHOD FOR A CLOSED QUEUEING SYSTEM WITH HYSTERETIC STRATEGY OF THE SERVICE TIME CHANGE

THE POTENTIALS METHOD FOR A CLOSED QUEUEING SYSTEM WITH HYSTERETIC STRATEGY OF THE SERVICE TIME CHANGE Joural of Applied Mathematics ad Computatioal Mechaics 5, 4(), 3-43 www.amcm.pcz.pl p-issn 99-9965 DOI:.75/amcm.5..4 e-issn 353-588 THE POTENTIALS METHOD FOR A CLOSED QUEUEING SYSTEM WITH HYSTERETIC STRATEGY

More information

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A

More information

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS R775 Philips Res. Repts 26,414-423, 1971' THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS by H. W. HANNEMAN Abstract Usig the law of propagatio of errors, approximated

More information

Series with Central Binomial Coefficients, Catalan Numbers, and Harmonic Numbers

Series with Central Binomial Coefficients, Catalan Numbers, and Harmonic Numbers 3 47 6 3 Joural of Iteger Sequeces, Vol. 5 (0), Article..7 Series with Cetral Biomial Coefficiets, Catala Numbers, ad Harmoic Numbers Khristo N. Boyadzhiev Departmet of Mathematics ad Statistics Ohio Norther

More information

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Discrete Orthogonal Moment Features Using Chebyshev Polynomials Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical

More information

Chapter 7 Isoperimetric problem

Chapter 7 Isoperimetric problem Chapter 7 Isoperimetric problem Recall that the isoperimetric problem (see the itroductio its coectio with ido s proble) is oe of the most classical problem of a shape optimizatio. It ca be formulated

More information

Council for Innovative Research

Council for Innovative Research ABSTRACT ON ABEL CONVERGENT SERIES OF FUNCTIONS ERDAL GÜL AND MEHMET ALBAYRAK Yildiz Techical Uiversity, Departmet of Mathematics, 34210 Eseler, Istabul egul34@gmail.com mehmetalbayrak12@gmail.com I this

More information

Random Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A.

Random Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A. Radom Walks o Discrete ad Cotiuous Circles by Jeffrey S. Rosethal School of Mathematics, Uiversity of Miesota, Mieapolis, MN, U.S.A. 55455 (Appeared i Joural of Applied Probability 30 (1993), 780 789.)

More information

A) is empty. B) is a finite set. C) can be a countably infinite set. D) can be an uncountable set.

A) is empty. B) is a finite set. C) can be a countably infinite set. D) can be an uncountable set. M.A./M.Sc. (Mathematics) Etrace Examiatio 016-17 Max Time: hours Max Marks: 150 Istructios: There are 50 questios. Every questio has four choices of which exactly oe is correct. For correct aswer, 3 marks

More information

Evaluation of Some Non-trivial Integrals from Finite Products and Sums

Evaluation of Some Non-trivial Integrals from Finite Products and Sums Turkish Joural of Aalysis umber Theory 6 Vol. o. 6 7-76 Available olie at http://pubs.sciepub.com/tjat//6/5 Sciece Educatio Publishig DOI:.69/tjat--6-5 Evaluatio of Some o-trivial Itegrals from Fiite Products

More information

The Riemann Zeta Function

The Riemann Zeta Function Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we

More information

Research Article A Unified Weight Formula for Calculating the Sample Variance from Weighted Successive Differences

Research Article A Unified Weight Formula for Calculating the Sample Variance from Weighted Successive Differences Discrete Dyamics i Nature ad Society Article ID 210761 4 pages http://dxdoiorg/101155/2014/210761 Research Article A Uified Weight Formula for Calculatig the Sample Variace from Weighted Successive Differeces

More information

TECHNIQUES OF INTEGRATION

TECHNIQUES OF INTEGRATION 7 TECHNIQUES OF INTEGRATION Simpso s Rule estimates itegrals b approimatig graphs with parabolas. Because of the Fudametal Theorem of Calculus, we ca itegrate a fuctio if we kow a atiderivative, that is,

More information

ECON 3150/4150, Spring term Lecture 3

ECON 3150/4150, Spring term Lecture 3 Itroductio Fidig the best fit by regressio Residuals ad R-sq Regressio ad causality Summary ad ext step ECON 3150/4150, Sprig term 2014. Lecture 3 Ragar Nymoe Uiversity of Oslo 21 Jauary 2014 1 / 30 Itroductio

More information

NAJČEŠĆE KORIŠĆENE BAZISNE FUNKCIJE U IMPULSNOM PROSTORU

NAJČEŠĆE KORIŠĆENE BAZISNE FUNKCIJE U IMPULSNOM PROSTORU UNIVERZITET U NIŠU PRIRODNO-MATEMATIČKI FAKULTET DEPARTMAN ZA FIZIKU NAJČEŠĆE KORIŠĆENE BAZISNE FUNKCIJE U IMPULSNOM PROSTORU M A S T E R R A D Metor: Dr Iva D. Mačev Studet: Marko J. Račić Niš, 0. godia

More information

Exact Solutions of the Generalized Benjamin Equation and (3 + 1)- Dimensional Gkp Equation by the Extended Tanh Method

Exact Solutions of the Generalized Benjamin Equation and (3 + 1)- Dimensional Gkp Equation by the Extended Tanh Method Available at http://pvamuedu/aam Appl Appl Math ISSN: 93-9466 Vol 7, Issue (Jue 0), pp 75 87 Applicatios ad Applied Mathematics: A Iteratioal Joural (AAM) Exact Solutios of the Geeralized Bejami Equatio

More information

Bernoulli numbers and the Euler-Maclaurin summation formula

Bernoulli numbers and the Euler-Maclaurin summation formula Physics 6A Witer 006 Beroulli umbers ad the Euler-Maclauri summatio formula I this ote, I shall motivate the origi of the Euler-Maclauri summatio formula. I will also explai why the coefficiets o the right

More information

Global Stability and Controllability of Switched Boolean Networks

Global Stability and Controllability of Switched Boolean Networks Global Stability ad Cotrollability of Switched Boolea Networks LI Haitao School of Cotrol Sciece ad Egieerig, Shadog Uiversity, Jia 250061, P R Chia E-mail: haitaoli09@gmailcom Abstract: This paper ivestigates

More information