Stationary processes and equilibrium states in non-symmetric neural networks

Size: px
Start display at page:

Download "Stationary processes and equilibrium states in non-symmetric neural networks"

Transcription

1 IVESTIGACIÓ REVISTA MEXICAA DE FÍSICA 48 (4) AGOSTO 2002 Statonary processes an equlbrum states n non-symmetrc neural networks A. Castellanos Departamento e Físca, Unversa e Sonora Apo. Post. 1626, Hermosllo 83000, Son. Méxco. acastell@fsca.uson.mx L. Vana Depto. e Físca Teórca, CCMC - UAM Apo. Post. 2681, Ensenaa B.C., Méxco. laura@ccmc.unam.mx Recbo el 9 e novembre e 2001; aceptao el 4 e abrl e 2002 Statonary processes an equlbrum states are scusse n fnte separable recurrent neural networks sequental ynamcs an away from saturaton. We escrbe thermal fluctuatons of the ynamcal orer parameters orgnate as fnte sze effects of orer O 1/2 by means of ther corresponng Fokker-Planck equaton, an fn ther tme epenent probablty strbuton. We ntrouce the concept of extene entropy of fluctuatons n orer to fn a general conton to characterze statonary states n eural etworks non symmetrc nteractons. Dvergence an rotatonal of the probablty current n the space of fluctuatons are also use to fferentate between statonary an equlbrum states. Beses, algebrac contons are foun to know when statonary states can exst. The results are llustrate by analyzng a neural network a macroscopc ynamcal fxe pont but not satsfyng etale balance at mcroscopc level. Keywors: Statstcal physcs; thermoynamcs an nonlnear ynamcal systems; neural networks; fuzzy logc; artfcal ntellgence; stochastc processes; probablty theory; stochastc processes an statstcs. Se scuten los procesos estaconaros y los estaos e equlbro en rees neuronales recurrentes, fntas y separables, con námca secuencal y muy lejos e la saturacón. Por meo e la corresponente ecuacón e Fokker-Planck, escrbmos las fluctuacones térmcas e la námca e los parámetros e oren, orgnaas como efectos e tamaño fnto e oren ( 1/2 ) y encontramos la strbucón e probabla epenente el tempo. Introucmos el concepto e entropía extena e las fluctuacones para encontrar una concón general que caracterce los estaos estaconaros en Rees euronales con nteraccones no smétrcas. Tambén se utlzan la vergenca y el rotaconal e la corrente e probabla en el espaco e fluctuacones para ferencar entre estaos estaconaros y e equlbro. Aemás, se encuentran las concones algebracas para saber cuano pueen exstr los estaos estaconaros. Los resultaos son lustraos meante el análss e una re neuronal con un punto fjo en la námca macroscópca pero que no satsface el balance etallao a nvel mcroscópco. Descrptores: Físca estaístca; termonámca y sstemas námcos no lneales; rees neuronales; lógca fusa; ntelgenca artfcal; procesos estocástcos; teoría e la probabla; procesos estocástcos y estaístca. PACS: 05.; Mh; Ey; r 1. Introucton Early works n statstcal physcs of eural etworks () were concentrate n the mean fel stuy of equlbrum propertes of nfnte systems symmetrc nteractons (for a revew see Ref. 1). However, one of the basc characterstcs of s ther hgh number of barrers an valleys n the energy space, whch makes the stuy of the ynamcs essental. In the case of non symmetrc nteractons, ynamcal technques are n fact the only tool avalable, as n ths kn of systems etale balance oes not hol an the usual statstcal concepts such as free energy cannot even be efne. On the other han, moels of the mean fel type, have prectons val for nfntely large systems ( ). Takng the thermoynamcal lmt has prove to be a useful approxmaton to stuy gases, lqus an sols, where the number of elements s of the orer of ; however, the ncluson of fnte sze effects s clearly mportant when conserng fnte systems lke the bran havng neurons workng n small groups whose sze ranges from a few thousans to 10 5 neurons. These fnte sze effects have been shown to play an mportant role n both symmetrc an non symmetrc 2, 3], even affectng, n some partcular cases, the macroscopc behavour of some systems 4]. Although conserng fnte sze effects an non symmetrc nteractons s not enough to obtan a goo escrpton of bologcal, t s clearly a frst an mportant step n ths recton; other characterstcs shoul be nclue nto the moels (see Ref. 5 for a goo revew) such as spatal structure 6 9, 5], as well as short an long-range connectvty 10, 11]. However these conseratons are out of the scope of ths paper. Untl now, most work n asymmetrc neural networks has been concentrate n the stuy of sequental retreval of embee patterns 12], an n the storage capacty of those networks 13], for nfnte systems. We beleve the stuy of the stochastc behavour resultng as a fnte sze effect n both symmetrc an non symmetrc wll contrbute towars a better unerstanng of fferent aspects of the ynamcs of an therefore t eserves specal attenton.

2 STATIOARYPROCESSESADEQUILIBRIUMSTATESIO-SYMMETRICEURALETWORKS 311 The objectve of ths work s to stuy the fference between statonary an equlbrum states from the analyss of the fluctuatons orgnate as fnte sze effects. Our formalsm escrbes the ynamcs of to frst non trval orer n the system sze 15, 16], by separatng systematcally the = from the fnte contrbutons n the Kramers- Moyal expanson 17] an keepng the frst non trval orer n (1/). In ths way t s obtane a Fokker-Planck Equaton (FPE) whch can be wrtten n terms of the fluctuatons of the ynamcal orer parameters, from whose soluton t can be constructe the probablty strbuton for these fluctuatons. Ths formalsm oes not requre symmetry of the nteracton matrx, as oppose to other methos 14]; so ynamcal flow agrams for ether symmetrc or non-symmetrc nteractons can be obtane by solvng numercally the etermnstc mean fel equatons for the fluctuaton moments. In the partcular case of symmetrc nteractons the approach to equlbrum s warrante by knowng Lyapunov functons. We know that any system n thermoynamcal equlbrum must obey etale balance. However, non symmetrc nteractons can never obey etale balance; ths s reflecte on the behavour as the lack of compensaton between postve an negatve ranom fluctuatons n the orer parameters, n a way that nuces a rotaton of the probablty current n the space of the fluctuaton varables. In ths sense, the vergence an rotatonal of the probablty current are useful tools to characterze the statonary an equlbrum states. In aton to ths geometrc pcture, we obtan the algebrac contons uner whch statonary states can exst. We also ntrouce the concept of extene entropy à la Boltzmann, for asymmetrc, to characterze the statonary states. The lmtng (maxmum) value of ths functon allows us to obtan a conton that has to be satsfe by a n a statonary state; we fn out that ths conton s met, n partcular, by any symmetrc n equlbrum, so we can see the equlbrum state as a partcular case of a statonary process. These results are llustrate by conserng a pecular system prevously analyze: a a macroscopc fxe pont but unable to satsfy etale balance at the mcroscopc level 16]. The paper s organze as follows: n Sec. 2, t s presente an overvew of the general formalsm (for etals see Refs. 15, 16); n Sec. 3, t s ntrouce the concept of extene entropy n orer to scuss statstcal equlbrum an statonary processes; n Sec. 4, the lack of etale balance an ts relaton to the rotaton of the probablty current are examne to obtan a geometrc pcture of ths phenomenon an the algebrac contons for the exstence of statonary states; fnally, n Sec. 5, t s scusse a a fxe pont but out etale balance as an example to llustrate the prevous sectons. Our work can be apple to stuy fnte sze effects n any kn of fnte, separable recurrent sequental ynamcs an away from saturaton, prove that transference functons are sgmoal. Ths nclues moels ether symmetrc or asymmetrc nteractons, an even moels where synapses are state epenent 18, 19]. The analyss of synchronous upatng 19, 20] s outse the lmts of ths theory. 2. Fnte sze effects. General formalsm Let us conser a system compose by a large, but fnte, number of nterconnecte neurons σ { 1, 1}, so the vector σ(t) = (σ 1 (t), σ 2 (t),..., σ (t)) escrbes the state of the system at a gven tme. Each varable σ (t) evolves n tme by escrbng sequental stochastc algnment of the spns uner the acton of an external fel gven by J j = 1 h (σ) = j J j σ j +θ, p ξ µ A µν ξj ν (1 δ j ), µν=1 J j s the (separable) synaptc strength of the connecton gong from neuron j to neuron, θ s a response threshol, an the auto nteractons J are exclue. These nteractons store a fnte number p of quenche bnary patterns {ξ µ } chosen at ranom µ = 1,..., p an the notaton ξ = (ξ 1,..., ξp ) { 1, 1}p, or alternatvely, ξ µ = (ξ µ 1,..., ξµ ) { 1, 1}. An, as we can see, nteractons o not nee to be symmetrcal. The mcro-ynamcs of the system s efne by a master equaton for the probablty ensty p t (σ), gven by t p t(σ) = {w (F σ)p t (F σ) w (σ)p t (σ)}, (1) w (σ) = σ tanh(βh (σ))], where w (σ) s the temperature (T = β 1 ) epenent transton probablty σ (t) σ (t), an F s an operator that flps the th spn. If we efne the probablty ensty for the pattern overlaps m as P t (m) σ p t(σ) δ m m(σ)] the usual efnton for the overlaps m(σ) = (m 1 (σ),..., m p (σ)), m µ = 1 ξ µ σ the master Eq. (1) can be rewrtten n terms of the macroscopc ynamcal varables m(σ) by nsertng δm m(σ)] nto Eq. (1), an makng a Taylor expanson n powers of the vector 2σ ξ /. If we expan the resultng Kramers-Moyal equaton n powers of (1/) an keep the two leang orers 23], we obtan a FPE val on fnte tme-scales (.e. not scalng ) for the stochastc vector m(σ). We can now wrte m(σ) as the sum of a etermnstc part m (t) whch correspons to the nfnte system, an the leang orer stochastc contrbuton q(t)/ vanshng as m = m (t) + 1 q(t) +...

3 312 A.CASTELLAOSADL.VIAA we can now separate the ynamcal equatons for these two orer parameters. In ths way, m (t) s the etermnstc soluton to the Louvlle equaton gven by F µ q; t] = K µ m (t)] + ν L µν m (t)]q ν (t), (4) t m (t) = ξ tanh β ξ Am (t)+θ] ξ,θ m (t), (2) where we efne ξ,θ as an average over an ensemble of characterze by fferent realzatons of patterns. We can now obtan a FPE for the remanng stochastc part q(σ) resultng from the fnte sze effects as follows: t P t(q) = µ q µ {P t (q) F µ q; t]} + µν 2 q µ q ν {P t (q) D µν q; t]}. (3) In ths equaton the flow term F µ q; t] s gven by { K µ x] = lm ξ µ tanh βξ Ax+θ] ξ,θ } 1 ξ µ tanh β ξ Ax+θ ], (5) L µν x] = δ µν β λ ξ µ ξ λ 1 tanh 2 βξ Ax+θ] ] ξ,θ A λν. (6) It s mportant to notce that K µ (t) makes a frozen correcton to the flux, that perssts even at T = 0 an epens on the specfc pattern confguraton {ξ µ }. On the other han, the ffuson term D µν t] of the FPE, Eq. (3), s foun to be symmetrc, an gven by D µν t] = δ µν e t lm 1 ξ µ ξν σ (0) tanh βξ Am(t)+θ ] t 0 s e s t ξ µ ξ ν tanh βξ Am (s)+θ] tanh βξ Am (t)+θ] ξ,θ. (7) It s mportant to notce that the flux s lnear n q(t) an the ffuson term s nepenent from q(t), therefore the process can be entfe as a tme epenent Ornsten-Uhlenbeck process, whose formal soluton s Gaussan: P t (q) = exp { 1 2 q q t] Ξ 1 (t) q q t ] } (8) (2π) p 2 et Ξ (t) Ξ, the tme epenent correlaton matrx for the fluctuatons Ξ µν = q µ q ν t q µ t q ν t, where < > t s the average at tme t over fferent evolutons, for a gven pattern confguraton {ξ µ }. Therefore, the process can be fully characterze by ts two frst statstcal moments, whch are foun to evolve n tme accorng to the followng etermnstc equatons t q t + F µ q; t] t = 0 (9) t Ξ(t) = L(t)Ξ(t) Ξ(t)LT (t) + 2D(t). (10) In orer to arrve to explct expressons for the flux n Eq. (4) an fusson n Eq. (7) terms, for any partcular system consere, we nee to choose nepenently rawn unbase pattern components ξ µ { 1, 1} ( equal probabltes) an nepenently rawn threshols θ (from some probablty strbuton W (θ)). Afterwars, by solvng equatons (9) an (10) n terms of the flow an fusson we can construct the probablty strbuton for the fluctuatons of the ynamcal orer parameters (8) aroun mean fel trajectores gven by Eq. (2). These stochastc fluctuatons epen on the actual realzaton of store patterns, through the frozen correcton term n Eq. (5) an ther sze s of orer O ( 1/2). 3. Statstcal equlbrum an statonary processes A physcal system s sa to be n a statonary state when the correlatons of the stochastc moments of q µ, become nepenent of tme translatons, that s f q µ (t 1 + τ) q µ (t 2 + τ)...q µ (t n + τ) = q µ (t 1 ) q µ (t 2 )...q µ (t n ) for any tme τ > 0, an µ = 1,...p 21]; ths conton s satsfe n partcular by any system n thermoynamcal equlbrum. However, thermoynamcal equlbrum requres atonally the presence of etale balance, whch means that at mcroscopcal level the probablty for a rect process

4 STATIOARYPROCESSESADEQUILIBRIUMSTATESIO-SYMMETRICEURALETWORKS 313 s equal to the probablty of the reverse process, so there s no net probablty current between any two mcroscopc states of the system. The equlbrum state has strong stablty propertes: any perturbaton actng on t wll e out fast; ths behavour s reflecte n the secon law of thermoynamcs, accorng to whch equlbrum s reache n an rreversble fashon an correspons to a maxmum entropy at constant energy, or mnmum free energy at constant temperature. However, all ths scheme s no longer useful for systems non-symmetrc nteractons, as n ths case t s not possble to efne an energy functon an etale balance no longer hols, so t s not possble to talk about equlbrum. However, t s observe that these kn of systems stll evolves from less probable to more probable states, an ths evoluton leas them towars more prvlege parts of the space of states. Therefore, for a system not satsfyng etale balance we can efne an extene entropy 21] as follows: S(t) = H(t) H(t) qp t (q) lnp t (q) (11) ( Boltzmann s constant k B 1). Ths extene entropy s an extensve quantty, whch, accorng to the H-theorem, ncreases monotoncally untl t reaches ts maxmum value, prove t exsts a statonary state; an ths result s val for any system satsfyng, or not, etale balance. The contrbuton to the extene entropy, ue to fnte sze fluctuatons, at tme t s obtane by substtutng (8) nto the prevous expresson (11). In ths way we fn S q (t) = ln P t (q) q;t = 1 ln {et Ξ (t)]} + C, 2 where C s a constant. Therefore, the ncrease n entropy between two nstants t 1 an t 2 such that t 1 < t 2, s gven by S = ln { } 1 et Ξ (t2 2 )] 0, et Ξ (t 1 )] whch mples et Ξ (t 2 )] et Ξ (t 1 )]. If we combne the Wronsk entty t ln {et B]} = T r B 1 t B] Eq. (10) escrbng the tme evoluton of the secon moments, we fn t {ln et Ξ (t)]} = 2 Tr D Ξ 1 L ], whch mples that n any relaxaton process t s satsfe Tr D Ξ 1] Tr L], where the equal sgn apples when the system reaches a statonary state. Therefore, Boltzmann s H-theorem tells us that f a statonary state exsts, t satsfes the expresson T r D Ξ 1] = T r L]. We can compare ths result to the stronger conton D Ξ 1 = L, satsfe by symmetrc systems n an equlbrum state (therefore, satsfyng etale balance) 16]. 4. on etale balance, rotatng probablty current an the exstence of statonary states In ths secton t s gven another way to characterze the fference between statonary an equlbrum states, n the framework of. We also fn the contons neee for the exstence of a statonary state. The frst objectve s reache f we analyze the probablty current n the space of fluctuaton varables by rewrtng the FPE as a contnuty equaton for the probablty flux J t (q) as follows: t P t(q) + J t (q) = 0, J t (q) = J mec t (q) + J s t (q), (12) where we have separate the probablty flux n two terms, one relate to the flux, an the other to the ffuson; for an Ornsten-Uhlenbeck process, lke ths one, these are gven by J mec t (q) = P t (q) F q, t], J s t (q) = P t (q)dξ 1 q(t) q(t) ]. (13) Van Kampen 22] has note that the frst term, resultng from the presence of an external fel, can be calle purely mechancal as t has the characterstcs of a Louvlle flux, whle the secon term s purely ffusve as t s relate to an rreversble ffusve process. From expresson n Eq. (12), t shoul be clear that when the vergence of the total probablty current J t (q) = J mec t (q) + J s t (q) s zero the system s n a statonary state. If we make q /t = 0, Ξ/t = 0 n Eqs. (9-10) for long tmes, we obtan F q] = 0 = q = L 1 K (14) L Ξ + Ξ L = 2D (15) The frst of Eqs. (14) shows that, n general, the mechancal probablty current vanshes. On the other han, Eq. (15) s known as the fluctuaton-sspaton theorem (FDT), an t shows that the correlaton Ξ µν may be n a statonary state as a consequence of the compensaton between ampng, represente by the convecton matrx L (6) an ffuson, represente by D (7). In ths lmt, we also have L = I β D A. Ths state s an equlbrum state f, atonally, the probablty current s rrotatonal, that s, f J µ /q ν = J ν /q µ, for all µ, ν. But ths happens only for A µν = A νµ, whch mples symmetrc nteractons, as expecte. On the contrary, rotatonal probablty current s present when asymmetrc nteractons appear an etale balance cannot exst. Beses, t P (q) = 0 mples that J t (q) = 0, so that, for symmetrc A, there wll be etale balance n the probablty current whenever there s a compensaton between the mechancal an the sspatve currents, that s: J µ = ( J mec µ + J s ) µ = 0. q µ q µ

5 314 A.CASTELLAOSADL.VIAA In partcular, the conton A µν = A νµ s satsfe by the use of the Hebb rule or a varaton of t where each pattern {ξ µ } s gven a fferent weght n the learnng prescrpton 24]. Therefore, vergence an rotatonal operatng on the probablty current, J t (q), gve us two contons to characterze statonary an equlbrum states, respectvely, n the framework of. To complete ths scusson from the formal pont of vew, we nee to present the algebrac contons for the exstence of statonary states. Obvously, the frst conton s that the convecton matrx shoul be nonsngular. But n aton, we nee to fn the formal soluton to the autocorrelaton matrx Ξ. We start from the FDT: L Ξ + Ξ L = 2 D an by choosng the egenvectors { λ n } of the matrx L L as the bass to represent our matrxes. Then we get L L = λ n λ n, an multplyng (15) by L on the left, an by L on the rght, we obtan λ n L L Ξ L λ m + λ n L Ξ L L λ m = λ n 2 L D L λ m. From Eq. (6) we know that f the nteracton matrx satsfes A µν = ± A νµ, then the convecton matrx L s also symmetrc or antsymmetrc, respectvely. In ths case we can use that λ n L L = ±λ n λ n, rearrangng an multplyng by L 1 from the rght we arrve to the soluton Ξ n,n = λ n L D L G n L 1 λ n G n = 2 p m=1 λ m λ n ±λ n + λ m. Then we fn, as expecte, that f the nteracton matrx A µν s symmetrc, then the soluton exsts. On the other han, for an antsymmetrc nteracton matrx t s a conton for havng a statonary state that all egenvalues {λ n } shoul be fferent. For other types of nteractons t s not possble to know, n general, f the soluton to the FDT exsts (15), so each partcular problem must be analyze separately, as we wll o t n the next secton by analyzng a system prevously stue 16]. 5. Stuy of a system a fxe pont but out etale balance In ths secton we wll stuy the ynamcs of a partcular that, even though ue to ts structure t s unable to satsfy etale balance, t has a stable fxe pont n the lmt:. In ths neural network, p bnary patterns {ξ µ } = {+1, 1}, ranomly chosen equal probablty for each of ther two values, have been store accorng to a (p p) asymmetrc nteracton matrx gven by A µν = δ µν +ε δ µ1 δ ν2, so ε can be seen as a parameter breakng the symmetry. As we wll corroborate later, n the nfnte verson t only matters the probablty strbuton P (ξ µ ), whle n the fnte case t s also relevant the partcular pattern realzaton {ξ µ } (measure by the overlap between pars of patterns). As we mentone before, the tme evoluton of the etermnstc part of the ynamcal orer parameter, m µ(t), can be obtane by numercal soluton of the Louvlle equaton, Eq. (2). If ε = 0 the characterstcs of the system are as follows: there s complete symmetry between all the Louvlle equatons for the parameters m µ (2), the nfnte system has 2p fxe ponts relate to pure states (only one orer parameter s fferent to zero) gven by m µ = ±m β, µ = 1,..., p whose values are temperature epenent an are gven by the soluton of m β = tanh (βm β ). Beses, there are other unesrable stable states relate to mxtures of the orer parameters (more than one of them are smultaneously fferent to zero), smaller basns of attracton 25]. The system evolves towars one or another of the fxe ponts, epenng on the ntal contons. In the case ε 0 the symmetry between patterns s lost, an the shape of the basns an the placng of the attractors epen on the value of ε. The convecton matrx s gven by L µν = δ µν β ξ µ ξ ν 1 tanh 2 β (ξ Am + θ) ] ξ,θ +β ε δ ν2 ξ µ ξ 1 1 tanh 2 β (ξ Am + θ) ] ξ,θ an ts etermnant s fferent from zero, so the soluton to q = L 1 K exsts. For smplcty, we wll analyze the statonary states for p = 2; as generalzaton to other values of p s straghtforwar an oes not a new nsght nto the physcs of the problem. In ths case, at T = 0 for 0 < ε < 1 there exst four attractors, two of them relate to pattern {ξ 1 } an two others havng a bg overlap pattern {ξ 2 } an a small, but stll macroscopc overlap pattern {ξ 1 }; for ε > 1 there exst only the two attractors relate to pattern 1. For hgher temperatures an analytcal expresson has been foun for ε c (T ) 16], the crtcal value of ε where stablty s lost. The FDT (15) can be summarze nto three relevant equatons to be wrtten as SΦ = Λ, Φ = (Ξ 11, Ξ 12, Ξ 22 ) an Λ = (D 11, D 12, D 22 ). For the uner conseraton, S s nonsngular whenever (L 11 ) 2 = (L 22 ) 2 L 12 L 21. Snce ths conton s met, we can conclue that the statonary state exsts. Then the FPE can be wrtten as J = 0 to establsh that compensaton between ampng an ffuson occurs accorng to the analyss wrtten before. It s just a matter of algebra to stuy the fnal sze effects close to the statonary state m = m β (1, 0). Accorng to Eq. (9), at the statonary state, m µ = m β δ µ1, the average flux F s equal to zero an the value of the average fluctuatons q epens on the explct confguraton of patterns an t s weghte by the overlap between patterns: R 12 = (1/ ) ξ1 ξ2. The fxe pont m s an stable attractor when the overlap between patterns s non negatve (R 12 0), whle for the case R 12 < 0, the state of the system escapes from ths basn of attracton an evolves towars m = m β ( 1, 0) 4]; ths

6 STATIOARYPROCESSESADEQUILIBRIUMSTATESIO-SYMMETRICEURALETWORKS 315 s an mportant fnte sze effect, snce escapng s mpossble n the nfnte case (R 0). The autocorrelaton matrx s gven by the next expresson ] 1 Ξ µν = H m (T ) {δ µν 1 + δ µ1 δ ν1 2 β2 ε 2 Hm 2 (T ) + 1 } 2 βεh m (T ) δ µ1 δ ν2 + δ µ2 δ ν1 ], n whch H m (T ) = 1 (m ) 2 1 β1 (m ) 2 ], H m (0) = 0, H m (1) =, so the contrbuton to the extene entropy, ue to fnte sze fluctuatons, s gven by S q ( ) = 1 2 ln {et Ξ ]} +C = H m (T ) ] 4 β2 ε 2 Hm 2 (T ). Ths expresson allows to see that the entropy assocate to the statonary state of the system not satsfyng etale balance, S est, can be wrtten as follows S est = S eq ln + 1 ] 4 β2 ε 2 Hm 2 (T ), where S eq = (1/2) ln H m (t)]+c. The rotaton ( of the) probablty current: J (q) = P (q) βε 1 m 2 β, prouces ths atonal entropy. 6. Dscusson an conclusons In general terms, n the scentfc lterature there exsts some confuson about the terms statonarty an equlbrum an sometmes they are even treate as synonymous. In ths work we scusse the conceptual fferences between these two terms from a formal pont of vew an prove a set of tools that can be use to llustrate ths fference n terms of the analyss of fnte sze effects. For ths purpose t was ntrouce the concept of Extene Entropy à la Boltzmann, an t was shown that n any process leang towars a statonary state t s satsfe the relatonshp Tr D Ξ 1] Tr L]. Ths conton s less restrctve than the conton satsfe by a process leang to an equlbrum state, namely D Ξ 1 L, whch can only be reache by symmetrc systems satsfyng etale balance. We also showe how the lack of etale balance s relate to the rotaton of the probablty current for the fnte sze fluctuatons q(t). For the case of antsymmetrc, whch are often use to retreve patterns sequentally 12], we emonstrate that they can reach a statonary state only f ts convecton matrx s nonsngular an the egenvalues of { λ n } are all fferent. We emonstrate our results by examnng the behavour of a pecular system whose nfnte verson has fxe ponts even though these states o not satsfy etale balance, an therefore, t cannot reach an equlbrum state. We emonstrate how, f we conser fnte sze effects, these fxe ponts transform nto the wer category of attractors of the ynamcs whch keep the system movng aroun but never reachng them, as the probablty strbuton of the fnte sze fluctuatons has a rotatonal fferent from zero. Our theory escrbes fnte sze effects aroun the orer parameter values whch are exact n the thermoynamcal lmt, an therefore help us to fn out how the fnte sze of the sample mofes the macroscopc behavour of the system. These fluctuatons, orgnate by the fnte sze of the system, are specfc for each partcular system an epen on the specfc set of store patterns. As we mentone before, ths theory can be apple n prncple, to stuy fnte sze effects n any kn of fnte, separable recurrent sequental ynamcs an away from saturaton. However, we shoul keep n mn that ths theory s correct only n the low storage regme as n the ervaton of the theory t s essental the assumpton that the overlap between any par of store patterns (chosen at ranom) goes to zero n the thermoynamcal lmt. In the case of systems parallel upatng a theory to escrbe fnte sze effects has not yet been evelope. Acknowlegments Ths work was partally supporte by grant DGAPA I (UAM). 1. A. C. C. Coolen an D. Sherrngton, n Mathematcal Approaches to eural etworks, e. Taylor JG, (orth-hollan, Amsteram, 1993) B. M. Forrest, J. Phys. A21 (1988) G. A. Kohrng, J. Phys. A: Math. Gen. 23 (1990) L. Vana, A. Castellanos an A. A. C. Coolen, Lect. otes n Comp. Sc (1999) E. Domany, J. L. van Hemmen, an K. Schulten, (e), Moels of eural etworks II, (Berln: Sprnger, 1994). 6. A. Cannng an E. Garner J. Phys. A: Math. Gen. 21 (1988) A. J. oest Europhys. Lett. 6 (1988) E. Domany, W. Knzel an R. Mer J. Phys. A: Math. Gen. 22 (1989) 2081.

7 316 A.CASTELLAOSADL.VIAA 9. A. C. C. Coolen an L. Vana J. Phys. A: Math. Gen S. Skantzos an A. C. C. Coolen, J. Phys. A: Math. Gen. 33 (2000) S. Skantzos an A. C. C. Coolen, J. Phys. A: Math. Gen. 33 (2000) H. shmor, T. akamura an M. Shno, Phys. Rev. A41 (1990) P. Peretto, J. Phys. France 49 (1988) J. M. G. Amaro, an J. F. Pérez, Europhys. Lett. 5 (1988) A. Castellanos, PhD Thess CICESE, Méxco (1998). 16. A. Castellanos, A.C.C. Coolen, L. Vana J. Phys. A31 (1998) A. C. C. Coolen an Th. W. Rujgrok, Phys. Rev. A. 38 (1988) I. J. Matus, an P. Pérez, Phys. Rev. A41 (1990) F. Zertuche, R. López-Peña an H. Waelbroeck, J. Phys. A27 (1994) O. Berner, Europhys. Lett. 16 (1991) G. Van Kampen, Can. J. Phys. 39 (1961) 551; Av. Chem. Phys. 34 (1976) G. Van Kampen, Stochastc Processes n Physcs an Chemstry (2n. e. orth-hollan 1992), p Accorng to Pawula s theorem, a consstent probablstc theory s obtane n only two cases a) f only the frst non trval orer n (1/) s kept, or b) f t s mantane the whole Kramers-Moyal expanson. 24. L. Vana, J. Phys. 49 (1988) L. Vana an A. C. C. Coolen, J. Phys. I 3 (1993) 777.

Analytical classical dynamics

Analytical classical dynamics Analytcal classcal ynamcs by Youun Hu Insttute of plasma physcs, Chnese Acaemy of Scences Emal: yhu@pp.cas.cn Abstract These notes were ntally wrtten when I rea tzpatrck s book[] an were later revse to

More information

Entropy Production in Nonequilibrium Systems Described by a Fokker-Planck Equation

Entropy Production in Nonequilibrium Systems Described by a Fokker-Planck Equation Brazlan Journal of Physcs, vol. 36, no. 4A, December, 2006 1285 Entropy Proucton n Nonequlbrum Systems Descrbe by a Fokker-Planck Equaton Tâna Tomé Insttuto e Físca Unversae e São Paulo Caxa postal 66318

More information

Large-Scale Data-Dependent Kernel Approximation Appendix

Large-Scale Data-Dependent Kernel Approximation Appendix Large-Scale Data-Depenent Kernel Approxmaton Appenx Ths appenx presents the atonal etal an proofs assocate wth the man paper [1]. 1 Introucton Let k : R p R p R be a postve efnte translaton nvarant functon

More information

Yukawa Potential and the Propagator Term

Yukawa Potential and the Propagator Term PHY304 Partcle Physcs 4 Dr C N Booth Yukawa Potental an the Propagator Term Conser the electrostatc potental about a charge pont partcle Ths s gven by φ = 0, e whch has the soluton φ = Ths escrbes the

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

The Noether theorem. Elisabet Edvardsson. Analytical mechanics - FYGB08 January, 2016

The Noether theorem. Elisabet Edvardsson. Analytical mechanics - FYGB08 January, 2016 The Noether theorem Elsabet Evarsson Analytcal mechancs - FYGB08 January, 2016 1 1 Introucton The Noether theorem concerns the connecton between a certan kn of symmetres an conservaton laws n physcs. It

More information

Conservative ensembles for nonequilibrium lattice-gas systems

Conservative ensembles for nonequilibrium lattice-gas systems Eur. Phys. J. B 64, 409 414 (2008) DOI: 10.1140/epb/e2008-00156-3 THE EUROPEAN PHYSICAL JOURNAL B Conservatve ensembles for nonequlbrum lattce-gas systems M.J. e Olvera a an T. Tomé Insttuto e Físca, Unversae

More information

New Liu Estimators for the Poisson Regression Model: Method and Application

New Liu Estimators for the Poisson Regression Model: Method and Application New Lu Estmators for the Posson Regresson Moel: Metho an Applcaton By Krstofer Månsson B. M. Golam Kbra, Pär Sölaner an Ghaz Shukur,3 Department of Economcs, Fnance an Statstcs, Jönköpng Unversty Jönköpng,

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function Commun. Theor. Phys. Bejng, Chna 49 008 pp. 97 30 c Chnese Physcal Socety Vol. 49, No., February 15, 008 Hgh-Orer Hamlton s Prncple an the Hamlton s Prncple of Hgh-Orer Lagrangan Functon ZHAO Hong-Xa an

More information

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes On a one-parameter famly of Roran arrays an the weght strbuton of MDS coes Paul Barry School of Scence Waterfor Insttute of Technology Irelan pbarry@wte Patrck Ftzpatrck Department of Mathematcs Unversty

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Entropy production in irreversible systems described by a Fokker-Planck equation

Entropy production in irreversible systems described by a Fokker-Planck equation Entropy proucton n rreversble systems escrbe by a Fokker-Planck equaton Tâna Tomé an Máro J. e Olvera Insttuto e Físca Unversae e São Paulo Caxa postal 66318 05314-970 São Paulo- SP, Brazl (Date: July

More information

GENERIC CONTINUOUS SPECTRUM FOR MULTI-DIMENSIONAL QUASIPERIODIC SCHRÖDINGER OPERATORS WITH ROUGH POTENTIALS

GENERIC CONTINUOUS SPECTRUM FOR MULTI-DIMENSIONAL QUASIPERIODIC SCHRÖDINGER OPERATORS WITH ROUGH POTENTIALS GENERIC CONTINUOUS SPECTRUM FOR MULTI-DIMENSIONAL QUASIPERIODIC SCHRÖDINGER OPERATORS WITH ROUGH POTENTIALS YANG FAN AND RUI HAN Abstract. We stuy the mult-mensonal operator (H xu) n = m n = um + f(t n

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015 Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecture 3 Contnuous Systems an Fels (Chapter 3) Where Are We Now? We ve fnshe all the essentals Fnal wll cover Lectures through Last two lectures: Classcal Fel Theory Start wth wave equatons

More information

Competitive Experimentation and Private Information

Competitive Experimentation and Private Information Compettve Expermentaton an Prvate Informaton Guseppe Moscarn an Francesco Squntan Omtte Analyss not Submtte for Publcaton Dervatons for te Gamma-Exponental Moel Dervaton of expecte azar rates. By Bayes

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

The characteristic landscape equation for an AR(2) landscape

The characteristic landscape equation for an AR(2) landscape Apple Mathematcs Letters 19 (2006) 260 266 www.elsever.com/locate/aml The characterstc lanscape equaton for an AR(2) lanscape B. Dmova, J.W. Barnes,E.Popova Grauate Program n Operatons Research an Inustral

More information

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set

More information

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.

More information

arxiv:quant-ph/ v1 6 Jun 2003

arxiv:quant-ph/ v1 6 Jun 2003 Quantum jumps an entropy proucton Henz Peter Breuer Fachberech Physk, Carl von Ossetzky Unverstät, D-6111 Olenburg, Germany an Physkalsches Insttut, Unverstät Freburg, D-79104 Freburg, Germany Date: November

More information

arxiv: v1 [quant-ph] 2 Jul 2007

arxiv: v1 [quant-ph] 2 Jul 2007 Non-Markovan quantum ynamcs an the metho of correlate projecton superoperators Henz-Peter Breuer 1, 1 Physkalsches Insttut, Unverstät Freburg, Hermann-Herer-Strasse 3, D-79104 Freburg, Germany (Date: February

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Chapter 2 Transformations and Expectations. , and define f

Chapter 2 Transformations and Expectations. , and define f Revew for the prevous lecture Defnton: support set of a ranom varable, the monotone functon; Theorem: How to obtan a cf, pf (or pmf) of functons of a ranom varable; Eamples: several eamples Chapter Transformatons

More information

ENGI9496 Lecture Notes Multiport Models in Mechanics

ENGI9496 Lecture Notes Multiport Models in Mechanics ENGI9496 Moellng an Smulaton of Dynamc Systems Mechancs an Mechansms ENGI9496 Lecture Notes Multport Moels n Mechancs (New text Secton 4..3; Secton 9.1 generalzes to 3D moton) Defntons Generalze coornates

More information

Limit Theorems for A Degenerate Fixed Point Equation

Limit Theorems for A Degenerate Fixed Point Equation Dscrete Mathematcs an Theoretcal Computer Scence DMTCS vol. subm.), by the authors, Lmt Theorems for A Degenerate Fxe Pont Equaton Mchael Drmota an Markus Kuba Insttute of Dscrete Mathematcs an Geometry,

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

CENTRAL LIMIT THEORY FOR THE NUMBER OF SEEDS IN A GROWTH MODEL IN d WITH INHOMOGENEOUS POISSON ARRIVALS

CENTRAL LIMIT THEORY FOR THE NUMBER OF SEEDS IN A GROWTH MODEL IN d WITH INHOMOGENEOUS POISSON ARRIVALS The Annals of Apple Probablty 1997, Vol. 7, No. 3, 82 814 CENTRAL LIMIT THEORY FOR THE NUMBER OF SEEDS IN A GROWTH MODEL IN WITH INHOMOGENEOUS POISSON ARRIVALS By S. N. Chu 1 an M. P. Qune Hong Kong Baptst

More information

Kinematics of Fluid Motion

Kinematics of Fluid Motion Knematcs of Flu Moton R. Shankar Subramanan Department of Chemcal an Bomolecular Engneerng Clarkson Unversty Knematcs s the stuy of moton wthout ealng wth the forces that affect moton. The scusson here

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Robust Dynamic Programming for Discounted Infinite-Horizon Markov Decision Processes with Uncertain Stationary Transition Matrice

Robust Dynamic Programming for Discounted Infinite-Horizon Markov Decision Processes with Uncertain Stationary Transition Matrice roceengs of the 2007 IEEE Symposum on Approxmate Dynamc rogrammng an Renforcement Learnng (ADRL 2007) Robust Dynamc rogrammng for Dscounte Infnte-Horzon Markov Decson rocesses wth Uncertan Statonary Transton

More information

ENTROPIC QUESTIONING

ENTROPIC QUESTIONING ENTROPIC QUESTIONING NACHUM. Introucton Goal. Pck the queston that contrbutes most to fnng a sutable prouct. Iea. Use an nformaton-theoretc measure. Bascs. Entropy (a non-negatve real number) measures

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

Numerical modeling of a non-linear viscous flow in order to determine how parameters in constitutive relations influence the entropy production

Numerical modeling of a non-linear viscous flow in order to determine how parameters in constitutive relations influence the entropy production Technsche Unverstät Berln Fakultät für Verkehrs- un Maschnensysteme, Insttut für Mechank Lehrstuhl für Kontnuumsmechank un Materaltheore, Prof. W.H. Müller Numercal moelng of a non-lnear vscous flow n

More information

The non-negativity of probabilities and the collapse of state

The non-negativity of probabilities and the collapse of state The non-negatvty of probabltes and the collapse of state Slobodan Prvanovć Insttute of Physcs, P.O. Box 57, 11080 Belgrade, Serba Abstract The dynamcal equaton, beng the combnaton of Schrödnger and Louvlle

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Non-negative Matrices and Distributed Control

Non-negative Matrices and Distributed Control Non-negatve Matrces an Dstrbute Control Yln Mo July 2, 2015 We moel a network compose of m agents as a graph G = {V, E}. V = {1, 2,..., m} s the set of vertces representng the agents. E V V s the set of

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Entropy generation in a chemical reaction

Entropy generation in a chemical reaction Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

International Journal of Pure and Applied Sciences and Technology

International Journal of Pure and Applied Sciences and Technology Int. J. Pure Appl. Sc. Technol., 4() (03), pp. 5-30 Internatonal Journal of Pure and Appled Scences and Technology ISSN 9-607 Avalable onlne at www.jopaasat.n Research Paper Schrödnger State Space Matrx

More information

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION?

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? ISAHP 001, Berne, Swtzerlan, August -4, 001 WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? Masaak SHINOHARA, Chkako MIYAKE an Kekch Ohsawa Department of Mathematcal Informaton Engneerng College

More information

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON PIOTR NAYAR AND TOMASZ TKOCZ Abstract We prove a menson-free tal comparson between the Euclean norms of sums of nepenent ranom vectors

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD SIMUATION OF WAVE POPAGATION IN AN HETEOGENEOUS EASTIC OD ogéro M Saldanha da Gama Unversdade do Estado do o de Janero ua Sào Francsco Xaver 54, sala 5 A 559-9, o de Janero, Brasl e-mal: rsgama@domancombr

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Supplementary material 1

Supplementary material 1 Supplementary materal 1 to Correlated connectvty and the dstrbuton of frng rates n the neocortex by Alexe Koulakov, Tomas Hromadka, and Anthony M. Zador The emergence of log-normal dstrbuton n neural nets.

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

MHD Kelvin-Helmholtz instability in non-hydrostatic equilibrium

MHD Kelvin-Helmholtz instability in non-hydrostatic equilibrium Journal of Physcs: Conference Seres MHD Kelvn-Helmholt nstablty n non-hyrostatc equlbrum To cte ths artcle: Y Laghouat et al 7 J. Phys.: Conf. Ser. 64 Ve the artcle onlne for upates an enhancements. Ths

More information

LECTURE 8-9: THE BAKER-CAMPBELL-HAUSDORFF FORMULA

LECTURE 8-9: THE BAKER-CAMPBELL-HAUSDORFF FORMULA LECTURE 8-9: THE BAKER-CAMPBELL-HAUSDORFF FORMULA As we have seen, 1. Taylor s expanson on Le group, Y ] a(y ). So f G s an abelan group, then c(g) : G G s the entty ap for all g G. As a consequence, a()

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,

More information

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave

More information

On Liu Estimators for the Logit Regression Model

On Liu Estimators for the Logit Regression Model CESIS Electronc Workng Paper Seres Paper No. 59 On Lu Estmators for the Logt Regresson Moel Krstofer Månsson B. M. Golam Kbra October 011 The Royal Insttute of technology Centre of Excellence for Scence

More information

π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m

π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m Homework Solutons Problem In solvng ths problem, we wll need to calculate some moments of the Gaussan dstrbuton. The brute-force method s to ntegrate by parts but there s a nce trck. The followng ntegrals

More information

p(z) = 1 a e z/a 1(z 0) yi a i x (1/a) exp y i a i x a i=1 n i=1 (y i a i x) inf 1 (y Ax) inf Ax y (1 ν) y if A (1 ν) = 0 otherwise

p(z) = 1 a e z/a 1(z 0) yi a i x (1/a) exp y i a i x a i=1 n i=1 (y i a i x) inf 1 (y Ax) inf Ax y (1 ν) y if A (1 ν) = 0 otherwise Dustn Lennon Math 582 Convex Optmzaton Problems from Boy, Chapter 7 Problem 7.1 Solve the MLE problem when the nose s exponentally strbute wth ensty p(z = 1 a e z/a 1(z 0 The MLE s gven by the followng:

More information

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Explicit bounds for the return probability of simple random walk

Explicit bounds for the return probability of simple random walk Explct bouns for the return probablty of smple ranom walk The runnng hea shoul be the same as the ttle.) Karen Ball Jacob Sterbenz Contact nformaton: Karen Ball IMA Unversty of Mnnesota 4 Ln Hall, 7 Church

More information

ON THE CURENT DENSITY AND OVERTENSION SIGNS II. THE CASE OF THE MULTI-ELECTRODIC INTERFACE

ON THE CURENT DENSITY AND OVERTENSION SIGNS II. THE CASE OF THE MULTI-ELECTRODIC INTERFACE ON HE CUREN DENSIY AND OVERENSION SIGNS II. HE CASE OF HE MULI-ELECRODIC INERFACE C. Mhalcuc an S. Lupu abstract: For a spontaneous electroe reacton the entropy proucton an the current ensty across the

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

Field and Wave Electromagnetic. Chapter.4

Field and Wave Electromagnetic. Chapter.4 Fel an Wave Electromagnetc Chapter.4 Soluton of electrostatc Problems Posson s s an Laplace s Equatons D = ρ E = E = V D = ε E : Two funamental equatons for electrostatc problem Where, V s scalar electrc

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1 CHE 54 A. Panagotopoulos STATSTCAL MECHACAL ESEMBLES MCROSCOPC AD MACROSCOPC ARABLES The central queston n Statstcal Mechancs can be phrased as follows: f partcles (atoms, molecules, electrons, nucle,

More information

Homological methods in Non-commutative Geometry Tokyo, 2007/2008 1

Homological methods in Non-commutative Geometry Tokyo, 2007/2008 1 Homologcal methos n Non-commutatve Geometry Tokyo, 2007/2008 1 Lecture 2. Secon bcomplex for cyclc homology. Connes fferental. Cyclc homology an the e Rham cohomology n the HKR case. Homology of small

More information

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem 6 th Worl Congress on Structural an Multscplnary Optmzaton Ro e Janero, 30 May - 03 June 2005, Brazl Topologcal Senstvty Analyss for Three-mensonal Lnear Elastcty Problem A.A. Novotny 1, R.A. Fejóo 1,

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

ESCI 341 Atmospheric Thermodynamics Lesson 10 The Physical Meaning of Entropy

ESCI 341 Atmospheric Thermodynamics Lesson 10 The Physical Meaning of Entropy ESCI 341 Atmospherc Thermodynamcs Lesson 10 The Physcal Meanng of Entropy References: An Introducton to Statstcal Thermodynamcs, T.L. Hll An Introducton to Thermodynamcs and Thermostatstcs, H.B. Callen

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Information Acquisition in Auctions: Sealed Bids vs. Open Bids

Information Acquisition in Auctions: Sealed Bids vs. Open Bids Informaton Acquston n Auctons: Seale Bs vs. Open Bs Ángel Hernano-Vecana Unversa e Alcante Uner Revson, Do Not Crculate Aprl 1, 2005 Abstract Ths paper stues the ncentves of a ber to acqure nformaton n

More information

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals

More information

Visualization of 2D Data By Rational Quadratic Functions

Visualization of 2D Data By Rational Quadratic Functions 7659 Englan UK Journal of Informaton an Computng cence Vol. No. 007 pp. 7-6 Vsualzaton of D Data By Ratonal Quaratc Functons Malk Zawwar Hussan + Nausheen Ayub Msbah Irsha Department of Mathematcs Unversty

More information

Hard Problems from Advanced Partial Differential Equations (18.306)

Hard Problems from Advanced Partial Differential Equations (18.306) Har Problems from Avance Partal Dfferental Equatons (18.306) Kenny Kamrn June 27, 2004 1. We are gven the PDE 2 Ψ = Ψ xx + Ψ yy = 0. We must fn solutons of the form Ψ = x γ f (ξ), where ξ x/y. We also

More information

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations Physcs 178/278 - Davd Klenfeld - Wnter 2019 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons Our goal to derve the form of the abstract quanttes n rate equatons, such as synaptc

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON PIOTR NAYAR AND TOMASZ TKOCZ Abstract We prove a menson-free tal comparson between the Euclean norms of sums of nepenent ranom vectors

More information

Chapter 7: Conservation of Energy

Chapter 7: Conservation of Energy Lecture 7: Conservaton o nergy Chapter 7: Conservaton o nergy Introucton I the quantty o a subject oes not change wth tme, t means that the quantty s conserve. The quantty o that subject remans constant

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information