Solving Differential Equations by the Laplace Transform and by Numerical Methods
|
|
- Harold Waters
- 6 years ago
- Views:
Transcription
1 36CH_PHCalter_TechMath_ //007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the Laplace tranform of a function by direct integration or by uing a table. Find the Laplace tranform of an expreion containing derivative, and ubtitute initial condition. Determine the invere of a Laplace tranform by completing the quare or by partial fraction. Solve firt-order and econd-order differential equation uing Laplace tranform. Solve electrical application uing Laplace tranform. Solve firt-order differential equation uing Euler method, the modified Euler method, or the Runge-Kutta method. Solve econd-order differential equation uing the modified Euler method or the Runge-Kutta method. The method for olving differential equation that we learned in Chapter 8 and 9 are often called claical method. In thi chapter we learn other powerful technique: the Laplace tranform and numerical method. The Laplace tranform enable u to tranform a differential equation into an algebraic one. We can then olve the algebraic equation and apply an invere tranform to obtain the olution to the differential equation. The Laplace tranform i good for finding particular olution of differential equation. For example, it enable u to olve initial value problem, that i, when the value of the function i known at t 0. Thu the equation that we deal with here are function of time, uch a electric circuit problem, rather than function of x. We firt ue the Laplace tranform to olve firt- and econd-order differential equation with contant coefficient. We then do ome application and compare our reult with thoe obtained by claical method in Chapter 8 and 9. We go on to olve differential equation uing numerical method. Thee method ue
2 36CH_PHCalter_TechMath_ //007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method ucceive approximation to give a olution in the form of number, rather than equation. One great advantage of numerical method i that they can be programmed on a computer. The Laplace Tranform of a Function Laplace Tranform of Some Simple Function The Laplace tranform i named after the French analyt, probabilit, atronomer, and phyicit Pierre Laplace (749 87). If y i ome function of time, o that y f(t), the Laplace tranform of that function i defined by the following improper integral: Laplace Tranform [f(t)] 0 f(t)et dt The tranformed expreion i a function of, which we call F(). Thu [ f(t)] F() We can write the tranform of an expreion by direct integration, that i, by writing an integral uing Eq. and then evaluating it. Glance back at Sec. 8 in Method of Integration if you have forgotten how to evaluate an improper integral. y Example : If y f(t), then [ f(t)] [] 0 ()et dt To evaluate the integral, we multiply by and compenate with /. y = e u [] 0 et ( dt) A et 0 Figure how a graph of y e u, which will help u to evaluate the limit. Note that a u approache infinity, e u approache zero, and a u approache zero, e u approache one. So we take e t equal to zero at the upper limit t and e t at the lower limit t 0. Thu 0 u [] (0 ) FIGURE Rule from our table of integral in Text Appendix C ay that af(x) dx a f(x) dx. It follow then that [af(t)] a[ f (t)] The Laplace tranform of a contant time a function i equal to the contant time the tranform of the function. Example : If [] /, a above, then [5] 5[] 5 Example 3: If y f(t) t, then [ f(t)] [t] 0 tet dt
3 36CH_PHCalter_TechMath_ //007 :8 PM Page 3 Section The Laplace Tranform of a Function 3 We evaluate the right ide uing Text Rule 37 with u t and a. [t] e t (t ) A 0 At the upper limit, e t approache zero, while t approache infinity. But e t i hrinking fater than t i growing, o their product approache zero. Thi i confirmed by Fig., which how the product t(e t ) firt increaing but then quickly approaching zero. So we get [t] e t )A (t 0 0 y y = te t t FIGURE Example 4: Find [in at]. Solution: Uing Eq., we find that [in at] 0 et in at dt Uing Text Rule 4 give t e [in at] ( in at a co at)a a 0 At the upper limit of infinity, e t approache zero, cauing the entire expreion to have a value of zero at that limit. At the lower limit of zero, e t approache one, in at approache zero, and co at approache one. Our expreion then become a [in at] 0 (a) a a Note that in each cae the tranformed expreion i a function of only. Tranform of a Sum If we have the um of everal term, f(t) a g(t) b h(t) then [ f(t)] [a g(t) b h(t) 0 ]et dt a 0 g(t)et dt b 0 h(t)et dt Thu [a g(t) b h(t) ] a[g(t)] b[h(t)]
4 36CH_PHCalter_TechMath_ //007 :8 PM Page 4 4 Solving Differential Equation by the Laplace Tranform and by Numerical Method The Laplace tranform of a um of everal term i the um of the tranform of each term. Thi property allow u to work term by term in writing the tranform of an expreion made up of um or difference. Example 5: If y 4 5t in 3t, then [ y] When we write the Laplace tranform of a function, we often ay that we are taking the tranform of the function. Tranform of a Derivative To olve differential equation, we mut be able to take the tranform of a derivative. Let f (t) be the derivative of ome function f (t). Then [ f (t)] f 0 (t)et dt Integrating by part, we let u e t and d f(t) dt. Sodu e t dt and f (t) dt f(t). Then [ f(t)] f(t)e A t f 0 0 (t)et dt f(t)e t A 0 [ f(t)] Thu: 0 f(0) [ f(t)] Tranform of a Derivative [ f (t)] [ f(t)] f(0) The tranform of a derivative of a function equal time the tranform of the function minu the function evaluated at t 0. Thu if we have a function y f (t), then [ y] [ y] y(0) Note that the tranform of a derivative contain f(0) or y(0), the value of the function at t 0. Thu to evaluate the Laplace tranform of a derivative, we mut know the initial condition. If two function are equal, their tranform are equal. Thu we can take the tranform of both ide of an equation without changing the meaning of the equation. Example 6: Take the Laplace tranform of both ide of the differential equation y 6t, with the initial condition that y(0) 5. Solution: y 6t [y] [6t] 6 [ y] y(0) or, ubtituting, [ y] 5 6
5 36CH_PHCalter_TechMath_ //007 :8 PM Page 5 Section The Laplace Tranform of a Function 5 ince y(0) 5. Later we will olve uch an equation for [y] and then find the invere of the tranform to obtain y. Similarly, we can find the tranform of the econd derivative ince the econd derivative of f(t) i the derivative of f (t). [ f (t)] [ f (t)] f (0) Subtituting, we get [ f (t)] {[ f(t)] f (0)} f (0) Tranform of the Second Derivative [ f (t)] [ f(t)] f(0) f (0) Another way of expreing thi equation, if y f(t), i [ y] [ y] y(0) y(0) Example 7: Tranform both ide of the econd-order differential equation y 3y 4y 5t if y(0) 6 and y(0) 7. Solution: y 3y 4y 5t [ y] 3[ y] 4[ y] [5t] 5 [y] y(0) y(0) 3{[y] y(0)} 4[y] Subtituting y(0) 6 and y(0) 7, we obtain [ y] 6 7 3[ y] 3(6) 4[ y] 5 or [y] 3[y] 4[y] 6 5 Tranform of an Integral Let u now find the tranform of the integral t f(t) dt. By our definition of the 0 Laplace tranform, Eq., t f(t) dt 0 0 t f(t) dt 0 et dt Integrating by part, we let o that Integrating yield u t f(t) dt and d 0 et dt du f(t) dt and e t / t f(t) dt 0 t (/)et f(t) dt A (/) f(t)et dt At the upper limit of infinity, the quantity in the bracket vanihe becaue e t goe to zero. At the lower limit of zero, the quantity in the bracket vanihe becaue the
6 36CH_PHCalter_TechMath_ //007 :8 PM Page 6 6 Solving Differential Equation by the Laplace Tranform and by Numerical Method upper and lower limit of the integral are equal. Further, the integral in the econd term i the definition (Eq. ) of the Laplace tranform of f(t). Thu: Tranform of an Integral t f(t) dt 0 (/) f(t) We will need thi equation later for olving electrical problem. Example 8: The voltage acro a capacitor of capacitance C i given by (C) t i dt 0 Find the tranform of thi voltage. Solution: Uing the tranform of an integral, we get (C) t i dt 0 (C)[i] Table of Laplace Tranform The tranform of ome common function are given in Table. Intead of tranforming a function tep by tep, we imply look it up in the table. TABLE Short table of Laplace tranform. Here n i a poitive integer. Tranform Number f(t) [ f (t)] F() f(t) F() e t f(t) dt f (t) [ f(t)] f(0) 3 f (t) [ f(t)] f(0) f(0) 4 a g(t) b h(t) a[g(t)] b[h(t)] 5 6 t 7 t n n! n t 8 (n )! n 9 e at a 0 e at a ( a) te at ( a) e at ( at) ( a) n 3 t n e at n! ( a) n 4 i n e at ( n )! n ( a) 5 e at e bt b a ( a)( b) 0
7 36CH_PHCalter_TechMath_ //007 :8 PM Page 7 Section The Laplace Tranform of a Function 7 TABLE Continued Tranform Number f(t) [ f (t)] F() 6 ae at be bt (a b) ( a) ( b) 7 in at a a 8 co at a 9 t in at a ( a ) a 0 t co at ( a ) a co at ( a ) 3 a at in at ( a ) 3 e at b in bt ( a ) b 4 e at a co bt ( a) b a 3 5 in at at co at ( a ) 6 in at at co at a ( a ) 7 co at at in at 3 ( a ) b 8 a (e at b at ) ( a) 9 t f(t) dt [f (t)] 0 Example 9: Find the Laplace tranform of t 3 e t by uing Table. Solution: Our function matche Tranform 3, with n 3 and a, o [t 3 e t 3! 6 ] ( ) 3 ( ) 4 Exercie The Laplace Tranform of a Function Tranform by Direct Integration Find the Laplace tranform of each function by direct integration.. f(t) 6. f(t) t 3. f(t) t 4. f(t) t 5. f(t) co 5t 6. f(t) e t in t Tranform by Table Ue Table to find the Laplace tranform of each function. 7. f(t) t 4 8. f(t) t 3 t 3t 4
8 36CH_PHCalter_TechMath_ //007 :8 PM Page 8 8 Solving Differential Equation by the Laplace Tranform and by Numerical Method 9. f(t) 3e t e t 0. f(t) 5te 3t. f(t) in t co 3t. f(t) 3 e 4t 3. f(t) 5e 3t co 5t 4. f(t) e 4t t 5. f(t) 5e t t in 3t 6. f(t) t 3 4t 3e t Tranform of Derivative Find the Laplace tranform of each expreion, and ubtitute the given initial condition. 7. y y, y(0) 8. 3y y, y(0) 3 9. y 4y, y(0) y 3y, y(0). y 3y y, y(0) and y(0) 3. y y y, y(0) and y(0) 0 3. y 3y y, y(0) and y(0) y y y, y(0) and y(0) Invere Tranform Before we can ue the Laplace tranform to olve differential equation, we mut be able to tranform a function of back to a function of t. The invere Laplace tranform i denoted by. Thu if [ f(t)] F(), then we have the following equation: We ll ee that finding the invere i often harder than finding the tranform itelf. Invere Laplace Tranform [F()] f(t) We ue Table to find the invere of ome Laplace tranform. We put our given expreion in a form that matche a right-hand entry of Table and then read the correponding entry at the left. Example 0: If F() 4/( 6), find f(t). Solution: We earch Table in the right-hand column for an expreion of imilar form and find Tranform 7, a [in at] a which matche our expreion if a 4. Thu f(t) in 4t Example : Find y if [y] 5/( 7) 4. Solution: From the table we find Tranform 3, [t n e at n! ] ( a) n
9 36CH_PHCalter_TechMath_ //007 :8 PM Page 9 Section Invere Tranform 9 In order to match our function, we mut have a 7 and n 3. The numerator then mut be 3! 3() 6. So we inert 6 in the numerator and compenate with 6 in the denominator. We then rewrite our function a [ y] ( 7) 4 which now i the ame a Tranform 3 with a coefficient of 5 6. The invere i then y 5 6 t3 e 7t Completing the Square To get our given function to match a right-hand table entry, we may have to do ome algebra. Sometime we mut complete the quare to make our function match one in the table. Example : Find y if [ y]. 4 0 Solution: Thi doe not now match any function in our table, but it will if we complete the quare on the denominator. 4 0 ( 4 ) 0 ( 4 4) 0 4 ( ) 4 The denominator i now of the ame form a in Tranform 4. However, to ue Tranform 4, the numerator mut be, not. So let u add 3 and ubtract 3 from the numerator, o that ( ) 3. Then ( ) 3 [y] 4 0 ( ) 4 ( ) 4 3 ( ) 4 3 ( ) ( ) 4 Now the firt expreion matche Tranform 4 and the econd matche Tranform 3. We find the invere of thee tranform to be Thi i no different from the method we ued earlier when we completed the quare. y e t co 4t 3 4 et in 4t e t co 4t 3 4 in 4t Partial Fraction We ued partial fraction in Sec. 5 in Method of Integration to make a given expreion match one lited in the table of integral. Now we ue partial fraction to make an expreion match one lited in our table of Laplace tranform. Example 3: Find y if [y] /( 8). Solution: We eparate [y] into partial fraction. [y] 8 A B 4 o A( ) B( 4) (A B) (A 4B) We could alo complete the quare here, but we would find that the reulting expreion would not match a table entry. Some trial-and-error work in algebra may be neceary to get the function to match a table entry.
10 36CH_PHCalter_TechMath_ //007 :8 PM Page 0 0 Solving Differential Equation by the Laplace Tranform and by Numerical Method from which A B 0 and A 4B. Solving imultaneouly give A and B, o [ y] 8 4 Uing Tranform 9 twice, we get y e 4t e t Exercie Invere Tranform Find the invere of each tranform ( 6) ( 9) ( 4) 5. ( ) ( ) ( )( ) ( ) ( ) ( 4)( ) 3 Solving Differential Equation by the Laplace Tranform To olve a differential equation with the Laplace tranform:. Take the tranform of each ide of the equation.. Solve for [y] F(). 3. Manipulate F() until it matche one or more table entrie. 4. Take the invere tranform to find y f(t). We tart with a very imple example. Example 4: Solve the firt-order differential equation y y if y(0) 0.
11 36CH_PHCalter_TechMath_ //007 :8 PM Page Section 3 Solving Differential Equation by the Laplace Tranform Solution:. We write the Laplace tranform of both ide and get [y] [y] [] [y] y(0) [y]. We ubtitute 0 for y(0) and olve for [y]. ( )[y] [y] ( ) 3. Our function will match Tranform 0 if we write it a [y] ( ) 4. Taking the invere give y ( e t ) In our next example we expre F() a three partial fraction in order to take the invere tranform. Example 5: Solve the firt-order differential equation y 3y 4 9t if y(0). Solution:. We take the Laplace tranform of both ide and get [y] [3y] [4] [9t] [y] y(0) 3[y] 4 9. We ubtitute for y(0) and olve for [y]. 9 ( 3)[y] [y] ( 3) ( 3) 3 3. We match our expreion with Tranform 8, 0, and 9. [y] ( 4. Taking the invere give 9 4 3) 3 3 ( 3) 3 y e 3t 3t 4 3 ( e3t ) e 3t 5 3 e3t 3t 3 We now ue the Laplace tranform to olve a econd-order differential equation. Example 6: Solve the econd-order differential equation y 4y 3 0 where y i a function of t and y and y are both zero at t 0.
12 36CH_PHCalter_TechMath_ //007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method Solution:. Taking the Laplace tranform of both ide, [y] [4y] [3] 0 [y] y(0) y(0) 4[y] 3 0 Subtituting y(0) 0 and y(0) 0 give. Solving for [y], we have 3. We match it to Tranform. 4. Taking the invere give [y] 0 0 4[y] 3 0 ( 4)[y] 3 [y] ( [y] ( 3 3 4) 4 3 4) 4 ( ) y 3 4 ( co t) We have een that the hardet part in olving a differential equation by Laplace tranform i often in doing the algebra neceary to get the given function to match a table entry. We mut often ue partial fraction or completing the quare, and ometime both, a in the following example. Example 7: Solve y 4y 5y t, where y(0) and y(0). Solution:. Taking the Laplace tranform of both ide and ubtituting initial condition, we find that [y] y(0) y(0) 4[ y] 4y(0) 5[y]. Solving for [y] give [ y]( 4 5) 6 6 [y] () ( 4 5) Taking for now jut the firt fraction in Equation (), which we ll call F (), we obtain F () A B ( 4 5) C D 4 5 Multiplying by ( 4 5) give A 3 4A 5A B 4B 5B C 3 D (A C) 3 (4A B D) (5A 4B) 5B
13 36CH_PHCalter_TechMath_ //007 :8 PM Page 3 Section 3 Solving Differential Equation by the Laplace Tranform 3 Equating coefficient give A 4 B 5 5 C 4 D 5 5 So F () 4/ 5 /5 4 /5 / We can ue Tranform 5 and 6 for 4/ and 5/, but there i no match for the third term in the parenthee. However, by completing the quare, we can get the denominator to match thoe in Tranform 3 and 4. Thu ( ) We now manipulate that third term into the form of Tranform 3 and (4) 4 5 ( ) / 4 ( ) (4) ( 3/ 4 ) (4) ( ) (3) ( ) Combining the third term with the firt and econd give 5 F () 5 4 (4) ( ) (3) ( ) Returning to the remaining fraction in Equation (), which we call F (), we complete the quare in the denominator and ue partial fraction to get 6 F () ( ) 6 ( ) ( ) 4 ( ) 4. Taking the invere tranform of F () and F (), we obtain y [5t 4 4e t co t 3e t in t] e t co t 4e t in t 5 which implifie to y (5t 4 9e t co t 03e t in t) 5 Exercie 3 Solving Differential Equation by the Laplace Tranform Firt-Order Equation Solve each differential equation by the Laplace tranform.. y 3y 0, y(0). y y, y(0) y y t, y(0) 0 4. y 5y e t, y(0) 5. 3y y t, y(0) 3 6. y 3y in t, y(0) 0 7. y y co t, y(0) y y 3t 3, y(0) 0
14 36CH_PHCalter_TechMath_ //007 :8 PM Page 4 4 Solving Differential Equation by the Laplace Tranform and by Numerical Method Second-Order Equation Solve each differential equation by the Laplace tranform. 9. y y 3y 0, y(0) 0 and y(0) 0. y y y, y(0) 0 and y(0) 0. y 3y 3, y(0) and y(0). y y, y(0) 0 and y(0) 3 3. y y 4t, y(0) 3 and y(0) 0 4. y y 5y t, y(0) and y(0) 5. y 4y 3y t, y(0) and y(0) 6. y 4y 8t 3, y(0) 0 and y(0) y y in t, y(0) and y(0) 3 8. y y 3, y(0) and y(0) 9. y y y e t, y(0) 0 and y(0) 0 0. y 3y co t, y(0) 0 and y(0). y y 3y te t, y(0) 0 and y(0) 0. 3y y y in 3t, y(0) 0 and y(0) 0 Motion in a Reiting Medium 3. A 5.0-kg object i dropped from ret through air whoe reiting force i equal to.85 time the object peed (in m/). Find it peed after.00. Exponential Growth and Decay 4. A certain teel ingot i 900 F and cool at a rate (in F/min) equal to.5 time it preent temperature ( F). Find it temperature after 5.00 min. 5. The rate of growth (bacteria/h) of a colony of bacteria i equal to.5 time the number preent at any intant. How many bacteria are there after 4 h if there are 5000 at firt? Mechanical Vibration 6. A weight that hang from a pring i pulled down.00 in. at t 0 and then i releaed from ret. The differential equation of motion i x 6.5x 5.8x 0. Write the equation for x a a function of time. 7. An alternating force i applied to a weight uch that the equation of motion i x 6.5x 45.3 co.5t. If and x are zero at t 0, write an equation for x a a function of time. 4 Electrical Application R Many of the differential equation for motion and electric circuit covered in Chapter 8 and 9 are nicely handled by the Laplace tranform. However, ince the Laplace tranform i ued mainly for electric circuit, we emphaize that application here. The method i illutrated by example of everal type of circuit with dc and ac ource. E i C Serie RC Circuit with dc Source FIGURE 3 RC circuit. Example 8: Capacitor Dicharging. A fully charged capacitor, Fig. 3, i dicharged by throwing the witch from poition to poition, at t 0. Write expreion for the current and the voltage acro the capacitor.
15 36CH_PHCalter_TechMath_ //007 :8 PM Page 5 Section 4 Electrical Application 5 Solution: If the voltage acro the capacitor i, then the voltage acro the reitor mut be, ince the um of the voltage around the loop mut be zero. Since the current through the reitor (/R) mut equal the current through the capacitor (C d/dt) we write or Taking the tranform of each term or Solving for [i] we get R C d dt 0 RC [i] 0 [i] 0 RC [i] [i] E RC E [i] /RC Uing tranform 9 to find the invere give Voltage acro a Dicharging Capacitor t /RC Ee 3 We next find the current by dividing the voltage by R. Current in a Dicharging Capacitor i E R e t/rc Example 9: Capacitor Charging. A fully dicharched capacitor, Fig. 4, i charged by throwing the witch from poition to poition at t 0. Write expreion for the current and voltage acro the capacitor. 4 Note that the reult of thi example and the two example to follow are the ame a we obtained by claical method for olving differential equation. Solution: Summing voltage around the loop give Drop acro reitor Drop acro capacitor Battery voltage () By Text Eq. 08, the voltage acro the capacitor i C t i dt, o or Ri C t i dt E 0 0 E R i E RC i t i dt EC 0 FIGURE 4
16 36CH_PHCalter_TechMath_ //007 :8 PM Page 6 6 Solving Differential Equation by the Laplace Tranform and by Numerical Method Taking the tranform of each term RC[i] [i] E C Solving for [i] and rearranging give [i] E R /RC Uing tranform 9 to find the invere tranform, we get Current in a Charging Capacitor i E R e trc 5 Returing to () to get the voltage acro the capacitor, Drop acro capacitor = Battery voltage Drop acro reitor E R E e t RC R or Voltage acro a Charging Capacitor E e t/rc 6 Serie RL Circuit with dc Source Example 0: Inductor dicharging. A fully charged inductor i dicharged by throwing the witch, Fig. 5, from poition to poition at t = 0. Write expreion for the current and the voltage acro the inductor. Ω Solution: Summing voltage around the loop we get 4 V i 4 H Battery voltage Drop acro reitor Drop acro inductor () E Ri L d i dt FIGURE 5 or E L R i i L Taking the tranform of each term E L R [i] [i] i(0) L Subtituting 0 for i(0) and olving for [i], [i] E R R/L ( R/L)
17 36CH_PHCalter_TechMath_ //007 :8 PM Page 7 Section 4 Electrical Application 7 Then by tranform 0, Current in a Charging Inductor i E R e Rt /L 7 Now uing () to find the voltage acro the inductor, Voltage acro inductor Battery voltage Voltage acro the reitor E Ri E R E R e Rt /L Voltage acro a Charging Inductor Ee Rt L 8 Serie RL Circuit with ac Source Example : A witch (Fig. 6) i cloed at t 0 when the applied voltage i zero and increaing. Find the current. Solution: Summing the voltage around the loop give Ri L d i dt Subtituting value and rearranging, we find that 4i 0.0 d i 80 in 400t dt Taking the tranform of each term yield 80(400) 4[i] 0.0[[i] i(0)] (400 ) Subtituting 0 for i(0) and olving for [i], we have [i] [i](0.0 4) 3,000 [ (400) ](0.0 4) 3, 000 ( 400) We now ue the method of Sec. 5 in Method of Integration for eparating into partial fraction a rational fraction that ha quadratic factor in the denominator.,600,000 [ (400) ] ( 00),600,000 [ (400) ]( 00) A B C (400) 00 4 Ω 80 in 400t V 0.0 H FIGURE 6 RL circuit with ac ource.
18 36CH_PHCalter_TechMath_ //007 :8 PM Page 8 8 Solving Differential Equation by the Laplace Tranform and by Numerical Method Multiplying by [ (400) ]( 00) give,600,000 (A B)( 00) C[ (400) ] which implifie to,600,000 (A C) (00A B) 00B 400 C Equating coefficient give A C 0 00A B 0 00B 400C,600,000 Solving imultaneouly for A, B, and C give A 8 B 600 C 8 Subtituting back, we get [i] ( 400) (400) 8 (400) 8 00 Taking the invere tranform uing Tranform 9, 7, and 8, we find that i 4 in 400t 8 co 400t 8e 00t We now combine the firt two term in the form I in(400t v), where I and v tan rad o i 8.94 in(400t.07) 8e 00t A A graph of thi wave, a well a the applied voltage wave, i hown in Fig. 7. i (A) 0 5 e (V) 00 i e t (m) Ω FIGURE V.5 H 4.75 F Serie RLC Circuit with dc Source FIGURE 8 RLC circuit. Example : A witch (Fig. 8) i cloed at t 0, and there i no initial charge on the capacitor. Write an expreion for the current.
19 36CH_PHCalter_TechMath_ //007 :8 PM Page 9 Section 4 Electrical Application 9 Solution:. The differential equation for thi circuit i Ri L d i dt C t i dt E 0. Tranforming each term and ubtituting value give 5[i].5{ [i] i(0)} 6 0 [ i] Setting i(0) equal to zero and olving for [i] yield [i] , We now try to match thi expreion with thoe in the table. Let u complete the quare in the denominator ,000 ( ) 40, ( 75) (34,000) ( 75) j (367) ( 75 j367)( 75 j367) after replacing by j and factoring the difference of two quare. So 50.3 [i] 50 40, ( 75 j367)( 75 j367) We can now find partial fraction ( 75 j367)( 75 j367) A 75 j367 Multiplying through by 75 j367 give 50.3 A B( 75 j367) 7 5 j j367 Setting 75 j367, we get A j j j367 j734 We find B by multiplying through by 75 j367 and then etting j36775, and we get B j0.065 (work not hown). Subtituting give u j j [i] 75 j j367 Thi till doen t match any table entry. Let u now combine the two fraction over a common denominator. j j [i] 75 j j367 j0.0685( 75 j367) j0.0685( 75 j367) ( 75 j367)( 75 j367) j ,000 ( 75) (367) B 75 j367 We ee that mot of the work when uing the Laplace tranform i in algebra, rather than in calculu.
20 36CH_PHCalter_TechMath_ //007 :8 PM Page 0 0 Solving Differential Equation by the Laplace Tranform and by Numerical Method 5. We finally have a match, with Tranform 3. Taking the invere tranform give i 37e 75t in 367t ma which agree with the reult found in Chapter 9, Example 9. Exercie 4 Electrical Application. The current in a certain RC circuit atifie the equation 7i 750i 5. If i i zero at t 0, how that i 4.8 e 6t ma. In an RL circuit, R 3750 and L 0.50 H. It i connected to a dc ource of 50 V at t 0. If i i zero at t 0, how that i 66.7 e 5,000t ma 3. The current in an RLC circuit atifie the equation 0.55i 48i 7350i 0 If i and i are zero at t 0, how that i 6.4 e 6t ma 4. In an RLC circuit, R 950, L 0.0 H, and C.55 mf. It i connected to a dc ource of 5 V at t 0. If i i zero at t 0, how that i 39.7e t ma 5. The current in a certain RLC circuit atifie the equation 3.5i 0.03i 680 t i dt If i i zero at t 0, how that i 69e 70.5t in 65t ma 6. In an RL circuit, R 450 and L 0.7 H. It i connected to a dc ource of 8.4 V at t 0. If i i zero at t 0, how that i 6.84 e 3,700t ma 7. In an RL circuit, R 8.37 and L 0.50 H. It i connected to a dc ource of 50.5 V at t 0. If i i zero at t 0, how that i 6.03 e 33.5t A 5 Numerical Solution of Firt-Order Differential Equation In Sec. 8, we howed how to ue Euler method to graphically and numerically olve a differential equation, but we mentioned there that it i not the mot accurate method available. It will, however, erve a a good introduction to other method. We will now how both the modified Euler method that include predictor-corrector tep and the Runge-Kutta method. We will apply thee method to both firt- and econd-order equation. Modified Euler Method with Predictor-Corrector Step The variou predictor-corrector method all have two main tep. Firt, the predictor tep trie to predict the next point on the curve from preceding value, and then the corrector tep trie to improve the predicted value. 0
21 36CH_PHCalter_TechMath_ //007 :8 PM Page Section 5 Numerical Solution of Firt-Order Differential Equation Here we modify Euler method of Sec. 8 a the predictor, uing the coordinate and lope at the preent point P to etimate the next point Q. The corrector tep then recompute Q uing the average of the lope at P and Q. Our iteration formula i then the following: Modified Euler Method x q x p x y q y p m p m q x 9 Example 3: Find an approximate olution to y x /y, with the boundary condition that y when x 3. Calculate y for x 3 to 0 in tep of. Solution: Uing data from before, the lope m at (3, ) wa 4.5, the econd point wa predicted to be (4, 6.5), and the lope m at (4, 6.5) wa.46. The recomputed ordinate at the econd point i then, by Eq. 9, y () 5.48 o our corrected econd point i (4, 5.48). The lope at thi point i We ue thi to predict y 3. m y(4, 5.48).99 y 3 (predicted) y m x () The lope at (5, 8.400) i m 3 y(5, 8.400).976 The corrected y 3 i then y 3 (corrected) () 8.49 The remaining value are given in Table. TABLE x Approximate y Exact y Error Note that the error at x 0 i about 0.6%, compared with 0.6% before. We have gained in accuracy at the cot of added computation.
22 36CH_PHCalter_TechMath_ //007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method Thi method i named after two German mathematician, Carl Runge (856 97) and Wilhelm Kutta ( ). Runge-Kutta Method The Runge-Kutta method, like the modified Euler, i a predictor-corrector method. Prediction: Starting at P (Fig. 9), we ue the lope m p to predict R, the lope m r to predict S, and the lope m to predict Q. y y q Q m y S m r y p 0 P x p x x r R m p x x q x x FIGURE 9 Runge-Kutta method. Correction: We then take a weighted average of the lope at P, Q, R, and S, giving twice the weight to the lope at the midpoint R and S a at the endpoint P and Q. m avg 6 (m p m r m m q ) We ue the average lope and the coordinate of P to find Q. y q y p m avg x x q x p x y q y p m avg x Runge- Kutta Method where m avg 6 (m p m r m m q ) and m p f(x p, y p ) m r f x p x, y p m p x m f x p x, y p m r x m q f(x p x, y p m x) 0 Example 4: Repeat Example 3 uing the Runge-Kutta method and a tep of.
23 36CH_PHCalter_TechMath_ //007 :8 PM Page 3 Section 5 Numerical Solution of Firt-Order Differential Equation 3 Solution: Starting at P(3, ), with m p 4.5, y r 4.5(0.5) 4.5 m r f(3.5, 4.5) ( ) y.88(0.5) 3.44 m f(3.5, 3.44) ( 3.5) y q 3.560() m q f(4, 5.560) (.88) (3.560).878 m avg y q (corrected) y p m avg x 3.377() So our econd point i (3, 5.377). The ret of the calculation i given in Table 3. TABLE 3 x Approximate y Exact y Error The error here at x 0 with the Runge-Kutta method i about 0.0%, compared with 0.6% for the modified Euler method. Comparion of the Three Method Table 4 how the error obtained when olving the ame differential equation by our three different method, each with a tep ize of. TABLE 4 x Euler Modified Euler Runge-Kutta Divergence A with mot numerical method, the computation can diverge, with our anwer getting further and further from the true value. The computation hown here can
24 36CH_PHCalter_TechMath_ //007 :8 PM Page 4 4 Solving Differential Equation by the Laplace Tranform and by Numerical Method diverge if either the tep ize x i too large or the range of x over which we compute i too large. Numerical method can alo diverge if the point ued in the computation are near a dicontinuity, a maximum or minimum point, or a point of inflection. When doing the computation on the computer (the only practical way), it i eay to change the tep ize. A good practice i to reduce the tep ize by a factor of 0 each time (, 0., 0.0,...) until the anwer do not change in a given decimal place. Further reduction of the tep ize will eventually reult in number too mall for the computer to handle, and unreliable reult will occur. Exercie 5 Numerical Solution of Firt-Order Differential Equation Your anwer may differ from thoe given for thi exercie, depending on the method ued and roundoff error. Solve each differential equation and find the value of y aked for. Ue either the modified Euler method or the Runge-Kutta method, with a tep ize of 0... y xy 3; y(4). Find y(5).. y 5y 3x; y(0) 3. Find y(). 3. x 4 y ye x xy 3 ; y(). Find y(3). 4. xy y 3 ln x; y() 3. Find y() x y 8.3x in y.99y 0; y() 3. Find y(). 6. y ln x 3.99(x y) 4.8y; y().76. Find y(3). 7. (x y)y.76x y 5.33y; y(0). Find y(). 8. y.97xe x.9x co y 0; y(). Find y(3). 9. xy xy 8y lnhxyh; y() 3. Find y(). 0. x y 8.85x ln y 3.7x 0; y() 3.9. Find y(). Computer. Write a program or ue a preadheet to olve a firt-order differential equation by Euler method. Try it on any of the equation above.. Write a program or ue a preadheet to olve a firt-order differential equation by the modified Euler method. Try it on any of the equation above. 3. Write a program or ue a preadheet to olve a firt-order differential equation by the Runge-Kutta method. Try it on any of the equation above. 4. If you have done one of the calculation above uing a preadheet that ha a graphing utility, ue that utility to make a graph of x veru y over the given range. 5. Some computer algebra ytem have built-in function for numerically olving a differential equation. In Derive, the command Euler will invoke Euler method, and the command RK will call for the Runge-Kutta method. In Maple, the command firteuler, impeuler, and rungekutta are ued for olution by Euler method, the modified Euler method, and the Runge-Kutta method. Mathematica ue NDSolve for the numerical olution of a differential equation. Conult your manual for pecific intruction on how to ue thee command, and ue one or more to olve any of the equation in thi exercie et. 6 Numerical Solution of Second-Order Differential Equation To olve a econd-order differential equation numerically, we firt tranform the given equation into two firt-order equation by ubtituting another variable, m, for y.
25 36CH_PHCalter_TechMath_ //007 :8 PM Page 5 Section 6 Numerical Solution of Second-Order Differential Equation 5 Example 5: By making the ubtitution y m, we can tranform the econdorder equation into the two firt-order equation and y 3y 4xy 0 y m m 3m 4xy f(x, y, m) A before, we tart with the given point P(x p, y p ), at which we mut alo know the lope m p. We proceed from P to the next point Q(x q, y q ) a follow:. Compute the average firt derivative, m avg, over the interval x.. Compute the average econd derivative m avg over that interval. 3. Find the coordinate x q and y q at Q by uing the equation and x q x p x y q y p m avg x 4. Find the firt derivative m q at Q by uing m q m p m avg x 5. Return to Step and repeat the equence for the next interval. We can find m avg and m avg in Step and either by the modified Euler method or by the Runge-Kutta method. The latter method, a before, require u to find value at the intermediate point R and S. The equation for econd-order equation are imilar to thoe we had for firt-order equation. At P: x p, y p, and m p are given: m p f(x p, y p, m p ). At R: x r x p x y r y p m p x m r m p m p x m r f(x r, y r, m r ) At S: x x p x y y p m r x m m p m r x m f(x, y, m ) At Q: x q x p x y q y p m x m q m p m x m q f(x q, y q, m q ) Example 6: Find an approximate olution to the equation y y y e x with boundary condition y e at (, 0). Find point from to 3 in tep of 0..
26 36CH_PHCalter_TechMath_ //007 :8 PM Page 6 6 Solving Differential Equation by the Laplace Tranform and by Numerical Method Solution: Replacing y by m and olving for m give m f(x, y, m) e x y m At P: The boundary value are x p, y p 0, and m p Then Then and So and m p f(x p, y p, m p ) e 0 (5.4366) At R: x r x p x y r y p m p x (0.05) 0.78 m r m p m p x (0.05) 6.5 m r e (6.5) At S: x.05 y y p m r x 0 6.5(0.05) 0.36 m m p m r x (0.05) m e (6.3360) At Q: x q x p x. y q y p m x (0.) m q m p m x (0.) m q e (7.437) 9.86 m p m r m m q m avg (6.5) (6.3360) m p m r m m q m avg (7.9477) (8.0707) y q y p m avg x (0.) m q m p m avg x (0.) 7.40 The computation i then repeated. The remaining value are given in Table 5, along with the exact y and exact lope, found by analytical olution of the given equation.
27 36CH_PHCalter_TechMath_ //007 :8 PM Page 7 Review Problem 7 TABLE 5 x y Exact y Slope Exact Slope Exercie 6 Numerical Solution of Second-Order Differential Equation Solve each equation by the Runge-Kutta method, and find the value of y at x. Take a tep ize of 0. unit.. y y xy 3, y(, ). y xy 3xy y, y(, ) 3. 3y y ye x xy 3, y(, ) 4. y xy y 5 ln x, y(, 0) 5. xy.4x y 5.3x in y 7.4y 0, y(, ) y y ln x.69(x y) 6.6y, y(, 3.7) x y (x y)y 8.6x y.83y, y(, ) 8. y y 0.533xye x 0.36x co y 0, y(, 73.4) y xy 59.xy 74.y lnhxyh, y(, ) 0. xy x y.4x ln y 6.x 0, y(, 8) 5. Computer. Write a program or ue a preadheet to olve a econd-order differential equation by the Runge-Kutta method. Try your program on any of the equation above. If your preadheet ha a graphing utility, ue it to graph x veru y over the given range. REVIEW PROBLEMS Find the Laplace tranform by direct integration.. f(t) t. f(t) 3t 3. f(t) co 3t
28 36CH_PHCalter_TechMath_ //007 :8 PM Page 8 8 Solving Differential Equation by the Laplace Tranform and by Numerical Method Find the Laplace tranform by uing Table. 4. f(t) t 3 3t 5. f(t) 3te t 6. f(t) t 3e t 7. f(t) 3e t t 8. f(t) 3t 3 5t 4e t 9. y 3y, y(0) 0. y 6y, y(0). 3y y, y(0) 0. y y 3y, y(0) and y(0) 3 3. y 3y 4y, y(0) and y(0) 3 Find the invere tranform F() 4 5. F() ( 3) 6 6. F() ( ) 7. F() ( 5) 3 8. F() ( 4 ) 5 9. F() F() ( 3) 3. F(). F() F() 4 4 Solve each differential equation by the Laplace tranform. 4. y y 0, y(0) 5. y 3y t, y(0) 6. 3y y t, y(0) 7. y y y, y(0) 0 and y(0) 0 8. y 3y y t, y(0) and y(0) 9. y y 4t, y(0) 3 and y(0) y 3y y e t, y(0) and y(0) 3. A.5-lb ball i dropped from ret through air whoe reiting force i equal to.75 time the ball peed in ft/. Write an expreion for the peed of the ball. 3. A weight that hang from a pring i pulled down cm at t 0 and releaed from ret. The differential equation of motion i x 3.x Write the equation for x a a function of time.
29 36CH_PHCalter_TechMath_ //007 :8 PM Page 9 Review Problem The current in a certain RC circuit atifie the equation 8.4i 49i 00. If i i zero at t 0, write an equation for i. 34. The current in an RLC circuit atifie the equation 5.45i 80i 9730i 0. If i and i are zero at t 0, write an equation for i. 35. In an RLC circuit, R 3750, L 0.50 H, and C.5 µf. The circuit i connected to a dc ource of 50 V at t 0. If i i zero at t 0, write an equation for the intantaneou current. 36. In an RL circuit, R 3750 and L 0.50 H. The circuit i connected to an ac ource of 8.4 in 84t V at t 0. If i i zero at t 0, write an equation for the intantaneou current. Solve each differential equation. Ue any numerical method with a tep ize of 0. unit. 37. y xy 4y, y(). Find y(3). 38. xy 3y 5xy, y() 5. Find y() y y xy ln y, y(, ) 8. Find y(). 40. yy 7.yy.8x in x.y 0, y(, 5) 8.. Find y(). 4. y y ln y 5.7(x y).84x, y(, ) 3. Find y(). 4. x y ( y)y 3.67x 5.8x, y(, ) 3. Find y(). 43. y xy 6e x 6 co y 0, y(, ) 5. Find y(). 44. y xy 8 ln y, y() 7. Find y() y.85x in y 7.34x 0, y() 4.3. Find y(3). 46. yy y 77.y 8.4x lnhxyh, y(3, 5) 3. Find y(4). Writing 47. Suppoe that your company get a new project requiring you to olve many differential equation. You ee thi a a chance to get that computer you ve alway wanted. Write a memo to your bo explaining how the computer can be ued to olve thoe differential equation. Team Project 48. Given y y e t 0 where y(0), olve thi DE by all of the method at your dipoal: (a) analytically (b) graphically (c) numerically (d) by the Laplace tranform (e) by CAS Find y when t.
Chapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationLaplace Transformation
Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationSolutions to homework #10
Solution to homework #0 Problem 7..3 Compute 6 e 3 t t t 8. The firt tep i to ue the linearity of the Laplace tranform to ditribute the tranform over the um and pull the contant factor outide the tranform.
More informationR L R L L sl C L 1 sc
2260 N. Cotter PRACTICE FINAL EXAM SOLUTION: Prob 3 3. (50 point) u(t) V i(t) L - R v(t) C - The initial energy tored in the circuit i zero. 500 Ω L 200 mh a. Chooe value of R and C to accomplih the following:
More information1. /25 2. /30 3. /25 4. /20 Total /100
Circuit Exam 2 Spring 206. /25 2. /30 3. /25 4. /20 Total /00 Name Pleae write your name at the top of every page! Note: ) If you are tuck on one part of the problem, chooe reaonable value on the following
More informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
More informationReading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions
Chapter 4 Laplace Tranform 4 Introduction Reading aignment: In thi chapter we will cover Section 4 45 4 Definition and the Laplace tranform of imple function Given f, a function of time, with value f(t
More informationReading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions
Chapter 4 Laplace Tranform 4 Introduction Reading aignment: In thi chapter we will cover Section 4 45 4 Definition and the Laplace tranform of imple function Given f, a function of time, with value f(t
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationSIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2005. Initial Condition Source 0 V battery witch flip at t 0 find i 3 (t) Component value:
More informationName: Solutions Exam 3
Intruction. Anwer each of the quetion on your own paper. Put your name on each page of your paper. Be ure to how your work o that partial credit can be adequately aeed. Credit will not be given for anwer
More informationModeling in the Frequency Domain
T W O Modeling in the Frequency Domain SOLUTIONS TO CASE STUDIES CHALLENGES Antenna Control: Tranfer Function Finding each tranfer function: Pot: V i θ i 0 π ; Pre-Amp: V p V i K; Power Amp: E a V p 50
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationSIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. R 4 := 100 kohm
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2003. Cacaded Op Amp [DC&L, problem 4.29] An ideal op amp ha an output impedance of zero,
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationECE382/ME482 Spring 2004 Homework 4 Solution November 14,
ECE382/ME482 Spring 2004 Homework 4 Solution November 14, 2005 1 Solution to HW4 AP4.3 Intead of a contant or tep reference input, we are given, in thi problem, a more complicated reference path, r(t)
More informationQuestion 1 Equivalent Circuits
MAE 40 inear ircuit Fall 2007 Final Intruction ) Thi exam i open book You may ue whatever written material you chooe, including your cla note and textbook You may ue a hand calculator with no communication
More informationThe Laplace Transform , Haynes Miller and Jeremy Orloff
The Laplace Tranform 8.3, Hayne Miller and Jeremy Orloff Laplace tranform baic: introduction An operator take a function a input and output another function. A tranform doe the ame thing with the added
More informationTMA4125 Matematikk 4N Spring 2016
Norwegian Univerity of Science and Technology Department of Mathematical Science TMA45 Matematikk 4N Spring 6 Solution to problem et 6 In general, unle ele i noted, if f i a function, then F = L(f denote
More informationMAE140 Linear Circuits Fall 2012 Final, December 13th
MAE40 Linear Circuit Fall 202 Final, December 3th Intruction. Thi exam i open book. You may ue whatever written material you chooe, including your cla note and textbook. You may ue a hand calculator with
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationEE Control Systems LECTURE 14
Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationLinear Motion, Speed & Velocity
Add Important Linear Motion, Speed & Velocity Page: 136 Linear Motion, Speed & Velocity NGSS Standard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objective: 3.A.1.1, 3.A.1.3 Knowledge/Undertanding
More informationME 375 EXAM #1 Tuesday February 21, 2006
ME 375 EXAM #1 Tueday February 1, 006 Diviion Adam 11:30 / Savran :30 (circle one) Name Intruction (1) Thi i a cloed book examination, but you are allowed one 8.5x11 crib heet. () You have one hour to
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationChapter 7: The Laplace Transform Part 1
Chapter 7: The Laplace Tranform Part 1 王奕翔 Department of Electrical Engineering National Taiwan Univerity ihwang@ntu.edu.tw November 26, 213 1 / 34 王奕翔 DE Lecture 1 Solving an initial value problem aociated
More informationSECTION x2 x > 0, t > 0, (8.19a)
SECTION 8.5 433 8.5 Application of aplace Tranform to Partial Differential Equation In Section 8.2 and 8.3 we illutrated the effective ue of aplace tranform in olving ordinary differential equation. The
More informationDIFFERENTIAL EQUATIONS Laplace Transforms. Paul Dawkins
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationPHYS 110B - HW #2 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS 11B - HW # Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed [1.] Problem 7. from Griffith A capacitor capacitance, C i charged to potential
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationThe Laplace Transform
Chapter 7 The Laplace Tranform 85 In thi chapter we will explore a method for olving linear differential equation with contant coefficient that i widely ued in electrical engineering. It involve the tranformation
More informationConstant Force: Projectile Motion
Contant Force: Projectile Motion Abtract In thi lab, you will launch an object with a pecific initial velocity (magnitude and direction) and determine the angle at which the range i a maximum. Other tak,
More informationECE-202 FINAL December 13, 2016 CIRCLE YOUR DIVISION
ECE-202 Final, Fall 16 1 ECE-202 FINAL December 13, 2016 Name: (Pleae print clearly.) Student Email: CIRCLE YOUR DIVISION DeCarlo- 8:30-9:30 Talavage-9:30-10:30 2021 2022 INSTRUCTIONS There are 35 multiple
More informationLAPLACE TRANSFORM REVIEW SOLUTIONS
LAPLACE TRANSFORM REVIEW SOLUTIONS. Find the Laplace tranform for the following function. If an image i given, firt write out the function and then take the tranform. a e t inh4t From #8 on the table:
More informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
More informationLecture 3. January 9, 2018
Lecture 3 January 9, 208 Some complex analyi Although you might have never taken a complex analyi coure, you perhap till know what a complex number i. It i a number of the form z = x + iy, where x and
More informationThe Laplace Transform (Intro)
4 The Laplace Tranform (Intro) The Laplace tranform i a mathematical tool baed on integration that ha a number of application It particular, it can implify the olving of many differential equation We will
More informationFUNDAMENTALS OF POWER SYSTEMS
1 FUNDAMENTALS OF POWER SYSTEMS 1 Chapter FUNDAMENTALS OF POWER SYSTEMS INTRODUCTION The three baic element of electrical engineering are reitor, inductor and capacitor. The reitor conume ohmic or diipative
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationSolving Radical Equations
10. Solving Radical Equation Eential Quetion How can you olve an equation that contain quare root? Analyzing a Free-Falling Object MODELING WITH MATHEMATICS To be proficient in math, you need to routinely
More informationChapter 13. Root Locus Introduction
Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will
More informationPHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS B - HW #6 Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed. Problem. from Griffith Show that the following, A µo ɛ o A V + A ρ ɛ o Eq..4 A
More informationCHE302 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHE3 ECTURE V APACE TRANSFORM AND TRANSFER FUNCTION Profeor Dae Ryook Yang Fall Dept. of Chemical and Biological Engineering Korea Univerity CHE3 Proce Dynamic and Control Korea Univerity 5- SOUTION OF
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informatione st t u(t 2) dt = lim t dt = T 2 2 e st = T e st lim + e st
Math 46, Profeor David Levermore Anwer to Quetion for Dicuion Friday, 7 October 7 Firt Set of Quetion ( Ue the definition of the Laplace tranform to compute Lf]( for the function f(t = u(t t, where u i
More informationEECS2200 Electric Circuits. RLC Circuit Natural and Step Responses
5--4 EECS Electric Circuit Chapter 6 R Circuit Natural and Step Repone Objective Determine the repone form of the circuit Natural repone parallel R circuit Natural repone erie R circuit Step repone of
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More informationDelhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:
Serial : LS_N_A_Network Theory_098 Delhi Noida Bhopal Hyderabad Jaipur Lucknow ndore Pune Bhubanewar Kolkata Patna Web: E-mail: info@madeeay.in Ph: 0-4546 CLASS TEST 08-9 NSTRUMENTATON ENGNEERNG Subject
More informationCHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHBE3 ECTURE V APACE TRANSFORM AND TRANSFER FUNCTION Profeor Dae Ryook Yang Spring 8 Dept. of Chemical and Biological Engineering 5- Road Map of the ecture V aplace Tranform and Tranfer function Definition
More informationBASIC INDUCTION MOTOR CONCEPTS
INDUCTION MOTOS An induction motor ha the ame phyical tator a a ynchronou machine, with a different rotor contruction. There are two different type of induction motor rotor which can be placed inide the
More informationLecture 9: Shor s Algorithm
Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationChapter 7: The Laplace Transform
Chapter 7: The Laplace Tranform 王奕翔 Department of Electrical Engineering National Taiwan Univerity ihwang@ntu.edu.tw November 2, 213 1 / 25 王奕翔 DE Lecture 1 Solving an initial value problem aociated with
More informationPulsed Magnet Crimping
Puled Magnet Crimping Fred Niell 4/5/00 1 Magnetic Crimping Magnetoforming i a metal fabrication technique that ha been in ue for everal decade. A large capacitor bank i ued to tore energy that i ued to
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationAdvanced methods for ODEs and DAEs
Lecture : Implicit Runge Kutta method Bojana Roić, 9. April 7 What you need to know before thi lecture numerical integration: Lecture from ODE iterative olver: Lecture 5-8 from ODE 9. April 7 Bojana Roić
More informationName: Solutions Exam 2
Intruction. Anwer each of the quetion on your own paper. Put your name on each page of your paper. Be ure to how your work o that partial credit can be adequately aeed. Credit will not be given for anwer
More informationMath 273 Solutions to Review Problems for Exam 1
Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c
More informationME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004
ME 375 FINAL EXAM SOLUTIONS Friday December 7, 004 Diviion Adam 0:30 / Yao :30 (circle one) Name Intruction () Thi i a cloed book eamination, but you are allowed three 8.5 crib heet. () You have two hour
More informationFunction and Impulse Response
Tranfer Function and Impule Repone Solution of Selected Unolved Example. Tranfer Function Q.8 Solution : The -domain network i hown in the Fig... Applying VL to the two loop, R R R I () I () L I () L V()
More informationMidterm Review - Part 1
Honor Phyic Fall, 2016 Midterm Review - Part 1 Name: Mr. Leonard Intruction: Complete the following workheet. SHOW ALL OF YOUR WORK. 1. Determine whether each tatement i True or Fale. If the tatement i
More informationHomework #7 Solution. Solutions: ΔP L Δω. Fig. 1
Homework #7 Solution Aignment:. through.6 Bergen & Vittal. M Solution: Modified Equation.6 becaue gen. peed not fed back * M (.0rad / MW ec)(00mw) rad /ec peed ( ) (60) 9.55r. p. m. 3600 ( 9.55) 3590.45r.
More informationMath 201 Lecture 17: Discontinuous and Periodic Functions
Math 2 Lecture 7: Dicontinuou and Periodic Function Feb. 5, 22 Many example here are taken from the textbook. he firt number in () refer to the problem number in the UA Cutom edition, the econd number
More information7-5-S-Laplace Transform Issues and Opportunities in Mathematica
7-5-S-Laplace Tranform Iue and Opportunitie in Mathematica The Laplace Tranform i a mathematical contruct that ha proven very ueful in both olving and undertanding differential equation. We introduce it
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationDesign of Digital Filters
Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function
More informationLecture 4 : Transform Properties and Interpretations. Continued to the Next (Higher) Level. 1. Example 1. Demo of the mult-by-t property.
Lecture 4 : Tranform Propertie and Interpretation Continued to the Next (Higher) Level 1. Example 1. Demo of the mult-by-t property. (i) Conider above graph of f (t) and g(t) = tf (t). Set K = 1. (ii)
More informationCumulative Review of Calculus
Cumulative Review of Calculu. Uing the limit definition of the lope of a tangent, determine the lope of the tangent to each curve at the given point. a. f 5,, 5 f,, f, f 5,,,. The poition, in metre, of
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationEE C128 / ME C134 Problem Set 1 Solution (Fall 2010) Wenjie Chen and Jansen Sheng, UC Berkeley
EE C28 / ME C34 Problem Set Solution (Fall 200) Wenjie Chen and Janen Sheng, UC Berkeley. (0 pt) BIBO tability The ytem h(t) = co(t)u(t) i not BIBO table. What i the region of convergence for H()? A bounded
More informationV = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr
0.1 Related Rate In many phyical ituation we have a relationhip between multiple quantitie, and we know the rate at which one of the quantitie i changing. Oftentime we can ue thi relationhip a a convenient
More informationECE Linear Circuit Analysis II
ECE 202 - Linear Circuit Analyi II Final Exam Solution December 9, 2008 Solution Breaking F into partial fraction, F 2 9 9 + + 35 9 ft δt + [ + 35e 9t ]ut A 9 Hence 3 i the correct anwer. Solution 2 ft
More informationDelhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:
Serial : LS_B_EC_Network Theory_0098 CLASS TEST (GATE) Delhi Noida Bhopal Hyderabad Jaipur Lucknow ndore Pune Bhubanewar Kolkata Patna Web: E-mail: info@madeeay.in Ph: 0-4546 CLASS TEST 08-9 ELECTRONCS
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More informationRoot Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0
Objective Root Locu Diagram Upon completion of thi chapter you will be able to: Plot the Root Locu for a given Tranfer Function by varying gain of the ytem, Analye the tability of the ytem from the root
More informationIII.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES
III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation
More informationLecture 6: Resonance II. Announcements
EES 5 Spring 4, Lecture 6 Lecture 6: Reonance II EES 5 Spring 4, Lecture 6 Announcement The lab tart thi week You mut how up for lab to tay enrolled in the coure. The firt lab i available on the web ite,
More informationNCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011
NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic
More informationSOLUTIONS FOR HOMEWORK SECTION 6.4 AND 6.5
SOLUTIONS FOR HOMEWORK SECTION 6.4 AND 6.5 Problem : For each of the following function do the following: (i) Write the function a a piecewie function and ketch it graph, (ii) Write the function a a combination
More informationUniform Acceleration Problems Chapter 2: Linear Motion
Name Date Period Uniform Acceleration Problem Chapter 2: Linear Motion INSTRUCTIONS: For thi homework, you will be drawing a coordinate axi (in math lingo: an x-y board ) to olve kinematic (motion) problem.
More informationName: Solutions Exam 2
Name: Solution Exam Intruction. Anwer each of the quetion on your own paper. Put your name on each page of your paper. Be ure to how your work o that partial credit can be adequately aeed. Credit will
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationME 375 FINAL EXAM Wednesday, May 6, 2009
ME 375 FINAL EXAM Wedneday, May 6, 9 Diviion Meckl :3 / Adam :3 (circle one) Name_ Intruction () Thi i a cloed book examination, but you are allowed three ingle-ided 8.5 crib heet. A calculator i NOT allowed.
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More information18.03SC Final Exam = x 2 y ( ) + x This problem concerns the differential equation. dy 2
803SC Final Exam Thi problem concern the differential equation dy = x y ( ) dx Let y = f (x) be the olution with f ( ) = 0 (a) Sketch the iocline for lope, 0, and, and ketch the direction field along them
More information4e st dt. 0 e st dt. lim. f (t)e st dt. f (t) e st dt + 0. f (t) e. e (2 s)t dt + 0. e (2 s)4 1 ] = 1 = 1. te st dt + = t s e st
Worked Solution Chapter : The Laplace Tranform 6 a F L4] 6 c F L f t] 4 4e t dt e t dt 4 e t 4 ] e t e 4 if > 6 e F L f t] 6 g Uing integration by part, f te t dt f t e t dt + e t dt + e t + 4 4 4 f te
More informationDepartment of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002
Department of Mechanical Engineering Maachuett Intitute of Technology 2.010 Modeling, Dynamic and Control III Spring 2002 SOLUTIONS: Problem Set # 10 Problem 1 Etimating tranfer function from Bode Plot.
More information4e st dt. 0 e st dt. lim. f (t)e st dt. f (t) e st dt + 0. e (2 s)t dt + 0. e (2 s)4 1 = 1. = t s e st
Worked Solution 8 Chapter : The Laplace Tranform 6 a F L] e t dt e t dt e t ] lim t e t e if > for > 6 c F L f t] f te t dt f t e t dt + e t dt + e t + f t e t dt e t dt ] e e ] 6 e F L f t] f te t dt
More information5.5 Application of Frequency Response: Signal Filters
44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order
More informationSuggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall
Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall
More informationRiemann s Functional Equation is Not Valid and its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.
Riemann Functional Equation i Not Valid and it Implication on the Riemann Hypothei By Armando M. Evangelita Jr. On November 4, 28 ABSTRACT Riemann functional equation wa formulated by Riemann that uppoedly
More informationMidterm Test Nov 10, 2010 Student Number:
Mathematic 265 Section: 03 Verion A Full Name: Midterm Tet Nov 0, 200 Student Number: Intruction: There are 6 page in thi tet (including thi cover page).. Caution: There may (or may not) be more than one
More information