Bucknell University Using ODE45 MATLAB Help
|
|
- Elaine Harvey
- 6 years ago
- Views:
Transcription
1 Bucknell University Using ODE45 MATLAB Help MATLAB's standard slver fr rdinary differential equatins (ODEs) is the functin de45. This functin implements a Runge-Kutta methd with a variable time step fr efficient cmputatin. de45 id designed t handle the fllwing general prblem dy = f(, t y) y( t ) = y [] dt where t is the independent variable (time, psitin, vlume) and y is a vectr f dependent variables (temperature, psitin, cncentratins) t be fund. The mathematical prblem is specified when the vectr f functins n the right-hand side f Eq. [], f(, t y), is set and the initial cnditins, y = y at time t, are specified. The ntes here apply t versins f MATLAB abve 5.0 and cver the basics f using the functin de45. Fr mre infrmatin n this and ther ODE slvers in MATLAB, see the n-line help. Cntents: Syntax fr de45... Integrating a single, first-rder equatin... 3 Getting the slutin at particular values f the independent variable... 4 Integrating a set f cupled first-rder equatins... 4 Integrating a secnd-rder initial-value prblem (IVP)... 7 Integrating an Nth-rder initial-value prblem... 8 Changing mdel parameters... 9 Integrating a secnd-rder bundary-value prblem (BVP)... Setting ptins in de45... Ging beynd de ENGR0 Using ODE45
2 Syntax fr de45 de45 may be invked frm the cmmand line via [t,y] = de45('fname', tspan, y0, pts) where fname name f the functin Mfile used t evaluate the right-hand-side functin in Eq. [] at a given value f the independent variable and dependent variable(s) (string). The functin definitin line usually has the frm functin dydt = fname(t,y) The utput variable (dydt) must be a vectr with the same size as y. Nte that the independent variable (t here) must be included in the input argument list even if it des nt explicitly appear in the expressins used t generate dydt. tspan y0 pts t y -element vectr defining the range f integratin ([t tf]) thugh variatins are pssible. vectr f initial cnditins fr the dependent variable. There shuld be as many initial cnditins as there are dependent variables. a MATLAB structure variable that allws yu t cntrl the details f cmputatin (if yu want t). This argument is ptinal and, if nt prvided, de45 will use default values (see the examples belw). Value f the independent variable at which the slutin array (y) is calculated. Nte that by default this will nt be a unifrmly distributed set f values. Values f the slutin t the prblem (array). Each clumn f y is a different dependent variable. The size f the array is length(t)-by-length(y0) Specific examples f using de45 nw fllw. Mfiles fr these examples are in the bdy f this dcument and shuld als be available in the flder that cntains this dcument. ENGR0 Using ODE45
3 Integrating a single, first-rder equatin The height f fluid in a tank (h(t)) whse utlet flw is dependent n the pressure head (height f fluid) inside the tank and whse inlet flw is a functin f time may be mdeled via the equatin dh dt = α() t β h h( 0 ) = h [] Find the slutin, h(t), fr 0< t < 30 if the fllwing values fr the parameters are given. Input flw: α( t) = sin( t) β = h = Step : Identify f(, t y) and write a MATLAB functin Mfile t evaluate it. In this case, we have time as the independent variable and the tank height as the (single) dependent variable. Thus, we have f(, t y) f(, t h) = α() t β h [3] Frm the given infrmatin fr this prblem, the required Mfile, named tankfill.m, is functin dhdt = tankfill(t,h) % RHS functin fr tank-fill prblem A = 0 + 4*sin(t); H = *sqrt(h); dhdt = A - H; % alpha(t) % beta*sqrt(h) % ef - tankfill.m Step : Use de45 t slve the prblem The initial cnditin has the height at fr t = 0 and we want t integrate until t = 30. The fllwing set f cmmands shw explicitly hw the slutin is put tgether. >> tspan = [0 30]; (integratin range) >> h0 = ; (initial cnditin, h(0)) >> [t,h] = de45('tankfill',tspan,h0); (slve the prblem) Step 3: Lk at the slutin The slutin can be viewed via the plt cmmand as in >> plt(t,h) The "curve" is a little chppy thugh it is accurate t the default relative tlerance (0.00). Nte that the places where the slutin is given are nt unifrmly spread ut. See the next sectin fr imprving appearances. ENGR0 Using ODE45 3
4 Getting the slutin at particular values f the independent variable de45 uses a variable-step-length algrithm t find the slutin fr a given ODE. Thus, de45 varies the size f the step f the independent variable in rder t meet the accuracy yu specify at any particular pint alng the slutin. If de45 can take "big" steps and still meet this accuracy, it will d s and will therefre mve quickly thrugh regins where the slutin des nt "change" greatly. In regins where the slutin changes mre rapidly, de45 will take "smaller" steps. While this strategy is gd frm an efficiency r speed pint f view, it means that the slutin des nt appear at a fixed set f values fr the independent variable (as a fixed-step methd wuld) and smetimes the slutin curves lk a little ragged. The simplest way t imprve n the density f slutin pints is t mdify the input tspan frm a -element vectr t an N-element vectr via smething like >> tspan = linspace(t,tf,500) ; and use this new versin in the input list t de45. Smther curves can als be generated by interplatin (spline interplatin usually wrks nicely). Fr example, if yu wanted a smther result frm the slutin fr the tank-fill prblem, yu might d the fllwing >> ti = linspace(tspan(),tspan(),300); (300 pints - yu culd use mre) >> hi = spline(t,h,ti); >> plt(t,h,,ti,hi); The interplated curve smthes ut the rugh edges caused by simply cnnecting the data pints (which is what plt des) and s makes the graph mre appealing, in a visual sense. Integrating a set f cupled first-rder equatins Chemical-kinetics prblems ften lead t sets f cupled, first-rder ODEs. Fr example, cnsider the reactin netwrk A B C [4] Assuming a first-rder reactin-rate expressin fr each transfrmatin, material balances fr each species lead t the fllwing set f ODEs: da = ka + kb dt db = ka kb kb 3 dt dc = kb 3 dt [5] with the initial cnditins, A( 0) = A, B( 0) = B, C( 0) = C. Since the equatins are cupled, yu cannt slve each ne separately and s must slve them simultaneusly. The system in Eq. [5] can be put in the standard frm fr de45 (Eq. []) by defining the vectrs y, y and f as ENGR0 Using ODE45 4
5 A A + ky ky y= B y( 0) = y = B f( t, y) = ky ( k+ k3) y C C ky 3 [6] Slving the system represented by Eq. [6] is a simple extensin f what was dne fr slving a single equatin. We'll demnstrate the slutin fr the fllwing situatin k = 5 k = k = A = B = C = 0 3 Step : Write a functin Mfile t evaluate the right-hand-side expressin The primary difference here, cmpared t the single-equatin case, is that the input variable y will be a vectr. The first element f y represents the cncentratin f species A at a time t, and the secnd and third elements representing the cncentratins f species B and C, respectively, at the same time, t. This rdering f variables is defined by Eq. [6]. There is n "right" rder t the variables but whatever rder yu d chse, use it cnsistently. We'll call the Mfile react.m. It lks like this: functin dydt = react(t,y) % Slve the kinetics example dydt = zers(size(y)); % Parameters - reactin-rate cnstants k = 5; k = ; k3 = ; A = y(); B = y(); C = y(3); We'll be explicit abut it here thugh yu can d the calculatins directly with the y-values. % Evaluate the RHS expressin dydt() = -k*a + k*b; dydt() = k*a - (k+k3)*b; dydt(3) = k3*b; % ef - react.m Nte that the input arguments must be t and y (in that rder) even thugh t is nt explicitly used in the functin. Step : Use de45 t slve the prblem N time interval is given s we'll pick ne (0 t 4) and see what the slutin lks like. If a lnger r shrter interval is needed, we can simply re-execute the functin with a new value fr the ending time. Fllwing the utline fr the single-equatin prblem, the call t de45 is, >> [t,y] = de45('react',[0 4],[ 0 0]); Nte that the initial cnditin is prvided directly in the call t de45. Yu culd als have defined a variable y0 prir t the call t de45 and used that variable as an input. ENGR0 Using ODE45 5
6 Take a mment t lk at the utputs. The number f pints at which the slutin is knwn is >> length(t) Als cnsider the shape f the utput variable y: >> size(y) Is the result as stated abve (i.e., is it length(t)-by-length(y0))? Step 3: Lk at the slutin If yu want t see the time-curse f all species, use the cmmand >> plt(t,y) The blue line will be the first clumn f y (species A). The green and red lines will be the secnd and third clumns f y (species B and C, respectively). If yu wanted t lk at nly ne species (fr example, species B), yu wuld give the cmmand >> plt(t,y(:,)) since the secnd clumn f y hlds the infrmatin n species B. ENGR0 Using ODE45 6
7 Integrating a secnd-rder initial-value prblem (IVP) A mass-spring-dashpt system can be mdeled via the fllwing secnd-rder ODE y + cy + ω y= g() t y( 0) = y, v( 0) = y ( 0 ) = v [7] In this mdel, c represents a retarding frce prprtinal t the velcity f the mass, ω is the natural frequency f the system and g(t) is the frcing (r input) functin. The initial cnditins are the initial psitin (y ) and initial velcity (v ). de45 is set up t handle nly first-rder equatins and s a methd is needed t cnvert this secndrder equatin int ne (r mre) first-rder equatins which are equivalent. The cnversin is accmplished thrugh a technique called "reductin f rder". We'll illustrate the slutin fr the particular set f cnditins c= 5 ω = y( 0) = v( 0) = 0 g() t = sin( t) Step: Define the cmpnents f a vectr p = [ p p ] T as fllws: p = y p = y [8] Step : Frm the first derivatives f each f the cmpnents f p Using the given differential equatin, we can write a system f first-rder equatins as p = y = p p = y = g() t cy ω y = gt () cp ω p [9] In writing the expressin fr the secnd cmpnent, we've used the gverning ODE (Eq. [7]). Step 3: Cast the prblem in the frmat needed t use de45. dp d p (, p) p = = f(, p) = () = p f t dt dt p gt cp p f (, t p) = t ω [0] Step 4: Cllect the initial cnditins int a single vectr ( 0) ( 0) p( 0) = p = p ( 0) = y ( 0) = y p y v [] Step 5: Apply de45 t slve the system f equatins The Mfile fr the RHS functin fr this prblem will be called spring.m. Here it is: functin pdt = spring(t,p) % Spring example prblem ENGR0 Using ODE45 7
8 % Parameters - damping cefficient and natural frequency c = 5; w = ; g = sin(t); % frcing functin pdt = zers(size(p)); pdt() = p(); pdt() = g - c*p() - (w^)*p(); % ef - spring.m The call t de45 is, fr a slutin interval f 0 t 0, >> p0 = [ 0]; (initial psitin and velcity) >> [t,p] = de45('spring',[0 0],p0); Step 6: Lk at the results If yu wanted t lk at nly the displacement, yu'd want t lk at the first clumn f p (see the definitin f p in the first step, Eq. [8]). Hence, yu wuld give the cmmand >> plt(t,p(:,)) An interesting plt fr these srts f prblems is the phase-plane plt, a plt f the velcity f the mass versus its psitin. This plt is easily created frm yur slutin via >> plt(p(:,),p(:,)) Phase-plane plts are useful in analyzing general features f dynamic systems. Integrating an Nth-rder initial-value prblem T use de45 t integrate an Nth-rder ODE, yu simply cntinue the prcess utlined in the sectin n integrating a nd-rder ODE. The first element f the vectr p is set t the dependent variable and then subsequent elements are defined as the derivatives f the dependent variable up t ne less than the rder f the equatin. Finally, the initial cnditins are cllected int ne vectr t give the frmat presented in Eq. []. Fr example, the 4th-rder equatin wuld generate the first-rder system a d 4 y b dy 3 c dy d dy + ey = 0 [] 4 3 dx dx dx dx p p d p p3 = dt p3 p4 p ( bp + cp + dp + ep )/ a [3] which, alng with an apprpriate set f initial cnditins wuld cmplete the set-up fr de45. ENGR0 Using ODE45 8
9 Changing mdel parameters In all the examples given abve, the parameter values were given specific variables in the Mfile used t evaluate the RHS functin (the mdel f the system). This is fine fr ne-sht cases and in instances where yu dn't anticipate a desire t change the parameters. Hwever, this situatin is nt fine where yu want t be able t change the parameters (e.g., change the damping cefficient in rder t see the result in a phase-plane plt). One apprach t changing parameters is t simply edit the file every time yu want t make a change. While having the advantage f simplicity, this apprach suffers frm inflexibility, especially as the number f parameters and as the frequency f the changes increase. T get ther parameters int the functin, yu need t use an expanded versin f the syntax fr de45, ne that allws ther infrmatin t be prvided t the derivative functin when de45 uses that functin. This is mst easily seen by an example (and by reading the n-line help n de45). Step : Write the Mfile fr the RHS functin s that it allws mre input variables. The parameter list starts after the dependent variable and after a required input called flag. Fr this example, we will re-write spring.m s that c and w are given their values via the functin definitin line. The altered Mfile is functin pdt = spring(t,p,flag,c,w) % Spring example prblem % Parameters c is the damping cefficient and % w is the natural frequency pdt = zers(size(p)); g = sin(t); % frcing functin pdt() = p(); pdt() = g - c*p() - (w^)*p(); % ef spring.m Step : Write a driver script that implements yur lgic and allws yu t set values f the parameters fr the prblem. The script created here will d the fllwing. Implement a lp that a. asks fr values f the parameters b. slves the ODE c. plts the phase-plane view f the slutin. Exits if n inputs are given Certainly mre sphisticated scripts are pssible but this has the essence f the idea. The script is called dspring.m and it is: %DOSPRING Interactive pltting f the phase-plane while % infinite lp C_SPRING = input('damping cefficient [c]: '); ENGR0 Using ODE45 9
10 end if isempty(c_spring) % hw t get ut break end W_SPRING = input('natural frequency [w]: '); [t,p] = de45('spring',[0 0],[ 0],[],C_SPRING,W_SPRING); plt(p(:,),p(:,)) title('phase-plane plt f results') xlabel('psitin') ylabel('velcity') % es - dspring.m Nte the additins t the call t de45. First, a placehlder fr the ptins input is inserted (an empty array) s that default ptins are used. Then, the parameters fr the mdel are prvided in the rder that they appear in the definitin line f the RHS functin. Try the script ut and mdify it (e.g., yu culd add the frequency and/r amplitude f the frcing functin as smething t be changed). ENGR0 Using ODE45 0
11 Integrating a secnd-rder bundary-value prblem (BVP) de45 was written t slve initial-value prblems (IVPs). Hence it cannt be (directly) used t slve, fr example, the fllwing prblem derived frm a mdel f heat-transfer in a rd: d y dy y = 0 y(0) = = 0 [4] dx dx x= since the value f the derivative at x = 0 is nt specified (it is knwn at x =, thugh). Equatin [4] is a bundary-value prblem (BVP) and is cmmn in mdels based n transprt phenmena (heat transfer, mass transfer and fluid mechanics). All is nt lst because ne way t slve a BVP is t pretend it is an IVP. T make up fr the lack f knwledge f the derivative at the initial pint, yu can guess a value, d the integratin and then check yurself by seeing hw clse yu are t meeting the cnditins at the ther end f the interval. When yu have guessed the right starting values, yu have the slutin t the prblem. This apprach is smetimes called the "shting methd" by analgy t the ballistics prblem f landing an artillery shell n a target by specifying nly it's set-up (pwder charge and angle f the barrel). Step : Set up the prblem s that de45 can slve it Using the apprach f turning a secnd rder equatin int a pair f cupled first-rder equatins, we have d p p dx p = p 0 = p( ) [5] v where v has been used t represent the (unknwn) value f the derivative at x = 0. The Mfile used t evaluate the RHS is as fllws functin dpdx = htrd(x,p) % Ht-rd prblem illustrating the shting methd dpdx = zers(size(p)); dpdx() = p(); dpdx() = p(); % ef - htrd.m Step : Guess a value f the initial slpe and integrate t x = The prblem will be iterative s it's nt likely that the first guess will be right. Frm the physics f the prblem, the end f the rd (at x = ) will be clder than the place we are starting frm (x = 0) and s we'll guess a negative value fr the initial slpe. >> v = -; >> [x,p] = de45('htrd',[0 ],[ v]); The value f the derivative at x = is the last value in the secnd clumn f p (why?). Thus, we can check the accuracy f the first guess via ENGR0 Using ODE45
12 >> p(length(x),) which I fund t be That s t lw (it shuld be zer). Step 3: Iterate until the bundary cnditin at x = is met Yu can use brute frce here if yu have nly ne prblem r yu culd finesse it by hking the whle thing up t fzer and have fzer d the guessing. Here are my brute-frce results: Value f v Slpe at x = The trend is bvius and s the initial slpe is arund (the exact value is -tanh() = ). Using fzer wuld be a gd alternative if this prblem were t be slved many times ver. Step 4: Lk at the results Even thugh we are guessing the initial slpe t slve the prblem, it is the slutin, y(x), that we are really interested in. This slutin is in the first clumn f p and may be viewed via >> plt(x,p(:,)) Setting ptins in de45 The input pts is a MATLAB structure variable that cn be used t cntrl the perfrmance f the varius ODE-slvers in MATLAB. The mst cmmn ptin that yu ll likely want t alter is the accuracy t which slutins are cmputed. T make this prcess easy, a pair f functins are available deset fr creating and changing ptins and deget fr displaying infrmatin n ptins. T see what the current settings are, try the cmmand >> deset Default values fr any setting are dented by the braces, {}. MATLAB uses tw accuracy measures fr slving ODEs the relative tlerance (RelTl in pts) and the abslute tlerance (AbsTl in pts). Each step in the integratin is taken s that it satisfies the cnditin Errr at step j max( RelTl y, AbsTl ) k where the subscript k ranges ver all the cmpnents f the slutin vectr at time step j. T alter the default settings, use cmmands such as >> ldopts = deset; >> newopts = deset(ldopts, RelTl,e-6) Infrmatin n the settings fr the ther ptins is available in the n-line help. jk k ENGR0 Using ODE45
13 Ging beynd de45 The slver de45 is nt the be-all and end-all f ODE-slvers. While de45 shuld be yur first chice fr integratin, there are prblems that the functin perfrms prly n r even fails n. In such cases, there are fallback slvers that are available. All these slvers use the same syntax as de45 (see page ) but have ptins fr handling mre difficult r sphisticated prblems. Here are sme suggestins fr handling nn-standard ODE prblems: If accuracy yu desire is nt btainable via de45, try the functin de3. This slver uses a variable rder methd that may be able t imprve ver what de45 des. If de45 is taking t lng t cmpute a slutin, yur prblem may be stiff (i.e., it invlves a system with a wide range f time cnstants). Try the functin de5s. If yur system f equatins has the frm d M y = f(, t y) dt where M is a (typically nn-singular) matrix, try the functin de5s. Yu ll find mre infrmatin n these functin in the n-line help and dcumentatin. Fr example, try the n-line functin reference (available thrugh the cmmand helpdesk) n any f the slvers nted abve. ENGR0 Using ODE45 3
Kinetic Model Completeness
5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins
More informationCS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007
CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is
More informationDifferentiation Applications 1: Related Rates
Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm
More informationActivity Guide Loops and Random Numbers
Unit 3 Lessn 7 Name(s) Perid Date Activity Guide Lps and Randm Numbers CS Cntent Lps are a relatively straightfrward idea in prgramming - yu want a certain chunk f cde t run repeatedly - but it takes a
More informationFall 2013 Physics 172 Recitation 3 Momentum and Springs
Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.
More informationSections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.
Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage
More information[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )
(Abut the final) [COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t m a k e s u r e y u a r e r e a d y ) The department writes the final exam s I dn't really knw what's n it and I can't very well
More informationRevision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax
.7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationMaterials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion
Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals f Diffusin Diffusin: Transprt in a slid, liquid, r gas driven by a cncentratin gradient (r, in the case f mass transprt, a chemical ptential
More informationAP Physics Kinematic Wrap Up
AP Physics Kinematic Wrap Up S what d yu need t knw abut this mtin in tw-dimensin stuff t get a gd scre n the ld AP Physics Test? First ff, here are the equatins that yu ll have t wrk with: v v at x x
More information20 Faraday s Law and Maxwell s Extension to Ampere s Law
Chapter 20 Faraday s Law and Maxwell s Extensin t Ampere s Law 20 Faraday s Law and Maxwell s Extensin t Ampere s Law Cnsider the case f a charged particle that is ming in the icinity f a ming bar magnet
More informationWe can see from the graph above that the intersection is, i.e., [ ).
MTH 111 Cllege Algebra Lecture Ntes July 2, 2014 Functin Arithmetic: With nt t much difficulty, we ntice that inputs f functins are numbers, and utputs f functins are numbers. S whatever we can d with
More informationNUMBERS, MATHEMATICS AND EQUATIONS
AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS NUMBERS, MATHEMATICS AND EQUATIONS An integral part t the understanding f ur physical wrld is the use f mathematical mdels which can be used t
More informationThis section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.
Sectin 3.2: Many f yu WILL need t watch the crrespnding vides fr this sectin n MyOpenMath! This sectin is primarily fcused n tls t aid us in finding rts/zers/ -intercepts f plynmials. Essentially, ur fcus
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 11: Mdeling with systems f ODEs In Petre Department f IT, Ab Akademi http://www.users.ab.fi/ipetre/cmpmd/ Mdeling with differential equatins Mdeling strategy Fcus
More informationFigure 1a. A planar mechanism.
ME 5 - Machine Design I Fall Semester 0 Name f Student Lab Sectin Number EXAM. OPEN BOOK AND CLOSED NOTES. Mnday, September rd, 0 Write n ne side nly f the paper prvided fr yur slutins. Where necessary,
More informationPhysics 2010 Motion with Constant Acceleration Experiment 1
. Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin
More informationCHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.
MATH 1342 Ch. 24 April 25 and 27, 2013 Page 1 f 5 CHAPTER 24: INFERENCE IN REGRESSION Chapters 4 and 5: Relatinships between tw quantitative variables. Be able t Make a graph (scatterplt) Summarize the
More informationChapter 2 GAUSS LAW Recommended Problems:
Chapter GAUSS LAW Recmmended Prblems: 1,4,5,6,7,9,11,13,15,18,19,1,7,9,31,35,37,39,41,43,45,47,49,51,55,57,61,6,69. LCTRIC FLUX lectric flux is a measure f the number f electric filed lines penetrating
More informationMODULE FOUR. This module addresses functions. SC Academic Elementary Algebra Standards:
MODULE FOUR This mdule addresses functins SC Academic Standards: EA-3.1 Classify a relatinship as being either a functin r nt a functin when given data as a table, set f rdered pairs, r graph. EA-3.2 Use
More informationPreparation work for A2 Mathematics [2017]
Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 2: Mdeling change. In Petre Department f IT, Åb Akademi http://users.ab.fi/ipetre/cmpmd/ Cntent f the lecture Basic paradigm f mdeling change Examples Linear dynamical
More information, which yields. where z1. and z2
The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin
More informationPlan o o. I(t) Divide problem into sub-problems Modify schematic and coordinate system (if needed) Write general equations
STAPLE Physics 201 Name Final Exam May 14, 2013 This is a clsed bk examinatin but during the exam yu may refer t a 5 x7 nte card with wrds f wisdm yu have written n it. There is extra scratch paper available.
More informationTuring Machines. Human-aware Robotics. 2017/10/17 & 19 Chapter 3.2 & 3.3 in Sipser Ø Announcement:
Turing Machines Human-aware Rbtics 2017/10/17 & 19 Chapter 3.2 & 3.3 in Sipser Ø Annuncement: q q q q Slides fr this lecture are here: http://www.public.asu.edu/~yzhan442/teaching/cse355/lectures/tm-ii.pdf
More informationAP Statistics Notes Unit Two: The Normal Distributions
AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).
More informationWRITING THE REPORT. Organizing the report. Title Page. Table of Contents
WRITING THE REPORT Organizing the reprt Mst reprts shuld be rganized in the fllwing manner. Smetime there is a valid reasn t include extra chapters in within the bdy f the reprt. 1. Title page 2. Executive
More informationMODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b
. REVIEW OF SOME BASIC ALGEBRA MODULE () Slving Equatins Yu shuld be able t slve fr x: a + b = c a d + e x + c and get x = e(ba +) b(c a) d(ba +) c Cmmn mistakes and strategies:. a b + c a b + a c, but
More informationReview Problems 3. Four FIR Filter Types
Review Prblems 3 Fur FIR Filter Types Fur types f FIR linear phase digital filters have cefficients h(n fr 0 n M. They are defined as fllws: Type I: h(n = h(m-n and M even. Type II: h(n = h(m-n and M dd.
More informationSurface and Contact Stress
Surface and Cntact Stress The cncept f the frce is fundamental t mechanics and many imprtant prblems can be cast in terms f frces nly, fr example the prblems cnsidered in Chapter. Hwever, mre sphisticated
More informationLecture 5: Equilibrium and Oscillations
Lecture 5: Equilibrium and Oscillatins Energy and Mtin Last time, we fund that fr a system with energy cnserved, v = ± E U m ( ) ( ) One result we see immediately is that there is n slutin fr velcity if
More informationOF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
U. S. FOREST SERVICE RESEARCH PAPER FPL 50 DECEMBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
More informationThermodynamics Partial Outline of Topics
Thermdynamics Partial Outline f Tpics I. The secnd law f thermdynamics addresses the issue f spntaneity and invlves a functin called entrpy (S): If a prcess is spntaneus, then Suniverse > 0 (2 nd Law!)
More informationCHAPTER 8b Static Equilibrium Units
CHAPTER 8b Static Equilibrium Units The Cnditins fr Equilibrium Slving Statics Prblems Stability and Balance Elasticity; Stress and Strain The Cnditins fr Equilibrium An bject with frces acting n it, but
More informationA Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture
Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu
More informationTrigonometric Ratios Unit 5 Tentative TEST date
1 U n i t 5 11U Date: Name: Trignmetric Ratis Unit 5 Tentative TEST date Big idea/learning Gals In this unit yu will extend yur knwledge f SOH CAH TOA t wrk with btuse and reflex angles. This extensin
More informationENSC Discrete Time Systems. Project Outline. Semester
ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding
More informationChE 471: LECTURE 4 Fall 2003
ChE 47: LECTURE 4 Fall 003 IDEL RECTORS One f the key gals f chemical reactin engineering is t quantify the relatinship between prductin rate, reactr size, reactin kinetics and selected perating cnditins.
More informationMaking and Experimenting with Voltaic Cells. I. Basic Concepts and Definitions (some ideas discussed in class are omitted here)
Making xperimenting with Vltaic Cells I. Basic Cncepts Definitins (sme ideas discussed in class are mitted here) A. Directin f electrn flw psitiveness f electrdes. If ne electrde is mre psitive than anther,
More informationENGI 4430 Parametric Vector Functions Page 2-01
ENGI 4430 Parametric Vectr Functins Page -01. Parametric Vectr Functins (cntinued) Any nn-zer vectr r can be decmpsed int its magnitude r and its directin: r rrˆ, where r r 0 Tangent Vectr: dx dy dz dr
More informationBuilding to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.
Building t Transfrmatins n Crdinate Axis Grade 5: Gemetry Graph pints n the crdinate plane t slve real-wrld and mathematical prblems. 5.G.1. Use a pair f perpendicular number lines, called axes, t define
More informationLead/Lag Compensator Frequency Domain Properties and Design Methods
Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin
More informationPreparation work for A2 Mathematics [2018]
Preparatin wrk fr A Mathematics [018] The wrk studied in Y1 will frm the fundatins n which will build upn in Year 13. It will nly be reviewed during Year 13, it will nt be retaught. This is t allw time
More informationPipetting 101 Developed by BSU CityLab
Discver the Micrbes Within: The Wlbachia Prject Pipetting 101 Develped by BSU CityLab Clr Cmparisns Pipetting Exercise #1 STUDENT OBJECTIVES Students will be able t: Chse the crrect size micrpipette fr
More information1 Course Notes in Introductory Physics Jeffrey Seguritan
Intrductin & Kinematics I Intrductin Quickie Cncepts Units SI is standard system f units used t measure physical quantities. Base units that we use: meter (m) is standard unit f length kilgram (kg) is
More informationAP Physics. Summer Assignment 2012 Date. Name. F m = = + What is due the first day of school? a. T. b. = ( )( ) =
P Physics Name Summer ssignment 0 Date I. The P curriculum is extensive!! This means we have t wrk at a fast pace. This summer hmewrk will allw us t start n new Physics subject matter immediately when
More informationLecture 7: Damped and Driven Oscillations
Lecture 7: Damped and Driven Oscillatins Last time, we fund fr underdamped scillatrs: βt x t = e A1 + A csω1t + i A1 A sinω1t A 1 and A are cmplex numbers, but ur answer must be real Implies that A 1 and
More informationALE 21. Gibbs Free Energy. At what temperature does the spontaneity of a reaction change?
Name Chem 163 Sectin: Team Number: ALE 21. Gibbs Free Energy (Reference: 20.3 Silberberg 5 th editin) At what temperature des the spntaneity f a reactin change? The Mdel: The Definitin f Free Energy S
More informationExaminer: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data
Benha University Cllege f Engineering at Banha Department f Mechanical Eng. First Year Mechanical Subject : Fluid Mechanics M111 Date:4/5/016 Questins Fr Final Crrective Examinatin Examiner: Dr. Mhamed
More informationThermodynamics and Equilibrium
Thermdynamics and Equilibrium Thermdynamics Thermdynamics is the study f the relatinship between heat and ther frms f energy in a chemical r physical prcess. We intrduced the thermdynamic prperty f enthalpy,
More informationKinematic transformation of mechanical behavior Neville Hogan
inematic transfrmatin f mechanical behavir Neville Hgan Generalized crdinates are fundamental If we assume that a linkage may accurately be described as a cllectin f linked rigid bdies, their generalized
More informationDead-beat controller design
J. Hetthéssy, A. Barta, R. Bars: Dead beat cntrller design Nvember, 4 Dead-beat cntrller design In sampled data cntrl systems the cntrller is realised by an intelligent device, typically by a PLC (Prgrammable
More informationWork, Energy, and Power
rk, Energy, and Pwer Physics 1 There are many different TYPES f Energy. Energy is expressed in JOULES (J 419J 4.19 1 calrie Energy can be expressed mre specifically by using the term ORK( rk The Scalar
More informationFebruary 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA
February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal
More informationCHM112 Lab Graphing with Excel Grading Rubric
Name CHM112 Lab Graphing with Excel Grading Rubric Criteria Pints pssible Pints earned Graphs crrectly pltted and adhere t all guidelines (including descriptive title, prperly frmatted axes, trendline
More informationLab 1 The Scientific Method
INTRODUCTION The fllwing labratry exercise is designed t give yu, the student, an pprtunity t explre unknwn systems, r universes, and hypthesize pssible rules which may gvern the behavir within them. Scientific
More informationSection 6-2: Simplex Method: Maximization with Problem Constraints of the Form ~
Sectin 6-2: Simplex Methd: Maximizatin with Prblem Cnstraints f the Frm ~ Nte: This methd was develped by Gerge B. Dantzig in 1947 while n assignment t the U.S. Department f the Air Frce. Definitin: Standard
More informationExperiment #3. Graphing with Excel
Experiment #3. Graphing with Excel Study the "Graphing with Excel" instructins that have been prvided. Additinal help with learning t use Excel can be fund n several web sites, including http://www.ncsu.edu/labwrite/res/gt/gt-
More informationPhysics 212. Lecture 12. Today's Concept: Magnetic Force on moving charges. Physics 212 Lecture 12, Slide 1
Physics 1 Lecture 1 Tday's Cncept: Magnetic Frce n mving charges F qv Physics 1 Lecture 1, Slide 1 Music Wh is the Artist? A) The Meters ) The Neville rthers C) Trmbne Shrty D) Michael Franti E) Radiatrs
More informationECEN 4872/5827 Lecture Notes
ECEN 4872/5827 Lecture Ntes Lecture #5 Objectives fr lecture #5: 1. Analysis f precisin current reference 2. Appraches fr evaluating tlerances 3. Temperature Cefficients evaluatin technique 4. Fundamentals
More informationMath 105: Review for Exam I - Solutions
1. Let f(x) = 3 + x + 5. Math 105: Review fr Exam I - Slutins (a) What is the natural dmain f f? [ 5, ), which means all reals greater than r equal t 5 (b) What is the range f f? [3, ), which means all
More information1 The limitations of Hartree Fock approximation
Chapter: Pst-Hartree Fck Methds - I The limitatins f Hartree Fck apprximatin The n electrn single determinant Hartree Fck wave functin is the variatinal best amng all pssible n electrn single determinants
More informationYou need to be able to define the following terms and answer basic questions about them:
CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f
More informationMore Tutorial at
Answer each questin in the space prvided; use back f page if extra space is needed. Answer questins s the grader can READILY understand yur wrk; nly wrk n the exam sheet will be cnsidered. Write answers,
More informationES201 - Examination 2 Winter Adams and Richards NAME BOX NUMBER
ES201 - Examinatin 2 Winter 2003-2004 Adams and Richards NAME BOX NUMBER Please Circle One : Richards (Perid 4) ES201-01 Adams (Perid 4) ES201-02 Adams (Perid 6) ES201-03 Prblem 1 ( 12 ) Prblem 2 ( 24
More informationMedium Scale Integrated (MSI) devices [Sections 2.9 and 2.10]
EECS 270, Winter 2017, Lecture 3 Page 1 f 6 Medium Scale Integrated (MSI) devices [Sectins 2.9 and 2.10] As we ve seen, it s smetimes nt reasnable t d all the design wrk at the gate-level smetimes we just
More informationPhysics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018
Michael Faraday lived in the Lndn area frm 1791 t 1867. He was 29 years ld when Hand Oersted, in 1820, accidentally discvered that electric current creates magnetic field. Thrugh empirical bservatin and
More informationExercise 3 Identification of parameters of the vibrating system with one degree of freedom
Exercise 3 Identificatin f parameters f the vibrating system with ne degree f freedm Gal T determine the value f the damping cefficient, the stiffness cefficient and the amplitude f the vibratin excitatin
More informationLecture 6: Phase Space and Damped Oscillations
Lecture 6: Phase Space and Damped Oscillatins Oscillatins in Multiple Dimensins The preius discussin was fine fr scillatin in a single dimensin In general, thugh, we want t deal with the situatin where:
More informationSemester 2 AP Chemistry Unit 12
Cmmn In Effect and Buffers PwerPint The cmmn in effect The shift in equilibrium caused by the additin f a cmpund having an in in cmmn with the disslved substance The presence f the excess ins frm the disslved
More information5 th grade Common Core Standards
5 th grade Cmmn Cre Standards In Grade 5, instructinal time shuld fcus n three critical areas: (1) develping fluency with additin and subtractin f fractins, and develping understanding f the multiplicatin
More informationI. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is
Length L>>a,b,c Phys 232 Lab 4 Ch 17 Electric Ptential Difference Materials: whitebards & pens, cmputers with VPythn, pwer supply & cables, multimeter, crkbard, thumbtacks, individual prbes and jined prbes,
More informationSPH3U1 Lesson 06 Kinematics
PROJECTILE MOTION LEARNING GOALS Students will: Describe the mtin f an bject thrwn at arbitrary angles thrugh the air. Describe the hrizntal and vertical mtins f a prjectile. Slve prjectile mtin prblems.
More informationDepartment of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets
Department f Ecnmics, University f alifrnia, Davis Ecn 200 Micr Thery Prfessr Giacm Bnann Insurance Markets nsider an individual wh has an initial wealth f. ith sme prbability p he faces a lss f x (0
More informationCHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India
CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce
More informationNGSS High School Physics Domain Model
NGSS High Schl Physics Dmain Mdel Mtin and Stability: Frces and Interactins HS-PS2-1: Students will be able t analyze data t supprt the claim that Newtn s secnd law f mtin describes the mathematical relatinship
More informationPhysics 101 Math Review. Solutions
Physics 0 Math eview Slutins . The fllwing are rdinary physics prblems. Place the answer in scientific ntatin when apprpriate and simplify the units (Scientific ntatin is used when it takes less time t
More information**DO NOT ONLY RELY ON THIS STUDY GUIDE!!!**
Tpics lists: UV-Vis Absrbance Spectrscpy Lab & ChemActivity 3-6 (nly thrugh 4) I. UV-Vis Absrbance Spectrscpy Lab Beer s law Relates cncentratin f a chemical species in a slutin and the absrbance f that
More informationFive Whys How To Do It Better
Five Whys Definitin. As explained in the previus article, we define rt cause as simply the uncvering f hw the current prblem came int being. Fr a simple causal chain, it is the entire chain. Fr a cmplex
More informationCHAPTER 2 Algebraic Expressions and Fundamental Operations
CHAPTER Algebraic Expressins and Fundamental Operatins OBJECTIVES: 1. Algebraic Expressins. Terms. Degree. Gruping 5. Additin 6. Subtractin 7. Multiplicatin 8. Divisin Algebraic Expressin An algebraic
More informationAdmin. MDP Search Trees. Optimal Quantities. Reinforcement Learning
Admin Reinfrcement Learning Cntent adapted frm Berkeley CS188 MDP Search Trees Each MDP state prjects an expectimax-like search tree Optimal Quantities The value (utility) f a state s: V*(s) = expected
More informationTHERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES
PREFERRED RELIABILITY PAGE 1 OF 5 PRACTICES PRACTICE NO. PT-TE-1409 THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC Practice: Perfrm all thermal envirnmental tests n electrnic spaceflight hardware in a flight-like
More information3. Design of Channels General Definition of some terms CHAPTER THREE
CHAPTER THREE. Design f Channels.. General The success f the irrigatin system depends n the design f the netwrk f canals. The canals may be excavated thrugh the difference types f sils such as alluvial
More informationA Correlation of. to the. South Carolina Academic Standards for Mathematics Precalculus
A Crrelatin f Suth Carlina Academic Standards fr Mathematics Precalculus INTRODUCTION This dcument demnstrates hw Precalculus (Blitzer), 4 th Editin 010, meets the indicatrs f the. Crrelatin page references
More informationA - LEVEL MATHEMATICS 2018/2019
A - LEVEL MATHEMATICS 2018/2019 STRUCTURE OF THE COURSE Yur maths A-Level Maths curse cvers Pure Mathematics, Mechanics and Statistics. Yu will be eamined at the end f the tw-year curse. The assessment
More informationRigid Body Dynamics (continued)
Last time: Rigid dy Dynamics (cntinued) Discussin f pint mass, rigid bdy as useful abstractins f reality Many-particle apprach t rigid bdy mdeling: Newtn s Secnd Law, Euler s Law Cntinuus bdy apprach t
More informationLyapunov Stability Stability of Equilibrium Points
Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),
More informationLim f (x) e. Find the largest possible domain and its discontinuity points. Why is it discontinuous at those points (if any)?
THESE ARE SAMPLE QUESTIONS FOR EACH OF THE STUDENT LEARNING OUTCOMES (SLO) SET FOR THIS COURSE. SLO 1: Understand and use the cncept f the limit f a functin i. Use prperties f limits and ther techniques,
More informationPHYS 314 HOMEWORK #3
PHYS 34 HOMEWORK #3 Due : 8 Feb. 07. A unifrm chain f mass M, lenth L and density λ (measured in k/m) hans s that its bttm link is just tuchin a scale. The chain is drpped frm rest nt the scale. What des
More informationx 1 Outline IAML: Logistic Regression Decision Boundaries Example Data
Outline IAML: Lgistic Regressin Charles Suttn and Victr Lavrenk Schl f Infrmatics Semester Lgistic functin Lgistic regressin Learning lgistic regressin Optimizatin The pwer f nn-linear basis functins Least-squares
More informationmaking triangle (ie same reference angle) ). This is a standard form that will allow us all to have the X= y=
Intrductin t Vectrs I 21 Intrductin t Vectrs I 22 I. Determine the hrizntal and vertical cmpnents f the resultant vectr by cunting n the grid. X= y= J. Draw a mangle with hrizntal and vertical cmpnents
More informationSection I5: Feedback in Operational Amplifiers
Sectin I5: eedback in Operatinal mplifiers s discussed earlier, practical p-amps hae a high gain under dc (zer frequency) cnditins and the gain decreases as frequency increases. This frequency dependence
More informationand the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:
Algrithm fr Estimating R and R - (David Sandwell, SIO, August 4, 2006) Azimith cmpressin invlves the alignment f successive eches t be fcused n a pint target Let s be the slw time alng the satellite track
More informationLHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers
LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the
More information37 Maxwell s Equations
37 Maxwell s quatins In this chapter, the plan is t summarize much f what we knw abut electricity and magnetism in a manner similar t the way in which James Clerk Maxwell summarized what was knwn abut
More informationENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS
ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS On cmpletin f this tutrial yu shuld be able t d the fllwing. Define viscsity
More informationDepartment of Electrical Engineering, University of Waterloo. Introduction
Sectin 4: Sequential Circuits Majr Tpics Types f sequential circuits Flip-flps Analysis f clcked sequential circuits Mre and Mealy machines Design f clcked sequential circuits State transitin design methd
More information