082-Engineering-Mathematics-Anthony-Croft-Robert-Davison-Martin-Hargreaves-James-Flint-Edisi-5-2017
610 Chapter 20 Ordinary differential equations IIy CCRi~LFigure20.6Thecircuitcould be modelled using v C andiasthe statevariables.each other. For example, for the electrical system of Figure 20.6, the choice of v Cand3v Cwouldleadtothetwostatevariablesbeingdependentoneachother.Avalidchoicewouldbe v Candiwhicharenotdirectlydependentoneachother.Thisisasecond-ordersystemand so only requires two statevariables tomodel itsbehaviour.Many engineering systems may have a high order and so require several differentialequations to model their behaviour. For this reason, a standard way of laying out theseequations has evolved to reduce the chance of making errors. There is also an addedadvantage in that the standard layout makes it easier to present the equations to a digitalcomputer for solution. Before introducing the standard form, an example will bepresented toillustratethe statevariable method.Engineeringapplication20.1State-spacemodelforaspring-mass-dampersystemConsiderthemechanicalsystemillustratedinFigure20.7.Amassrestsonafrictionlesssurfaceandisconnectedtoafixedwallbymeansofanidealspringandanidealdamper. Aforce, f (t), isapplied tothe mass and the position ofthe mass is w.KBMwf(t)Figure20.7Asecond-order mechanicalsystem.By considering the forces acting on the mass,M, it is possible to devise a differentialequation that models the behaviour of the system. The force produced by thespring isKw (K is the spring stiffness) and opposes the forward motion. The forceproduced by the damper isB dw , where B is the damping coefficient, and this alsodtopposestheforwardmotion.Sincethemassisconstant,Newton’ssecondlawofmotionstates that the net force on the mass is equal to the product of the mass and itsacceleration. Thus wefindf(t)−B dw −Kw=M d2 w(20.8)dt dt 2Note that this is a second-order differential equation and so the system is a secondordersystem.Inordertoobtainastate-spacemodelthestatevariableshavetobechosen.Thereareseveralpossiblechoices.Anobviousoneistheposition, w,ofthemass.Asecond
20.4 State-space models 611statevariableisrequiredasthesystemissecondorder.Inthiscasethevelocityofthemasswillbechosen.Thevelocityofthemassisnotdirectlydependentonitspositionand so the two variables are independent. Another possible choice would have beendw−w.Thismayseemaclumsychoicebutforcertainproblemssuchchoicesmaydtlead tosimplifications inthe statevariable equations.Itiscustomarytousethesymbolsx 1,x 2,x 3,...torepresentvariablesforreasonsthat will become clear shortly. So,x 1=w(20.9a)x 2= dw(20.9b)dtBecause of the particular choice of state variables, it is easy to obtain the first ofthefirst-orderdifferentialequations--thusillustratingtheneedforexperiencewhenchoosing statevariables.Differentiating Equation (20.9a) givesdx 1= dw =xdt dt 2This isthe first stateequation although itisusually written asẋ 1=x 2where ẋ 1denotes dx 1. The second of the first-order equations is obtained by rearrangingEquationdt(20.8):d 2 w= − K dt 2 M w − B dw+ 1 f (t) (20.10)M dt MHowever, differentiating Equation (20.9b) we getd 2 wdt 2 =ẋ 2Then, usingEquation (20.10) weobtainẋ 2=− K M x 1 − B M x 2 + 1 M f(t)Finally, itisusual toarrange the stateequations inaparticular way:ẋ 1= + x 2ẋ 2=− K M x 1 − B M x 2 + 1 M f(t)w=x 1Note that it is conventional to relate the output of the system to the state variables.Assume that for this system the required output variable is the position of the mass,that is w.Itisstraightforward torewrite these equations inmatrix form:) ( )( ) ( )(ẋ1 0 1 x1 0=+ f(t)ẋ2 −K/M −B/M x 21/Mw = ( 1 0 ) ( )x 1x 2
- Page 579 and 580: 19.5 Second-order linear equations
- Page 581 and 582: 19.6 Constant coefficient equations
- Page 583 and 584: 19.6 Constant coefficient equations
- Page 585 and 586: 19.6 Constant coefficient equations
- Page 587 and 588: 19.6 Constant coefficient equations
- Page 589 and 590: 19.6 Constant coefficient equations
- Page 591 and 592: thatis,β = 332The required particu
- Page 593 and 594: 19.6 Constant coefficient equations
- Page 595 and 596: 19.6 Constant coefficient equations
- Page 597 and 598: 19.6 Constant coefficient equations
- Page 599 and 600: 19.6 Constant coefficient equations
- Page 601 and 602: 19.6 Constant coefficient equations
- Page 603 and 604: 19.6 Constant coefficient equations
- Page 605 and 606: 19.7 Series solution of differentia
- Page 607 and 608: 19.8 Bessel’s equation and Bessel
- Page 609 and 610: 19.8 Bessel’s equation and Bessel
- Page 611 and 612: fromwhicha 5=− 124 a 3 = 1192 a 1
- Page 613 and 614: n 5 2n 5 3n5 419.8 Bessel’s equat
- Page 615 and 616: 19.8 Bessel’s equation and Bessel
- Page 617 and 618: 19.8 Bessel’s equation and Bessel
- Page 619 and 620: 19.8 Bessel’s equation and Bessel
- Page 621 and 622: Review exercises 19 601quencies, kn
- Page 623 and 624: 20 OrdinarydifferentialequationsIIC
- Page 625 and 626: 20.2 Analogue simulation 605y iy i
- Page 627 and 628: 20.3 Higher order equations 607High
- Page 629: 20.4 State-space models 609Solution
- Page 633 and 634: 20.4 State-space models 613R a L aT
- Page 635 and 636: 20.5 Numerical methods 615Combining
- Page 637 and 638: 20.6 Euler’s method 617and we can
- Page 639 and 640: 20.6 Euler’s method 619Whenx = 1,
- Page 641 and 642: 20.7 Improved Euler method 621Table
- Page 643 and 644: 20.8 Runge--Kutta method of order 4
- Page 645 and 646: 20.8.1 Higherorderequations20.8 Run
- Page 647 and 648: 21 TheLaplacetransformContents 21.1
- Page 649 and 650: 21.3 Laplace transforms of some com
- Page 651 and 652: 21.4 Properties of the Laplace tran
- Page 653 and 654: (b) Usethe first shift theorem with
- Page 655 and 656: 21.5 Laplace transform of derivativ
- Page 657 and 658: 21.5 Laplace transform of derivativ
- Page 659 and 660: 21.6 Inverse Laplace transforms 639
- Page 661 and 662: 21.7 Using partial fractions to fin
- Page 663 and 664: 21.8 Finding the inverse Laplace tr
- Page 665 and 666: 21.8 Finding the inverse Laplace tr
- Page 667 and 668: 21.9 The convolution theorem 64721.
- Page 669 and 670: 21.9 Solving linear constant coeffi
- Page 671 and 672: 21.10 Solving linear constant coeff
- Page 673 and 674: 21.10 Solving linear constant coeff
- Page 675 and 676: X(s)= s3 −3s 2 −4s+6s 2 (s 2
- Page 677 and 678: 21.10 Solving linear constant coeff
- Page 679 and 680: 21.11 Transfer functions 659(c) x
20.4 State-space models 611
statevariableisrequiredasthesystemissecondorder.Inthiscasethevelocityofthe
masswillbechosen.Thevelocityofthemassisnotdirectlydependentonitsposition
and so the two variables are independent. Another possible choice would have been
dw
−w.Thismayseemaclumsychoicebutforcertainproblemssuchchoicesmay
dt
lead tosimplifications inthe statevariable equations.
Itiscustomarytousethesymbolsx 1
,x 2
,x 3
,...torepresentvariablesforreasons
that will become clear shortly. So,
x 1
=w
(20.9a)
x 2
= dw
(20.9b)
dt
Because of the particular choice of state variables, it is easy to obtain the first of
thefirst-orderdifferentialequations--thusillustratingtheneedforexperiencewhen
choosing statevariables.
Differentiating Equation (20.9a) gives
dx 1
= dw =x
dt dt 2
This isthe first stateequation although itisusually written as
ẋ 1
=x 2
where ẋ 1
denotes dx 1
. The second of the first-order equations is obtained by rearranging
Equation
dt
(20.8):
d 2 w
= − K dt 2 M w − B dw
+ 1 f (t) (20.10)
M dt M
However, differentiating Equation (20.9b) we get
d 2 w
dt 2 =ẋ 2
Then, usingEquation (20.10) weobtain
ẋ 2
=− K M x 1 − B M x 2 + 1 M f(t)
Finally, itisusual toarrange the stateequations inaparticular way:
ẋ 1
= + x 2
ẋ 2
=− K M x 1 − B M x 2 + 1 M f(t)
w=x 1
Note that it is conventional to relate the output of the system to the state variables.
Assume that for this system the required output variable is the position of the mass,
that is w.
Itisstraightforward torewrite these equations inmatrix form:
) ( )( ) ( )
(ẋ1 0 1 x1 0
=
+ f(t)
ẋ2 −K/M −B/M x 2
1/M
w = ( 1 0 ) ( )
x 1
x 2