Seismologia ja maan rakenne A Seismology and structure of the Earth
|
|
- Gyles Manning
- 5 years ago
- Views:
Transcription
1 Seismologia ja maan rakenne A Seismology and structure of the Earth Prac-cals 1: Seismic data analysis: poles and zeros, convolu-on and deconvolu-on Main materials: 1) New Manual of Seismological Observatory Prac-ce (NMSOP- 2): hjp:// bib.telegrafenberg.de/publizieren/vertrieb/nmsop/ hjp://geoladm.geol.queensu.ca Chapter 5. Erhardt Wielandt. Seismic sensors and their calibra-on Chapter 11. Peter Bormann, Klaus Klinge and Siegfried Wendt. Data Analysis and Seismogram Interpreta-on 2) Havskow, J., Alguacil, G Instrumenta-on in earthquake seismology FOR POLES AND ZEROS YOU CAN FIND MATERIALS IN CONTROL SYSTEMS IN INTERNET, ALSO IN YOUTUBE: hjps://
2 The basic requirements in analog and digital rou-ne observatory prac-ce i.e., to: recognize the occurrence of an earthquake in a record; iden-fy and annotate the seismic phases; determine onset -me and polarity correctly; measure the maximum ground amplitude and related period; calculate slowness and azimuth; determine source parameters such as the hypocenter, origin -me, magnitude, source mechanism, etc.. In modern digital observatory prac-ce these procedures are implemented in computer programs.
3 Seismic observables Period ranges (order of magnitudes) Sound s Earthquakes s (surface waves, body waves) Eigenmodes of the Earth 1000 s Coseismic deforma-on 1 s 1000 s Postseismic deforma-on s Seismic explora-on s Laboratory signals s s - > What are the consequences for sampling intervals, data volumes, etc.? Modern Seismology Data processing and inversion 3
4 Different factors/sub- systems (without seismic noise) which influence a seismic record (yellow boxes) and the informa-on that can be derived from record analysis (blue boxes)
5 SEISMOMETERS Output voltage is propordonal to ground velocity in certain frequency bandwidth. requirements ü Very low noise ü ResoluDon ü Broad bandwidth Weight and size MagneDc resistance Price
6 BASIC DEFINITIONS 1) Signals and systems A signal is a descrip-on of how one parameter varies with another parameter. For instance, voltage changing over -me in an electronic circuit, or brightness varying with distance in an image. A system is any process that produces an output signal in response to an input signal Con-nuous systems input and output con-nuous signals, such as in analog electronics. Discrete systems input and output discrete signals, such as computer programs that manipulate the values stored in arrays. x(t) x[n] Continuous System Discrete System y(t) y[n] Con-nuous signals are usually represented with parentheses, while discrete signals use brackets. All signals use lower case lejers, reserving the upper case for the frequency domain. Unless there is a bejer name available, the input signal is called: x(t) or x[n], while the output is called: y(t) or y[n]
7 The Laplace transformadon and transfer funcdon of a linear system A signal that has a definite beginning in -me (such as the seismic waves from an earthquake) can be decomposed into exponen-ally growing, sta-onary, or exponen-ally decaying sinusoidal signals with the Laplace integral transforma-on: st ( s) f ( t) e dt 0 F. f ( t) 1 2 j j j F( s) e st ds s=σ + jω Generally, the transfer func-on H(s) of an linear system is the complex func-on for which G( s) H ( s) F( s) with F(s) and G(s) represen-ng the Laplace transforms of the input and output signals
8 GENERAL PRESENTATION OF A TRANSFER FUNCTION OF A LINEAR SYSTEM: where B and A and are polynomials in s, M is the order of the numerator polynomial, b m is the m- th coefficient of the numerator polynomial, N is the order of the denominator polynomial, and a n is the n- th coefficient of the denominator polynomial. The zeros of the system are roots of the numerator polynomial: TF - à 0 The poles of the system are roots of the denominator polynomial: TF - à inf
9 It is oren convenient to factor the polynomials in the numerator and denominator, and to write the transfer func-on in terms of those factors: on in terms of those factors: H(s) = N(s) D(s) = K (s z 1)(s z 2 )...(s z m 1 )(s z m ) (s p 1 )(s p 2 )...(s p n 1 )(s p n ), where the numerator and denominator polynomials, N(s) and D(s), have real coefficients defined by the system s differential equation and K = b m /a n. As written in Eq. (2) the z i s are the roots of the equation N(s) =0, (3) and are defined to be the system zeros, and the p i s are the roots of the equation D(s) =0, (4) and are defined to be the system poles. In Eq. (2) the factors in the numerator and denominator are written so that when s = z i the numerator N(s) = 0 and the transfer function vanishes, that is lim H(s) =0. s z i
10 Pole- zero plot IMAGINARY REAL
11 System poles define stability of a system: system is stable when f(t)à 0 if tà inf The specifica-on of the form of components of the homogeneous response from the system pole loca-ons on the pole- zero plot: EXAMPLE: H(s)=1/(s+2) à s=- 2, if we take inverce Laplace transform, then f(t)= e - 2 H(s)=1/(s- 2) e 2: The à specification s=2, then off(t)=e the form 2 of components of the homogeneous response from the s
12 e 2: The specification of the form of components of the homogeneous response from the s 1. A real pole p i = σ in the left-half of the s-plane defines an exponentially decaying component, Ce σt, in the homogeneous response. The rate of the decay is determined by the pole location; poles far from the origin in the left-half plane correspond to components that decay rapidly, while poles near the origin correspond to slowly decaying components. 2. A pole at the origin p i = 0 defines a component that is constant in amplitude and defined by the initial conditions. 3. A real pole in the right-half plane corresponds to an exponentially increasing component Ce σt in the homogeneous response; thus defining the system to be unstable. 4. A complex conjugate pole pair σ ± jω in the left-half of the s-plane combine to generate a response component that is a decaying sinusoid of the form Ae σt sin (ωt + φ) where A and φ are determined by the initial conditions. The rate of decay is specified by σ; thefrequency of oscillation is determined by ω. 5. An imaginary pole pair, that is a pole pair lying on the imaginary axis, ±jω generates an oscillatory component with a constant amplitude determined by the initial conditions. 6. A complex pole pair in the right half plane generates an exponentially increasing component.
13 Example A linear system is described by the differential equation Find the system poles and zeros. Solution: d 2 y dt 2 +5dy dt +6y =2du dt +1. From the differential equation the transfer function is H(s) = which may be written in factored form 2s +1 s 2 +5s +6. (5) H(s) = 1 s +1/2 2 (s + 3)(s + 2) = 1 s ( 1/2) 2 (s ( 3))(s ( 2)). (6) The system therefore has a single real zero at s = 1/2, and a pair of real poles at s = 3 and s = 2.
14 Seismic instruments: The standard iner-al seismometer with damping Spring force Damping force It is convenient to define a resonant angular frequency and a dumping parameter as follows: These subs-tu-ons give: This equa-on shows that the Earth accelera-on may be recovered by measuring the displacement of the mass and its -me deriva-ves.
15 EXERCICE&1:&write&the&expression&for&transfer&func8on& Example A linear system is described by the differential equation Find the system poles and zeros. Solution: d 2 y dt 2 +5dy dt +6y =2du dt +1. From the differential equation the transfer function is H(s) = which may be written in factored form 2s +1 s 2 +5s +6. (5) H(s) = 1 s +1/2 2 (s + 3)(s + 2) = 1 s ( 1/2) 2 (s ( 3))(s ( 2)). (6) The system therefore has a single real zero at s = 1/2, and a pair of real poles at s = 3 and s = 2.
16 Exercice 1: A second- order system has a pair of complex conjugate poles a s = 2±j3 and a single zero at the origin of the s- plane. Find the transfer func-on and use the pole- zero plot
17 Solution: From the problem description H(s) = s K (s ( 2+j3))(s ( 2 j3)) = s K s 2 +4s + 13 The pole-zero plot is shown in Fig. 6. From the figure the transfer function is
18 Poles and zeros must either be real or symmetric to the real axis, as men-oned above. When the numerator polynomial is s m, then s = 0 is an m- fold zero of the transfer func-on, and the system is a high- pass filter of order m. (Zeros at nonzero frequency do normally not appear in the transfer func-on of broadband seismographs because, if they occur mathema-cally, their effect must prac-cally be cancelled by nearby poles; otherwise the response would not be called broadband.) Adding a zero at s = 0 produces the transfer func-on Hz(s) = sh(s), and the step response of this system is purely the deriva-ve of the step response of the original system. Depending on the order n of the denominator and accordingly on the number of poles, the response may be flat at high frequencies (n = m), or the system may act as a low- pass filter there (n > m). The case n < m can occur only as an approxima-on in a limited bandwidth because no prac-cal system can have an unlimited gain at high frequencies.
19 Seismic instruments: The standard iner-al seismometer We obtained: and: When h=1, the system is said to be cri-cally damped. Seismometers generally perform op-mally at values of damping close to cri-cal. A polarity reversal at high frequencies.
20 EXERSICE 2: Plot poles and zeros for broadband seismometers Seismograph Zeros Poles STS2 (0.0, 0.0) (0.0, 0.0) (0.0, 0.0) STS1(GRF) (0.0, 0.0) (0.0, 0.0) (0.0, 0.0) STS1(VBB)) (0.0, 0.0) (0.0, 0.0) (0.0, 0.0) (-3.674E-2, E-3) (-3.674E-2, 3.675E-3) ( , ) ( , ) ( , 0.0) ( , 4.574) ( , ) (-7.006, ) (-7.006, ) ( , ) ( , ) ( E-02, E-02) ( E-02, E-02)
21 on in terms of those factors: H(s) = N(s) D(s) = K (s z 1)(s z 2 )...(s z m 1 )(s z m ) (s p 1 )(s p 2 )...(s p n 1 )(s p n ), The transfer function may be evaluated for any value of s = σ + jω, and in general, when s is complex the function H(s) itself is complex. It is common to express the complex value of the transfer function in polar form as a magnitude and an angle: with a magnitude H(s) and an angle φ(s) given by R {} I {} H(s) = H(s) e jφ(s), (17) H(s) = R {H(s)} 2 + I {H(s)} 2, (18) ( ) I {H(s)} φ(s) = tan 1 (19) R {H(s)}
22 The frequency response may be written in terms of the system poles and zeros by substituting jω for s directly into the factored form of the transfer function: H(jω) =K (jω z 1)(jω z 2 )...(jω z m 1 )(jω z m ) (jω p 1 )(jω p 2 )...(jω p n 1 )(jω p n ). (30)
23 EXERCISE 3: Convolu-on and deconvolu-on in seismic signal processing
24 ConvoluDon W Input displacement Output displacement Amplitude Input Output Convolu-on (Kanasewich 1981) is a mathema-cal opera-on defining the change of shape of a waveform resul-ng from its passage through a filter. Thus, for example, a seismic pulse generated by an explosion is altered in shape by filtering effects, both in the ground and in the recording system, so that the seismogram (the filtered output) differs significantly from the ini-al seismic pulse (the input). Time Fig. 2.9 The principle of filtering illustrated by the perturbation of a suspended weight system. As a simple example of filtering, consider a weight suspended from the end of a ver-cal spring. If the top of the spring is perturbed by a sharp up- and- down movement (the input), the mo-on of the weight (the filtered output) is a series of damped oscilla-ons out of phase with the ini-al perturba-on.
25 The effect of a filter is described mathema-cally by a convolu-on opera-on such that, if the input signal g(t) to the filter is convolved with the impulse response f(t) of the filter, known as the convolu-on operator, the filtered output y(t) is obtained (asterisk denotes a convolu-on operator): yt ()= gt ()* f() t Spike input Output = Impulse response Filter
26 The mathema-cal implementa-on of convolu-on involves -me inversion (or folding) of one of the func-ons and its progressive sliding past the other func-on, the individual terms in the convolved output being de- rived by summa-on of the cross- mul-plica-on products over the overlapping parts of the two func-ons. In general, if g i (i = 1, 2,..., m) is an input func-on and f j ( j = 1, 2,..., n) is a convolu-on operator, then the convolu-on output func-on y k is given by y = g f k = 12,,..., m+ n-1 k m  i= 1 i k-i ( )
27 In Fig the individual steps in the convolu-on process are shown for two digital func-ons, a double spike func-on given by g i = g 1,g 2,g 3 = 2,0,1 and an impulse response func-on given by f i = f 1, f 2, f 3, f 4 = 4, 3, 2, 1, where the numbers refer to discrete amplitude values at the sampling points of the two func-ons. From Fig it can be seen that the convolved output y i = y 1, y 2, y 3, y 4, y 5, y 6 = 8, 6, 8, 5, 2, 1. Note that the convolved output is longer than the input waveforms; if the func- -ons to be convolved have lengths of m and n, the convolved output has a length of (m + n - 1) Cross-products Sum = = = = = 2 Fig A method of calculating the convolution of two digital functions = 1
28 Deconvolu-on or inverse filtering (Kanasewich 1981) is a process that counteracts a previous convolu-on (or filtering) ac-on. Consider the convolu-on opera-on given in equa-on yt gt f t ()= () () * y(t) is the filtered output derived by passing the input waveform g(t) through a filter of impulse response f(t). Knowing y(t) and f(t), the recovery of g(t) represents a de- convolu-on opera-on. In the seismic case, y(t) is the seismic record resul-ng from the passage of a seismic wave g(t) through a por-on of the Earth, which acts as a filter with an impulse response f(t).
29 Reflectivity and convolution The seismic wave is sensitive to the sequence of impedance contrasts The reflectivity series (R) We input a source wavelet (W) which is reflected at each impedance contrast The seismogram recorded at the surface (S) is the convolution of the two S = W * R
30 Convolution Reflectivity series 1 ½ ½ Reflec-vity series describes the geological medium with reflec-vity contrasts Source wavelet ½ -½ 1 This is the seismic wave propaga-ng through the medium Output 0 Recorded waveform TASK: Calculate and plot the recorded seismic waveform using convolu-on operator
31 Spiking deconvolution Recorded waveform 1-1 ¾ -½ Deconvolution operator ¼ 1 1 Output 0 1 Recovered reflectivity series TASK: recover the reflec-vity series using spiking deconvolu-on operator
32 The Fourier transformadon, spectra, amplitude and phase response of a linear system ~ F( ) f ( t) e j t dt f 1 ~ j ( ) ( ) t t F e d, 2 The signal is here assumed to have a finite energy so that the integrals converge. The condition that no signal is present at negative times can be dropped in this case. The Fourier j t e transformation decomposes the signal into purely harmonic (sinusoidal) waves. The direct and inverse Fourier transformation are also known as a harmonic analysis and synthesis. We may consider the Fourier transformation as a special version of the Laplace transformation for real frequencies, i.e. for s = jω. In fact, F (ω ) = F ( jω ), i.e. the Fourier transform for real angular frequencies ω is identical to the Laplace transform for imaginary s = jω. In the case of s=jω the function H(ω) is called the complex frequency response of the system. The absolute value of H (ω ) is called the amplitude response, and the phase of H (ω ) the phase response of the system. The distinction between H(ω) and H(s) is essential when systems are characterized by their poles and zeros. Usually, poles and zeros are given for the Laplace transform.
33
34
Review of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 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
More information1 The frequency response of the basic mechanical oscillator
Seismograph systems The frequency response of the basic mechanical oscillator Most seismographic systems are based on a simple mechanical oscillator really just a mass suspended by a spring with some method
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 informationDisplacement at very low frequencies produces very low accelerations since:
SEISMOLOGY The ability to do earthquake location and calculate magnitude immediately brings us into two basic requirement of instrumentation: Keeping accurate time and determining the frequency dependent
More informationLecture 27 Frequency Response 2
Lecture 27 Frequency Response 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/12 1 Application of Ideal Filters Suppose we can generate a square wave with a fundamental period
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationDynamic circuits: Frequency domain analysis
Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution
More informationChapter 6: The Laplace Transform. Chih-Wei Liu
Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace
More informationAdvanced Analog Building Blocks. Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc
Advanced Analog Building Blocks Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc 1 Topics 1. S domain and Laplace Transform Zeros and Poles 2. Basic and Advanced current
More informationThe Laplace Transform
The Laplace Transform Syllabus ECE 316, Spring 2015 Final Grades Homework (6 problems per week): 25% Exams (midterm and final): 50% (25:25) Random Quiz: 25% Textbook M. Roberts, Signals and Systems, 2nd
More informationSIGNAL PROCESSING. B14 Option 4 lectures. Stephen Roberts
SIGNAL PROCESSING B14 Option 4 lectures Stephen Roberts Recommended texts Lynn. An introduction to the analysis and processing of signals. Macmillan. Oppenhein & Shafer. Digital signal processing. Prentice
More informationIntroduction & Laplace Transforms Lectures 1 & 2
Introduction & Lectures 1 & 2, Professor Department of Electrical and Computer Engineering Colorado State University Fall 2016 Control System Definition of a Control System Group of components that collectively
More informationTransform Representation of Signals
C H A P T E R 3 Transform Representation of Signals and LTI Systems As you have seen in your prior studies of signals and systems, and as emphasized in the review in Chapter 2, transforms play a central
More informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationANALOG AND DIGITAL SIGNAL PROCESSING CHAPTER 3 : LINEAR SYSTEM RESPONSE (GENERAL CASE)
3. Linear System Response (general case) 3. INTRODUCTION In chapter 2, we determined that : a) If the system is linear (or operate in a linear domain) b) If the input signal can be assumed as periodic
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationEE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models
EE/ME/AE324: Dynamical Systems Chapter 7: Transform Solutions of Linear Models The Laplace Transform Converts systems or signals from the real time domain, e.g., functions of the real variable t, to the
More informationThe Laplace Transform
The Laplace Transform Generalizing the Fourier Transform The CTFT expresses a time-domain signal as a linear combination of complex sinusoids of the form e jωt. In the generalization of the CTFT to the
More information2.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
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationQuestion Paper Code : AEC11T02
Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationLecture 28 Continuous-Time Fourier Transform 2
Lecture 28 Continuous-Time Fourier Transform 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/14 1 Limit of the Fourier Series Rewrite (11.9) and (11.10) as As, the fundamental
More informationExplanations and Discussion of Some Laplace Methods: Transfer Functions and Frequency Response. Y(s) = b 1
Engs 22 p. 1 Explanations Discussion of Some Laplace Methods: Transfer Functions Frequency Response I. Anatomy of Differential Equations in the S-Domain Parts of the s-domain solution. We will consider
More informationUnit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace
Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationFrequency Response. Re ve jφ e jωt ( ) where v is the amplitude and φ is the phase of the sinusoidal signal v(t). ve jφ
27 Frequency Response Before starting, review phasor analysis, Bode plots... Key concept: small-signal models for amplifiers are linear and therefore, cosines and sines are solutions of the linear differential
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationDESIGN OF CMOS ANALOG INTEGRATED CIRCUITS
DESIGN OF CMOS ANALOG INEGRAED CIRCUIS Franco Maloberti Integrated Microsistems Laboratory University of Pavia Discrete ime Signal Processing F. Maloberti: Design of CMOS Analog Integrated Circuits Discrete
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationDynamic Behavior. Chapter 5
1 Dynamic Behavior In analyzing process dynamic and process control systems, it is important to know how the process responds to changes in the process inputs. A number of standard types of input changes
More informationDr. Ian R. Manchester
Dr 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
More informationCircuit Analysis Using Fourier and Laplace Transforms
EE2015: Electrical Circuits and Networks Nagendra Krishnapura https://wwweeiitmacin/ nagendra/ Department of Electrical Engineering Indian Institute of Technology, Madras Chennai, 600036, India July-November
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
More informationDetailed Solutions to Exercises
Detailed Solutions to Exercises Digital Signal Processing Mikael Swartling Nedelko Grbic rev. 205 Department of Electrical and Information Technology Lund University Detailed solution to problem E3.4 A
More informationLaplace Transform Part 1: Introduction (I&N Chap 13)
Laplace Transform Part 1: Introduction (I&N Chap 13) Definition of the L.T. L.T. of Singularity Functions L.T. Pairs Properties of the L.T. Inverse L.T. Convolution IVT(initial value theorem) & FVT (final
More informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
More informationLECTURE 12 Sections Introduction to the Fourier series of periodic signals
Signals and Systems I Wednesday, February 11, 29 LECURE 12 Sections 3.1-3.3 Introduction to the Fourier series of periodic signals Chapter 3: Fourier Series of periodic signals 3. Introduction 3.1 Historical
More informationStability. X(s) Y(s) = (s + 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = +2 and s = -2. Y(s) 8X(s) G 1 G 2
Stability 8X(s) X(s) Y(s) = (s 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = 2 and s = -2 If all poles are in region where s < 0, system is stable in Fourier language s = jω G 0 - x3 x7 Y(s)
More informationR-L-C Circuits and Resonant Circuits
P517/617 Lec4, P1 R-L-C Circuits and Resonant Circuits Consider the following RLC series circuit What's R? Simplest way to solve for is to use voltage divider equation in complex notation. X L X C in 0
More informationSome of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e
Transform methods Some of the different forms of a signal, obtained by transformations, are shown in the figure. X(s) X(t) L - L F - F jw s s jw X(jw) X*(t) F - F X*(jw) jwt e z jwt z e X(nT) Z - Z X(z)
More informationTransform Solutions to LTI Systems Part 3
Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised
More informationDiscrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture
Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},
More informationTransient Response of a Second-Order System
Transient Response of a Second-Order System ECEN 830 Spring 01 1. Introduction In connection with this experiment, you are selecting the gains in your feedback loop to obtain a well-behaved closed-loop
More informationControl of Manufacturing Processes
Control of Manufacturing Processes Subject 2.830 Spring 2004 Lecture #18 Basic Control Loop Analysis" April 15, 2004 Revisit Temperature Control Problem τ dy dt + y = u τ = time constant = gain y ss =
More information27. The pole diagram and the Laplace transform
124 27. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More information3.2 Complex Sinusoids and Frequency Response of LTI Systems
3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex
More informationQ. 1 Q. 25 carry one mark each.
Q. Q. 5 carry one mark each. Q. Consider a system of linear equations: x y 3z =, x 3y 4z =, and x 4y 6 z = k. The value of k for which the system has infinitely many solutions is. Q. A function 3 = is
More informationLaplace Transforms and use in Automatic Control
Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationESS Finite Impulse Response Filters and the Z-transform
9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of
More informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More informationHow to manipulate Frequencies in Discrete-time Domain? Two Main Approaches
How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches Difference Equations (an LTI system) x[n]: input, y[n]: output That is, building a system that maes use of the current and previous
More informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
More information9.5 The Transfer Function
Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 3 Signals & Systems Prof. Mark Fowler Note Set #9 C-T Systems: Laplace Transform Transfer Function Reading Assignment: Section 6.5 of Kamen and Heck /7 Course Flow Diagram The arrows here show conceptual
More information20. The pole diagram and the Laplace transform
95 0. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More informationEarthquake Detec-on and Loca-on using Cascadia Ini-a-ve data. Susan L. Bilek and Emily Morton Earth and Environmental Science Dept.
Earthquake Detec-on and Loca-on using Cascadia Ini-a-ve data Susan L. Bilek and Emily Morton Earth and Environmental Science Dept. New Mexico Tech Mo-va-on Use unique Cascadia Ini-a-ve dataset to explore
More informationEA2.3 - Electronics 2 1
In the previous lecture, I talked about the idea of complex frequency s, where s = σ + jω. Using such concept of complex frequency allows us to analyse signals and systems with better generality. In this
More informationOverview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros)
Overview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros) J. McNames Portland State University ECE 222 Bode Plots Ver.
More informationLinear systems, small signals, and integrators
Linear systems, small signals, and integrators CNS WS05 Class Giacomo Indiveri Institute of Neuroinformatics University ETH Zurich Zurich, December 2005 Outline 1 Linear Systems Crash Course Linear Time-Invariant
More informationLecture 4: R-L-C Circuits and Resonant Circuits
Lecture 4: R-L-C Circuits and Resonant Circuits RLC series circuit: What's V R? Simplest way to solve for V is to use voltage divider equation in complex notation: V X L X C V R = in R R + X C + X L L
More informationSignals and Systems: Introduction
Dependent variable Signals and Systems: Introduction What is a signal? Signals may describe a wide variety of physical phenomena. The information in a signal is contained in a pattern of variations of
More information6.003: Signals and Systems Lecture 18 April 13, 2010
6.003: Signals and Systems Lecture 8 April 3, 200 6.003: Signals and Systems Discrete-Time Frequency Representations Signals and/or Systems Two perspectives: feedback and control (focus on systems) X +
More informationLinear Systems Theory
ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace
More informationCHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationBasic Electronics. Introductory Lecture Course for. Technology and Instrumentation in Particle Physics Chicago, Illinois June 9-14, 2011
Basic Electronics Introductory Lecture Course for Technology and Instrumentation in Particle Physics 2011 Chicago, Illinois June 9-14, 2011 Presented By Gary Drake Argonne National Laboratory Session 2
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
More informationChapter 3: Second Order ODE 3.8 Elements of Particle Dy
Chapter 3: Second Order ODE 3.8 Elements of Particle Dynamics 3 March 2018 Objective The objective of this section is to explain that any second degree linear ODE represents the motion of a particle. This
More informationSeismometer. Stefanie Donner
Seismometer Stefanie Donner Some history instrument for measuring the seasonal winds and the movements of the Earth Zhang Heng, China 132 AD #2 Some history 1755 the Great Lisbon earthquake Subsequent
More informationModule 4. Related web links and videos. 1. FT and ZT
Module 4 Laplace transforms, ROC, rational systems, Z transform, properties of LT and ZT, rational functions, system properties from ROC, inverse transforms Related web links and videos Sl no Web link
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationControl of Manufacturing Processes
Control of Manufacturing Processes Subject 2.830 Spring 2004 Lecture #19 Position Control and Root Locus Analysis" April 22, 2004 The Position Servo Problem, reference position NC Control Robots Injection
More informationCHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME
CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME Shri Mata Vaishno Devi University, (SMVDU), 2013 Page 13 CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME When characterizing or modeling a random variable, estimates
More informationCHAPTER 5: LAPLACE TRANSFORMS
CHAPTER 5: LAPLACE TRANSFORMS SAMANTHA RAMIREZ PREVIEW QUESTIONS What are some commonly recurring functions in dynamic systems and their Laplace transforms? How can Laplace transforms be used to solve
More informationEE Homework 13 - Solutions
EE3054 - Homework 3 - Solutions. (a) The Laplace transform of e t u(t) is s+. The pole of the Laplace transform is at which lies in the left half plane. Hence, the Fourier transform is simply the Laplace
More informationDEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010
[E2.5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010 EEE/ISE PART II MEng. BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: 2:00 hours There are FOUR
More informationThe Laplace Transform
The Laplace Transform Introduction There are two common approaches to the developing and understanding the Laplace transform It can be viewed as a generalization of the CTFT to include some signals with
More informationAbsolute strain determination from a calibrated seismic field experiment
Absolute strain determination Absolute strain determination from a calibrated seismic field experiment David W. Eaton, Adam Pidlisecky, Robert J. Ferguson and Kevin W. Hall ABSTRACT The concepts of displacement
More informationCast of Characters. Some Symbols, Functions, and Variables Used in the Book
Page 1 of 6 Cast of Characters Some s, Functions, and Variables Used in the Book Digital Signal Processing and the Microcontroller by Dale Grover and John R. Deller ISBN 0-13-081348-6 Prentice Hall, 1998
More informationSIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB
INTRODUCTION SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB Laplace transform pairs are very useful tools for solving ordinary differential equations. Most
More information8 sin 3 V. For the circuit given, determine the voltage v for all time t. Assume that no energy is stored in the circuit before t = 0.
For the circuit given, determine the voltage v for all time t. Assume that no energy is stored in the circuit before t = 0. Spring 2015, Exam #5, Problem #1 4t Answer: e tut 8 sin 3 V 1 For the circuit
More informationHomework 4. May An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt
Homework 4 May 2017 1. An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt Determine the impulse response of the system. Rewriting as y(t) = t e (t
More informationECE-314 Fall 2012 Review Questions for Midterm Examination II
ECE-314 Fall 2012 Review Questions for Midterm Examination II First, make sure you study all the problems and their solutions from homework sets 4-7. Then work on the following additional problems. Problem
More informationCore Concepts Review. Orthogonality of Complex Sinusoids Consider two (possibly non-harmonic) complex sinusoids
Overview of Continuous-Time Fourier Transform Topics Definition Compare & contrast with Laplace transform Conditions for existence Relationship to LTI systems Examples Ideal lowpass filters Relationship
More informationPoles and Zeros in z-plane
M58 Mixed Signal Processors page of 6 Poles and Zeros in z-plane z-plane Response of discrete-time system (i.e. digital filter at a particular frequency ω is determined by the distance between its poles
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationLecture 9 Infinite Impulse Response Filters
Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9
More informationForced Mechanical Vibrations
Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State
More informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
More information