System simulation using Matlab, state plane plots
|
|
- Madison Porter
- 5 years ago
- Views:
Transcription
1 System simulation using Matlab, state plane plots his lab is mainly concerned with making state plane (also referred to as phase plane ) plots for various linear and nonlinear systems with two states he second part of this document presents the derivation of the differential equations which describe the dynamics of the inverted pendulum Problem ake for example the system described by the differential equation Denoting x = y, x = y we now write x = x x = α( x ) x x he equilibrium point is x =, x = he linear model of the system around this equilibrium point is x x x = α x y+ α( y ) y + y = he following Matlab code simulates the system and draws the resulting trajectories in the phase plane (state plane) function main clear all;close all;clc; nt=; % nt is the number of trajectories for k=:nt, x=[(rand(,)-5*[ ]')*]; % each one of the two initial states is %chosen randomly in the interval [-,+] =;tspan=[ ]; [t,x]= ode(@system,tspan,x); plot(x(),x(),'or',x(:,),x(:,)); end title('phase plane plot'); xlabel('x_'); ylabel('x_') function xdot=system(t,x) x=[x() x()]'; alpha=; xdot=[x(); -alpha*(x()^-)*x()-x()]; % xdot=[x(); - x()+alpha*x()]; % this is to be uncommented when the %simulation of the linear model of the system is desired a Run the simulation for various numbers of trajectories Are the trajectories moving toward or away from the equilibrium point? If they are moving toward (away of) the equilibrium point then it is a stable (unstable) equilibrium b Run the simulation for various values of α Notice how this parameter changes the shape of the trajectories
2 c Notice that the trajectories are converging to a closed orbit in the state space hat is a limit cycle Since the trajectories are converging towards it this means that it is a stable limit cycle Problem Consider the system described by the differential equation x = x y = x y a Calculate the equilibrium points bdetermine the linear model of the system around those equilibrium points c Simulate the two linear approximations of the nonlinear system and draw the phase plane plots One equilibrium point is stable while the other one is a saddle point he following figures present the results in the two cases Phase plane plot 5 x x Figure State plane trajectories for the linear model around the stable equilibrium (,) Notice that for the linear model of the system, valid around this equilibrium point, this equilibrium is moved to the point (,) Phase plane plot 5 x x Figure State plane trajectories for the linear model around the stable equilibrium (,-)
3 dsimulate the nonlinear system and draw the phase plane plot he next figure presents such a result Phase plane plot x x Figure State plane trajectories for the nonlinear system Notice how in small regions around the two equilibrium points the phase plot for the nonlinear system resembles the ones obtained at part c he state trajectories go towards the equilibrium (,) from all directions he equilibrium (,-) is avoided on the direction of x while the trajectories are attracted to it on the direction of x he Matlab code which was used to obtain the previous figures is the following: function main clear all;close all;clc; nt=6; % nt is the number of trajectories for k=:nt, x=[(rand()-5)*; (rand())*-]; % each one of the two initial states is chosen randomly in a given interval =;tspan=[ ]; [t,x]= ode(@system,tspan,x); plot(x(),x(),'or',x(:,),x(:,)); end title('phase plane plot'); xlabel('x_'); ylabel('x_') function xdot=system(t,x) x=[x() x()]'; %xdot=[-x(); -x()^-x()^]; xdot=[- ; ]*x; %saddle point equilibrium %xdot=[- ; -]*x; %stable equilibrium
4 Problem Consider the system described by the differential equations x = x+ xx x = x xx a Calculate the equilibrium points b Determine the linear model of the system around those equilibrium points c Simulate the two linear approximations of the nonlinear system and draw the phase plane plots d Simulate the nonlinear system and draw the phase plane plot Problem 4 Simulate the following systems and draw the phase plane plot x = y( + x y ) a y = x( + y x ) x = y b y = x r = r( r) c θ = sin θ + ( r) In the case of part b you will see that the trajectories in the state plane are circles In the case of part c, plot the state trajectories in the plane (x,y) where x= rcos θ, y = r You will see that the equilibrium point in the origin is unstable and that all the state trajectories are attracted to a circular trajectory with radius equal to one (a stable limit cycle) 4
5 INVERED PENDULUM he inverted pendulum on a cart is representative of a class of systems that includes stabilization of a rocket during launch he position of the cart is p, the angle of the rod is θ, the force input to the cart is f, the cart mass is M, the mass of the bob is m, and the length of the rod is L he coordinates of the bob are (p,z ) p m θ L z g f M p In the following we derive the differential equations which describe the dynamics of the inverted pendulum system use Lagrange's equation he kinetic energy of the cart is K = Mp he kinetic energy of the bob is Inverted Pendulum K = m( p + z ) where p = p + L, z = L so that p = p + Lθ, z = Lθ herefore, the total kinetic energy is K = K + K = Mp + m( p + p θ L + L θ ) he potential energy is due to the bob and is U = mgz = mgl he Lagrangian is (note that on the left hand side of this equation L denotes the Lagrangian function while L in the right hand side is the constant length of the pendulum arm) L = K U = ( M + m) p + mlp θ + ml θ mgl 5
6 he generalized coordinates are selected as q = [ q ] q = [ p θ ] so that Lagrange's equations are d L L = f dt p p d L L = dt θ θ Substituting for L and performing the partial differentiation yields ( M + m) p + mlθ mlθ = f mlp + ml θ mgl = hese dynamical equations must now be placed into state-space form o accomplish this, write the Lagrange equation in terms of matrices as M + m ml p mlθ + f = ml ml θ mgl his is a mechanical system in typical Lagrangian form, with the inertia matrix multiplying the acceleration vector he term θ ml is a centripetal term and mgl is a gravity term Invert the inertia matrix and simplify to obtain mg mlθ f p = m cos θ ( M + m) ( M + m) g + mlθ + f θ = ml cos θ ( M + m) L Now, the state may be defined as x = [ x x x x ] [ ] = p p θ θ 4 and the input as u= f hen the nonlinear state equation may be written as x mg sin x cos x mlx4 sin x u m cos x ( M + m) x = = x4 ( M + m) g sin x + mlx4 sin x cos x + u cos x ml cos x ( M + m) L f ( x, u) Given these nonlinear state equations, it is very easy to simulate the inverted pendulum behavior on a computer We now want to linearize this and obtain the linear state equations he equilibrium point of interest is x= p p θ θ = [ ], where the rod is upright Instead of using the linearization method, for simplicity in this case we will use approximations which are valid near the origin sin x x,cos x In addition, all squared state components are very small and so set equal to zero his yields the linear state equation 6
7 x mgx u + mg x M M x M = + u = Ax + Bu x = 4 ( M + m) g ( M + m) gx u ML ML ML he output equation depends on the measurements taken, which depends on the sensors available Assuming measurements of cart position and rod angle, the output equation is p y = = x = h( x, u) θ he cart position may be measured by placing an optical encoder on one of the wheels, and the rod angle by placing an encoder at the rod pivot point It is difficult to measure the velocities p, θ, but this might be achieved by placing tachometers on a wheel and at the rod pivot point hen, the output equation will change Given the linear state-space equations, a controller can be designed to keep the rod upright hough the controller is designed using the linear state equations, the performance of the controller should be simulated in a closed-loop system using the full nonlinear dynamics x = f ( x, u) Problem 5 Write the Matlab code for the simulation of the inverted pendulum Present the results of the simulation, considering the initial state [ ] x = and u =, using a visual representation of the inverted pendulum in a similar fashion as in the case of the robotic arm from the previous laboratory session Note that this will be part of a future mini-project which will require simulation of the inverted pendulum system in the case when the control input is calculated using an optimal state feedback control law 7
Lagrange s Equations of Motion and the Generalized Inertia
Lagrange s Equations of Motion and the Generalized Inertia The Generalized Inertia Consider the kinetic energy for a n degree of freedom mechanical system with coordinates q, q 2,... q n. If the system
More informationSimple and Physical Pendulums Challenge Problem Solutions
Simple and Physical Pendulums Challenge Problem Solutions Problem 1 Solutions: For this problem, the answers to parts a) through d) will rely on an analysis of the pendulum motion. There are two conventional
More informationMatlab-Based Tools for Analysis and Control of Inverted Pendula Systems
Matlab-Based Tools for Analysis and Control of Inverted Pendula Systems Slávka Jadlovská, Ján Sarnovský Dept. of Cybernetics and Artificial Intelligence, FEI TU of Košice, Slovak Republic sjadlovska@gmail.com,
More informationThe Inverted Pendulum
Lab 1 The Inverted Pendulum Lab Objective: We will set up the LQR optimal control problem for the inverted pendulum and compute the solution numerically. Think back to your childhood days when, for entertainment
More informationEE Homework 3 Due Date: 03 / 30 / Spring 2015
EE 476 - Homework 3 Due Date: 03 / 30 / 2015 Spring 2015 Exercise 1 (10 points). Consider the problem of two pulleys and a mass discussed in class. We solved a version of the problem where the mass was
More informationPhys 7221, Fall 2006: Midterm exam
Phys 7221, Fall 2006: Midterm exam October 20, 2006 Problem 1 (40 pts) Consider a spherical pendulum, a mass m attached to a rod of length l, as a constrained system with r = l, as shown in the figure.
More informationApplication of Neural Networks for Control of Inverted Pendulum
Application of Neural Networks for Control of Inverted Pendulum VALERI MLADENOV Department of Theoretical Electrical Engineering Technical University of Sofia Sofia, Kliment Ohridski blvd. 8; BULARIA valerim@tu-sofia.bg
More informationLab 6a: Pole Placement for the Inverted Pendulum
Lab 6a: Pole Placement for the Inverted Pendulum Idiot. Above her head was the only stable place in the cosmos, the only refuge from the damnation of the Panta Rei, and she guessed it was the Pendulum
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 informationTorque and Simple Harmonic Motion
Torque and Simple Harmonic Motion Recall: Fixed Axis Rotation Angle variable Angular velocity Angular acceleration Mass element Radius of orbit Kinematics!! " d# / dt! " d 2 # / dt 2!m i Moment of inertia
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 informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationEE 16B Midterm 2, March 21, Name: SID #: Discussion Section and TA: Lab Section and TA: Name of left neighbor: Name of right neighbor:
EE 16B Midterm 2, March 21, 2017 Name: SID #: Discussion Section and TA: Lab Section and TA: Name of left neighbor: Name of right neighbor: Important Instructions: Show your work. An answer without explanation
More informationThe dynamics of a Mobile Inverted Pendulum (MIP)
The dynamics of a Mobile Inverted Pendulum (MIP) 1 Introduction Saam Ostovari, Nick Morozovsky, Thomas Bewley UCSD Coordinated Robotics Lab In this document, a Mobile Inverted Pendulum (MIP) is a robotic
More informationProblem Solving Session 10 Simple Harmonic Oscillator Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics 8.01 Problem Solving Session 10 Simple Harmonic Oscillator Solutions W13D3-0 Group Problem Gravitational Simple Harmonic Oscillator Two identical
More informationClassical Mechanics. FIG. 1. Figure for (a), (b) and (c). FIG. 2. Figure for (d) and (e).
Classical Mechanics 1. Consider a cylindrically symmetric object with a total mass M and a finite radius R from the axis of symmetry as in the FIG. 1. FIG. 1. Figure for (a), (b) and (c). (a) Show that
More informationLab 4 Numerical simulation of a crane
Lab 4 Numerical simulation of a crane Agenda Time 10 min Item Review agenda Introduce the crane problem 95 min Lab activity I ll try to give you a 5- minute warning before the end of the lab period to
More informationInverted Pendulum. Objectives
Inverted Pendulum Objectives The objective of this lab is to experiment with the stabilization of an unstable system. The inverted pendulum problem is taken as an example and the animation program gives
More informationSRV02-Series Rotary Experiment # 7. Rotary Inverted Pendulum. Student Handout
SRV02-Series Rotary Experiment # 7 Rotary Inverted Pendulum Student Handout SRV02-Series Rotary Experiment # 7 Rotary Inverted Pendulum Student Handout 1. Objectives The objective in this experiment is
More informationDynamic Modeling of Rotary Double Inverted Pendulum Using Classical Mechanics
ISBN 978-93-84468-- Proceedings of 5 International Conference on Future Computational echnologies (ICFC'5) Singapore, March 9-3, 5, pp. 96-3 Dynamic Modeling of Rotary Double Inverted Pendulum Using Classical
More informationDesign and Comparison of Different Controllers to Stabilize a Rotary Inverted Pendulum
ISSN (Online): 347-3878, Impact Factor (5): 3.79 Design and Comparison of Different Controllers to Stabilize a Rotary Inverted Pendulum Kambhampati Tejaswi, Alluri Amarendra, Ganta Ramesh 3 M.Tech, Department
More information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More informationFigure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody
Lecture 27. THE COMPOUND PENDULUM Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody diagram The term compound is used to distinguish the present
More informationSTATE VARIABLE (SV) SYSTEMS
Copyright F.L. Lewis 999 All rights reserved Updated:Tuesday, August 05, 008 STATE VARIABLE (SV) SYSTEMS A natural description for dynamical systems is the nonlinear state-space or state variable (SV)
More informationHarmonic Oscillator - Model Systems
3_Model Systems HarmonicOscillators.nb Chapter 3 Harmonic Oscillator - Model Systems 3.1 Mass on a spring in a gravitation field a 0.5 3.1.1 Force Method The two forces on the mass are due to the spring,
More informationThe Acrobot and Cart-Pole
C H A P T E R 3 The Acrobot and Cart-Pole 3.1 INTRODUCTION A great deal of work in the control of underactuated systems has been done in the context of low-dimensional model systems. These model systems
More informationRotational Kinetic Energy
Lecture 17, Chapter 10: Rotational Energy and Angular Momentum 1 Rotational Kinetic Energy Consider a rigid body rotating with an angular velocity ω about an axis. Clearly every point in the rigid body
More informationSolution to phys101-t112-final Exam
Solution to phys101-t112-final Exam Q1. An 800-N man stands halfway up a 5.0-m long ladder of negligible weight. The base of the ladder is.0m from the wall as shown in Figure 1. Assuming that the wall-ladder
More informationAnalysis of Linear State Space Models
Max f (lb) Max y (in) 1 ME313 Homework #7 Analysis of Linear State Space Models Last Updated November 6, 213 Repeat the car-crash problem from HW#6. Use the Matlab function lsim to perform the simulation.
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationPhysics Mechanics. Lecture 32 Oscillations II
Physics 170 - Mechanics Lecture 32 Oscillations II Gravitational Potential Energy A plot of the gravitational potential energy U g looks like this: Energy Conservation Total mechanical energy of an object
More informationProblem 1. Mathematics of rotations
Problem 1. Mathematics of rotations (a) Show by algebraic means (i.e. no pictures) that the relationship between ω and is: φ, ψ, θ Feel free to use computer algebra. ω X = φ sin θ sin ψ + θ cos ψ (1) ω
More informationMEAM 510 Fall 2011 Bruce D. Kothmann
Balancing g Robot Control MEAM 510 Fall 2011 Bruce D. Kothmann Agenda Bruce s Controls Resume Simple Mechanics (Statics & Dynamics) of the Balancing Robot Basic Ideas About Feedback & Stability Effects
More informationTwo-Dimensional Rotational Kinematics
Two-Dimensional Rotational Kinematics Rigid Bodies A rigid body is an extended object in which the distance between any two points in the object is constant in time. Springs or human bodies are non-rigid
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More informationDifferent Modelling and Controlling Technique For Stabilization Of Inverted Pendulam
International Journal of Scientific & Engineering Research Volume 4, Issue 2, February-23 Different Modelling and Controlling Technique For Stabilization Of Inverted Pendulam K.CHAKRABORTY,, R.R. MUKHERJEE,
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationPHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010
Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationMCE 366 System Dynamics, Spring Problem Set 2. Solutions to Set 2
MCE 366 System Dynamics, Spring 2012 Problem Set 2 Reading: Chapter 2, Sections 2.3 and 2.4, Chapter 3, Sections 3.1 and 3.2 Problems: 2.22, 2.24, 2.26, 2.31, 3.4(a, b, d), 3.5 Solutions to Set 2 2.22
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 1 The concept of feedback The role of a controller What is a state? The concept of model in systems engineering: a
More informationTHE REACTION WHEEL PENDULUM
THE REACTION WHEEL PENDULUM By Ana Navarro Yu-Han Sun Final Report for ECE 486, Control Systems, Fall 2013 TA: Dan Soberal 16 December 2013 Thursday 3-6pm Contents 1. Introduction... 1 1.1 Sensors (Encoders)...
More informationMEAM 510 Fall 2012 Bruce D. Kothmann
Balancing g Robot Control MEAM 510 Fall 2012 Bruce D. Kothmann Agenda Bruce s Controls Resume Simple Mechanics (Statics & Dynamics) of the Balancing Robot Basic Ideas About Feedback & Stability Effects
More informationEl péndulo invertido: un banco de pruebas para el control no lineal. XXV Jornadas de Automática
El péndulo invertido: un banco de pruebas para el control no lineal Javier Aracil and Francisco Gordillo Escuela Superior de Ingenieros Universidad de Sevilla XXV Jornadas de Automática Ciudad Real, 8-1
More informationAP Physics. Harmonic Motion. Multiple Choice. Test E
AP Physics Harmonic Motion Multiple Choice Test E A 0.10-Kg block is attached to a spring, initially unstretched, of force constant k = 40 N m as shown below. The block is released from rest at t = 0 sec.
More informationStatic Equilibrium, Gravitation, Periodic Motion
This test covers static equilibrium, universal gravitation, and simple harmonic motion, with some problems requiring a knowledge of basic calculus. Part I. Multiple Choice 1. 60 A B 10 kg A mass of 10
More informationECE-320: Linear Control Systems Homework 8. 1) For one of the rectilinear systems in lab, I found the following state variable representations:
ECE-30: Linear Control Systems Homework 8 Due: Thursday May 6, 00 at the beginning of class ) For one of the rectilinear systems in lab, I found the following state variable representations: 0 0 q q+ 74.805.6469
More informationA Normal Form for Energy Shaping: Application to the Furuta Pendulum
Proc 4st IEEE Conf Decision and Control, A Normal Form for Energy Shaping: Application to the Furuta Pendulum Sujit Nair and Naomi Ehrich Leonard Department of Mechanical and Aerospace Engineering Princeton
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationROBOTICS Laboratory Problem 02
ROBOTICS 2015-2016 Laboratory Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with or without friction along the horizontal
More informationVersion 001 circular and gravitation holland (2383) 1
Version 00 circular and gravitation holland (383) This print-out should have 9 questions. Multiple-choice questions may continue on the next column or page find all choices before answering. AP B 993 MC
More informationNonlinear Oscillators: Free Response
20 Nonlinear Oscillators: Free Response Tools Used in Lab 20 Pendulums To the Instructor: This lab is just an introduction to the nonlinear phase portraits, but the connection between phase portraits and
More informationis conserved, calculating E both at θ = 0 and θ = π/2 we find that this happens for a value ω = ω given by: 2g
UNIVERSITETET I STAVANGER Institutt for matematikk og naturvitenskap Suggested solutions, FYS 500 Classical Mechanics Theory 2016 fall Set 5 for 23. September 2016 Problem 27: A string can only support
More informationClicker Quiz. a) 25.4 b) 37.9 c) 45.0 d) 57.1 e) 65.2
Clicker Quiz Assume that the rock is launched with an angle of θ = 45. With what angle with respect to the horizontal does the rock strike the ground in front of the castle? v 0 = 14.2 m/s v f = 18.5 m/s
More informationSwinging-Up and Stabilization Control Based on Natural Frequency for Pendulum Systems
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrC. Swinging-Up and Stabilization Control Based on Natural Frequency for Pendulum Systems Noriko Matsuda, Masaki Izutsu,
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationTime Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors
Time Response of Dynamic Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Multi-dimensional trajectories Numerical integration Linear and nonlinear systems Linearization
More informationTopic # Feedback Control Systems
Topic #19 16.31 Feedback Control Systems Stengel Chapter 6 Question: how well do the large gain and phase margins discussed for LQR map over to DOFB using LQR and LQE (called LQG)? Fall 2010 16.30/31 19
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More informationkx m x B N 1 C L, M Mg θ
.004 MODELING DYNAMICS AND CONTROL II Spring 00 Solutions to Problem Set No. 7 Problem 1. Pendulum mounted on elastic support. This problem is an execise in the application of momentum principles. Two
More informationInteractions Between Two Non-Stationary Pendulums
Interactions Between Two Non-Stationary Pendulums Alexander Rich Harvey Mudd College 3 December 2013 Abstract Should two pendulums on a frictionless cart synchronize? This experiment measures the angular
More informationUnforced Oscillations
Unforced Oscillations Simple Harmonic Motion Hooke s Law Newton s Second Law Method of Force Competition Visualization of Harmonic Motion Phase-Amplitude Conversion The Simple Pendulum and The Linearized
More informationFALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Saturday, 14 December 2013, 1PM to 4 PM, AT 1003
FALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Saturday, 14 December 2013, 1PM to 4 PM, AT 1003 NAME: STUDENT ID: INSTRUCTION 1. This exam booklet has 14 pages. Make sure none are missing 2. There is
More informationReal-Time Implementation of a LQR-Based Controller for the Stabilization of a Double Inverted Pendulum
Proceedings of the International MultiConference of Engineers and Computer Scientists 017 Vol I,, March 15-17, 017, Hong Kong Real-Time Implementation of a LQR-Based Controller for the Stabilization of
More informationINPUT-STATE LINEARIZATION OF A ROTARY INVERTED PENDULUM
0 Asian Journal of Control Vol 6 No pp 0-5 March 004 Brief Paper INPU-SAE LINEARIZAION OF A ROARY INVERED PENDULUM Chih-Keng Chen Chih-Jer Lin and Liang-Chun Yao ABSRAC he aim of this paper is to design
More informationNumerics and Control of PDEs Lecture 1. IFCAM IISc Bangalore
1/1 Numerics and Control of PDEs Lecture 1 IFCAM IISc Bangalore July 22 August 2, 2013 Introduction to feedback stabilization Stabilizability of F.D.S. Mythily R., Praveen C., Jean-Pierre R. 2/1 Q1. Controllability.
More informationθ + mgl θ = 0 or θ + ω 2 θ = 0 (2) ω 2 = I θ = mgl sinθ (1) + Ml 2 I = I CM mgl Kater s Pendulum The Compound Pendulum
Kater s Pendulum The Compound Pendulum A compound pendulum is the term that generally refers to an arbitrary lamina that is allowed to oscillate about a point located some distance from the lamina s center
More informationTable of Contents 1.0 OBJECTIVE APPARATUS PROCEDURE LAB PREP WORK POLE-PLACEMENT DESIGN
LAB 4 ENGI 38: ADVANCED CONTROLS -------------------------------------------------------------------- Lab Partners: (Alphabetically) Figliomeni, Dan Malyshev, Andrey McGrath, Adam TO: WARREN PAJU ELECTRICAL
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationLinear control of inverted pendulum
Linear control of inverted pendulum Deep Ray, Ritesh Kumar, Praveen. C, Mythily Ramaswamy, J.-P. Raymond IFCAM Summer School on Numerics and Control of PDE 22 July - 2 August 213 IISc, Bangalore http://praveen.cfdlab.net/teaching/control213
More information. d. v A v B. e. none of these.
General Physics I Exam 3 - Chs. 7,8,9 - Momentum, Rotation, Equilibrium Oct. 28, 2009 Name Rec. Instr. Rec. Time For full credit, make your work clear to the grader. Show the formulas you use, the essential
More informationA. B. C. D. E. v x. ΣF x
Q4.3 The graph to the right shows the velocity of an object as a function of time. Which of the graphs below best shows the net force versus time for this object? 0 v x t ΣF x ΣF x ΣF x ΣF x ΣF x 0 t 0
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationApr 14, Calculus with Algebra and Trigonometry II Lecture 20More physics Aprapplications
Calculus with Algebra and Trigonometry II Lecture 20 More physics applications Apr 14, 2015 14, 2015 1 / 14 Motion in two dimensions A particle s motion can be described by specifying how the coordinates
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationpaths 1, 2 and 3 respectively in the gravitational field of a point mass m,
58. particles of mass m is moving in a circular path of constant radius r such that its centripetal acceleration a c is varying with time t as a c = k 2 rt 2 where k is a constant. The power delivered
More informationFeedback Optimal Control for Inverted Pendulum Problem by Using the Generating Function Technique
(IJACSA) International Journal o Advanced Computer Science Applications Vol. 5 No. 11 14 Feedback Optimal Control or Inverted Pendulum Problem b Using the Generating Function echnique Han R. Dwidar Astronom
More informationNeural Network Control of an Inverted Pendulum on a Cart
Neural Network Control of an Inverted Pendulum on a Cart VALERI MLADENOV, GEORGI TSENOV, LAMBROS EKONOMOU, NICHOLAS HARKIOLAKIS, PANAGIOTIS KARAMPELAS Department of Theoretical Electrical Engineering Technical
More informationReverse Order Swing-up Control of Serial Double Inverted Pendulums
Reverse Order Swing-up Control of Serial Double Inverted Pendulums T.Henmi, M.Deng, A.Inoue, N.Ueki and Y.Hirashima Okayama University, 3-1-1, Tsushima-Naka, Okayama, Japan inoue@suri.sys.okayama-u.ac.jp
More informationChapter 14. Oscillations. Oscillations Introductory Terminology Simple Harmonic Motion:
Chapter 14 Oscillations Oscillations Introductory Terminology Simple Harmonic Motion: Kinematics Energy Examples of Simple Harmonic Oscillators Damped and Forced Oscillations. Resonance. Periodic Motion
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationSWINGING UP A PENDULUM BY ENERGY CONTROL
Paper presented at IFAC 13th World Congress, San Francisco, California, 1996 SWINGING UP A PENDULUM BY ENERGY CONTROL K. J. Åström and K. Furuta Department of Automatic Control Lund Institute of Technology,
More informationFundamentals Physics. Chapter 15 Oscillations
Fundamentals Physics Tenth Edition Halliday Chapter 15 Oscillations 15-1 Simple Harmonic Motion (1 of 20) Learning Objectives 15.01 Distinguish simple harmonic motion from other types of periodic motion.
More informationLagging Pendulum. Martin Ga
02 Lagging Pendulum Martin Ga 1 Task A pendulum consists of a strong thread and a bob. When the pivot of the pendulum starts moving along a horizontal circumference, the bob starts tracing a circle which
More informationROTATIONAL DOUBLE INVERTED PENDULUM
ROTATIONAL DOUBLE INVERTED PENDULUM Thesis Submitted to The School of Engineering of the UNIVERSITY OF DAYTON In Partial Fulfillment of the Requirements for The Degree of Master of Science in Electrical
More information= o + t = ot + ½ t 2 = o + 2
Chapters 8-9 Rotational Kinematics and Dynamics Rotational motion Rotational motion refers to the motion of an object or system that spins about an axis. The axis of rotation is the line about which the
More informationMSMS Matlab Problem 02
MSMS 2014-2015 Matlab Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with friction along the horizontal rail; the cart supports
More information2.003 Quiz #1 Review
2.003J Spring 2011: Dynamics and Control I Quiz #1 Review Massachusetts Institute of Technology March 5th, 2011 Department of Mechanical Engineering March 6th, 2011 1 Reference Frames 2.003 Quiz #1 Review
More informationKINETIC ENERGY SHAPING IN THE INVERTED PENDULUM
KINETIC ENERGY SHAPING IN THE INVERTED PENDULUM J. Aracil J.A. Acosta F. Gordillo Escuela Superior de Ingenieros Universidad de Sevilla Camino de los Descubrimientos s/n 49 - Sevilla, Spain email:{aracil,
More informationConservation of Energy
Name Period Date Conservation of Energy Driving Questions How does the energy of a cart poised at the top of a hill compare to its energy at the bottom of the hill? Background Gravitational potential energy
More informationMechatronic System Case Study: Rotary Inverted Pendulum Dynamic System Investigation
Mechatronic System Case Study: Rotary Inverted Pendulum Dynamic System Investigation Dr. Kevin Craig Greenheck Chair in Engineering Design & Professor of Mechanical Engineering Marquette University K.
More information14.4 Energy in Simple Harmonic Motion 14.5 Pendulum Motion.notebook January 25, 2018
The interplay between kinetic and potential energy is very important for understanding simple harmonic motion. Section 14.4 Energy in Simple Harmonic Motion For a mass on a spring, when the object is at
More informationPHYS 1114, Lecture 33, April 10 Contents:
PHYS 1114, Lecture 33, April 10 Contents: 1 This class is o cially cancelled, and has been replaced by the common exam Tuesday, April 11, 5:30 PM. A review and Q&A session is scheduled instead during class
More informationMEM04: Rotary Inverted Pendulum
MEM4: Rotary Inverted Pendulum Interdisciplinary Automatic Controls Laboratory - ME/ECE/CHE 389 April 8, 7 Contents Overview. Configure ELVIS and DC Motor................................ Goals..............................................3
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationPHY 141 Midterm 1 Oct 2, 2014 Version A
PHY 141 Midterm 1 Oct 2, 2014 Version A Put FULL NAME, ID#, and EXAM VERSION on the front cover of the BLUE BOOKLET! To avoid problems in grading: do all problems in order, write legibly, and show all
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationLecture 10. Example: Friction and Motion
Lecture 10 Goals: Exploit Newton s 3 rd Law in problems with friction Employ Newton s Laws in 2D problems with circular motion Assignment: HW5, (Chapter 7, due 2/24, Wednesday) For Tuesday: Finish reading
More informationCourse Name: AP Physics. Team Names: Jon Collins. Velocity Acceleration Displacement
Course Name: AP Physics Team Names: Jon Collins 1 st 9 weeks Objectives Vocabulary 1. NEWTONIAN MECHANICS and lab skills: Kinematics (including vectors, vector algebra, components of vectors, coordinate
More information