Collective Inference on Markov Models for Modeling Bird Migration
|
|
- Evan Hancock
- 5 years ago
- Views:
Transcription
1 Collective Inference on Markov Models for Modeling Bird Migration Daniel Sheldon Cornell University M. A. Saleh Elmohamed Cornell University Dexter Kozen Cornell University Abstract We investigate a family of inference problems on Markov models, where many sample paths are drawn from a Markov chain and partial information is revealed to an observer who attempts to reconstrct the sample paths. We present algorithms and hardness reslts for several variants of this problem which arise by revealing different information to the observer and imposing different reqirements for the reconstrction of sample paths. Or algorithms are analogos to the classical Viterbi algorithm for Hidden Markov Models, which finds the single most probable sample path given a seqence of observations. Or work is motivated by an important application in ecology: inferring bird migration paths from a large database of observations. 1 Introdction Hidden Markov Models (HMMs) assme a generative model for seqential data whereby a seqence of states (or sample path) is drawn from a Markov chain in a hidden experiment. Each state generates an otpt symbol from alphabet Σ, and these otpt symbols constitte the data or observations. A classical problem, solved by the Viterbi algorithm, is to find the most probable sample path given the observations for a given Markov model. We call this the single path problem; it is well sited to labeling or tagging a single seqence of data. For example, HMMs have been sccessflly applied in speech recognition [1], natral langage processing [2], and biological seqencing [3]. We introdce two generalizations of the single path problem for performing collective inference on Markov models, motivated by an effort to model bird migration patterns sing a large database of static observations. The ebird database hosted by the Cornell Lab of Ornithology contains millions of bird observations from throghot North America, reported by the general pblic sing the ebird web application. 1 Observations report location, date, species and nmber of birds observed. The ebird data set is very rich; the hman eye can easily discern migration patterns from animations showing the observations as they nfold over time on a map of North America. 2 However, the ebird data are static, and they do not explicitly record movement, only the distribtions at different points in time. Conclsions abot migration patterns are made by the hman observer. Or goal is to bild a mathematical framework to infer dynamic migration models from the static ebird data. Qantitative migration models are of great scientific and practical import: for example, this problem arose ot of an interdisciplinary project at Cornell University to model the possible spread of avian inflenza in North America throgh wild bird migration. The migratory behavior for a species of birds can be modeled as a single generative process that independently governs how individal birds fly between locations, giving rise to the following inference problem: a hidden experiment simltaneosly draws many independent sample paths from
2 a Markov chain, and the observations reveal aggregate information abot the collection of sample paths at each time step. For example, the ebird data estimate the geographical distribtion of a species on sccessive days, bt do not track individal birds. We discss two problems within this framework. In the mltiple path problem, we assme that M independent sample paths are simltaneosly drawn from a Markov model, and the observations reveal the nmber of paths that otpt symbol α at time t, for each α and t. The observer seeks the most likely collection of paths given the observations. The fractional path problem is a frther generalization in which paths are divisible entities. The observations reveal the fraction of paths that otpt symbol α at time t, and the observer s job is to find the most likely (in a sense to be defined later) weighted collection of paths given the observations. Conceptally, the fractional path problem can be derived from the mltiple path problem by letting M go to infinity; or it has a probabilistic interpretation in terms of distribtions over paths. After discssing some preliminaries in section 2, sections 3 and 4 present algorithms for the mltiple and fractional path problems, respectively, sing network flow techniqes on the trellis graph of the Markov model. The mltiple path problem in its most general form is NP-hard, bt can be solved as an integer program. The special case when otpt symbols niqely identify their associated states can be solved efficiently as a flow problem; althogh the single path problem is trivial in this case, the mltiple and fractional path problems remain interesting. The fractional path problem can be solved by linear programming. We also introdce a practical extension to the fractional paths problem, inclding slack variables allowing the soltion to deviate slightly from potentially noisy observations. In section 5, we demonstrate or techniqes with visalizations for the migration of Archilochs colbris, the Rby-throated Hmmingbird, devoting some attention to a challenging problem we have neglected so far: estimating species distribtions from ebird observations. We briefly mention some related work. Carana et al. [4] and Phillips et al. [5] sed machine learning techniqes to model bird distribtions from observations and environmental featres. For problems on seqential data, many variants of HMMs have been proposed [3], and recently, conditional random fields (CRFs) have become a poplar alternative [6]. Roth and Yih [7] present an integer programming inference framework for CRFs that is similar to or problem formlations. 2 Preliminaries 2.1 Data Model and Notation A Markov model (V, p, Σ, σ) is a Markov chain with state set V and transition probabilities p(, v) for all, v V. Each state generates a niqe otpt symbol from alphabet Σ, given by the mapping σ : V Σ. Althogh some presentations allow each state to otpt mltiple symbols with different emission probabilities, we lose no generality assming that each state emits a niqe symbol to encode a model where state v otpt mltiple symbols, we simply dplicate v for each symbol and encode the emission probabilities into the transitions. Of corse, σ need not be one-to-one. It is sefl to think of σ as a partition of the states, letting V α = σ 1 (α) be the set of all states that otpt α. We assme each model has a distingished start state s and otpt symbol start. Let Y = V T be the set of all possible sample paths of length T. We represent a path y Y as a row vector y = (y 1,..., y T ), and a collection of M paths as the M T matrix Y = (y it ), with each row y i representing an independent sample path. The transition probabilities indce a distribtion λ on Y, where λ(y) = T 1 t=1 p(y t, y t+1 ). We will also consider arbitrary distribtions π over Y, letting Y = (Y 1,..., Y T ) denote a random path from π. Then, for example, we write Pr π [Y t = ] to be the probability nder π that the tth state is, and E π [f(y )] to be the expected vale of f(y ) for any fnction f of a random path Y drawn from π. Note that Y (boldface) denotes a matrix of M paths, while Y denotes a random path. 2.2 The Trellis Graph and Viterbi as Shortest Path To develop or flow-based algorithms, it is instrctive to bild pon a shortest-path interpretation of the Viterbi algorithm [7]. In an instance of the single path problem we are given a model (V, p, Σ, σ) and observations α 1,..., α T, and we seek the most probable path y given the observations. We call path y feasible if σ(y t ) = α t for all t; then we wish to maximize λ(y) over feasible y. The problem 2
3 p(, ) c(, ) s p(s, ) v V 0 V 0 V 0 p(, w) s c(s, ) v V 0 V 0 V 0 c(, w) w V 1 V 1 V 1 w V 1 V 1 V 1 start 0 Observations 0 1 start 0 Observations 0 1 (a) (b) Figre 1: Trellis graph for Markov model with states {s,, v, w} and alphabet {start, 0, 1}. States and v otpt the symbol 0, and state w otpts the symbol 1. (a) The bold path is feasible for the specified observations, with probability p(s, )p(, )p(, w). (b) Infeasible edges have been removed (indicated by light dashed lines), and probabilities changed to costs. The bold path has cost c(s, ) + c(, ) + c(, w). is conveniently illstrated sing the trellis graph of the Markov model (Figre 1). Here, the states are replicated for each time step, and edges connect a state at time t to its possible sccessors at time t + 1, labeled with the transition probability. A feasible path mst pass throgh partition V αt at step t, so we can prne all edges incident on other partitions, leaving only feasible paths. By defining the cost of an edge as c(, v) = log p(, v), and letting the path cost c(y) be the sm of its edge costs, straightforward algebra shows that arg max y λ(y) = arg min y c(y), i.e., the path of maximm probability becomes the path of minimm cost nder this transformation. Ths the Viterbi algorithm finds the shortest feasible path in the trellis sing edge lengths c(, v). 3 Mltiple Path Problem In the mltiple path problem, M sample paths are drawn from the model and the observations reveal the nmber of paths N t (α) that otpt α at time t, for all α and t; or, eqivalently, the mltiset A t of otpt symbols at time t. The objective is to find the most probable collection Y that is feasible, meaning it prodces mltisets A 1,..., A T. The probability λ(y) is jst the prodct of the path-wise probabilities: M M λ(y) = λ(y i ) = i=1 Then the formal specification of this problem is T 1 i=1 t=1 p(y i,t, y i,t+1 ). (1) max Y λ(y) sbject to {i : y i,t V α } = N t (α) for all α, t. (2) 3.1 Redction to the Single Path Problem A naive approach to the mltiple path problem redces it to the single path problem by creating a new Markov model on state set V M where state v 1,..., v M encodes an entire tple of original states, and the transition probabilities are given by the prodct of the element-wise transition probabilities: p( 1,..., M, v 1,..., v M ) = M p( i, v i ). A state from the prodct space V M corresponds to an entire colmn of the matrix Y, and by changing the order of mltiplication in (1), we see that the probability of a path in the new model is eqal to the probability of the entire collection of paths in the old model. To complete the redction, we form a new alphabet ˆΣ whose symbols represent mltisets of size M on Σ. Then the soltion to (2) can be fond by rnning the Viterbi algorithm to find the most likely seqence of states from V M that prodce otpt symbols (mltisets) A 1,..., A T. The rnning time is polynomial in V M and ˆΣ, bt exponential in M. i=1 3
4 3.2 Graph Flow Formlation Can we do better than the naive approach? Viewing the cost of a path as the cost of roting one nit of flow along that path in the trellis, a minimm cost collection of M paths is eqivalent to a minimm cost flow of M nits throgh the trellis given M paths, we can rote one nit along each to get a flow, and we can decompose any flow of M nits into paths each carrying a single nit of flow. Ths we can write the optimization problem in (2) as the following flow integer program, with additional constraints that the flow paths generate the correct observations. The decision variable x t v indicates the flow traveling from to v at time t; or, the nmber of sample paths that transition from to v at time t. (IP) min,v,t c(, v)x t v s.t. x t v = x t+1 vw for all v, t, (3) w x t v = N t (α) for all α, t, (4) V α,v V x t v N for all, v, t. The flow conservation constraints (3) are standard: the flow into v at time t is eqal to the flow leaving v at time t + 1. The observation constraints (4) specify that N t (α) nits of flow leave partition V α at time t. These also imply that exactly M nits of flow pass throgh each level of the trellis, by smming over all α, x t v = x t v = N t (α) = M.,v α V α,v V α Withot the observation constraints, IP wold be an instance of the minimm-cost flow problem [8], which is solvable in polynomial time by a variety of algorithms [9]. However, we cannot hope to encode the observation constraints into the flow framework, de to the following reslt. Lemma 1. The mltiple path problem is NP-hard. The proof of Lemma 1 is by redction from SET COVER, and is omitted. One may se a general prpose integer program solver to solve IP directly; this may be efficient in some cases despite the lack of polynomial time performance garantees. In the following sections we discss alternatives that are efficiently solvable. 3.3 An Efficient Special Case In the special case when σ is one-to-one, the otpt symbols niqely identify their generating states, so we may assme that Σ = V, and the otpt symbol is always the name of the crrent state. To see how the problem IP simplifies, we now have V = {} for all, so each partition consists of a single state, and the observations completely specify the flow throgh each node in the trellis: x t v = N t () for all, t. (4 ) v Sbstitting the new observation constraints (4 ) for time t+1 into the RHS of the flow conservation constraints (3) for time t yield the following replacements: x t v = N t+1 (v) for all v, t. (3 ) This gives an eqivalent set of constraints, each of which refers only to variables x t v for a single t. Hence the problem can be decomposed into T 1 disjoint sbproblems for t = 1,..., T 1. The tth sbproblem IP t is given in Figre 2(a), and illstrated on the trellis in Figre 2(b). State at time t has a spply of N t () nits of flow coming from the previos step, and we mst rote N t+1 (v) nits of flow to state v at time t+1, so we place a demand of N t+1 (v) at the corresponding node. Then the problem redces to finding a minimm cost roting of the spply from time t to meet the demand at time t + 1, solved separately for all t = 1,..., T 1. The problem IP t an instance of the transportation problem [10], a special case of the minimm-cost flow problem. There are a variety of efficient algorithms to solve both problems [8,9], or one may se a general prpose linear program (LP) solver; any basic soltion to the LP relaxation of IP t is garanteed to be integral [8]. 4
5 (IP t ) min c(, v)x t v,v s.t. x t v = N t+1 (v) v, (3 ) x t v = N t (), (4 ) v x t v N, v. (a) Nt( ) Spply v 1 c(v1, v1) v 2 v 3 v 1 v 2 v 3 t t + 1 (b) Nt+1( ) Demand Figre 2: (a) The definition of sbproblem IP t. (b) Illstration on the trellis. 4 Fractional Paths Problem In the fractional paths problem, a path is a divisible entity. The observations specify q t (α), the fraction of paths that otpt α at time t, and the observer chooses π(y) fractional nits of each path y, totaling one nit, sch that q t (α) nits otpt α at time t. The objective is to maximize y Y λ(y)π(y). Pt another way, π is a distribtion over paths sch that Pr π [Y t V α ] = q t (α), i.e., q t specifies the marginal distribtion over symbols at time t. By taking the logarithm, an eqivalent objective is to maximize E π [log λ(y )], so we seek the distribtion π that maximizes the expected log-probability of a path Y drawn from π. Conceptally, the fractional paths problem arises by letting M in the mltiple paths problem and normalizing to let q t (α) = N t (α)/m specify the fraction of paths that otpt α at time t. Operationally, the fractional paths problem is modeled by the LP relaxation of IP, which rotes one splittable nit of flow throgh the trellis. (RELAX) min,v,t c(, v)x t v s.t. x t v = x t+1 vw for all v, t, w x t v = q t (α) for all α, t, (5) V α v V x t v 0 for all, v, t. It is easy to see that a single-nit flow x corresponds to a probability distribtion π. Given π, let x t v = Pr π [Y t =, Y t+1 = v]; then x is a flow becase the probability a path enters v at time t is eqal to the probability it leaves v at time t + 1. Conversely, given x, any path decomposition assigning flow π(y) to each y Y is a probability distribtion becase the total flow is one. In general, the decomposition is not niqe, bt any choice yields a distribtion π with the same objective vale. Frthermore, nder this correspondence, x satisfies the marginal constraints (5) if and only if π has the correct marginals: x t v = Pr [Y t =, Y t+1 = v] = Pr [Y t = ] = Pr [Y t V α ]. v V v V V α V α V α Finally, we can rewrite the objective fnction in terms of paths: π(y)c(y) = E π [c(y )] = E π [ log λ(y )].,v,t c(, v)x t v = y Y By switching signs and changing from minimization to maximization, we see that RELAX solves the fractional paths problem. This problem is very similar to maximm entropy or minimm cross entropy modeling, bt the details are slightly different: sch a model wold typically find the distribtion π with the correct marginals that minimizes the cross entropy or Kllback-Leibler divergence [11] between λ and π, which, after removing a constant term, redces to minimizing E λ [ log π(y )]. Like IP, the RELAX problem also decomposes into sbproblems in the case when σ is one-to-one, bt this simplification is incompatible with the slack variables introdced in the following section. 5
6 4.1 Incorporating Slack In or application, the marginal distribtions q t ( ) are themselves estimates, and it is sefl to allow the LP to deviate slightly from these marginals to find a better overall soltion. To accomplish this, we add slack variables δ t into the marginal constraints (5), and charge for the slack in the objective fnction. The new marginals constraints are x t v = q t (α) + δα t for all α, t, (5 ) V α v V and we add the term α,t γt α δ t α into the objective fnction to charge for the slack, sing a standard LP trick [8] to model the absolte vale term. The slack costs γ t α can be tailored to individal inpt vales; for example, yo may want to charge more to deviate from a confident estimate. This will depend on the specific application. We also add the necessary constraints to ensre that the new marginals q t(α) = q t (α) + δ t α form a valid probability distribtion for all t. 5 Demonstration In this section, we demonstrate or techniqes by sing the fractional paths problem to help visalize the migration of Archilochs colbris, the Rby-throated Hmmingbird, a common bird whose range is relatively well covered by ebird observations. We work in discretized space and time, dividing the map into grid cells and the year into weeks. In addition to the Markov model for transitions between locations (grid cells) in sccessive weeks, we reqire estimates q t ( ) for the weekly distribtions of hmmingbirds across locations. Since the actal ebird observations are highly non-niform in space and time, estimating weekly distribtions reqires significant inference for locations with few or no observations. Althogh or focs is on migration inference, this is an important part of the overall problem. In the appendix, we otline one approach based on on harmonic energy minimzation [12], bt or methods are compatible with any techniqe that prodces weekly distribtions q t () and slack costs γ t. Finally, althogh or final observations are distribtions over states (locations) and not otpt symbols i.e., σ is one-to-one we cannot se the simplification from section 3.3 becase we incorporate slack into the model. 5.1 ebird Data Lanched in 2002, ebird is a citizen science project rn by the Cornell Lab of Ornithology, leveraging the data gathering power of the pblic. On the ebird website, birdwatchers sbmit checklists of birds they observe, indicating a cont for each species, along with the location, date, time and additional information. Or data set consists of the 428,648 complete checklists from 1995 throgh 2007, meaning the reporter listed all species observed. This means we can infer a cont of zero, or a negative observation, for any species not listed. Using a USGS land cover map, we divide North America into grid cells that are roghly 225 km on a side. All years of data are aggregated into one, and the year is divided into weeks so t = 1,..., 52 represents the week of the year. 5.2 Migration Inference Given weekly distribtions q t () and slack costs γ t (see the appendix), it remains to specify the Markov model. We se a simple Gassian model favoring short flights, letting p(, v) exp( d(, v) 2 /σ 2 ), where d(, v) measres the distance between grid cell centers. This corresponds to a sqared distance cost fnction. To redce problem size, we omitted variables x t v from the LP when d(, v) > 1350 km, effectively setting p(, v) = 0. We also fond it sefl to impose pper bonds δ t 2q t () on the slack variables so no single vale cold increase by more than a factor of 3. Or final LP, which was solved sing the MOSEK optimization toolbox, had 78,521 constraints and 3,031,116 variables. Figre 3 displays the migration paths or model inferred for the for weeks indicated. The top row shows the distribtion and paths inferred by the model; grid cells colored in darker shades of red have more birds (higher vales for q t()). Ble arrows indicate flight paths (x t v) between the week shown and the following week, with line width in proportion to x t v. In the bottom row, the raw data is given for comparison. Yellow dots indicate negative observations; ble sqares indicate 6
7 Week 10 Week 20 Week 30 Week 40 3/5 5/14 7/28 10/01 Figre 3: Rby-throated Hmmingbird migration. See text for description. Observations... Sink t 1 t t+1 Figre 4: Range of Rby-throated Hmmingbird. From BNA [15]. Figre 5: Illstration of local connectivity for lattice point t (appendix). positive observations, with size proportional to cont. Locations with both positive and negative observations appear dark gray. The inferred distribtions and paths are consitent with both seasonal ranges (Figre 4), and what experts crrently know abot migration rotes. For example, in the smmary paragraph on migration from the Archilochs colbris species accont in Birds of North America, Robinson et al. write Many fly across Glf of Mexico, bt many also follow coastal rote. Rotes may differ for north- and sothbond birds. Acknowledgements We are gratefl to Daniel Fink, Wesley Hochachka and Steve Kelling from the Cornell Lab of Ornithology for sefl discssions. This work was spported in part by ONR Grant N and by NSF grant CCF The views and conclsions herein are those of the athors and do not necessarily represent the official policies or endorsements of these organizations or the US Government. 7
8 References [1] L. R. Rabiner. A ttorial on hidden Markov models and selected applications in speech recognition. Proceedings of the IEEE, 77(2): , [2] E. Charniak. Statistical techniqes for natral langage parsing. AI Magazine, 18(4):33 44, [3] R. Drbin, S. Eddy, A. Krogh, and G. Mitchison. Biological seqence analysis: Probabilistic models of proteins and ncleic acids. Cambridge University Press, [4] R. Carana, M. Elhawary, A. Mnson, M. Riedewald, D. Sorokina, D. Fink, W. M. Hochachka, and S. Kelling. Mining citizen science data to predict prevalence of wild bird species. In SIGKDD, [5] S. J. Phillips, M. Ddík, and R. E. Schapire. A maximm entropy approach to species distribtion modeling. In ICML, [6] J. Lafferty, A. McCallm, and F. Pereira. Conditional random fields: Probabilistic models for segmenting and labeling seqence data. ICML, [7] D. Roth and W. Yih. Integer linear programming inference for conditional random fields. ICML, [8] V. Chvátal. Linear Programming. W.H. Freeman, New York, NY, [9] A. V. Goldberg, S. A. Plotkin, and E. Tardos. Combinatorial algorithms for the generalized circlation problem. Math. Oper. Res., 16(2): , [10] G. B. Dantzig. Application of the simplex method to a transportation problem. In T. C. Koopmans, editor, Activity Analysis of Prodction and Allocation, volme 13 of Cowles Commission for Research in Economics, pages Wiley, [11] J. Shore and R. Johnson. Properties of cross-entropy minimization. IEEE Trans. on Information Theory, 27: , [12] X. Zh, Z. Ghahramani, and J. Lafferty. Semi-spervised learning sing Gassian fields and harmonic fnctions. In ICML, [13] P. G. Doyle and J. L. Snell. Random walks and electric networks. abstract?id=oai:arxiv.org:math/ , [14] D. Aldos and J. Fill. Reversible Markov Chains and Random Walks on Graphs. Monograph in Preparation, [15] T. R. Robinson, R. R. Sargent, and M. B. Sargent. Rby-throated Hmmingbird (Archilochs colbris). In A. Poole and F. Gill, editors, The Birds of North America, nmber 204. The Academy of Natral Sciences, Philadelphia, and The American Ornithologists Union, Washington, D.C., A Estimating Weekly Distribtions from ebird Or goal is to estimate q t (), the fraction of birds in grid cell dring week t. Given enogh observations, we can estimate q t () sing the average nmber of birds conted per checklist, a qantity we call the rate r t (). However, even for a bird with good ebird coverage, there are cells with few or no observations dring some weeks. To fill these gaps, we se the harmonic energy minimization techniqe [12] to determine vales for empty cells based on cells that are nearby in space and time. To start, we bild a 3-dimensional lattice on the points t, where t represents cell dring week t. Edges are weighted, with weights representing similarity between points. The connections for point t are illstrated in Figre 5. We exclde the edge between cells t, v t when the line connecting the centers of and v is more than half water. All edges have weight 1, except for edges connecting t to t 1 and t+1, which have weight 1/4, to prioritize spatial similarity over temporal similarity. Given the lattice, harmonic energy minimzation learns a fnction f on the lattice; the idea is to match r t () on points with sfficient data and find vales for other points according to the similarity strctre. To this end, we designate some bondary points for which the vale of f is fixed, while other points are interior points. The vale of f at an interior point t is determined by the expected vale of the following random experiment. Perform a random walk starting from t, following otgoing edges with probability proportional to their weight. When a bondary point is reached, terminate the walk and accept the bondary vale. In this way, the vales at interior points are a weighted average of nearby bondary points, where nearness is interpreted as the absportion probability in an absorbing random walk. We derive a measre of confidence in the vale f( t ) from the same experiment: let h( t ) be the expected nmber of steps for the random walk from t 8
9 to hit the bondary (the hitting time [14] of the bondary set). When h( t ) is small, t is close to the bondary and we are more confident in f( t ). Rather than choosing a threshold on the nmber of observations reqired to be a bondary point, we create a soft bondary by designating all lattice points t as interior points, and adding one bondary node connected to t for each observation, with vale eqal to the nmber of birds observed. As t gains more observations, the probability of reaching an observation in the first step approaches one, and f( t ) approaches r t (), the average of the observations. As a conservative measre, each node is also connected to a sink with bondary vale 0, to prevent vales from propagating over very long distances. We compte h and f iteratively [13]. Since f( t ) approximates the rate r t (), we mltiply by the land mass of cell to get a (relative) estimate ˆq t () for the nmber of birds in cell at time t. Finally, we normalize ˆq for each t, taking q t () = ˆq t ()/ ˆq t(). For slack costs, we set γ t = γ 0 /h( t ) to be inversely proportional to bondary hitting time, with γ chosen in conjnction with the transition costs in section 5.2 so the average cost for a nit of slack is the same as moving 40 pixels. 9
Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationClassify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.
Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports
More information1. Tractable and Intractable Computational Problems So far in the course we have seen many problems that have polynomial-time solutions; that is, on
. Tractable and Intractable Comptational Problems So far in the corse we have seen many problems that have polynomial-time soltions; that is, on a problem instance of size n, the rnning time T (n) = O(n
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationNetwork Coding for Multiple Unicasts: An Approach based on Linear Optimization
Network Coding for Mltiple Unicasts: An Approach based on Linear Optimization Danail Traskov, Niranjan Ratnakar, Desmond S. Ln, Ralf Koetter, and Mriel Médard Abstract In this paper we consider the application
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationStability of Model Predictive Control using Markov Chain Monte Carlo Optimisation
Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationAssignment Fall 2014
Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationChapter 4 Supervised learning:
Chapter 4 Spervised learning: Mltilayer Networks II Madaline Other Feedforward Networks Mltiple adalines of a sort as hidden nodes Weight change follows minimm distrbance principle Adaptive mlti-layer
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationModelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationA New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation
A ew Approach to Direct eqential imlation that Acconts for the Proportional ffect: Direct ognormal imlation John Manchk, Oy eangthong and Clayton Detsch Department of Civil & nvironmental ngineering University
More informationUniversal Scheme for Optimal Search and Stop
Universal Scheme for Optimal Search and Stop Sirin Nitinawarat Qalcomm Technologies, Inc. 5775 Morehose Drive San Diego, CA 92121, USA Email: sirin.nitinawarat@gmail.com Vengopal V. Veeravalli Coordinated
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
More informationSensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters
Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed
More informationThe Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks
The Cryptanalysis of a New Pblic-Key Cryptosystem based on Modlar Knapsacks Yeow Meng Chee Antoine Jox National Compter Systems DMI-GRECC Center for Information Technology 45 re d Ulm 73 Science Park Drive,
More informationTheoretical and Experimental Implementation of DC Motor Nonlinear Controllers
Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationBayes and Naïve Bayes Classifiers CS434
Bayes and Naïve Bayes Classifiers CS434 In this lectre 1. Review some basic probability concepts 2. Introdce a sefl probabilistic rle - Bayes rle 3. Introdce the learning algorithm based on Bayes rle (ths
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationIntrodction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc
A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper
More informationEssentials of optimal control theory in ECON 4140
Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as
More informationJoint Transfer of Energy and Information in a Two-hop Relay Channel
Joint Transfer of Energy and Information in a Two-hop Relay Channel Ali H. Abdollahi Bafghi, Mahtab Mirmohseni, and Mohammad Reza Aref Information Systems and Secrity Lab (ISSL Department of Electrical
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationLecture 8: September 26
10-704: Information Processing and Learning Fall 2016 Lectrer: Aarti Singh Lectre 8: September 26 Note: These notes are based on scribed notes from Spring15 offering of this corse. LaTeX template cortesy
More informationA Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model
entropy Article A Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model Hyenkyn Woo School of Liberal Arts, Korea University of Technology and Edcation, Cheonan
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationCHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy
CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK Wassim Joini and Christophe Moy SUPELEC, IETR, SCEE, Avene de la Bolaie, CS 47601, 5576 Cesson Sévigné, France. INSERM U96 - IFR140-
More informationReflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by
Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length
More informationWorst-case analysis of the LPT algorithm for single processor scheduling with time restrictions
OR Spectrm 06 38:53 540 DOI 0.007/s009-06-043-5 REGULAR ARTICLE Worst-case analysis of the LPT algorithm for single processor schedling with time restrictions Oliver ran Fan Chng Ron Graham Received: Janary
More informationStep-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter
WCE 007 Jly - 4 007 London UK Step-Size onds Analysis of the Generalized Mltidelay Adaptive Filter Jnghsi Lee and Hs Chang Hang Abstract In this paper we analyze the bonds of the fixed common step-size
More informationPrediction of Transmission Distortion for Wireless Video Communication: Analysis
Prediction of Transmission Distortion for Wireless Video Commnication: Analysis Zhifeng Chen and Dapeng W Department of Electrical and Compter Engineering, University of Florida, Gainesville, Florida 326
More informationModel Discrimination of Polynomial Systems via Stochastic Inputs
Model Discrimination of Polynomial Systems via Stochastic Inpts D. Georgiev and E. Klavins Abstract Systems biologists are often faced with competing models for a given experimental system. Unfortnately,
More informationMulti-Voltage Floorplan Design with Optimal Voltage Assignment
Mlti-Voltage Floorplan Design with Optimal Voltage Assignment ABSTRACT Qian Zaichen Department of CSE The Chinese University of Hong Kong Shatin,N.T., Hong Kong zcqian@cse.chk.ed.hk In this paper, we stdy
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationEXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.
.4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre
More informationVectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2
MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,
More informationE ect Of Quadrant Bow On Delta Undulator Phase Errors
LCLS-TN-15-1 E ect Of Qadrant Bow On Delta Undlator Phase Errors Zachary Wolf SLAC Febrary 18, 015 Abstract The Delta ndlator qadrants are tned individally and are then assembled to make the tned ndlator.
More informationDecoder Error Probability of MRD Codes
Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehighed Zhiyan Yan Department of Electrical
More informationMathematical Analysis of Nipah Virus Infections Using Optimal Control Theory
Jornal of Applied Mathematics and Physics, 06, 4, 099- Pblished Online Jne 06 in SciRes. http://www.scirp.org/jornal/jamp http://dx.doi.org/0.436/jamp.06.464 Mathematical Analysis of Nipah Virs nfections
More informationOn the tree cover number of a graph
On the tree cover nmber of a graph Chassidy Bozeman Minerva Catral Brendan Cook Oscar E. González Carolyn Reinhart Abstract Given a graph G, the tree cover nmber of the graph, denoted T (G), is the minimm
More informationMethods of Design-Oriented Analysis The GFT: A Final Solution for Feedback Systems
http://www.ardem.com/free_downloads.asp v.1, 5/11/05 Methods of Design-Oriented Analysis The GFT: A Final Soltion for Feedback Systems R. David Middlebrook Are yo an analog or mixed signal design engineer
More informationOn relative errors of floating-point operations: optimal bounds and applications
On relative errors of floating-point operations: optimal bonds and applications Clade-Pierre Jeannerod, Siegfried M. Rmp To cite this version: Clade-Pierre Jeannerod, Siegfried M. Rmp. On relative errors
More information3 2D Elastostatic Problems in Cartesian Coordinates
D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates
More informationA Note on Johnson, Minkoff and Phillips Algorithm for the Prize-Collecting Steiner Tree Problem
A Note on Johnson, Minkoff and Phillips Algorithm for the Prize-Collecting Steiner Tree Problem Palo Feofiloff Cristina G. Fernandes Carlos E. Ferreira José Coelho de Pina September 04 Abstract The primal-dal
More informationMove Blocking Strategies in Receding Horizon Control
Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) MAE 5 - inite Element Analysis Several slides from this set are adapted from B.S. Altan, Michigan Technological University EA Procedre for
More informationQUANTILE ESTIMATION IN SUCCESSIVE SAMPLING
Jornal of the Korean Statistical Society 2007, 36: 4, pp 543 556 QUANTILE ESTIMATION IN SUCCESSIVE SAMPLING Hosila P. Singh 1, Ritesh Tailor 2, Sarjinder Singh 3 and Jong-Min Kim 4 Abstract In sccessive
More informationWe automate the bivariate change-of-variables technique for bivariate continuous random variables with
INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,
More informationConvergence analysis of ant colony learning
Delft University of Technology Delft Center for Systems and Control Technical report 11-012 Convergence analysis of ant colony learning J van Ast R Babška and B De Schtter If yo want to cite this report
More informationRESOLUTION OF INDECOMPOSABLE INTEGRAL FLOWS ON A SIGNED GRAPH
RESOLUTION OF INDECOMPOSABLE INTEGRAL FLOWS ON A SIGNED GRAPH BEIFANG CHEN, JUE WANG, AND THOMAS ZASLAVSKY Abstract. It is well-known that each nonnegative integral flow of a directed graph can be decomposed
More informationIII. Demonstration of a seismometer response with amplitude and phase responses at:
GG5330, Spring semester 006 Assignment #1, Seismometry and Grond Motions De 30 Janary 006. 1. Calibration Of A Seismometer Using Java: A really nifty se of Java is now available for demonstrating the seismic
More information1. State-Space Linear Systems 2. Block Diagrams 3. Exercises
LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space
More informationEXPT. 5 DETERMINATION OF pk a OF AN INDICATOR USING SPECTROPHOTOMETRY
EXPT. 5 DETERMITIO OF pk a OF IDICTOR USIG SPECTROPHOTOMETRY Strctre 5.1 Introdction Objectives 5.2 Principle 5.3 Spectrophotometric Determination of pka Vale of Indicator 5.4 Reqirements 5.5 Soltions
More informationFEA Solution Procedure
EA Soltion rocedre (demonstrated with a -D bar element problem) MAE - inite Element Analysis Many slides from this set are originally from B.S. Altan, Michigan Technological U. EA rocedre for Static Analysis.
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationSimulation investigation of the Z-source NPC inverter
octoral school of energy- and geo-technology Janary 5 20, 2007. Kressaare, Estonia Simlation investigation of the Z-sorce NPC inverter Ryszard Strzelecki, Natalia Strzelecka Gdynia Maritime University,
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More informationRELIABILITY ASPECTS OF PROPORTIONAL MEAN RESIDUAL LIFE MODEL USING QUANTILE FUNC- TIONS
RELIABILITY ASPECTS OF PROPORTIONAL MEAN RESIDUAL LIFE MODEL USING QUANTILE FUNC- TIONS Athors: N.UNNIKRISHNAN NAIR Department of Statistics, Cochin University of Science Technology, Cochin, Kerala, INDIA
More informationScheduling parallel jobs to minimize the makespan
J Sched (006) 9:433 45 DOI 10.1007/s10951-006-8497-6 Schedling parallel jobs to minimize the makespan Berit Johannes C Springer Science + Bsiness Media, LLC 006 Abstract We consider the NP-hard problem
More informationQueueing analysis of service deferrals for load management in power systems
Qeeing analysis of service deferrals for load management in power systems Andrés Ferragt and Fernando Paganini Universidad ORT Urgay Abstract With the advent of renewable sorces and Smart- Grid deployments,
More informationElements of Coordinate System Transformations
B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and
More information4.2 First-Order Logic
64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore
More informationCONTENTS. INTRODUCTION MEQ curriculum objectives for vectors (8% of year). page 2 What is a vector? What is a scalar? page 3, 4
CONTENTS INTRODUCTION MEQ crriclm objectives for vectors (8% of year). page 2 What is a vector? What is a scalar? page 3, 4 VECTOR CONCEPTS FROM GEOMETRIC AND ALGEBRAIC PERSPECTIVES page 1 Representation
More informationREINFORCEMENT LEARNING AND OPTIMAL ADAPTIVE CONTROL
Lewis c11.tex V1-10/19/2011 4:10pm Page 461 11 REINFORCEMENT LEARNING AND OPTIMAL ADAPTIVE CONTROL In this book we have presented a variety of methods for the analysis and design of optimal control systems.
More informationA Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units
A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationDecoder Error Probability of MRD Codes
Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehigh.ed Zhiyan Yan Department of Electrical
More informationMath 116 First Midterm October 14, 2009
Math 116 First Midterm October 14, 9 Name: EXAM SOLUTIONS Instrctor: Section: 1. Do not open this exam ntil yo are told to do so.. This exam has 1 pages inclding this cover. There are 9 problems. Note
More informationRoy Aleksan Centre d Etudes Nucleaires, Saclay, DAPNIA/SPP, F Gif-sur-Yvette, CEDEX, France
CSUHEP -1 DAPNIA--95 LAL--83 Measring the Weak Phase γ in Color Allowed B DKπ Decays Roy Aleksan Centre d Etdes Ncleaires, Saclay, DAPNIA/SPP, F-91191 Gif-sr-Yvette, CEDEX, France Troels C. Petersen LAL
More informationNonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator
Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational
More informationDecision Oriented Bayesian Design of Experiments
Decision Oriented Bayesian Design of Experiments Farminder S. Anand*, Jay H. Lee**, Matthew J. Realff*** *School of Chemical & Biomoleclar Engineering Georgia Institte of echnology, Atlanta, GA 3332 USA
More informationConstructive Root Bound for k-ary Rational Input Numbers
Constrctive Root Bond for k-ary Rational Inpt Nmbers Sylvain Pion, Chee Yap To cite this version: Sylvain Pion, Chee Yap. Constrctive Root Bond for k-ary Rational Inpt Nmbers. 19th Annal ACM Symposim on
More informationHome Range Formation in Wolves Due to Scent Marking
Blletin of Mathematical Biology () 64, 61 84 doi:1.16/blm.1.73 Available online at http://www.idealibrary.com on Home Range Formation in Wolves De to Scent Marking BRIAN K. BRISCOE, MARK A. LEWIS AND STEPHEN
More informationAdaptive Fault-tolerant Control with Control Allocation for Flight Systems with Severe Actuator Failures and Input Saturation
213 American Control Conference (ACC) Washington, DC, USA, Jne 17-19, 213 Adaptive Falt-tolerant Control with Control Allocation for Flight Systems with Severe Actator Failres and Inpt Satration Man Wang,
More informationKonyalioglu, Serpil. Konyalioglu, A.Cihan. Ipek, A.Sabri. Isik, Ahmet
The Role of Visalization Approach on Stdent s Conceptal Learning Konyaliogl, Serpil Department of Secondary Science and Mathematics Edcation, K.K. Edcation Faclty, Atatürk University, 25240- Erzrm-Trkey;
More informationAn Investigation into Estimating Type B Degrees of Freedom
An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information
More informationPath-SGD: Path-Normalized Optimization in Deep Neural Networks
Path-SGD: Path-Normalized Optimization in Deep Neral Networks Behnam Neyshabr, Rslan Salakhtdinov and Nathan Srebro Toyota Technological Institte at Chicago Department of Compter Science, University of
More information1 Undiscounted Problem (Deterministic)
Lectre 9: Linear Qadratic Control Problems 1 Undisconted Problem (Deterministic) Choose ( t ) 0 to Minimize (x trx t + tq t ) t=0 sbject to x t+1 = Ax t + B t, x 0 given. x t is an n-vector state, t a
More informationB-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Quadratic Optimization Problems in Continuous and Binary Variables
B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Qadratic Optimization Problems in Continos and Binary Variables Naohiko Arima, Snyong Kim and Masakaz Kojima October 2012,
More informationIntegration of Basic Functions. Session 7 : 9/23 1
Integration o Basic Fnctions Session 7 : 9/3 Antiderivation Integration Deinition: Taking the antiderivative, or integral, o some nction F(), reslts in the nction () i ()F() Pt simply: i yo take the integral
More informationLINEAR COMBINATIONS AND SUBSPACES
CS131 Part II, Linear Algebra and Matrices CS131 Mathematics for Compter Scientists II Note 5 LINEAR COMBINATIONS AND SUBSPACES Linear combinations. In R 2 the vector (5, 3) can be written in the form
More informationKrauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.
Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore
More informationPREDICTABILITY OF SOLID STATE ZENER REFERENCES
PREDICTABILITY OF SOLID STATE ZENER REFERENCES David Deaver Flke Corporation PO Box 99 Everett, WA 986 45-446-6434 David.Deaver@Flke.com Abstract - With the advent of ISO/IEC 175 and the growth in laboratory
More informationFrequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1
Freqency Estimation, Mltiple Stationary Nonsinsoidal Resonances With Trend 1 G. Larry Bretthorst Department of Chemistry, Washington University, St. Lois MO 6313 Abstract. In this paper, we address the
More informationFREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS
7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,
More information3.1 The Basic Two-Level Model - The Formulas
CHAPTER 3 3 THE BASIC MULTILEVEL MODEL AND EXTENSIONS In the previos Chapter we introdced a nmber of models and we cleared ot the advantages of Mltilevel Models in the analysis of hierarchically nested
More informationStudy on the Mathematic Model of Product Modular System Orienting the Modular Design
Natre and Science, 2(, 2004, Zhong, et al, Stdy on the Mathematic Model Stdy on the Mathematic Model of Prodct Modlar Orienting the Modlar Design Shisheng Zhong 1, Jiang Li 1, Jin Li 2, Lin Lin 1 (1. College
More information