Lecture 2 – Linear Systems

Total Page:16

File Type:pdf, Size:1020Kb

Lecture 2 – Linear Systems Lecture 2 – Linear Systems This lecture: EE263 material recap + some controls motivation • Continuous time (physics) • Linear state space model • Transfer functions • Black-box models; frequency domain analysis • Linearization EE392m - Spring 2005 Control Engineering 2-1 Gorinevsky Modeling and Analysis Input signal System Model Output signal (control) [internal states] (observation) state evolution This lecture considers • Linear models. More detail on modeling in Lecture 7 • Simulation: computing state evolution and output signal • Stability: does the solution diverge after some time? • Approximate linear models EE392m - Spring 2005 Control Engineering 2-2 Gorinevsky Linear Models • Model is a mathematical representations of a system – Models allow simulating the system – Models can be used for conceptual analysis – Models are never exact • Linear models – Have simple structure – Can be analyzed using powerful mathematical tools – Can be matched against real data using known procedures – Many simple physics models are linear – They are just models, not the real systems EE392m - Spring 2005 Control Engineering 2-3 Gorinevsky State space model • Generic state space model Example: F16 Longitudinal Model – is described by ODEs from Aircraft Control and Simulation – e.g., physics-based system by Stevens and Lewis model θ α u = δ – state vector x e V – observation vector y x – control vector u δ dx y = 57.3θ z e = f (x,u,t) state evolution dt = = ⋅ -2 + α − θ − + y g(x,t) observation V& -1.93 10 V 8.82 32.2 0.58q 0.17δe α = ⋅ -4 − α + − ⋅ -3δ x1 - velocity V [ft/sec] & -2.54 10 V 1.02 0.91q 2.15 10 e x2 - angle of attack α [rad] θ& = q x3 - pitch angle θ [rad] q = 2.95 ⋅10-12V + 0.82α −1.08q − 0.18δ x4 - pitch rate q [rad/sec] & e δ e - elevator deflection [deg] EE392m - Spring 2005 Control Engineering 2-4 Gorinevsky Linear state space model • Linear Time Invariant (LTI) state space model: dx = Ax + Bu state evolution dt y = Cx observations • Can be integrated analytically or numerically (simulation) • Can be well analyzed: stability, response Example: F16 Longitudinal Model ⎡−1.93⋅10−2 8.82 − 32.2 − 0.48⎤ ⎡ 0.17 ⎤ ⎢ −4 ⎥ ⎢ −3 ⎥ dx − 2.54 ⋅10 −1.02 0 0.91 − 2.15⋅10 = ⎢ ⎥ ⋅ x + ⎢ ⎥ ⋅ u dt ⎢ 0 0 0 1 ⎥ ⎢ 0 ⎥ ⎢ −12 ⎥ ⎢ − ⎥ ⎣ 2.95⋅10 0.82 0 −1.08⎦ ⎣ 0.18 ⎦ y = []0 0 57.3 0 ⋅ x EE392m - Spring 2005 Control Engineering 2-5 Gorinevsky Integrating linear autonomous system dx Example: = Ax >>A = [-1.93e-2 8.82 -32.2 -0.58; dt -2.54e-4 -1.02 0 0.91; 0 0 0 1; • Matrix exponential 2.95e-12 0.82 0 -1.08]; – Can be computed in Matlab >> expm(A) as expm(A) ans = 0.9802 3.1208 -31.8830 -10.2849 – The definition corresponds -0.0002 0.5004 0.0031 0.3604 -0.0000 0.2216 1.0002 0.6709 to integrating the ODE by -0.0001 0.3239 0.0007 0.4773 Euler method x(T + ∆t) ≈ x(T ) + ∆tAx(t) = ()I + ∆tA x(T ) x(T + n∆t) = ()I + ∆tA n x(T ) → exp( An∆t)x(T ) ∆ exp( At) = lim ()I + ∆tA t / t ∆t→0 EE392m - Spring 2005 Control Engineering 2-6 Gorinevsky Integrating linear autonomous system 3 dx = Ax, x(0) = x 2 dt 0 1 0 5 10 15 20 • Initial condition response 0.02 = 0 x(t) exp( At)x0 -0.02 -0.04 0 5 10 15 20 0.02 Example: 0 %Take A, B, C from the F16 example -0.02 -0.04 >> k = 0.02 -0.06 >> G = A + k*B*C; -0.08 >> x0 = [0.1700; 0 5 10 15 20 -0.0022; 0.05 0; 0 -0.1800]; -0.05 >> -0.1 for j = 1:length(t); -0.15 x(:,j)=expm(G*t(j))*x0; end; 0 5 10 15 20 EE392m - Spring 2005 Control EngineeringTIME 2-7 Gorinevsky Eigenvalues and Stability Example: % take A from % the F16 example • Consider eigenvalues of the matrix A >> eig(A) ans = λ = λ − = -1.9125 { j } eig( A); det( I A) 0 -0.1519 + 0.1143i -0.1519 - 0.1143i 0.0970 • Suppose A has all different and nonzero eigenvalues, then = ⋅ λ ⋅ −1 A V diag{ j } V λ − exp( At) = V ⋅ diag{e jt } ⋅V 1 λ < • The system solution is exponentially stable if Re j 0 • If A has eigenvalues with multiplicity more than 1, things are a bit more complicated: Jordan blocks, polynomials in t λ < • Still the condition of exponentially stability is Re j 0 EE392m - Spring 2005 Control Engineering 2-8 Gorinevsky Input-output models • Black-box models – describe system P as an operator u P y input signal output signal x internal state (hidden) • Historically (50 years ago) – Black-box models EE – State-space models ME, AA EE392m - Spring 2005 Control Engineering 2-9 Gorinevsky Linear System (input-output) • Linearity ⋅ ⎯⎯→P ⋅ ⋅ ⎯⎯→P ⋅ u1( ) y1( ) u2 ( ) y2 ( ) ⋅ + ⋅ ⎯⎯→P ⋅ + ⋅ au1( ) bu2 ( ) ay1( ) by2 ( ) • Linear Time-Invariant systems - LTI u(⋅ − T ) ⎯⎯→P y(⋅ − T ) u y ⎯⎯→P t t EE392m - Spring 2005 Control Engineering 2-10 Gorinevsky Convolution representation • Convolution integral t signal y(t) = ∫ h(t −τ )u(τ )dτ processing y = h * u −∞ notation • Impulse response u(t) = δ (t) ⇒ y(t) = h(t) • Step response: u = 1 for t > 0 t t d g(t) = ∫ h(t −τ )dτ =∫ h(τ )dτ h(t) = g(t) 0 0 dt EE392m - Spring 2005 Control Engineering 2-11 Gorinevsky Impulse Response for State Space Model dx Example: = + Ax Bu state evolution >> A = [-0.3130 56.7000 0; dt -0.0139 -0.4260 0; y = Cx observation 0 56.7000 0]; >> B = [0.232; 0.0203; 0]; >> C = [0, 0, 1]; • Impulse response for the state x 1 x(∆t) ≈ ∆tAx(0) + B∆tu 0.5 { 0 1 -0.5 x(∆t) ≈ B∆tu = B 0 10 20 0.04 x(t) = exp( At)B 0.02 0 -0.02 0 10 20 • System impulse response 1 h(t) = C exp( At)B 0.5 0 0 10 20 TIME EE392m - Spring 2005 Control Engineering 2-12 Gorinevsky Formal transfer function • Rational transfer function = IIR (Infinite Impulse Response) model – Broad class of input-output linear models d m y d m−1 y d n u d n−1u a + a + + a + y = b + b + + b + u 1 dt m 2 dt m−1 K m 1 1 dt n 2 dt n−1 K n 1 d • Differentiation operator → s dt • Formal transfer function – rational function of s N (s) N (s) = a s m + + a s + a y = H (s) ⋅ u = ⋅ u 1 K m m+1 D(s) = n + + + D(s) b1s K bn s bn+1 • For a causal system m ≤ n EE392m - Spring 2005 Control Engineering 2-13 Gorinevsky Poles, Impulse Response N (s) N (s) = 0 - zeros y = ⋅ u D(s) D(s) = 0 - poles = − M 1 ⋅ ⋅ − M K • Expand D(s) (s p1 ) K (s pK ) N (s) N (s) Then y = 1 ⋅ u + + K ⋅ u − M 1 K − M K (s p1 ) (s p K ) • Quasi-polynomial impulse response – see a textbook − − M 1 1 p1t M K 1 pK t h(t) = (c t + + c + )e + + (c t + + c + )e 1,1 K 1,M 1 1 K K ,1 K K ,M K 1 2 d y = 1 ⋅ •Example: = u Transfer function: y 2 u dt 2 s = + Impulse response: h(t) v0t x0 EE392m - Spring 2005 Control Engineering 2-14 Gorinevsky Transfer Function for State Space • Formal transfer function for a state space model − = + y = ()sI − A 1 B ⋅ u d → sx Ax Bu s − dt y = Cx H (s) = C()sI − A 1 B • Characteristic polynomial N (s) = det()sI − A = 0 Poles det()sI − A = 0 eigenvalues • Poles are the same as eigenvalues of the state-space matrix A λ • For stability we need Re pk = Re k < 0 EE392m - Spring 2005 Control Engineering 2-15 Gorinevsky Laplace transform ∞ • Laplace integral transform: x(t) → xˆ(s) = ∫ x(t)e st dt 0 • Laplace transform of the convolution integral yields ∞ Imag s yˆ(s) = H (s)uˆ(s) H (s) = ∫ h(t)est dt 0 • Transfer function: H (s) Re s – function of complex variable s – analytical in a right half-plane Re s ≥ a – for a stable system a ≤ 0 N(s) dx –for an IIR model H(s) = → sxˆ(s) D(s) dt EE392m - Spring 2005 Control Engineering 2-16 Gorinevsky Frequency decomposition • Sinusoids are eigenfunctions of an LTI system: y = H (s)u LTI Plant iωt iωt ω e H (i )e ω d ω ω s()ei t → ei t = iωei t dt • Frequency domain analysis ω ω i kt → = ω i kt ∑uk e y ∑uk H(i k )e EE392m - Spring 2005 Control Engineering 2-17 Gorinevsky Frequency domain description • Frequency domain analysis 1 − ω 1 − ω u = u~(ω)e i t dω ⇒ y = H(iω)u~(ω)e i t dω 2π ∫ 2π ∫ 142 43 ~y (ω) ω Packet ei t Packet eiωt u of u~(ω) H(iω) of ~y(ω) y sinusoids sinusoids • Fourier transform – numerical analysis • Laplace transform – complex analysis u(t) = 0, for t < 0 ∞ ∞ u~(ω) = u(t)eiωt dω = ⎛ u(t)est dω⎞ = uˆ(iω) −∞ ⎜ −∞ ⎟ ∫ ⎝∫ ⎠s=iω EE392m - Spring 2005 Control Engineering 2-18 Gorinevsky Continuous systems in frequency domain ∞ ~x(ω) = ∫ x(t)eiωt dt • Fourier transform −∞ [−∞,∞] → [−∞,∞] ∞ 1 − ω x(t) = ~x(ω)e i t dω • Inverse Fourier transform π ∫ 2 −∞ t y(t) = ∫ h(t −τ )u(τ )dt • I/O impulse response model −∞ ∞ H (s) = ∫ h(t)est dt • Transfer function 0 ~ ω = ω ~ ω y( ) H (i )u( ) • System frequency response EE392m - Spring 2005 Control Engineering 2-19 Gorinevsky Frequency domain description M (ω) = H(iω) • Bode plots: ϕ(ω) = arg H(iω) u = eiωt Bode Diagram 15 iωt = ω y H(i )e 10 • Example: 5 1 0 Magnitude (dB) H(s) = -5 s − 0.7 0 -45 •|H| is often measured -90 in dB -135 – [dB] = 20 log M (deg) Phase 10 -180 -2 -1 0 10 10 10 Frequency (rad/sec) EE392m - Spring 2005 Control Engineering 2-20 Gorinevsky Model Approximation • Model structure – physics, computational • Determine parameters from data • Step/impulse responses are close the input/output models are close • Example – fit step response • Linearization of nonlinear model EE392m - Spring 2005 Control Engineering 2-21 Gorinevsky Black-box model from data • Linear black-box model can be determined from the data, e.g., step response data, or frequency response STEP RESPONSE 1 0.8 0.6 0.4 HEAT FLUX 0.2 0 0 20 40 60 80 100 TIME • Example problem: fit an IIR model of a given order • This is called model identification • Considered in more detail in Lecture 8 EE392m - Spring 2005 Control Engineering 2-22 Gorinevsky Linear PDE models • Include functions of spatial variables – electromagnetic fields – mass and heat transfer – fluid dynamics – structural deformations • Example: sideways heat equation x ∂ ∂2 T = T k 2 ∂t ∂x Tinside=u Toutside=0 T (0) = u; T (1)
Recommended publications
  • 2.161 Signal Processing: Continuous and Discrete Fall 2008
    MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts Institute of Technology Department of Mechanical Engineering 2.161 Signal Processing - Continuous and Discrete Fall Term 2008 1 Lecture 1 Reading: • Class handout: The Dirac Delta and Unit-Step Functions 1 Introduction to Signal Processing In this class we will primarily deal with processing time-based functions, but the methods will also be applicable to spatial functions, for example image processing. We will deal with (a) Signal processing of continuous waveforms f(t), using continuous LTI systems (filters). a LTI dy nam ical system input ou t pu t f(t) Continuous Signal y(t) Processor and (b) Discrete-time (digital) signal processing of data sequences {fn} that might be samples of real continuous experimental data, such as recorded through an analog-digital converter (ADC), or implicitly discrete in nature. a LTI dis crete algorithm inp u t se q u e n c e ou t pu t seq u e n c e f Di screte Si gnal y { n } { n } Pr ocessor Some typical applications that we look at will include (a) Data analysis, for example estimation of spectral characteristics, delay estimation in echolocation systems, extraction of signal statistics. (b) Signal enhancement. Suppose a waveform has been contaminated by additive “noise”, for example 60Hz interference from the ac supply in the laboratory. 1copyright c D.Rowell 2008 1–1 inp u t ou t p u t ft + ( ) Fi lte r y( t) + n( t) ad d it ive no is e The task is to design a filter that will minimize the effect of the interference while not destroying information from the experiment.
    [Show full text]
  • Lecture 2 LSI Systems and Convolution in 1D
    Lecture 2 LSI systems and convolution in 1D 2.1 Learning Objectives Understand the concept of a linear system. • Understand the concept of a shift-invariant system. • Recognize that systems that are both linear and shift-invariant can be described • simply as a convolution with a certain “impulse response” function. 2.2 Linear Systems Figure 2.1: Example system H, which takes in an input signal f(x)andproducesan output g(x). AsystemH takes an input signal x and produces an output signal y,whichwecan express as H : f(x) g(x). This very general definition encompasses a vast array of di↵erent possible systems.! A subset of systems of particular relevance to signal processing are linear systems, where if f (x) g (x)andf (x) g (x), then: 1 ! 1 2 ! 2 a f + a f a g + a g 1 · 1 2 · 2 ! 1 · 1 2 · 2 for any input signals f1 and f2 and scalars a1 and a2. In other words, for a linear system, if we feed it a linear combination of inputs, we get the corresponding linear combination of outputs. Question: Why do we care about linear systems? 13 14 LECTURE 2. LSI SYSTEMS AND CONVOLUTION IN 1D Question: Can you think of examples of linear systems? Question: Can you think of examples of non-linear systems? 2.3 Shift-Invariant Systems Another subset of systems we care about are shift-invariant systems, where if f1 g1 and f (x)=f (x x )(ie:f is a shifted version of f ), then: ! 2 1 − 0 2 1 f (x) g (x)=g (x x ) 2 ! 2 1 − 0 for any input signals f1 and shift x0.
    [Show full text]
  • Modelling and Control Methods with Applications to Mechanical Waves
    Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 1174 Modelling and Control Methods with Applications to Mechanical Waves HANS NORLANDER ACTA UNIVERSITATIS UPSALIENSIS ISSN 1651-6214 ISBN 978-91-554-9023-2 UPPSALA urn:nbn:se:uu:diva-229793 2014 Dissertation presented at Uppsala University to be publicly examined in room 2446, ITC, Lägerhyddsvägen 2, Uppsala, Friday, 17 October 2014 at 10:15 for the degree of Doctor of Philosophy. The examination will be conducted in English. Faculty examiner: Prof. Jonas Sjöberg (Chalmers). Abstract Norlander, H. 2014. Modelling and Control Methods with Applications to Mechanical Waves. Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 1174. 72 pp. Uppsala: Acta Universitatis Upsaliensis. ISBN 978-91-554-9023-2. Models, modelling and control design play important parts in automatic control. The contributions in this thesis concern topics in all three of these concepts. The poles are of fundamental importance when analyzing the behaviour of a system, and pole placement is an intuitive and natural approach for control design. A novel parameterization for state feedback gains for pole placement in the linear multiple input case is presented and analyzed. It is shown that when the open and closed loop poles are disjunct, every state feedback gain can be parameterized. Other properties are also investigated. Hammerstein models have a static non-linearity on the input. A method for exact compensation of such non-linearities, combined with introduction of integral action, is presented. Instead of inversion of the non-linearity the method utilizes differentiation, which in many cases is simpler.
    [Show full text]
  • Design Analysis Method for Multidisciplinary Complex Product Using Sysml
    MATEC Web of Conferences 139, 00014 (2017) DOI: 10.1051/matecconf/201713900014 ICMITE 2017 Design Analysis Method for Multidisciplinary Complex Product using SysML Jihong Liu1,*, Shude Wang1, and Chao Fu1 1 School of Mechanical Engineering and Automation, Beihang University, 100191 Beijing, China Abstract. In the design of multidisciplinary complex products, model-based systems engineering methods are widely used. However, the methodologies only contain only modeling order and simple analysis steps, and lack integrated design analysis methods supporting the whole process. In order to solve the problem, a conceptual design analysis method with integrating modern design methods has been proposed. First, based on the requirement analysis of the quantization matrix, the user’s needs are quantitatively evaluated and translated to system requirements. Then, by the function decomposition of the function knowledge base, the total function is semi-automatically decomposed into the predefined atomic function. The function is matched into predefined structure through the behaviour layer using function-structure mapping based on the interface matching. Finally based on design structure matrix (DSM), the structure reorganization is completed. The process of analysis is implemented with SysML, and illustrated through an aircraft air conditioning system for the system validation. 1 Introduction decomposition and function modeling, function structure mapping and evaluation, and structure reorganization. In the process of complex product design, system The whole process of the analysis method is aimed at engineering is a kind of development methods which providing designers with an effective analysis of ideas covers a wide range of applications across from system and theoretical support to reduce the bad design and requirements analysis, function decomposition to improve design efficiency.
    [Show full text]
  • Introduction to Aircraft Aeroelasticity and Loads
    JWBK209-FM-I JWBK209-Wright November 14, 2007 2:58 Char Count= 0 Introduction to Aircraft Aeroelasticity and Loads Jan R. Wright University of Manchester and J2W Consulting Ltd, UK Jonathan E. Cooper University of Liverpool, UK iii JWBK209-FM-I JWBK209-Wright November 14, 2007 2:58 Char Count= 0 Introduction to Aircraft Aeroelasticity and Loads i JWBK209-FM-I JWBK209-Wright November 14, 2007 2:58 Char Count= 0 ii JWBK209-FM-I JWBK209-Wright November 14, 2007 2:58 Char Count= 0 Introduction to Aircraft Aeroelasticity and Loads Jan R. Wright University of Manchester and J2W Consulting Ltd, UK Jonathan E. Cooper University of Liverpool, UK iii JWBK209-FM-I JWBK209-Wright November 14, 2007 2:58 Char Count= 0 Copyright C 2007 John Wiley & Sons Ltd, The Atrium, Southern Gate, Chichester, West Sussex PO19 8SQ, England Telephone (+44) 1243 779777 Email (for orders and customer service enquiries): [email protected] Visit our Home Page on www.wileyeurope.com or www.wiley.com All Rights Reserved. No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form or by any means, electronic, mechanical, photocopying, recording, scanning or otherwise, except under the terms of the Copyright, Designs and Patents Act 1988 or under the terms of a licence issued by the Copyright Licensing Agency Ltd, 90 Tottenham Court Road, London W1T 4LP, UK, without the permission in writing of the Publisher. Requests to the Publisher should be addressed to the Permissions Department, John Wiley & Sons Ltd, The Atrium, Southern Gate, Chichester, West Sussex PO19 8SQ, England, or emailed to [email protected], or faxed to (+44) 1243 770620.
    [Show full text]
  • Linear System Theory
    Chapter 2 Linear System Theory In this course, we will be dealing primarily with linear systems, a special class of sys- tems for which a great deal is known. During the first half of the twentieth century, linear systems were analyzed using frequency domain (e.g., Laplace and z-transform) based approaches in an effort to deal with issues such as noise and bandwidth issues in communication systems. While they provided a great deal of intuition and were suf- ficient to establish some fundamental results, frequency domain approaches presented various drawbacks when control scientists started studying more complicated systems (containing multiple inputs and outputs, nonlinearities, noise, and so forth). Starting in the 1950’s (around the time of the space race), control engineers and scientists started turning to state-space models of control systems in order to address some of these issues. These time-domain approaches are able to effectively represent concepts such as the in- ternal state of the system, and also present a method to introduce optimality conditions into the controller design procedure. We will be using the state-space (or “modern”) approach to control almost exclusively in this course, and the purpose of this chapter is to review some of the essential concepts in this area. 2.1 Discrete-Time Signals Given a field F,asignal is a mapping from a set of numbers to F; in other words, signals are simply functions of the set of numbers that they operate on. More specifically: • A discrete-time signal f is a mapping from the set of integers Z to F, and is denoted by f[k]fork ∈ Z.Eachinstantk is also known as a time-step.
    [Show full text]
  • Digital Image Processing Lectures 3 & 4
    2-D Signals and Systems 2-D Fourier Transform Digital Image Processing Lectures 3 & 4 M.R. Azimi, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2017 M.R. Azimi Digital Image Processing 2-D Signals and Systems 2-D Fourier Transform Definitions and Extensions: 2-D Signals: A continuous image is represented by a function of two variables e.g. x(u; v) where (u; v) are called spatial coordinates and x is the intensity. A sampled image is represented by x(m; n). If pixel intensity is also quantized (digital images) then each pixel is represented by B bits (typically B = 8 bits/pixel). 2-D Delta Functions: They are separable 2-D functions i.e. 1 (u; v) = (0; 0) Dirac: δ(u; v) = δ(u) δ(v) = 0 Otherwise 1 (m; n) = (0; 0) Kronecker: δ(m; n) = δ(m) δ(n) = 0 Otherwise M.R. Azimi Digital Image Processing 2-D Signals and Systems 2-D Fourier Transform Properties: For 2-D Dirac delta: 1 1- R R x(u0; v0)δ(u − u0; v − v0)du0dv0 = x(u; v) −∞ 2- R R δ(u; v)du dv = 1 8 > 0 − For 2-D Kronecker delta: 1 1 1- x(m; n) = P P x(m0; n0)δ(m − m0; n − n0) m0=−∞ n0=−∞ 1 2- P P δ(m; n) = 1 m;n=−∞ Periodic Signals: Consider an image x(m; n) which satisfies x(m; n + N) = x(m; n) x(m + M; n) = x(m; n) This signal is said to be doubly periodic with horizontal and vertical periods M and N, respectively.
    [Show full text]
  • CHAPTER TWO - Static Aeroelasticity – Unswept Wing Structural Loads and Performance 21 2.1 Background
    Static aeroelasticity – structural loads and performance CHAPTER TWO - Static Aeroelasticity – Unswept wing structural loads and performance 21 2.1 Background ........................................................................................................................... 21 2.1.2 Scope and purpose ....................................................................................................................... 21 2.1.2 The structures enterprise and its relation to aeroelasticity ............................................................ 22 2.1.3 The evolution of aircraft wing structures-form follows function ................................................ 24 2.2 Analytical modeling............................................................................................................... 30 2.2.1 The typical section, the flying door and Rayleigh-Ritz idealizations ................................................ 31 2.2.2 – Functional diagrams and operators – modeling the aeroelastic feedback process ....................... 33 2.3 Matrix structural analysis – stiffness matrices and strain energy .......................................... 34 2.4 An example - Construction of a structural stiffness matrix – the shear center concept ........ 38 2.5 Subsonic aerodynamics - fundamentals ................................................................................ 40 2.5.1 Reference points – the center of pressure..................................................................................... 44 2.5.2 A different
    [Show full text]
  • Understanding the Rational Function Model: Methods and Applications
    UNDERSTANDING THE RATIONAL FUNCTION MODEL: METHODS AND APPLICATIONS Yong Hu, Vincent Tao, Arie Croitoru GeoICT Lab, York University, 4700 Keele Street, Toronto M3J 1P3 - {yhu, tao, ariec}@yorku.ca KEY WORDS: Photogrammetry, Remote Sensing, Sensor Model, High-resolution, Satellite Imagery ABSTRACT: The physical and generalized sensor models are two widely used imaging geometry models in the photogrammetry and remote sensing. Utilizing the rational function model (RFM) to replace physical sensor models in photogrammetric mapping is becoming a standard way for economical and fast mapping from high-resolution images. The RFM is accepted for imagery exploitation since high accuracies have been achieved in all stages of the photogrammetric process just as performed by rigorous sensor models. Thus it is likely to become a passkey in complex sensor modeling. Nowadays, commercial off-the-shelf (COTS) digital photogrammetric workstations have incorporated the RFM and related techniques. Following the increasing number of RFM related publications in recent years, this paper reviews the methods and key applications reported mainly over the past five years, and summarizes the essential progresses and address the future research directions in this field. These methods include the RFM solution, the terrain- independent and terrain-dependent computational scenarios, the direct and indirect RFM refinement methods, the photogrammetric exploitation techniques, and photogrammetric interoperability for cross sensor/platform imagery integration. Finally, several open questions regarding some aspects worth of further study are addressed. 1. INTRODUCTION commercial companies, such as Space Imaging (the first high- resolution satellite imagery vendor), to adopt the RFM scheme A sensor model describes the geometric relationship between in order to deliver the imaging geometry model has also the object space and the image space, or vice visa.
    [Show full text]
  • Reliability Block Diagram (Rbd)
    RELIABILITY BLOCK DIAGRAM (RBD) 1. System performance 2. RBD definition 3. Structure function Master ISMP Master ISMP - Castanier 4. Reliability function of a non repairable complex system 5. Reliability metrics for a repairable complex system 59 System performance Problem: Individual definition of the failure modes and components What about the global performance of the « system »? Master ISMP Master ISMP - Castanier 60 Reliability Block Diagram definition Definition A Reliability Block Diagram is a success-oriented graph which illustrates from a logical perspective how the different functional blocks ensure the global mission/function of the system. The structure of the reliability block diagram is mathematically described through structure functions which allow to assess some reliability measures of the (complex) Master ISMP Master ISMP - Castanier system. Function model EThe valve can be closed S and stop the flow SDV1 61 Reliability Block Diagram definition Decomposition of a complex system in functional blocks Component 1 8 9 4 4 2 7 8 10 Master ISMP Master ISMP - Castanier 5 6 3 9 10 Functional block 1 62 Reliability Block Diagram definition The classical architectures « Serie » structure A system which is operating if and only of all of its components are operating is called a serie structure/system. « Parallel » structure A system which is operating if at least one of its components is Master ISMP Master ISMP - Castanier operating is called a parallel structure/system. « oo » structure A system which is operating if at least of its component are operating is called a over structure/system 63 Reliability Block Diagram definition Example Let consider 2 independent safety-related valves V1 and V2 which are physically installed in series.
    [Show full text]
  • A Model-Driven Data-Analysis Architecture Enabling Reuse and Insight in Open Data
    A model-driven data-analysis architecture enabling reuse and insight in open data Robin Hoogervorst July 2018 Master's Thesis Master of Computer Science Specialization Software Technology University of Twente Faculty of Electrical Engineering, Mathematics and Computer Science Supervising committee: dr. Luis Ferreira Pires dr.ir. Maurice van Keulen prof.dr.ir. Arend Rensink Abstract The last years have shown an increase in publicly available data, named open data. Organisations can use open data to enhance data analysis, but tradi- tional data solutions are not suitable for data sources not controlled by the organisation. Hence, each external source needs a specific solution to solve accessing its data, interpreting it and provide possibilities verification. Lack of proper standards and tooling prohibits generalization of these solutions. Structuring metadata allows structure and semantics of these datasets to be described. When this structure is properly designed, these metadata can be used to specify queries in an abstract manner, and translated these to dataset its storage platform. This work uses Model-Driven Engineering to design a metamodel able to represent the structure different open data sets as metadata. In addition, a function metamodel is designed and used to define operations in terms of these metadata. Transformations are defined using these functions to generate executable code, able to execute the required data operations. Other transformations apply the same operations to the metadata model, allowing parallel transformation of metadata and data, keeping them synchronized. The definition of these metamodels, as well as their transformations are used to develop a prototype application framework able to load external datasets and apply operations to the data and metadata simultaneously.
    [Show full text]
  • Integration Definition for Function Modeling (IDEF0)
    NIST U.S. DEPARTMENT OF COMMERCE PUBLICATIONS £ Technology Administration National Institute of Standards and Technology FIPS PUB 183 FEDERAL INFORMATION PROCESSING STANDARDS PUBLICATION INTEGRATION DEFINITION FOR FUNCTION MODELING (IDEFO) » Category: Software Standard SUBCATEGORY: MODELING TECHNIQUES 1993 December 21 183 PUB FIPS JK- 45C .AS A3 //I S3 IS 93 FIPS PUB 183 FEDERAL INFORMATION PROCESSING STANDARDS PUBLICATION INTEGRATION DEFINITION FOR FUNCTION MODELING (IDEFO) Category: Software Standard Subcategory: Modeling Techniques Computer Systems Laboratory National Institute of Standards and Technology Gaithersburg, MD 20899 Issued December 21, 1993 U.S. Department of Commerce Ronald H. Brown, Secretary Technology Administration Mary L. Good, Under Secretary for Technology National Institute of Standards and Technology Arati Prabhakar, Director Foreword The Federal Information Processing Standards Publication Series of the National Institute of Standards and Technology (NIST) is the official publication relating to standards and guidelines adopted and promulgated under the provisions of Section 111 (d) of the Federal Property and Administrative Services Act of 1949 as amended by the Computer Security Act of 1987, Public Law 100-235. These mandates have given the Secretary of Commerce and NIST important responsibilities for improving the utilization and management of computer and related telecommunications systems in the Federal Government. The NIST, through its Computer Systems Laboratory, provides leadership, technical guidance,
    [Show full text]