Chapter 7: Application Issues

Size: px
Start display at page:

Download "Chapter 7: Application Issues"

Transcription

1 Chapter 7: Applcaton Issues hs chapter wll brefly summarze several of the ssues that arse n mplementaton of the Carpal Wrst. he specfc ssues nvolved n puttng ths wrst nto producton are both task and manpulator dependent and are therefore unque to the desred applcaton of the overall manpulator system. However general approaches to modelng and solvng task-dependent operatons are necessary. hs chapter wll focus on manpulator systems topcs that are general n nature and that can be addressed usng the knematc and dynamc wrst model developed n ths dssertaton. he followng four topcs wll be dscussed n ths chapter: ) Methods for provdng full orentatonal control on the Wrst structure. 2) Redundancy resoluton through optmal performance requrements n dexterty mechancal advantage and dynamc response. 3) Path plannng or trajectory synthess 4) rajectory synthess based on nput torque requrements or other dynamc constrants. Whle not exhaustve these topcs represent a varety of the ssues crtcal to manpulator applcaton. Performance of a manpulator through advanced control results from accurate modelng of the robot archtecture ts envronment and task. he topcs demonstrate how the Wrst model can be ncorporated wth ts ntended task to provde the performance requred n applcaton. 7. Full Orentatonal Capablty of the Carpal Wrst he Carpal Wrst s a three degree-of-freedom devce provdng two orentatonal degrees of freedom ptch and yaw and a thrd translatonal degree of freedom called plunge. hese three degrees of freedom result from the knematc constrants of ts parallel archtecture and are fully controlled through the three nput jont angles. A typcal manpulator conssts of two sectons or structures one the arm (called regonal structure) whch s used to generally poston the tool and the other the wrst (called the local structure) whch provdes the requred orentaton of the tool. In applcaton the Carpal Wrst can be mounted on the end of a robotc arm to become the local structure. he most general robotc manpulator provdes full pose control of ts tool a sx degree of freedom (dof) task. Here the arm of a sx degree-of-freedom robot should provde full three dof poston capablty whle the wrst should provde complete orentaton capablty also a three-dof task. When mounted on a typcal three dof manpulator arm the added Carpal Wrst would make the robot a sx degree of freedom devce. he Carpal Wrst however lacks full orentatonal control snce t cannot provde roll moton between ts basal and dstal plates. Dependng on the choce of arm capabltes and applcaton needs a number of possble arrangements could exst: he frst possble arrangement s to mount the Wrst on an arm that provdes two postonal and one orentatonal degrees of freedom. hs would result n a sx dof manpulator that provdes full pose control of the tool. However there are several factors that may lmt or prevent ths form of applcaton. Frst the type of arm specfed n ths system s non-standard snce current robotc systems expect roll orentaton to be performed by the wrst structure rather 69

2 Stephen L. Canfeld Chapter 7: Applcaton Issues 70 than the arm. Also ths arrangement requres the Carpal Wrst to provde one of the postonng degrees of freedom for the manpulator. he plunge or translaton moton of the Carpal Wrst s lmted n range of travel however and would not satsfy the postonng requrements of a typcal manpulator. A second arrangement could consst of the Wrst mounted on a standard three degree of freedom postonng arm. hs arrangement results n a sx dof manpulator that s not capable of provdng tool roll but wth one redundant degree of postonal freedom.. hs arrangement has many advantages that would encourage ts use n producton. It s smlar to the typcal manpulator systems that use a three dof postonng arm. Also ths arrangement does not alter the Carpal Wrst structure. Fnally the redundant degree of freedom can be used to ncrease manpulator performance through a refned control algorthm. he drawback of ths arrangement s obvously ts lack of ablty to provde tool roll. However many applcatons do not requre orentaton of the tool about ts pontng drecton. hese applcatons or tasks are called axssymmetrc and nclude processes such as arc weldng spot weldng deburrng applyng adhesves polshng grndng drllng spray-pantng scannng many pck-and-place jobs and mountng and fxturng tasks. In applcatons that do not requre orentaton about the tool axs (or orentaton about any axs fxed relatve to the tool) ths arrangement provdes the best combnaton of arm and Wrst capabltes. he thrd arrangement proposed s to mount the Wrst on a standard three dof postonng arm and provde an extra roll degree of freedom to the manpulator. hs results n a system wth full pose control (6 dof) plus a redundant degree of freedom (for a total of 7 dof). However complcatons arse n provdng the addtonal dof n such a way that t does not lmt the advantages provded by the Carpal Wrst. Addng the roll degree of freedom can be accomplshed n a number of ways. he strongest and least complcated approach s to provde a rotatonal degree of freedom located between the manpulator arm and the Wrst. hs would rotate the basal plate normal as: θ roll B Z (7.) where θ roll s the nput roll angle and B Z s the z axs of the basal frame the frame of the Wrst base. In order to provde pure roll about the dstal plate extra orentaton control must be provded by the Wrst such that only the component of θ roll about the D Z axs s seen by the tool and the D X and D Y components are removed. Gven an nput θroll at the base the tool would be rotated by: D D B θ = R 3 θ 3 (7.2) roll B j roll where: D R 3 s the rotaton between the basal and dstal plates expressed n tensor notaton and B j D roll θ are the components of orentaton at the dstal plate for = -3. θ he negatve = 2 components of D roll are ncluded n the path plannng requrements and are sent through the nverse knematcs of the Wrst. he advantage of ths form of roll control s

3 Stephen L. Canfeld Chapter 7: Applcaton Issues 7 prmarly the locaton of the roll axs. Locatng t at the base of the Wrst mnmzes the sze and weght constrants that are placed on the Wrst components resultng n a stronger jont and lghter Wrst. he prmary dsadvantage s that the arm archtecture must be modfed an opton that s not attractve to companes who would adapt the Carpal Wrst to current robotc arms. o avod modfyng the arm requrements the roll axs can be ntegrated wth the Carpal Wrst archtecture. hs approach seems to be the most acceptable to robot manufacturers for applcatons requrng the full sx dof pose control snce t provdes a wrst wth full orentatonal control. A conceptual desgn of ths four-degree-of-freedom Wrst has been prepared n response to requests from several robot manufacturers. he conceptual desgn s shown on the Wrst n Fg. 7. and n an exploded vew n Fg hs arrangement provdes a compact hghstrength hgh-resoluton roll axs that mantans the hollow tunnel characterstcs of the Carpal Wrst. Knematcally the conceptual desgn conssts of addng a roll axs on the dstal plate drected along the dstal plate normal ( D Z). he forward and nverse knematcs of ths four degree-of-freedom arrangement have been presented n the postonal knematc analyss and demonstrate no addtonal dffculty n control. he conceptual desgn for the added roll dof conssts of four basc elements: a motor a drve-shaft a gear reducer and a bearng support structure. A hollow-shaft pancake-type DC servo motor s mounted on the bottom of the dstal plate. A pancake-type 000: harmonc-drve gear reducer s mounted on the opposng sde of the dstal plate along wth a flange-type bearng retaner that carres all axal and radal loads. A hollow drve shaft runnng through the dstal plate connects the motor to gear reducer. Fgure 7.: Conceptual Roll Axs Desgn

4 Stephen L. Canfeld Chapter 7: Applcaton Issues 72 Bearng Flange Dstal Plate Harmonc Drve Reducer Hollow Drve- Shaft Hollow-Shaft DC Servo Fgure 7.2: Exploded Vew of Roll dof Addng a roll axs to the Wrst creates a unque knematc devce n that the structure s no longer parallel but a hybrd parallel-seral devce. Hybrd manpulator devces have been dscussed n the lterature prmarly n devces that employ n-parallel actuaton. Snce the advantages of the Carpal Wrst stem from ts parallel archtecture addng an axs n seres would create a potental lmtng factor. However ths hybrd arrangement avods the major dsadvantage of seral manpulators; accumulatng moment loads on jont axes due to lnks farther along the chan. In the proposed Wrst hybrd arrangement the roll axs s placed at the end of the manpulator and thus does not see moment loads from any other lnk. A fnal approach to provdng full pose control wth the Wrst s to requre the tool to provde the roll moton. hs would be logcal n many applcatons for example a drllng tool or a rotary drvng tool. In effect ths approach s smlar to addng the roll on the dstal plate. 7.2 Redundancy Resoluton When matng the Carpal Wrst wth a typcal three degree-of-freedom postonng arm one redundant degree of freedom wll result due to the Wrst s plunge capablty. In ths case the redundancy s apparent from physcal consderatons of the manpulator system dof s whch can be thought of as postonng and orentng components. Mathematcally redundances are demonstrated n the system Jacoban matrx J where the rows of J represent output degrees of freedom (a maxmum of sx) of the manpulator and the columns represent system nputs or controlled jonts. hus a square n x n Jacoban of full rank represents a non-redundant system wth n degrees of freedom whle a non-square n x m wth m>n ( short ) Jacoban represents a

5 Stephen L. Canfeld Chapter 7: Applcaton Issues 73 redundant manpulator system wth (m-n) redundances. For a manpulator that does not requre roll the Carpal Wrst can be attached to a three dof arm and the resultng Jacoban would be a 5 x 6 matrx demonstratng 6-5 = redundant degree of freedom. he redundancy resultng n ths system s of nterest snce t provdes a free degree of freedom (.e. one that s not rgdly defned from the knematc path specfcatons) that can be used to mprove the knematc and dynamc performance of the system. he general approach for takng advantage of ths free dof s to specfy some performance crtera that provdes an addtonal constrant on the system at the knematc level thus resolvng the redundancy. hs wll be demonstrated on three levels for a Wrst-ntegrated manpulator system: One approach wll optmze on dexterty characterstcs at the nstantaneous knematcs level another wll provde a desred mechancal advantage or force control output also at the velocty level. he fnal approach wll optmze on dynamc performance at the acceleraton level Redundancy Resoluton Lterature Revew Usng manpulator system redundances to enhance knematc and dynamc performance has been studed by many researchers. One of the most common uses for redundances n robotc systems s obstacle-avodance. hs area has been studed by a number of researchers ncludng Macejewsk and Klen 985 Krcansk and Vukobratovc 986 Nearchou and Aspragathos 996 and McLean and Cameron 996. Whtney 969 demonstrated mprovements n rate resoluton based on the pseudo-nverse of the Jacoban matrx. hs approach has been refned and mproved by several other researchers ncludng Hollerbach and Suh 985 Nakamuraa and Hanafusa 987 Klen et al. 995 and Roberts and Macejewsk 996. Usng the system redundancy to mnmze the requred actuator nput energy was demonstrated by Vukobratovc and Krcansk 984 whle usng the system redundancy to mnmze nstantaneous jont torques was studed by Hollerbach and Suh 985. Usng redundancy to optmze the mechancal advantage of the manpulator was demonstrated by Dubrey and Luh 986 and Park et al Fnally Km and Rastegar 997 used redundancy n avodng actuatng torques that may excte the natural modes of vbratons of the manpulator system Optmzng on Dexterty In ths secton an approach for redundancy resoluton of the Carpal Wrst based on dexterty crtera s demonstrated. he dexterty crtera used n ths technque s formulated from the dexterty defntons presented by Soper et al. 997 Canfeld et al. 997 and the nstantaneous knematc analyss developed for the Carpal Wrst n Canfeld et al. 996 and the materal presented n Chap. 4. he dexterty defnton presented n ths work s based on the dea of the relatve stretchng between the nput and output velocty vectors defned by the Jacoban matrx: v D = = (7.3) ω µ

6 Stephen L. Canfeld Chapter 7: Applcaton Issues 74 where D s the dexterty µ s the mechancal advantage and the nput and output veloctes are related through the Jacoban as: v = Jω. (7.4) Expandng D: and then solvng for ω from Eq. 7.4 Substtutng ths result n Eq. 7.5 gves: D= D= ( v v) 2( ) ω = J 2 ω ω (7.5) v (7.6) ( ( ) ) v v v J J v. (7.7) hs equaton gves a general equaton for the dexterty of a manpulator as a functon of workspace poston n a desred drecton. Note that t assumes a square full rank (nonredundant non-sngular) Jacoban. When m>n a non-sngular Jacoban requres the rank of J to equal n. In order for the output velocty to be consstent n Eq. 7.4 t must be n the column space of J. When ths system of equatons s consstent the nput velocty s gven by the pseudo-nverse of J as: ω =( J J) J v he dexterty D s then found by substtutng Eq. 7.8 nto Eq. 7.5: D= [( ) ] ( ) v v v J J J J J J v (7.8). (7.9) Lookng at Eq. 7.4 a polar decomposton can be used to decompose the Jacoban nto two matrx operators: a matrx R whch stretches and a matrx Q whch rotates the nput velocty vector. hs s demonstrated for the m x n Jacoban for two cases (Soper et al. 997): m n: J = QR (7.0) where R s postve defnte R = ( J J) m < n: where R s postve defnte R = ( JJ ) 2 2 and Q s orthonormal Q= JR. J = RQ (7.) and Q s orthonormal Q= R J. he matrx R s of prmary nterest for the dexterty defnton because t contans the relatve stretchng nformaton between the nput and output veloctes. Consder the egenvalue decomposton of R: R = L L (7.2)

7 Stephen L. Canfeld Chapter 7: Applcaton Issues 75 where Σ s the dagonal matrx of the egenvalues of R. he elements of Σ represent the extremes n stretchng v relatve to ω and wll be called the prncpal dextertes. he assocated egenvectors wll be called the prncpal drectons. for m n: QRω = v (7.3) defnng y as y = Q v gves: Rω = y (7.4) and Eq. 7.3 becomes: D = y ω (7.5) whch s the rato of maxmum to mnmum prncpal dextertes. he drectons assocated wth these prncpal dextertes y prncpal are the egenvectors of R (u )rotated by Q: y prncpal = Qu (7.6) Fnally for the case where m<n the same approach s appled: RQω = v (7.7) and defnng y = Qω gves: Ry = v (7.8) and Eq. 7.3 becomes D = v y (7.9) whch s agan the rato of maxmum to mnmum prncpal dextertes. Now the drectons of the prncpal dextertes are the drectons of the egenvectors u. he polar decomposton for the case m<n allows the dexterty to be calculated as gven n Eqs. 7.5 and 7.9: D= ( ( ) ) v v v R R v. (7.20) Wth the dexterty equatons developed dexterty can be optmzed to resolve the redundant plunge degree of freedom n the Wrst. he approach wll proceed as follows. Frst the optmal dexterty wll be quantfed for example as the largest stretchng possble between the output and nput veloctes whch s gven n the maxmum prncpal dexterty. he above dexterty defntons wll then be used to fnd the maxmum prncpal dexterty and ts drecton. Wth ths Eq. 7.8 determnes the plunge velocty component of the vector v whch s then mposed on the knematcs to provde the addtonal constrant needed for redundancy resoluton.

8 Stephen L. Canfeld Chapter 7: Applcaton Issues Optmzng for Mechancal Advantage Usng mechancal advantage as the crtera n elmnatng redundancy follows the procedure for maxmzng dexterty. From the dexterty defnton proposed by Soper et al. 997 and repeated n Eq. 7.3 above the mechancal advantage s nversely related to dexterty. herefore the greatest mechancal advantage wll occur n the drecton of mnmum prncple dexterty. Equaton 7.8 s used agan to determne the plunge velocty component of v whch s then added as a constrant equaton n the knematcs provdng fully defned jont parameter specfcatons Optmzng for Dynamc Performance Optmzng the manpulator dynamc performance provdes a thrd method of redundancy resoluton. hs approach conssts of expressng the jont space coordnate motons as trajectores. hese trajectores are a functon of tme and coeffcents known as trajectory parameters. he trajectory parameters are determned from the knematc and dynamc constrants of the user-specfed tool-task and path. For jont space trajectores contanng k trajectory parameters one degree of redundancy adds an addtonal k parameters that can be used to satsfy path/task requrements. o optmze on dynamc performance there must be a suffcent number of total parameters to exceed the knematc moton requrements. he remanng parameters are avalable to satsfy dynamc requrements. hs s demonstrated n full detal n the followng sectons. 7.3 Path Plannng In operaton a manpulator typcally moves through a specfed path or trajectory. he trajectory defnes the tme hstory of poston velocty and acceleraton of the manpulator. he trajectory exsts n both tool space and jont space (whch are related through the manpulator knematcs) but s generally specfed n tool space. he tool space specfcatons are commonly user specfed n terms of goal poses wth some added path crtera. For example a smple pckand-place task has an ntal pck pose specfed as a goal a fnal place pose specfed as a goal and path crtera whch may nclude the speed of travel general obstacle avodance and tmng events. A general approach to path plannng appled to the Carpal Wrst s presented that can be extended to nclude the arm knematcs and thus the entre manpulator. hs approach s descrbed n the followng eght steps: ) he desred output path s gven n tool-space coordnates. hs ncludes the path crtera event tmng etc. 2) he output path s descrbed wth a fnte set of precson postons called va ponts (Crag 989). 3) he precson postons or va ponts are descrbed n jont space as: ( ϑ θ θ ) ( α β ) 2 3 = f p (7.2)

9 Stephen L. Canfeld Chapter 7: Applcaton Issues 77 where the jont angles θ -θ 3 are determned from α β and p (the tool-space Wrst coordnates) through the functon f whch represents the nverse knematcs presented n Chap. 2. 4) Wth the precson ponts specfed n jont space trajectores for the jonts between the precson postons can be descrbed usng trajectory curves to meet the specfc trajectory crtera. 5) he trajectory curves are wrtten as a functon of ther trajectory parameters and tme: n θ = k j f j j= () t (7.22) where n represents the degree order or number of harmoncs dependng on the type of path chosen. 6) Moton crtera are then gven as parameters on the trajectory posed as equatons. he type of moton and degree of crtera determne the type and degree of the trajectory path and ts parameters. 7) he jont trajectory equatons (Eq. 7.22) are substtuted nto the moton crtera equatons and the resultng system of equatons s solved for the trajectory parameters. 8) Wth the trajectory parameters the jont trajectores are known and are ready to be sent to the control system to result n the desred tool path. For example a specfed path and set of precson postons can be gven requrng contnuous poston velocty and acceleraton values between all va ponts wth specfc end condtons. hs crtera can be expressed as: θ( t fk ) = θ θ k ( t fk ) = θ θ k ( t fk ) = θ k (7.23) θ( t fk ) = θk θ ( t fk ) = θ k + θ( t fk + ) = θ k + he jont space angles can be suffcently descrbed usng a ffth-order polynomal: θ () t = a5 t + a4 t + a3 t + a2 t + a t + a0 (7.24) Solvng for the k th nterval gves:

10 Stephen L. Canfeld Chapter 7: Applcaton Issues 78 a a a a a a = θ = θ k k = θ k 2 ( ) t ( ) 20θ 20θ 8θ + 2θ 3θ θ = 2t 2 k+ k o k+ k k+ k k+ k+ 3 k+ ( ) t ( ) 30θ 30θ + 4θ + 6θ 3θ 2θ = 4 2t 2 k o k+ k+ k k+ k k+ k+ k+ ( ) t ( ) 2θ 2θ 6θ + 6θ θ θ = 5 2t 2 k+ k o k+ k k+ k k+ k+ k+ t t t. (7.25) 7.4 Path Plannng Based on Input orque Requrements / Constrants rajectory synthess has been demonstrated n meetng the knematc requrements of the desred tool path. he general procedure conssted of descrbng the jont space varables as a functons of tme and trajectory parameters. hen the knematc requrements were cast as constrant equatons and used to solve the trajectory parameters thus descrbng the jont space trajectores. In ths dscusson ths general approach wll be extended to nclude dynamc / nput actuaton requrements as well as the knematc requrements placed on the jont space trajectores. Summarzng the results from Chap. 5 the equatons of moton of the Carpal Wrst manpulator can be wrtten as: = A θθ + B θθθ + Cθ (7.26) where: M s the th nput motor moment () ( ) () M A(θ) s the mass matrx representng nertal components of the wrst and tool B( θθθ ) represents the Corols and centrfugal nerta terms and C(θ) represents the gravtatonal terms. Let the jont space trajectores be descrbed as: n general or for example as a polynomal: ( ) r k θ = f k a k t k= 0 r k θ = a k t k= 0 (7.27) (7.28)

11 Stephen L. Canfeld Chapter 7: Applcaton Issues 79 where: θ are the jont space coordnates cast n the form of a polynomal (n ths example) a k k=0-r are the r+ trajectory parameters and t s tme. Let s be the number of knematc constrants or requrements gven by the specfed tool path. hus s trajectory parameters are requred to satsfy the s specfed end condtons. Snce a total of r+ trajectory parameters are avalable n Eq or 7.28 above the extra (r+ - s) trajectory parameters can be used to satsfy dynamc / motor torque requrements. he torque requrements can be represented as: m ( ) ( ) M = f p t ; l = r+ k (7.29) l j j j= where M l s the h torque constrant for l = -> (r+-k) of the th motor ( j ) f j p t represents the m functons makng up the th order torque constrant equaton and are gven along wth the tool-path specfcatons. All that remans s solvng the r+ trajectory parameters a k. he soluton s relatvely straghtforward. It conssts of takng the frst and second dervatves of Eqs (or 7.28): k k ( ka k t ) f k ( a k t ) θ = r k= 0 k 2 k k k [( kk ( ) ak t ) f k ( ak t ) ( kak t ) f k ( ak t )] θ = + r k= 0 n general or the dervatves of the polynomal example: θ = r k= 0 k k θ = kk a t r k= 0 ka t ( ) k k 2 (7.30) (7.3) and then substtutng these nto Eq he moments n Eq are also replaced wth the motor torque requrements Eq Fnally the resultng system of equatons s solved for the (r+-k) trajectory parameters gvng the jont space trajectores. 7.5 Optmal Path Plannng Based on Input orque Requrements he procedure demonstrated above can be altered to generate an optmzaton algorthm. he value of ths les n meetng general motor torque requrements that are not easly satsfed wth a fnte number of constrant equatons. hus the procedure above s altered at Eq to wrte the motor requrements as a functon to be optmzed followed by a gven set of constrant

12 Stephen L. Canfeld Chapter 7: Applcaton Issues 80 equatons. Optmzaton follows yeldng a set of trajectory parameters to best meet the desred motor torque profle whle satsfyng knematc moton constrants. 7.6 Summary of Applcaton Issues hs chapter has addressed several of the ssues that wll arse n applcaton of the Carpal Wrst. hs chapter s not ntended to cover all mplementaton ssues Rather ts purpose s to demonstrate the applcaton of the Carpal Wrst model knematcs and dynamcs n resolvng applcaton conflcts and n mprovng performance. A fourth degree of freedom can be added to the Carpal Wrst to provde full orentatonal capabltes he solutons presented however wll ncrease costs and are thus ntended only n applcatons requrng the addtonal degree of freedom. Redundancy resoluton technques were demonstrated at both the nstantaneous knematc and the dynamc level. Results from the applcaton of one of these specfc technques can be seen n Fg. 4.8 a mnmum dexterty plot of the Carpal Wrst. Secton and able 4.3 demonstrate the sgnfcant mprovement n Wrst dexterty due n part to the redundant degree of freedom. hs mprovement n dexterty over conventonal robotc wrsts (or the Carpal Wrst operated n fxed-plunge mode) s provded at the cost of one addtonal actuator. rajectory Synthess approaches for the Carpal Wrst and assocated arm were presented. Path plannng s requred n all manpulator systems. he technques gven n ths secton demonstrate how the knematcs and dynamc equatons of moton can be derved to optmze Wrst and manpulator performance under varous operatng crtera.

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Iterative General Dynamic Model for Serial-Link Manipulators

Iterative General Dynamic Model for Serial-Link Manipulators EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Advanced Mechanical Elements

Advanced Mechanical Elements May 3, 08 Advanced Mechancal Elements (Lecture 7) Knematc analyss and moton control of underactuated mechansms wth elastc elements - Moton control of underactuated mechansms constraned by elastc elements

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

coordinates. Then, the position vectors are described by

coordinates. Then, the position vectors are described by Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics) CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

The Analysis of Coriolis Effect on a Robot Manipulator

The Analysis of Coriolis Effect on a Robot Manipulator Internatonal Journal of Innovatons n Engneerng and echnology (IJIE) he Analyss of Corols Effect on a Robot Manpulator Pratap P homas Assstant Professor Department of Mechancal Engneerng K G Reddy college

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists * How Strong Are Weak Patents? Joseph Farrell and Carl Shapro Supplementary Materal Lcensng Probablstc Patents to Cournot Olgopolsts * September 007 We study here the specal case n whch downstream competton

More information

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

More information

Lecture 8 Modal Analysis

Lecture 8 Modal Analysis Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Physics 111: Mechanics Lecture 11

Physics 111: Mechanics Lecture 11 Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton

More information

Pivot-Wheel Drive Crab with a Twist! Clem McKown Team November-2009 (eq 1 edited 29-March-2010)

Pivot-Wheel Drive Crab with a Twist! Clem McKown Team November-2009 (eq 1 edited 29-March-2010) Pvot-Wheel Drve Crab wth a Twst! Clem McKown Team 1640 13-November-2009 (eq 1 edted 29-March-2010) 4-Wheel Independent Pvot-Wheel Drve descrbes a 4wd drve-tran n whch each of the (4) wheels are ndependently

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15 NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

More information

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

1 Derivation of Point-to-Plane Minimization

1 Derivation of Point-to-Plane Minimization 1 Dervaton of Pont-to-Plane Mnmzaton Consder the Chen-Medon (pont-to-plane) framework for ICP. Assume we have a collecton of ponts (p, q ) wth normals n. We want to determne the optmal rotaton and translaton

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

CHAPTER III Neural Networks as Associative Memory

CHAPTER III Neural Networks as Associative Memory CHAPTER III Neural Networs as Assocatve Memory Introducton One of the prmary functons of the bran s assocatve memory. We assocate the faces wth names, letters wth sounds, or we can recognze the people

More information

1 GSW Iterative Techniques for y = Ax

1 GSW Iterative Techniques for y = Ax 1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn

More information

PHYS 705: Classical Mechanics. Canonical Transformation II

PHYS 705: Classical Mechanics. Canonical Transformation II 1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m

More information

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

More information

AP Physics 1 & 2 Summer Assignment

AP Physics 1 & 2 Summer Assignment AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013 Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique Outlne and Readng Dynamc Programmng The General Technque ( 5.3.2) -1 Knapsac Problem ( 5.3.3) Matrx Chan-Product ( 5.3.1) Dynamc Programmng verson 1.4 1 Dynamc Programmng verson 1.4 2 Dynamc Programmng

More information