Linear Program for Partially Observable Markov Decision Processes. MS&E 339B June 9th, 2004 Erick Delage

Size: px
Start display at page:

Download "Linear Program for Partially Observable Markov Decision Processes. MS&E 339B June 9th, 2004 Erick Delage"

Transcription

1 Linea Pogam fo Patiall Obsevable Makov Decision Pocesses MS&E 339B June 9th 2004 Eick Delage

2 Intoduction Patiall Obsevable Makov Decision Pocesses Etension of the Makov Decision Pocess to a wold with uncetaint about the cuent state Double objective decisions: Confim state vs. Reduce epected immediate cost Fomulation as a continuous state space Makov Decision Pocess Methods to appoimate the Value function Value function appoimation using Linea Pogamming Etension of LP appoach to continuous state-space MDP Application to POMDP Implementation fo small poblems

3 Patiall Obsevable Makov Decision Pocess Descibed b: a set of states a set of actions a set of obsevations a set of tansition pobabilities a set of obsevation pobabilities a set of costs a discount facto an initial belief the hbid DP Value function S { s1 s2... ss} A { a1 a2... aa} Z { z1 z2... zz} si a sj p sj si a O zi a sj p zi sj R gs:s A α 01 [ ] p0 s J p min{ Ε [ α g k µ p µ p U 0 ~p0 k 0 k k k+ 1 ]}

4 POMDP 2

5 POMDP 3 In bief: MDP with uncetaint about the state MDP is a specific case of POMDP Applications: Robot Navigation A obot gathes infomation about its position and decides of movements to get to destination A obot gathes infomation about a moving objects position and decides of pope movements to get close to it. Dialogue Modeling A machine needs to ask the ight questions to lean what is a use s intentions Limits of Scalabilit: Cuse of dimensionalit : state space is the S -1 dimension pobabilit simple Cuse of histo : the value function gets moe comple with depth of hoizon

6 POMDP as a Continuous state Makov Decision Pocess POMDP can be fomulated : a conve state space B b s p s a set of actions A { a1 a2... aa} a set of tansition pobabilities bi a bj p bjbi a a cost function gb: B R a discount facto α 01 [ ] an initial belief b0 Belief Space: Discete pobabilit distibution ove the set of envionment state S S -dimensional pobabilit simple which vetices ae beliefs whee the state is known with 100% cetaint Cost function on the belief space is the epected immediate cost g b b a Ε[ gs s a] b s gs s a s S Continuous state-space Bellman s equation J b min{ Ε[ g b µ b + αj b]} π

7 Beliefs ansition pobabilities bi a bj p bjbi a p zibi a if bj s p s a s~b i s Using Bae s ule: p za p zs p s a s S p a s S p zs p s a p zs a O z a s p s a p s s a p s s S p z s a p zs a p s a p s a z p za p za

8 Some Solutions - Heuistics Most likel state Dnamic Pogam Noubakhsh et al J Ε 0 min k s s α g sk µ sk π k0 Polic will be satisfied with making most pobable the lowest cost state In fact tue immediate cost is which penalizes uncetaint Bi-Objective algoithm Cetaint vs. Rewad Cassanda et al Limits the entop level on the wa to the lowest cost state H p s N k 1 p slog2 p s aps sstem in states that ae suounded b high uncetaint J gb b Ε[ gs s] π p s s p s Jsagma p s p s agmin Ε[ H p sa z] if H p s > κ π p s a othewise

9 Some Solutions - Appoimations Gid Based Appoimaton Hauskecht 2000 N discete points in the belief space the value function evaluated at these points an intepolation function fo all othe points Vˆ B ˆ { bˆˆ 1b2...ˆ bn} Vˆ b: Bˆ R V b Int b Vˆ 1 b min{ gb b a + E[ Int b Vˆ k]} b Bˆ k+ α π b PCA & Eponential Famil PCA compession of state space Ro 2003 Assumes tajectoies fo most lage eal wold POMDP lie nea a low dimensional manifold embedded in the belief space Finds an m-dimension manifold m<<n ove which to poject the belief space Discetise the m-dimension manifold Use Gid based Appoimation ove the manifold

10 Linea Pogam Appoimation Appoach J min Ε[ g µ + αj ] µ U J min{ g µ + αε[ J ]} In the case of POMDP if we assume Z {z1z2 } is a finite set then: owads fomulating the Linea Pogam: Minimize Subject to µ U J min { g µ + α p µ J d } J µ U min { g µ + α µ U c J g µ + p Y µ p Y µ µ J } J J µ µ

11 LP Appoimation Appoach 2 Minimize Subject to ] [ 2 1 J k d c p g Y µ µ µ µ c d c d c d c d c. 2 1 Σ Y Y Y p p p µ µ µ µ µ µ µ

12 LP Appoimation Appoach 3 Minimize Subject to c [ Σ. ] g µ µ µ Computing the integals Could be computed analticall Becomes moe pactical to use c as a distibution to appoimate the integals Sampling the constaints Discete ALP concepts suggests to use c as the andom distibution of samples Numbe of constaints equied popotional to numbe of featue functions Choosing the featue functions Integal must make sense ove the pobabilit simple Since value function is concave intuition suggests using a quadatic function

13 Choosing c C is a measue of whee the LP should concentate its effots Poblem dependent fo best esults In POMDP not all beliefs ae visited C should eflect this

14 Implementation Small poblems Sstem had ~10 states ~10-D belief space Featue functions: 1000 linea constaints + C distibution of beliefs when behaving andoml Pocessing time: less than 1 minute [ 2 2L ] 0 1 s 2 s S s S S

15 Conclusion Appoimate Linea Pogam etends easil to the continuous state space Simple paametes to deal with : featue functions & c Most of pocessing done b the linea pogam solve Computation should gow lineal with numbe of featue functions Numbe of featue functions should gow linea with numbe of states S

16 Futue Wok Appl Continuous ALP to lage scale POMDPs Compae pefomances with othe techniques Confim bounds on compleit gowth Stud bounds on optimal polic pefomance Stud acceptable class of featue functions

17 Refeences Cassanda A. R. Kaelbling L. and Kuien J. A Acting unde uncetaint: Discete Baesian models fo mobile-obot navigation. In Poceedings of the IEEE/RSJ Intenational Confeence on Intelligent Robotic Sstems IROS. De Faias D.P. Van Ro B he Linea Pogamming Appoach to Appoimate Dnamic Pogamming. In Opeations Reseach 51:6. Hauskecht M Value-function appoimations fo patiall obsevable Makov decision pocesses. Jounal of Atificial Intelligence Reseach 13: Hauskecht M. and Kveton B Linea pogam appoimations fo factoed continuous-state Makov decision pocesses. In Neual Infomation Pocessing Sstems. Pineau J. Godon G. hun S Point-based value iteation: An antime algoithm fo POMDPs. Intenational Joint Confeence on Atificial Intelligence Noubakhsh I. Powes R. and Bichfield S.1995 Devish an office-navigation obot. AI Magazine 162: Ro N Finding Appoimate POMDP Solutions hough Belief Compession. Ph.D. thesis Canegie Mellon Univesit Pittsbugh PA.

A Markov Decision Approach for the Computation of Testability of RTL Constructs

A Markov Decision Approach for the Computation of Testability of RTL Constructs A Makov Decision Appoach fo the Computation of Testability of RTL Constucts José Miguel Fenandes Abstact In the analysis of digital cicuits, to study testability estimation measues, dissipated powe and

More information

1 Explicit Explore or Exploit (E 3 ) Algorithm

1 Explicit Explore or Exploit (E 3 ) Algorithm 2.997 Decision-Making in Lage-Scale Systems Mach 3 MIT, Sping 2004 Handout #2 Lectue Note 9 Explicit Exploe o Exploit (E 3 ) Algoithm Last lectue, we studied the Q-leaning algoithm: [ ] Q t+ (x t, a t

More information

Bayesian Congestion Control over a Markovian Network Bandwidth Process

Bayesian Congestion Control over a Markovian Network Bandwidth Process Bayesian Congestion Contol ove a Makovian Netwok Bandwidth Pocess Paisa Mansouifad,, Bhaska Kishnamachai, Taa Javidi Ming Hsieh Depatment of Electical Engineeing, Univesity of Southen Califonia, Los Angeles,

More information

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India

More information

Temporal-Difference Learning

Temporal-Difference Learning .997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation

More information

Value Prediction with FA. Chapter 8: Generalization and Function Approximation. Adapt Supervised Learning Algorithms. Backups as Training Examples [ ]

Value Prediction with FA. Chapter 8: Generalization and Function Approximation. Adapt Supervised Learning Algorithms. Backups as Training Examples [ ] Chapte 8: Genealization and Function Appoximation Objectives of this chapte:! Look at how expeience with a limited pat of the state set be used to poduce good behavio ove a much lage pat.! Oveview of function

More information

Chapter 8: Generalization and Function Approximation

Chapter 8: Generalization and Function Approximation Chapte 8: Genealization and Function Appoximation Objectives of this chapte: Look at how expeience with a limited pat of the state set be used to poduce good behavio ove a much lage pat. Oveview of function

More information

Goodness-of-fit for composite hypotheses.

Goodness-of-fit for composite hypotheses. Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test

More information

MONTE CARLO STUDY OF PARTICLE TRANSPORT PROBLEM IN AIR POLLUTION. R. J. Papancheva, T. V. Gurov, I. T. Dimov

MONTE CARLO STUDY OF PARTICLE TRANSPORT PROBLEM IN AIR POLLUTION. R. J. Papancheva, T. V. Gurov, I. T. Dimov Pliska Stud. Math. Bulga. 14 (23), 17 116 STUDIA MATHEMATICA BULGARICA MOTE CARLO STUDY OF PARTICLE TRASPORT PROBLEM I AIR POLLUTIO R. J. Papancheva, T. V. Guov, I. T. Dimov Abstact. The actual tanspot

More information

International Journal of Mathematical Archive-3(12), 2012, Available online through ISSN

International Journal of Mathematical Archive-3(12), 2012, Available online through  ISSN Intenational Jounal of Mathematical Achive-3(), 0, 480-4805 Available online though www.ijma.info ISSN 9 504 STATISTICAL QUALITY CONTROL OF MULTI-ITEM EOQ MOEL WITH VARYING LEAING TIME VIA LAGRANGE METHO

More information

4/18/2005. Statistical Learning Theory

4/18/2005. Statistical Learning Theory Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse

More information

Bayesian Analysis of Topp-Leone Distribution under Different Loss Functions and Different Priors

Bayesian Analysis of Topp-Leone Distribution under Different Loss Functions and Different Priors J. tat. Appl. Po. Lett. 3, No. 3, 9-8 (6) 9 http://dx.doi.og/.8576/jsapl/33 Bayesian Analysis of Topp-Leone Distibution unde Diffeent Loss Functions and Diffeent Pios Hummaa ultan * and. P. Ahmad Depatment

More information

Information Retrieval Advanced IR models. Luca Bondi

Information Retrieval Advanced IR models. Luca Bondi Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the

More information

INFORMATION GEOMETRY OF PROPAGATION ALGORITHMS AND APPROXIMATE INFERENCE

INFORMATION GEOMETRY OF PROPAGATION ALGORITHMS AND APPROXIMATE INFERENCE 2nd Intenational Symposium on Infomation Geomety and its Applications Decembe 12-16, 2005, Tokyo Pages 000 000 INFORMATION GEOMETRY OF PROPAGATION ALGORITHMS AND APPROXIMATE INFERENCE SHIRO IKEDA 1 intoduction

More information

Regularization. Stephen Scott and Vinod Variyam. Introduction. Outline. Machine. Learning. Problems. Measuring. Performance.

Regularization. Stephen Scott and Vinod Variyam. Introduction. Outline. Machine. Learning. Problems. Measuring. Performance. leaning can geneally be distilled to an optimization poblem Choose a classifie (function, hypothesis) fom a set of functions that minimizes an objective function Clealy we want pat of this function to

More information

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two

More information

Prediction of Motion Trajectories Based on Markov Chains

Prediction of Motion Trajectories Based on Markov Chains 2011 Intenational Confeence on Compute Science and Infomation Technology (ICCSIT 2011) IPCSIT vol. 51 (2012) (2012) IACSIT Pess, Singapoe DOI: 10.7763/IPCSIT.2012.V51.50 Pediction of Motion Tajectoies

More information

Control Chart Analysis of E k /M/1 Queueing Model

Control Chart Analysis of E k /M/1 Queueing Model Intenational OPEN ACCESS Jounal Of Moden Engineeing Reseach (IJMER Contol Chat Analysis of E /M/1 Queueing Model T.Poongodi 1, D. (Ms. S. Muthulashmi 1, (Assistant Pofesso, Faculty of Engineeing, Pofesso,

More information

Spatial Distribution of Two-Agent Clusters for Efficient Navigation

Spatial Distribution of Two-Agent Clusters for Efficient Navigation Spatial Distibution of Two-Agent Clustes fo Efficient Navigation Jan De Mot Eic Feon Laboatoy fo Infomation and Decision Systems Massachusetts Institute of Technology Cambidge MA 02139 USA {jdemotfeon}@mit.edu

More information

Scalable Reinforcement Learning via Trajectory Optimization and Approximate Gaussian Process Regression

Scalable Reinforcement Learning via Trajectory Optimization and Approximate Gaussian Process Regression Scalable Reinfocement Leaning via Tajectoy Optimization and Appoximate Gaussian Pocess Regession Yunpeng Pan,, Xinyan Yan,, Evangelos Theodoou,, and Byon Boots, Institute fo Robotics and Intelligent Machines,

More information

Adaptive Checkpointing in Dynamic Grids for Uncertain Job Durations

Adaptive Checkpointing in Dynamic Grids for Uncertain Job Durations Adaptive Checkpointing in Dynamic Gids fo Uncetain Job Duations Maia Chtepen, Bat Dhoedt, Filip De Tuck, Piet Demeeste NTEC-BBT, Ghent Univesity, Sint-Pietesnieuwstaat 41, Ghent, Belgium {maia.chtepen,

More information

GENLOG Multinomial Loglinear and Logit Models

GENLOG Multinomial Loglinear and Logit Models GENLOG Multinomial Loglinea and Logit Models Notation This chapte descibes the algoithms used to calculate maximum-likelihood estimates fo the multinomial loglinea model and the multinomial logit model.

More information

Performance Loss Bounds for Approximate Value Iteration with State Aggregation

Performance Loss Bounds for Approximate Value Iteration with State Aggregation MATHEMATICS OF OPERATIONS RESEARCH Vol. 31, No. 2, May 2006, pp. 234 244 issn 0364-765X eissn 1526-5471 06 3102 0234 infoms doi 10.1287/moo.1060.0188 2006 INFORMS Pefomance Loss Bounds fo Appoximate Value

More information

Using Laplace Transform to Evaluate Improper Integrals Chii-Huei Yu

Using Laplace Transform to Evaluate Improper Integrals Chii-Huei Yu Available at https://edupediapublicationsog/jounals Volume 3 Issue 4 Febuay 216 Using Laplace Tansfom to Evaluate Impope Integals Chii-Huei Yu Depatment of Infomation Technology, Nan Jeon Univesity of

More information

A Converse to Low-Rank Matrix Completion

A Converse to Low-Rank Matrix Completion A Convese to Low-Rank Matix Completion Daniel L. Pimentel-Alacón, Robet D. Nowak Univesity of Wisconsin-Madison Abstact In many pactical applications, one is given a subset Ω of the enties in a d N data

More information

Fall 2014 Randomized Algorithms Oct 8, Lecture 3

Fall 2014 Randomized Algorithms Oct 8, Lecture 3 Fall 204 Randomized Algoithms Oct 8, 204 Lectue 3 Pof. Fiedich Eisenband Scibes: Floian Tamè In this lectue we will be concened with linea pogamming, in paticula Clakson s Las Vegas algoithm []. The main

More information

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018 Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining

More information

Solving Some Definite Integrals Using Parseval s Theorem

Solving Some Definite Integrals Using Parseval s Theorem Ameican Jounal of Numeical Analysis 4 Vol. No. 6-64 Available online at http://pubs.sciepub.com/ajna///5 Science and Education Publishing DOI:.69/ajna---5 Solving Some Definite Integals Using Paseval s

More information

3.1 Random variables

3.1 Random variables 3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated

More information

Multiple Criteria Secretary Problem: A New Approach

Multiple Criteria Secretary Problem: A New Approach J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and

More information

Simulation of a 2-link Brachiating Robot with Open-Loop Controllers

Simulation of a 2-link Brachiating Robot with Open-Loop Controllers Simulation of a -link Bachiating Robot with Open-oop Contolles David Uffod Nothwesten Univesit June 009 . Poject Oveview The goal of this poject was to wite a complete simulation of a -link swinging obot

More information

Analytical Expressions for Positioning Uncertainty Propagation in Networks of Robots

Analytical Expressions for Positioning Uncertainty Propagation in Networks of Robots Analytical Expessions fo Positioning Uncetainty Popagation in Netwoks of Robots Ioannis M Rekleitis and Stegios I Roumeliotis Abstact In this pape we pesent an analysis of the positioning uncetainty incease

More information

Computers and Mathematics with Applications

Computers and Mathematics with Applications Computes and Mathematics with Applications 58 (009) 9 7 Contents lists available at ScienceDiect Computes and Mathematics with Applications jounal homepage: www.elsevie.com/locate/camwa Bi-citeia single

More information

AQI: Advanced Quantum Information Lecture 2 (Module 4): Order finding and factoring algorithms February 20, 2013

AQI: Advanced Quantum Information Lecture 2 (Module 4): Order finding and factoring algorithms February 20, 2013 AQI: Advanced Quantum Infomation Lectue 2 (Module 4): Ode finding and factoing algoithms Febuay 20, 203 Lectue: D. Mak Tame (email: m.tame@impeial.ac.uk) Intoduction In the last lectue we looked at the

More information

Mobility of atoms and diffusion. Einstein relation.

Mobility of atoms and diffusion. Einstein relation. Mobility of atoms and diffusion. Einstein elation. In M simulation we can descibe the mobility of atoms though the mean squae displacement that can be calculated as N 1 MS ( t ( i ( t i ( 0 N The MS contains

More information

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012 Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,

More information

Quantum Fourier Transform

Quantum Fourier Transform Chapte 5 Quantum Fouie Tansfom Many poblems in physics and mathematics ae solved by tansfoming a poblem into some othe poblem with a known solution. Some notable examples ae Laplace tansfom, Legende tansfom,

More information

Supporting information

Supporting information Electonic Supplementay Mateial (ESI) fo Physical Chemisty Chemical Physics. This jounal is the Owne Societies 18 Suppoting infomation Nonstoichiometic oxides as a continuous homologous seies: linea fee-enegy

More information

The Substring Search Problem

The Substring Search Problem The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is

More information

Estimation of entropy rate and Rényi entropy rate for Markov chains

Estimation of entropy rate and Rényi entropy rate for Markov chains 06 IEEE Intenational Symposium on Infomation Theoy Estimation of entopy ate and Rényi entopy ate fo Makov chains Sudeep Kamath, Segio Vedú Depatment of Electical Engineeing Pinceton Univesity e-mail: sukamath,

More information

A Simple Model of Communication APIs Application to Dynamic Partial-order Reduction

A Simple Model of Communication APIs Application to Dynamic Partial-order Reduction Simple Model of Communication PIs pplication to Dynamic Patial-ode Reduction Cistian Rosa Stephan Mez Matin Quinson VOCS 2010 22/09/2010 1 / 18 Motivation Distibuted lgoithms ae had to get ight: lack of

More information

New problems in universal algebraic geometry illustrated by boolean equations

New problems in universal algebraic geometry illustrated by boolean equations New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic

More information

Robust optimal policies for Markov decision processes with safety-threshold constraints

Robust optimal policies for Markov decision processes with safety-threshold constraints Robust optimal policies fo Makov decision pocesses with safety-theshold constaints Rayna Dimitova, Jie Fu, Ufuk Topcu Abstact We study the synthesis of obust optimal contol policies fo Makov decision pocesses

More information

Description of TEAM Workshop Problem 33.b: Experimental Validation of Electric Local Force Formulations.

Description of TEAM Workshop Problem 33.b: Experimental Validation of Electric Local Force Formulations. 1 Desciption of TEAM Wokshop Poblem 33.b: Expeimental Validation of Electic Local Foce Fomulations. Olivie Baé, Pascal Bochet, Membe, IEEE Abstact Unde extenal stesses, the esponse of any mateial is a

More information

6 Matrix Concentration Bounds

6 Matrix Concentration Bounds 6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom

More information

Research Article A Case Study on Air Combat Decision Using Approximated Dynamic Programming

Research Article A Case Study on Air Combat Decision Using Approximated Dynamic Programming Mathematical Poblems in Engineeing, Aticle ID 84, pages http://dx.doi.og/.55/4/84 Reseach Aticle A Case Study on Ai Combat Decision Using Appoximated Dynamic Pogamming Yaofei Ma, Xiaole Ma, and Xiao Song

More information

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations MATH 415, WEEK 3: Paamete-Dependence and Bifucations 1 A Note on Paamete Dependence We should pause to make a bief note about the ole played in the study of dynamical systems by the system s paametes.

More information

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in

More information

Conjugate Gradient Methods. Michael Bader. Summer term 2012

Conjugate Gradient Methods. Michael Bader. Summer term 2012 Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent

More information

Reliability analysis examples

Reliability analysis examples Reliability analysis examples Engineeing Risk Analysis Goup, Technische Univesität München. Acisst., 80 Munich, Gemany. May 8, 08 Intoduction In the context of eliability analysis and ae event estimation,

More information

Directed Regression. Benjamin Van Roy Stanford University Stanford, CA Abstract

Directed Regression. Benjamin Van Roy Stanford University Stanford, CA Abstract Diected Regession Yi-hao Kao Stanfod Univesity Stanfod, CA 94305 yihaoao@stanfod.edu Benjamin Van Roy Stanfod Univesity Stanfod, CA 94305 bv@stanfod.edu Xiang Yan Stanfod Univesity Stanfod, CA 94305 xyan@stanfod.edu

More information

J. Electrical Systems 1-3 (2005): Regular paper

J. Electrical Systems 1-3 (2005): Regular paper K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato

More information

Bounds on the performance of back-to-front airplane boarding policies

Bounds on the performance of back-to-front airplane boarding policies Bounds on the pefomance of bac-to-font aiplane boading policies Eitan Bachmat Michael Elin Abstact We povide bounds on the pefomance of bac-to-font aiplane boading policies. In paticula, we show that no

More information

Conservative Averaging Method and its Application for One Heat Conduction Problem

Conservative Averaging Method and its Application for One Heat Conduction Problem Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem

More information

Rejection Based Face Detection

Rejection Based Face Detection Rejection Based Face Detection Michael Elad* Scientific Computing and Computational Mathematics Stanfod Univesit The Compute Vision Wokshop Mach 18 th, 00 * Collaboation with Y. Hel-O and R. Keshet 1 Pat

More information

MONTE CARLO SIMULATION OF FLUID FLOW

MONTE CARLO SIMULATION OF FLUID FLOW MONTE CARLO SIMULATION OF FLUID FLOW M. Ragheb 3/7/3 INTRODUCTION We conside the situation of Fee Molecula Collisionless and Reflective Flow. Collisionless flows occu in the field of aefied gas dynamics.

More information

16 Modeling a Language by a Markov Process

16 Modeling a Language by a Markov Process K. Pommeening, Language Statistics 80 16 Modeling a Language by a Makov Pocess Fo deiving theoetical esults a common model of language is the intepetation of texts as esults of Makov pocesses. This model

More information

Identification of the degradation of railway ballast under a concrete sleeper

Identification of the degradation of railway ballast under a concrete sleeper Identification of the degadation of ailway ballast unde a concete sleepe Qin Hu 1) and Heung Fai Lam ) 1), ) Depatment of Civil and Achitectual Engineeing, City Univesity of Hong Kong, Hong Kong SAR, China.

More information

Anyone who can contemplate quantum mechanics without getting dizzy hasn t understood it. --Niels Bohr. Lecture 17, p 1

Anyone who can contemplate quantum mechanics without getting dizzy hasn t understood it. --Niels Bohr. Lecture 17, p 1 Anyone who can contemplate quantum mechanics without getting dizzy hasn t undestood it. --Niels Boh Lectue 17, p 1 Special (Optional) Lectue Quantum Infomation One of the most moden applications of QM

More information

1D2G - Numerical solution of the neutron diffusion equation

1D2G - Numerical solution of the neutron diffusion equation DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both

More information

Internet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks

Internet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate

More information

On Utilization of K-Means for Determination

On Utilization of K-Means for Determination Jounal of Softwae Engineeing and Applications, 07, 0, 605-64 http://www.scip.og/jounal/jsea ISSN Online: 945-34 ISSN Pint: 945-36 On Utilization of K-Means fo Detemination of q-paamete fo Tsallis-Entopy-Maximized-FCM

More information

A New Approach to General Relativity

A New Approach to General Relativity Apeion, Vol. 14, No. 3, July 7 7 A New Appoach to Geneal Relativity Ali Rıza Şahin Gaziosmanpaşa, Istanbul Tukey E-mail: aizasahin@gmail.com Hee we pesent a new point of view fo geneal elativity and/o

More information

Scientific Computing II

Scientific Computing II Scientific Computing II Conjugate Gadient Methods Michael Bade Summe 2014 Conjugate Gadient Methods, Summe 2014 1 Families of Iteative Solves elaxation methods: Jacobi-, Gauss-Seidel-Relaxation,... Ove-Relaxation-Methods

More information

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events Intenational Jounal of Contempoay Mathematical Sciences Vol. 12, 2017, no. 5, 243-253 HIKARI Ltd, www.m-hikai.com https://doi.og/10.12988/ijcms.2017.7728 Hypothesis Test and Confidence Inteval fo the Negative

More information

PAPER 39 STOCHASTIC NETWORKS

PAPER 39 STOCHASTIC NETWORKS MATHEMATICAL TRIPOS Pat III Tuesday, 2 June, 2015 1:30 pm to 4:30 pm PAPER 39 STOCHASTIC NETWORKS Attempt no moe than FOUR questions. Thee ae FIVE questions in total. The questions cay equal weight. STATIONERY

More information

Research Article Robust Evaluation for Transportation Network Capacity under Demand Uncertainty

Research Article Robust Evaluation for Transportation Network Capacity under Demand Uncertainty Hindawi Jounal of Advanced Tanspotation Volume 2017, Aticle ID 9814909, 11 pages https://doi.og/10.1155/2017/9814909 Reseach Aticle Robust Evaluation fo Tanspotation Netwok Capacity unde Demand Uncetainty

More information

On Computing Optimal (Q, r) Replenishment Policies under Quantity Discounts

On Computing Optimal (Q, r) Replenishment Policies under Quantity Discounts Annals of Opeations Reseach manuscipt No. will be inseted by the edito) On Computing Optimal, ) Replenishment Policies unde uantity Discounts The all - units and incemental discount cases Michael N. Katehakis

More information

A Compositional Approach to Reactive Games under Temporal Logic Specifications

A Compositional Approach to Reactive Games under Temporal Logic Specifications A Compositional Appoach to Reactive Games unde Tempoal Logic Specifications Abhishek Ninad Kulkani and Jie Fu Abstact We study the poblem of compositional synthesis of contolles fo eactive games with linea

More information

Bifurcation Analysis for the Delay Logistic Equation with Two Delays

Bifurcation Analysis for the Delay Logistic Equation with Two Delays IOSR Jounal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 5 Ve. IV (Sep. - Oct. 05), PP 53-58 www.iosjounals.og Bifucation Analysis fo the Delay Logistic Equation with Two Delays

More information

Modeling and Calculation of Optical Amplification in One Dimensional Case of Laser Medium Using Finite Difference Time Domain Method

Modeling and Calculation of Optical Amplification in One Dimensional Case of Laser Medium Using Finite Difference Time Domain Method Jounal of Physics: Confeence Seies PAPER OPEN ACCESS Modeling and Calculation of Optical Amplification in One Dimensional Case of Lase Medium Using Finite Diffeence Time Domain Method To cite this aticle:

More information

MULTILAYER PERCEPTRONS

MULTILAYER PERCEPTRONS Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR

More information

Grouped data clustering using a fast mixture-model-based algorithm

Grouped data clustering using a fast mixture-model-based algorithm Poceedings of the 29 IEEE Intenational Confeence on Systems, Man, and Cybenetics San Antonio, TX, USA - Octobe 29 Gouped data clusteing using a fast mixtue-model-based algoithm Allou SAMÉ Laboatoie des

More information

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc. Numeical solution of diffusion mass tansfe model in adsoption systems Pof., D.Sc. Agenda Mass Tansfe Mechanisms Diffusion Mass Tansfe Models Solving Diffusion Mass Tansfe Models Paamete Estimation 2 Mass

More information

10/04/18. P [P(x)] 1 negl(n).

10/04/18. P [P(x)] 1 negl(n). Mastemath, Sping 208 Into to Lattice lgs & Cypto Lectue 0 0/04/8 Lectues: D. Dadush, L. Ducas Scibe: K. de Boe Intoduction In this lectue, we will teat two main pats. Duing the fist pat we continue the

More information

Nuclear Medicine Physics 02 Oct. 2007

Nuclear Medicine Physics 02 Oct. 2007 Nuclea Medicine Physics Oct. 7 Counting Statistics and Eo Popagation Nuclea Medicine Physics Lectues Imaging Reseach Laboatoy, Radiology Dept. Lay MacDonald 1//7 Statistics (Summaized in One Slide) Type

More information

Explosive Contagion in Networks (Supplementary Information)

Explosive Contagion in Networks (Supplementary Information) Eplosive Contagion in Netwoks (Supplementay Infomation) Jesús Gómez-Gadeñes,, Laua Loteo, Segei N. Taaskin, and Fancisco J. Péez-Reche Institute fo Biocomputation and Physics of Comple Systems (BIFI),

More information

Lecture 8 - Gauss s Law

Lecture 8 - Gauss s Law Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.

More information

763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012

763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012 763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012 1. Continuous Random Walk Conside a continuous one-dimensional andom walk. Let w(s i ds i be the pobability that the length of the i th displacement

More information

FUSE Fusion Utility Sequence Estimator

FUSE Fusion Utility Sequence Estimator FUSE Fusion Utility Sequence Estimato Belu V. Dasaathy Dynetics, Inc. P. O. Box 5500 Huntsville, AL 3584-5500 belu.d@dynetics.com Sean D. Townsend Dynetics, Inc. P. O. Box 5500 Huntsville, AL 3584-5500

More information

A hybrid cross entropy algorithm for solving dynamic transit network design problem

A hybrid cross entropy algorithm for solving dynamic transit network design problem A hybid coss entopy algoithm fo solving dynamic tansit netok design poblem TAI-YU MA Tanspot Economics Laboatoy Univesity Lyon 2 - CNRS Lyon, Fance This pape poposes a hybid multiagent leaning algoithm

More information

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical

More information

Modelling Vortex Fields in Metal Smelting Furnaces

Modelling Vortex Fields in Metal Smelting Furnaces Int. Jnl. of Multiphysics Volume 4 Numbe 4 00 35 Modelling Votex Fields in Metal Smelting Funaces Oleg Kazak, Oleksand Semko Depatment Of Physics, Donetsk National Univesity, Donetsk, Ukaine Email: olegkazak@yandex.u;

More information

Green s function Monte Carlo algorithms for elliptic problems

Green s function Monte Carlo algorithms for elliptic problems Mathematics and Computes in Simulation 63 (2003) 587 604 Geen s function Monte Calo algoithms fo elliptic poblems I.T. Dimov, R.Y. Papancheva Cental Laboatoy fo Paallel Pocessing, Depatment of Paallel

More information

Energy Savings Achievable in Connection Preserving Energy Saving Algorithms

Energy Savings Achievable in Connection Preserving Energy Saving Algorithms Enegy Savings Achievable in Connection Peseving Enegy Saving Algoithms Seh Chun Ng School of Electical and Infomation Engineeing Univesity of Sydney National ICT Austalia Limited Sydney, Austalia Email:

More information

An extended target tracking method with random finite set observations

An extended target tracking method with random finite set observations 4th Intenational Confeence on Infomation Fusion Chicago Illinois USA July 5-8 0 An extended taget tacing method with andom finite set obsevations Hongyan Zhu Chongzhao Han Chen Li Dept. of Electonic &

More information

ASTR415: Problem Set #6

ASTR415: Problem Set #6 ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal

More information

arxiv: v2 [cs.ai] 14 Aug 2014

arxiv: v2 [cs.ai] 14 Aug 2014 A Paallel Algoithm fo Exact Bayesian Stuctue Discovey in Bayesian Netwoks axiv:1408.1664v2 [cs.ai] 14 Aug 2014 Yetian Chen Jin Tian Depatment of Compute Science Iowa State Univesity Ames, IA 50011, USA

More information

OSCILLATIONS AND GRAVITATION

OSCILLATIONS AND GRAVITATION 1. SIMPLE HARMONIC MOTION Simple hamonic motion is any motion that is equivalent to a single component of unifom cicula motion. In this situation the velocity is always geatest in the middle of the motion,

More information

A Stochastic EOQ Policy of Cold-Drink-For a Retailer

A Stochastic EOQ Policy of Cold-Drink-For a Retailer Vietnam Jounal of Mathematics 33:4 005) 437 44 9LHWQDP -RXUQDO RI 0$7+0$7,&6 9$67 A Stochastic EOQ Policy of Cold-Dink-Fo a Retaile Shib Sanka Sana 1 Kipasindhu Chaudhui 1 Depatment of Math., Bhanga Mahavidyalaya

More information

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS Pogess In Electomagnetics Reseach, PIER 73, 93 105, 2007 COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS T.-X. Song, Y.-H. Liu, and J.-M. Xiong School of Mechanical Engineeing

More information

CSCE 478/878 Lecture 4: Experimental Design and Analysis. Stephen Scott. 3 Building a tree on the training set Introduction. Outline.

CSCE 478/878 Lecture 4: Experimental Design and Analysis. Stephen Scott. 3 Building a tree on the training set Introduction. Outline. In Homewok, you ae (supposedly) Choosing a data set 2 Extacting a test set of size > 3 3 Building a tee on the taining set 4 Testing on the test set 5 Repoting the accuacy (Adapted fom Ethem Alpaydin and

More information

Duality between Statical and Kinematical Engineering Systems

Duality between Statical and Kinematical Engineering Systems Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.

More information

Markscheme May 2017 Calculus Higher level Paper 3

Markscheme May 2017 Calculus Higher level Paper 3 M7/5/MATHL/HP3/ENG/TZ0/SE/M Makscheme May 07 Calculus Highe level Pape 3 pages M7/5/MATHL/HP3/ENG/TZ0/SE/M This makscheme is the popety of the Intenational Baccalaueate and must not be epoduced o distibuted

More information

15.081J/6.251J Introduction to Mathematical Programming. Lecture 6: The Simplex Method II

15.081J/6.251J Introduction to Mathematical Programming. Lecture 6: The Simplex Method II 15081J/6251J Intoduction to Mathematical Pogamming ectue 6: The Simplex Method II 1 Outline Revised Simplex method Slide 1 The full tableau implementation Anticycling 2 Revised Simplex Initial data: A,

More information

HINDCASTING OF WIND AND WAVE CLIMATE OF SEAS AROUND RUSSIA

HINDCASTING OF WIND AND WAVE CLIMATE OF SEAS AROUND RUSSIA Leonid J. Lopatoukhin, Alexande V. Boukhanovsky, Ecateina S. Chenysheva, Segey V. Ivanov HINDCASTING OF WIND AND WAVE CLIMATE OF SEAS AROUND RUSSIA St. Petesbug State Univesity, Institute fo High Pefomance

More information

Machine Learning and Rendering

Machine Learning and Rendering East Building, Balloom BC nvidia.com/siggaph2018 Machine Leaning and Rendeing Alex Kelle, Diecto of Reseach Machine Leaning and Rendeing Couse web page at https://sites.google.com/site/mlandendeing/ 14:00

More information

CBE Transport Phenomena I Final Exam. December 19, 2013

CBE Transport Phenomena I Final Exam. December 19, 2013 CBE 30355 Tanspot Phenomena I Final Exam Decembe 9, 203 Closed Books and Notes Poblem. (20 points) Scaling analysis of bounday laye flows. A popula method fo measuing instantaneous wall shea stesses in

More information

Homework Set 3 Physics 319 Classical Mechanics

Homework Set 3 Physics 319 Classical Mechanics Homewok Set 3 Phsics 319 lassical Mechanics Poblem 5.13 a) To fin the equilibium position (whee thee is no foce) set the eivative of the potential to zeo U 1 R U0 R U 0 at R R b) If R is much smalle than

More information

LEARNING PARAMETRIC CLOSED-LOOP POLICIES

LEARNING PARAMETRIC CLOSED-LOOP POLICIES LARNING PARAMTRIC CLOSD-LOOP POLICIS FOR MARKOV POTNTIAL GAMS Segio Valcacel Macua PROWLR.io Cambidge, UK segio@powle.io Javie Zazo, Santiago Zazo Infomation Pocessing and Telecommunications Cente Univesidad

More information