Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
|
|
- Tracey Hill
- 6 years ago
- Views:
Transcription
1 Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system is a set of eqations (differential eqations) that represents the dynamics of the system sing physics laws. The model permits to stdy system transients and steady state performance. Model compleity As model becomes more detailed it also can become more accrate. Model accracy needed for control system design is normally simpler than the model sed for system simlation. Simpler models: -ignore some physical phenomena -approimate linearly nonlinear characteristics -se lmped parameters approimation of distribted parameters systems. For the design of a control system: -an initial simplified model is bilt for conceptal design -a more model is sed for controller design and parameters determination 3-
2 Linear Systems For linear systems the principle of sperposition is valid and the response to a comple inpt can be calclated by smming p the responses to its components. Linear Time Invariant (LTI ) Systems verss Linear Time Varying Systems Linear Time Invariant (LTI ) Systems = systems: - represented by lmped components- - described by linear differential eqations - parameters of the eqations are time invariant. Systems with parameters that vary in time are called linear time varying systems. Eamples: a car in motion or a rocket in flight have weight that diminishes as the fel is consmed. 3-
3 3. Transfer Fnction L{otpt} Transfer Fnction = G(s) = L{inpt} for zero initial conditions. Inpt (s) Transfer Fnction=G(s) Otpt y(s) G (s) = y(s) (s) The transfer fnction of a system represents the link between the inpt to the system to the otpt of the system. The transfer fnction of a system G(s) is a comple fnction that describes system dynamics in s-domains opposed t the differential eqations that describe system dynamics in time domain. The transfer fnction is independent of the inpt to the system and does not provide any information concerning the internal strctre of the system. Same transfer fnction can represent different systems. The transfer fnction permits to calclate the otpt or response for varios inpts. The transfer fnction can be calclated analytically starting from the physics eqations or can be determined eperimentally by measring the otpt to varios known inpts to the system. 3-3
4 Eample: Car sspension model m d Y o /dt m m k b Y o Y i k(y o -Y i ) B (dy o /dt - dy i /dt) d Y m dt o dyodt dyidt b dt k(y o Y ) = 0 i Laplace transform for zero initial conditions gives ms Y (s) bs(y (s) Y (s)) k(y (s) Y (s)) = 0 or o o ms Y (s) bsy (s) ky (s) = bsy (s) ky (s) o The transfer fnction is o i o i o i i Y o (s) Y (s) i = bs k ms bs k 3-4
5 Implse Response The Laplace transform of an implse fnction δ(t) is given by L { δ (t)} = The otpt of a system de to an implse inpt (s)= δ(s) = is y (s) = G(s) (s) = G(s) The implse response of a system is identical to the transfer fnction of that system. The inverse Laplace transform of the implse response G(s) L {G(s)} = g(t) The transfer fnction that contains complete information abot the dynamic characteristics of a system can be obtained by applying an implse inpt (t)= δ(t) and measring system response y(t) which in this case is identical to g(t). The transfer fnction will then be G(s)=L{g(t)} 3-5
6 3.3 Block Diagrams Block diagrams are a graphical representation of the system model. The blocks represent physical or fnctional components of the system. Each block has inscribed the transfer fnction of that component the relate the otpt of the component to its inpt. Block Diagrams Block diagrams consist of. blocks. smmation jnctions 3. paths 4. branching points. Block Eample: car sspension system Y i (s) bsk ms bsk Y o (s) Blocks represent physical or fnctional components in the system. In the block is inscribed the transfer fnction of that component of the system 3-6
7 . Smmation Jnction B A C - D=AB-C A B Σ D=AB-C -C A smming jnction reslts in the addition or sbtraction of inpt signals for a single otpt. 3. Path X(s) Signal X(s) flows along the directed path 4. Branching Point X(s) X(s) X(s) At the branching point a signal splits into two signals of the same vale. 3-7
8 3.4 Block Diagram of a Closed Loop System In a closed loop control system also called feedback control system the otpt variable y(s) is measred as y m (s) sbtracted from its desired vale y d (s) to calclate the error e(s)= y d (s) - y m (s) y d (s) _ e (s)= y d (s)- y m (s) G(s) y y m(s) H(s) G(s) is feed forward transfer fnction (of the controller actator and the system) G(s) = y(s) / e(s) H(s) is feedback transfer fnction (of the sensor) H(s) = y m (s) / y(s) G(s) H(s) = open loop transfer fnction G(s) H(s) = [y m (s) / y(s)][ y(s) / e(s)]= y m (s) / e(s) y d (s) / y(s) is closed loop transfer fnction obtained from the above eqations by eliminating e(s) and y m (s) y(s) = G(s) e(s) = G(s) [y d (s) - y m (s)]= G(s) [y d (s) H(s)y (s)] y(s) G(s) H(s) y(s) =G(s) y d (s) y(s)/y d (s)= G(s) / [ G(s) H(s)] This eqation gives the single block eqivalent of the above closed loop system y d (s) G(s) G(s)H(s) y(s) 3-8
9 The transfer fnction represents the closed loop system dynamics with comple fnctions. The otpt y(s) is given by y(s)= G(s) / [ G(s) H(s)] y d (s) and depends on the closed loop transfer fnction and the desired vale of the otpt y d (s) called also the inpt (to the closed loop system). The following positive feedback block diagram y d (s) G(s) y (s) y m(s) H(s) is eqivalent to a negative feedback one if H(s) is replaced by H(s) y d (s) - G(s) y (s) y m(s) -H(s) which is eqivalent to the block y d (s) G(s) G(s)H(s) y(s) 3-9
10 Closed loop Control System Advantage Ideally i.e it is reqired that y(s)/y d (s) G(s) / [ G(s) H(s)] or / [/ G(s) H(s)] which is achieved for a very high vale of G(s) G(s) G(s) represents the controller actator and the system and G(s) and the high feed forward gain can be achieved by a very high vale of the controller transfer fnction. In the stability stdy it will be seen that there is a limit to sch high vale de to stability constraints. 3-0
11 Error Error of a closed loop control system can be obtained from: e(s)=y d (s)-y m (s)= y d (s)-h(s) y (s)=y d (s){- H(s)y(s)/y d (s)}= y d (s) {- H(s)G(s) / [ G(s) H(s)] }= {- H(s)G(s) / [ G(s) H(s)] }y d (s)= { / [ G(s) H(s)] }y d (s) y d (s) G(s)H(s) e(s) Steady State Error Final vale theorem gives the steady state error of a system e( ) i.e. error when t By taking the limit for s 0 of e( )=lim t e(t)= lim s 0 {se(s)}= lim s 0 { / [ G(s) H(s)] }y d (s)= { / [ G(0) H(0)] }y d (0) The Components of Closed Loop Control System -The system to be controlled S (s) -The transdcer (sensor) H(s) to measre the otpts y(s) and feed them back as y m (s) -the comparator to calclate the error e(s) between the desired vale of the otpt y d (s) and otpt measrement y m (s) -the controller C(s) that ses this error signal e\(s) and generates the command c (s) to the system -the actator M(s) sch that 3-
12 G(s) = C(s) M(s) S(s) y d (s) _ e (s) C(s) CONTROLLER c power spply M(s) (s) ACTUATOR S(s) SYSTEM y (s) y m(s) H(s) Pertrbations Pertrbations p (s) to a control system generally act on the system G(s) and occr as an additional power inpt that adds to the power inpt (s) to the system. y d (s) _ e (s) C(s) c power spply M(s) (s) p(s) S(s) y (s) y m(s) H(s) Sperposition principle is applied for linear systems and consider each inpt independently and the otpts corresponding to each inpt alone can be added to give the total otpt. The principle of sperposition is applied as follows -the otpt of the system y (s) is calclated as a reslt of the inpt command y d (s) for p(s)=0. 3-
13 y G(s) (s) = G(s)H(s) y (s) d -the otpt of the system y (s) is calclated as a reslt of the p(s) for y d (s)=0. y (s)=s(s){c(s)/[c(s)m(s)s(s)h(s)]}p(s) or y C(s) (s) = p(s) G(s)H(s) - add the two otpts y(s)= y (s) y (s) (s) = y (s) y (s) G(s) = G(s)H(s) y y (s) d C(s) G(s)y (s) C(s)p(s) p(s) = d G(s)H(s) G(s)H(s) The effect of the pertrbation p(s) is cancelled for G(s)H(s)>> and M(s) >>. For G(s)H(s)>> G(s)H(s) G(s)H(s) and C(s)/[G(s)H(s)]= C(s)/[C(s)M(s)S(s)H(s)] C(s)/[C(s)M(s)S(s)H(s)] /[ M(s)S(s)H(s)] sch that for y(s) = y (s) y(s) = y (s) p(s) H(s) d M(s)S(s)H(s) For M(s) >> /[ M(s)S(s)H(s)] 0 and y(s) = y (s) y(s) = y (s) H(s) d i.e. the otpt y(s) is not affected by p(s) and G(s). H(s)= will reslt in the ideal y(s)= y d (s) 3-3
14 Block Diagram Redction Block diagram for actal systems can contain a large nmber of blocks and block diagram redction is reqired to redce it to a single eqivalent block sing the following two rles:. The prodct of the transfer fnctions in the feedforward direction shold remain the same. The prodct of the transfer fnctions arond the loop shold remain the same Eample: G (s) - G (s) H (s) G 3 (s) H (s) The top negative feedback loop cannot be redced nless the smming point is moved ahead G (s). - G (s) G (s) a(s) G 3 (s) H (s) We will determine the new nknown feedback transfer fnction a(s) sing the above two rles.. The first rle is actally satisfied as both block diagrams have the same prodct of the transfer fnctions in the feedforward direction G (s) G (s) G 3 (s) 3-4
15 . The prodct of the transfer fnctions arond the loop for the initial block diagram is G (s) [G (s) G 3 (s)/( G (s)g 3 (s) H (s)]= G (s) G (s) G 3 (s)/[( G (s)g 3 (s) H (s)] and for the modified block diagram G (s) G (s) G 3 (s)/( G (s)g (s)g 3 (s) a(s)] In order to be the same given the same nmerators the denominators have to be eqal G (s)g 3 (s) H (s)= G (s)g (s)g 3 (s) a(s) or G (s)g 3 (s) H (s)= G (s)g (s)g 3 (s) a(s) or H (s)= G (s) a(s) which gives the soltion for the nknown new feedback transfer fnction a(s) = H (s)/g (s) Previos block diagram is eqivalent to H (s)/g (s) - G (s) G (s) G 3 (s) H (s) as the switch of the smming points satisfy the two rles. Now the bottom positive feedback loop can be redced to G (s)g (s)/[-g (s)g (s)h (s)] sch that 3-5
16 - G (s)g (s)/[-g (s)g (s)h (s)] H (s)/g (s) G 3 (s) and final redction gives the final transfer fnction G (s)g (s) G 3 (s)/[-g (s)g (s)h (s)]/{ G (s)g (s) G 3 (s)/[- G (s)g (s)h (s) H (s)/g (s)] (See Eample 3- and A-3- to A-3-5) 3-6
17 Controllers Controllers are the part of the overall feedback control system that are the main focs of control engineering. desired otpt or error C(s) CONTROLLER otpt measrement controller otpt H(s) power spply M(s) ACTUATOR inpt S(s) SYSTEM otpt of the system y d (s) _ e (s) C(s) power spply c (s) M(s) (s) p(s) S(s) y (s) y m(s) H(s) Classifications of Controllers By the type of implementation: -pnematic controller -hydralic controller -analog electronic controllers -digital controllers Modern implementation is as digital controllers. By control law: On-off Proportional P-control PD-control PID-control In Control I and II only the last three are stdied. 3-7
18 P-control The transfer fnction is c (s) = k e(s) p where k p is the proportional gain The block diagram is e(s) k p c (s) The P-controller corresponds to an operational amplifier with an adjstable gain k p. PD Control The transfer fnction is c (s) e(s) = k p k s d = k p ( T s) d where k d is the derivative gain and T d = k d / k p The block diagram is e(s) k p k p s c (s) 3-8
19 PID Control The PDI controller transfer fnction of is c (s) e(s) = k p k s d k i s = k p ( T s d ) T s i where k i is the integral gain and /T i = k i / k p The block diagram is e(s) k p k p sk i /s c (s) PID controllers are freqently sed in applications. 3-9
20 3.5 State Space Representation State space based control approaches are developed in Rssia in 960s. This is called Modern Control Theory or State space Control Theory or Time Domain Approach (as opposed to Conventional Control Theory of Classic Control Theory or Freqency Domain Approach The major limitation of the conventional freqency domain models is the reqirement that the system be linear and Single Inpt Single Otpt (SISO). Definitions State Variables: are the minimm set of variables that niqely define the state of a dynamic system at any instant of time. States: are the vales of the set of state variables. For known crrent state vales and of the inpt the state of the system in the ftre time can be calclated. State Vector: is the vector of state variables. State space: is the hyperspace in which the state of the system takes vales. State Space Model consists of: a) State eqations for a nonlinear system 3-0
21 3- t) k... n... ( n f n t) k... n... ( f t) k... n... ( f = = = & & & b) otpt eqations for a nonlinear system t) k... n... ( m g m y... t) k... n... ( g y t) k... n... ( g y = = = Matri state space model is defined for -state vector [n ] = (t).. (t) (t) (t) n -otpt vector [m ] = (t) y.. (t) y (t) y y(t) m
22 -inpt vector [k ] (t) = k (t) (t).. (t) - right hand side of state eqation f (... n f (... n f ( t) =.. f... k... ( n n k -right hand side of otpt eqation t) t) t) k g ( t) f ( f ( = f... n... n. ( n n k.... k... t) t) t) k the above state and otpt eqations can be written in a compact form as matri eqations & (t) = f ( t) y(t) = g( t) 3-
23 Linearization of Nonlinear Eqations The operation of a system for a short time dration is arond an operating point with small variations abot the eqilibrim point. In this case is sitable to model the operation by approimations of the nonlinear system sing local linear approimations. Sch a linearized LTI model is sed for controllers design. Linearizing state and otpt nonlinear eqations abot a given vale of the state vector gives Linear Time Variant eqations & (t) = A(t) (t) B(t) (t) y(t) = C(t) (t) D(t) (t) A(t) is state matri B(t) is the inpt matri C(t) is the otpt matri D(t) is the direct transmission matri. A Linear Time Invariant (LTI) system is given by & (t) = A(t) B(t) y(t) = C(t) D(t) where A B C and D are matrices of constant vales. Eample: A mass-damper-spring vertical M-B-K system is sbject to a force f(t) K B M Y(t) f(t) The free body diagram is 3-3
24 KY(t) BdY(t)\dt M f(t) Y(t) The eqation of motion is given by Newton second law d Y(t) dy(t) M = B KY(t) f (t) dt dt Inertia force is M d Y(t)/dt and gives the following free body diagram KY(t) M d Y(t)/dt BdY(t)\dt M f(t) Y(t) D Alembert principle gives the force balance eqation d Y(t) dy(t) M B KY(t) f (t) = 0 dt dt Let s define the following set of state variables (t) and (t) (t)=y(t) (t)=dy(t)/dt= d (t)/dt and d (t) /dt= d Y(t)/dt =(/M)[-B dy(t)/dt-ky(t)f(t)] The above second order eqation of motion can be replaced by two first order state eqations sing state variables (t) and (t) and inpt variable (t)=f(t) 3-4
25 d (t)/dt= (t) d (t) /dt= -(K/M) (t) - (B/M) (t) (/M) (t)] and system otpt eqation y(t) = (t) in case that the otpt is the position (or y(t) = (t) if the otpt is the velocity (t)=dy(t)/dt of the mass M) Matri form of state dynamics and otpt eqations is & (t) & (t) y(t) = or = [ 0] 0 - K/M (t) (t) & (t) = A(t) B(t) y(t) = C(t) D(t) where (t) = (t) (t) - B/M (t) (t) 0 / M (t) A = 0 - K/M - B/M 0 B = / M C = D=0 [ 0] 3-5
26 Relationship between Transfer Fnctions and State-Space Eqations System transfer fnction y(s) G (s) = (s) has to be obtained from the matri form of state and otpt eqations by elimination the other variable the vector (s). This reqires Laplace transform of these matri eqations for zero initial conditions s(s) = A(s) B(s) y(s) = C(s) D(s) or (si A) (s) = B(s) y(s) = C(s) D(s) where I is a by identity matri. The elimination is carried ot by solving first eqation for (s) ( s) = (si A) B (s) and replacing (s) in the second eqation y ( s) = C(sI A) B (s) D(s) or y ( s) = [C(sI A) B D] (s) For single inpt y( s) = y(s) and single otpt ( s) = (s) sch that y (s) = [C(sI A) B D](s) This eqation gives the relationship between the transfer fnction and the matrices A B C and D of the state space representation 3-6
27 y(s) G (s) = = [C(sI A) B D] (s) Given that the inverse of the matri si-a is given by the matri of adjoints si-a divided by the determinant si-a (si-a) - = si-a / si-a Assming D=0 G(s) =C si-a B/ si-a Sch that the poles of the transfer fnction G(s) are given by si-a =0 that is actally the eqation that permits the calclation of the vector of eigenvales λ of the sqare matri A λ I-A =0 This indicates that the characteristic eqation of the system is given by (si-a). That means that the eigenvales of A is identical to the poles of G(s) in the Laplace domain. The poles of the transfer fnction G(s) can conseqently be calclated as eigenvales of the matri A. Eample: Consider the M-B-K system analysed before. State space and otpt eqations for the system were given by & (t) (t) 0 = 0 (t) & (t) (t) - K/M - B/M / M (t) y(t) = [ 0] (t) To obtain the transfer fnction from the state space eqations see (See Eample 3-4 from Ogata tetbook) The reslt is G(s)= /( Ms bsk) The poles of G(s) are given by Ms BsK=0 Which is the determinant of the matri A A =0 3-7
28 State Space representation of nth-order systems of Linear Differential Eqations a) The case when the forcing fnction does not involve derivative terms b) When the inpt involves derivative terms in the forcing fnction (See Ch. 3-5 from Ogata tetbook) (See Eample 3-5 in Ogata tetbook) 3-8
Linear System Theory (Fall 2011): Homework 1. Solutions
Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationState Space Models Basic Concepts
Chapter 2 State Space Models Basic Concepts Related reading in Bay: Chapter Section Sbsection 1 (Models of Linear Systems) 1.1 1.1.1 1.1.2 1.1.3 1.1.5 1.2 1.2.1 1.2.2 1.3 In this Chapter we provide some
More informationControl Systems
6.5 Control Systems Last Time: Introdction Motivation Corse Overview Project Math. Descriptions of Systems ~ Review Classification of Systems Linear Systems LTI Systems The notion of state and state variables
More informationClassify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.
Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More information1. State-Space Linear Systems 2. Block Diagrams 3. Exercises
LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space
More informationHybrid modelling and model reduction for control & optimisation
Hybrid modelling and model redction for control & optimisation based on research done by RWTH-Aachen and TU Delft presented by Johan Grievink Models for control and optimiation market and environmental
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More information2E1252 Control Theory and Practice
2E1252 Control Theory and Practice Lectre 11: Actator satration and anti wind-p Learning aims After this lectre, yo shold nderstand how satration can case controller states to wind p know how to modify
More informationLinear and Nonlinear Model Predictive Control of Quadruple Tank Process
Linear and Nonlinear Model Predictive Control of Qadrple Tank Process P.Srinivasarao Research scholar Dr.M.G.R.University Chennai, India P.Sbbaiah, PhD. Prof of Dhanalaxmi college of Engineering Thambaram
More informationDevelopment of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model
Development of Second Order Pls Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Lemma D. Tfa*, M. Ramasamy*, Sachin C. Patwardhan **, M. Shhaimi* *Chemical Engineering Department, Universiti
More informationSareban: Evaluation of Three Common Algorithms for Structure Active Control
Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood
More informationTheoretical and Experimental Implementation of DC Motor Nonlinear Controllers
Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationA fundamental inverse problem in geosciences
A fndamental inverse problem in geosciences Predict the vales of a spatial random field (SRF) sing a set of observed vales of the same and/or other SRFs. y i L i ( ) + v, i,..., n i ( P)? L i () : linear
More informationControl Performance Monitoring of State-Dependent Nonlinear Processes
Control Performance Monitoring of State-Dependent Nonlinear Processes Lis F. Recalde*, Hong Ye Wind Energy and Control Centre, Department of Electronic and Electrical Engineering, University of Strathclyde,
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationUncertainties of measurement
Uncertainties of measrement Laboratory tas A temperatre sensor is connected as a voltage divider according to the schematic diagram on Fig.. The temperatre sensor is a thermistor type B5764K [] with nominal
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationAN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS
AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS Fang-Ming Y, Hng-Yan Chng* and Chen-Ning Hang Department of Electrical Engineering National Central University, Chngli,
More informationAN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS
AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS Fang-Ming Y, Hng-Yan Chng* and Chen-Ning Hang Department of Electrical Engineering National Central University, Chngli,
More informationEstimating models of inverse systems
Estimating models of inverse systems Ylva Jng and Martin Enqvist Linköping University Post Print N.B.: When citing this work, cite the original article. Original Pblication: Ylva Jng and Martin Enqvist,
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationNonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator
Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationDr Ian R. Manchester
Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign
More informationIII. Demonstration of a seismometer response with amplitude and phase responses at:
GG5330, Spring semester 006 Assignment #1, Seismometry and Grond Motions De 30 Janary 006. 1. Calibration Of A Seismometer Using Java: A really nifty se of Java is now available for demonstrating the seismic
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationDesign and Data Acquisition for Thermal Conductivity Matric Suction Sensors
68 TRANSPORTATION RSARCH RCORD 1432 Design and Data Acqisition for Thermal Condctivity Matric Sction Sensors J. K.-M. GAN, D. G. FRDLUND, A. XING, AND W.-X. LI The principles behind sing the thermal condctivity
More informationSystem identification of buildings equipped with closed-loop control devices
System identification of bildings eqipped with closed-loop control devices Akira Mita a, Masako Kamibayashi b a Keio University, 3-14-1 Hiyoshi, Kohok-k, Yokohama 223-8522, Japan b East Japan Railway Company
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationWorkshop on Understanding and Evaluating Radioanalytical Measurement Uncertainty November 2007
1833-3 Workshop on Understanding and Evalating Radioanalytical Measrement Uncertainty 5-16 November 007 Applied Statistics: Basic statistical terms and concepts Sabrina BARBIZZI APAT - Agenzia per la Protezione
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More informationCOMPARISON OF MODEL INTEGRATION APPROACHES FOR FLEXIBLE AIRCRAFT FLIGHT DYNAMICS MODELLING
COMPARISON OF MODEL INTEGRATION APPROACHES FOR FLEXIBLE AIRCRAFT FLIGHT DYNAMICS MODELLING Christian Reschke and Gertjan Looye DLR German Aerospace Center, Institte of Robotics and Mechatronics D-82234
More informationACTUATION AND SIMULATION OF A MINISYSTEM WITH FLEXURE HINGES
The 4th International Conference Comptational Mechanics and Virtal Engineering COMEC 2011 20-22 OCTOBER 2011, Brasov, Romania ACTUATION AND SIMULATION OF A MINISYSTEM WITH FLEXURE HINGES D. NOVEANU 1,
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationGradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How
1 Gradient Projection Anti-windp Scheme on Constrained Planar LTI Systems Jstin Teo and Jonathan P. How Technical Report ACL1 1 Aerospace Controls Laboratory Department of Aeronatics and Astronatics Massachsetts
More informationSimulation Based Analysis of Two Different Control Strategies for PMSM
International Jornal of Engineering Trends and Technology (IJETT) - Volme4Isse4- April 23 Simlation Based Analysis of Two Different Control Strategies for PMSM Lindita Dhamo #, Aida Spahi #2 # Department
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) MAE 5 - inite Element Analysis Several slides from this set are adapted from B.S. Altan, Michigan Technological University EA Procedre for
More information3.4-Miscellaneous Equations
.-Miscellaneos Eqations Factoring Higher Degree Polynomials: Many higher degree polynomials can be solved by factoring. Of particlar vale is the method of factoring by groping, however all types of factoring
More informationINPUT-OUTPUT APPROACH NUMERICAL EXAMPLES
INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES EXERCISE s consider the linear dnamical sstem of order 2 with transfer fnction with Determine the gain 2 (H) of the inpt-otpt operator H associated with this sstem.
More informationThe Oscillatory Stable Regime of Nonlinear Systems, with two time constants
6th WSES International Conference on CIRCUITS SYSTEMS ELECTRONICSCONTROL & SIGNL PROCESSING Cairo Egpt Dec 9-3 7 5 The Oscillator Stable Regime of Nonlinear Sstems with two time constants NUŢU VSILE *
More informationSimulation investigation of the Z-source NPC inverter
octoral school of energy- and geo-technology Janary 5 20, 2007. Kressaare, Estonia Simlation investigation of the Z-sorce NPC inverter Ryszard Strzelecki, Natalia Strzelecka Gdynia Maritime University,
More informationElectron Phase Slip in an Undulator with Dipole Field and BPM Errors
CS-T--14 October 3, Electron Phase Slip in an Undlator with Dipole Field and BPM Errors Pal Emma SAC ABSTRACT A statistical analysis of a corrected electron trajectory throgh a planar ndlator is sed to
More informationState variable feedback
State variable feedbak We have previosly disssed systems desribed by the linear state-spae eqations Ax B y Cx n m with xt () R the internal state, t () R the ontrol inpt, and yt () R the measred otpt.
More informationThe Real Stabilizability Radius of the Multi-Link Inverted Pendulum
Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract
More informationFRÉCHET KERNELS AND THE ADJOINT METHOD
PART II FRÉCHET KERNES AND THE ADJOINT METHOD 1. Setp of the tomographic problem: Why gradients? 2. The adjoint method 3. Practical 4. Special topics (sorce imaging and time reversal) Setp of the tomographic
More informationReflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by
Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More informationFEA Solution Procedure
EA Soltion rocedre (demonstrated with a -D bar element problem) MAE - inite Element Analysis Many slides from this set are originally from B.S. Altan, Michigan Technological U. EA rocedre for Static Analysis.
More informationCOMPARATIVE STUDY OF ROBUST CONTROL TECHNIQUES FOR OTTO-CYCLE MOTOR CONTROL
ACM Symposim Series in Mechatronics - Vol. - pp.76-85 Copyright c 28 by ACM COMPARATIVE STUDY OF ROUST CONTROL TECHNIQUES FOR OTTO-CYCLE MOTOR CONTROL Marcos Salazar Francisco, marcos.salazar@member.isa.org
More informationMove Blocking Strategies in Receding Horizon Control
Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it
More informationIJSER. =η (3) = 1 INTRODUCTION DESCRIPTION OF THE DRIVE
International Jornal of Scientific & Engineering Research, Volme 5, Isse 4, April-014 8 Low Cost Speed Sensor less PWM Inverter Fed Intion Motor Drive C.Saravanan 1, Dr.M.A.Panneerselvam Sr.Assistant Professor
More informationModel Predictive Control Lecture VIa: Impulse Response Models
Moel Preictive Control Lectre VIa: Implse Response Moels Niet S. Kaisare Department of Chemical Engineering Inian Institte of Technolog Maras Ingreients of Moel Preictive Control Dnamic Moel Ftre preictions
More informationPredictive Control- Exercise Session 4 Adaptive Control: Self Tuning Regulators and Model Reference Adaptive Systems
Predictive Control- Exercise Session 4 Adaptive Control: Self Tning Reglators and Model Reference Adaptive Systems 1. Indirect Self Tning Reglator: Consider the system where G(s)=G 1 (s)g 2 (s) G 1 (s)=
More informationReducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation
JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationDynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining
Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining P. Hgenin*, B. Srinivasan*, F. Altpeter**, R. Longchamp* * Laboratoire d Atomatiqe,
More informationMATHEMATICAL MODELING OF CONTROL SYSTEMS
1 MATHEMATICAL MODELING OF CONTROL SYSTEMS Sep-14 Dr. Mohammed Morsy Outline Introduction Transfer function and impulse response function Laplace Transform Review Automatic control systems Signal Flow
More informationConditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane
Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching
More informationControl of a Power Assisted Lifting Device
Proceedings of the RAAD 212 21th International Workshop on Robotics in Alpe-Adria-Danbe Region eptember 1-13, 212, Napoli, Italy Control of a Power Assisted Lifting Device Dimeas otios a, Kostompardis
More informationA Model-Free Adaptive Control of Pulsed GTAW
A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationEXPT. 5 DETERMINATION OF pk a OF AN INDICATOR USING SPECTROPHOTOMETRY
EXPT. 5 DETERMITIO OF pk a OF IDICTOR USIG SPECTROPHOTOMETRY Strctre 5.1 Introdction Objectives 5.2 Principle 5.3 Spectrophotometric Determination of pka Vale of Indicator 5.4 Reqirements 5.5 Soltions
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationLeast-squares collocation with covariance-matching constraints
J Geodesy DOI 1.17/s19-7-133-5 ORIGINAL ARTICLE Least-sqares collocation with covariance-matching constraints Christopher Kotsakis Received: 18 Agst 26 / Accepted: 28 December 26 Springer-Verlag 27 Abstract
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationEXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.
.4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationAn extremum seeking approach to parameterised loop-shaping control design
Preprints of the 19th World Congress The International Federation of Atomatic Control An extremm seeking approach to parameterised loop-shaping control design Chih Feng Lee Sei Zhen Khong Erik Frisk Mattias
More informationAdjoint-Based Sensitivity Analysis for Computational Fluid Dynamics
Adjoint-Based Sensitivity Analysis for Comptational Flid Dynamics Dimitri J. Mavriplis Department of Mecanical Engineering niversity of Wyoming Laramie, WY Motivation Comptational flid dynamics analysis
More informationCreating a Sliding Mode in a Motion Control System by Adopting a Dynamic Defuzzification Strategy in an Adaptive Neuro Fuzzy Inference System
Creating a Sliding Mode in a Motion Control System by Adopting a Dynamic Defzzification Strategy in an Adaptive Nero Fzzy Inference System M. Onder Efe Bogazici University, Electrical and Electronic Engineering
More information3 2D Elastostatic Problems in Cartesian Coordinates
D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates
More informationMathematical Models of Physiological Responses to Exercise
Mathematical Models of Physiological Responses to Exercise Somayeh Sojodi, Benjamin Recht, and John C. Doyle Abstract This paper is concerned with the identification of mathematical models for physiological
More informationMethods of Design-Oriented Analysis The GFT: A Final Solution for Feedback Systems
http://www.ardem.com/free_downloads.asp v.1, 5/11/05 Methods of Design-Oriented Analysis The GFT: A Final Soltion for Feedback Systems R. David Middlebrook Are yo an analog or mixed signal design engineer
More informationSafe Manual Control of the Furuta Pendulum
Safe Manal Control of the Frta Pendlm Johan Åkesson, Karl Johan Åström Department of Atomatic Control, Lnd Institte of Technology (LTH) Box 8, Lnd, Sweden PSfrag {jakesson,kja}@control.lth.se replacements
More informationOn the circuit complexity of the standard and the Karatsuba methods of multiplying integers
On the circit complexity of the standard and the Karatsba methods of mltiplying integers arxiv:1602.02362v1 [cs.ds] 7 Feb 2016 Igor S. Sergeev The goal of the present paper is to obtain accrate estimates
More informationMULTIOBJECTIVE OPTIMAL STRUCTURAL CONTROL OF THE NOTRE DAME BUILDING MODEL BENCHMARK *
Thi d d i h F M k 4 0 4 Earthqake Engineering and Strctral Dynamics, in press. MULTIOBJECTIVE OPTIMAL STRUCTURAL CONTROL OF THE NOTRE DAME BUILDING MODEL BENCHMARK * E. A. JOHNSON, P. G. VOULGARIS, AND
More informationMaterial Transport with Air Jet
Material Transport with Air Jet Dr. István Patkó Bdapest Tech Doberdó út 6, H-1034 Bdapest, Hngary patko@bmf.h Abstract: In the field of indstry, there are only a very few examples of material transport
More informationDEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1. C. Silvestre A. Pascoal
DEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1 C. Silvestre A. Pascoal Institto Sperior Técnico Institte for Systems and Robotics Torre Norte - Piso 8, Av. Rovisco
More informationPhysicsAndMathsTutor.com
. Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the
More informationAlireza Mousavi Brunel University
Alireza Mousavi Brunel University 1 » Control Process» Control Systems Design & Analysis 2 Open-Loop Control: Is normally a simple switch on and switch off process, for example a light in a room is switched
More informationOnline Solution of State Dependent Riccati Equation for Nonlinear System Stabilization
> REPLACE American HIS Control LINE Conference WIH YOUR PAPER IDENIFICAION NUMBER (DOUBLE-CLICK HERE O EDI) FrC3. Marriott Waterfront, Baltimore, MD, USA Jne 3-Jly, Online Soltion of State Dependent Riccati
More informationm = Average Rate of Change (Secant Slope) Example:
Average Rate o Change Secant Slope Deinition: The average change secant slope o a nction over a particlar interval [a, b] or [a, ]. Eample: What is the average rate o change o the nction over the interval
More informationBumpless transfer for discrete-time switched systems
29 American Control Conference Hyatt Regency Riverfront, St. Lois, O, USA Jne 1-12, 29 WeB12.3 Bmpless transfer for discrete-time switched systems I. ALLOCI, L. HETEL, J. DAAFOUZ, ember, IEEE, C. IUNG,
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationSTABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION
8 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION Čeliovsý S.*, Hospodář P.** *CTU Prage, Faclty
More informationA Macroscopic Traffic Data Assimilation Framework Based on Fourier-Galerkin Method and Minimax Estimation
A Macroscopic Traffic Data Assimilation Framework Based on Forier-Galerkin Method and Minima Estimation Tigran T. Tchrakian and Sergiy Zhk Abstract In this paper, we propose a new framework for macroscopic
More informationReduction of over-determined systems of differential equations
Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical
More informationAdaptive Fault-tolerant Control with Control Allocation for Flight Systems with Severe Actuator Failures and Input Saturation
213 American Control Conference (ACC) Washington, DC, USA, Jne 17-19, 213 Adaptive Falt-tolerant Control with Control Allocation for Flight Systems with Severe Actator Failres and Inpt Satration Man Wang,
More informationMEG 741 Energy and Variational Methods in Mechanics I
MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)
More informationComputational Geosciences 2 (1998) 1, 23-36
A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), 23-36 Abstract.
More informationMULTIPLE REGRESSION WITH TWO EXPLANATORY VARIABLES: EXAMPLE. ASVABC + u
MULTIPLE REGRESSION MULTIPLE REGRESSION WITH TWO EXPLANATORY VARIABLES: EXAMPLE EARNINGS α + HGC + ASVABC + α EARNINGS ASVABC HGC This seqence provides a geometrical interpretation of a mltiple regression
More information3.3 Operations With Vectors, Linear Combinations
Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram
More informationDirect Modeling and Robust Control of a Servo-pneumatic System. Ph.D. Thesis
Bdapest University of Technology and Economics Department of Mechatronics, Optics and Engineering Informatics Direct Modeling and Robst Control of a Servo-pnematic System Ph.D. Thesis Károly Széll Spervisor:
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
More informationData-Efficient Control Policy Search using Residual Dynamics Learning
Data-Efficient Control Policy Search sing Residal Dynamics Learning Matteo Saveriano 1, Ychao Yin 1, Pietro Falco 1 and Donghei Lee 1,2 Abstract In this work, we propose a model-based and data efficient
More information