A DISTRIBUTED QUANTUM EVOLUTIONARY ALGORITHM WITH A NEW CYCLING OPERATOR AND ITS APPLICATION IN FRACTAL IMAGE COMPRESSION

Size: px
Start display at page:

Download "A DISTRIBUTED QUANTUM EVOLUTIONARY ALGORITHM WITH A NEW CYCLING OPERATOR AND ITS APPLICATION IN FRACTAL IMAGE COMPRESSION"

Transcription

1 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 A DISTRIBUTED QUANTUM EVOLUTIONARY ALGORITHM WITH A NEW CYCLING OPERATOR AND ITS APPLICATION IN FRACTAL IMAGE COMPRESSION Ali Nodehi, Hosein Mohamadi and Mohamad Tayarani 3 Deparmen of Compuer Engineering, Gorgan Branch, Islamic Azad Universiy alinodehi@gorganiau.ac.ir Deparmen of Compuer Engineering, Azadshahr Branch, Islamic Azad Universiy h_mohamadi@homail.com 3 Deparmen of Compuer Engineering, Souhampon Universiy, Unied Kingdom ayarani@ieee.org ABSTRACT QUANTUM EVOLUTIONARY ALGORITHM (QEA) IS A NOVEL OPTIMIZATION ALGORITHM, PROPOSED FOR COMBINATORIAL PROBLEMS LIKE KNAPSACK AND TRAP PROBLEMS. WHILE FRACTAL IMAGE COMPRESSION IS IN THE CLASS OF NP-HARD PROBLEMS AND QEA IS HIGHLY SUITABLE FOR THE CLASS OF COMBINATORIAL PROBLEMS, QEA IS NOT WIDELY USED IN FRACTAL IMAGE COMPRESSION YET. IN ORDER TO IMPROVE THE PERFORMANCE OF FRACTAL IMAGE COMPRESSION ALGORITHMS, THIS PAPER PROPOSES A DISTRIBUTED QEA WITH A NOVEL OPERATOR CALLED CYCLING QUANTUM EVOLUTIONARY ALGORITHM. IN STANDARD QEA THE DIVERSITY IN THE POPULATION DECREASES ACROSS THE GENERATIONS. DECREASING THE DIVERSITY OF THE POPULATION DECREASES THE EXPLORATION PERFORMANCE OF THE ALGORITHM AND CAUSES THE ALGORITHM TRAPPING IN THE LOCAL OPTIMA. IN THE PROPOSED ALGORITHM, THERE ARE SOME SUBPOPULATIONS SEARCHING THE SEARCH SPACE. AFTER THE SUBPOPULATIONS ARE TRAPPED IN A LOCAL OPTIMUM, THE BEST OBSERVED POSSIBLE SOLUTIONS IN THE SUBPOPULATIONS ARE EXCHANGED IN A CYCLIC MANNER. THE PROPOSED ALGORITHM IS USED IN FRACTAL IMAGE COMPRESSION AND EXPERIMENTAL RESULTS ON SEVERAL IMAGES SHOW BETTER PERFORMANCE FOR THE PROPOSED ALGORITHM THAN GENETIC ALGORITHMS AND QEA. IN COMPARISON WITH CONVENTIONAL FRACTAL IMAGE COMPRESSION, THE PROPOSED ALGORITHM FINDS A SUITABLE SOLUTION WITH MUCH LESS COMPUTATIONAL COMPLEXITY KEYWORDS QUANTUM EVOLUTIONARY ALGORITHM, OPTIMIZATION ALGORITHM, FRACTAL IMAGE COMPRESSION, GENETIC ALGORITHM. INTRODUCTION Fracals exis in naure widely. Fracal images conain few amoun of informaion, bu posses high-level of visual complexiy []. Fracal image compression is poenially a grea coding scheme since i feaures a high compression raio and good rerieved image qualiy. In recen years, many researchers have sudied and improved he fracal image encoding and have goen a lo of achievemens. In 988, he fracal image compression was firsly proposed and uilized a number of affine mappings o denoe he original image. Those mappings are ieraed convergen and heir limi is very close o he original image. In 99, Jacquin inroduced an auomaic DOI : 0.5/ijaia

2 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 encoding algorihm for he firs ime ha called baseline fracal image compression or BFC [ 6]. This mehod breaks he original image ino sub-blocks and needs o find he bes mached subblocks according o self-similariies in he image. BFC makes he fracal image encoding become a very hopeful echnique o improve he sorage schemes applied in he consumer elecronics. Alhough BFC is very charming, a grea deal of ime cos during encoding limis i o widely pracical applicaions. In order o solve his problem some improved approaches have been presened. In [6-7] auhors proposed classificaion mehods base on he feaure of domain blocks. [8] proposed a kind of neighborhood maching mehod based on spaial correlaion which makes use of he informaion of mached range blocks and effecively reduced he encoding ime. Fracal Image Compression has araced much ineres due o is high compression raio and qualiy of compressed image. Fracal Image Compression problem is in he class of NP-Hard problems. Quanum Evoluionary Algorihm is a novel opimizaion algorihm proposed for class of combinaorial problems [9]. The probabilisic represenaion of possible soluions in QEA helps he algorihm escaping from local opima. While QEA is highly suiable for class of NP- Hard problems, i is no widely used in Fracal Image Compression. Several works have focused on improving he performance of QEA in solving knapsack and some oher problems. In order o se he parameers of QEA for pracical applicaions, [0] proposes some guidelines for parameer seings. In anoher work, by leading immune conceps in Arificial Immune Sysems, [] proposes a novel algorihm ha is called he immune quanum-inspired evoluionary algorihm. In a more recen work, [] proposes a novel muli-universe parallel Immune QEA ha uses a learning mechanism. The proposed mehod mainains he populaion diversiy along wih beer convergence speed. A muli-agen QEA is inroduced in [3] and he muli-agen operaors, compeiion and cooperaion are performed o. The roaion gae is also considered as leverage o improve he performance of QEA. In [4] a new mehod of calculaing roaion angle of quanum gae is proposed ha causes a rapid convergence for QEA. In order o preserve he diversiy in populaion and empower he search abiliy of QEA, [5] proposes a novel diversiy preservaion operaor for QEA. Reference [6] proposes a sinusoid sized populaion QEA ha makes a radeoff beween exploraion and exploiaion. Several works ry o improve he algorihm of fracal image compression using Geneic algorihm. The low speed of fracal image compression blocks is way o pracical applicaion. In [7] a geneic algorihm approach is used o improve he speed of searching in fracal image compression. A new mehod for geneic fracal image compression based on an eliis model proposed in [8]. In he proposed approach he search space for finding he bes self similariy is grealy decreased. Reference [9] makes an improvemen on he fracal image coding algorihm by applying geneic algorihm. Many researches increase he speed of fracal image compression bu he qualiy of he image will decrease. In [0] he speed of fracal image compression is improved wihou significan loss of image qualiy. In order o preven he premaure convergence of GA in fracal image compression a new approach is proposed in [], which conrols he parameers of GA adapively. A schema geneic algorihm for fracal image compression is proposed in [] o find he bes self similariy in fracal image compression. In evoluionary algorihms he diversiy in he populaion is decreased as generaions advance. When he algorihm reaches is maure sae, he diversiy in he populaion is decreased o is lowes quaniy. In order o deal wih his problem, his paper proposes a disribued QEA wih a 64

3 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 cycling exchange of bes observed possible soluions. The proposed algorihm is applied on fracal image compression problem and shows beer performance for he proposed algorihm han convenional QEA and GA. This paper is organized as follows, in Secion, he heoreical basics of fracal image compression is presened. In secion 3, QEA and is represenaion is presened, Secion 4 proposes he novel disribued QEA wih he cycling operaor, in Secion 5 he proposed algorihm is applied on fracal image compression. Several experimenal resuls are performed o in Secion 6 o es he proposed algorihm, and finally Secion 7 concludes he paper.. THEORETICAL BASICS OF FRACTAL IMAGE COMPRESSION The fracal image compression is based on he local self-similariy propery and PIFS. The relaed definiions and heorems are saed as follows [5-6]: Definiion.. Le X be a meric space wih meric d x. A map w:x X is Lipschiz wih Lipschiz facor s, if here exiss a real value s such ha dx( w( x), w( y)) sdx( x, y), 0 s, x, y X () We also say ha w is conracive wih conraciviy s. Definiion.. Le X be a meric space and d x be is meric. For a poin x X and a nonempy se A X, le us firs define he disance of x o A by d ( x, A) = infd ( x, a) () y a A Then he Hausdorff disance beween A and B is defined for any nonempy ses A, Where x B X as d( A, B ) = Max( d ( A, B ), d ( B, A )) (3) d h h h ( A, B ) = sup d ( a, B ) (4) a A Le X be he se of N N gray level images. The meric is defined as he usual Euclidean disance by regarding he elemens in X as vecors of dimension N N. Le I be a given image belonging o X. The goal is o find a se of ransformaions { w, w,..., wn }, each of which is a resriced funcion and saisfies (), such ha he given image I is an approximae aracor. The se w, w,..., w } is called PIFS. The following heorem is an imporan fac for PIFS. { n Theorem.. Consider a PIFS w, w,..., wn wih w i : X X for all i. Le W =Uw i. Then here exiss a unique poin A X such ha for any poin B X Y n A = W( A ) = Lim w ( B ) (5) n The poin A in (5) is called he fixed poin or he aracor of he mapping W. Nex, he famous Collage heorem will be inroduced. Theorem.. Le { w, w,..., wn } be a PIFS wih conraciviy facor s. Le B be any nonempy compac se in X. Then we have d( A, B ) ( s ) d( B, W( B )) (6) where A is he aracor of he mapping W and d is he Hausdorff meric (3). 65

4 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 Le d ( I, W( I )) ε where e is a very small posiive real number. By he Collage heorem, one can obain ha ε d ( A, I ) (7) ε From Eq. (7), one can see ha afer a large number of ieraions, an aracor A is generaed which is sufficienly close o he given image I. For pracical implemenaion, le I be a given gray level image. The domain pool D is defined as he se of all possible blocks of size 6 6 of he image f, which makes up (56-6+) (56-6+) =5808 blocks. The range pool R is defined o be he se of all nonoverlapping blocks of size 8 8, which makes up (56/8) (56/8) = 04 blocks. For each block v from he range pool, he fracal ransformaion is consruced by searching all elemens in he domain pool D he mos similar block. Le u denoe a sub-sampled domain block which is of he same size as v. The similariy of u and v is measured using Mean Square Error (MSE) defined by MSE= J 0 I 0 [ u( i, j ) v( i, j )] The fracal ransformaion allows he dihedral ransformaion of he domain blocks, i.e., he eigh orienaions of he blocks generaed by roaing he blocks counerclockwise a angles 0, 90, 80, and 70 degrees and flipping wih respec o he line y = x, respecively. Thus, for a given block from he range pool, here are = 464,648 MSE compuaions o obain he mos similar block from he domain pool. Thus, in oal, one needs ,648 = 475,799,55 MSE compuaions o encode he whole image using his full search compression mehod. For a given range block v, he fracal ransformaion also allows he adjusmen of he conras p and he brighness q on he subsample domain block u. The similariy is measured by he quaniy d = p. uk + q v, where u k, 0 k 7,are he eigh orienaions of u. By simple opimizaion mehod, p and q can be compued direcly as (8) And [ N u, v u, v, ] p = (9) [ N u, u u, ] q = N.[ v. p u. ] (0) Where N=64. The posiion of he mos similar domain block, he conras p, he brighness q, and he orienaion k consiue he fracal code of he given range block v. In pracice, for image, 6 bis are required o represen he posiion of he domain block. Finally, as v runs over all 04 range blocks in he range pool R, he encoding process is compleed. To decode, one firs makes up he 04 affine ransformaions from he compression codes and chooses any iniial image. Nex, one performs he 04 affine ransformaions on he image o obain a new image, and hen proceeds recursively. According o Theorems. and., he sequence of images will converge. The sopping crierion of he recursion is designed according 66

5 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 o user s applicaion and he final image is he rerieved image of fracal coding. 3. QUANTUM EVOLUTIONARY ALGORITHMS QEA is inspired from he principles of quanum compuaion, and is superposiion of saes is based on q-bis, he smalles uni of informaion sored in a wo-sae quanum compuer [8,0]. A q-bi could be eiher in sae 0 or, or in any superposiion of he wo as described below: ψ = α 0 + β () Where α and β are complex number, which denoe he corresponding sae appearance probabiliy, following below consrain: α + β = () This probabilisic represenaion implies ha if here is a sysem of m q-bis, he sysem can m represen saes simulaneously. A each observaion, a q-bis quanum sae collapses o a single sae as deermined by is corresponding probabiliies. 3. Represenaion QEA uses a novel represenaion based on he above concep of q-bis. Consider ih individual in h generaion defined as an m-q-bi as below: q i α i αi αii αim = ( 3) βi βi βii βim Where α + β =, j=,,,m, m is he number of q-bis, i.e., he sring lengh of he q-bi ij ij individual, i=,,,n, n is he number of possible soluion in populaion and is generaion number of he evoluion. Since a q-bi is a probabilisic represenaion, any superposiion of saes is simulaneously represened. If here is, for insance, a hree-q-bis(m = 3) individual such as (4): 3 q = i (4) 3 3 Or alernaively, he possible saes of he individual can be represened as: i q = (5) 6 3 Noe ha he square of above numbers are rue probabiliies, i.e. he above resul means ha he probabiliies o represen he sae 000, 00, 00, 00 are /4, /8, /4 and / respecively. Consequenly, he hree-q-bis sysem of (4) has all eigh saes informaion a he same ime. Evoluionary compuing wih he q-bi represenaion has a beer characerisic of diversiy han classical approaches since i can represen superposiion of saes. Only one q-bi individual such as (4) is enough o represen eigh saes, whereas in classical represenaion eigh individuals are needed. Addiionally, along wih he convergence of he quanum individuals, he diversiy will gradually fade away and he algorihm converges

6 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 3. QEA Srucure T In he iniializaion sep of QEA, [ αij βij ] of all q0i are iniialized wih. This implies ha 0 each q-bi individual q i represens he linear superposiion of all possible saes wih equal probabiliy. The nex sep makes a se of binary insans; x i by observing Q ( ) = { q, q,..., q n } saes, where X ( ) = { x, x,..., xi,..., xn } a generaion is a random insan of q-bi populaion. Each binary insan, x i of lengh m, is formed by selecing each bi using he probabiliy of q-bi, eiher α ij or β ij of q i. Each insan x i is evaluaed o give some n measure of is finess. The iniial bes soluion b = max{ f ( xi )} is hen seleced and sored from i= among he binary insans of X(). Then, in updae Q(), quanum gaes U updae his se of q-bi individuals Q() as discussed below. This process is repeaed in a while loop unil convergence is achieved. The appropriae quanum gae is usually designed in accordance wih problems under consideraion. The pseudo-code of QEA algorihm is defined as [4]: Procedure QEA begin =0. iniialize Q(0).. make X(0) by observing he saes of Q(0). 3. evaluae X(0). 4. Sore X(0) ino B(0). Sore he bes soluion among X(0) ino b. 5. while no erminaion condiion do begin =+ 6. make X() by observing he saes of Q(-) 7. evaluae X() 8. updae Q() using Q-gaes 9. sore he bes soluions among B(-) and X() ino B() 0. sore he bes soluion among B() ino b. if global migraion condiion hen migrae b o B() globally. else if local migraion condiion hen migrae b k in B() o B() locally end end QEA has a populaion of quanum individuals ( ) q, q,..., q } Q = { n, where is generaion sep and n is he size of populaion. The QEA procedure is described as: 0 αij. In he iniializaion sep all q-bis 0 and ij i =,,..., n and j =,,..., m are iniialized wih.. I means he probabiliy of observing "0" and "" for all q-bis is equal.. In his sep he binary soluions observing Q(0). Observing xij from qubi x ij X (0) = { x 0 0 0, x,..., xn T ij ij ] [ α β 0 if R(0,) < α ij = oherwise } a generaion = 0 are creaing by is performed as below: Where R (, ), is a uniform random number generaor. 3. All soluions in X() are evaluaed wih finess funcion. 4. Sore X(0) ino B(0). Selec bes soluion among X(0) and sore i o b. (6) 68

7 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 5. The while loop is running unil erminaion condiion is saisfied. Terminaion condiion can be considered as maximum generaion condiion or convergence condiion. 6. Observing X() from Q(-). 7. Evaluae X() by finess funcion 8. Updae Q() 9, 0. Sore he bes soluions among B(-) and X() o B(). If he fies soluion among B() is fier han b hen sore he bes soluion ino b.,. If global migraion condiion is saisfied copy b o all he soluions in B(). If local migraion condiion is saisfied replace some of soluions in B() wih bes one of hem. 3.3 Quanum Gaes Assignmen The common muaion is a random disurbance of each individual, promoing exploraion while also slowing convergence. Here, he quanum bi represenaion can be simply inerpreed as a biased muaion operaor. Therefore, he curren bes individual can be used o seer he direcion of his muaion operaor, which will speed up he convergence. The evoluionary process of quanum individual is compleed hrough he sep of updae Q(). A crossover operaor, quanum roaion gae, is described below. Specifically, a q-bi individual q i is updaed by using he roaion gae U(θ) in his algorihm. The jh q-bi value of ih quanum individual in T generaion [ α β is updaed as: ij ij ] + α ij cos( θ ) sin( θ ) α ij = + (7) β ij sin( θ ) cos( θ ) β ij Where θ is roaion angle and conrols he speed of convergence and deermined from Table I. Reference [8] shows ha hese values for θ have beer performance. TABLE I. Lookup able OF θ. xi b i f ( x) f ( b) θ 0 0 false rue 0 0 false 0.0π 0 rue 0 0 false 0.0π 0 rue 0 false 0 rue 0 4. DISTRIBUTED QEA WITH CYCLING OPERATOR In QEA when he algorihm converges, all he q-individuals reach similar values and heir quanum bis converge o eiher or 0. This means ha here is lile if any exploraion during laer sages of he algorihm s convergence and ha all of he individuals search abou a single opimum before converging o i. In oher words, he algorihm canno explore he oher pars of 69

8 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 he search space well, and all of he individuals are exploiing a single opimum. This secion proposes a novel Disribued QEA wih a cycling operaor o improve he diversiy of he populaion and so he performance of he QEA. In he proposed algorihm he populaion is divided ino some dependen subpopulaions searching he search space. Since each subpopulaion searches he search space independenly, each of hem reaches a differen local opimum. Here convergence esimaion is used o deermine if he subpopulaion is rapped in a local opimum. The converged subpopulaions are considered as he rapped subpopulaions and he cycling operaor will be performed on hem. In he cycling operaor, he bes observed of he subpopulaions will be exchanged. When he bes observed possible soluions are exchanged, he q-individuals of each subpopulaion will move oward anoher local opimum. During his movemen, he q-individuals have a chance o mee new local opima and i increases he chance of possible soluions in finding global opimum. The algorihm of he proposed algorihm is briefly shown in figure. The propose algorihm is similar o excep in seps 4 and 6. In his sep he diversiy preserving operaor is performed on all he q-individuals. In sep 4 he bes possible soluion in each subpopulaion is found and he bes possible soluion hroughou he search is sored in Bl. All he q-individuals in each subpopulaion search around Bl. The disribued naure of he proposed algorihm prevens he whole populaion from rapping in local opima. Alhough here is a chance for each of he subpopulaion rapping in local opima, he saisical naure of QEA guaranees he subpopulaions searching disinc local opima. Afer rapping in local opima, he cycling operaor is applied on he subpopulaions. Here convergence esimaion is used o deermine if he subpopulaion is rapped in a local opimum. The q-individual q i converges if i saisfies he below consrain [3]: m i, k m k = α > γ Where, m is he size of q-individuals. This consrain is saisfied only when he average disance beween he q-bis of he q-individual or [ β ] T = [ 0 ] T q i (8) and he converged saes ( [ α β ] T = [ 0 ] T α ) is less han γ. In sep 6, when he subpopulaions are converged, he bes observed possible soluions are exchanged among he subpopulaion in a cyclic manner: - B S = B - For l= o S- 3- B i = B i+ Where B is he bes observed possible soluion of he h subpopulaion and S is he number of subpopulaions. 70

9 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 5. DISTRIBUTED QEA FOR FRACTAL IMAGE COMPRESSION The proposed algorihm for each range block, searches among all he domain pool o find he bes domain block and he bes ransformaion. The coding mehod for he q-individuals in he proposed mehod is as figure. P x P y P Figure. The coding of he addresses for domain blocks. P x is he horizonal posiion of domain block, P y is he verical posiion of domain block and P is he ransformaion. In he proposed approach each q-individual, has hree pars: P x shows he horizonal posiion of domain block, P y shows he verical posiion of he domain block and P shows he ransformaion. The ransformaions are he 8 ordinary ransformaions: roae 0, 90, 80, 70, flip verically, horizonally, flip relaive o 45, and relaive o 35. Each par of each soluion is a real number and is convered o an ineger number before evaluaion process. Figure. The algorihm of he propose algorihm. 7

10 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 6. EXPERIMENTAL RESULTS This secion experimens he proposed algorihm and compares he proposed algorihm wih he performance of GA in fracal image compression. The proposed algorihm is examined on images Lena, Pepper and Baboon wih he size of and gray scale. The size of range blocks is considered as 8 8 and he size of domain blocks is considered as 6 6. In order o compare he qualiy of resuls, he PSNR es is performed: Where m n is he size of image. PSNR= 0 log n m m n i= j= 55 ( f ( i, j) g( i, j) ) The crossover rae in GA is 0.8 and he probabiliy of muaion is for each allele. he parameers of QEA is se o Table I. The size of subpopulaions is se o 5 for all he experimens. The convergence crieria, γ is considered as Table II shows he experimenal resuls on he proposed algorihm and GA. The number of ieraions for GA, QEA and he proposed algorihm for all he experimens is 00. According o able II he proposed algorihm improves he performance of fracal image compression for all he experimenal resuls. (9) 7

11 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 TABLE II. Experimenal resuls on Lena, Pepper and Baboon Picure Lena Pepper Baboon Mehod Populaion MSE Size compuaions PSNR Full Search - 59,474, ,44, QEA 5 5,0, ,096, ,07, ,44, PROPOSED 5 5,0, ALGORITHM 0 4,096, ,07, ,44, GA 5 5,0, ,096, ,07, Full Search - 59,474, ,44, QEA 5 5,0, ,096, ,07, PROPOSED ALGORITHM GA 30 6,44, ,0, ,096, ,07, ,44, ,0, ,096, ,07, Full Search - 59,474, ,44, QEA 5 5,0, ,096, ,07, PROPOSED ALGORITHM GA 30 6,44, ,0, ,096, ,07, ,44, ,0, ,096, ,07,

12 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 (a) (b) (c) (d) (e) (f) Figure 3. (a) Original image Lena of size 56 56,(b) iniial image for decoder, (c) full search mehod, (d) proposed mehod, (e) GA mehod, (e) QEA mehod 7. CONCLUSION Quanum Evoluionary Algorihms are novel algorihm proposed for combinaorial opimizaion problems like knapsack problem. Since fracal image compression is in he class of NP-Hard problems, QEA is highly suiable for his problem bu is no properly applied o fracal image compression ye. This paper proposes a Disribued QEA wih a cycling Operaor for QEA in fracal image compression. In evoluionary algorihms during he search process he diversiy in he populaion is decreased, and he algorihm is rapped in he local opima. This is rue for QEA which uses a probabilisic represenaion for he individuals. Here we propose a new mehod o help he algorihm escaping from he local opima. In he proposed algorihm he populaion is divided ino some subpopulaions and a cycling operaor exchanges he bes observed possible soluions among he subpopulaions. The experimenal resuls on Lena, Pepper, and Baboon picure show an improvemen on fracal image compression. 74

13 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 ACKNOWLEDGEMENTS The auhors would like o reurn hanks he financial suppor he Islamic Azad Universiy, Gorgan Branch, Iran. REFERENCES [] Peigen H.O., Jurgens H. and Saupe D., Chao and Fracals: New Froniers of Science, Springer Verlag, New York, 99. [] Barnsley, M.F., Fracal Everywhere, Academic Press, California, 988. [3] Barnsley, M.F. and Demko, S., Ieraed funcion sysems and he global consrucion of fracals. Proceedings of he Royal Sociey, PP 43-75,985. [4] Barnsley M.F., Elon J.H. and Hardin D.P, Recurren ieraed funcion sysems Consr., Approx. 3 3, 989. [5] Jacquin A.E, coding based on a fracal heory of ieraed conracive image ransformaions, IEEE Trans. Image Process,PP 8-30, 99. [6] Fisher Y, Fracal Image Compression: Theory and Applicaion, Springer Verlag, Berlin, 995. [7] Wang Z., Zhang D. and Yu Y., Hybrid image coding based on parial fracal mapping, Signal Process. Image Commun, pp , 000. [8] Truong T.K, Kung C.M and Jeng J.H., Hsieh M.L., Fas fracal image compression using spaial correlaion, Chaos Soliions Frac, Chaos Soliions Frac, pp , 004. [9] K. Han, J. Kim, Quanum-inspired evoluionary algorihm for a class of combinaorial opimizaion, IEEE Trans. on Evoluionary Compu, vol., no. 6, 00. [0] Kuk-Hyun Han Jong-Hwan Kim "On seing he parameers of quanum-inspired evoluionary algorihm for pracical applicaion", IEEE Congress on Evoluionary Compuaion, 003. [] Ying Li Yanning Zhang Rongchun Zhao Licheng Jiao, "The immune quanum-inspired evoluionary algorihm ", Inernaional Conference on Sysems, Man and Cyberneics, 004 IEEE. [] Xiaoming You, Sheng Liu and Dianxun Shuai "On Parallel Immune Quanum Evoluionary Algorihm Based on Learning Mechanism and Is Convergence", Springer, Advanced in Naural Compuaion [3] Chaoyong Qin Jianguo Zheng and Jiyu Lai, "A Muliagen Quanum Evoluionary Algorihm for Global Numerical Opimizaion", Springer, Life Sysem Modeling and simulaion, Vol 4689, pp , 007. [4] Rui Zhang Hui Gao, "Improved Quanum Evoluionary Algorihm for Combinaorial Opimizaion Problem", IEEE Inernaional Conference on Machine Learning and Cyberneics, 007. [5] M.-H. Tayarani- N. and M.- R. Akbarzadeh-T A Cellular Srucure and Diversiy Preserving operaor in Quanum Evoluionary Algorihms. IEEE World Conference on Compuaional Inelligence 008. [6] M.-H. Tayarani- N. and M.- R. Akbarzadeh-T A Sinusoid Size Ring Srucure Quanum Evoluionary Algorihm. IEEE Inernaional Conference on Cyberneics and Inelligen Sysems Roboics, Auomaion and Mechanics 008. [7] Xiao Chen, Guang-Xi Zhu, and Yaoing Zhu, Fracal image coding mehod based on geneic algorihms Inernaional Symposium on Mulispecral Image Processing 998. [8] Mira, S.K. Murhy, C.A. Kundu, M.K Technique for fracal image compression using geneic algorihm, IEEE Trans on Image Processing Vol 7 no 4 pp [9] Lu Xun Yu Zhongqiu The applicaion of GA in fracal image compression, 3rd IEEE World Congress on Inelligen Conrol and Auomaion, 000. [0] Gafour, A. Faraoun, K. Lehireche, A Geneic fracal image compression. ACS/IEEE Inernaional Conference on Compuer Sysems and Applicaions, 003. [] Lifeng Xi, Liangbin Zhang, A Sudy of Fracal Image Compression Based on an Improved Geneic Algorihm, Inernaional Journal of Nonlinear Science Vol 3, No, pp 6-4m

14 Inernaional Journal of Arificial Inelligence & Applicaions (IJAIA), Vol.3, No., January 0 [] Ming-Sheng Wu, Jyh-Horng Jeng, and Jer-Guang Hsieh, Schema geneic algorihm for fracal image compression Elsevier Journal of Engineering Applicaions of Arificial Inelligence Vol 0, Issue 4, pp [3] K. Han, J. Kim, Quanum-Inspired Evoluionary Algorihms wih a New Terminaion Crierion, H Gae, and Two-Phase Scheme, IEEE Trans. on Evoluionary Compuaion, vol. 8, no.,

Particle Swarm Optimization

Particle Swarm Optimization Paricle Swarm Opimizaion Speaker: Jeng-Shyang Pan Deparmen of Elecronic Engineering, Kaohsiung Universiy of Applied Science, Taiwan Email: jspan@cc.kuas.edu.w 7/26/2004 ppso 1 Wha is he Paricle Swarm Opimizaion

More information

Particle Swarm Optimization Combining Diversification and Intensification for Nonlinear Integer Programming Problems

Particle Swarm Optimization Combining Diversification and Intensification for Nonlinear Integer Programming Problems Paricle Swarm Opimizaion Combining Diversificaion and Inensificaion for Nonlinear Ineger Programming Problems Takeshi Masui, Masaoshi Sakawa, Kosuke Kao and Koichi Masumoo Hiroshima Universiy 1-4-1, Kagamiyama,

More information

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints IJCSI Inernaional Journal of Compuer Science Issues, Vol 9, Issue 1, No 1, January 2012 wwwijcsiorg 18 Applying Geneic Algorihms for Invenory Lo-Sizing Problem wih Supplier Selecion under Sorage Capaciy

More information

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Space

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Space Inernaional Journal of Indusrial and Manufacuring Engineering Applying Geneic Algorihms for Invenory Lo-Sizing Problem wih Supplier Selecion under Sorage Space Vichai Rungreunganaun and Chirawa Woarawichai

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

MODULE - 9 LECTURE NOTES 2 GENETIC ALGORITHMS

MODULE - 9 LECTURE NOTES 2 GENETIC ALGORITHMS 1 MODULE - 9 LECTURE NOTES 2 GENETIC ALGORITHMS INTRODUCTION Mos real world opimizaion problems involve complexiies like discree, coninuous or mixed variables, muliple conflicing objecives, non-lineariy,

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

A Framework for Efficient Document Ranking Using Order and Non Order Based Fitness Function

A Framework for Efficient Document Ranking Using Order and Non Order Based Fitness Function A Framework for Efficien Documen Ranking Using Order and Non Order Based Finess Funcion Hazra Imran, Adii Sharan Absrac One cenral problem of informaion rerieval is o deermine he relevance of documens

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

Technical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.

Technical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models. Technical Repor Doc ID: TR--203 06-March-203 (Las revision: 23-Februar-206) On formulaing quadraic funcions in opimizaion models. Auhor: Erling D. Andersen Convex quadraic consrains quie frequenl appear

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence Supplemen for Sochasic Convex Opimizaion: Faser Local Growh Implies Faser Global Convergence Yi Xu Qihang Lin ianbao Yang Proof of heorem heorem Suppose Assumpion holds and F (w) obeys he LGC (6) Given

More information

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi

More information

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen

More information

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems Single-Pass-Based Heurisic Algorihms for Group Flexible Flow-shop Scheduling Problems PEI-YING HUANG, TZUNG-PEI HONG 2 and CHENG-YAN KAO, 3 Deparmen of Compuer Science and Informaion Engineering Naional

More information

Chapter Floating Point Representation

Chapter Floating Point Representation Chaper 01.05 Floaing Poin Represenaion Afer reading his chaper, you should be able o: 1. conver a base- number o a binary floaing poin represenaion,. conver a binary floaing poin number o is equivalen

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Longest Common Prefixes

Longest Common Prefixes Longes Common Prefixes The sandard ordering for srings is he lexicographical order. I is induced by an order over he alphabe. We will use he same symbols (,

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

4.1 Other Interpretations of Ridge Regression

4.1 Other Interpretations of Ridge Regression CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning

More information

Logic in computer science

Logic in computer science Logic in compuer science Logic plays an imporan role in compuer science Logic is ofen called he calculus of compuer science Logic plays a similar role in compuer science o ha played by calculus in he physical

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions

More information

Retrieval Models. Boolean and Vector Space Retrieval Models. Common Preprocessing Steps. Boolean Model. Boolean Retrieval Model

Retrieval Models. Boolean and Vector Space Retrieval Models. Common Preprocessing Steps. Boolean Model. Boolean Retrieval Model 1 Boolean and Vecor Space Rerieval Models Many slides in his secion are adaped from Prof. Joydeep Ghosh (UT ECE) who in urn adaped hem from Prof. Dik Lee (Univ. of Science and Tech, Hong Kong) Rerieval

More information

BU Macro BU Macro Fall 2008, Lecture 4

BU Macro BU Macro Fall 2008, Lecture 4 Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec

More information

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Západočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France

Západočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France ADAPTIVE SIGNAL PROCESSING USING MAXIMUM ENTROPY ON THE MEAN METHOD AND MONTE CARLO ANALYSIS Pavla Holejšovsá, Ing. *), Z. Peroua, Ing. **), J.-F. Bercher, Prof. Assis. ***) Západočesá Univerzia v Plzni,

More information

Single and Double Pendulum Models

Single and Double Pendulum Models Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves Rapid Terminaion Evaluaion for Recursive Subdivision of Bezier Curves Thomas F. Hain School of Compuer and Informaion Sciences, Universiy of Souh Alabama, Mobile, AL, U.S.A. Absrac Bézier curve flaening

More information

Orientation. Connections between network coding and stochastic network theory. Outline. Bruce Hajek. Multicast with lost packets

Orientation. Connections between network coding and stochastic network theory. Outline. Bruce Hajek. Multicast with lost packets Connecions beween nework coding and sochasic nework heory Bruce Hajek Orienaion On Thursday, Ralf Koeer discussed nework coding: coding wihin he nework Absrac: Randomly generaed coded informaion blocks

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

Task-based Configuration Optimization of Modular and Reconfigurable Robots using a Multi-solution Inverse Kinematics Solver

Task-based Configuration Optimization of Modular and Reconfigurable Robots using a Multi-solution Inverse Kinematics Solver CARV 2007 Task-based Configuraion Opimizaion of Modular and Reconfigurable Robos using a Muli-soluion Inverse Kinemaics Solver S. Tabandeh 1, C. Clark 2, W. Melek 3 Absrac: Modular and Reconfigurable Robos

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1 SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision

More information

The Effect of Introducing a Tag Node in Solving Scheduling Problems Using Edge Histogram Based Sampling Algorithms

The Effect of Introducing a Tag Node in Solving Scheduling Problems Using Edge Histogram Based Sampling Algorithms Proceedings of he 2003 Congress on Evoluionary Compuaion (CEC' 03) The Effec of Inroducing a Tag Node in Solving Scheduling Problems Using Edge Hisogram Based Sampling Algorihms Shigeyoshi Tsusui Deparmen

More information

Presentation Overview

Presentation Overview Acion Refinemen in Reinforcemen Learning by Probabiliy Smoohing By Thomas G. Dieerich & Didac Busques Speaer: Kai Xu Presenaion Overview Bacground The Probabiliy Smoohing Mehod Experimenal Sudy of Acion

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

References are appeared in the last slide. Last update: (1393/08/19)

References are appeared in the last slide. Last update: (1393/08/19) SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

More information

The Rosenblatt s LMS algorithm for Perceptron (1958) is built around a linear neuron (a neuron with a linear

The Rosenblatt s LMS algorithm for Perceptron (1958) is built around a linear neuron (a neuron with a linear In The name of God Lecure4: Percepron and AALIE r. Majid MjidGhoshunih Inroducion The Rosenbla s LMS algorihm for Percepron 958 is buil around a linear neuron a neuron ih a linear acivaion funcion. Hoever,

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

Ensamble methods: Bagging and Boosting

Ensamble methods: Bagging and Boosting Lecure 21 Ensamble mehods: Bagging and Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Ensemble mehods Mixure of expers Muliple base models (classifiers, regressors), each covers a differen par

More information

A Hop Constrained Min-Sum Arborescence with Outage Costs

A Hop Constrained Min-Sum Arborescence with Outage Costs A Hop Consrained Min-Sum Arborescence wih Ouage Coss Rakesh Kawara Minnesoa Sae Universiy, Mankao, MN 56001 Email: Kawara@mnsu.edu Absrac The hop consrained min-sum arborescence wih ouage coss problem

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED 0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable

More information

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information

Removing Useless Productions of a Context Free Grammar through Petri Net

Removing Useless Productions of a Context Free Grammar through Petri Net Journal of Compuer Science 3 (7): 494-498, 2007 ISSN 1549-3636 2007 Science Publicaions Removing Useless Producions of a Conex Free Grammar hrough Peri Ne Mansoor Al-A'ali and Ali A Khan Deparmen of Compuer

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Keywords Digital Infinite-Impulse Response (IIR) filter, Digital Finite-Impulse Response (FIR) filter, DE, exploratory move

Keywords Digital Infinite-Impulse Response (IIR) filter, Digital Finite-Impulse Response (FIR) filter, DE, exploratory move Volume 5, Issue 7, July 2015 ISSN: 2277 128X Inernaional Journal of Advanced Research in Compuer Science and Sofware Engineering Research Paper Available online a: www.ijarcsse.com A Hybrid Differenial

More information

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal? EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of

More information

Optimal Design of LQR Weighting Matrices based on Intelligent Optimization Methods

Optimal Design of LQR Weighting Matrices based on Intelligent Optimization Methods Opimal Design of LQR Weighing Marices based on Inelligen Opimizaion Mehods Inernaional Journal of Inelligen Informaion Processing, Volume, Number, March Opimal Design of LQR Weighing Marices based on Inelligen

More information

Entanglement and complexity of many-body wavefunctions

Entanglement and complexity of many-body wavefunctions Enanglemen and complexiy of many-body wavefuncions Frank Versraee, Universiy of Vienna Norber Schuch, Calech Ignacio Cirac, Max Planck Insiue for Quanum Opics Tobias Osborne, Univ. Hannover Overview Compuaional

More information

A Global Convergence Proof for a Class of Genetic Algorithms

A Global Convergence Proof for a Class of Genetic Algorithms A Global Convergence Proof for a Class of Geneic Algorihms Richard F. Harl * Absrac: In his paper a varian of a geneic algorihm is invesigaed which combines he advanages of geneic algorihms and of simulaed

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

A Forward-Backward Splitting Method with Component-wise Lazy Evaluation for Online Structured Convex Optimization

A Forward-Backward Splitting Method with Component-wise Lazy Evaluation for Online Structured Convex Optimization A Forward-Backward Spliing Mehod wih Componen-wise Lazy Evaluaion for Online Srucured Convex Opimizaion Yukihiro Togari and Nobuo Yamashia March 28, 2016 Absrac: We consider large-scale opimizaion problems

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

Ensemble Confidence Estimates Posterior Probability

Ensemble Confidence Estimates Posterior Probability Ensemble Esimaes Poserior Probabiliy Michael Muhlbaier, Aposolos Topalis, and Robi Polikar Rowan Universiy, Elecrical and Compuer Engineering, Mullica Hill Rd., Glassboro, NJ 88, USA {muhlba6, opali5}@sudens.rowan.edu

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel, Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

A Video Vehicle Detection Algorithm Based on Improved Adaboost Algorithm Weiguang Liu and Qian Zhang*

A Video Vehicle Detection Algorithm Based on Improved Adaboost Algorithm Weiguang Liu and Qian Zhang* A Video Vehicle Deecion Algorihm Based on Improved Adaboos Algorihm Weiguang Liu and Qian Zhang* Zhongyuan Universiy of Technology, Zhengzhou 450000, China lwg66123@163.com, 2817343431@qq.com *The corresponding

More information

Sequential Importance Resampling (SIR) Particle Filter

Sequential Importance Resampling (SIR) Particle Filter Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle

More information

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011 2229-12 School and Workshop on Marke Microsrucure: Design, Efficiency and Saisical Regulariies 21-25 March 2011 Some mahemaical properies of order book models Frederic ABERGEL Ecole Cenrale Paris Grande

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course OMP: Arificial Inelligence Fundamenals Lecure 0 Very Brief Overview Lecurer: Email: Xiao-Jun Zeng x.zeng@mancheser.ac.uk Overview This course will focus mainly on probabilisic mehods in AI We shall presen

More information

Content-Based Shape Retrieval Using Different Shape Descriptors: A Comparative Study Dengsheng Zhang and Guojun Lu

Content-Based Shape Retrieval Using Different Shape Descriptors: A Comparative Study Dengsheng Zhang and Guojun Lu Conen-Based Shape Rerieval Using Differen Shape Descripors: A Comparaive Sudy Dengsheng Zhang and Guojun Lu Gippsland School of Compuing and Informaion Technology Monash Universiy Churchill, Vicoria 3842

More information

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when

More information

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne

More information