How to create Matlab Script and Simulink Model for Designing a Pole Placement Controller
|
|
- Lindsey Arnold
- 5 years ago
- Views:
Transcription
1 How to reate Matlab Sript an Simulink Moel for Designing a Pole Plaement Controller
2 Pole Plaement Controller Design H tot Y( Kr B ( R ( AC ( ( BD ( ( P ( AC ( ( BD ( (
3 Polplaement Equation ( ( ( ( [egree] ( ( ( ( ( B P K D C n n n n n n n P D B C A r n n n n a b b a p Choose egree: Set:
4 Pole Plaement Example Bertil Thomas, Moern reglerteknik y( k Y,8Y H ( B(,8y U Y,5 k A(,8,5u,5U,5 k,5,5,8,5u,5u k { Z :} (x Desire: ubble pole in,5
5 Pole Plaement Example,5,5 ( ( n p P,8 (,5,5 ( A B Desire: ubble pole in,5 ( ( D n n C n n a b,5 (,5 (,5 (,8 ( : P BD AC
6 Pole Plaement Example,5,8 (,5,85,8 (,5,5,8,8 (,5 (,5 (,8 (,5,5,8,5,8,5,8,5,5,8 ACBD P
7 Pole Plaement Example,8,5,5,8,5,5,8,8,5,5,8,8,5,5,36,3,5» X[,.5; -.8,.5]\[-.8,.5] X »
8 Pole Plaement Example P B,5,5,5 K r P( B(,5,5,5,5 E U R K E C r Y D
9 Pole Plaement Example E U R K E C r Y D E,5 R,5Y U E,363 U E,363U Differens equations: U (,363 E e,5r, 5y k k k u k ek k,363u
10
11 Integrating Blok To avoi resiual error if we have a isturbane V we nee an integrating ontroller. A ( A( ( * During alulations the integrator is plae in A(, here alle A * (. ( easiest
12 With Integrating Blok Searhe P( has ubble pole,5 an extra pole in origo.,5 (,5 ( ( n p P (,8 ( (,8 (,5,5 ( * A A B ( ( D n n C n n a b *,5 (,5 (,5 ( (,8 ( : P BD A C
13 With Integrating Blok 3,5,5 (,8,5,5,8 (,8,5,8 ( (,5 (,5 ( (,8 (,5,8,5 (,8,57,5,5,8,5,5,5,8 (,8,8,5,5,8 (» X[,.5,;-.8,.5,.5;.8,,.5]\[.8,-.57,] X »
14 With Integrating Blok P B,5,5,5 K r P( B(,5,5,5,5 C * D( C( (,876,66 ( K r,5,6,38,88,63 (,38,68,5,5,8 Integrator now in C-blok! * C B A,38 D
15 Controller Differene Equations E U R K E * C r Y D E,5 R,88Y,63Y U,6 E,38 U,6U,38U E U E,6U,38U Differene Equations: ek,5rk,88 yk,63 yk u,6 k ek uk, 38uk
16
17 Matlab Sript for Lab 3 Pole plaement requires extensive alulations every time you want to try new poles. This is something you o not have time to o uring the lab you have to automate alulations using a Matlab sript. When solving the preparation tasks, evelop a MATLAB sript to be able to quikly reo the alulations uring the lab with the values of the real proess.
18 Matlab Sript for Lab 3 % poleplae.m % % Preparation for IE34 Lab 3: Pole plaement ontroller. % % E U GP % R -->Kr-->O--->/C(--->B(/A(----> Y % Htot(: -^ % % D(< % % Proess transfer funtion GP has two time onstants Start the sript with a omment (%. That omment is printe when writing help poleplae in the Matlab omman winow goo pratie to show what the sript oes.
19 Matlab Sript for Lab 3 Prior to the lab, we assume that proess gain an time onstants have the following values: % Proess transfer funtion GP [m] % has two time onstants Kp 3 % Change proess gain if neee T6 % Change first time onstant if neee T % Change seon time onstant if neee s tf('s'; GP Kp / ( T*s / ( T*s [GpNum,GpDen]tfata(GP,'v' These values must be hange to atual proess values at the lab [V]
20 Assume values prior to lab. KP GP { KP 3 T 6 T } ( Ts( Ts s s s s ( 6 ( Perform the system ientifiation at lab. Note! Use sript: [a,p]ap(q,n
21 Matlab Sript for Lab 3 It is ommon pratie to insert results from the omman winow as ommans in the sript! Kp 3 % Change proess gain if neee T 6 % Change first time onstant if neee T % Change seon time onstant if neee s tf('s'; GP Kp / ( T*s / ( T*s [GpNum,GpDen]tfata(GP,'v' % Default transfer funtion: % 3 % % 6 s^ 7 s % GpNum 3 % GpDen 6 7
22 Matlab Sript for Lab 3 For pole plaement alulations, we nee the proess isretetime moel. Matlab an print the moel in positive ( n or negative ( -n form. % Disretie the proess isp('now isretiing' h.6; % sampling interval HP(GP,h; fprintf('\nhp Negative notation' HP.variable '^-' fprintf('\nhp Positive notation' HP.variable '' % Disrete transfer funtion: % HP Negative notation % Transfer funtion: %.47 ^ ^- % % ^ ^- % Sampling time (seons:.6 % % HP Positive notation % Transfer funtion: % % % ^ % Sampling time (seons:.6
23 Matlab Sript for Lab 3 We nee proess transfer funtion in negative form, split into enominator (A an numerator (B. % HPB/A [B,A]tfata(HP,'v' % HP Negative notation % Transfer funtion: %.47 ^ ^- % % ^ ^- % B [.4.39] % A [ ]
24 Matlab Sript for Lab 3 A an B are both of seon egee an for now we use the following preliminary poles to reate the pole plaement polynomial:,7 ±,i an. % 3 poles suggestion: p.7.i; % Change if neee p.7-.i; % Change if neee p3; P ( AC ( ( BD ( ( % Prompt user for pole values % p input('first pole (abi: '; % p input('seon pole (a-bi: '; % isplay('thir pole at origo i'; % Poleplaement polynomial P poly([p,p,p3] % P [ ] When trying many ifferent pole values it might be onvenient to have the sript prompt the user!
25 Matlab Sript for Lab 3 Controller polynomials C an D: % Controller polynomials C an D % C _*^- % D *^- AC ( ( BD ( ( P ( Parameters oh _ are alulate by ientifiation of the ontroller polynomial s equivalent terms. This algebra an not be performe with Matlab (when the Maple moule is missing but may be one with paper an penil or with Mathematia. When the algebra is mae, Matlab an solve the equation system numerially. You must present the algebrai alulations at the lab!
26 Parameter names in Matlab % P [ ] [P(4 P(3 P( P(] P(4*^-3 P(4*^- P(4*^- P(4*^- % HPB/A [B,A]tfata(HP,'v' % B [.4.39] [B(3 B( B(] B(3*^- B(*^- B(*^- % A [ ] [A(3 A( A(] A(3*^- A(*^- A(*^- % Controller polynomials C an D % C [ _] _*^- % D [ ] *^- Ientify: AC ( ( BD ( ( P (
27 Matlab Sript for Lab 3 Desribe your algebrai alulations at the lab! This Matlab an solve left [ B( A( B(3 B( A(3 B(3 ]; right [P(-A(, P(3-A(3, P(4]' ; left \ right; _ ans(, _ ans(, _ ans(3 % _.48 % _ % _ Now we have the parameters _ alulate!
28 Matlab Sript for Lab 3 K r P( B( % KrP(/B( Kr polyval(p, / polyval(b, % Kr Now we have all the parameters!
29 Matlab Sript for Lab 3 % Controller bloks D [ ] D_blok tf( D,[ ],h; fprintf('\nblok D positive notation' D_blok.variable '' % D % % Blok D positive notation % Transfer funtion: % % % % Sampling time (seons:.6
30 Matlab Sript for Lab 3 % Controller bloks C [ _] C_blok tf([ ],C,h ; fprintf('\nblok C positive notation' C_blok.variable '' % C..48 % % Blok C positive notation % Transfer funtion: % % %.48 % Sampling time (seons:.6
31 Matlab Sript for Lab 3 % Total transfer funtion % Kr*B Htot % Htot % A*C B*D % Transfer funtion of ontrolle system: fprintf('\ntotal system, positive notation' Htot tf( Kr*B, onv(a,c onv(b,d,h ; Htot.variable '' Kr B ( AC ( ( BD ( ( % Transfer funtion: % % % ^3 -.4 ^ e-7
32 Matlab Sript for Lab 3 Commans to plot system info. % Print system info. figure('name','step Response of Proess','numbertitle','off' step(gp % step response of proess opt stepdataoptions('stepamplitue',.5, 'InputOffset',.5; figure('name','step Response of Feebak Loop','numbertitle','off' step(htot, opt % step response of feebak loop figure('name','margins of Feebak Loop','numbertitle','off' margin(htot % Gain an phase margins of feebak loop [Y,T]step(Htot, opt; stepinfo(y,t,'settlingtimethreshol',.5
33 Matlab Sript for Lab 3 How o you hek that the sript is orret an that the poles are orretly plae? [P,Z]pmap(Htot % P %.7.i %.7-.i % Z Poles are as esire, the sript is orret!
34
35 Simulink Moel
36 Export transfer funtions to Simulink Simulink an simulate ontinuous-time an isrete-time systems together. For the proess it is onvenient to use the ontinuous-time moel. [GpNum,GpDen]tfata(GP,'v' GpNum an GpDen polynomials are importe in Simulink moel. William Sanqvist
37 Export transfer funtions to Simulink % Controller bloks D [ ] D_blok tf( D,[ ],h; fprintf('\nblok D positive notation' D_blok.variable '' % D % % Blok D positive notation % Transfer funtion: % % % % Sampling time (seons:.6 D polynomial is importe in Simulink moel.
38 Export transfer funtions to Simulink C [ _] C_blok tf([ ],C,h ; fprintf('\nblok C positive notation' C_blok.variable '' % C..48 % % Blok C positive notation % Transfer funtion: % % %.48 % Sampling time (seons:.6 C polynomial is importe in Simulink moel.
39 Export transfer funtions to Simulink % KrP(/B( P( B( Kr polyval(p, / polyval(b, % Kr6.664 K r Kr value is importe in Simulink moel.
40
41
42 Plot Simulation Result in Matlab The following Matlab sript reates one figure for ontrol input an another for proess output. The plots are upate when a new simulation is mae. % stepplot.m % Plots output (y an ontrol input (u. figure('name','output, y','numbertitle','off' plot(t,y linkata on figure('name','control Input,u','numbertitle','off' stairs(t,u linkata on
43 Simulation Plots Simulation result when referene value is hange from.5v to 3V after 5s. Proess output y(k is ontinous, while ontrol input u(k is isrete.
44
45 Matlab sript for C-oe (Not part of the ourse. Matlab an write to a file! For example, the pole plaement ontroller an be printe as a C funtion, whih an be inlue when ompiling a C program. fifopen('pp.','w'; fprintf(fi,'float PolePlaeController(float r, float y\n'; fprintf(fi,'{\n'; fprintf(fi,' stati float pre_u ;\n'; fprintf(fi,' stati float pre_y ;\n'; fprintf(fi,' float e;\n'; fprintf(fi,' float u;\n\n'; fprintf(fi,' e %6.3f*r - %6.3f*y - %6.3f*pre_y;\n',Kr,_,_; fprintf(fi,' u e - %6.3f*pre_u;\n\n',_; fprintf(fi,' pre_u u;\n'; fprintf(fi, pre_y y;\n\n'; fprintf(fi,' return u;\n}\n'; flose(fi;
46 C-funtion (not part of the Course
How to Create Matlab Script and Simulink Model for Designing a Pole Placement Controller
How to Create Matlab Script and Simulink Model for Designing a Pole Placement Controller 1(37) Pole Placement Controller Design H tot Y ( z) Kr B ( z) R( z) A( z) C( z) B( z) D( z) P( z) A( z) C( z) B(
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationMcCreight s Suffix Tree Construction Algorithm. Milko Izamski B.Sc. Informatics Instructor: Barbara König
1. Introution MCreight s Suffix Tree Constrution Algorithm Milko Izamski B.S. Informatis Instrutor: Barbara König The main goal of MCreight s algorithm is to buil a suffix tree in linear time. This is
More informationECE-320 Linear Control Systems. Winter 2013, Exam 1. No calculators or computers allowed, you may leave your answers as fractions.
ECE-320 Linear Control Systems Winter 2013, Exam 1 No alulators or omputers allowed, you may leave your answers as frations. All problems are worth 3 points unless noted otherwise. Total /100 1 Problems
More informationSampler-B. Secondary Mathematics Assessment. Sampler 521-B
Sampler-B Seonary Mathematis Assessment Sampler 51-B Instrutions for Stuents Desription This sample test inlues 15 Selete Response an 5 Construte Response questions. Eah Selete Response has a value of
More informationUTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationEE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas
EE 370L Controls Laboratory Laboratory Exercise #7 Root Locus Department of Electrical an Computer Engineering University of Nevaa, at Las Vegas 1. Learning Objectives To emonstrate the concept of error
More informationEE451/551: Digital Control. Relationship Between s and z Planes. The Relationship Between s and z Planes 11/10/2011
/0/0 EE45/55: Digital Control Chater 6: Digital Control System Design he Relationshi Between s and Planes As noted reviously: s j e e e e r s j where r e and If an analog system has oles at: s n jn a jd
More informationThe simulation analysis of the bridge rectifier continuous operation in AC circuit
Computer Appliations in Eletrial Engineering Vol. 4 6 DOI 8/j.8-448.6. The simulation analysis of the bridge retifier ontinuous operation in AC iruit Mirosław Wiślik, Paweł Strząbała Kiele University of
More informationSome facts you should know that would be convenient when evaluating a limit:
Some fats you should know that would be onvenient when evaluating a it: When evaluating a it of fration of two funtions, f(x) x a g(x) If f and g are both ontinuous inside an open interval that ontains
More informationController Design Based on Transient Response Criteria. Chapter 12 1
Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors
More informationChapter 15 Equilibrium. Reversible Reactions & Equilibrium. Reversible Reactions & Equilibrium. Reversible Reactions & Equilibrium 2/3/2014
Amount of reatant/produt //01 quilibrium in Chemial Reations Lets look bak at our hypothetial reation from the kinetis hapter. A + B C Chapter 15 quilibrium [A] Why doesn t the onentration of A ever go
More information1. Description of Finite Heat Release Function
ME 40 Day 27 Desription of Finite Heat elease Funtion - SI Engines Differential Equations to Moel Cyle Software Implementation in EES Questions that an be answere. Desription of Finite Heat elease Funtion
More informationEE451/551: Digital Control. Chapter 7: State Space Representations
EE45/55: Digital Control Chapter 7: State Spae Representations State Variables Definition 7.: The system state is a minimal set of { variables x ( t ), i =, 2,, n needed together with the i input ut, and
More informationarxiv: v1 [physics.gen-ph] 5 Jan 2018
The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationForce Reconstruction for Nonlinear Structures in Time Domain
Fore Reonstrution for Nonlinear Strutures in ime Domain Jie Liu 1, Bing Li 2, Meng Li 3, an Huihui Miao 4 1,2,3,4 State Key Laboratory for Manufaturing Systems Engineering, Xi an Jiaotong niversity, Xi
More informationSpeed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion
Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3
More information1. A dependent variable is also known as a(n). a. explanatory variable b. control variable c. predictor variable d. response variable ANSWER:
1. A epenent variale is also known as a(n). a. explanatory variale. ontrol variale. preitor variale. response variale FEEDBACK: A epenent variale is known as a response variale. Definition of the Simple
More informationChapter 15 Equilibrium. Reversible Reactions & Equilibrium. Reversible Reactions & Equilibrium. Reversible Reactions & Equilibrium 5/27/2014
Amount of reatant/produt 5/7/01 quilibrium in Chemial Reations Lets look bak at our hypothetial reation from the kinetis hapter. A + B C Chapter 15 quilibrium [A] Why doesn t the onentration of A ever
More informationProblem set 6 for the course Theoretical Optics Sample Solutions
Karlsruher Institut für Tehnologie KIT) Institut für theoretishe Festkörperphysik SS01 Prof. Dr. G. Shön, Dr. R. Frank 15.06.01 http://www.tfp.kit.eu/stuium-lehre.php Tutorial: Group 1, Name: Group, Group
More information6.4 Dividing Polynomials: Long Division and Synthetic Division
6 CHAPTER 6 Rational Epressions 6. Whih of the following are equivalent to? y a., b. # y. y, y 6. Whih of the following are equivalent to 5? a a. 5, b. a 5, 5. # a a 6. In your own words, eplain one method
More informationMethods Higher Tier Practice Paper Unit 1 Markscheme
Methods Higher Tier Pratie Paper Unit Marksheme GCSE MATHEMATICS LINKED PAIR METHODS FOUNDATION NOTES ON MARKING PRINCIPLES Types of mark M marks: method marks A marks: auray marks B marks: unonditional
More informationAlgebra I. Slide 1 / 216. Slide 2 / 216. Slide 3 / 216. Polynomials
Slide 1 / 216 Slide 2 / 216 lgebra I Polynomials 2015-11-02 www.njctl.org Table of ontents efinitions of Monomials, Polynomials and egrees dding and Subtracting Polynomials Multiplying a Polynomial by
More informationmn 3 17x 2 81y 4 z Algebra I Definitions of Monomials, Polynomials and Degrees 32,457 Slide 1 / 216 Slide 2 / 216 Slide 3 / 216 Slide 4 / 216
Slide 1 / 216 Slide 2 / 216 lgebra I Polynomials 2015-11-02 www.njctl.org Slide 3 / 216 Table of ontents efinitions of Monomials, Polynomials and egrees dding and Subtracting Polynomials Multiplying a
More informationMCH T 111 Handout Triangle Review Page 1 of 3
Hnout Tringle Review Pge of 3 In the stuy of sttis, it is importnt tht you e le to solve lgeri equtions n tringle prolems using trigonometry. The following is review of trigonometry sis. Right Tringle:
More informationGeometry of Transformations of Random Variables
Geometry of Transformations of Random Variables Univariate distributions We are interested in the problem of finding the distribution of Y = h(x) when the transformation h is one-to-one so that there is
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationFundamental Theorem of Calculus
Chater 6 Fundamental Theorem of Calulus 6. Definition (Nie funtions.) I will say that a real valued funtion f defined on an interval [a, b] is a nie funtion on [a, b], if f is ontinuous on [a, b] and integrable
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationCommon Mistakes & How to avoid them Class X - Math. Unit: Algebra. Types of Question Common Mistakes Points to be emphasised. points.
Common Mistakes & How to avoid them Class X - Math Unit: Algera Chapter: Pair of Linear Equations in Two Variales Types of Question Common Mistakes Points to e emphasised Solving the system of (i) Error
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationMODELS FOR VARIABLE RECRUITMENT (continued)
ODL FOR VARIABL RCRUITNT (ontinue) The other moel ommonly use to relate reruitment strength with the size of the parental spawning population is a moel evelope by Beverton an Holt (957, etion 6), whih
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationSampler-A. Secondary Mathematics Assessment. Sampler 521-A
Sampler-A Seondary Mathematis Assessment Sampler 521-A Instrutions for Students Desription This sample test inludes 14 Seleted Response and 4 Construted Response questions. Eah Seleted Response has a
More informationYTÜ Mechanical Engineering Department
YTÜ Mechanical Engineering Department Lecture of Special Laboratory of Machine Theory, System Dynamics and Control Division Coupled Tank 1 Level Control with using Feedforward PI Controller Lab Report
More information6665/01 Edexcel GCE Core Mathematics C3 Silver Level S4
Paper Referene(s) 6665/0 Edexel GCE Core Mathematis C Silver Level S4 Time: hour 0 minutes Materials required for examination papers Mathematial Formulae (Green) Items inluded with question Nil Candidates
More informationCanimals. borrowed, with thanks, from Malaspina University College/Kwantlen University College
Canimals borrowed, with thanks, from Malaspina University College/Kwantlen University College http://ommons.wikimedia.org/wiki/file:ursus_maritimus_steve_amstrup.jpg Purpose Investigate the rate of heat
More informationFREQUENCY DOMAIN FEEDFORWARD COMPENSATION. F.J. Pérez Castelo and R. Ferreiro Garcia
FREQUENCY DOMAIN FEEDFORWARD COMPENSATION F.J. Pérez Castelo and R. Ferreiro Garia Dept. Ingeniería Industrial. Universidad de La Coruña javierp@ud.es, Phone: 98 7.Fax: -98-7 ferreiro@ud.es, Phone: 98
More information6 Dynamic Optimization in Continuous Time
6 Dynami Optimization in Continuous Time 6.1 Dynami programming in ontinuous time Consider the problem Z T max e rt u (k,, t) dt (1) (t) T s.t. k ú = f (k,, t) (2) k () = k, (3) with k (T )= k (ase 1),
More informationSELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc.
SELF-ERECTING, ROTARY MOTION INVERTED PENDULUM Quanser Consulting Inc. 1.0 SYSTEM DESCRIPTION The sel-erecting rotary motion inverte penulum consists of a rotary servo motor system (SRV-02) which rives
More informationControl Systems. Control Systems Design Lead-Lag Compensator.
Design Lead-Lag Compensator hibum@seoulteh.a.kr Outline Lead ompensator design in frequeny domain Lead ompensator design steps. Example on lead ompensator design. Frequeny Domain Design Frequeny response
More informationModelling and Simulation. Study Support. Zora Jančíková
VYSOKÁ ŠKOLA BÁŇSKÁ TECHNICKÁ UNIVERZITA OSTRAVA FAKULTA METALURGIE A MATERIÁLOVÉHO INŽENÝRSTVÍ Modelling and Simulation Study Support Zora Jančíková Ostrava 5 Title: Modelling and Simulation Code: 638-3/
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More informationAnnouncements. Office Hours Swap: OH schedule has been updated to reflect this.
SA Solving Announements Offie Hours Swap: Zavain has offie hours from 4-6PM toay in builing 460, room 040A. Rose has offie hours tonight from 7-9PM in Gates B26B. Keith has offie hours hursay from 2-4PM
More informationFinal Exam: know your section, bring your ID!
Chapter 15: Equilibrium Part 1 Read: BLB 15.1 3 HW: BLB 15:13,14, 21 Supplemental 15:1 4 Know: Chemial Equilibrium Catalysts Equilibrium Constant Equilibrium onstant expression Homogeneous/Heterogeneous
More informationCHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS. Professor Dae Ryook Yang
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 208 Dept. of Chemial and Biologial Engineering 0- Road Map of the Leture X Stability of losed-loop ontrol system
More informationImproved Extended Kalman Filter for Parameters Identification
Improved Extended Kalman Filter for Parameters Identifiation Peipei Zhao. Zhenyu Wang. ao Lu. Yulin Lu Institute of Disaster Prevention, Langfang,Hebei Provine, China SUMMARY: he extended Kalman filter
More informationCSE 5311 Notes 18: NP-Completeness
SE 53 Notes 8: NP-ompleteness (Last upate 7//3 8:3 PM) ELEMENTRY ONEPTS Satisfiability: ( p q) ( p q ) ( p q) ( p q ) Is there an assignment? (Deision Problem) Similar to ebugging a logi iruit - Is there
More informationChapter 2: One-dimensional Steady State Conduction
1 Chapter : One-imensional Steay State Conution.1 Eamples of One-imensional Conution Eample.1: Plate with Energy Generation an Variable Conutivity Sine k is variable it must remain insie the ifferentiation
More informationInstytut Fizyki Doświadczalnej Wydział Matematyki, Fizyki i Informatyki UNIWERSYTET GDAŃSKI
Instytut Fizyki Doświadzalnej Wydział Matematyki, Fizyki i Informatyki UNIWERSYTET GDAŃSKI I. Bakground theory. 1. Laser doppler anemometry LDA (Doppler model): a) LDA priniples; b) single and multi-hannel
More informationPN Code Tracking Loops
Wireless Information Transmission System Lab. PN Coe Traking Loops Institute of Communiations Engineering National Sun Yat-sen University Introution Coe synhronization is generally arrie out in two steps
More informationSimplify each expression. 1. 6t + 13t 19t 2. 5g + 34g 39g 3. 7k - 15k 8k 4. 2b b 11b n 2-7n 2 3n x 2 - x 2 7x 2
9-. Plan Objetives To desribe polynomials To add and subtrat polynomials Examples Degree of a Monomial Classifying Polynomials Adding Polynomials Subtrating Polynomials 9- What You ll Learn To desribe
More information18 Numerical Integration of Functions
Slightly moifie //9, /8/6 Firstly written at Marh 5 8 Numerial Integration of Funtions Introution Romberg Integration Gauss Quarature Aaptive Quarature Case Stuy: Root-Mean-Square Current DM869/Computational
More informationProcess engineers are often faced with the task of
Fluids and Solids Handling Eliminate Iteration from Flow Problems John D. Barry Middough, In. This artile introdues a novel approah to solving flow and pipe-sizing problems based on two new dimensionless
More informationMATLAB/Simulink Based Study of Different Approaches Using Mathematical Model of Differential Equations
I.J. Intelligent Systems and Appliations, 04, 05, -4 Published Online April 04 in MECS (http://www.mes-press.org/) DOI: 0.585/ijisa.04.05.0 MATLAB/Simulink Based Study of Different Approahes Using Mathematial
More informationComputer Engineering 4TL4: Digital Signal Processing (Fall 2003) Solutions to Final Exam
Computer Engineering TL: Digital Signal Proessing (Fall 3) Solutions to Final Exam The step response ynof a ausal, stable LTI system is: n [ ] = [ yn ] un, [ ] where un [ ] is the unit step funtion a Find
More informationState variable feedback
State variable feedbak We have previosly disssed systems desribed by the linear state-spae eqations Ax B y Cx n m with xt () R the internal state, t () R the ontrol inpt, and yt () R the measred otpt.
More informationChapter 13, Chemical Equilibrium
Chapter 13, Chemial Equilibrium You may have gotten the impression that when 2 reatants mix, the ensuing rxn goes to ompletion. In other words, reatants are onverted ompletely to produts. We will now learn
More informationAn Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances
An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory
More information2 How far? Equilibrium Answers
How far? Equilibrium Answers ratie: pages 37 39 1 Answer is D. Only a hange in temperature harges the value of the equilibrium onstant. Answer is D. [B] /[A] so [B] [A] or [B] [A] 1/ 3 Answer is B. Amounts
More informationExercise 3: Quadratic sequences
Exerise 3: s Problem 1: Determine whether eah of the following sequenes is: a linear sequene; a quadrati sequene; or neither.. 3. 4. 5. 6. 7. 8. 8;17;3;53;80; 3 p ;6 p ;9 p ;1 p ;15 p ; 1;,5;5;8,5;13;
More informationCommunicating Special Relativity Theory s Mathematical Inconsistencies
Communiating Speial Relatiity Theory s Mathematial Inonsistenies Steen B Bryant Primitie Logi, In, 704 Sansome Street, San Franiso, California 94111 Stee.Bryant@RelatiityChallenge.Com Einstein s Speial
More informationKAMILLA OLIVER AND HELMUT PRODINGER
THE CONTINUED RACTION EXPANSION O GAUSS HYPERGEOMETRIC UNCTION AND A NEW APPLICATION TO THE TANGENT UNCTION KAMILLA OLIVER AND HELMUT PRODINGER Abstrat Starting from a formula for tan(nx in the elebrated
More information18.06 Problem Set 4 Due Wednesday, Oct. 11, 2006 at 4:00 p.m. in 2-106
8. Problem Set Due Wenesy, Ot., t : p.m. in - Problem Mony / Consier the eight vetors 5, 5, 5,..., () List ll of the one-element, linerly epenent sets forme from these. (b) Wht re the two-element, linerly
More informationChapter 15: Chemical Equilibrium
Chapter 5: Chemial Equilibrium ahoot!. At eq, the rate of the forward reation is the rate of the reverse reation. equal to, slower than, faster than, the reverse of. Selet the statement that BEST desribes
More informationExperiment I Electric Force
Experiment I Electric Force Twenty-five hunre years ago, the Greek philosopher Thales foun that amber, the harene sap from a tree, attracte light objects when rubbe. Only twenty-four hunre years later,
More informationFalling Mass and Improved Euler
Lesson 3: Falling Mass and Improved Euler 3.1 Applied Problem. Consider a mass which is falling through a viscous medium. One example is a falling rock, and a second example is a person in parachute. These
More informationFinal Review. A Puzzle... Special Relativity. Direction of the Force. Moving at the Speed of Light
Final Review A Puzzle... Diretion of the Fore A point harge q is loated a fixed height h above an infinite horizontal onduting plane. Another point harge q is loated a height z (with z > h) above the plane.
More informationMasonry Beams. Ultimate Limit States: Flexure and Shear
Masonry Beams 4:30 PM 6:30 PM Bennett Banting Ultimate Limit States: Flexure and Shear Leture Outline 1. Overview (5) 2. Design for Flexure a) Tension Reinforement (40) b) Compression Reinforement (20)
More informationExperiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.01 Purpose of the Experiment: Experiment 03: Work and Energy In this experiment you allow a art to roll down an inlined ramp and run into
More informationTAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS
MISN-0-4 TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS f(x ± ) = f(x) ± f ' (x) + f '' (x) 2 ±... 1! 2! = 1.000 ± 0.100 + 0.005 ±... TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS by Peter Signell 1.
More informationIn this assignment you will build a simulation of the presynaptic terminal.
9.16 Problem Set #2 In this assignment you will buil a simulation of the presynapti terminal. The simulation an be broken own into three parts: simulation of the arriving ation potential (base on the Hogkin-Huxley
More informationFactorising FACTORISING.
Ftorising FACTORISING www.mthletis.om.u Ftorising FACTORISING Ftorising is the opposite of expning. It is the proess of putting expressions into rkets rther thn expning them out. In this setion you will
More informationPhysicsAndMathsTutor.com 1
PhysisAndMathsTutor.om. (a (i beam splitter [or semi-silvered mirror] (ii a ompensator [or a glass blok] allows for the thikness of the (semi-silvered mirror to obtain equal optial path lengths in the
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationVariable Impedance Control with an Artificial Muscle Manipulator Using Instantaneous Force and MR Brake
1 IEEE/RSJ International Conferene on Intelligent Robots an Systems (IROS) November -7, 1. Tokyo, Japan ariable Impeane Control with an Artifiial Musle Manipulator Using Instantaneous Fore an MR Brake
More informationExtended Spectral Nonlinear Conjugate Gradient methods for solving unconstrained problems
International Journal of All Researh Euation an Sientifi Methos IJARESM ISSN: 55-6 Volume Issue 5 May-0 Extene Spetral Nonlinear Conjuate Graient methos for solvin unonstraine problems Dr Basim A Hassan
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationOptimal torque control of permanent magnet synchronous machines using magnetic equivalent circuits
This oument ontains a post-print version of the paper Optimal torque ontrol of permanent magnet synhronous mahines using magneti equivalent iruits authore by W. Kemmetmüller, D. Faustner, an A. Kugi an
More informationTo work algebraically with exponential functions, we need to use the laws of exponents. You should
Prealulus: Exponential and Logisti Funtions Conepts: Exponential Funtions, the base e, logisti funtions, properties. Laws of Exponents memorize these laws. To work algebraially with exponential funtions,
More informationThe Lorenz Transform
The Lorenz Transform Flameno Chuk Keyser Part I The Einstein/Bergmann deriation of the Lorentz Transform I follow the deriation of the Lorentz Transform, following Peter S Bergmann in Introdution to the
More informationPlanar sheath and presheath
5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.
More informationTo derive the other Pythagorean Identities, divide the entire equation by + = θ = sin. sinθ cosθ tanθ = 1
Syllabus Objetives: 3.3 The student will simplify trigonometri expressions and prove trigonometri identities (fundamental identities). 3.4 The student will solve trigonometri equations with and without
More informationToday we re going to talk about time delays, affectionately referred to as e -st in the Laplace domain. Why e -st? Recall: time.
6.302 Feedba Systems Reitation 23: ime Delays oday we re going to tal about time delays, affetionately referred to as e -st in the Laplae domain. Why e -st? Reall: time h(t) delay h(t-) Laplae transform:
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationZero-Knowledge Protocols
he People Zero-Knowlege Protools 2 he wars Prover (Peggy) Claim I Verifier (Vi) S Seret Deision 2 {true, false} zero-knowlege protool allows Peggy to Convine Vi that her laim is true an that she knows
More informationPhysics 6C. Special Relativity. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB
Physis 6C Speial Relatiity Two Main Ideas The Postulates of Speial Relatiity Light traels at the same speed in all inertial referene frames. Laws of physis yield idential results in all inertial referene
More informationREVIEW QUESTIONS Chapter 15
hemistry 10 ANSWER EY REVIEW QUESTIONS hapter 15 1. A mixture of 0.10 mol of NO, 0.050 mol of H and 0.10 mol of HO is plaed in a 1.0-L flask and allowed to reah equilibrium as shown below: NO (g) + H (g)
More informationDepartment of Electrical and Computer Engineering ECED4601 Digital Control System Lab3 Digital State Space Model
Department of Electrical and Computer Engineering ECED46 Digital Control System Lab3 Digital State Space Model Objectives. To learn some MATLAB commands that deals with the discrete time systems.. To give
More informationLecture 14 - Using the MATLAB Control System Toolbox and Simulink Friday, February 8, 2013
Today s Objectives ENGR 105: Feedback Control Design Winter 2013 Lecture 14 - Using the MATLAB Control System Toolbox and Simulink Friday, February 8, 2013 1. introduce the MATLAB Control System Toolbox
More informationUNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS. V. N. Matveev and O. V. Matvejev
UNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS V. N. Matveev and O. V. Matvejev Joint-Stok Company Sinerta Savanoriu pr., 159, Vilnius, LT-315, Lithuania E-mail: matwad@mail.ru Abstrat
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationA population of 50 flies is expected to double every week, leading to a function of the x
4 Setion 4.3 Logarithmi Funtions population of 50 flies is epeted to doule every week, leading to a funtion of the form f ( ) 50(), where represents the numer of weeks that have passed. When will this
More information1. Which two values of temperature are equivalent to the nearest degree when measured on the Kelvin and on the
. Whih two values of teperature are equivalent to the nearest degree when easured on the Kelvin and on the Celsius sales of teperature? Kelvin sale Celsius sale A. 40 33 B. 273 00 C. 33 40 D. 373 0 2.
More informationInverse Functions. Review from Last Time: The Derivative of y = ln x. [ln. Last time we saw that
Inverse Functions Review from Last Time: The Derivative of y = ln Last time we saw that THEOREM 22.0.. The natural log function is ifferentiable an More generally, the chain rule version is ln ) =. ln
More informationA MATLAB Method of Lines Template for Evolution Equations
A MATLAB Metho of Lines Template for Evolution Equations H.S. Lee a, C.J. Matthews a, R.D. Braok a, G.C. Saner b an F. Ganola a a Faulty of Environmental Sienes, Griffith University, Nathan, QLD, 4111
More information