3 Modeling and Simulation
|
|
- Howard Armstrong
- 5 years ago
- Views:
Transcription
1 ME 132, Spring 2005, UC Berkeley, A Packard 13 3 Modeling and Simulation 31 Systems of 1st order, coupled differential equations Consider an input/output system, with m inputs (denoted by d), q outputs (denoted by e), and governed by a set of n, 1st order, coupled differential equations, of the form ẋ 1 (t) ẋ 2 (t) ẋ n (t) e 1 (t) e 2 (t) e q (t) = f 1 (x 1 (t), x 2 (t),, x n (t), d 1 (t), d 2 (t),, d m (t), t) f 2 (x 1 (t), x 2 (t),, x n (t), d 1 (t), d 2 (t),, d m (t), t) f n (x 1 (t), x 2 (t),, x n (t), d 1 (t), d 2 (t),, d m (t), t) h 1 (x 1 (t), x 2 (t),, x n (t), d 1 (t), d 2 (t),, d m (t), t) h 2 (x 1 (t), x 2 (t),, x n (t), d 1 (t), d 2 (t),, d m (t), t) h q (x 1 (t), x 2 (t),, x n (t), d 1 (t), d 2 (t),, d m (t), t) where the functions f i, h i are given functions of the n variables x 1, x 2,, x n, the m variables d 1, d 2,, d m and also explicit functions of t For shorthand, we write (79) as (8) ẋ(t) = f (x(t), d(t), t) e(t) = h (x(t), d(t), t) (9) Given an initial condition vector x 0, and a forcing function d(t) for t t 0, we wish to solve for the solutions x 1 (t) e 1 (t) x 2 (t) e 2 (t) x(t) = x n (t), e(t) = on the interval [t 0, t F ], given the initial condition and the input forcing function d( ) x (t 0 ) = x 0 e q (t) simulink solves for this using numerical integration techniques, such as 4th and 5th order Runge-Kutta formulae You can learn more about numerical integration by taking Math 128 We will not discuss this important topic in detail in this class However, to understand in a very crude manner how ODE solvers work, consider the Euler method of solution If the functions f and u are reasonably well behaved in x and t, then the solution exists, is continuous, and differentiable at all points Hence, it is reasonable that a Taylor series for x at a given time t will be predictive of the values of x around t
2 ME 132, Spring 2005, UC Berkeley, A Packard 14 If we do a Taylor s expansion on a function x, and ignore the higher order terms, we get an approximation formula x(t + δ) x(t) + δẋ(t) = x(t) + δf(x(t), d(t), t) Roughly, the smaller δ is, the closer the left-hand-side is to the actual value of x(t + δ) Euler s method propogates a solution to (79) by using this approximation repeatedly for a fixed δ, called the stepsize Hence, Euler s method gives that for any integer k >= 0, the solution to (79) approximately satisfies x((k + 1)δ) = x(kδ) + δ f(x(kδ), d(kδ), kδ) Writing out the first 4 time steps (ie, t = 0, δ, 2δ, 3δ, 4δ) gives x(δ) = x(0) + δf (x(0), d(0), 0) x(2δ) = x(δ) + δf (x(δ), d(δ), δ) x(3δ) = x(2δ) + δf (x(2δ), d(2δ), 2δ) x(4δ) = x(3δ) + δf (x(3δ), d(3δ), 3δ) (10) and so on So, as long as you have a subroutine that can evaluate f(x, d, t), given x and t, you can quickly propogate an approximate solution simply by calling the subroutine once for every timestep Computing the output, e(t) simply involves evaluating the function h(x(t), d(t), t) at the solution points In the Runge-Kutta method, a more sophisticated approximation is made, which results in more computations (4 function evaluations of f for every time step), but much greater accuracy In effect, more terms of the Taylor series are used, involving matrices of partial derivatives, and even their derivatives, df dx, d2 f dx 2, d3 f dx 3 but without actually requiring explicit knowledge of these derivatives of the function f 32 Remarks about Integration Options in simulink The Simulation SimulationParameters Solver page is used to set additional optional properties, including integration step-size options, and may need to be used to obtain smooth plots Additional options are
3 ME 132, Spring 2005, UC Berkeley, A Packard 15 Term RelTol AbsTol MaxStep Meaning Relative error tolerance, default 1e-3, probably leave it alone, though if instructions below don t work, try making it a bit smaller Absolute error tolerance, default 1e-6, probably leave it alone, though if instructions below don t work, try making it a bit smaller maximum step size I believe the default is (StopTime- StartTime)/50 In general, make it smaller than (StopTime-StartTime)/50 if your plots are jagged 33 Problems 1 A sealed-box loudspeaker is shown below From Amplifier Pole Piece Voice Coil Vent Magnet Back Plate Top Plate Spider Voice Coil Former Basket Subwoofer Electrical Leads Dust Cap Woofer Cone Surround Signal From Accelerometer Sealed Enclosure Ignore the wire marked accelerometer mounted on cone We will develop a model for accoustic radiation from a sealed-box loudspeaker The equations are Speaker Cone: force balance, Voice-coil motor (a DC motor): Magnetic flux/length Factor: m z(t) = F vc (t) + F k (t) + F d (t) + F b + F env (t) V in (t) LI(t) RI(t) B l (t)ż(t) = 0 F vc (t) = B l (t)i(t) B l (t) = BL BL 1 z(t) 4
4 ME 132, Spring 2005, UC Berkeley, A Packard 16 Suspension/Surround: Sealed Enclosure: F k (t) = K 0 z(t) K 1 z 2 (t) K 2 z 3 (t) F d (t) = R S ż(t) F b (t) = P 0 A Environment (Baffled Half-Space) where x(t) is 6 1, and A e := ( 1 + Az(t) V 0 ) γ ẋ(t) = A e x(t) + B e ż(t) F env (t) = AP 0 C e x(t) D e ż(t) B e := C e := [ ] D e = 4662 This is an approximate model of the impedance seen by the face of the loudspeaker as it radiates into an infinite half-space The values for the various constants are
5 ME 132, Spring 2005, UC Berkeley, A Packard 17 Symbol Value A meters 2 V meters 3 P Pa m 0117 kg L H R 3 Ω BL Tesla meters BL meters 4 K N/meter K 1 0 K N/meter 3 R S 128 N sec/meter (a) Build a Simulink model for the system Use the Subsystem grouping capability to manage the complexity of the diagram Each of the equations above should represent a different subsystem Make sure you think through the question what is the input(s) and what is the output(s)? for each subsystem (b) Write a function that has two input arguments (denoted V and Ω) and two output arguments, z max,pos and z max,neg The functional relationship is defined as follows: Suppose V in (t) = V sin Ω 2πt, where t is in seconds, and V in volts In other words, the input is a sin-wave at a frequency of Ω Hertz Simulate the loudspeaker behavior for about 25/Ω seconds The displacement z of the cone will become nearly sinusoidal by the end of the simulation Let z max,pos and z max,neg be the maximum positive and negative values of the displacement in the last full cycle (ie, the last 1/Ω seconds of the simulation, which we will approximately decree as the steady-state response to the input ) (c) Using bisection, determine the value (approximately) of V so that the steadystate maximum (positive) excursion of z is 0007meters if the frequency of the excitation is Ω = 25 What is the steady-state minimum (negative) excursion of z?
Contents. 0. Introduction Loudspeaker s Impedance Why is this useful? How good is it? Bibliography...
Contents 0. Introduction... 2... 2 1. Loudspeaker s Impedance... 2 2. Why is this useful?... 10 3. How good is it?... 13 4. Bibliography... 14 0. Introduction In this article I present a way to model a
More informationNonlinear Losses in Electro-acoustical Transducers Wolfgang Klippel, Daniel Knobloch
The Association of Loudspeaker Manufacturers & Acoustics International (ALMA) Nonlinear Losses in Electro-acoustical Transducers Wolfgang Klippel, Daniel Knobloch Institute of Acoustics and Speech Communication
More information24, B = 59 24, A = 55
Math 128a - Homework 8 - Due May 2 1) Problem 8.4.4 (Page 555) Solution: As discussed in the text, the fourth-order Adams-Bashforth formula is a formula of the type x n+1 = x n + h[af n + Bf n 1 + Cf n
More informationPRO 5W-8. Professional series. SPL vs Freq. Specification. Materials of Construction
Professional series PRO 5W-8 Small in size, large in performance. Featuring an advanced X5 pulp cone, high temperature copper wire voice coil, and truncated cast frame chassis, the Pro 5W can handle your
More informationNumerical Algorithms for ODEs/DAEs (Transient Analysis)
Numerical Algorithms for ODEs/DAEs (Transient Analysis) Slide 1 Solving Differential Equation Systems d q ( x(t)) + f (x(t)) + b(t) = 0 dt DAEs: many types of solutions useful DC steady state: state no
More informationPhysically Based Modeling Differential Equation Basics
Physically Based Modeling Differential Equation Basics Andrew Witkin and David Baraff Pixar Animation Studios Please note: This document is 2001 by Andrew Witkin and David Baraff. This chapter may be freely
More informationFunctional differentiation
Functional differentiation March 22, 2016 1 Functions vs. functionals What distinguishes a functional such as the action S [x (t] from a function f (x (t, is that f (x (t is a number for each value of
More informationPhysics 106a/196a Problem Set 7 Due Dec 2, 2005
Physics 06a/96a Problem Set 7 Due Dec, 005 Version 3, Nov 7, 005 In this set we finish up the SHO and study coupled oscillations/normal modes and waves. Problems,, and 3 are for 06a students only, 4, 5,
More informationPhysically Based Modeling: Principles and Practice Differential Equation Basics
Physically Based Modeling: Principles and Practice Differential Equation Basics Andrew Witkin and David Baraff Robotics Institute Carnegie Mellon University Please note: This document is 1997 by Andrew
More informationChecking the Radioactive Decay Euler Algorithm
Lecture 2: Checking Numerical Results Review of the first example: radioactive decay The radioactive decay equation dn/dt = N τ has a well known solution in terms of the initial number of nuclei present
More informationKlippel Non-Linear Test Results. LSI (Large Signal Identification) Model #: PS Introduction. Large Signal Modeling. Nonlinear Characteristics
Klippel Non-Linear Test Results LSI (Large Signal Identification) Model #: PS180-8 Introduction Large Signal Modeling At higher amplitudes, loudspeakers produce substantial distortion in the output signal,
More informationKlippel Non-Linear Test Results LSI (Large Signal Identification) Driver Name: ND90-8. Introduction. Nonlinear Parameters. Large Signal Modeling
Klippel Non-Linear Test Results LSI (Large Signal Identification) Driver Name: ND90-8 Introduction Large Signal Modeling At higher amplitudes, loudspeakers produce substantial distortion in the output
More informationCHAPTER 5 SIMULATION AND TEST SETUP FOR FAULT ANALYSIS
47 CHAPTER 5 SIMULATION AND TEST SETUP FOR FAULT ANALYSIS 5.1 INTRODUCTION This chapter describes the simulation model and experimental set up used for the fault analysis. For the simulation set up, the
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationδ Substituting into the differential equation gives: x i+1 x i δ f(t i,x i ) (3)
Solving Differential Equations Numerically Differential equations are ubiquitous in Physics since the laws of nature often take on a simple form when expressed in terms of infinitesimal changes of the
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationPh 22.1 Return of the ODEs: higher-order methods
Ph 22.1 Return of the ODEs: higher-order methods -v20130111- Introduction This week we are going to build on the experience that you gathered in the Ph20, and program more advanced (and accurate!) solvers
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. For each of the following linear systems, determine whether the origin is asymptotically stable and, if so, plot the step response and frequency response for the system. If there are multiple
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationf = 1 T 6 a.c. (Alternating Current) Circuits Most signals of interest in electronics are periodic : they repeat regularly as a function of time.
Analogue Electronics (Aero).66 66 Analogue Electronics (Aero) 6.66 6 a.c. (Alternating Current) Circuits Most signals of interest in electronics are periodic : they repeat regularly as a function of time.
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationNew Trends in Modeling and Identification of Loudspeaker with Nonlinear Distortion
New Trends in Modeling and Identification of Loudspeaker with Nonlinear Distortion Pascal Brunet and Bahram Shafai Department of Electrical and Computer Engineering, Northeastern University, Boston, MA,
More informationUsing Simulink to analyze 2 degrees of freedom system
Using Simulink to analyze 2 degrees of freedom system Nasser M. Abbasi Spring 29 page compiled on June 29, 25 at 4:2pm Abstract A two degrees of freedom system consisting of two masses connected by springs
More informationCHARGE TO MASS RATIO FOR THE ELECTRON
CHARGE TO MASS RATIO FOR THE ELECTRON OBJECTIVE: To measure the ratio of the charge of an electron to its mass. METHOD: A stream of electrons is accelerated by having them "fall" through a measured potential
More informationImpedance and Loudspeaker Parameter Measurement
ECEN 2260 Circuits/Electronics 2 Spring 2007 2-26-07 P. Mathys Impedance and Loudspeaker Parameter Measurement 1 Impedance Measurement Many elements from which electrical circuits are built are two-terminal
More informationAccurate Determination of Loudspeaker Parameters using Audio Analyzer Type 2012 and Laser Velocity Transducer Type 3544
7/6-'89 Application Note Accurate Determination of Loudspeaker Parameters using Audio Analyzer Type 202 and Laser Velocity Transducer Type 3544 by Søren Jønsson, Brüel & Kjær, Denmark A method to determine
More informationEE 16B Final, December 13, Name: SID #:
EE 16B Final, December 13, 2016 Name: SID #: Important Instructions: Show your work. An answer without explanation is not acceptable and does not guarantee any credit. Only the front pages will be scanned
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationSolved Problems. Electric Circuits & Components. 1-1 Write the KVL equation for the circuit shown.
Solved Problems Electric Circuits & Components 1-1 Write the KVL equation for the circuit shown. 1-2 Write the KCL equation for the principal node shown. 1-2A In the DC circuit given in Fig. 1, find (i)
More informationKnud Thorborg Scan-Speak, Videbæk, Denmark,
Knud Thorborg Scan-Speak, Videbæk, Denmark, kt@scan-speak.dk Traditional and Advanced Models for the Dynamic Loudspeaker The traditional equivalent circuit for a loudspeaker, based on the so-called Thiele-Small
More informationLab 1: Dynamic Simulation Using Simulink and Matlab
Lab 1: Dynamic Simulation Using Simulink and Matlab Objectives In this lab you will learn how to use a program called Simulink to simulate dynamic systems. Simulink runs under Matlab and uses block diagrams
More informationCurrent in a Magnetic Field Learning Outcomes. Force on a Current-Carrying Conductor
1 Current in a Magnetic Field Learning Outcomes Discuss the force on a current-carrying conductor in a magnetic field. Demonstrate this force. Solve problems about this force. Discuss applications of this
More informationECE2262 Electric Circuits. Chapter 6: Capacitance and Inductance
ECE2262 Electric Circuits Chapter 6: Capacitance and Inductance Capacitors Inductors Capacitor and Inductor Combinations Op-Amp Integrator and Op-Amp Differentiator 1 CAPACITANCE AND INDUCTANCE Introduces
More informationQuiescent Steady State (DC) Analysis The Newton-Raphson Method
Quiescent Steady State (DC) Analysis The Newton-Raphson Method J. Roychowdhury, University of California at Berkeley Slide 1 Solving the System's DAEs DAEs: many types of solutions useful DC steady state:
More informationMOTORS AND GENERATORS
DO PHYSCS ONLNE MOTORS AND GENERATORS view 1 Charge q Q [coulomb C] view 2 Current i [ampere A] view 3 Potential difference v V [volt V] Electric ield E [V.m -1 N.C -1 ] view 4 Resistance R [ohm ] view
More informationA First Course on Kinetics and Reaction Engineering Supplemental Unit S5. Solving Initial Value Differential Equations
Supplemental Unit S5. Solving Initial Value Differential Equations Defining the Problem This supplemental unit describes how to solve a set of initial value ordinary differential equations (ODEs) numerically.
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
More informationCoupled Electrical Oscillators Physics Advanced Physics Lab - Summer 2018 Don Heiman, Northeastern University, 5/24/2018
Coupled Electrical Oscillators Physics 3600 - Advanced Physics Lab - Summer 08 Don Heiman, Northeastern University, 5/4/08 I. INTRODUCTION The objectives of this experiment are: () explore the properties
More informationMechatronic design optimisation of subwoofer system. P6-project spring 2010 Group MCE6-621
Mechatronic design optimisation of subwoofer system P6-project spring 2010 Group MCE6-621 Title: Mechatronic design optimisation of subwoofer system Semester: 6th semester 2010 Semester theme: Mechatronic
More informationAPPLICATION OF THE FINITE ELEMENT METHOD TO MODEL THE NONLINEAR VOICE COIL MOTION PRODUCED BY A LOUDSPEAKER MAGNET ASSEMBLY.
APPLICATION OF THE FINITE ELEMENT METHOD TO MODEL THE NONLINEAR VOICE COIL MOTION PRODUCED BY A LOUDSPEAKER MAGNET ASSEMBLY. Mark Dodd Celestion International & KEF Audio (UK) Ltd. 1. INTRODUCTION Moving
More information4.7 Magnetism and electromagnetism
4.7 Magnetism and electromagnetism Electromagnetic effects are used in a wide variety of devices. Engineers make use of the fact that a magnet moving in a coil can produce electric current and also that
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationDepartment of Architecture & Civil Engineering ( ) 2 2a. L = 65 2 ρπa4 L. + asinα = 3aθ 2. ( ) = a 1 cos( θ ρπa4 L.
MODE ANSWER age: 1 QUESTION Mass of tube = ρπ 3a ( ) ( a) Moment of inertia of tube = ρπ 3 Mass of bar = ρπa Moment of inertia of bar = = 5ρπa ( 3a) 4 a ( ) 4 ρπ ( a )4 = 1 3 ρπa4 Horizontal displacement
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationSection 6.0 : Advanced Transmission Line Modeling Techniques By Martin J King, 07/05/02 Copyright 2002 by Martin J. King. All Rights Reserved.
Section 6.0 : Advanced Transmission Line Modeling Techniques In Section 4.0, two simple models were derived for a closed end and an open end transmission line. In both of these models the driver is located
More informationEuler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability.
Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability. REVIEW: We start with the differential equation dy(t) dt = f (t, y(t)) (1.1) y(0) = y 0 This equation can be
More informationDifferential Equations
Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationLecture 17: Ordinary Differential Equation II. First Order (continued)
Lecture 17: Ordinary Differential Equation II. First Order (continued) 1. Key points Maple commands dsolve dsolve[interactive] dsolve(numeric) 2. Linear first order ODE: y' = q(x) - p(x) y In general,
More information7 PDES, FOURIER SERIES AND TRANSFORMS Last Updated: July 16, 2012
Problem List 7.1 Fourier series expansion of absolute value function 7.2 Fourier series expansion of step function 7.3 Orthogonality of functions 7.4 Fourier transform of sinusoidal pulse 7.5 Introducing
More informationOften, in this class, we will analyze a closed-loop feedback control system, and end up with an equation of the form
ME 32, Spring 25, UC Berkeley, A. Packard 55 7 Review of SLODEs Throughout this section, if y denotes a function (of time, say), then y [k or y (k) denotes the k th derivative of the function y, y [k =
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
More informationYell if you have any questions
Class 31: Outline Hour 1: Concept Review / Overview PRS Questions possible exam questions Hour : Sample Exam Yell if you have any questions P31 1 Exam 3 Topics Faraday s Law Self Inductance Energy Stored
More informationFor those of you who are taking Calculus AB concurrently with AP Physics, I have developed a
AP Physics C: Mechanics Greetings, For those of you who are taking Calculus AB concurrently with AP Physics, I have developed a brief introduction to Calculus that gives you an operational knowledge of
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (35 points) Leslie Leroy Irvin bails out of an airplane at the altitude of 16,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s
More informationFile: KappaPro15SmallVentedBox500Watts.bb6
KappaPro15A, Small Vented Box, For Mid/High use only. By McJerry, Eminence Speaker LLC Displacement limited to 3 watts if run down to 1 Hz. Displacement and thermally limited to 5 watts if run down to
More informationEE C128 / ME C134 Midterm Fall 2014
EE C128 / ME C134 Midterm Fall 2014 October 16, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket calculator
More informationMath 240: Spring/Mass Systems II
Math 240: Spring/Mass Systems II Ryan Blair University of Pennsylvania Monday, March 26, 2012 Ryan Blair (U Penn) Math 240: Spring/Mass Systems II Monday, March 26, 2012 1 / 12 Outline 1 Today s Goals
More informationENGG4420 LECTURE 7. CHAPTER 1 BY RADU MURESAN Page 1. September :29 PM
CHAPTER 1 BY RADU MURESAN Page 1 ENGG4420 LECTURE 7 September 21 10 2:29 PM MODELS OF ELECTRIC CIRCUITS Electric circuits contain sources of electric voltage and current and other electronic elements such
More informationCoordinate systems and vectors in three spatial dimensions
PHYS2796 Introduction to Modern Physics (Spring 2015) Notes on Mathematics Prerequisites Jim Napolitano, Department of Physics, Temple University January 7, 2015 This is a brief summary of material on
More informationPHYS 1441 Section 001 Lecture #23 Monday, Dec. 4, 2017
PHYS 1441 Section 1 Lecture #3 Monday, Dec. 4, 17 Chapter 3: Inductance Mutual and Self Inductance Energy Stored in Magnetic Field Alternating Current and AC Circuits AC Circuit W/ LRC Chapter 31: Maxwell
More informationET3-7: Modelling II(V) Electrical, Mechanical and Thermal Systems
ET3-7: Modelling II(V) Electrical, Mechanical and Thermal Systems Agenda of the Day 1. Resume of lesson I 2. Basic system models. 3. Models of basic electrical system elements 4. Application of Matlab/Simulink
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationAP Physics C - E & M
Slide 1 / 27 Slide 2 / 27 AP Physics C - E & M Current, Resistance & Electromotive Force 2015-12-05 www.njctl.org Slide 3 / 27 Electric Current Electric Current is defined as the movement of charge from
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationNumerical Analysis Exam with Solutions
Numerical Analysis Exam with Solutions Richard T. Bumby Fall 000 June 13, 001 You are expected to have books, notes and calculators available, but computers of telephones are not to be used during the
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5. Numerical Methods Gaëtan Kerschen Space Structures & Systems Lab (S3L) Why Different Propagators? Analytic propagation: Better understanding of the perturbing forces. Useful
More informationAssume that you have made n different measurements of a quantity x. Usually the results of these measurements will vary; call them x 1
#1 $ http://www.physics.fsu.edu/users/ng/courses/phy2048c/lab/appendixi/app1.htm Appendix I: Estimates for the Reliability of Measurements In any measurement there is always some error or uncertainty in
More information3.4 Introduction to power series
3.4 Introduction to power series Definition 3.4.. A polynomial in the variable x is an expression of the form n a i x i = a 0 + a x + a 2 x 2 + + a n x n + a n x n i=0 or a n x n + a n x n + + a 2 x 2
More informationChapter 21. Magnetic Forces and Magnetic Fields
Chapter 21 Magnetic Forces and Magnetic Fields 21.1 Magnetic Fields The needle of a compass is permanent magnet that has a north magnetic pole (N) at one end and a south magnetic pole (S) at the other.
More informationMulti Acoustic Prediction Program (MAPP tm ) Recent Results Perrin S. Meyer and John D. Meyer
Multi Acoustic Prediction Program (MAPP tm ) Recent Results Perrin S. Meyer and John D. Meyer Meyer Sound Laboratories Inc., Berkeley, California, USA Presented at the Institute of Acoustics (UK), Reproduced
More informationElectrodynamics Exam 3 and Final Exam Sample Exam Problems Dr. Colton, Fall 2016
Electrodynamics Exam 3 and Final Exam Sample Exam Problems Dr. Colton, Fall 016 Multiple choice conceptual questions 1. An infinitely long, straight wire carrying current passes through the center of a
More informationECE2262 Electric Circuits. Chapter 6: Capacitance and Inductance
ECE2262 Electric Circuits Chapter 6: Capacitance and Inductance Capacitors Inductors Capacitor and Inductor Combinations 1 CAPACITANCE AND INDUCTANCE Introduces two passive, energy storing devices: Capacitors
More informationMotors and Generators
Motors and Generators Joe Khachan and Chris Stewart Introduction This part of the syllabus draws upon many ideas and concepts from the old HSC syllabus, perhaps with a bit more of a historical twist thrown
More informationOrdinary Differential Equations
Ordinary Differential Equations In this lecture, we will look at different options for coding simple differential equations. Start by considering bicycle riding as an example. Why does a bicycle move forward?
More informationExam 3 Topics. Displacement Current Poynting Vector. Faraday s Law Self Inductance. Circuits. Energy Stored in Inductor/Magnetic Field
Exam 3 Topics Faraday s Law Self Inductance Energy Stored in Inductor/Magnetic Field Circuits LR Circuits Undriven (R)LC Circuits Driven RLC Circuits Displacement Current Poynting Vector NO: B Materials,
More informationCopyright 2014 by Martin J. King. All Rights Reserved. Configuration : Extended Range Driver w/ Bass Driver Mounted on an Open Baffle
Open Baffle w/ 3 Drivers - Acoustic Response 12/3/1 Software : by Martin J. King e-mail MJKing57@aol.com Copyright 21 by Martin J. King. All Rights Reserved. Configuration : Extended Range Driver w/ Bass
More informationOpen Baffle w/ 2 Drivers - Acoustic Response 12/30/14. Copyright 2014 by Martin J. King. All Rights Reserved.
Open Baffle w/ 2 Drivers - Acoustic Response 12/3/14 Software : by Martin J. King e-mail MJKing57@aol.com Copyright 214 by Martin J. King. All Rights Reserved. Configuration : Extended Range Driver w/
More informationMATH 31B: MIDTERM 2 REVIEW. sin 2 x = 1 cos(2x) dx = x 2 sin(2x) 4. + C = x 2. dx = x sin(2x) + C = x sin x cos x
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate sin x and cos x. Solution: Recall the identities cos x = + cos(x) Using these formulas gives cos(x) sin x =. Trigonometric Integrals = x sin(x) sin x = cos(x)
More informationInitial Value Problems
Numerical Analysis, lecture 13: Initial Value Problems (textbook sections 10.1-4, 10.7) differential equations standard form existence & uniqueness y 0 y 2 solution methods x 0 x 1 h h x 2 y1 Euler, Heun,
More informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationK&J Magnetics Mendocino Motor Worksheet
K&J Magnetics Mendocino Motor Worksheet This report summarizes some theoretical work we did to answer these questions: What solar panels should we use? What gauge of magnet wire should we use? How many
More informationIntroduction to Series and Sequences Math 121 Calculus II Spring 2015
Introduction to Series and Sequences Math Calculus II Spring 05 The goal. The main purpose of our study of series and sequences is to understand power series. A power series is like a polynomial of infinite
More informationGeneral Physics I. Lecture 14: Sinusoidal Waves. Prof. WAN, Xin ( 万歆 )
General Physics I Lecture 14: Sinusoidal Waves Prof. WAN, Xin ( 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Motivation When analyzing a linear medium that is, one in which the restoring force
More informationSection 7.0 : Design of a Back Loaded Exponential Horn By Martin J. King, 08/17/04 Copyright 2004 by Martin J. King. All Rights Reserved.
Section 7. : Design of a Back Loaded Exponential Horn By Martin J. King, 8/7/4 Copyright 4 by Martin J. King. All Rights Reserved. Section 7. : Design of a Back Loaded Exponential Horn For over a year,
More informationMath 252 Fall 2002 Supplement on Euler s Method
Math 5 Fall 00 Supplement on Euler s Method Introduction. The textbook seems overly enthusiastic about Euler s method. These notes aim to present a more realistic treatment of the value of the method and
More informationFeedback Basics. David M. Auslander Mechanical Engineering University of California at Berkeley. copyright 1998, D.M. Auslander
Feedback Basics David M. Auslander Mechanical Engineering University of California at Berkeley copyright 1998, D.M. Auslander 1 I. Feedback Control Context 2 What is Feedback Control? Measure desired behavior
More informationNonlinear Modeling of the Heat Transfer in Loudspeakers
Nonlinear Modeling of the Heat Transfer in Loudspeakers Wolfgang Klippel Klippel GmbH Dresden, 01277, Germany www.klippel.de ABSTRACT Traditional modeling describes the heat flow in loudspeakers by an
More informationCompare the measured behavior with desired behavior, and
Chapter Introduction The objective of this text is to introduce methods for design and analysis of classical feedback control systems. The basic principle of feedback is to: Use a sensor to measure the
More informationDON T PANIC! If you get stuck, take a deep breath and go on to the next question. Come back to the question you left if you have time at the end.
Math 307, Midterm 2 Winter 2013 Name: Instructions. DON T PANIC! If you get stuck, take a deep breath and go on to the next question. Come back to the question you left if you have time at the end. There
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationMATH 105: PRACTICE PROBLEMS FOR CHAPTER 3: SPRING 2010
MATH 105: PRACTICE PROBLEMS FOR CHAPTER 3: SPRING 010 INSTRUCTOR: STEVEN MILLER (SJM1@WILLIAMS.EDU Question 1 : Compute the partial derivatives of order 1 order for: (1 f(x, y, z = e x+y cos(x sin(y. Solution:
More informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationCHAPTER 3 : MATHEMATICAL MODELLING PRINCIPLES
CHAPTER 3 : MATHEMATICAL MODELLING PRINCIPLES When I complete this chapter, I want to be able to do the following. Formulate dynamic models based on fundamental balances Solve simple first-order linear
More informationME 132, Dynamic Systems and Feedback. Class Notes. Spring Instructor: Prof. A Packard
ME 132, Dynamic Systems and Feedback Class Notes by Andrew Packard, Kameshwar Poolla & Roberto Horowitz Spring 2005 Instructor: Prof. A Packard Department of Mechanical Engineering University of California
More information