Control of industrial robots. Control of the interaction

Size: px
Start display at page:

Download "Control of industrial robots. Control of the interaction"

Transcription

1 Control o industrial robots Control o the interation Pro. Paolo Roo (paolo.roo@polimi.it) Politenio di Milano Dipartimento di Elettronia, Inormazione e Bioingegneria

2 Introdution So ar we have assumed that the robot is ontrolled in a ree environment However the robot has to interat with the environment either to manipulate it, or to avoid to ollide with it, or to interat with other devies or humans In order to eetively interat with the environment, the robot must be endowed with eternal (eteroeptive) sensors, whih make the robot ontroller aware o the situation Common sensors used to measure suh interation are ore sensors and vision sensors In the ollowing we will study the ontrol systems based on ore measurements. Control o industrial robots Control o the interation Paolo Roo

3 Interation with the environment A robot ommonly interats with the working environment, manipulating objets and perorming operations on suraes (think o grinding, deburring, part assembly). In modern robotis, situations are ommon where the robot physially interats with the operator (think o telemanipulation or rehabilitation). Control o industrial robots Control o the interation Paolo Roo

4 Interation with the environment We are now dealing with the modelling o the environment and the use o eternal sensors: eternal sensors environment task planner trajetory planner ontrol robot internal sensors Control o industrial robots Control o the interation Paolo Roo

5 Passive ontrol o ompliane A irst way to address the interation with the environment is to endow the manipulator with devies that ailitate the eeution o the task in a passive way: RCC Remote Center o Compliane The RCC is used in assembly tasks (peg-in-a-hole) It is typially plaed between the robot s wrist and the gripper The RCC lets the gripper assembly move in the plane perpendiular to the peg s ais. This allows the peg to rotate With the RCC, the ores generated by any misalignment an be ompensated in a passive way Control o industrial robots Control o the interation Paolo Roo

6 Ative ontrol o the interation During more general interation tasks, the use o a purely positional ontrol strategy (the same adopted in ree motion) may lead to problems due to positioning errors and unertainties in task planning, related to a inomplete knowledge o the environment. In many appliations ontrolling the ontat ore is as important as ontrolling the speed o the robot alongthe presribed diretiont. Control o industrial robots Control o the interation Paolo Roo

7 Ative ontrol o the interation Eamples where an ative ontrol o the interation is desirable inlude: mehanial mahining (deburring, surae inishing, polishing, assembly,...) tele-manipulation ooperation o multi-manipulator systems deterous robot hands physial human-robot interation For an ative ontrol o interation it is neessary to use devies or ore (and moment) measurement. Control o industrial robots Control o the interation Paolo Roo

8 Fore sensors Fore/torque sensors are devies that return the measurements o three omponents o ore and three omponents o moment with respet to a loal rame. Control o industrial robots Control o the interation Paolo Roo

9 Internal behaviour Sensors are based on strain gauges, i.e. devies that an measure a strain. The strain indues a variation in the resistane o a iruit. Strain gauges are suitably mounted, in suh a way as to return inormation rom whih three omponents o ore and moment an be omputed. Control o industrial robots Control o the interation Paolo Roo

10 Statis o the manipulator When studying the interation o the robot with the environment, the problem arises how to inlude the eet o a ore applied at the end-eetor in the model. Suppose then that a ore is applied at the endeetor. What are the joint torques that keep the system in equilibrium???? Control o industrial robots Control o the interation Paolo Roo

11 Systems o ores As a matter o at a system o ores may at at the end-eetor. How an this system be redued? From lassial mehanis One ore applied in a point (the resultant) One torque It ollows that we will always be able to represent a system o ores with two vetors: A ore vetor applied in a point O (the resultant) A moment vetor (the moment o the original system o ores with respet to the ulrum O) Control o industrial robots Control o the interation Paolo Roo

12 Statis o the manipulator Let us then onsider a manipulator subjet to a set o ores at the end eetor, having a resultant and moment, w.r.t. the origin o the end eetor rame, µ. We then deine vetor: h = µ The virtual work o the ores/torques at the joints is: T δw τ = τ δq T T δw = δp + µ ωdt = J µ h Let us all τ the vetor o torques (or ores) applied to the joints. where δq is the elementary displaement o the joint oordinates. The virtual work o the ores at the end-eetor is: T T ( q) δq + J ( q) δq = h J( q) δq Control o industrial robots Control o the interation Paolo Roo T P where δp is the virtual displaement o the origin o the end-eetor rame, onsistent with the dierential kinematis o the robot, and ω is the angular veloity o the end eetor rame. O

13 Statis o the manipulator T δw τ = τ δq T δw = J µ h P T T ( q) δq + J ( q) δq = h J( q) δq, O J ( q) = J J P O ( q) ( ) q position Jaobian orientation Jaobian Based on the priniple o virtual work the manipulator is in stati equilibrium i and only i: δw τ = δw h τ = J T ( q)h Notie that the same matri, the Jaobian matri, whih deines the dierential kinematis o the manipulator, is also involved in the statis o the system (kinetostati duality). Control o industrial robots Control o the interation Paolo Roo

14 Taonomy Fore based interation ontrol strategies an be divided into two ategories: Impedane/admittane ontrol The goal is to assign a presribed dynami relation between interation ores and position errors We want that the manipulator, in ase o unplanned interations, behaves like a generalized mass-spring-damper system Hybrid position and ore ontrol We separate diretions whih are only position ontrolled rom diretions whih are only ore/torque ontrolled In the diretions onstrained by the environment we want that a ore/torque with a speiied value be established Furthermore we will talk about: Epliit ontrol The ore (or the impedane) is diretly assigned, ating on the ontrol variables (joint torques) Impliit ontrol The ore (or the impedane) is indiretly assigned, ating on the set-points o the position ontrol loops Control o industrial robots Control o the interation Paolo Roo

15 Mehanial impedane The low o power among two physial systems an always be deined as the produt o two onjugate quantities: a (generalized) eort and a (generalized) low In the eletrial ramework, the low is the urrent, while the eort is the voltage. In the mehanial ramework, instead, the low is the (linear or angular) veloity, the eort is the ore (or the moment) The mehanial impedane is then deined as dynamial relation that is established between ore and veloity (or displaement) or a mehanial system. The admittane is just the reiproal o the impedane The impedane ontrol aims at making the manipulator, position ontrolled and in interation with the environment, assume a desired mehanial impedane, like a generalized mass-spring-damper system Control o industrial robots Control o the interation Paolo Roo

16 Mehanial impedane Eamples o impedane/admittane ontrol or physial human-robot interation: Control o industrial robots Control o the interation Paolo Roo

17 A 1 d.o.. system In order to understand how an impedane ontroller operates, let us irst onsider the simple eample o a body with mass M interating with the environment, whih eerts a ore : p u M The model o the system is simply: Ma = u + Notie that the model o the interation with the environment is not speiied: we just know that there is an eternal ore ating on the system Control o industrial robots Control o the interation Paolo Roo

18 A 1 d.o.. system I we assume a ontrol proportional to position and veloity: u = k1p k2v We obtain: Ma k v + k p = k 1 k 2 p M The ontrolled system reats to the eternal ore like a mass-spring-damper The spring and damper ators are tunable, the mass is the physial one o the system. Control o industrial robots Control o the interation Paolo Roo

19 A 1 d.o.. system I we want to hange the mass, too, we need to measure the ore: u M M ( k p k v + ) = 1 2 d We obtain: M d k 1 M d a k v + k p = k 2 p All the parameters o the mass, spring and damper system are now tunable. Control o industrial robots Control o the interation Paolo Roo

20 Etension to the whole manipulator We now want to etend the impedane ontrol to the whole manipulator. The dynami model in interation with the environment is given by the equation: ( ) q T + C( q, q ) q + g( q) = τ J ( q)h B q where h is the vetor o the ontat ores (and moments) applied by the end eetor on the environment. Let us onsider now the adoption o an inverse dynamis ontrol law: τ ( q) y + C( q, q ) q g( q) = B + Substituting we obtain: nonlinear oupling term q = y B 1 T ( q) J ( q)h Control o industrial robots Control o the interation Paolo Roo

21 Inverse dynamis ontrol Assume now or y the epression introdued in the inverse dynami ontrol in the operational spae: y By substituting the epression or y we obtain the equation o the dynamis in losed loop: ~ where = d, while: B A 1 = J A ( ) q q ( q) + K ~ + K ~ J ( q, ) d T 1 ( q) J ( q) B( q) J ( q) D ~ K P ~ 1 K D ~ + + = B A = A A P ( q) h A is the inertia matri in the operational spae and: A generalized ores perorming work on ẋ h A = T T A ( ) h, J = T A ( ) J A T J h = J T A h A Control o industrial robots Control o the interation Paolo Roo

22 Impedane ontrol with ore measurement ~ K P ~ 1 K D ~ + + = B A ( q) h A We have obtained an impedane relation whih is oupled and only partially assignable: a ore applied along a Cartesian diretion yields motion in other diretions as well. I ore measurements are available, we an set: τ T ( ) y + C( q, q ) q + g( q) + J ( q)h = B q y = J 1 A ( ) q q h 1 ( q) M M + D ~ + K ~ M J ( q, ) d d d d d d A A where M d, D d, K d are diagonal, positive deinite matries. We obtain: M ~ + D ~ + K ~ = d d d h A i.e. a ompletely deoupled system. Control o industrial robots Control o the interation Paolo Roo

23 Impedane ontrol The equation: M ~ + D ~ + K ~ = d d d h A deines a mehanial impedane relation between the ores/moments deined at the right hand side o the equations and the position/orientation errors in the operational spae. This impedane has the orm o a mehanial system with mass matri M d, damping D d, and stiness K d. The dynami behaviour along eah diretion o the operational spae an thus be speiied: a ore or a moment along a ertain Cartesian diretion yields a motion in suh diretion, whih is related to the ore and moment by a ully ontrollable dynamis. But what do we eatly mean or impedane in the operational spae? Control o industrial robots Control o the interation Paolo Roo

24 Translational impedane The irst three degrees o reedom in the operational spae are the position oordinate o the end-eetor. In this ase the deinition o the impedane is rather intuitive: p d : desired position p e : atual position p = p d p e The impedane relation is: M p + D p + K p = p p p We thus apply to eah degree o reedom what we have already seen or the single mass. Control o industrial robots Control o the interation Paolo Roo

25 Rotational impedane The remaining three degrees o reedom in the operational spae onsist in a minimal representation (Euler angles) o the end-eetor orientation. In this ase the deinition o the impedane is less intuitive: φ d : desired orientation φ e : atual orientation φ = φ d φ e The impedane relation is: M φ + D φ + K φ = φ φ φ Τ T ( φ )µ e Due to the presene o matri T, this relation is subjet to representation singularities. Better deinition o the rotational impedane have been proposed in the literature. Control o industrial robots Control o the interation Paolo Roo

26 Impliit impedane ontrol (admittane ontrol) When a positional ontrol is already applied to the mehanial system, it might be onvenient to lose the impedane ontrol outside the position ontrol. Consider again the ase o single d.o.. system, position ontrolled: d R u M The system an be modeled with the blok diagram: u d R(s) G(s) 1 ( s) 2 G = Ms Control o industrial robots Control o the interation Paolo Roo

27 Impliit impedane ontrol (admittane ontrol) Consider the ollowing blok diagram, where the ore applied by the environment is measured and iltered through a system with transer untion G d (s), to give a variation o the reerene d : d + + R(s) u + + G(s) G d (s) impedane ontroller Elaborating the blok diagram: = G s 1+ R s ( ) ( ) G( s) + ( G ( s) ) d d ( ) G( s) ( ) G( s) R s 1+ R s Also alled admittane ontroller (it aepts a ore and delivers a displaement) Control o industrial robots Control o the interation Paolo Roo

28 Impliit impedane ontrol (admittane ontrol) In a band o requenies well inside the bandwidth o the position ontrol system, we an assume the ore disturbane well rejeted and the position reerene orretly traked. Thus: d Gd ( s) = G ( s) 1~ d I we impose: G d ( s) = M d s D d s + K d M ~ d + Dd ~ + Kd ~ = we assign to the manipulator the desired impedane behaviour. The advantage o this sheme is that the position ontrol is kept unhanged: the impedane ontrol modiies the set-points o the position ontrol, whih is usually the only openness allowed by an industrial motion ontrol system. Control o industrial robots Control o the interation Paolo Roo

29 Admittane ontrol or the whole manipulator The same sheme as with the single mass an be used, applied to eah Cartesian diretion o the end-eetor: p p d p d + R u, y, z The approah an be generalized to aount also or moments ating at the end-eetor and onsequent variations to the desired orientation o the end-eetor F Admittane ilter Control o industrial robots Control o the interation Paolo Roo

30 Admittane ontrol used or manual guidane Control o industrial robots Control o the interation Paolo Roo

31 Natural onstraints and artiiial onstraints An alternative way to perorm an interation task is to assign reerene values o the ores and o the positions, onsistent with the geometry o the environment. Some diretions an be in at subjeted to onstraints on the position, other ones to onstraints on the ores that an be applied. In general, it is not possible to impose both the ore, and the position, along the same diretion. We make the ollowing distintion: Natural onstraints: Artiiial onstraints: they are imposed by the environment along eah degree o reedom o the task; they depend on the geometry o the task. onstraints imposed by the ontrol system, related to the task eeution strategy. For task desription we reer to a onstraint rame, y, z. Control o industrial robots Control o the interation Paolo Roo

32 Sliding on a lat surae Natural Constraints p ω ω z y y µ z Artiiial Constraints µ z µ p p ω y y z This piture is taken rom the tetbook: B. Siiliano, L. Siavio, L. Villani, G. Oriolo: Robotis: Modelling, Planning and Control, 3rd Ed. Springer, 2009 Control o industrial robots Control o the interation Paolo Roo

33 Insertion o a peg in a hole This piture is taken rom the tetbook: B. Siiliano, L. Siavio, L. Villani, G. Oriolo: Robotis: Modelling, Planning and Control, 3rd Ed. Springer, 2009 Natural Constraints p p ω ω y y z µ z Artiiial Constraints. Usually we set p z 0 and all the other artiiial onstraints to zero. µ µ y p ω y z z Control o industrial robots Control o the interation Paolo Roo

34 Rotation o a rank Oriented like the rank Natural Constraints p p ω z Artiiial Constraints z µ Idle knob ω y y µ z µ p ω y y z This piture is taken rom the tetbook: B. Siiliano, L. Siavio, L. Villani, G. Oriolo: Robotis: Modelling, Planning and Control, 3rd Ed. Springer, 2009 Normally artiiial onstraints on ores and moments are assumed to be zero. Control o industrial robots Control o the interation Paolo Roo

35 Seletion matri We assume that the onstraints an always be epressed along mutually orthogonal diretions o rame C. Thus on eah diretion either a natural onstraint on position, or a natural onstraint on ore, eists. In order to speiy the onstraints we use a seletion matri Σ diagonal and with null or unitary elements: Natural onstraints Artiiial onstraints n ( I ) h h n Σv = v Σ = ( I Σ) v = v Σh = h a Σ = a The seletion matri etrats diretions where a natural onstraint in veloity applies and thus an artiiial onstraint in ore an be enored Control o industrial robots Control o the interation Paolo Roo

36 Hybrid ore/position ontrol Hybrid ore/position ontrol is a ontrol strategy that uses the artiiial onstraints to speiy the goals o the ontrol system, in terms o both ore and position. The 6 diretions o the operational spae are divided into position/veloity ontrolled diretions and ore/moment ontrolled diretions: Σ h Ad Σ INTERACTION CONTROL y F + y NONLINEAR u MANIPULATOR h A v d I Σ MOTION CONTROL y v + COMPENSATIONS AND ENVIRONMENT v All quantities reerred to the rame C I Σ Control o industrial robots Control o the interation Paolo Roo

37 Soures o inonsisteny Hybrid position/ore ontrol is based on a nominal model o the interation. Inonsisteny may however our in the measurements, due e.g. to: 1. rition at the ontat (a ore is deteted in a nominally ree diretion) 2. ompliane in the robot struture and/or at the ontat (a displaement is deteted in a diretion whih is nominally onstrained in motion) 3. unertainty in the environment geometry at the ontat The irst two soures o inonsisteny are automatially iltered out through the seletion matries. The third soure o inonsisteny an be mitigated by real time estimation proess. Control o industrial robots Control o the interation Paolo Roo

38 A 1 d.o.. system In order to have a better idea on how the ontrol system operates along the ore ontrolled diretions, let s onsider the simple eample o a body o mass M in interation with the environment: u M We assume a purely elasti model o the environment. Blok diagram: u + G(s) K u G (s) G K ( s) = 2 Ms + K Control o industrial robots Control o the interation Paolo Roo

39 A 1 d.o.. system u G (s) G 2 K n ( s) = ω = ωn = Ms + K s + ω, n K M The system presents a unitary gain and resonant, high requeny (or a sti ontat) dynamis. Notie that = u or K. The best way to ontrol this system is through an integral ontroller: d + u R R (s) G (s) ( s) = K s i ω K i This way the bandwidth o the ore ontrol is diretly assigned with the integrator, leaving outside the bandwidth the high requeny dynamis o the interation with the environment. Control o industrial robots Control o the interation Paolo Roo

40 Impliit ore ontrol As with impedane ontrol, when a position ontrol is already ative, it might be onvenient to lose the ore ontrol around the position ontrol loop. Consider or simpliity only the ase o a, position ontrolled, single degree o reedom: d R u M The system an be modelled through the ollowing blok diagram: d + R(s) u + G(s) K The advantage is that the position ontrol is let unhanged. Control o industrial robots Control o the interation Paolo Roo

41 Impliit ore ontrol Solving the blok diagram: ( s) ( s) = d 1+ G s + R( s) G( s) K ( ) K R( s) G( s) Under the assumption o an ininitely sti ontat (K ) we obtain =0 and u: d R(s) u 1 The position ontrol loop gets open! d I a PID position ontroller is used: ( ) R s we an adopt or the ore ontrol the sheme: + R (s) C(s) d 2 KD s + KPs + K = s u R(s) 1 I ( ) C s = K Partially ompensates the eet o the PID ontroller D s K P s + K I Control o industrial robots Control o the interation Paolo Roo

42 Impliit ore ontrol Then, ater ompensating the eet o the PID ontroller: d + R (s) 1 s We an take R as a PI ontroller: ( s) R = Kp + K s i Tuning o this ontroller is straightorward: ω K p L sk p ( s) = 2 s + K i K i K p Control o industrial robots Control o the interation Paolo Roo

43 Industrial appliation o ore ontrol ABB RobotWare Mahining Fore Control (mahining o blades): Control o industrial robots Control o the interation Paolo Roo

Control of industrial robots

Control of industrial robots Control o industrial robots Control with eternal sensors Pro. Paolo Roo (paolo.roo@polimi.it) Politenio di Milano Dipartimento di Elettronia, Inormazione e Bioingegneria Introdution So ar we have assumed

More information

f 2 f n where m is the total mass of the object. Expression (6a) is plotted in Figure 8 for several values of damping ( ).

f 2 f n where m is the total mass of the object. Expression (6a) is plotted in Figure 8 for several values of damping ( ). F o F o / k A = = 6 k 1 + 1 + n r n n n RESONANCE It is seen in Figure 7 that displaement and stress levels tend to build up greatly when the oring requeny oinides with the natural requeny, the buildup

More information

Control Theory association of mathematics and engineering

Control 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 information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Finite Formulation of Electromagnetic Field

Finite Formulation of Electromagnetic Field Finite Formulation o Eletromagneti Field Enzo TONTI Dept.Civil Engin., Univ. o Trieste, Piazzale Europa 1, 34127 Trieste, Italia. e-mail: tonti@univ.trieste.it Otober 16, 2000 Abstrat This paper shows

More information

Force/Vision Based Active Damping Control of Contact Transition in Dynamic Environments

Force/Vision Based Active Damping Control of Contact Transition in Dynamic Environments Fore/Vision Based Ative Damping Control of Contat Transition in Dynami Environments Tomas Olsson, Rolf Johansson, Anders Robertsson Department of Automati Control, Lund Inst. of Teh., Lund University,

More information

Part G-4: Sample Exams

Part G-4: Sample Exams Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

Contact State Estimation using Multiple Model Estimation and Hidden Markov Models

Contact State Estimation using Multiple Model Estimation and Hidden Markov Models Thomas J. Debus Pierre E. Dupont Department of Aerospae and Mehanial Engineering Boston University Boston, MA 02445, USA Robert D. Howe Division of Engineering and Applied Sienes Harvard University Cambridge,

More information

Motor Sizing Application Note

Motor Sizing Application Note PAE-TILOGY Linear Motors 70 Mill orest d. Webster, TX 77598 (8) 6-7750 ax (8) 6-7760 www.trilogysystems.om E-mail emn_support_trilogy@parker.om Motor Sizing Appliation Note By Jak Marsh Introdution Linear

More information

Relativistic Dynamics

Relativistic Dynamics Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable

More information

Department of Mechanical Engineering

Department of Mechanical Engineering Department o Mehanial Engineering AMEE41 / ATO4 Aerodynamis Instrutor: Marios M. Fyrillas Email: eng.m@it.a.y Homework Assignment #4 QESTION 1 Consider the boundary layer low on a lat plate o width b (shown

More information

Practice Exam 2 Solutions

Practice Exam 2 Solutions MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department o Physis Physis 801T Fall Term 004 Problem 1: stati equilibrium Pratie Exam Solutions You are able to hold out your arm in an outstrethed horizontal position

More information

SEISMIC ANALYSIS OF SPHERICAL TANKS INCLUDING FLUID-STRUCTURE-SOIL INTERACTION

SEISMIC ANALYSIS OF SPHERICAL TANKS INCLUDING FLUID-STRUCTURE-SOIL INTERACTION 3 th World Conferene on Earthquake Engineering Vanouver, B.C., Canada August -6, 2004 aper o. 84 SEISMIC AALYSIS OF SHERICAL TAKS ICLUDIG FLUID-STRUCTURE-SOIL ITERACTIO T.L. Karavasilis, D.C. Rizos 2,

More information

Development of moving sound source localization system

Development of moving sound source localization system Eletroni Journal «Tehnial Aoustis» http://www.ejta.org 2006, 11 Doh-Hyoung Kim 1, Youngjin Park 2 Korea Advaned Institute o Siene and Tehnology, ME3076, 373-1, Guseong-Dong, Yuseong-Gu, Daejon, 305-701,

More information

v = fy c u = fx c z c The Pinhole Camera Model Camera Projection Models

v = fy c u = fx c z c The Pinhole Camera Model Camera Projection Models The Pinhole Camera Model Camera Projetion Models We will introdue dierent amera projetion models that relate the loation o an image point to the oordinates o the orresponding 3D points. The projetion models

More information

Advanced Computational Fluid Dynamics AA215A Lecture 4

Advanced Computational Fluid Dynamics AA215A Lecture 4 Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas

More information

Speed Regulation of a Small BLDC Motor using Genetic-Based Proportional Control

Speed Regulation of a Small BLDC Motor using Genetic-Based Proportional Control World Aademy of Siene, Engineering and Tehnology 47 8 Speed Regulation of a Small BLDC Motor using Geneti-Based Proportional Control S. Poonsawat, and T. Kulworawanihpong Abstrat This paper presents the

More information

Relativity in Classical Physics

Relativity in Classical Physics Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of

More information

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept

More information

Comparison of Alternative Equivalent Circuits of Induction Motor with Real Machine Data

Comparison of Alternative Equivalent Circuits of Induction Motor with Real Machine Data Comparison of Alternative Equivalent Ciruits of Indution Motor with Real Mahine Data J. radna, J. auer, S. Fligl and V. Hlinovsky Abstrat The algorithms based on separated ontrol of the motor flux and

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion

Speed-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 information

Measuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach

Measuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La

More information

A model for measurement of the states in a coupled-dot qubit

A model for measurement of the states in a coupled-dot qubit A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:

More information

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,

More information

The Unified Geometrical Theory of Fields and Particles

The Unified Geometrical Theory of Fields and Particles Applied Mathematis, 014, 5, 347-351 Published Online February 014 (http://www.sirp.org/journal/am) http://dx.doi.org/10.436/am.014.53036 The Unified Geometrial Theory of Fields and Partiles Amagh Nduka

More information

Kinematics and inverse dynamics analysis for a general 3-PRS spatial parallel mechanism Yangmin Li and Qingsong Xu

Kinematics and inverse dynamics analysis for a general 3-PRS spatial parallel mechanism Yangmin Li and Qingsong Xu Robotia (2005) volume 23, pp. 219 229. 2005 Cambridge University Press doi:10.1017/s0263574704000797 Printed in the United Kingdom Kinematis and inverse dynamis analysis for a general 3-PRS spatial parallel

More information

Doppler effect of the rupture process of the great M W 7.9 Wenchuan earthquake

Doppler effect of the rupture process of the great M W 7.9 Wenchuan earthquake Earthq Si (1)3: 535 539 535 539 Doi: 1.17/s11589-1-75-4 Doppler eet o the rupture proess o the great M W 7.9 Wenhuan earthquake Ge Jin 1,, Youai Tang 1 Shiyong Zhou 1 and Yongshun John Chen 1 1 Institute

More information

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J.

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J. A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE W. J. Haering* Senior Projet Engineer General Motors Corporation Warren, Mihigan R.

More information

Effect of Droplet Distortion on the Drag Coefficient in Accelerated Flows

Effect of Droplet Distortion on the Drag Coefficient in Accelerated Flows ILASS Amerias, 19 th Annual Conerene on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Eet o Droplet Distortion on the Drag Coeiient in Aelerated Flows Shaoping Quan, S. Gopalakrishnan

More information

Maximum Entropy and Exponential Families

Maximum 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 information

PREDICTION OF CONCRETE COMPRESSIVE STRENGTH

PREDICTION OF CONCRETE COMPRESSIVE STRENGTH PREDICTION OF CONCRETE COMPRESSIVE STRENGTH Dunja Mikuli (1), Ivan Gabrijel (1) and Bojan Milovanovi (1) (1) Faulty o Civil Engineering, University o Zagreb, Croatia Abstrat A ompressive strength o onrete

More information

Acoustic Attenuation Performance of Helicoidal Resonator Due to Distance Change from Different Cross-sectional Elements of Cylindrical Ducts

Acoustic Attenuation Performance of Helicoidal Resonator Due to Distance Change from Different Cross-sectional Elements of Cylindrical Ducts Exerpt rom the Proeedings o the COMSOL Conerene 1 Paris Aousti Attenuation Perormane o Helioidal Resonator Due to Distane Change rom Dierent Cross-setional Elements o Cylindrial Duts Wojieh ŁAPKA* Division

More information

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry 9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived

More information

Conformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices

Conformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices AAS 03-190 Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX 77843-3141 Abstrat This

More information

A novel Infrared Thermography (IRT) based experimental technique for distributed temperature measurements in hot gas flows

A novel Infrared Thermography (IRT) based experimental technique for distributed temperature measurements in hot gas flows QIRT 1 1 http://dx.doi.org/1.1611/qirt.1.45 th International Conerene on Quantitative InraRed Thermography July 7-3, 1, Québe (Canada A novel Inrared Thermography (IRT based experimental tehnique or distributed

More information

arxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006

arxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006 The simplest derivation of the Lorentz transformation J.-M. Lévy Laboratoire de Physique Nuléaire et de Hautes Energies, CNRS - IN2P3 - Universités Paris VI et Paris VII, Paris. Email: jmlevy@in2p3.fr

More information

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are

More information

SURFACE WAVES OF NON-RAYLEIGH TYPE

SURFACE WAVES OF NON-RAYLEIGH TYPE SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh

More information

Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems

Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems Rudranarayan M. Mukherjee, Kurt S. Anderson Computational Dynamis Laboratory Department of Mehanial Aerospae

More information

Normative and descriptive approaches to multiattribute decision making

Normative and descriptive approaches to multiattribute decision making De. 009, Volume 8, No. (Serial No.78) China-USA Business Review, ISSN 57-54, USA Normative and desriptive approahes to multiattribute deision making Milan Terek (Department of Statistis, University of

More information

Aharonov-Bohm effect. Dan Solomon.

Aharonov-Bohm effect. Dan Solomon. Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that

More information

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo

More information

Finite Element Modelling of Vibro-Acoustic Systems for Active Noise Reduction

Finite Element Modelling of Vibro-Acoustic Systems for Active Noise Reduction ECHNISCHE MECHANIK, Band 25, Heft 3-4, (25), 241-247 Manuskripteingang: 11. Dezember 25 Finite Element Modelling of Vibro-Aousti Systems for Ative Noise Redution J. Lefèvre, U. Gabbert Over the past years

More information

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential

More information

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia

More information

INFORMATION CONCERNING MATERIALS TO BE USED IN THE DESIGN

INFORMATION CONCERNING MATERIALS TO BE USED IN THE DESIGN TITLE 5 DESIGN CHAPTER 8 INFORMATION CONCERNING MATERIALS TO BE USED IN THE DESIGN Artile 38. Charateristis o steel or reinorements 38.1 General The harateristis o the steel used or the design desribed

More information

13 Fitzroy Street London W1T 4BQ Telephone: +44 (0) Facsimile: +44 (0)

13 Fitzroy Street London W1T 4BQ Telephone: +44 (0) Facsimile: +44 (0) Oasys AdSe Theory 13 Fitzroy Street London W1T 4BQ Telephone: +44 (0) 0 7755 330 Fasimile: +44 (0) 0 7755 370 Central Square Forth Street Newastle Upon Tyne NE1 3PL Telephone: +44 (0) 191 38 7559 Fasimile:

More information

Finite-time stabilization of chaotic gyros based on a homogeneous supertwisting-like algorithm

Finite-time stabilization of chaotic gyros based on a homogeneous supertwisting-like algorithm OP Conferene Series: Materials Siene Engineering PAPER OPEN ACCESS Finite-time stabilization of haoti gyros based on a homogeneous supertwisting-like algorithm To ite this artile: Pitha Khamsuwan et al

More information

The simulation analysis of the bridge rectifier continuous operation in AC circuit

The 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 information

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite. Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need

More information

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

23.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 information

A note on a variational formulation of electrodynamics

A note on a variational formulation of electrodynamics Proeedings of the XV International Workshop on Geometry and Physis Puerto de la Cruz, Tenerife, Canary Islands, Spain September 11 16, 006 Publ. de la RSME, Vol. 11 (007), 314 31 A note on a variational

More information

Directional Coupler. 4-port Network

Directional Coupler. 4-port Network Diretional Coupler 4-port Network 3 4 A diretional oupler is a 4-port network exhibiting: All ports mathed on the referene load (i.e. S =S =S 33 =S 44 =0) Two pair of ports unoupled (i.e. the orresponding

More information

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances

An 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 information

Wave Propagation through Random Media

Wave Propagation through Random Media Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene

More information

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR Bulletin of the Transilvania University of Braşov Series II: Forestry Wood Industry Agriultural Food Engineering Vol. 9 (58) No. - 06 IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED

More information

Dynamic Programming and Multi Objective Linear Programming approaches

Dynamic Programming and Multi Objective Linear Programming approaches 0 NSP Vol. () (0), 9 Dynami Programming and Multi Objetive Linear Programming approahes P. K. De and Amita Bhinher Department o Mathematis, National Institute o Tehnology SILCHAR - 788 00 (Assam) India

More information

Grasp Planning: How to Choose a Suitable Task Wrench Space

Grasp Planning: How to Choose a Suitable Task Wrench Space Grasp Planning: How to Choose a Suitable Task Wrenh Spae Ch. Borst, M. Fisher and G. Hirzinger German Aerospae Center - DLR Institute for Robotis and Mehatronis 8223 Wessling, Germany Email: [Christoph.Borst,

More information

PIEZOELECTRIC materials are widely used for structural

PIEZOELECTRIC materials are widely used for structural Parametri Analysis o the Vibration Control o Sandwih Beams Through Shear-Based Piezoeletri Atuation M. A. TRINDADE,* A. BENJEDDOU AND R. OHAYON Strutural Mehanis and Coupled Systems Laboratory, Conservatoire

More information

Simplified Buckling Analysis of Skeletal Structures

Simplified Buckling Analysis of Skeletal Structures Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation

More information

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines General Closed-form Analytial Expressions of Air-gap Indutanes for Surfaemounted Permanent Magnet and Indution Mahines Ronghai Qu, Member, IEEE Eletroni & Photoni Systems Tehnologies General Eletri Company

More information

Particle-wave symmetry in Quantum Mechanics And Special Relativity Theory

Particle-wave symmetry in Quantum Mechanics And Special Relativity Theory Partile-wave symmetry in Quantum Mehanis And Speial Relativity Theory Author one: XiaoLin Li,Chongqing,China,hidebrain@hotmail.om Corresponding author: XiaoLin Li, Chongqing,China,hidebrain@hotmail.om

More information

Special Relativity Entirely New Explanation

Special Relativity Entirely New Explanation 8-1-15 Speial Relatiity Entirely New Eplanation Mourii Shahter mourii@gmail.om mourii@walla.o.il ISRAEL, HOLON 54-54855 Introdution In this paper I orret a minor error in Einstein's theory of Speial Relatiity,

More information

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading 2th International LS-DYNA Users Conferene Constitutive Modeling(2) xperimental Investigation and F Analysis of Fiber Woven Layered Composites under Dynami Loading Pavel A. Mossakovsky, Fedor K. Antonov,

More information

STRUCTURAL AND BEHAVIORAL OPTIMIZATION OF THE NONLINEAR HILL MODEL

STRUCTURAL AND BEHAVIORAL OPTIMIZATION OF THE NONLINEAR HILL MODEL THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 12, Number 3/2011, pp. 213 220 STRUCTURAL AND BEHAVIORAL OPTIMIZATION OF THE NONLINEAR HILL MODEL Tudor

More information

Sensitivity Analysis in Markov Networks

Sensitivity Analysis in Markov Networks Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores

More information

Bayesian Optimization Under Uncertainty

Bayesian Optimization Under Uncertainty Bayesian Optimization Under Unertainty Justin J. Beland University o Toronto justin.beland@mail.utoronto.a Prasanth B. Nair University o Toronto pbn@utias.utoronto.a Abstrat We onsider the problem o robust

More information

Vibration Control of Smart Structure Using Sliding Mode Control with Observer

Vibration Control of Smart Structure Using Sliding Mode Control with Observer JOURNAL OF COMPUERS, VOL. 7, NO., FEBRUARY Vibration Control of Smart Struture Using Sliding Mode Control with Observer Junfeng Hu Shool of Mehanial & Eletrial Engineering, Jiangxi University of Siene

More information

A Big Bang-Big Crunch Optimization Based Approach for Interval Type-2 Fuzzy PID Controller Design

A Big Bang-Big Crunch Optimization Based Approach for Interval Type-2 Fuzzy PID Controller Design A Big Bang-Big Crunh Optimization Based Approah or Interval Type-2 Fuzzy PID Controller Design Tuan Kumbasar Control Engineering Department Istanbul Tehnial University Istanbul, Turkey kumbasart@itu.edu.tr

More information

The homopolar generator: an analytical example

The homopolar generator: an analytical example The homopolar generator: an analytial example Hendrik van Hees August 7, 214 1 Introdution It is surprising that the homopolar generator, invented in one of Faraday s ingenious experiments in 1831, still

More information

Influence of transverse cracks on the onset of delamination: application to L-angle specimens. F. Laurin*, A. Mavel, P. Nuñez, E.

Influence of transverse cracks on the onset of delamination: application to L-angle specimens. F. Laurin*, A. Mavel, P. Nuñez, E. Influene of transverse raks on the onset of delamination: appliation to L-angle speimens F. Laurin*, A. Mavel, P. Nuñez, E. Auguste Composite strutures subjeted to 3D loading Wings Strutures under 3D loadings

More information

REGULATION AND INPUT DISTURBANCE SUPPRESSION FOR PORT-CONTROLLED HAMILTONIAN SYSTEMS. Luca Gentili Arjan van der Schaft

REGULATION AND INPUT DISTURBANCE SUPPRESSION FOR PORT-CONTROLLED HAMILTONIAN SYSTEMS. Luca Gentili Arjan van der Schaft REGULATION AND INPUT DISTURBANE SUPPRESSION FOR PORTONTROLLED HAMILTONIAN SYSTEMS Lua Gentili Arjan van der Shaft ASYDEIS University of Bologna viale Risorgimento 4136 Bologna Italy email: lgentili@deisuniboit

More information

University of Groningen

University of Groningen University of Groningen Port Hamiltonian Formulation of Infinite Dimensional Systems II. Boundary Control by Interonnetion Mahelli, Alessandro; van der Shaft, Abraham; Melhiorri, Claudio Published in:

More information

Determination of the reaction order

Determination of the reaction order 5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order

More information

Active Magnetic Bearings for Frictionless Rotating Machineries

Active Magnetic Bearings for Frictionless Rotating Machineries Ative Magneti Bearings for Fritionless Rotating Mahineries Joga Dharma Setiawan Abstrat Ative magneti bearing (AMB systems an support a rotor without physial ontat and enable users to preisely ontrol rotor

More information

Fig 1: Variables in constant (1+1)D acceleration. speed of time. p-velocity & c-time. velocities (e.g. v/c) & times (e.g.

Fig 1: Variables in constant (1+1)D acceleration. speed of time. p-velocity & c-time. velocities (e.g. v/c) & times (e.g. Proper veloity and frame-invariant aeleration in speial relativity P. Fraundorf Department of Physis & Astronomy University of Missouri-StL, St. Louis MO (November, 99) We examine here a possible endpoint

More information

ETNA Kent State University

ETNA Kent State University % { Eletroni Transations on Numerial Analysis. Volume 21, pp. 20-27, 2005. Copyright 2005,. ISSN 1068-9613. QR FACTORIZATIONS USING A RESTRICTED SET OF ROTATIONS DIANNE P. O LEARY AND STEPHEN S. BULLOCK

More information

Inference-based Ambiguity Management in Decentralized Decision-Making: Decentralized Diagnosis of Discrete Event Systems

Inference-based Ambiguity Management in Decentralized Decision-Making: Decentralized Diagnosis of Discrete Event Systems Inerene-based Ambiguity Management in Deentralized Deision-Making: Deentralized Diagnosis o Disrete Event Systems Ratnesh Kumar and Shigemasa Takai Abstrat The task o deentralized deision-making involves

More information

CSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM

CSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM CSC2515 Winter 2015 Introdu3on to Mahine Learning Leture 5: Clustering, mixture models, and EM All leture slides will be available as.pdf on the ourse website: http://www.s.toronto.edu/~urtasun/ourses/csc2515/

More information

Beams on Elastic Foundation

Beams on Elastic Foundation Professor Terje Haukaas University of British Columbia, Vanouver www.inrisk.ub.a Beams on Elasti Foundation Beams on elasti foundation, suh as that in Figure 1, appear in building foundations, floating

More information

DO PHYSICS ONLINE. SPECIAL RELATIVITY Frames of Reference

DO PHYSICS ONLINE. SPECIAL RELATIVITY Frames of Reference DO PHYSICS ONLINE SPACE SPECIAL RELATIVITY Frames of Referene Spae travel Apollo 11 spaeraft: Earth Moon v ~ 40x10 3 km.h -1 Voyager spaeraft: v ~ 60x10 3 km.h -1 (no sling shot effet) Ulysses spaeraft:

More information

A multiscale description of failure in granular materials

A multiscale description of failure in granular materials A multisale desription of failure in granular materials Nejib Hadda, François Niot, Lu Sibille, Farhang Radjai, Antoinette Tordesillas et al. Citation: AIP Conf. Pro. 154, 585 (013); doi: 10.1063/1.4811999

More information

An Evolutionary Algorithm for Constrained Optimization

An Evolutionary Algorithm for Constrained Optimization An Evolutionary Algorim or Constrained Optimization Tapabrata Ray, Tai Kang and Seow Kian Chye Center or Advaned Numerial Engineering Simulations Shool o ehanial and Prodution Engineering Nanyang Tehnologial

More information

Properties of Quarks

Properties of Quarks PHY04 Partile Physis 9 Dr C N Booth Properties of Quarks In the earlier part of this ourse, we have disussed three families of leptons but prinipally onentrated on one doublet of quarks, the u and d. We

More information

QUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1

QUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1 QUANTUM MECHANICS II PHYS 57 Solutions to Problem Set #. The hamiltonian for a lassial harmoni osillator an be written in many different forms, suh as use ω = k/m H = p m + kx H = P + Q hω a. Find a anonial

More information

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW P. М. Меdnis Novosibirs State Pedagogial University, Chair of the General and Theoretial Physis, Russia, 636, Novosibirs,Viljujsy, 8 e-mail: pmednis@inbox.ru

More information

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple

More information

Some Properties on Nano Topology Induced by Graphs

Some Properties on Nano Topology Induced by Graphs AASCIT Journal of anosiene 2017; 3(4): 19-23 http://wwwaasitorg/journal/nanosiene ISS: 2381-1234 (Print); ISS: 2381-1242 (Online) Some Properties on ano Topology Indued by Graphs Arafa asef 1 Abd El Fattah

More information

IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES

IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES Arthur P. Berkhoff University of Twente, Faulty of Eletrial Engineering, P.O. Box 217, 7 AE Enshede, The Netherlands email: a.p.berkhoff@el.utwente.nl

More information

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena Page 1 of 10 Physial Laws, Absolutes, Relative Absolutes and Relativisti Time Phenomena Antonio Ruggeri modexp@iafria.om Sine in the field of knowledge we deal with absolutes, there are absolute laws that

More information

3 Tidal systems modelling: ASMITA model

3 Tidal systems modelling: ASMITA model 3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level

More information

TENSOR FORM OF SPECIAL RELATIVITY

TENSOR FORM OF SPECIAL RELATIVITY TENSOR FORM OF SPECIAL RELATIVITY We begin by realling that the fundamental priniple of Speial Relativity is that all physial laws must look the same to all inertial observers. This is easiest done by

More information

On Chaotic Cartesian Product of Graphs and Their Retractions

On Chaotic Cartesian Product of Graphs and Their Retractions ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Siene Vol.5(03) No.,pp.86-90 On Chaoti Cartesian Produt of Graphs and Their Retrations M. Abu-Saleem Department of Mathematis,

More information

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution. arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat

More information

Combined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2

Combined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2 Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque

More information

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer Proeedings of the 4th IIAE International Conferene on Industrial Appliation Engineering 216 Stabilization of the Preision Positioning Stage Working in the Vauum Environment by Using the Disturbane Observer

More information

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum

More information

PHY 396 T: SUSY Solutions for problem set #12.

PHY 396 T: SUSY Solutions for problem set #12. PHY 396 T: SUSY Solutions or problem set #. Problem a: In priniple the non-perturbative superpotential o the theory may depend on the dual quark and antiquark ields q and q as well as the singlets Φ but

More information