A Guidance Law for a Mobile Robot for Coverage Applications: A Limited Information Approach
|
|
- Quentin Norman
- 6 years ago
- Views:
Transcription
1 Proceedings of the Third International Conference on Adances in Control and Optimization of Dynamical Systems, Indian Institute of Technology anpur, India, March 1-15, 014 Frcsh1.6 A Guidance Law for a Mobile Robot for Coerage Applications: A Limited Information Approach Twinle Tripathy Arpita Sinha Interdisciplinary Programme in Systems and Control Engineering, IIT Bombay, Mumbai , India twinle.tripathy@sc.iitb.ac.in, asinha@sc.iitb.ac.in Abstract: In this paper a guidance law has been proposed to sole applications related to coerage of a pre-defined area using a mobile robot. In many situations, it may not be possible to measure all the information due to practical limitations or from computational point of iew. So, we hae used only range information to sole the aboe specified problem. The target point is considered to be stationary in this particular case. Now, for any initial condition for which the elocity ector of the robot is perpendicular to the line of sight with respect to the stationary target, the autonomous agent displays interesting behaiours which hae been analysed theoretically. All the analytical results obtained here hae been alidated through MATLAB simulations also. eywords: Coerage problems, minimal information, autonomous mobile robots, unicycle model, root locus. 1. INTRODUCTION Oer the past few decades the problem of coerage of an area by a mobile robot has been a ery challenging problem. It finds seeral applications these days. Most common applications would be acuum cleaning, painting, humanitarian demining, foraging, farm irrigation and so on. Some of the more crucial applications would be monitor the area around a stationary point lie a building. The coerage problem refers to the problem of guiding an autonomous agent to coer an unnown or predefined area. Seeral techniues hae been proposed to achiee coerage, which ary depending on the application. Sureillance based systems, discussed in Dhillon et al. 00), Lai et al. 007) and Wang et al. 008), use wireless sensor networs to sole the coerage problem. The sensor nodes are capable of behaing autonomously. In Dhillon et al. 00) and Wang et al. 008), the coerage achieed is aerage, but more emphasis is laid on the placement of sensors for effectie coerage; Dhillon et al. 00) hae based their wor on fixed sensor nodes. Lai et al. 007) discuss sureillance systems based on wireless sensor networs. In Lai et al. 007), the deployed sensors are diided into disjoint subsets of sensors, or sensor coers, such that each sensor coer can coer all targets and wor in turns. A sureillance system based on wireless sensor networs has been discussed by Yan et al. 00). The system proides a differentiated sureillance system which deploys sensors nodes in appropriate areas depending on the degree of sureillance reuired. But, the wireless sensor networs suffer from the problem of energy reuirements which increase with the number of sensor nodes used in the networ. There are seeral algorithms which are used to coer unnown areas. For example, a lot of wor has been done on map based coerage techniues see Wong et al. 00), Zelinsy et al. 199), Stachniss et al. 00) and Rutishauser et al. 009) and references therein) which wor by diiding the area into cells and then assigning alues to the cells based on the presence of obstacles or free space. Wong et al. 00) used topological maps to sole the coerage problem and their results are erified by simulation tests which show that oer 99% of the surface is coered. The topological maps, represent the enironment as graphs where landmars are nodes, and edges represent the connectiity between the landmars. The algorithm achiees coerage with single mobile agent. Zelinsy et al. 199) also use maps for coerage but with a different approach. The solution ensures complete coerage by the use of a distance transform path planning methodology. The solution was simulated and implemented on an autonomous robot called Yamabico to gie satisfactory performance. Stachniss et al. 00) introduced the concept of coerage maps where each cell of a gien grid corresponds to the amount of a cell which is coered by an obstacle; the coerage maps are improement of occupancy grids which are based on the assumption that each cell is either occupied or free. The model presented in the paper allows updation of the coerage maps upon input obtained from sensors. Rutishauser et al. Rutishauser et al. 009) sole the problem of collaboratie coerage using a swarm of networed miniature robots again by the use of grid based methods. For a multi-agent system, the problem of coerage, addressed by Bataia et al. 009), has been soled with high accuracy for a semi- structured enironment. The approach has the adantage of an operator driing the outline of a desired coerage area as input to a coerage generation algorithm. Acar et al. 001) hae 47
2 Tx T, y T ) R max ρr T R min 1 ρ)r VTx V T, y V T ) α θ Fig. 1. Region to be coered Fig.. Coerage obtained by an Autonomous Agent soled the coerage problem by the decomposition of the enironment, by the use of oronoi diagrams. In this paper we hae proposed a guidance strategy that uses only range information. The transmission of minimal information reduces the bandwidth reuirement. This leads to energy conseration. Priliminary results on the guidance strategy hae been done by Tripathy et al. 01). In the wor by Tripathy et al. 01), we hae shown that the coerage strategy is stable for eery initial condition. The problem discussed in this paper is formulated as follows. The autonomous agent is modeled as an unicycle. We assume, the region to be coered, is circular or annular around the target as shown in Fig. 1. The agent uses only distance information to coer the annular region. The annular region is specified by R min and R max. It is shown that R min and R max depend on the initial conditions. Using this, we hae deried: The trajectory of the agent under the initial condition for which the elocity ector is perpendicular to the line of sight with respect to the target. The initial conditions necessary to achiee coerage of an entire circular dis of any radius. This paper is organized as follows. In Section and, we hae analysed the problem for a particular situation. Simulations results are presented in Section 4 and Section 5 concludes the paper and discusses the areas on which further wor can be done in future.. COVERAGE STRATEGY The paper addresses the coerage problem as shown in Fig. 1. Let the target be at x T, y T ) and is stationary at any instant t > 0, position of agent is gien by xt), yt)), Ax, y) Fig.. Engagement geometry Reference as shown in Fig.. The inematics of agent is as shown below: ẋt) = cos αt) 1) ẏt) = sin αt) ) αt) = ut) ) where is the linear elocity, which is a constant for our case, αt) is the heading angle and ut) is the control input. In Fig., rt) is the line-of-sightlos) distance and θt) is the LOS angle. Since only range measurement is feasible, agent nows only instantaneous rt) and not θt). Let V T denote a irtual point, called the irtual target which is at a distance ρrt) from the target along the LOS, where ρ 0, 1). In the control law, the lateral acceleration of the autonomous agent is made proportional to the distance between the irtual target and the agent, gien by x V T xt)) + y V T yt)) ) 1/ = 1 ρ)rt). So, the algorithm is defined in terms of a M t) as: a M t) = h1 ρ)rt) 4) αt) = a M t) = rt) 5) where h is the constant of proportionality between a M and the distance, 1 ρ)r, and = h1 ρ). The parameter h is a system gain. The concept of the irtual target increases the flexibility of the algorithm as we get another parameter ρ to be aried as needed. The algorithm generates arious coerage patterns surrounding the stationary target. One of the patterns generated is shown in Fig.. In the next section, the control algorithm is analysed for area coerage applications. Consider the Fig., the LOS between the agent and the target can be charaterised as follows: r = ṙ = cosα θ) 6) θ = r θ = sinα θ) 7) Let us define φ = α θ. By combining 5)-7) and soling the integration, we establish a relationship between r and φ as, r sin φ + r = r 0 sin φ 0 + r 0 8) In the aboe euation, r 0 and φ 0 are the initial conditions. We define a ariable which contains all the terms corresponding to the initial conditions, such that, 48
3 r Bφ = π) point A and thereafter it decreases as r < 0 in the path A D C. Hence, R max occurs at Aφ = π ) and R min occurs at Cφ = π ). Cφ = π ) Fig. 4. θ, r ) plot C B r D A Fig. 5. Engagement trajectory 1 in θ, r ) space θ Dφ = 0) r0 sinφ 0 ) = θ Aφ = π ) r B C A Fig. 6. Engagement trajectory in θ, r ) space ) + r0 9) r r sin φ + = 0 10) By setting 6) to zero, we obsere that r can tae its maximum and minimum alues that is R min and R max when φ euals either π or π. Next we loo into how the inematics aries in the θ, r ) plane. Upon suaring and adding 6) and 7), we find that the instantaneous θ and r lie on a circle in the θ, r ) plane as shown in Fig. 4. By the use of 6) and 7), we can calculate the alue of φ corresponding to each set of instantaneous r and θ. Point A corresponds to θ, r ) =, 0), thus, φ = π. Similarly, we get that at B, φ = π, at C, φ = π and at D, φ = 0. Hence, the clocwise moement about the circle shown in Fig. 4 gies the direction of increasing alues φ. Combining 5) - 7), we get, dφ dt = + r ) 11) r The rate of change of φ with respect to time aries depending on the sign of. Considering the two cases separately, we hae: Case 1: > 0 The aboe euation indicates that for > 0 we get φ > 0. This corresponds to clocwise moement in the θ, r ) plane as shown in Fig. 5. Now, r > 0 in the path C B A which means instantaneous r increases in this path upto θ Case : 0 r and θ are continuous functions of φ. By re-arranging 10), we get sin φ = r ) r. Hence, from 7) for all alues of 0, θ 0. This indicates that the engagement geometry can exist only in the first and fourth uadrants in the θ, r ) plane. So, the only feasible inflexion point is φ = π. Hence, considering the continuity of the θ and r functions, we conclude that the θ, r ) plot oscillates about φ = π as shown in Fig. 6. Now, r > 0 in the path A B A which means instantaneous r increases in this path till the second time it reaches point A. Thereafter, it decreases as r < 0 in the path A C A till it comes to the point A for the second time. Hence, in this case, both R max and R min occur at Aφ = π ). Thus, by using the coerage control algorithm gien in 5), we get two types of engagement trajectories in r, θ ) plane, depending on or the initial conditions r 0, φ 0 ). By differentiating 6) and 7), we get: r = θ φ 1) θ = r φ 1) At the points B and C shown in Fig. 6, r and θ become zero simultaneously, which implies φ = 0. Using this condition in 11), we get the following result for the instantaneous alue of r: r = ) 1 14) where = ). Substituting 14) in 10) and further soling, we get an expression of instantaneous φ as, sin φ = 4 ) 1 15) En. 15) gies a bound on = ) as sine function can tae alues between -1 to 1. Hence, the minimum alue of can be obtained by setting sin φ = 1, which gies = =. Let us denote the lower alue of by: br =. Thus, coerage becomes feasible for negatie alues of only if br 0. Next we address two problems: If we are gien an initial condition, what will be the coerage area and if we are gien an area to be coered, what should be the initial conditions to achiee the coerage? Gien an initial condition r 0, φ 0 ), we hae the following result: Theorem 1. Tripathy et al. 01)) Gien a set of initial conditions r 0, φ 0 ) for an agent with inematics 1)-) and a control law gien by 5), coerage will be achieed for a region gien by the maximum radius R max and minimum radius R min which are gien as: R max = az 1 + bz 16) R min = az + bz 4 17) 49
4 Theorem. Tripathy et al. 01)) For r 0 = and φ 0 = π, the trajectory of an agent with inematics 1)- ) and a control law gien by 5), is a circle of radius. Ax, y) α φ = π θ Tx T, y T ) Reference The detailed proof of this theorem has been coered in the paper by Tripathy et al. 01). Proposition. If φ 0 = π and r 0 [0, ], then R min = r 0. Fig. 7. Agent and target initial engagement for φ 0 = π where 1 + i a =, if br 0 1, if 0 < < 1 i b =, if br 0 1, if 0 < < z 1 = ) z = ) 1 4 z if br 0 z = ) if 0 < < z 1 if br 0 z 4 = ) if 0 < <. The detailed proof of the theorem is gien in the wor by Tripathy et al. 01).. AREA OF COVERAGE In this section we hae considered the special case when the initial conditions are selected such that φ 0 = π. In physical sense, this is the configuration when that the initial heading angle is perpendicular to the LOS between the agent and the stationary target as shown in fig. 7). Critical Initial Position: Critical initial position refers to that alue of r 0, φ 0 ) for which R max and R min are eual. This maes the agent moe on a circular path around the stationary target for all time. Also, it is interesting to note that for these initial coordinates, we get R max = R min = r 0. Proof. Substituting φ 0 = π in 9), we find out that for r 0 [0, ] 0 and otherwise > 0. Let us consider the case r 0 [0, ]. From Section, we now that for 0, R max and R min occur at φ = π. Now, we substitute φ 0 = π in 8) and sole it to get the following roots: r 1 = r 0 18) r = r 0 + r0 + 19) 4 One of the roots is R max and the other one is R min. The negatie roots hae been ignored. r 1 r = r 0 1 r ) From 0) it can be shown that r 1 r 0 for r 0. Hence, R min = r 0 for r 0 [0, ]. Proposition 4. If φ 0 = π and r 0 [, ], then R max = r 0. Proof. In the proof gien for Proposition, we hae seen that for r 0 [0, ] 0. Hence, in this region both R max and R min occur at φ 0 = π. So, substituting the alue of φ 0 in 8), we get the roots 18) and 19). Again, we hae ignored the negatie root. Using 0), it can be shown that r 1 r 0 for r 0 [, ]. Hence, R max = r 0 for r 0 [, ]. Theorem 5. The agent coers a maximum radius of r 0 if it starts at r 0 with heading φ 0 = π. Proof. For r 0 [, ], it follows directly from Proposition 4 that R max = r 0. From Section, we now that R max always occurs at φ 0 = π irrespectie of the sign of. Substituting the alue of φ 0 in 8), we get the positie roots as, r 1 = r 0 1) ) 1 r = ) ) Using simple calculations, it can be shown that r 1 > r for eery r 0 >. Hence, we conclude that for r 0, R max = r
5 π and r0 [0, ], the trajectory of the agent is bound by a circle of radius. Corollary 6. Gien φ0 = Proof. In Proposition 4, we found out that for r0 p [, ] Rmax = r0. Hence, Rmax becomes maximum when r0 is maximum. Hence,in the gien range the maximum alue of Rmax is. For the case when p r0 [0, ], the expression for Rmax is gien by 19). The first deriatie of the expression is gien by 1 1 r0 r0 + 1 ). The first deriatie is always negatie, so Rmax decreases with respect to r0. Hence, Rmax is maximum for r0 = 0. Using 19), the maximum alue of Fig. 8. Trajectory plot: r0 = 45, φ0 = Fig. 9. r ersus t : r0 = 45, φ0 = 110π 180. Since, for the gien range of r0 the maximum alue of Rmax has an upper bound of, the trajectory is also bounded by a circle of radius. Theorem 7. To coer an entire region of radius, the Rmax is eual to agent must hae an initial condition gien by d, π ) or π 0, ) and = d. Proof. In Corollary 6, we deried that the maximum alue of Rmax is eual to which occurs at r0 = 0 and. Using 18) and 19), we calculated that for r0 = 0 and, the alue of Rmin is eual to zero. This proes that the agent coers the entire circular dis of radius π π for the initial configurations, ) and 0, ). 110π 180 Let us consider a circle of radius d which is desired to be coered. The alue of the controller gain should be chosen such that d =. Hence, as explained in the aboe paragraph, for = d the agent should start at π π d, ) and 0, ) in order to coer a circular area of radius d. 4. SIMULATION RESULTS In this section, we alidate the results obtained in Section through MATLAB simulations. The simulations are performed with the parameters h = 0.01, ρ = 0. and = 5. The control algorithm deeloped in 5) leads to formation of patterns around the target ensuring coerage of the area surrounding the target. 4.1 Critical Initial Position p For initial conditions, π ), the trajectory of the agent is circular with the target as the centre as gien by Theorem. This result has been erified by simulation as shown in fig. 10) 4. Range of Rmax and Rmin In Proposition, we deried p that Rmin aries linearly with respect to r0 for r0 [0, ]. This can be erified by the fig. 11 which has been obtained using MATLAB Fig. 10. Agent trajectory for critical initial position simulation. Similarly in Proposition 4, we found outthat p Rmax aries linearly with respect to r0 for r0 [, ]. This result has been shown in fig CONCLUSION In the paper, we hae analysed the control law 5) under initial conditions such that the elocity ector is perpendicular to the line of sight with respect to the stationary target. For this initial heading, the agent exhibits seeral interesting behaiours. For the critical initial initial condition gien in Theorem, the agent circles around the target. Hence, the agent can continuously monitor the target using only range information. The radius of the circular trajectory depends on the controller gain and speed and, hence, can be specified by the user. In Theorem 7, it has been shown that by adjusting the controller gain, 441
6 Fig. 11. Variation of Rmin with respect to r0 Fig. 1. Variation of Rmax with respect to r0 Acar, E. U., Choset, H., and Atar, P. N., Complete Sensor-based Coerage with Extended-range Detectors: A Hierarchical Decomposition in Terms of Critical Points and Voronoi Diagrams, Proceedings of the IEEW/RSJ International Conference on Intelligent Robots and Systems, , 001 Dhillon, S. S. and Charabarty,. Sensor Placement for Effectie Coerage and Sureillance in Distributed Sensor Networs Proceedings of IEEE Wireless Communications and Networing Conference, , 00. Lai, C. C., Ting, C. C., and o, R. S. An Effectie Genetic Algorithm to Improe Wireless Sensor Networ Lifetime for Large-Scale Sureillance Applications IEEE Congress on Eolutionary Computation, 5158, 007. Wang, W., Sriniasan, V., Wang, B., and Chua,. C. Coerage for Target Localization in Wireless Sensor Networs IEEE Transactions on Wireless Communications, Vol. 7, No., , 008. Yan, T., He, T., and Stanoic, J. A. Differentiated Sureillance for Sensor Networs First ACM Conference on Embedded Networed Sensor Systems, 51-6, 00. Tripathy, T. and Sinha, A. Guidance of an autonomous agent for coerage applications using range only measurement Proceedings of AIAA, Guidance, Naigation & Control Conference, 01, the agent can be made to coer an entire dis. The wor can be extended to find out the effect of the controller gain and speed. Also, more wor can be done to mae the controller robust. Since, the agent uses only onerange) parameter s data for the control law, the noise associated with the system reduces. The guidance law also uses minimal information for coerage. This leads to cost saing. It also minimises the energy consumption which is one of the major issues of today! REFERENCES Wong, S. C. and MacDonald, B. A. A topological coerage algorithm for mobile robots Proceedings of the IEEW/RSJ International Conference on Intelligent Robots and Systems,, , 00. Zelinsy, A., Jaris, R. A., Byrne, J. C., and Yuta, S. Planning paths of complete coerage of an unstructured enironment by a mobile robots Proceedings of International Conference on Adanced Robotics, 5-58, 199. Stachniss, C. and Burgard, W. Mapping and Exploration with Mobile Robots using Coerage Maps Proceedings of the IEEW/RSJ International Conference on Intelligent Robots and Systems, Vol. 1, , 00. Rutishauser, S., Correll, N., and Martinoli, A. Collaboratie coerage using a swarm of networed miniature robots Robotics and Autonomous Systems, Vol. 57, 51755, 009. Bataia, P. H., Roth, S. A., and Singh, S. Autonomous Coerage Operations In Semi-Structured Outdoor Enironments Proceedings of the IEEW/RSJ International Conference on Intelligent Robots and Systems, ,
Generating Patterns with a Unicycle
Generating Patterns with a Unicycle Twinkle Tripathy and Arpita Sinha Abstract The objective of the paper is to make a unicycle trace hypotrochoid like patterns using a simple range based control scheme.
More informationUNDERSTAND MOTION IN ONE AND TWO DIMENSIONS
SUBAREA I. COMPETENCY 1.0 UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS MECHANICS Skill 1.1 Calculating displacement, aerage elocity, instantaneous elocity, and acceleration in a gien frame of reference
More informationChapter 2: 1D Kinematics Tuesday January 13th
Chapter : D Kinematics Tuesday January 3th Motion in a straight line (D Kinematics) Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Short summary Constant acceleration a special
More informationCollective circular motion of multi-vehicle systems with sensory limitations
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seille, Spain, December 12-15, 2005 MoB032 Collectie circular motion of multi-ehicle systems with
More informationDO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.
DO PHYSICS ONLINE DISPLACEMENT VELOCITY ACCELERATION The objects that make up space are in motion, we moe, soccer balls moe, the Earth moes, electrons moe, - - -. Motion implies change. The study of the
More informationChapter 1: Kinematics of Particles
Chapter 1: Kinematics of Particles 1.1 INTRODUCTION Mechanics the state of rest of motion of bodies subjected to the action of forces Static equilibrium of a body that is either at rest or moes with constant
More informationCollective circular motion of multi-vehicle systems with sensory limitations
Collectie circular motion of multi-ehicle systems with sensory limitations Nicola Ceccarelli, Mauro Di Marco, Andrea Garulli, Antonello Giannitrapani Dipartimento di Ingegneria dell Informazione Uniersità
More information0 a 3 a 2 a 3 0 a 1 a 2 a 1 0
Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.
More informationFeb 6, 2013 PHYSICS I Lecture 5
95.141 Feb 6, 213 PHYSICS I Lecture 5 Course website: faculty.uml.edu/pchowdhury/95.141/ www.masteringphysics.com Course: UML95141SPRING213 Lecture Capture h"p://echo36.uml.edu/chowdhury213/physics1spring.html
More informationA possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k
A possible mechanism to explain wae-particle duality L D HOWE No current affiliation PACS Numbers: 0.50.-r, 03.65.-w, 05.60.-k Abstract The relationship between light speed energy and the kinetic energy
More informationWould you risk your live driving drunk? Intro
Martha Casquete Would you risk your lie driing drunk? Intro Motion Position and displacement Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Constant acceleration: A special
More informationDynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates
ACC 9, Submitted St. Louis, MO Dynamic Vehicle Routing with Moing Demands Part II: High speed demands or low arrial rates Stephen L. Smith Shaunak D. Bopardikar Francesco Bullo João P. Hespanha Abstract
More informationPHYS 1443 Section 004 Lecture #4 Thursday, Sept. 4, 2014
PHYS 1443 Section 004 Lecture #4 Thursday, Sept. 4, 014 One Dimensional Motion Motion under constant acceleration One dimensional Kinematic Equations How do we sole kinematic problems? Falling motions
More information4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion.
4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion. We now hae deeloped a ector model that allows the ready isualization
More informationA Geometric Review of Linear Algebra
A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. The order of presentation is unconentional, with emphasis on geometric intuition rather than
More informationDisplacement, Time, Velocity
Lecture. Chapter : Motion along a Straight Line Displacement, Time, Velocity 3/6/05 One-Dimensional Motion The area of physics that we focus on is called mechanics: the study of the relationships between
More information4-vectors. Chapter Definition of 4-vectors
Chapter 12 4-ectors Copyright 2004 by Daid Morin, morin@physics.harard.edu We now come to a ery powerful concept in relatiity, namely that of 4-ectors. Although it is possible to derie eerything in special
More informationVISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION
VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION Predict Obsere Explain Exercise 1 Take an A4 sheet of paper and a heay object (cricket ball, basketball, brick, book, etc). Predict what will
More informationChapter 2 Motion Along a Straight Line
Chapter Motion Along a Straight Line In this chapter we will study how objects moe along a straight line The following parameters will be defined: (1) Displacement () Aerage elocity (3) Aerage speed (4)
More informationChapter 11 Collision Theory
Chapter Collision Theory Introduction. Center o Mass Reerence Frame Consider two particles o masses m and m interacting ia some orce. Figure. Center o Mass o a system o two interacting particles Choose
More informationPhysics 1: Mechanics
Physics 1: Mechanics Đào Ngọc Hạnh Tâm Office: A1.53, Email: dnhtam@hcmiu.edu.n HCMIU, Vietnam National Uniersity Acknowledgment: Most of these slides are supported by Prof. Phan Bao Ngoc credits (3 teaching
More informationDynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates
Dynamic Vehicle Routing with Moing Demands Part II: High speed demands or low arrial rates Stephen L. Smith Shaunak D. Bopardikar Francesco Bullo João P. Hespanha Abstract In the companion paper we introduced
More informationMOTION OF FALLING OBJECTS WITH RESISTANCE
DOING PHYSICS WIH MALAB MECHANICS MOION OF FALLING OBJECS WIH RESISANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECORY FOR MALAB SCRIPS mec_fr_mg_b.m Computation
More informationMotion in Two and Three Dimensions
PH 1-1D Spring 013 Motion in Two and Three Dimensions Lectures 5,6,7 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More information21.60 Worksheet 8 - preparation problems - question 1:
Dynamics 190 1.60 Worksheet 8 - preparation problems - question 1: A particle of mass m moes under the influence of a conseratie central force F (r) =g(r)r where r = xˆx + yŷ + zẑ and r = x + y + z. A.
More informationA Geometric Review of Linear Algebra
A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. I assume the reader is familiar with basic (i.e., high school) algebra and trigonometry.
More informationMagnetic Fields Part 3: Electromagnetic Induction
Magnetic Fields Part 3: Electromagnetic Induction Last modified: 15/12/2017 Contents Links Electromagnetic Induction Induced EMF Induced Current Induction & Magnetic Flux Magnetic Flux Change in Flux Faraday
More informationState-space Modelling of Hysteresis-based Control Schemes
European Control Conference (ECC) July 7-9,, Zürich, Switzerland. State-space Modelling of Hysteresis-based Control Schemes Soumya Kundu Ian A. Hiskens Abstract The paper deelops a state-space model for
More informationUNIVERSITY OF TRENTO ITERATIVE MULTI SCALING-ENHANCED INEXACT NEWTON- METHOD FOR MICROWAVE IMAGING. G. Oliveri, G. Bozza, A. Massa, and M.
UNIVERSITY OF TRENTO DIPARTIMENTO DI INGEGNERIA E SCIENZA DELL INFORMAZIONE 3823 Poo Trento (Italy), Via Sommarie 4 http://www.disi.unitn.it ITERATIVE MULTI SCALING-ENHANCED INEXACT NEWTON- METHOD FOR
More informationThe Kinetic Theory of Gases
978-1-107-1788-3 Classical and Quantum Thermal Physics The Kinetic Theory of Gases CHAPTER 1 1.0 Kinetic Theory, Classical and Quantum Thermodynamics Two important components of the unierse are: the matter
More informationPurpose of the experiment
Impulse and Momentum PES 116 Adanced Physics Lab I Purpose of the experiment Measure a cart s momentum change and compare to the impulse it receies. Compare aerage and peak forces in impulses. To put the
More informationLecture #8-6 Waves and Sound 1. Mechanical Waves We have already considered simple harmonic motion, which is an example of periodic motion in time.
Lecture #8-6 Waes and Sound 1. Mechanical Waes We hae already considered simple harmonic motion, which is an example of periodic motion in time. The position of the body is changing with time as a sinusoidal
More informationE : Ground-penetrating radar (GPR)
Geophysics 3 March 009 E : Ground-penetrating radar (GPR) The EM methods in section D use low frequency signals that trael in the Earth by diffusion. These methods can image resistiity of the Earth on
More informationResearch and Demonstration of the Refraction Problems
American Journal of Applied Mathematics 8; 6(): 87-96 http://www.sciencepublishinggroup.com/j/ajam doi:.648/j.ajam.86.9 ISSN: 33-43 (Print); ISSN: 33-6X (Online) Research and Demonstration of the Refraction
More informationKinematics on oblique axes
Bolina 1 Kinematics on oblique axes arxi:physics/01111951 [physics.ed-ph] 27 No 2001 Oscar Bolina Departamento de Física-Matemática Uniersidade de São Paulo Caixa Postal 66318 São Paulo 05315-970 Brasil
More informationSemi-Empirical Formulas of Drag/Lift Coefficients for High Speed Rigid Body Manoeuvring in Water Column
Semi-Empirical Formulas of Drag/Lift Coefficients for High Speed Rigid Body Manoeuring in Water Column Peter C. Chu and Chenwu Fan Naal Ocean Analysis and Prediction Laboratory Naal Postgraduate School,
More informationA-level Mathematics. MM03 Mark scheme June Version 1.0: Final
-leel Mathematics MM0 Mark scheme 660 June 0 Version.0: Final Mark schemes are prepared by the Lead ssessment Writer and considered, together with the releant questions, by a panel of subject teachers.
More informationCooperative Target Capturing with Multiple Heterogeneous Vehicles
Cooperative Capturing with Multiple Heterogeneous Vehicles Sangeeta Daingade and Arpita Sinha Abstract This paper presents a Nonlinear Cyclic Pursuit based target capturing strategy for a group of autonomous
More informationPHYS 1441 Section 002 Lecture #6
PHYS 1441 Section 00 Lecture #6 Monday, Feb. 4, 008 Examples for 1-Dim kinematic equations Free Fall Motion in Two Dimensions Maximum ranges and heights Today s homework is homework #3, due 9pm, Monday,
More informationNoise constrained least mean absolute third algorithm
Noise constrained least mean absolute third algorithm Sihai GUAN 1 Zhi LI 1 Abstract: he learning speed of an adaptie algorithm can be improed by properly constraining the cost function of the adaptie
More informationPhysics 4A Solutions to Chapter 4 Homework
Physics 4A Solutions to Chapter 4 Homework Chapter 4 Questions: 4, 1, 1 Exercises & Problems: 5, 11, 3, 7, 8, 58, 67, 77, 87, 11 Answers to Questions: Q 4-4 (a) all tie (b) 1 and tie (the rocket is shot
More informationGeostrophy & Thermal wind
Lecture 10 Geostrophy & Thermal wind 10.1 f and β planes These are planes that are tangent to the earth (taken to be spherical) at a point of interest. The z ais is perpendicular to the plane (anti-parallel
More informationReview of Matrices and Vectors 1/45
Reiew of Matrices and Vectors /45 /45 Definition of Vector: A collection of comple or real numbers, generally put in a column [ ] T "! Transpose + + + b a b a b b a a " " " b a b a Definition of Vector
More informationNonlinear Lateral Control of Vision Driven Autonomous Vehicles*
Paper Machine Intelligence & Robotic Control, Vol. 5, No. 3, 87 93 (2003) Nonlinear ateral Control of Vision Drien Autonomous Vehicles* Miguel Ángel Sotelo Abstract: This paper presents the results of
More informationarxiv: v1 [math.pr] 25 Jun 2015
Random walks colliding before getting trapped Louigi Addario-Berry Roberto I. Olieira Yual Peres Perla Sousi ari:156.78451 math.pr 25 Jun 215 Abstract Let P be the transition matrix of a finite, irreducible
More informationPosition in the xy plane y position x position
Robust Control of an Underactuated Surface Vessel with Thruster Dynamics K. Y. Pettersen and O. Egeland Department of Engineering Cybernetics Norwegian Uniersity of Science and Technology N- Trondheim,
More informationMotors and Generators
Physics Motors and Generators New Reised Edition Brian Shadwick Contents Use the table of contents to record your progress through this book. As you complete each topic, write the date completed, then
More informationSpace Probe and Relative Motion of Orbiting Bodies
Space robe and Relatie Motion of Orbiting Bodies Eugene I. Butiko Saint etersburg State Uniersity, Saint etersburg, Russia E-mail: e.butiko@phys.spbu.ru bstract. Seeral possibilities to launch a space
More informationVariational Calculus & Variational Principles In Physics
Variational Calculus & Variational Principles In Physics Michael A. Carchidi Noember 3, 9 1. A Typical Problem In The Calculus Let R be the set of real numbers and let f : R R be a differentiable function
More informationPhysics 4C Spring 2016 Test 3
Physics 4C Spring 016 Test 3 Name: June 1, 016 Please show your work! Answers are not complete without clear reasoning. When asked for an expression, you must gie your answer in terms of the ariables gien
More informationarxiv: v1 [cs.ni] 13 Aug 2013
Robotic Message Ferrying for Wireless Networks using Coarse-Grained Backpressure Control Shangxing Wang Dept. of Electrical Engineering, Uniersity of Southern California, Los Angeles, CA Email: shangxiw@usc.edu
More informationAlgebra Based Physics. Motion in One Dimension. 1D Kinematics Graphing Free Fall 2016.notebook. August 30, Table of Contents: Kinematics
Table of Contents: Kinematics Algebra Based Physics Kinematics in One Dimension 06 03 www.njctl.org Motion in One Dimension Aerage Speed Position and Reference Frame Displacement Aerage Velocity Instantaneous
More informationThe Budgeted Minimum Cost Flow Problem with Unit Upgrading Cost
The Budgeted Minimum Cost Flow Problem with Unit Upgrading Cost Christina Büsing Sarah Kirchner Arie Koster Annika Thome October 6, 2015 Abstract The budgeted minimum cost flow problem (BMCF(K)) with unit
More informationInformation Propagation Speed in Bidirectional Vehicular Delay Tolerant Networks
Information Propagation Speed in Bidirectional Vehicular Delay Tolerant Networks Emmanuel Baccelli, Philippe Jacquet, Bernard Mans, Georgios Rodolakis INRIA, France, {firstname.lastname}@inria.fr Macquarie
More informationMath 425 Lecture 1: Vectors in R 3, R n
Math 425 Lecture 1: Vectors in R 3, R n Motiating Questions, Problems 1. Find the coordinates of a regular tetrahedron with center at the origin and sides of length 1. 2. What is the angle between the
More information(a)!! d = 17 m [W 63 S]!! d opposite. (b)!! d = 79 cm [E 56 N] = 79 cm [W 56 S] (c)!! d = 44 km [S 27 E] = 44 km [N 27 W] metres. 3.
Chapter Reiew, pages 90 95 Knowledge 1. (b). (d) 3. (b) 4. (a) 5. (b) 6. (c) 7. (c) 8. (a) 9. (a) 10. False. A diagram with a scale of 1 cm : 10 cm means that 1 cm on the diagram represents 10 cm in real
More informationLesson 6: Apparent weight, Radial acceleration (sections 4:9-5.2)
Beore we start the new material we will do another Newton s second law problem. A bloc is being pulled by a rope as shown in the picture. The coeicient o static riction is 0.7 and the coeicient o inetic
More informationDesign and Control of Novel Tri-rotor UAV
Design and Control of Noel Tri-rotor UAV Mohamed Kara Mohamed School of Electrical and Electronic Engineering The Uniersity of Manchester Manchester, UK, M 9PL Email: Mohamed.KaraMohamed@postgrad.manchester.ac.uk
More informationDIFFERENTIAL DRAG SPACECRAFT RENDEZVOUS USING AN ADAPTIVE LYAPUNOV CONTROL STRATEGY
(Preprint) IAA-AAS-DyCoSS1-XX-XX DIFFERENIAL DRAG SPACECRAF RENDEZVOUS USING AN ADAPIVE LYAPUNOV CONROL SRAEGY Daid Pérez * and Riccardo Beilacqua I. INRODUCION his paper introduces a noel Lyapuno-based
More informationParameters Identification of Equivalent Circuit Diagrams for Li-Ion Batteries
Parameters Identification of Equialent Circuit Diagrams for Li-Ion eries Ahmad ahmoun, Helmuth Biechl Uniersity of Applied ciences Kempten Ahmad.ahmoun@stud.fh-empten.de, biechl@fh-empten.de Abstract-eries
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More informationv v Downloaded 01/11/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
The pseudo-analytical method: application of pseudo-laplacians to acoustic and acoustic anisotropic wae propagation John T. Etgen* and Serre Brandsberg-Dahl Summary We generalize the pseudo-spectral method
More informationProbabilistic Engineering Design
Probabilistic Engineering Design Chapter One Introduction 1 Introduction This chapter discusses the basics of probabilistic engineering design. Its tutorial-style is designed to allow the reader to quickly
More informationCases of integrability corresponding to the motion of a pendulum on the two-dimensional plane
Cases of integrability corresponding to the motion of a pendulum on the two-dimensional plane MAXIM V. SHAMOLIN Lomonoso Moscow State Uniersity Institute of Mechanics Michurinskii Ae.,, 99 Moscow RUSSIAN
More informationIntroduction to Thermodynamic Cycles Part 1 1 st Law of Thermodynamics and Gas Power Cycles
Introduction to Thermodynamic Cycles Part 1 1 st Law of Thermodynamics and Gas Power Cycles by James Doane, PhD, PE Contents 1.0 Course Oeriew... 4.0 Basic Concepts of Thermodynamics... 4.1 Temperature
More informationUnit 11: Vectors in the Plane
135 Unit 11: Vectors in the Plane Vectors in the Plane The term ector is used to indicate a quantity (such as force or elocity) that has both length and direction. For instance, suppose a particle moes
More informationF = q v B. F = q E + q v B. = q v B F B. F = q vbsinφ. Right Hand Rule. Lorentz. The Magnetic Force. More on Magnetic Force DEMO: 6B-02.
Lorentz = q E + q Right Hand Rule Direction of is perpendicular to plane containing &. If q is positie, has the same sign as x. If q is negatie, has the opposite sign of x. = q = q sinφ is neer parallel
More informationOPTIMIZATION OF FLOWS AND ANALYSIS OF EQUILIBRIA IN TELECOMMUNICATION NETWORKS
OPTIMIZATION OF FLOWS AND ANALYSIS OF EQUILIBRIA IN TELECOMMUNICATION NETWORKS Ciro D Apice (a, Rosanna Manzo (b, Luigi Rarità (c (a, (b, (c Department of Information Engineering and Applied Mathematics,
More informationThree-dimensional Guidance Law for Formation Flight of UAV
ICCAS25 Three-dimensional Guidance aw for Formation Flight of UAV Byoung-Mun Min * and Min-Jea Tahk ** * Department of Aerospace Engineering KAIST Daejeon Korea (Tel : 82-42-869-3789; E-mail: bmmin@fdcl.kaist.ac.kr)
More informationA. unchanged increased B. unchanged unchanged C. increased increased D. increased unchanged
IB PHYSICS Name: DEVIL PHYSICS Period: Date: BADDEST CLASS ON CAMPUS CHAPTER B TEST REVIEW. A rocket is fired ertically. At its highest point, it explodes. Which one of the following describes what happens
More informationReversal in time order of interactive events: Collision of inclined rods
Reersal in time order of interactie eents: Collision of inclined rods Published in The European Journal of Physics Eur. J. Phys. 27 819-824 http://www.iop.org/ej/abstract/0143-0807/27/4/013 Chandru Iyer
More informationLABORATORY VI. ROTATIONAL DYNAMICS
LABORATORY VI. ROTATIONAL DYNAMICS So far this semester, you hae been asked to think of objects as point particles rather than extended bodies. This assumption is useful and sometimes good enough. Howeer,
More informationTheory of Network Communication
Theory of Network Communication Dr. Christian Scheideler Johns Hopkins Uniersity, Sept. 9 2002 Introduction The goal of this lecture is to gie an introduction to the state of the art in the theory of network
More informationMAGNETIC EFFECTS OF CURRENT-3
MAGNETIC EFFECTS OF CURRENT-3 [Motion of a charged particle in Magnetic field] Force On a Charged Particle in Magnetic Field If a particle carrying a positie charge q and moing with elocity enters a magnetic
More informationAlgorithms and Data Structures 2014 Exercises and Solutions Week 14
lgorithms and ata tructures 0 xercises and s Week Linear programming. onsider the following linear program. maximize x y 0 x + y 7 x x 0 y 0 x + 3y Plot the feasible region and identify the optimal solution.
More informationSimulations of Space Probes and their Motions Relative to the Host Orbital Station
Simulations of Space robes and their Motions Relatie to the Host Orbital Station Eugene I. Butiko St. etersburg State Uniersity, St. etersburg, Russia Abstract Launching a space probe from an orbiting
More informationMATHEMATICAL MODELLING AND IDENTIFICATION OF THE FLOW DYNAMICS IN
MATHEMATICAL MOELLING AN IENTIFICATION OF THE FLOW YNAMICS IN MOLTEN GLASS FURNACES Jan Studzinski Systems Research Institute of Polish Academy of Sciences Newelska 6-447 Warsaw, Poland E-mail: studzins@ibspan.waw.pl
More informationPhysics 111. Help sessions meet Sunday, 6:30-7:30 pm in CLIR Wednesday, 8-9 pm in NSC 098/099
ics Announcements day, ember 7, 2007 Ch 2: graphing - elocity s time graphs - acceleration s time graphs motion diagrams - acceleration Free Fall Kinematic Equations Structured Approach to Problem Soling
More informationMotion in Two and Three Dimensions
PH 1-A Fall 014 Motion in Two and Three Dimensions Lectures 4,5 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More informationIntroduction to the Hardy-Littlewood Circle method Rational Points seminar -Tomer Schlank
Introduction to the Hardy-Littlewood Circle method Rational Points seminar -Tomer Schlank by Asaf Katz, Hebrew Uniersity January 5, 200 Some History and General Remarks The circle method was inented by
More informationA nonstandard cubic equation
MATH-Jan-05-0 A nonstandard cubic euation J S Markoitch PO Box West Brattleboro, VT 050 Dated: January, 05 A nonstandard cubic euation is shown to hae an unusually econoical solution, this solution incorporates
More informationGeneral Lorentz Boost Transformations, Acting on Some Important Physical Quantities
General Lorentz Boost Transformations, Acting on Some Important Physical Quantities We are interested in transforming measurements made in a reference frame O into measurements of the same quantities as
More informationThe number of marks is given in brackets [ ] at the end of each question or part question. The total number of marks for this paper is 72.
ADVANCED GCE UNIT 76/ MATHEMATICS (MEI Mechanics MONDAY MAY 7 Additional materials: Answer booklet (8 pages Graph paper MEI Examination Formulae and Tables (MF Morning Time: hour minutes INSTRUCTIONS TO
More informationOptimal Joint Detection and Estimation in Linear Models
Optimal Joint Detection and Estimation in Linear Models Jianshu Chen, Yue Zhao, Andrea Goldsmith, and H. Vincent Poor Abstract The problem of optimal joint detection and estimation in linear models with
More informationChapter 14 Thermal Physics: A Microscopic View
Chapter 14 Thermal Physics: A Microscopic View The main focus of this chapter is the application of some of the basic principles we learned earlier to thermal physics. This will gie us some important insights
More informationOblivious Routing in Wireless networks. Costas Busch
Obliious Routing in Wireless networks Costas Busch Rensselaer Polytechnic Institute Joint work with: Malik Magdon-Ismail and Jing Xi 1 Outline of Presentation Introduction Network Model Obliious Algorithm
More informationP2.8 THE MECHANICS OF FALLING HAILSTONES AND HAILSWATHS. Kevin Vermeesch* and Ernest Agee Purdue University, West Lafayette, IN. where b is given as
P.8 THE MECHANICS OF FALLING HAILSTONES AND HAILSWATHS Kein Vermeesch* and Ernest Agee Purdue Uniersity, West Lafayette, IN. INTRODUCTION * This study has focused on the mechanics of falling hailstones
More informationTHE FIFTH DIMENSION EQUATIONS
JP Journal of Mathematical Sciences Volume 7 Issues 1 & 013 Pages 41-46 013 Ishaan Publishing House This paper is aailable online at http://www.iphsci.com THE FIFTH DIMENSION EQUATIONS Niittytie 1B16 03100
More informationDYNAMICS. Kinematics of Particles VECTOR MECHANICS FOR ENGINEERS: Tenth Edition CHAPTER
Tenth E CHAPTER 11 VECTOR MECHANICS FOR ENGINEERS: DYNAMICS Ferdinand P. Beer E. Russell Johnston, Jr. Phillip J. Cornwell Lecture Notes: Brian P. Self California Polytechnic State Uniersity Kinematics
More informationTools for Investigation of Dynamics of DC-DC Converters within Matlab/Simulink
Tools for Inestigation of Dynamics of DD onerters within Matlab/Simulink Riga Technical Uniersity, Riga, Latia Email: pikulin03@inbox.l Dmitry Pikulin Abstract: In this paper the study of complex phenomenon
More informationPhysics 107 HOMEWORK ASSIGNMENT #9b
Physics 07 HOMEWORK SSIGNMENT #9b Cutnell & Johnson, 7 th edition Chapter : Problems 5, 58, 66, 67, 00 5 Concept Simulation. reiews the concept that plays the central role in this problem. (a) The olume
More informationSpring 2000 HIGHER STILL. Physics. Student Materials Advanced Higher. Summary Notes Unit 3 Wave Phenomena. Physics (AH): Mechanics - Student Materials
Spring 2000 HIGHER STILL Physics Student Materials Adanced Higher Summary Notes Unit 3 Wae Phenomena Physics (AH): Mechanics - Student Materials WAVE PHENOMENA The Content Statements for this unit are
More informationSIMULATIONS OF CHARACTERISTICS OF TUNED LIQUID COLUMN DAMPER USING AN ELLIPTICAL FLOW PATH ESTIMATION METHOD
October -7, 008, Beijing, China SIMULATIONS OF CHARACTERISTICS OF TUNED LIQUID COLUMN DAMPER USING AN ELLIPTICAL FLOW PATH ESTIMATION METHOD P. Chaiiriyawong, S. Limkatanyu and T. Pinkaew 3 Lecturer, Dept.
More informationLinear Momentum and Collisions Conservation of linear momentum
Unit 4 Linear omentum and Collisions 4.. Conseration of linear momentum 4. Collisions 4.3 Impulse 4.4 Coefficient of restitution (e) 4.. Conseration of linear momentum m m u u m = u = u m Before Collision
More informationTo Feel a Force Chapter 13
o Feel a Force Chapter 13 Chapter 13:. A. Consequences of Newton's 2nd Law for a human being he forces acting on a human being are balanced in most situations, resulting in a net external force of zero.
More informationVelocity, Acceleration and Equations of Motion in the Elliptical Coordinate System
Aailable online at www.scholarsresearchlibrary.com Archies of Physics Research, 018, 9 (): 10-16 (http://scholarsresearchlibrary.com/archie.html) ISSN 0976-0970 CODEN (USA): APRRC7 Velocity, Acceleration
More informationLecture 12! Center of mass! Uniform circular motion!
Lecture 1 Center of mass Uniform circular motion Today s Topics: Center of mass Uniform circular motion Centripetal acceleration and force Banked cures Define the center of mass The center of mass is a
More informationEE243 Advanced Electromagnetic Theory Lec # 15 Plasmons. Reading: These PPT notes and Harrington on corrugated surface 4.8.
EE Applied EM Fall 6, Neureuther Lecture #5 Ver //6 EE43 Adanced Electromagnetic Theory Lec # 5 Plasmons Fixed Slides 6, (Group and Phase Velocity for Guided Waes) (Sign Reersals for Propagation Direction
More informationMA 266 Review Topics - Exam # 2 (updated)
MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential
More informationCIRCULAR MOTION EXERCISE 1 1. d = rate of change of angle
CICULA MOTION EXECISE. d = rate of change of angle as they both complete angle in same time.. c m mg N r m N mg r Since r A r B N A N B. a Force is always perpendicular to displacement work done = 0 4.
More information