ECM Control Engineering Dr Mustafa M Aziz (2013) SYSTEM RESPONSE

Size: px
Start display at page:

Download "ECM Control Engineering Dr Mustafa M Aziz (2013) SYSTEM RESPONSE"

Transcription

1 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) SYSTEM RESPONSE. Itroductio. Repoe Aalyi of Firt-Order Sytem 3. Secod-Order Sytem 4. Siuoidal Repoe of the Sytem 5. Bode Diagram 6. Baic Fact About Egieerig Sytem. Itroductio The order of a ytem i defied a beig the highet power of derivative i the differetial equatio, or beig the highet power of i the deomiator of the trafer fuctio. A firt-order ytem oly ha to the power oe i the deomiator, while a ecod-order ytem ha the highet power of i the deomiator beig two. Type of the iput fuctio (or tet iput igal) commoly ued are: Impule fuctio: I the time domai, u(t) cδ(t). I the domai, U() c. Step fuctio: I the time domai, u(t) c. I the domai, U() c/. Ramp fuctio: I the time domai, u(t) ct. I the domai, U() c/. Siuoidal fuctio: I the time domai, u(t) ci(ωt). I the domai, U() cω/( +ω ). where c i a cotat i all the above. With thee tet igal, mathematical ad experimetal aalye of cotrol ytem ca be carried out eaily ice the igal are very imple fuctio of time. Which of thee typical igal to ue for aalyig ytem characteritic may be determied by the form of the iput that the ytem will be ubjected to mot frequetly uder ormal operatio. If the iput to a cotrol ytem are gradually chagig fuctio of time, the a ramp fuctio of time may be a good tet igal. Similarly, if a ytem i ubjected to udde diturbace, a tep fuctio of time may be a good tet igal, ad for a ytem ubjected to a hock iput, a pule or a impule fuctio may be bet. Exercie: What are the order of the ytem decribed by the followig trafer fuctio: a) G() m + b + k b) G() RC + c) G() LC + RC +

2 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) The time repoe of a cotrol ytem coit of two part: the traiet repoe ad the teadytate repoe. The traiet repoe i defied a the part of the time repoe which goe from the iitial tate to the fial tate ad reduce to zero a time become very large. The teady-tate repoe i defied a the behaviour of the ytem a t approache ifiity after the traiet have died out. Thu the ytem repoe may be writte a: y t (t) + y (t) where y t (t) deote the traiet repoe, ad y (t) deote the teady-tate repoe.. Repoe Aalyi of Firt-Order Sytem May ytem are approximately firt-order. The importat feature i that the torage of ma, mometum ad eergy ca be captured by oe parameter. Example of firt-order ytem are velocity of a car o the road, cotrol of the velocity of a rotatig ytem, electric ytem where eergy torage i eetially i oe capacitor or oe iductor, icompreible fluid flow i a pipe, level cotrol of a tak, preure cotrol i a ga tak, temperature i a body with eetially uiform temperature ditributio (e.g. team filled veel). Next we will preet everal example to how how to obtai the dyamic equatio of firt-order ytem. Example : Mechaical ytem m i the ma, u(t) i the exteral force, i the velocity ad b i the frictio coefficiet. By Newto law, we have the followig differetial equatio: u(t) m d m + b u(t) dt b Example : Electrical ytem R i the reitace, C i the capacitace, u(t) i the iput voltage ad i the output voltage. By Kirchhoff law: R u(t) Ri(t) + ad i(t) Cd/dt u(t) C Thu d RC + u(t) dt A geeral form of a firt-order ytem ca be repreeted by the block diagram. R() + _ /T Y() R() T + Y()

3 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) Tak: Write the ytem output or repoe to iput uch a the uit-tep, uit-ramp, ad uitimpule fuctio, repectively. The iitial coditio are aumed to be zero. Draw the repoe curve. T i the time cotat of the ytem... Uit-tep repoe of firt-order ytem R() /, ad therefore the uit-tep repoe i: Y() (T + ) T Expadig Y() ito partial fractio: Y() T + + / T Take the ivere Laplace traform: - e -t/t, t. The olutio ha two part: a teady-tate repoe: y (t), ad a traiet repoe: t / T y (t) e, which decay to zero a t. t.9.8 Slope /T Uit-tep repoe, T T r t / T T The lope of the taget lie at t i /T. Pole locatio i the plae: -/T. At t T, y(t) e T i called the time cotat, ad it i the time it take for the tep repoe to rie to 63.% of it fial value. y(t).865; y(3t).95; y(4t).98; y(5t) It ca be ee that for t 4T, the repoe remai withi % of the fial value; thi time i kow a the ettlig time, T. 3

4 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) The rie time, T r, i defied a the time for the waveform to go from % to 9% of it fial value. The teady-tate error i the error after the traiet repoe ha decayed leavig oly the cotiuou repoe. The error igal: e(t) r(t) e -t/t e -t/t A t approache ifiity, e -t/t approache zero ad the teady-tate error i: e e( ) lim r(t) t [ ] The larger the time cotat T i, the lower the ytem repoe i. It i oted that the traiet repoe domiate the repoe of the ytem at time immediately after the iput i applied ad ca make igificat cotributio to the ytem repoe whe the time cotat i large. Y() Exercie: A RC circuit ha the followig trafer fuctio: R() + 4 For a tep iput r(t) V, what i the time take for the output of the RC circuit to reach 95% of it teady-tate repoe? Y() 5 Exercie: A ytem ha trafer fuctio: R() + 5 Fid the time cotat, T, the ettlig time, T, ad the rie time, T r for a uit-tep iput. 4

5 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3).. Uit-ramp repoe of firt-order ytem The iput, r(t) t for t. 6 Uit-ramp repoe, T Laplace traform: R() /. 5 The output traform: () Y T + 4 Steady-tate error Expadig Y() ito partial fractio: T T Y() + T + 3 r(t) t Takig the ivere Laplace traform: t - T + Te -t/t, t. Steady-tate error: e e( ) lim[ r(t) ] t T t / T.3. Uit-impule repoe of firt-order ytem The uit-impule iput, r(t) δ(t), t t δ( t) t. /T Uit-impule repoe, T Laplace traform: R()..8 The output traform: Y() T + Takig the ivere Laplace traform: e -t/t /T, t. Note that the impule iput yield the trafer fuctio of the ytem a output t / T 5

6 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3. Secod-Order Sytem Example : Mechaical ytem For the mechaical ytem how i the figure, m i the ma, k i the prig cotat, b i the frictio coefficiet, u(t) i the exteral force ad i the diplacemet. From Newto ecod law force ma: d d + b + k u(t) dt dt m k b m u(t) Example : Electrical ytem: RLC circuit Uig Kirchhoff law: di(t) u (t) Ri(t) + L + dt d where i (t) L dt Hece: d LC dt + u(t) R d RC + u(t) dt L C A geeral form of a ecod order ytem i: d + ζω + ω dt dt d kω u(t) Trafer fuctio: Y() R() kω + ζω + ω k: the gai of the ytem ζ: the dampig ratio of the ytem ω : the (udamped) atural frequecy of the ytem Solutio (root, or pole of the ytem) of the characteritic equatio are: Three cae: ζω ω ζ ad ζω + ω ζ ζ, critically damped cae ζ >, overdamped cae < ζ <, uderdamped cae We will tudy the above cae whe k for implicity. 6

7 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3.. Step Repoe of Secod-Order Sytem 3... Critically damped cae (ζ ) Two equal pole: -ζω ω For a uit-tep iput R() /, the output i: Y() ( + ω ) ω Expadig Y() ito partial fractio: Y() + ω ( + ω ) + t Takig the ivere Laplace traform: e ( t) Steady-tate error: e( ) Uit-tep repoe of d-order ytem (critically damped cae) -plae jω ω 3 ω -ζω σ Exercie: A ytem ha the followig trafer fuctio: Y() R() t (ec) What i the tate of dampig of the ytem whe it i ubjected to a uit-tep iput? Determie the atural frequecy of the ytem. 7

8 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3... Overdamped cae (ζ > ) We ca write the trafer fuctio of a ecod-order ytem by factorig the deomiator a: Y() ω ω R() ( + ζω + ω ζ )( + ζω ω ζ ) ( )( ) Two pole i the -plae: ζω ω ζ ad ζω + ω ζ For a uit-tep iput, the output i: ω Y() ( )( Takig the ivere Laplace traform yield the time repoe (prove thi time repoe a a exercie): t t e e Whe ζ i much greater tha uity, i.e. ζ >>, the >> ad the term ivolvig i the time repoe will decay fater tha the term ivolvig. The term ivolvig ca therefore be eglected ad the ytem become firt-order decided maily by the pole : Y() / >> ω ω ω / R() ( )( ) ( / )( ) ) Sice ω, the trafer fuctio become: Y() R() The uit-tep time repoe i: ( )t e, t Uit-tep repoe of d order ytem (overdamped cae), ζ, ω.9.8 e t -plae jω ω ( ζ ω ( ζ + ζ ζ ) ) σ e t e t ζ>>, ca be eglected. t (ec) 8

9 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) Uit-tep repoe of d order ytem (overdamped cae), ζ., ω e t e t e t t (ec) ζ, caot be eglected Uderdamped cae ( < ζ < ) Trafer fuctio: Y() R() ( + ζω ω + jω )( + ζω d jω d ) where ω d ω ζ i called the damped atural frequecy. The two pole are: ζω jω ζ ad ζω + ω j ζ -ζω +jω d jω ω d -plae ζ θ ta ζ co(θ) ζ θ -ζω ω σ For uit-tep iput R() /, the output i: -ζω -jω d -ω d Y() + ζω + ζω + ω + ζω ( + ζω ) + ω d ζω ( + ζω ) + ω d Takig the ivere Laplace traform uig the table of Laplace traform yield: - + ζω ζωt e co( ω ( + ζω ) + ωd - dt) ζω ζωt e i( ω ( + ζω ) + ωd d t) 9

10 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) The time repoe i: e where ζ θ ta. ζ ζω t e ζω t ζ co( ω d i( ω t) + d t + θ) ζ ζ i( ω d t) Whe ζ, the repoe become udamped ad ocillatio cotiue idefiitely at frequecy ω. The time repoe i thi cae become: co(ω t) The atural udamped frequecy, ω, i the frequecy of ocillatio of the ytem without dampig Uit-tep repoe of d order ytem (uderdamped cae), ω 3 ζ ζ. ζ.5 ζ t (ec)

11 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3.. Traiet ytem pecificatio.4 Uit-tep repoe of a d order ytem, ζ.9, ω 3. Maximum overhoot T r T p T t (ec) 3... Maximum overhoot The maximum amout by which the ytem output repoe proceed beyod the deired repoe. Let y max deote the maximum value of, ad y y( ) the teady-tate value of, the the maximum overhoot of i defied a: maximum overhoot y max - y The maximum overhoot i ofte repreeted by a percetage of the fial value of the tep repoe: y max y ζπ PO ( percet overhoot) % exp y ζ 3... Peak time, T p The time required for the repoe to reach the firt peak of the overhoot: π Tp ω ζ Rie time, T r The time required for the tep repoe to rie from % to 9% of it fial value for critical ad overdamped cae, ad from % to % for uderdamped cae. For the uderdamped cae: T r π θ ω d π ta ω ( ζ ζ / ζ)

12 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) Settlig time, T The time required for the tep repoe to ettle withi a certai percetage of it fial value. A frequetly ued figure i % i which cae the ettlig time i approximately: 4 T ζω 3 For varyig withi 5% of the fial value, the ettig time i: T ζω Note that NOT all thee pecificatio ecearily apply to ay give cae. For example, for a overdamped ytem, the term peak time ad maximum overhoot do ot apply. The traiet behaviour of a ecod-order ytem ca be decribed by: the wifte of the repoe, a repreeted by T r ad T p the cloee of the repoe to the deired repoe, a repreeted by PO ad T. From the deig requiremet, the wifter ad cloer, the better. However, whe ω i fixed, mall T r ad T p require a mall ζ, while mall T ad PO require a large ζ. Note that ( ζπ / ζ ) PO exp, T π/[ ω ζ ], T [ π ta ( ζ / ζ)]/[ ω ζ ]. p Thee lead to coflictig requiremet. A compromie mut be obtaied ometime. Exercie: A ecod-order ytem i uderdamped with a dampig ratio of.4 ad a atural frequecy of Hz. Fid: a) the trafer fuctio b) the time repoe whe it i ubjected to a uit-tep iput c) the percetage overhoot with uch a iput d) the rie time r

13 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) Exercie: Give the trafer fuctio: fid T p, PO, T ad T r. G () Ramp repoe of a ecod-order ytem The Laplace traform of a uit-ramp iput i R() / The output i: kω Y() ( + ζω + ω ) Agai, there are three cae: ζ, critically damped cae ζ >, overdamped cae < ζ <, uderdamped cae 3

14 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3.4. Impule repoe of a ecod-order ytem The Laplace traform of a uit-impule iput i R(). The output traform i therefore equal to the trafer fuctio of the ytem, i.e.: kω Y() + ζω + ω Fact: the uit-impule fuctio i the time derivative of the uit-tep fuctio. Therefore, the impule repoe of a LTI ytem ca be foud from the time derivative of the tep repoe for a give dampig. Takig the example of a critically damped ytem where ζ, the uit-tep repoe i give by: t e ( + t) The uit-impule repoe i therefore: d dt ω te t Note, agai, that the impule iput give the trafer fuctio of the ytem. 4

15 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 4. Siuoidal Repoe of the Sytem Although tep repoe are commoly ued i both imulatio ad experimetal tet, it i alo commo to udertake frequecy repoe tet o the ytem. The frequecy repoe of a ytem i defied a the teady-tate repoe of the ytem to a iuoidal iput igal. u(t) U o i(ωt) y (t) U o G(jω) i(ωt+φ) U o G(jω) U o φ t The liear, time-ivariat ytem G() ubjected to a iuoidal iput of amplitude U o ad frequecy ω decribe by: u(t) U o i(ωt) will, at teady tate, have a iuoidal output of the ame frequecy a the iput but, geerally, with differet amplitude ad phae give by: y (t) U o G(jω) i(ωt + φ(ω)) where U o G(jω) i the amplitude of the output ie wave: G(jω ) {Re[G( jω)]} + {Im[G( jω ad φ(ω) i the phae hift i radia or degree give by: Im[G( jω)] φ( ω) ta G( jω) Re[G( jω)] The iuoidal trafer fuctio of ay liear ytem i obtaied by ubtitutig jω for i the trafer fuctio of the ytem. )]} Proof: Coider a ytem decribed by: Y() G() U() u(t) U() G() y() Y() The iput u(t) i a ie wave with ad amplitude U o ad frequecy ω: u(t) U o i(ωt) U oω The Laplace traform of u(t) i: U() + ω 5

16 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) With zero iitial coditio, the Laplace traform of the output i: U oω Y() G() + ω A partial fractio expaio of a geeral ytem (aumig the pole of G() are ditict) yield: * c c c c c Y() + + L p + p + p + jω jω Partial fractio term from G() * where c are cotat ad c ad c are a complex cojugate pair that ca be obtaied uig the cover-up rule: G( jω)u o * G(jω)U o c c j j Sice G(jω) i a complex quatity, it ca be writte i the form: G(jω) G(jω) e jφ where G(jω) i the magitude ad φ i the phae give repectively by: Im[G( jω)] G(jω ) {Re[G( jω)]} + {Im[G( jω)]} φ( ω) ta G( jω) Re[G( jω)] Similarly: G(-jω) G(-jω) e -jφ G(jω) e -jφ Therefore: c G( jω) e j jφ U o c * G(jω) e j jφ U o The time repoe that correpod to Y() i: p t pt c e + c e + L + c pt jωt * j e + ce + c e ω pt If all the pole of the ytem repreet a table behaviour, the atural uforced repoe ( c e decay to zero at t ) will die out evetually ad therefore the teady-tate repoe of the ytem will be due olely to the iuoidal term which i caued by the iuoidal excitatio, i.e. jωt * jωt y (t) lim c e + c e Subtitutig for c ad t * c ad otig that i(x) (e jx e -jx )/(j), give the teady-tate output: y (t) U o G(jω) i(ωt + φ) t Advatage of frequecy domai aalyi: The trafer fuctio of complicated compoet ca be determied experimetally by frequecy repoe tet (without derivig their mathematical model) uig available igal geerator ad precie meauremet equipmet (e.g. pectrum aalyer). 6

17 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) The amplitude ad phae of the frequecy repoe ca be ued to predict both time-domai traiet ad teady-tate ytem performace. Sytem may be deiged to achieve traiet ad teady-tate requiremet uig frequecy repoe aalyi, ad uch aalyi ad deig may be exteded to certai oliear cotrol ytem. Exercie: For the iuoidal iput u(t) i(t) applied to the ytem: determie the teady-tate output of the ytem. G(), +. Liear Simulatio Reult Amplitude Time (ec.) 7

18 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) There are two commoly ued repreetatio of iuoidal trafer fuctio: ) Nyquit or polar plot, ad ) Bode diagram We hall focu o the more popular Bode aalyi ad how how we ca ue MATLAB to produce thee plot. 5. Bode Diagram (Aymptotic Approximatio) The frequecy repoe of a ytem G() ca be decribed a: G() jω G( jω) Re[G( jω)] + Im[G( jω)] R( ω) + jx( ω) A Bode diagram coit of two graph: - plot of the logarithm of the magitude of the a iuoidal trafer fuctio, G(jω) - plot of the phae agle, φ(ω) both are plotted agait the frequecy ω (rad/) o a logarithmic (bae ) cale. Logarithmic magitude (alo called gai) of G(jω), M log G(jω) (uit i decibel, db) Phae: φ( ω) ta X( ω) R( ω) Advatage of Bode diagram: Bode plot of ytem i erie imply add, which i quite coveiet. For example, coider the trafer fuctio: b m ( z)( z ) L( z m ) G() ( p )( p ) L( p ) The magitude of the frequecy repoe of the ytem i give by: b m ( z) ( z ) L ( z m ) G(jω) ( p ) ( p ) L ( p ) Takig the logarithm yield: jω log G(jω) log b m + log ( - z ) + log ( z ) + - log ( - p ) - log ( p ) - jω Kowig the repoe of each term, the algebraic um would give the total repoe i db. Bode' importat phae-gai relatiohip i plotted uig aymptotic approximatio o a logarithmic cale, which mea a wider rage of ytem behaviour, from low to high frequecie, ca be diplayed o a igle plot. Dyamic compeator (cotroller) deig ca be baed etirely o Bode plot. 8

19 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) Example: Fid the Bode plot for the followig RC filter with: G(jω ) RC( jω) + jωt + R Solutio: u(t) C Magitude: G(jω ) + ω T / T / T + ω Logarithmic Gai (db) log ω / G(jω) log (/ T) log (/ T + ) At low frequecie (ω ): gai (db) log (/ T) log (/ T) db At high frequecie (ω ): gai (db) log (/ T) log ( ) ω If we plot the logarithmic gai (db) agait log ( ω ), the the above equatio become the traight lie: y c + mx where y log G(jω ), c itercept, m lope, ad x log ( ω). Thu the logarithmic gai reduce by db (egative lope) per factor of (decade) icreae i frequecy ω (db/decade). The traitio betwee the high ad low frequecy aymptote i foud by equatig the low ad high frequecy limit thi i kow a the corer or cut-off frequecy: ω c T Sice G jωt (jω) R( ω) + jx( ω) + jωt + ω T + ω T X R Phae: φ( ω) ta ta ( ωt) (rad) At low frequecie (ω ): φ() ta ( ) rad At high frequecie (ω ): φ( ) ta ( ) π / At corer frequecy (ω /T): φ(/ T) ta ( ) π / 4 rad rad 9

20 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 5 Aymptotic curve of G(jω) Aymptote Corer frequecy Aymptote Aymptotic curve of G(jω) Exact curve db -5 Exact curve φ -π/4 Aymptote /T./T /T /T ω -π/./t./t /T /T ω Exercie: Sketch the Bode plot for the followig trafer fuctio: G() 4 +

21 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3 Bode Diagram Phae (deg) Magitude (db) ω (rad/ec) 5.. Bode plot with MATLAB The MATLAB commad bode(sys) compute the logarithmic gai ad phae agle of the frequecy repoe of the LTI SYStf(um,de), where um ad de are the umerator ad deomiator coefficiet of the ytem, repectively. For example, to plot the Bode diagram how for the trafer fuctio of the previou exercie, we eter o the MATLAB commad lie: um [ ]; de [4 ]; SYS tf(um,de); bode(sys) or bode(um,de)

22 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 3 Bode Diagram 5 Phae (deg); Magitude (db) Frequecy (rad/ec) 6. Baic Fact About Egieerig Sytem ENGINEERING SYSTEMS OFTEN ARE REQUIRED TO OPERATE IN A STEADY OUTPUT CONDITION with the ytem deiged o the bai of achievig the bet output from the give raw material or power available. The deig of the ytem to produce the deired teady-tate i a problem i it ow right but ot the ubject matter of thi coure. OPERATING CONDITIONS ARE OFTEN CHANGED BY OPERATOR OR COMPUTER INTERVENTION uch a chage i power upplied by a power tatio to the atioal grid due to chage i the overall atioal power coumptio. ENGINEERING SYSTMES HAVE DYNAMICAL BEHAVIOURS they do ot jut produce the deired output temperature, preure, cocetratio, voltage, curret, frequecy, poitio, velocity, acceleratio, force, torque, flow, level, cocetratio, reactio rate,.. etc eve if that output variable i required to be cotat. Thi ca be due to the effect of diturbace (kow or ukow), phyical effect withi the ytem or huma itervetio/iterferece uch a chage required i operatio coditio. ENGINEERING SYSTEMS REQUIRE CONTROL to couteract the dyamic behaviour preferred by the ytem ad replace them by acceptable dyamic repoe, the proce mut be augmeted by a cotrol ytem icorporatig feature of output meauremet (eor), iput variable chage (actuator) ad a (dyamical) data proceig device that procee both eor data ad deired output pecificatio to geerate a deired plat iput igal to the actuator. Thi ytem almot alway ha a FEEDBACK tructure reflectig the fact that it output i ued to create the deired iput i real time. CONTROL SYSTEMS REQUIRE DESIGNING otwithtadig the impreio left by the popular TV ciece programme Tomorrow World, cotrol ytem are deiged rather tha beig available (without thought required from the uer) i a form that imply eed to be hooked

23 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) up to the plat. Eve thoe cotroller that are available commercially require a elemet of deig either off-lie or o-lie ad hece a effective egieer require a appreciatio of the deig proce ad deig tool curretly available. CONTROLLERS DEPEND ON THE DETAILED DYNAMICAL CHARACTERISTICS OF THE PROCESS TO BE CONTROLLED If that were ot the cae the it would oly be eceary to have oe cotrol ytem o ale. Practical experiece ha how that thi i ot feaible it i foud that it i eceary to have appropriate data ad ome udertadig of the proce to be cotrolled. Thi typically take the form of EITHER ) a mathematical model expreed i differetial equatio, trafer fuctio or tate pace form, AND/OR ) data o the behaviour of the plat output() i repoe to kow iput (uch a tep or iuoid) from which (a) the deired parameter for the model (obtaied i ()) ca be etimated or (b) if a phyical model i ot available a model ca be cotructed from a curve fittig or a, o-called, idetificatio procedure. THE DEGREE OF CONTROL OBTAINED DEPENDES ON THE AVAILABILITY OF SUITABLE MEASUREMENTS OF SYSTEM BEHAVIOUR ad the more accurate ad exteive thee meauremet are, the better cotrol will be. YOU CANNOT ALWAYS MEASURE WHAT YOU WANT TO MEASURE if you caot meaure the output variable of iteret (due to extreme phyical coditio of peed, temperature or preure.. etc, it i eceary to create a itelliget device which oberve the (available) meauremet ad ue them to create a ueful etimate of the (uavailable) output. A a example, how i it poible to cotrol the temperature i the cetre of a furace whe the temperature eor are placed o the exteral wall? CONTROLLERS CAN DO AMAZING THINS A cotrol requiremet are a varied a the applicatio ad eed for ew product, eve the virtually impoible ha bee aked for i the earch for the itelliget cotroller. Thi require the developmet of a abtract way of thikig but ha amazig coequece e.g. (a) the developmet of cotrol elemet capable of obervig ad accurately etimatig variable that caot be meaured (b) the developmet of cotrol ytem capable of adaptig to ew ituatio ad learig from experiece. 3

24 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) TUTORIAL PROBLEM SHEET 4. Idetify the gai ad time cotat of the followig firt-order trafer fuctio: 3 G() +. A RC circuit ha the followig trafer fuctio: Y() G() R() + 4 For a tep iput u(t) V for t : a) What i the teady-tate repoe of the circuit? b) What i the time take for the output of the RC circuit to reach 95% of it teady-tate repoe? c) Check your reult with MATLAB. 3. The geeral form of a firt-order ytem i decribed by: d T + ku(t) dt where T i the time cotat, k i the gai, u(t) ad are the iput ad output of the ytem repectively. The uit-tep repoe of a firt-order ytem i how i the followig figure. Determie the parameter k ad T from thi figure..5 Step Repoe Amplitude Time (ec.) 4

25 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 4. Coider the firt-order ytem, Y() R() T + Obtai the uit-tep repoe curve for T.,.5,., 5. ad. repectively, with MATLAB. 5. Coider the firt-order ytem, Y() k R() + Obtai the uit-tep repoe curve for k.,.5,., 5., ad. repectively, with MATLAB. 6. A geeral ecod-order ytem ha the form: Y() kω R() + ζω + ω What are the value of k, ζ, ad ω for the followig ytem: Y() 3 R() A ecod-order ytem i decribed by the differetial equatio: d d u(t) dt dt a) Write dow the trafer fuctio Y()/U() of the ytem, where U() ad Y() are the Laplace traform of u(t) ad, repectively. b) Obtai the dampig ratio ζ ad the atural frequecy ω of the ytem. c) Calculate the rie time ad percet overhoot of the ytem. d) Evaluate for a uit-tep iput u(t). e) Check your awer of the above with MATLAB. 8. For the cotrol ytem how by the block diagram, the umerical value of J kg-m ad B N-m/(rad/ec). R() + _ K + J + B _ / Y() K a) Fid the trafer fuctio Y()/R(). b) Determie the value of the gai K ad velocity feedback cotat K o that the maximum overhoot i the uit-tep repoe i. ad peak time i ec. c) With thee value of K ad K, obtai the rie time ad ettlig time. d) Obtai the repoe for the uit-tep iput r(t). e) Check the above calculatio with MATLAB. 5

26 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3) 9. Whe the ecod-order ytem Y() R() T K + + K i ubjected to a uit-tep iput, the ytem output repod a how i the followig figure. Determie K ad T from the repoe curve. Step Repoe.4.3. Amplitude Time (ec.). A iuoidal iput u(t) i(t) i applied to a ytem with trafer fuctio: Y() U() ( + ) Determie the teady-tate output, y (t), of the ytem.. The figure below how a block diagram of a pace vehicle attitude cotrol ytem where R ad Y are the Laplace traform of the referece (or deired) ad actual attitude agle repectively. Determie the value of K P ad K D to yield a ettlig time of.5 ecod ad % overhoot i the cloe-loop ytem for a uit-tep iput. R() + K P + Y() K D 6

27 ECM5 - Cotrol Egieerig Dr Mutafa M Aziz (3). Coider the ecod-order ytem Y() R() + ζ + Obtai the uit-impule repoe curve for ζ.,.3,.5,.7,., ad 4. repectively, with MATLAB. 3. Coider the ecod-order ytem Y() R() kω + ζω + ω Aumig that ω, k, obtai the uit-impule repoe curve for ζ.,.3,.5,.7,., ad 4. repectively, with MATLAB. 7

ELEC 372 LECTURE NOTES, WEEK 4 Dr. Amir G. Aghdam Concordia University

ELEC 372 LECTURE NOTES, WEEK 4 Dr. Amir G. Aghdam Concordia University ELEC 37 LECTURE NOTES, WEE 4 Dr Amir G Aghdam Cocordia Uiverity Part of thee ote are adapted from the material i the followig referece: Moder Cotrol Sytem by Richard C Dorf ad Robert H Bihop, Pretice Hall

More information

State space systems analysis

State space systems analysis State pace ytem aalyi Repreetatio of a ytem i tate-pace (tate-pace model of a ytem To itroduce the tate pace formalim let u tart with a eample i which the ytem i dicuio i a imple electrical circuit with

More information

Introduction to Control Systems

Introduction to Control Systems Itroductio to Cotrol Sytem CLASSIFICATION OF MATHEMATICAL MODELS Icreaig Eae of Aalyi Static Icreaig Realim Dyamic Determiitic Stochatic Lumped Parameter Ditributed Parameter Liear Noliear Cotat Coefficiet

More information

CONTROL SYSTEMS. Chapter 7 : Bode Plot. 40dB/dec 1.0. db/dec so resultant slope will be 20 db/dec and this is due to the factor s

CONTROL SYSTEMS. Chapter 7 : Bode Plot. 40dB/dec 1.0. db/dec so resultant slope will be 20 db/dec and this is due to the factor s CONTROL SYSTEMS Chapter 7 : Bode Plot GATE Objective & Numerical Type Solutio Quetio 6 [Practice Book] [GATE EE 999 IIT-Bombay : 5 Mark] The aymptotic Bode plot of the miimum phae ope-loop trafer fuctio

More information

Time Response. First Order Systems. Time Constant, T c We call 1/a the time constant of the response. Chapter 4 Time Response

Time Response. First Order Systems. Time Constant, T c We call 1/a the time constant of the response. Chapter 4 Time Response Time Repoe Chapter 4 Time Repoe Itroductio The output repoe of a ytem i the um of two repoe: the forced repoe ad the atural repoe. Although may techique, uch a olvig a differetial equatio or takig the

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 30 Sigal & Sytem Prof. Mark Fowler Note Set #8 C-T Sytem: Laplace Traform Solvig Differetial Equatio Readig Aigmet: Sectio 6.4 of Kame ad Heck / Coure Flow Diagram The arrow here how coceptual flow

More information

Automatic Control Systems

Automatic Control Systems Automatic Cotrol Sytem Lecture-5 Time Domai Aalyi of Orer Sytem Emam Fathy Departmet of Electrical a Cotrol Egieerig email: emfmz@yahoo.com Itrouctio Compare to the implicity of a firt-orer ytem, a eco-orer

More information

ECE 422 Power System Operations & Planning 6 Small Signal Stability. Spring 2015 Instructor: Kai Sun

ECE 422 Power System Operations & Planning 6 Small Signal Stability. Spring 2015 Instructor: Kai Sun ECE 4 Power Sytem Operatio & Plaig 6 Small Sigal Stability Sprig 15 Itructor: Kai Su 1 Referece Saadat Chapter 11.4 EPRI Tutorial Chapter 8 Power Ocillatio Kudur Chapter 1 Power Ocillatio The power ytem

More information

ECEN620: Network Theory Broadband Circuit Design Fall 2014

ECEN620: Network Theory Broadband Circuit Design Fall 2014 ECE60: etwork Theory Broadbad Circuit Deig Fall 04 Lecture 3: PLL Aalyi Sam Palermo Aalog & Mixed-Sigal Ceter Texa A&M Uiverity Ageda & Readig PLL Overview & Applicatio PLL Liear Model Phae & Frequecy

More information

The Performance of Feedback Control Systems

The Performance of Feedback Control Systems The Performace of Feedbac Cotrol Sytem Objective:. Secify the meaure of erformace time-domai the firt te i the deig roce Percet overhoot / Settlig time T / Time to rie / Steady-tate error e. ut igal uch

More information

CONTROL ENGINEERING LABORATORY

CONTROL ENGINEERING LABORATORY Uiverity of Techology Departmet of Electrical Egieerig Cotrol Egieerig Lab. CONTROL ENGINEERING LABORATORY By Dr. Abdul. Rh. Humed M.Sc. Quay Salim Tawfeeq M.Sc. Nihad Mohammed Amee M.Sc. Waleed H. Habeeb

More information

Lecture 30: Frequency Response of Second-Order Systems

Lecture 30: Frequency Response of Second-Order Systems Lecture 3: Frequecy Repoe of Secod-Order Sytem UHTXHQF\ 5HVSRQVH RI 6HFRQGUGHU 6\VWHPV A geeral ecod-order ytem ha a trafer fuctio of the form b + b + b H (. (9.4 a + a + a It ca be table, utable, caual

More information

System Control. Lesson #19a. BME 333 Biomedical Signals and Systems - J.Schesser

System Control. Lesson #19a. BME 333 Biomedical Signals and Systems - J.Schesser Sytem Cotrol Leo #9a 76 Sytem Cotrol Baic roblem Say you have a ytem which you ca ot alter but it repoe i ot optimal Example Motor cotrol for exokeleto Robotic cotrol roblem that ca occur Utable Traiet

More information

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,

More information

Fig. 1: Streamline coordinates

Fig. 1: Streamline coordinates 1 Equatio of Motio i Streamlie Coordiate Ai A. Soi, MIT 2.25 Advaced Fluid Mechaic Euler equatio expree the relatiohip betwee the velocity ad the preure field i ivicid flow. Writte i term of treamlie coordiate,

More information

Brief Review of Linear System Theory

Brief Review of Linear System Theory Brief Review of Liear Sytem heory he followig iformatio i typically covered i a coure o liear ytem theory. At ISU, EE 577 i oe uch coure ad i highly recommeded for power ytem egieerig tudet. We have developed

More information

Last time: Completed solution to the optimum linear filter in real-time operation

Last time: Completed solution to the optimum linear filter in real-time operation 6.3 tochatic Etimatio ad Cotrol, Fall 4 ecture at time: Completed olutio to the oimum liear filter i real-time operatio emi-free cofiguratio: t D( p) F( p) i( p) dte dp e π F( ) F( ) ( ) F( p) ( p) 4444443

More information

ELEC 372 LECTURE NOTES, WEEK 1 Dr. Amir G. Aghdam Concordia University

ELEC 372 LECTURE NOTES, WEEK 1 Dr. Amir G. Aghdam Concordia University EEC 37 ECTURE NOTES, WEEK Dr Amir G Aghdam Cocordia Uiverity Part of thee ote are adapted from the material i the followig referece: Moder Cotrol Sytem by Richard C Dorf ad Robert H Bihop, Pretice Hall

More information

u t u 0 ( 7) Intuitively, the maximum principles can be explained by the following observation. Recall

u t u 0 ( 7) Intuitively, the maximum principles can be explained by the following observation. Recall Oct. Heat Equatio M aximum priciple I thi lecture we will dicu the maximum priciple ad uiquee of olutio for the heat equatio.. Maximum priciple. The heat equatio alo ejoy maximum priciple a the Laplace

More information

EE 508 Lecture 6. Dead Networks Scaling, Normalization and Transformations

EE 508 Lecture 6. Dead Networks Scaling, Normalization and Transformations EE 508 Lecture 6 Dead Network Scalig, Normalizatio ad Traformatio Filter Cocept ad Termiology 2-d order polyomial characterizatio Biquadratic Factorizatio Op Amp Modelig Stability ad Itability Roll-off

More information

Adaptive control design for a Mimo chemical reactor

Adaptive control design for a Mimo chemical reactor Automatio, Cotrol ad Itelliget Sytem 013; 1(3): 64-70 Publihed olie July 10, 013 (http://www.ciecepublihiggroup.com/j/aci) doi: 10.11648/j.aci.0130103.15 Adaptive cotrol deig for a Mimo chemical reactor

More information

TESTS OF SIGNIFICANCE

TESTS OF SIGNIFICANCE TESTS OF SIGNIFICANCE Seema Jaggi I.A.S.R.I., Library Aveue, New Delhi eema@iari.re.i I applied ivetigatio, oe i ofte itereted i comparig ome characteritic (uch a the mea, the variace or a meaure of aociatio

More information

2.004 Dynamics and Control II Spring 2008

2.004 Dynamics and Control II Spring 2008 MIT OpeCourseWare http://ocw.mit.edu 2.004 Dyamics ad Cotrol II Sprig 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts Istitute of Techology

More information

Performance-Based Plastic Design (PBPD) Procedure

Performance-Based Plastic Design (PBPD) Procedure Performace-Baed Platic Deig (PBPD) Procedure 3. Geeral A outlie of the tep-by-tep, Performace-Baed Platic Deig (PBPD) procedure follow, with detail to be dicued i ubequet ectio i thi chapter ad theoretical

More information

High-Speed Serial Interface Circuits and Systems. Lect. 4 Phase-Locked Loop (PLL) Type 1 (Chap. 8 in Razavi)

High-Speed Serial Interface Circuits and Systems. Lect. 4 Phase-Locked Loop (PLL) Type 1 (Chap. 8 in Razavi) High-Speed Serial Iterface Circuit ad Sytem Lect. 4 Phae-Locked Loop (PLL) Type 1 (Chap. 8 i Razavi) PLL Phae lockig loop A (egative-feedback) cotrol ytem that geerate a output igal whoe phae (ad frequecy)

More information

EE 508 Lecture 6. Scaling, Normalization and Transformation

EE 508 Lecture 6. Scaling, Normalization and Transformation EE 508 Lecture 6 Scalig, Normalizatio ad Traformatio Review from Lat Time Dead Network X IN T X OUT T X OUT N T = D D The dead etwork of ay liear circuit i obtaied by ettig ALL idepedet ource to zero.

More information

Isolated Word Recogniser

Isolated Word Recogniser Lecture 5 Iolated Word Recogitio Hidde Markov Model of peech State traitio ad aligmet probabilitie Searchig all poible aligmet Dyamic Programmig Viterbi Aligmet Iolated Word Recogitio 8. Iolated Word Recogier

More information

THE CONCEPT OF THE ROOT LOCUS. H(s) THE CONCEPT OF THE ROOT LOCUS

THE CONCEPT OF THE ROOT LOCUS. H(s) THE CONCEPT OF THE ROOT LOCUS So far i the tudie of cotrol yte the role of the characteritic equatio polyoial i deteriig the behavior of the yte ha bee highlighted. The root of that polyoial are the pole of the cotrol yte, ad their

More information

LECTURE 13 SIMULTANEOUS EQUATIONS

LECTURE 13 SIMULTANEOUS EQUATIONS NOVEMBER 5, 26 Demad-upply ytem LETURE 3 SIMULTNEOUS EQUTIONS I thi lecture, we dicu edogeeity problem that arie due to imultaeity, i.e. the left-had ide variable ad ome of the right-had ide variable are

More information

Heat Equation: Maximum Principles

Heat Equation: Maximum Principles Heat Equatio: Maximum Priciple Nov. 9, 0 I thi lecture we will dicu the maximum priciple ad uiquee of olutio for the heat equatio.. Maximum priciple. The heat equatio alo ejoy maximum priciple a the Laplace

More information

x z Increasing the size of the sample increases the power (reduces the probability of a Type II error) when the significance level remains fixed.

x z Increasing the size of the sample increases the power (reduces the probability of a Type II error) when the significance level remains fixed. ] z-tet for the mea, μ If the P-value i a mall or maller tha a pecified value, the data are tatitically igificat at igificace level. Sigificace tet for the hypothei H 0: = 0 cocerig the ukow mea of a populatio

More information

Last time: Ground rules for filtering and control system design

Last time: Ground rules for filtering and control system design 6.3 Stochatic Etimatio ad Cotrol, Fall 004 Lecture 7 Lat time: Groud rule for filterig ad cotrol ytem deig Gral ytem Sytem parameter are cotaied i w( t ad w ( t. Deired output i grated by takig the igal

More information

Erick L. Oberstar Fall 2001 Project: Sidelobe Canceller & GSC 1. Advanced Digital Signal Processing Sidelobe Canceller (Beam Former)

Erick L. Oberstar Fall 2001 Project: Sidelobe Canceller & GSC 1. Advanced Digital Signal Processing Sidelobe Canceller (Beam Former) Erick L. Obertar Fall 001 Project: Sidelobe Caceller & GSC 1 Advaced Digital Sigal Proceig Sidelobe Caceller (Beam Former) Erick L. Obertar 001 Erick L. Obertar Fall 001 Project: Sidelobe Caceller & GSC

More information

Lecture 11. Course Review. (The Big Picture) G. Hovland Input-Output Limitations (Skogestad Ch. 3) Discrete. Time Domain

Lecture 11. Course Review. (The Big Picture) G. Hovland Input-Output Limitations (Skogestad Ch. 3) Discrete. Time Domain MER4 Advaced Cotrol Lecture Coure Review (he ig Picture MER4 ADVANCED CONROL EMEER, 4 G. Hovlad 4 Mai heme of MER4 Frequecy Domai Aalyi (Nie Chapter Phae ad Gai Margi Iput-Output Limitatio (kogetad Ch.

More information

Assignment 1 - Solutions. ECSE 420 Parallel Computing Fall November 2, 2014

Assignment 1 - Solutions. ECSE 420 Parallel Computing Fall November 2, 2014 Aigmet - Solutio ECSE 420 Parallel Computig Fall 204 ovember 2, 204. (%) Decribe briefly the followig term, expoe their caue, ad work-aroud the idutry ha udertake to overcome their coequece: (i) Memory

More information

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION I liear regreio, we coider the frequecy ditributio of oe variable (Y) at each of everal level of a ecod variable (X). Y i kow a the depedet variable.

More information

Dynamic Response of Linear Systems

Dynamic Response of Linear Systems Dyamic Respose of Liear Systems Liear System Respose Superpositio Priciple Resposes to Specific Iputs Dyamic Respose of st Order Systems Characteristic Equatio - Free Respose Stable st Order System Respose

More information

Gain-scheduling of Acceleration Estimator for Low-velocity Measurement with Encoders

Gain-scheduling of Acceleration Estimator for Low-velocity Measurement with Encoders Jue -5, INTEX, Gyeoggi-Do, orea Gai-chedulig of Acceleratio Etimator for Low-velocity Meauremet with Ecoder Seug-Woo So, Sag-Hu Lee ad Jog-Sug Hur Electro-Mechaical Reearch Ititute, Hyudai Heavy Idutrie

More information

Chapter (a) ζ. ω. 5 2 (a) Type 0 (b) Type 0 (c) Type 1 (d) Type 2 (e) Type 3 (f) Type 3. (g) type 2 (h) type (a) K G s.

Chapter (a) ζ. ω. 5 2 (a) Type 0 (b) Type 0 (c) Type 1 (d) Type 2 (e) Type 3 (f) Type 3. (g) type 2 (h) type (a) K G s. Chapter 5 5 1 (a) ζ. ω 0 707 rad / ec (b) 0 ζ 0. 707 ω rad / ec (c) ζ 0. 5 1 ω 5 rad / ec (d) 0. 5 ζ 0. 707 ω 0. 5 rad / ec 5 (a) Type 0 (b) Type 0 (c) Type 1 (d) Type (e) Type 3 (f) Type 3 (g) type (h)

More information

EE Control Systems

EE Control Systems Copyright FL Lewis 7 All rights reserved Updated: Moday, November 1, 7 EE 4314 - Cotrol Systems Bode Plot Performace Specificatios The Bode Plot was developed by Hedrik Wade Bode i 1938 while he worked

More information

Dr. Seeler Department of Mechanical Engineering Fall 2009 Lafayette College ME 479: Control Systems and Mechatronics Design and Analysis

Dr. Seeler Department of Mechanical Engineering Fall 2009 Lafayette College ME 479: Control Systems and Mechatronics Design and Analysis Dr. Seeler Departmet of Mechaical Egieerig Fall 009 Lafayette College ME 479: Cotrol Systems ad Mechatroics Desig ad Aalysis Lab 0: Review of the First ad Secod Order Step Resposes The followig remarks

More information

MAHALAKSHMI ENGINEERING COLLEGE TIRUCHIRAPALLI

MAHALAKSHMI ENGINEERING COLLEGE TIRUCHIRAPALLI MAHALASHMI ENGINEERING COLLEGE TIRUCHIRAPALLI-63. QUESTION BAN DEPARTMENT: ECE SUBJECT CODE: EC55 SEMESTER - III SUBJECT NAME: CONTROL SYSTEMS UNIT- SYSTEMS AND THEIR REPRESENTATION. Write Mao gai formula.

More information

1it is said to be overdamped. When 1, the roots of

1it is said to be overdamped. When 1, the roots of Homework 3 AERE573 Fall 8 Due /8(M) SOLUTIO PROBLEM (4pt) Coider a D order uderdamped tem trafer fuctio H( ) ratio The deomiator i the tem characteritic polomial P( ) (a)(5pt) Ue the quadratic formula,

More information

ME 410 MECHANICAL ENGINEERING SYSTEMS LABORATORY REGRESSION ANALYSIS

ME 410 MECHANICAL ENGINEERING SYSTEMS LABORATORY REGRESSION ANALYSIS ME 40 MECHANICAL ENGINEERING REGRESSION ANALYSIS Regreio problem deal with the relatiohip betwee the frequec ditributio of oe (depedet) variable ad aother (idepedet) variable() which i (are) held fied

More information

SOLUTION: The 95% confidence interval for the population mean µ is x ± t 0.025; 49

SOLUTION: The 95% confidence interval for the population mean µ is x ± t 0.025; 49 C22.0103 Sprig 2011 Homework 7 olutio 1. Baed o a ample of 50 x-value havig mea 35.36 ad tadard deviatio 4.26, fid a 95% cofidece iterval for the populatio mea. SOLUTION: The 95% cofidece iterval for the

More information

Professor: Mihnea UDREA DIGITAL SIGNAL PROCESSING. Grading: Web: MOODLE. 1. Introduction. General information

Professor: Mihnea UDREA DIGITAL SIGNAL PROCESSING. Grading: Web:   MOODLE. 1. Introduction. General information Geeral iformatio DIGITL SIGL PROCESSIG Profeor: ihea UDRE B29 mihea@comm.pub.ro Gradig: Laboratory: 5% Proect: 5% Tet: 2% ial exam : 5% Coure quiz: ±% Web: www.electroica.pub.ro OODLE 2 alog igal proceig

More information

Statistical Inference Procedures

Statistical Inference Procedures Statitical Iferece Procedure Cofidece Iterval Hypothei Tet Statitical iferece produce awer to pecific quetio about the populatio of iteret baed o the iformatio i a ample. Iferece procedure mut iclude a

More information

Appendix: The Laplace Transform

Appendix: The Laplace Transform Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio

More information

University of California at Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences

University of California at Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences A Uiversity of Califoria at Berkeley College of Egieerig Departmet of Electrical Egieerig ad Computer Scieces U N I V E R S T H E I T Y O F LE T TH E R E B E LI G H T C A L I F O R N 8 6 8 I A EECS : Sigals

More information

Statistics and Chemical Measurements: Quantifying Uncertainty. Normal or Gaussian Distribution The Bell Curve

Statistics and Chemical Measurements: Quantifying Uncertainty. Normal or Gaussian Distribution The Bell Curve Statitic ad Chemical Meauremet: Quatifyig Ucertaity The bottom lie: Do we trut our reult? Should we (or ayoe ele)? Why? What i Quality Aurace? What i Quality Cotrol? Normal or Gauia Ditributio The Bell

More information

Comments on Discussion Sheet 18 and Worksheet 18 ( ) An Introduction to Hypothesis Testing

Comments on Discussion Sheet 18 and Worksheet 18 ( ) An Introduction to Hypothesis Testing Commet o Dicuio Sheet 18 ad Workheet 18 ( 9.5-9.7) A Itroductio to Hypothei Tetig Dicuio Sheet 18 A Itroductio to Hypothei Tetig We have tudied cofidece iterval for a while ow. Thee are method that allow

More information

Capacitors and PN Junctions. Lecture 8: Prof. Niknejad. Department of EECS University of California, Berkeley. EECS 105 Fall 2003, Lecture 8

Capacitors and PN Junctions. Lecture 8: Prof. Niknejad. Department of EECS University of California, Berkeley. EECS 105 Fall 2003, Lecture 8 CS 15 Fall 23, Lecture 8 Lecture 8: Capacitor ad PN Juctio Prof. Nikejad Lecture Outlie Review of lectrotatic IC MIM Capacitor No-Liear Capacitor PN Juctio Thermal quilibrium lectrotatic Review 1 lectric

More information

ECE 522 Power Systems Analysis II 3.2 Small Signal Stability

ECE 522 Power Systems Analysis II 3.2 Small Signal Stability ECE 5 Power Sytem Aalyi II 3. Small Sigal Stability Sprig 18 Itructor: Kai Su 1 Cotet 3..1 Small igal tability overview ad aalyi method 3.. Small igal tability ehacemet Referece: EPRI Dyamic Tutorial Chapter

More information

Chapter 2 Feedback Control Theory Continued

Chapter 2 Feedback Control Theory Continued Chapter Feedback Cotrol Theor Cotiued. Itroductio I the previous chapter, the respose characteristic of simple first ad secod order trasfer fuctios were studied. It was show that first order trasfer fuctio,

More information

Control of a Linear Permanent Magnet Synchronous Motor using Multiple Reference Frame Theory

Control of a Linear Permanent Magnet Synchronous Motor using Multiple Reference Frame Theory Cotrol of a Liear Permaet Maget Sychroou Motor uig Multiple Referece Frame Theory Jawad Faiz ad B. RezaeiAlam Departmet of Electrical ad Computer Egieerig Faculty of Egieerig Uiverity of Tehra, Tehra,

More information

MEM 255 Introduction to Control Systems: Analyzing Dynamic Response

MEM 255 Introduction to Control Systems: Analyzing Dynamic Response MEM 55 Itroductio to Cotrol Systems: Aalyzig Dyamic Respose Harry G. Kwaty Departmet of Mechaical Egieerig & Mechaics Drexel Uiversity Outlie Time domai ad frequecy domai A secod order system Via partial

More information

Sinusoidal Steady-state Analysis

Sinusoidal Steady-state Analysis Siusoidal Steady-state Aalysis Complex umber reviews Phasors ad ordiary differetial equatios Complete respose ad siusoidal steady-state respose Cocepts of impedace ad admittace Siusoidal steady-state aalysis

More information

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B. Amme 3500 : System Dyamics ad Cotrol Root Locus Course Outlie Week Date Cotet Assigmet Notes Mar Itroductio 8 Mar Frequecy Domai Modellig 3 5 Mar Trasiet Performace ad the s-plae 4 Mar Block Diagrams Assig

More information

a 1 = 1 a a a a n n s f() s = Σ log a 1 + a a n log n sup log a n+1 + a n+2 + a n+3 log n sup () s = an /n s s = + t i

a 1 = 1 a a a a n n s f() s = Σ log a 1 + a a n log n sup log a n+1 + a n+2 + a n+3 log n sup () s = an /n s s = + t i 0 Dirichlet Serie & Logarithmic Power Serie. Defiitio & Theorem Defiitio.. (Ordiary Dirichlet Serie) Whe,a,,3, are complex umber, we call the followig Ordiary Dirichlet Serie. f() a a a a 3 3 a 4 4 Note

More information

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 6 - Small Signal Stability

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 6 - Small Signal Stability ECE 4/5 Power System Operatios & Plaig/Power Systems Aalysis II : 6 - Small Sigal Stability Sprig 014 Istructor: Kai Su 1 Refereces Kudur s Chapter 1 Saadat s Chapter 11.4 EPRI Tutorial s Chapter 8 Power

More information

CHAPTER 6. Confidence Intervals. 6.1 (a) y = 1269; s = 145; n = 8. The standard error of the mean is = s n = = 51.3 ng/gm.

CHAPTER 6. Confidence Intervals. 6.1 (a) y = 1269; s = 145; n = 8. The standard error of the mean is = s n = = 51.3 ng/gm. } CHAPTER 6 Cofidece Iterval 6.1 (a) y = 1269; = 145; = 8. The tadard error of the mea i SE ȳ = = 145 8 = 51.3 g/gm. (b) y = 1269; = 145; = 30. The tadard error of the mea i ȳ = 145 = 26.5 g/gm. 30 6.2

More information

Introduction to Wavelets and Their Applications

Introduction to Wavelets and Their Applications Itroductio to Wavelet ad Their Applicatio a T. Bialaiewicz Uiverity of Colorado at Dever ad Health Sciece Ceter Wavelet are a powerful tool for tudyig time-frequecy behavior of fiite-eergy igal. Advatage

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

Chapter #5 EEE Control Systems

Chapter #5 EEE Control Systems Sprig EEE Chpter #5 EEE Cotrol Sytem Deig Bed o Root Locu Chpter / Sprig EEE Deig Bed Root Locu Led Cotrol (equivlet to PD cotrol) Ued whe the tedy tte propertie of the ytem re ok but there i poor performce,

More information

Chapter 9. Key Ideas Hypothesis Test (Two Populations)

Chapter 9. Key Ideas Hypothesis Test (Two Populations) Chapter 9 Key Idea Hypothei Tet (Two Populatio) Sectio 9-: Overview I Chapter 8, dicuio cetered aroud hypothei tet for the proportio, mea, ad tadard deviatio/variace of a igle populatio. However, ofte

More information

Chapter 7, Solution 1C.

Chapter 7, Solution 1C. hapter 7, Solutio 1. he velocity of the fluid relative to the immered olid body ufficietly far away from a body i called the free-tream velocity,. he uptream or approach velocity i the velocity of the

More information

Chapter 4 : Laplace Transform

Chapter 4 : Laplace Transform 4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic

More information

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to: 2.003 Egieerig Dyamics Problem Set 9--Solutio Problem 1 Fid the equatio of motio for the system show with respect to: a) Zero sprig force positio. Draw the appropriate free body diagram. b) Static equilibrium

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Hidden Markov Model Parameters

Hidden Markov Model Parameters .PPT 5/04/00 Lecture 6 HMM Traiig Traiig Hidde Markov Model Iitial model etimate Viterbi traiig Baum-Welch traiig 8.7.PPT 5/04/00 8.8 Hidde Markov Model Parameter c c c 3 a a a 3 t t t 3 c a t A Hidde

More information

Introduction to Signals and Systems, Part V: Lecture Summary

Introduction to Signals and Systems, Part V: Lecture Summary EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive

More information

We will look for series solutions to (1) around (at most) regular singular points, which without

We will look for series solutions to (1) around (at most) regular singular points, which without ENM 511 J. L. Baai April, 1 Frobeiu Solutio to a d order ODE ear a regular igular poit Coider the ODE y 16 + P16 y 16 + Q1616 y (1) We will look for erie olutio to (1) aroud (at mot) regular igular poit,

More information

Class 07 Time domain analysis Part II 2 nd order systems

Class 07 Time domain analysis Part II 2 nd order systems Class 07 Time domai aalysis Part II d order systems Time domai aalysis d order systems iput S output Secod order systems of the type α G(s) as + bs + c Time domai aalysis d order systems iput S α as +

More information

100(1 α)% confidence interval: ( x z ( sample size needed to construct a 100(1 α)% confidence interval with a margin of error of w:

100(1 α)% confidence interval: ( x z ( sample size needed to construct a 100(1 α)% confidence interval with a margin of error of w: Stat 400, ectio 7. Large Sample Cofidece Iterval ote by Tim Pilachowki a Large-Sample Two-ided Cofidece Iterval for a Populatio Mea ectio 7.1 redux The poit etimate for a populatio mea µ will be a ample

More information

EE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis

EE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis EE/ME/AE34: Dynamical Sytem Chapter 8: Tranfer Function Analyi The Sytem Tranfer Function Conider the ytem decribed by the nth-order I/O eqn.: ( n) ( n 1) ( m) y + a y + + a y = b u + + bu n 1 0 m 0 Taking

More information

STABILITY OF THE ACTIVE VIBRATION CONTROL OF CANTILEVER BEAMS

STABILITY OF THE ACTIVE VIBRATION CONTROL OF CANTILEVER BEAMS Iteratioal Coferece o Vibratio Problem September 9-,, Liboa, Portugal STBILITY OF THE CTIVE VIBRTIO COTROL OF CTILEVER BEMS J. Tůma, P. Šuráe, M. Mahdal VSB Techical Uierity of Otraa Czech Republic Outlie.

More information

20. CONFIDENCE INTERVALS FOR THE MEAN, UNKNOWN VARIANCE

20. CONFIDENCE INTERVALS FOR THE MEAN, UNKNOWN VARIANCE 20. CONFIDENCE INTERVALS FOR THE MEAN, UNKNOWN VARIANCE If the populatio tadard deviatio σ i ukow, a it uually will be i practice, we will have to etimate it by the ample tadard deviatio. Sice σ i ukow,

More information

Weak formulation and Lagrange equations of motion

Weak formulation and Lagrange equations of motion Chapter 4 Weak formulatio ad Lagrage equatio of motio A mot commo approach to tudy tructural dyamic i the ue of the Lagrage equatio of motio. Thee are obtaied i thi chapter tartig from the Cauchy equatio

More information

18.05 Problem Set 9, Spring 2014 Solutions

18.05 Problem Set 9, Spring 2014 Solutions 18.05 Problem Set 9, Sprig 2014 Solutio Problem 1. (10 pt.) (a) We have x biomial(, θ), o E(X) =θ ad Var(X) = θ(1 θ). The rule-of-thumb variace i jut 4. So the ditributio beig plotted are biomial(250,

More information

STRONG DEVIATION THEOREMS FOR THE SEQUENCE OF CONTINUOUS RANDOM VARIABLES AND THE APPROACH OF LAPLACE TRANSFORM

STRONG DEVIATION THEOREMS FOR THE SEQUENCE OF CONTINUOUS RANDOM VARIABLES AND THE APPROACH OF LAPLACE TRANSFORM Joural of Statitic: Advace i Theory ad Applicatio Volume, Number, 9, Page 35-47 STRONG DEVIATION THEORES FOR THE SEQUENCE OF CONTINUOUS RANDO VARIABLES AND THE APPROACH OF LAPLACE TRANSFOR School of athematic

More information

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web:     Ph: Serial :. PT_EE_A+C_Control Sytem_798 Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubanewar olkata Patna Web: E-mail: info@madeeay.in Ph: -4546 CLASS TEST 8-9 ELECTRICAL ENGINEERING Subject

More information

Dorf, R.C., Wan, Z., Johnson, D.E. Laplace Transform The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000

Dorf, R.C., Wan, Z., Johnson, D.E. Laplace Transform The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000 Dorf, R.C., Wa, Z., Joho, D.E. Laplace Traform The Electrical Egieerig Hadbook Ed. Richard C. Dorf Boca Rato: CRC Pre LLC, 6 Laplace Traform Richard C. Dorf Uiverity of Califoria, Davi Zhe Wa Uiverity

More information

EULER-MACLAURIN SUM FORMULA AND ITS GENERALIZATIONS AND APPLICATIONS

EULER-MACLAURIN SUM FORMULA AND ITS GENERALIZATIONS AND APPLICATIONS EULER-MACLAURI SUM FORMULA AD ITS GEERALIZATIOS AD APPLICATIOS Joe Javier Garcia Moreta Graduate tudet of Phyic at the UPV/EHU (Uiverity of Baque coutry) I Solid State Phyic Addre: Practicate Ada y Grijalba

More information

Explicit scheme. Fully implicit scheme Notes. Fully implicit scheme Notes. Fully implicit scheme Notes. Notes

Explicit scheme. Fully implicit scheme Notes. Fully implicit scheme Notes. Fully implicit scheme Notes. Notes Explicit cheme So far coidered a fully explicit cheme to umerically olve the diffuio equatio: T + = ( )T + (T+ + T ) () with = κ ( x) Oly table for < / Thi cheme i ometime referred to a FTCS (forward time

More information

STUDY OF SUBSYNCHRONOUS RESONANCE AND ANALYSIS OF SSR

STUDY OF SUBSYNCHRONOUS RESONANCE AND ANALYSIS OF SSR ISSN: 39-8753 Vol., Iue 7, July 03 STUDY OF SUBSYNCHRONOUS RESONANCE AND ANALYSIS OF SSR Rajit M. Zede, Suryakat H. Pawar P.. Studet, Departmet of Electrical Egieerig, ovt. College of Egg.Karad, Maharahtra,

More information

Dynamic System Response

Dynamic System Response Solutio of Liear, Costat-Coefficiet, Ordiary Differetial Equatios Classical Operator Method Laplace Trasform Method Laplace Trasform Properties 1 st -Order Dyamic System Time ad Frequecy Respose d -Order

More information

STUDENT S t-distribution AND CONFIDENCE INTERVALS OF THE MEAN ( )

STUDENT S t-distribution AND CONFIDENCE INTERVALS OF THE MEAN ( ) STUDENT S t-distribution AND CONFIDENCE INTERVALS OF THE MEAN Suppoe that we have a ample of meaured value x1, x, x3,, x of a igle uow quatity. Aumig that the meauremet are draw from a ormal ditributio

More information

ME 375 FINAL EXAM Wednesday, May 6, 2009

ME 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 information

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004 METR4200 Advanced Control Lecture 4 Chapter Nie Controller Deign via Frequency Repone G. Hovland 2004 Deign Goal Tranient repone via imple gain adjutment Cacade compenator to improve teady-tate error Cacade

More information

Numerical Solution of Coupled System of Nonlinear Partial Differential Equations Using Laplace-Adomian Decomposition Method

Numerical Solution of Coupled System of Nonlinear Partial Differential Equations Using Laplace-Adomian Decomposition Method I S S N 3 4 7-9 V o l u m e N u m b e r 0 8 J o u r a l o f A d v a c e i M a t h e m a t i c Numerical Solutio of Coupled Sytem of Noliear Partial Differetial Equatio Uig Laplace-Adomia Decompoitio Method

More information

Chapter 5. Root Locus Techniques

Chapter 5. Root Locus Techniques Chapter 5 Rt Lcu Techique Itrducti Sytem perfrmace ad tability dt determied dby cled-lp l ple Typical cled-lp feedback ctrl ytem G Ope-lp TF KG H Zer -, - Ple 0, -, - K Lcati f ple eaily fud Variati f

More information

ON THE SCALE PARAMETER OF EXPONENTIAL DISTRIBUTION

ON THE SCALE PARAMETER OF EXPONENTIAL DISTRIBUTION Review of the Air Force Academy No. (34)/7 ON THE SCALE PARAMETER OF EXPONENTIAL DISTRIBUTION Aca Ileaa LUPAŞ Military Techical Academy, Bucharet, Romaia (lua_a@yahoo.com) DOI:.96/84-938.7.5..6 Abtract:

More information

Digital Signal Processing, Fall 2010

Digital Signal Processing, Fall 2010 Digital Sigal Proeig, Fall 2 Leture 3: Samplig ad reotrutio, traform aalyi of LTI ytem tem Zheg-ua Ta Departmet of Eletroi Sytem Aalborg Uiverity, Demar t@e.aau.d Coure at a glae MM Direte-time igal ad

More information

8.6 Order-Recursive LS s[n]

8.6 Order-Recursive LS s[n] 8.6 Order-Recurive LS [] Motivate ti idea wit Curve Fittig Give data: 0,,,..., - [0], [],..., [-] Wat to fit a polyomial to data.., but wic oe i te rigt model?! Cotat! Quadratic! Liear! Cubic, Etc. ry

More information

EE 205 Dr. A. Zidouri. Electric Circuits II. Frequency Selective Circuits (Filters) Bode Plots: Complex Poles and Zeros.

EE 205 Dr. A. Zidouri. Electric Circuits II. Frequency Selective Circuits (Filters) Bode Plots: Complex Poles and Zeros. EE 5 Dr. A. Zidouri Electric Circuits II Frequecy Selective Circuits (Filters) Bode Plots: Complex Poles ad Zeros Lecture #4 - - EE 5 Dr. A. Zidouri The material to be covered i this lecture is as follows:

More information

Analysis of Analytical and Numerical Methods of Epidemic Models

Analysis of Analytical and Numerical Methods of Epidemic Models Iteratioal Joural of Egieerig Reearc ad Geeral Sciece Volue, Iue, Noveber-Deceber, 05 ISSN 09-70 Aalyi of Aalytical ad Nuerical Metod of Epideic Model Pooa Kuari Aitat Profeor, Departet of Mateatic Magad

More information

Zeta-reciprocal Extended reciprocal zeta function and an alternate formulation of the Riemann hypothesis By M. Aslam Chaudhry

Zeta-reciprocal Extended reciprocal zeta function and an alternate formulation of the Riemann hypothesis By M. Aslam Chaudhry Zeta-reciprocal Eteded reciprocal zeta fuctio ad a alterate formulatio of the Riema hypothei By. Alam Chaudhry Departmet of athematical Sciece, Kig Fahd Uiverity of Petroleum ad ieral Dhahra 36, Saudi

More information

ME 375 FINAL EXAM Friday, May 6, 2005

ME 375 FINAL EXAM Friday, May 6, 2005 ME 375 FINAL EXAM Friay, May 6, 005 Divisio: Kig 11:30 / Cuigham :30 (circle oe) Name: Istructios (1) This is a close book examiatio, but you are allowe three 8.5 11 crib sheets. () You have two hours

More information

Applied Mathematical Sciences, Vol. 9, 2015, no. 3, HIKARI Ltd,

Applied Mathematical Sciences, Vol. 9, 2015, no. 3, HIKARI Ltd, Applied Mathematical Sciece Vol 9 5 o 3 7 - HIKARI Ltd wwwm-hiaricom http://dxdoiorg/988/am54884 O Poitive Defiite Solutio of the Noliear Matrix Equatio * A A I Saa'a A Zarea* Mathematical Sciece Departmet

More information

Statistics Problem Set - modified July 25, _. d Q w. i n

Statistics Problem Set - modified July 25, _. d Q w. i n Statitic Problem Set - modified July 5, 04 x i x i i x i _ x x _ t d Q w F x x t pooled calculated pooled. f d x x t calculated / /.. f d Kow cocept of Gauia Curve Sytematic Error Idetermiate Error t-tet

More information