Ordered and disordered dynamics in random networks
|
|
- Paula Debra Hill
- 5 years ago
- Views:
Transcription
1 EUROPHYSICS LETTERS 5 March 998 Eurohys. Lett., 4 (6), (998) Ordered and disordered dynamics in random networks L. Glass ( )andc. Hill 2,3 Deartment of Physiology, McGill University 3655 Drummond Street, Montreal, Quebec, Canada H3G Y6 2 Deartment of Physics, McGill University - Montreal, Quebec, Canada 3 Deartment of Physics, Cornell University - Ithaca, NY, USA (received 20 October 997; acceted in final form 29 January 998) PACS b Theory and models of chaotic systems. PACS q Lattice theory and statistics; Ising roblems. PACS Jb Muscle contraction, nerve conduction, synatic transmission, memorization, and other neurohysiological rocesses (excluding ercetion rocesses and seech). Abstract. Random Boolean networks that model genetic networks show transitions between ordered and disordered dynamics as a function of the number of inuts er element, K, and the robability,, that the truth table for a given element will have a bias for being, in the limit as the number of elements N. We analyze transitions between ordered and disordered dynamics in randomly constructed ordinary differential equation analogues of the random Boolean networks. These networks show a transition from order to chaos for finite N. Qualitative features of the dynamics in a given network can be redicted based on the comutation of the mean dimension of the subsace admitting outflows during the integration of the equations. Genetic networks have been modeled by random Boolean networks in which time is discrete and each element comutes a Boolean function based on the values of inuts to that element []. Since the number of human genes is of the order of 00,000, and each gene is idealized as either on () or off (0), the state sace for the human gene activity is huge. An order-disorder transition has been described for random Boolean networks in the limit that the number of variables, N, as a function of the number of inuts er variable, K, and the robability,, that the truth table for a given element will have a bias for being [2]-[8]. The order-disorder boundary is given by K c = 2 c ( c ), () where K c and c reresent the values of K and on the boundary [2], [4], [6]. Kauffman has argued that for a network to be biologically meaningful, it should have relatively few attractors, ( ) glass@cnd.mcgill.ca c EDP Sciences
2 600 EUROPHYSICS LETTERS and the cycle length of attractors should be comaratively short []. In real biological systems there are not clocking devices to generate synchronous udating and theoretical models are more aroriately formulated as continuous differential equations [9]-[4]. Here we resent numerical evidence for an order-disorder transition in differential equation analogues of the discrete switching networks. We also resent a robabilistic model of the dynamics, in which we show that () also alies to the continuous equations. First consider a logical network consisting of N binary variables, X i =0,, i =,...,N. Since we consider the networks as models of genes, X i reresents the activity of gene i. In other contexts, logical variables may reresents sins or voltages. The network is udated by means of the dynamical equation X i (j +)=Λ i (X i (j),x i2 (j),...,x ik (j)), i =,...,N, (2) where Λ i (X i (j),x i2 (j),...,x ik (j)) {0,}and K is the number of inuts. More comactly, we have X(j +)=Λ(X(j)), i =,...,N. (3) Thus, for any state X(j), Λ is a truth table determining X(j +). The logical structure of eq. (2) can be catured by a differential equation [9]. To a continuous variable x i (t), we associate a discrete variable X i (t), X i (t) =0 if x i (t)<0; otherwise X i (t) =. (4) For any logical network, we define an analogous differential equation, dx i dt = x i + λ i (X i (j),x i2 (j),...,x ik (j)), i =,...,N, (5) where λ i (X i (j),x i2 (j),...,x ik (j)) is a scalar whose sign is negative (ositive) if the corresonding logical variable Λ i (X i (j),x i2 (j),...,x ik (j)) is 0 (). For each variable, the temoral evolution is governed by a first-order iecewise linear differential equation. Let {t,t 2,...,t k } denote the switch times when any variable of the network crosses 0. The solution of eq. (5) for each variable x i for t j <t<t j+ is x i (t) =x i (t j )e (t tj) +λ i (X i (j),x i2 (j),...,x ik (j))( e (t tj) ). (6) Deending on the articular network, eq. (5) can dislay steady states, limit cycles, quasieriodicity or chaos [9]-[3]. Although the origin of chaos in one articular fourdimensional network has been analyzed [3], no general methods have yet been develoed to determine whether chaotic dynamics exists in any given equation without integrating it. Chaos is the usual behavior in these networks for N = 64, 9 K 25, =[4]. In carrying out the comutations, we have made several additional assumtions concerning the structure of the equations. i) There is no self-inut or recirocal inut. This means that i cannot be an inut to itself, and if i is an inut to k, k is not also an inut to i. These assumtions eliminate the ossibility of asymtotic aroach to stable steady states in the neighborhood of threshold axes x i = 0 [9], []. If these assumtions are not made, the techniques for analysis of the dynamics based on the symbolic dynamics described below fail. ii) We set λ i in eq. (5) to be in the range ± or ±. iii) In the truth tables we assume that for a given value of, half the variables of the network have entries that are biased towards and half are biased towards 0. Although the iecewise-linear nature of eq. (5) facilitates the seed and accuracy of numerical integration over other integration techniques, numerical studies nevertheless grow
3 l. glass et al.: ordered and disordered dynamics etc. 60 (a) (b) (c) x(t) 0 x(t) 0 x(t) t t t Fig.. Time series for a tyical variable in a network with N = 50, K =8. (a)=, (b) =0.8, (c) =0.85. In anels (a) and (b), the time series satisfy our tests for chaos and in (c) the dynamics are eriodic. large raidly. Consequently, we focus on the order-disorder transition for limited regions of arameter sace. For each network, we investigated one initial condition and classified the dynamics after a transient of switch times. The dynamics can be classified by using symbolic dynamics, keeing track of the variable that switches as a function of time. Provided the network does not reach a steady state during the course of the integration, we generate a sequence of integers denoting the label of the variable that switches at each consecutive switch time. We search for eriodicities in this sequence of integers, searching for eriodicities u to length Given the restrictions on the equations mentioned above, if the symbolic sequence is eriodic, in the differential equation there is a stable limit cycle oscillation [9]. If the sequence is not eriodic, then in the differential equation, there can either be quasieriodic dynamics or chaotic dynamics. To identify quasieriodic dynamics, we need to kee track of the exact switching times of all the variables of the network. However, ilot investigations of several hundred randomly constructed networks failed to identify quasieriodic dynamics. Consequently, to simlify the analysis we lum quasieriodic rhythms together with chaotic networks, where we exect that the incidence of quasieriodic dynamics is negligible (less than %). Figure illustrates tyical dynamics for a single variable in a network with N = 50, K =8. The same connection matrix was assumed for three values, =,0.8, In fig. (a) and (b), the symbolic transition sequences corresonding to the time series do not show eriodicities. The Lyaunov number, which can be evaluated numerically using techniques in Ott [4], [5], is ositive for anels (a) and (b) and negative for anel (c). Based on the above, the dynamics are deterministic chaos in anels (a) and (b) and eriodic in anel (c). Figure 2 shows the numbers of networks dislaying steady states, limit cycles, and chaos for K =8andN= 50 anel (a), N = 00 anel (b), and N = 200 anel (c), for 50 different randomly generated networks for [,.0] incremented in stes of 2. As increases the number of chaotic networks decreases and the number of networks dislaying steady states increases. The values of associated with limit cycles are centered in an increasingly narrow range as N increases. Since the total number of orthants of hase sace that are visited during the simulations is a tiny fraction of the 2 N orthants of hase sace, we cannot exclude the ossibility that networks classified as chaotic will eventually reach a stable steady state or limit cycle if the integration times are extended. We now characterize the dynamics in eq. (5). Let j designate the time interval (t j,t j+ ), in eq. (6). At any given time in the interval j, the state in the continuous equation is maed to the logical state X(j). The distance between two logical states is the number of variables in which the activities differ. We call the distance between X(j) andx(j+ ), determined
4 602 EUROPHYSICS LETTERS from eq. (2), the outflow dimension, h(j), of X(j). The outflow dimension is a measure of the number of different variables that have the otential to cross 0 at the next switch time of the network [9]. If h(j) = 0, then the system will aroach a steady state. Using methods similar to those develoed in earlier work by Derrida and others [2], [5], [8], we comute the mean value of h. Denote the number of variables in logical state α, and with truth table entry β in time interval j by N αβ (j), where α, β {0,}.IfN 0 (j) =N 0 (j) =0, there is a fixed oint. If this is not the case there will eventually be a transition to a new state. The robability that there will be a transition of a variable from state 0 to state is P 0 = N 0 (j)/(n 0 (j)+n 0 (j)), and the robability that there will be a transition of a variable from state to state 0 is P 0 = N 0 (j)/(n 0 (j)+n 0 (j)). The exected number of inuts from the switching variable to variables in state (α, β) is ρn αβ,whereρ=k/n. Now assume that there is a transition of a variable from 0 at switch time t j+. We adot a robabilistic aroach to determine N αβ (j + ). Consider first the value of N 00 (j + ). This may be different from N 00 (j) if inuts from the variable that changed its state have inuts to variables in N 00 that leads to a change to N 0, or inuts to variables in N 0 that leads to changes to state N 00. Thus, we find that N 00 (j +)=N 00 (j) ρ( )N 00 (j)+ρn 0 (j). If there is a transition of a variable from 0attimet j, then the above exression would be changed to N 00 (j +)=N 00 (j)+ ρ( )N 00 (j)+ρn 0 (j). In similar fashion, by weighting the resective transition robabilities, we are led to the following system of equations for N αβ (j +): N 00 (j +)=N 00 (j)+p 0 ρ( )N 00 (j)+ρn 0 (j), N 0 (j +)=N 0 (j) P 0 + ρ( )N 00 (j) ρn 0 (j), N 0 (j +)=N 0 (j) P 0 ρ( )N 0 (j)+ρn (j), (7) N (j +)=N (j)+p 0 +ρ( )N 0 (j) ρn (j). At steady state, we have N αβ (j +)=N αβ (j) =Nαβ, where the asterisk reresents the steady-state value. Substituting this relation in eq. (7), we find that N0 = N 0. Substituting fraction of networks.0 (a) (b) (c) Fig. 2. Variation in the number of networks dislaying steady states (circles), limit cycles (squares), and chaos (diamonds) for K = 8, as a function of, (a)n= 50, (b) N = 00, (c) N = 200. A single initial condition was selected for each of 50 different networks for each value of.
5 l. glass et al.: ordered and disordered dynamics etc. 603 h*/n (a) fraction of networks.0 (b) 2 3 h* Fig. 3. (a) Normalized mean outflow dimension, h /N, as a function of for N = 50 (triangles), N = 00 (asterisks), N = 200 (lus signs) comared with the theoretical result from eq. (9). For each network, the value of h was averaged over 3000 switch times following a transient of switch times. (b) Fraction of chaotic networks as a function of h for N = 50 (triangles), N = 00 (asterisks), N = 200 (lus signs). this result in eq. (7), we find 2 = ρ( )N 00 + ρn0, 2 = ρ( )N 0 ρn. (8) Using the conservation condition, N00 +2N 0 +N = N, we have three simultaneous equations. Solving these equations and recalling that at steady state the mean outflow dimension, h,is given by N0 + N 0 =2N 0, we comute h N = K ( +2K 2K2 ). (9) In fig. 3(a), we comare the theoretical estimates from eq. (9) (solid curve) with the numerical comutations for N = 50, 00, 200. The theoretical estimate lies consistently above the numerically comuted values, but shows a similar functional deendence on. Whenh =0, we recover eq. (). In fig. 3(b), we lot the fraction of chaotic networks as a function of h. The transition occurs aroximately in the range 0 <h <20 for all values of N considered. The results in fig. 3 rovide a challenge for further theoretical analysis. These results have connections with the extensive studies carried out on the discrete time and discrete state sace switching networks [2]-[8]. In contrast to the earlier work, in which all finite networks must eventually cycle in the limit t and which do not therefore admit deterministic chaos, in eq. (5) deterministic chaos is ossible [3], [4]. At the moment, there are no general techniques to assert deterministic chaos in any given network, and it is ossible that eventually networks identified as chaotic here will reach limit cycles or steady states. Nevertheless, the extremely long transient behavior would aear to render these architectures imrobable for the highly constrained dynamics in real biological systems. In the continuous equations, the critical line () defines the line at which h = 0 and almost all networks with c dislay steady states. Thus, the continuous equations show transitions in dynamics even for finite N. Asthevalueofis varied from, there is a transition from chaotic dynamics to steady states, with an intervening zone of eriodic dynamics that becomes increasingly narrow as N increases.
6 604 EUROPHYSICS LETTERS *** This research has been suorted by grants from NSERC and FCAR. We thank T. Mestl, B. Sawhill, and S. Kauffman for helful conversations. The authors thank the Santa Fe Institute for hositality and for the use of comutational facilities. REFERENCES [] Kauffman S. A., J. Theor. Biol., 22 (969) 437; Kauffman S. A., Origins of Order: Self- Organization and Selection in Evolution (Oxford University Press, Oxford) 993. [2] Derrida B. and Pomeau Y., Eurohys. Lett.,, (986) 45; Derrida B. and Pomeau Y., Eurohys. Lett., 2 (986) 739. [3] Derrida B., Philos. Mag., 56 (987) 97. [4] Flyvbjerg H., J. Phys. A, 2 (988) L955. [5] Weisbuch G., Comlex Systems Dynamics (Addison-Wesley, Redwood City) 99. [6] Bastolla U. and Parisi G., Physica D, 98 (996) ; J. Theor. Biol., 87 (997) 7. [7] Bhattacharjya A. and Liang S., Phys. Rev. Lett., 77 (996) 644; Physica D, 95 (997) 29. [8] Luque B. and Solé R.V.,Phys. Rev. E, 55 (997) 257. [9] Glass L., J. Chem. Phys., 63 (975) 325; Glass L. and Pasternack J., J. Math. Biol., 6 (978) 207. [0] Thomas R. and Àri R. D., Biological Feedback (CRC Press, Boca Raton) 990. [] Mestl T., Plahte E. and Omholt S. W., J. Theor. Biol., 76 (995) 29; Mestl T., Plahte E. and Omholt S. W., Dyn. Stab. Systems, 2 (995) 79. [2] Bagley R. J. and Glass L., J. Theor. Biol., 83 (996) 269. [3] Mestl T., Lemay C. and Glass L., Physica D, 98 (996) 35. [4] Mestl T., Bagley R. J. and Glass L., Phys. Rev. Lett., 79 (997) 653. [5] The Lyaunov number can be exlicitly evaluated numerically using techniques in Ott E., Chaos in Dynamical Systems (Cambridge) 993, Chat. 4.
PHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationRobust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System
Vol.7, No.7 (4),.37-38 htt://dx.doi.org/.457/ica.4.7.7.3 Robust Predictive Control of Inut Constraints and Interference Suression for Semi-Trailer System Zhao, Yang Electronic and Information Technology
More informationarxiv: v2 [cond-mat.stat-mech] 10 Feb 2012
Ergodicity Breaking and Parametric Resonances in Systems with Long-Range Interactions Fernanda P. da C. Benetti, Tarcísio N. Teles, Renato Pakter, and Yan Levin Instituto de Física, Universidade Federal
More informationLyapunov exponents in random Boolean networks
Physica A 284 (2000) 33 45 www.elsevier.com/locate/physa Lyapunov exponents in random Boolean networks Bartolo Luque a;, Ricard V. Sole b;c a Centro de Astrobiolog a (CAB), Ciencias del Espacio, INTA,
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationarxiv:cond-mat/ v2 25 Sep 2002
Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,
More informationPeriod-two cycles in a feedforward layered neural network model with symmetric sequence processing
PHYSICAL REVIEW E 75, 4197 27 Period-two cycles in a feedforward layered neural network model with symmetric sequence rocessing F. L. Metz and W. K. Theumann Instituto de Física, Universidade Federal do
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationOutline. Markov Chains and Markov Models. Outline. Markov Chains. Markov Chains Definitions Huizhen Yu
and Markov Models Huizhen Yu janey.yu@cs.helsinki.fi Det. Comuter Science, Univ. of Helsinki Some Proerties of Probabilistic Models, Sring, 200 Huizhen Yu (U.H.) and Markov Models Jan. 2 / 32 Huizhen Yu
More informationControlling chaos in random Boolean networks
EUROPHYSICS LETTERS 20 March 1997 Europhys. Lett., 37 (9), pp. 597-602 (1997) Controlling chaos in random Boolean networks B. Luque and R. V. Solé Complex Systems Research Group, Departament de Fisica
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationProof: We follow thearoach develoed in [4]. We adot a useful but non-intuitive notion of time; a bin with z balls at time t receives its next ball at
A Scaling Result for Exlosive Processes M. Mitzenmacher Λ J. Sencer We consider the following balls and bins model, as described in [, 4]. Balls are sequentially thrown into bins so that the robability
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationBeam-Beam Stability in Electron-Positron Storage Rings
Beam-Beam Stability in Electron-Positron Storage Rings Bjoern S. Schmekel, Joseh T. Rogers Cornell University, Deartment of Physics, Ithaca, New York 4853, USA Abstract At the interaction oint of a storage
More informationChaotic Dynamics in an Electronic Model of a Genetic Network
Journal of Statistical Physics, Vol. 121, Nos. 5/6, December 2005 ( 2005) DOI: 10.1007/s10955-005-7009-y Chaotic Dynamics in an Electronic Model of a Genetic Network Leon Glass, 1 Theodore J. Perkins,
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More informationA Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split
A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Jul 1998
arxiv:cond-mat/98773v1 [cond-mat.stat-mech] 5 Jul 1998 Floy modes and the free energy: Rigidity and connectivity ercolation on Bethe Lattices P.M. Duxbury, D.J. Jacobs, M.F. Thore Deartment of Physics
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationCalculation of MTTF values with Markov Models for Safety Instrumented Systems
7th WEA International Conference on APPLIE COMPUTE CIENCE, Venice, Italy, November -3, 7 3 Calculation of MTTF values with Markov Models for afety Instrumented ystems BÖCÖK J., UGLJEA E., MACHMU. University
More informationA Simple Weight Decay Can Improve. Abstract. It has been observed in numerical simulations that a weight decay can improve
In Advances in Neural Information Processing Systems 4, J.E. Moody, S.J. Hanson and R.P. Limann, eds. Morgan Kaumann Publishers, San Mateo CA, 1995,. 950{957. A Simle Weight Decay Can Imrove Generalization
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationCharacterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations
Characterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations PINAR KORKMAZ, BILGE E. S. AKGUL and KRISHNA V. PALEM Georgia Institute of
More informationThe non-stochastic multi-armed bandit problem
Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at
More informationTemperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures
Indian Journal of Pure & Alied Physics Vol. 44, December 2006,. 953-958 Temerature, current and doing deendence of non-ideality factor for n and nn unch-through structures Khurshed Ahmad Shah & S S Islam
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationImproved Capacity Bounds for the Binary Energy Harvesting Channel
Imroved Caacity Bounds for the Binary Energy Harvesting Channel Kaya Tutuncuoglu 1, Omur Ozel 2, Aylin Yener 1, and Sennur Ulukus 2 1 Deartment of Electrical Engineering, The Pennsylvania State University,
More informationDetermining Momentum and Energy Corrections for g1c Using Kinematic Fitting
CLAS-NOTE 4-17 Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting Mike Williams, Doug Alegate and Curtis A. Meyer Carnegie Mellon University June 7, 24 Abstract We have used the
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationScaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering
More informationModeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation
6.3 Modeling and Estimation of Full-Chi Leaage Current Considering Within-Die Correlation Khaled R. eloue, Navid Azizi, Farid N. Najm Deartment of ECE, University of Toronto,Toronto, Ontario, Canada {haled,nazizi,najm}@eecg.utoronto.ca
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationUniversal Finite Memory Coding of Binary Sequences
Deartment of Electrical Engineering Systems Universal Finite Memory Coding of Binary Sequences Thesis submitted towards the degree of Master of Science in Electrical and Electronic Engineering in Tel-Aviv
More informationVariable Selection and Model Building
LINEAR REGRESSION ANALYSIS MODULE XIII Lecture - 38 Variable Selection and Model Building Dr. Shalabh Deartment of Mathematics and Statistics Indian Institute of Technology Kanur Evaluation of subset regression
More informationOnline Appendix to Accompany AComparisonof Traditional and Open-Access Appointment Scheduling Policies
Online Aendix to Accomany AComarisonof Traditional and Oen-Access Aointment Scheduling Policies Lawrence W. Robinson Johnson Graduate School of Management Cornell University Ithaca, NY 14853-6201 lwr2@cornell.edu
More informationwhere x i is the ith coordinate of x R N. 1. Show that the following upper bound holds for the growth function of H:
Mehryar Mohri Foundations of Machine Learning Courant Institute of Mathematical Sciences Homework assignment 2 October 25, 2017 Due: November 08, 2017 A. Growth function Growth function of stum functions.
More informationAn Adaptive Three-bus Power System Equivalent for Estimating Voltage Stability Margin from Synchronized Phasor Measurements
An Adative Three-bus Power System Equivalent for Estimating oltage Stability argin from Synchronized Phasor easurements Fengkai Hu, Kai Sun University of Tennessee Knoxville, TN, USA fengkaihu@utk.edu
More informationParticipation Factors. However, it does not give the influence of each state on the mode.
Particiation Factors he mode shae, as indicated by the right eigenvector, gives the relative hase of each state in a articular mode. However, it does not give the influence of each state on the mode. We
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationAI*IA 2003 Fusion of Multiple Pattern Classifiers PART III
AI*IA 23 Fusion of Multile Pattern Classifiers PART III AI*IA 23 Tutorial on Fusion of Multile Pattern Classifiers by F. Roli 49 Methods for fusing multile classifiers Methods for fusing multile classifiers
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
More informationA Note on the Positive Nonoscillatory Solutions of the Difference Equation
Int. Journal of Math. Analysis, Vol. 4, 1, no. 36, 1787-1798 A Note on the Positive Nonoscillatory Solutions of the Difference Equation x n+1 = α c ix n i + x n k c ix n i ) Vu Van Khuong 1 and Mai Nam
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationOn the capacity of the general trapdoor channel with feedback
On the caacity of the general tradoor channel with feedback Jui Wu and Achilleas Anastasooulos Electrical Engineering and Comuter Science Deartment University of Michigan Ann Arbor, MI, 48109-1 email:
More informationPretest (Optional) Use as an additional pacing tool to guide instruction. August 21
Trimester 1 Pretest (Otional) Use as an additional acing tool to guide instruction. August 21 Beyond the Basic Facts In Trimester 1, Grade 8 focus on multilication. Daily Unit 1: Rational vs. Irrational
More informationRe-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences
Re-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences J.A. Vallejos & S.M. McKinnon Queen s University, Kingston, ON, Canada ABSTRACT: Re-entry rotocols are
More informationComputer arithmetic. Intensive Computation. Annalisa Massini 2017/2018
Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More informationOil Temperature Control System PID Controller Algorithm Analysis Research on Sliding Gear Reducer
Key Engineering Materials Online: 2014-08-11 SSN: 1662-9795, Vol. 621, 357-364 doi:10.4028/www.scientific.net/kem.621.357 2014 rans ech Publications, Switzerland Oil emerature Control System PD Controller
More information+++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic +++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic
Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic Daniel Leitner, Johannes Krof, Günther Zauner, TU Vienna, Austria, dleitner@osiris.tuwien.ac.at Yuri Karov, Yuri Senichenkov, Yuri
More information%(*)= E A i* eiujt > (!) 3=~N/2
CHAPTER 58 Estimating Incident and Reflected Wave Fields Using an Arbitrary Number of Wave Gauges J.A. Zelt* A.M. ASCE and James E. Skjelbreia t A.M. ASCE 1 Abstract A method based on linear wave theory
More informationCorrespondence Between Fractal-Wavelet. Transforms and Iterated Function Systems. With Grey Level Maps. F. Mendivil and E.R.
1 Corresondence Between Fractal-Wavelet Transforms and Iterated Function Systems With Grey Level Mas F. Mendivil and E.R. Vrscay Deartment of Alied Mathematics Faculty of Mathematics University of Waterloo
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationAge of Information: Whittle Index for Scheduling Stochastic Arrivals
Age of Information: Whittle Index for Scheduling Stochastic Arrivals Yu-Pin Hsu Deartment of Communication Engineering National Taiei University yuinhsu@mail.ntu.edu.tw arxiv:80.03422v2 [math.oc] 7 Ar
More informationWave Drift Force in a Two-Layer Fluid of Finite Depth
Wave Drift Force in a Two-Layer Fluid of Finite Deth Masashi Kashiwagi Research Institute for Alied Mechanics, Kyushu University, Jaan Abstract Based on the momentum and energy conservation rinciles, a
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationResearch Article An iterative Algorithm for Hemicontractive Mappings in Banach Spaces
Abstract and Alied Analysis Volume 2012, Article ID 264103, 11 ages doi:10.1155/2012/264103 Research Article An iterative Algorithm for Hemicontractive Maings in Banach Saces Youli Yu, 1 Zhitao Wu, 2 and
More informationA Model for Randomly Correlated Deposition
A Model for Randomly Correlated Deosition B. Karadjov and A. Proykova Faculty of Physics, University of Sofia, 5 J. Bourchier Blvd. Sofia-116, Bulgaria ana@hys.uni-sofia.bg Abstract: A simle, discrete,
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationAbstract. We perform a perturbative QCD analysis of the quark transverse momentum
BIHEP-TH-96-09 The erturbative ion-hoton transition form factors with transverse momentum corrections Fu-Guang Cao, Tao Huang, and Bo-Qiang Ma CCAST (World Laboratory), P.O.Box 8730, Beijing 100080, China
More informationTHE 3-DOF helicopter system is a benchmark laboratory
Vol:8, No:8, 14 LQR Based PID Controller Design for 3-DOF Helicoter System Santosh Kr. Choudhary International Science Index, Electrical and Information Engineering Vol:8, No:8, 14 waset.org/publication/9999411
More informationOn Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada
On Using FASTEM2 for the Secial Sensor Microwave Imager (SSM/I) March 15, 2001 Godelieve Deblonde Meteorological Service of Canada 1 1. Introduction Fastem2 is a fast model (multile-linear regression model)
More informationAvailability and Maintainability. Piero Baraldi
Availability and Maintainability 1 Introduction: reliability and availability System tyes Non maintained systems: they cannot be reaired after a failure (a telecommunication satellite, a F1 engine, a vessel
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationarxiv: v1 [math-ph] 10 Sep 2010
Calculation of generalized ubbell rectangular source integral Jonathan Murley and Nasser Saad arxiv:9.7v [math-h] Se Deartment of Mathematics and Statistics, University of Prince Edward Island Charlottetown,
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
More informationAlgorithms for Air Traffic Flow Management under Stochastic Environments
Algorithms for Air Traffic Flow Management under Stochastic Environments Arnab Nilim and Laurent El Ghaoui Abstract A major ortion of the delay in the Air Traffic Management Systems (ATMS) in US arises
More informationFactors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier
Australian Journal of Basic and Alied Sciences, 5(12): 2010-2020, 2011 ISSN 1991-8178 Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doed Fiber Amlifier
More informationarxiv: v1 [cond-mat.stat-mech] 1 Dec 2017
Suscetibility Proagation by Using Diagonal Consistency Muneki Yasuda Graduate School of Science Engineering, Yamagata University Kazuyuki Tanaka Graduate School of Information Sciences, Tohoku University
More informationMaximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming
Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August
More informationLearning Sequence Motif Models Using Gibbs Sampling
Learning Sequence Motif Models Using Gibbs Samling BMI/CS 776 www.biostat.wisc.edu/bmi776/ Sring 2018 Anthony Gitter gitter@biostat.wisc.edu These slides excluding third-arty material are licensed under
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationLecture 6. 2 Recurrence/transience, harmonic functions and martingales
Lecture 6 Classification of states We have shown that all states of an irreducible countable state Markov chain must of the same tye. This gives rise to the following classification. Definition. [Classification
More informationJohn Weatherwax. Analysis of Parallel Depth First Search Algorithms
Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel
More informationModel checking, verification of CTL. One must verify or expel... doubts, and convert them into the certainty of YES [Thomas Carlyle]
Chater 5 Model checking, verification of CTL One must verify or exel... doubts, and convert them into the certainty of YES or NO. [Thomas Carlyle] 5. The verification setting Page 66 We introduce linear
More informationVISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES
VISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES V. V. Novikov ), K.W. Wojciechowski ) ) Odessa National Polytechnical University, Shevchenko Prosekt, 6544 Odessa, Ukraine; e-mail: novikov@te.net.ua
More informationA PEAK FACTOR FOR PREDICTING NON-GAUSSIAN PEAK RESULTANT RESPONSE OF WIND-EXCITED TALL BUILDINGS
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-1, 009, Taiei, Taiwan A PEAK FACTOR FOR PREDICTING NON-GAUSSIAN PEAK RESULTANT RESPONSE OF WIND-EXCITED TALL BUILDINGS M.F. Huang 1,
More informationEstimation of Separable Representations in Psychophysical Experiments
Estimation of Searable Reresentations in Psychohysical Exeriments Michele Bernasconi (mbernasconi@eco.uninsubria.it) Christine Choirat (cchoirat@eco.uninsubria.it) Raffaello Seri (rseri@eco.uninsubria.it)
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.
More informationDETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS
Proceedings of DETC 03 ASME 003 Design Engineering Technical Conferences and Comuters and Information in Engineering Conference Chicago, Illinois USA, Setember -6, 003 DETC003/DAC-48760 AN EFFICIENT ALGORITHM
More informationDomain Dynamics in a Ferroelastic Epilayer on a Paraelastic Substrate
Y. F. Gao Z. Suo Mechanical and Aerosace Engineering Deartment and Princeton Materials Institute, Princeton University, Princeton, NJ 08544 Domain Dynamics in a Ferroelastic Eilayer on a Paraelastic Substrate
More informationOne step ahead prediction using Fuzzy Boolean Neural Networks 1
One ste ahead rediction using Fuzzy Boolean eural etworks 1 José A. B. Tomé IESC-ID, IST Rua Alves Redol, 9 1000 Lisboa jose.tome@inesc-id.t João Paulo Carvalho IESC-ID, IST Rua Alves Redol, 9 1000 Lisboa
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationAnalysis of execution time for parallel algorithm to dertmine if it is worth the effort to code and debug in parallel
Performance Analysis Introduction Analysis of execution time for arallel algorithm to dertmine if it is worth the effort to code and debug in arallel Understanding barriers to high erformance and redict
More informationIntroduction to MVC. least common denominator of all non-identical-zero minors of all order of G(s). Example: The minor of order 2: 1 2 ( s 1)
Introduction to MVC Definition---Proerness and strictly roerness A system G(s) is roer if all its elements { gij ( s)} are roer, and strictly roer if all its elements are strictly roer. Definition---Causal
More informationDistributed K-means over Compressed Binary Data
1 Distributed K-means over Comressed Binary Data Elsa DUPRAZ Telecom Bretagne; UMR CNRS 6285 Lab-STICC, Brest, France arxiv:1701.03403v1 [cs.it] 12 Jan 2017 Abstract We consider a networ of binary-valued
More informationMeshless Methods for Scientific Computing Final Project
Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking
More informationarxiv: v1 [physics.data-an] 20 Sep 2011
Strategy of Cometition between Two Grous based on a Contrarian Oinion Model Qian Li 1, Lidia A. Braunstein, 2,1, Shlomo Havlin, 3, and H. Eugene Stanley 1 1 Deartment of Physics and Center for Polymer
More information