fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /10/10

Size: px
Start display at page:

Download "fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /10/10"

Transcription

1 2 Peri Nes Peer Marwedel TU Dormund, Informaik /0/0 Grahics: Alexandra Nole, Gesine Marwedel, 2003

2 Generalizaion of daa flow: Comuaional grahs Examle: Peri nes Inroduced in 962 by Carl Adam Peri in his PhD hesis. Focus on modeling causal deendencies; no global synchronizaion assumed (message assing only). Key elemens: Condiions Eiher me or no me. Evens May ake lace if cerain condiions are me. Flow relaion Relaes condiions and evens. Condiions, evens and he flow relaion form a biarie grah (grah wih wo kinds of nodes). 2,

3 Examle: Synchronizaion a single rack rail segmen Precondiions 2,

4 Playing he oken game 2, 2008 use normal view mode! - 4 -

5 Conflic for resource rack 2,

6 More comlex examle () Thalys rains beween Cologne, Amserdam, Brussels and Paris. [h:// 2,

7 More comlex examle (2) Slighly simlified: Synchronizaion a Brussels and Paris, using saions Gare du Nord and Gare de Lyon a Paris s 2, 2008 use normal view mode! - 7 -

8 Condiion/even nes Def.: N=(C,E,F) is called a ne, iff he following holds 2. C and E are disjoin ses 3. F (C E) (E C); is binary relaion, ( flow relaion ) Def.: Le N be a ne and le x (C E). x := {y y F x} is called he se of recondiions. x := {y x F y} is called he se of oscondiions. Examle: x x x 2,

9 Loos and ure nes Def.: Le (c,e) C E. (c,e) is called a loo iff cfe efc. Def.: Ne N=(C,E,F) is called ure, if F does no conain any loos. 2,

10 Simle nes Def.: A ne is called simle if no wo ransiions and 2 have he same ses of inu and ouu laces. Examle (no a simle ne): Def.: Simle nes wih no isolaed elemens meeing some addiional resricions are called condiion/even nes (C/E nes). 2,

11 Place/ransiion nes Def.: (P, T, F, K, W, M 0 ) is called a lace/ransiion ne iff 2. N=(P,T,F) is a ne wih laces P and ransiions T 3. K: P (N 0 {ω}) \{0} denoes he caaciy of laces (ω symbolizes infinie caaciy) 4. W: F (N 0 \{0}) denoes he weigh of grah edges 5. M 0 : P N 0 {ω} reresens he iniial marking of laces W (Segmen of some ne) M 0 defauls: K = ω W = 2,

12 Comuing changes of markings Firing ransiions generae new markings on each of he laces according o he following rules: 2,

13 Transiion is acivaed iff Acivaed ransiions Acivaed ransiions can ake lace or fire, bu don have o. We never alk abou ime in he conex of Peri nes. The order in which acivaed ransiions fire, is no fixed (i is non-deerminisic). 2,

14 Shorhand for changes of markings Slide 2: Le ( ) W (, )if \ + W (, )if \ W (, ) + W (, )if 0 = P: M () = M()+ () M = M+ +: vecor add 2,

15 Marix N describing all changes of markings W (, )if \ + W (, )if \ ( ) = W (, ) + W (, )if 0 Def.: Marix N of ne N is a maing N: P T Z (inegers) such ha T: N(,)=() Comonen in column and row indicaes he change of he marking of lace if ransiion akes lace. For ure nes, (N, M 0 ) is a comlee reresenaion of a ne. 2,

16 - 6-2, 2008 Examle: N = s

17 Place - invarians Sandardized echnique for roving roeries of sysem models For any ransiion j T we are looking for ses R P of laces for which he accumulaed marking is consan: Examle: R j ( ) = 0 2,

18 Characerisic Vecor R j ( ) = 0 Le: c R ( ) = if 0 if R R 0 = ( ) = j R P j ( ) c R ( ) = j c R Scalar roduc 2,

19 Condiion for lace invarians R j ( ) = ( ) c ( ) = c = j j R P R 0 Accumulaed marking consan for all ransiions if n... c c R R = = 0 Equivalen o N T c R = 0 where N T is he ransosed of N 2,

20 - 20-2, 2008 More deailed view of comuaions = ) (... ) ( ) ( ) ( )... (... ) ( )... ( ) ( )... ( n R R R n m m n n c c c Sysem of linear equaions. Soluion vecors mus consis of zeros and ones. Equaions wih mulile unknowns ha mus be inegers called diohanic ( Greek mahemaician Diohanos, ~300 B.C.). Diohanic linear equaion sysem more comlex o solve han sandard sysem of linear equaions (acually NP-hard)) Differen echniques for solving equaion sysem (manual,..)

21 Alicaion o Thalys examle N T c R = 0, wih N T = Soluions? Educaed guessing: rows =0 linear deendency among rows rank = 0- = 9 Dimension of soluion sace = 3 - rank = 4 4 comonens of (6,, 2, 3) of c R are indeenden se one of hese o and ohers o 0 o obain a basis for he soluion sace 2,

22 s basis Se one of comonens (6,, 2, 3) o, ohers o 0. s basis b : b (s 6 )=, b (s )=0, b (s 2 )=0, b (s 3 )=0 0 (s 0 ) b (s 0 ) + 0 (s ) b (s ) + 0 (s 3 ) b (s 3 ) = 0 b (s 0 ) = 0 9 (s 9 ) b (s 9 ) + 9 (s 0 ) b (s 0 ) = 0 b (s 9 )=0 b = (,,,,,,0,0,0,0,0,0,0) All comonens {0, } c R = b 2,

23 Inerreaion of he s invarian c R, = ( ) Characerisic vecor describes laces for Cologne rain. We roved ha: he number of rains along he ah remains consan. C R, s 2,

24 2 nd basis Se one of comonens (6,, 2, 3) o, ohers o 0. 2nd basis b 2 : b 2 (s 6 )=0, b 2 (s )=, b 2 (s 2 )=0, b 2 (s 3 )=0 0 (s 0 ) b 2 (s 0 ) + 0 (s ) b 2 (s ) + 0 (s 3 ) b 2 (s 3 ) = 0 b 2 (s 0 ) = 9 (s 9 ) b 2 (s 9 ) + 9 (s 0 ) b 2 (s 0 ) = 0 b 2 (s 9 )= b 2 = (0,-,-,-,0,0,0,0,,,,0,0) b = (,,,,,,0,0,0,0,0,0,0) 2, 2008 b 2 no a characerisic vecor, bu c R,2 =b +b 2 is c R,2 = (,0,0,0,,,0,0,,,,0,0-24 -

25 Inerreaion of he 2 nd invarian c R, 2 = (,0,0,0,,,0,0,,,,0,0) We roved ha: None of he Amserdam rains ges los (nice o know ). C R,2 s 2,

26 Seing b 3 (s 2 ) o and b 4 (s 3 ) o leads o an addiional 2 invarians c R, = c R, 2 = c R, = 3 ( ) ( ) ( ) C R,2 C R,3 C R, c R, 4 = ( ) s We roved ha: he number of rains serving Amserdam, Cologne and Paris remains consan. he number of rain drivers remains consan. C R,4 2,

27 Alicaions Modeling of resources; modeling of muual exclusion; modeling of synchronizaion. 2,

28 Predicae/ransiion nes Goal: comac reresenaion of comlex sysems. Key changes: Tokens are becoming individuals; Transiions enabled if funcions a incoming edges rue; Individuals generaed by firing ransiions defined hrough funcions Changes can be exlained by folding and unfolding C/E nes, semanics can be defined by C/E nes. 2,

29 Examle: Dining hilosohers roblem n> hilosohers siing a a round able; n forks, n laes wih saghei; hilosohers eiher hinking or eaing saghei (using lef and righ fork). 2 forks needed! How o model conflic for forks? How o guaranee avoiding sarvaion? 2,

30 Condiion/even ne model of he dining hilosohers roblem Le x {..3} x : x is hinking e x : x is eaing f x : fork x is available Model quie clumsy. Difficul o exend o more hilosohers. 2, 2008 Normal view mode!

31 Predicae/ransiion model of he dining hilosohers roblem () Le x be one of he hilosohers, le l(x) be he lef soon of x, le r(x) be he righ soon of x. 2 3 f f2 f3 Tokens: individuals. Semanics can be defined by relacing ne by equivalen condiion/even ne. 2,

32 Predicae/ransiion model of he dining hilosohers roblem (2) x 2 3 x v l(x) r(x) f f f2 f3 l(x) r(x) u x e x Model can be exended o arbirary numbers of eole. 2, 2008 use normal view mode!

33 Evaluaion Pros: Aroriae for disribued alicaions, Well-known heory for formally roving roeries, Iniially a quie bizarre oic, bu now acceed due o increasing number of disribued alicaions. Cons (for he nes resened) : roblems wih modeling iming, no rogramming elemens, no hierarchy. Exensions: Enormous amouns of effors on removing limiaions. 2, 2008 back o full screen mode

34 Summary Peri nes: focus on causal relaionshis Condiion/even nes Single oken er lace Place/ransiion nes Mulile okens er lace Predicae/ransiion nes Tokens become individuals Dining hilosohers used as an examle Exensions required o ge around limiaions 2,

fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /11/09

fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /11/09 2 Peri Nes Peer Marwedel TU Dormund, Informaik 2 2009//09 Grahics: Alexandra Nole, Gesine Marwedel, 2003 Models of comuaion considered in his course Communicaion/ local comuaions Undefined comonens Communicaing

More information

Petri Nets. Peter Marwedel TU Dortmund, Informatik /05/13 These slides use Microsoft clip arts. Microsoft copyright restrictions apply.

Petri Nets. Peter Marwedel TU Dortmund, Informatik /05/13 These slides use Microsoft clip arts. Microsoft copyright restrictions apply. 2 Peri Nes Peer Marwedel TU Dormund, Informaik 2 Grahics: Alexandra Nole, Gesine Marwedel, 2003 20/05/3 These slides use Microsof cli ars. Microsof coyrigh resricions aly. Models of comuaion considered

More information

Petri Nets. Peter Marwedel TU Dortmund, Informatik 年 10 月 31 日. technische universität dortmund. fakultät für informatik informatik 12

Petri Nets. Peter Marwedel TU Dortmund, Informatik 年 10 月 31 日. technische universität dortmund. fakultät für informatik informatik 12 Sringer, 2 2 Peri Nes Peer Marwedel TU Dormund, Informaik 2 22 年 月 3 日 These slides use Microsof cli ars. Microsof coyrigh resricions aly. Models of comuaion considered in his course Communicaion/ local

More information

Embedded Systems 4. Petri nets. Introduced in 1962 by Carl Adam Petri in his PhD thesis. Different Types of Petri nets known

Embedded Systems 4. Petri nets. Introduced in 1962 by Carl Adam Petri in his PhD thesis. Different Types of Petri nets known Embedded Sysems 4 - - Peri nes Inroduced in 962 by Carl Adam Peri in his PhD hesis. Differen Tyes of Peri nes known Condiion/even nes Place/ransiion nes Predicae/ransiion nes Hierachical Peri nes, - 2

More information

Embedded Systems 5 BF - ES - 1 -

Embedded Systems 5 BF - ES - 1 - Embedded Sysems 5 - - REVIEW: Peri Nes Def.: N=C,E,F) is called a Peri ne, iff he following holds. C and E are disjoin ses 2. F C E) E C); is binary relaion, flow relaion ) Def.: Le N be a ne and le x

More information

Embedded Systems CS - ES

Embedded Systems CS - ES Embedded Sysems - - Overview of embedded sysems design REVIEW - 2 - REVIEW - 3 - REVIEW - 4 - REVIEW - 5 - Scheduling rocesses in ES: Differences in goals REVIEW In classical OS, qualiy of scheduling is

More information

Embedded Systems 5. Midterm, Thursday December 18, 2008, Final, Thursday February 12, 2009, 16-19

Embedded Systems 5. Midterm, Thursday December 18, 2008, Final, Thursday February 12, 2009, 16-19 Embedded Sysems 5 - - Exam Daes / egisraion Miderm, Thursday December 8, 8, 6-8 Final, Thursday February, 9, 6-9 egisraion hrough HISPOS oen in arox. week If HISPOS no alicable Non-CS, Erasmus, ec send

More information

Embedded Systems CS - ES

Embedded Systems CS - ES Embedded Sysems - - REVIEW Peri Nes - 2 - Comuing changes of markings REVIEW Firing ransiions generae new markings on each of he laces according o he following rules: When a ransiion fires from a marking

More information

fakultät für informatik informatik 12 technische universität dortmund Petri nets Peter Marwedel Informatik 12 TU Dortmund Germany

fakultät für informatik informatik 12 technische universität dortmund Petri nets Peter Marwedel Informatik 12 TU Dortmund Germany 12 Petri nets Peter Marwedel Informatik 12 TU Dortmund Germany Introduction Introduced in 1962 by Carl Adam Petri in his PhD thesis. Focus on modeling causal dependencies; no global synchronization assumed

More information

I. Introduction to place/transition nets. Place/Transition Nets I. Example: a vending machine. Example: a vending machine

I. Introduction to place/transition nets. Place/Transition Nets I. Example: a vending machine. Example: a vending machine Inroducory Tuorial I. Inroducion o place/ransiion nes Place/Transiion Nes I Prepared by: Jörg Desel, Caholic Universiy in Eichsä and Karsen Schmid, Humbold-Universiä zu Berlin Speaker: Wolfgang Reisig,

More information

Chapter 3 Common Families of Distributions

Chapter 3 Common Families of Distributions Chaer 3 Common Families of Disribuions Secion 31 - Inroducion Purose of his Chaer: Caalog many of common saisical disribuions (families of disribuions ha are indeed by one or more arameers) Wha we should

More information

Computing with diode model

Computing with diode model ECE 570 Session 5 C 752E Comuer Aided Engineering for negraed Circuis Comuing wih diode model Objecie: nroduce conces in numerical circui analsis Ouline: 1. Model of an examle circui wih a diode 2. Ouline

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES FORECASS GENERAING FOR ARCH-GARCH PROCESSES USING HE MALAB PROCEDURES Dušan Marček, Insiue of Comuer Science, Faculy of Philosohy and Science, he Silesian Universiy Oava he Faculy of Managemen Science

More information

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011)

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011) ame ID number Mah 8- Sring 5 PRACTICE FIAL EXAM (modified from Mah 8 final exam, Aril 9, ) This exam is closed-book and closed-noe You may use a scienific calculaor, bu no one which is caable of grahing

More information

Curves and Surfaces. Chap. 8 Intro. to Computer Graphics, Spring 2009, Y. G. Shin

Curves and Surfaces. Chap. 8 Intro. to Computer Graphics, Spring 2009, Y. G. Shin Curves and Surfaces Cha. 8 Inro. o Comuer Grahics, Sring 9, Y. G. Shin Reresenaion of Curves and Surfaces Key words surface modeling arameric surface coninuiy, conrol oins basis funcions Bezier curve B-sline

More information

Synthesis of Concurrent Programs Based on Supervisory Control

Synthesis of Concurrent Programs Based on Supervisory Control Synhesis of Concurren Programs Based on Suervisory Conrol Marian V. Iordache School of Engineering and Eng. Tech. LeTourneau Universiy Longview, TX 75607, USA E-mail: MarianIordache@leu.edu Panos J. Ansaklis

More information

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2 7 Parameric equaions This chaer will show ou how o skech curves using heir arameric equaions conver arameric equaions o Caresian equaions find oins of inersecion of curves and lines using arameric equaions

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT Zajac, Budny ON DETERMINATION O SOME CHARACTERISTICS O SEMI MARKOV PROCESS OR DIERENT DISTRIBUTIONS O R&RATA # 2(3 ar 2 (Vol. 2 29, June ON DETERMINATION O SOME CHARACTERISTICS O SEMI-MARKOV PROCESS OR

More information

G. =, etc.

G. =, etc. Maerial Models υ υ3 0 0 0 υ υ 3 0 0 0 υ3 υ3 0 0 0 = 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 3 l (9..4) he subscris denoe he maerial axes, i.e., υ = υ and = (9..5) i j xi xj ii xi Since l is symmeric υ υ =, ec.

More information

Decentralized Control of Petri Nets

Decentralized Control of Petri Nets Decenralized Conrol of Peri Nes Marian V. Iordache and Panos J. Ansaklis Dearmen of Elecrical Engineering Universiy of Nore Dame Nore Dame, IN 6556 iordache.@nd.edu June, 00 M.V. Iordache and P.J. Ansaklis,

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito 1 Inermediae Macroeconomics: Mid-erm exam May 30 h, 2016 Makoo Saio Try he following hree roblems, and submi your answer in handwrien A4 aers. You are execed o dro your aers ino he mailbox assigned for

More information

Description of the MS-Regress R package (Rmetrics)

Description of the MS-Regress R package (Rmetrics) Descriion of he MS-Regress R ackage (Rmerics) Auhor: Marcelo Perlin PhD Suden / ICMA Reading Universiy Email: marceloerlin@gmail.com / m.erlin@icmacenre.ac.uk The urose of his documen is o show he general

More information

Reduction of the Supervisor Design Problem with Firing Vector Constraints

Reduction of the Supervisor Design Problem with Firing Vector Constraints wih Firing Vecor Consrains Marian V. Iordache School of Engineering and Eng. Tech. LeTourneau Universiy Longview, TX 75607-700 Panos J. Ansaklis Dearmen of Elecrical Engineering Universiy of Nore Dame

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

Embedded Systems 6 REVIEW. Place/transition nets. defaults: K = ω W = 1

Embedded Systems 6 REVIEW. Place/transition nets. defaults: K = ω W = 1 Embedded Systems 6-1 - Place/transition nets REVIEW Def.: (P, T, F, K, W, M 0 ) is called a place/transition net (P/T net) iff 1. N=(P,T,F) is a net with places p P and transitions t T 2. K: P (N 0 {ω})

More information

10. State Space Methods

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

More information

2. Nonlinear Conservation Law Equations

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

More information

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise: Mah 225-4 Week 3 April 9-3 EP 7.6 - convoluions; 6.-6.2 - eigenvalues, eigenvecors and diagonalizabiliy; 7. - sysems of differenial equaions. Mon Apr 9 EP 7.6 Convoluions and Laplace ransforms. Announcemens:

More information

EXERCISES FOR SECTION 1.5

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

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

APPLIED ECONOMETRIC TIME SERIES (2nd edition)

APPLIED ECONOMETRIC TIME SERIES (2nd edition) INSTRUCTOR S RESOURCE GUIDE TO ACCOMPANY APPLIED ECONOMETRIC TIME SERIES (2nd ediion) Waler Enders Universiy of Alabama Preared by Pin Chung American Reinsurance Comany and Iowa Sae Universiy Waler Enders

More information

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at BOX-JENKINS MODEL NOAION he Box-Jenkins ARMA(,q) model is denoed b he equaion + + L+ + a θ a L θ a 0 q q. () he auoregressive (AR) ar of he model is + L+ while he moving average (MA) ar of he model is

More information

Announcements: Warm-up Exercise:

Announcements: Warm-up Exercise: Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

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

More information

Decentralized Control of Petri Nets

Decentralized Control of Petri Nets Marian Iordache, Panos J. Ansaklis, Decenralized Conrol of Peri Nes, Proceedings of he Worksho on Discree Even Sysems Conrol a he h Inernaional Conference on Alicaion and Theory of Peri Nes (ATPN 00),.-58,

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate. Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since

More information

Study and Evaluation of Performances of the Digital Multimeter

Study and Evaluation of Performances of the Digital Multimeter Sudy and Evaluaion of Performances of he Digial Mulimeer. Inroducion Mulimeers allow measuring he mos common elecrical quaniies: volage, curren, resisance. They evolved from simle ones based on elecromechanical

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

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

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

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

20. Applications of the Genetic-Drift Model

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

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively: XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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

More information

On a problem of Graham By E. ERDŐS and E. SZEMERÉDI (Budapest) GRAHAM stated the following conjecture : Let p be a prime and a 1,..., ap p non-zero re

On a problem of Graham By E. ERDŐS and E. SZEMERÉDI (Budapest) GRAHAM stated the following conjecture : Let p be a prime and a 1,..., ap p non-zero re On a roblem of Graham By E. ERDŐS and E. SZEMERÉDI (Budaes) GRAHAM saed he following conjecure : Le be a rime and a 1,..., a non-zero residues (mod ). Assume ha if ' a i a i, ei=0 or 1 (no all e i=0) is

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

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

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

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y Review - Quiz # 1 (1) Solving Special Tpes of Firs Order Equaions I. Separable Equaions (SE). d = f() g() Mehod of Soluion : 1 g() d = f() (The soluions ma be given implicil b he above formula. Remember,

More information

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1.

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1. Timed Circuis Asynchronous Circui Design Chris J. Myers Lecure 7: Timed Circuis Chaper 7 Previous mehods only use limied knowledge of delays. Very robus sysems, bu exremely conservaive. Large funcional

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

Reducing Probabilistic Timed Petri Nets for Asynchronous Architectural Analysis

Reducing Probabilistic Timed Petri Nets for Asynchronous Architectural Analysis Reducing Probabilisic Timed Peri Nes for Asynchronous Archiecural Analysis Sangyun Kim, Sunan Tugsinavisu Universiy of Souhern California Los Angeles, CA 90089, USA {sangyun,ugsinav}@eiger.usc.edu Peer

More information

Wavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves

Wavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves Wavele Mehods for Time Series Analysis Par I: Inroducion o Waveles and Wavele Transforms waveles are analysis ools for ime series and images as a subjec, waveles are relaively new (983 o presen) a synhesis

More information

1 Review of Zero-Sum Games

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

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Appendix 14.1 The optimal control problem and its solution using

Appendix 14.1 The optimal control problem and its solution using 1 Appendix 14.1 he opimal conrol problem and is soluion using he maximum principle NOE: Many occurrences of f, x, u, and in his file (in equaions or as whole words in ex) are purposefully in bold in order

More information

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

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

More information

A Method for Deadlock Prevention in Discrete Event Systems Using Petri Nets

A Method for Deadlock Prevention in Discrete Event Systems Using Petri Nets A Mehod for Deadlock Prevenion in Discree Even Sysems Using Peri Nes Technical Reor of he ISIS Grou a he Universiy of Nore Dame ISIS-99-006 July, 999 Marian V. Iordache John O. Moody Panos J. Ansaklis

More information

A Method for Deadlock Prevention in Discrete Event Systems Using Petri Nets

A Method for Deadlock Prevention in Discrete Event Systems Using Petri Nets M.V. Iordache, J.O. Moody and P.J. Ansaklis, "A Mehod for Deadlock Prevenion in Discree Even Sysems Using Peri Nes," Technical Reor isis-99-006, De. of Elecrical Engr., Univ. of Nore Dame, July 999. A

More information

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag. Saionary Processes Sricly saionary - The whole join disribuion is indeenden of he dae a which i is measured and deends only on he lag. - E y ) is a finie consan. ( - V y ) is a finie consan. ( ( y, y s

More information

Online Convex Optimization Example And Follow-The-Leader

Online Convex Optimization Example And Follow-The-Leader CSE599s, Spring 2014, Online Learning Lecure 2-04/03/2014 Online Convex Opimizaion Example And Follow-The-Leader Lecurer: Brendan McMahan Scribe: Sephen Joe Jonany 1 Review of Online Convex Opimizaion

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Non-uniform circular motion *

Non-uniform circular motion * OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform

More information

f(t) dt, x > 0, is the best value and it is the norm of the

f(t) dt, x > 0, is the best value and it is the norm of the MATEMATIQKI VESNIK 66, 1 (214), 19 32 March 214 originalni nauqni rad research aer GENERALIZED HAUSDORFF OPERATORS ON WEIGHTED HERZ SPACES Kuang Jichang Absrac. In his aer, we inroduce new generalized

More information

Transform Techniques. Moment Generating Function

Transform Techniques. Moment Generating Function Transform Techniques A convenien way of finding he momens of a random variable is he momen generaing funcion (MGF). Oher ransform echniques are characerisic funcion, z-ransform, and Laplace ransform. Momen

More information

Rise-Time Distortion of Signal without Carrying Signal

Rise-Time Distortion of Signal without Carrying Signal Journal of Physics: Conference Series PAPER OPEN ACCESS Rise-Time Disorion of Signal wihou Carrying Signal To cie his aricle: N S Bukhman 6 J. Phys.: Conf. Ser. 738 8 View he aricle online for udaes and

More information

Decentralized Control of Petri Nets

Decentralized Control of Petri Nets Marian Iordache, Panos J.Ansaklis, Decenralized onrol of Peri Nes, ISIS Technical Reor ISIS-00-005, Ocober 00. Decenralized onrol of Peri Nes Technical Reor of he ISIS Grou a he Universiy of Nore Dame

More information

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he

More information

EE123 Digital Signal Processing

EE123 Digital Signal Processing Discree Transforms (Finie) EE3 Digial Signal Processing Lecure 9 DFT is only one ou of a LARGE class of ransforms Used for: Analysis Comression Denoising Deecion Recogniion Aroximaion (Sarse) Sarse reresenaion

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum Theory Probab. Al. Vol. 45, No.,, (5-36) Research Reor No. 45, 999, De. Theore. Sais. Aarhus Soing Brownian Moion wihou Aniciaion as Close as Possible o is Ulimae Maximum S. E. GRAVERSEN 3, G. PESKIR 3,

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

Notes on Kalman Filtering

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

More information

The following report makes use of the process from Chapter 2 in Dr. Cumming s thesis.

The following report makes use of the process from Chapter 2 in Dr. Cumming s thesis. Zaleski 1 Joseph Zaleski Mah 451H Final Repor Conformal Mapping Mehods and ZST Hele Shaw Flow Inroducion The Hele Shaw problem has been sudied using linear sabiliy analysis and numerical mehods, bu a novel

More information

Math 2214 Sol Test 2B Spring 2015

Math 2214 Sol Test 2B Spring 2015 Mah 14 Sol Tes B Sring 015 roblem 1: An objec weighing ounds sreches a verical sring 8 fee beond i naural lengh before coming o res a equilibrium The objec is ushed u 6 fee from i s equilibrium osiion

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

Chapter 2 Summary. Carnegie Learning

Chapter 2 Summary. Carnegie Learning Chaper Summary Key Terms inducion (.1) deducion (.1) counerexample (.1) condiional saemen (.1) proposiional form (.1) proposiional variables (.1) hypohesis (.1) conclusion (.1) ruh value (.1) ruh able

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

Math 315: Linear Algebra Solutions to Assignment 6

Math 315: Linear Algebra Solutions to Assignment 6 Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen

More information

Innova Junior College H2 Mathematics JC2 Preliminary Examinations Paper 2 Solutions 0 (*)

Innova Junior College H2 Mathematics JC2 Preliminary Examinations Paper 2 Solutions 0 (*) Soluion 3 x 4x3 x 3 x 0 4x3 x 4x3 x 4x3 x 4x3 x x 3x 3 4x3 x Innova Junior College H Mahemaics JC Preliminary Examinaions Paper Soluions 3x 3 4x 3x 0 4x 3 4x 3 0 (*) 0 0 + + + - 3 3 4 3 3 3 3 Hence x or

More information

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Two Coupled Oscillators / Normal Modes

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

More information

Solutions of Ecological models by Homotopy-Perturbation Method (HPM)

Solutions of Ecological models by Homotopy-Perturbation Method (HPM) Ques Journals Journal of Research in Alied Mahemaics Volume ~ Issue 0 (0) : 7- ISSN(Online) : 9-07 ISSN (Prin):9-075 www.quesjournals.org Research Paer Soluions of Ecological models by Homooy-Perurbaion

More information