Path planning with PH G2 splines in R2
|
|
- Giles Kennedy
- 5 years ago
- Views:
Transcription
1 Path panning with PH G2 spines in R2 Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru To cite this version: Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru. Path panning with PH G2 spines in R2. IEEE 1st Internationa Conference on Systems and Computer Science, Aug 2012, Lie, France <ha > HAL Id: ha Submitted on 17 Jan 2013 HAL is a muti-discipinary open access archive for the deposit and dissemination of scientific research documents, whether they are pubished or not. The documents may come from teaching and research institutions in France or abroad, or from pubic or private research centers. L archive ouverte puridiscipinaire HAL, est destinée au dépôt et à a diffusion de documents scientifiques de niveau recherche, pubiés ou non, émanant des étabissements d enseignement et de recherche français ou étrangers, des aboratoires pubics ou privés.
2 Path panning with PH G 2 spines in R 2 Laurent Gajny, Richard Béarée, Éric Nyiri, Oivier Gibaru Arts et Métiers ParisTech, LSIS - UMR CNRS 7296, 8 Bouevard Louis XIV, Lie Cedex, France INRIA Lie-Nord-Europe, 40 avenue Haey Vieneuve d Ascq, France Abstract In this artice, we justify the use of parametric panar Pythagorean Hodograph spine curves in path panning. The eegant properties of such spines enabe us to design an efficient interpoator agorithm, more precise than the cassica Tayor interpoators and faster than an interpoator based on arc ength computations. Index Terms Path panning, Pythagorean-hodograph, Spines. I. INTRODUCTION Designing a good geometric support from given points is a fundamenta issue in path panning for industria machines Machine-too, manipuator, robot... Traditionnay, the methods used consist in minimizing the L p norm of the second derivative of poynomia sines which interpoate the given points, for p N See for exampe [1] for the genera method and [9] for its appication. In this artice, we are interested in a subcass of poynomia spines caed Pythagorean Hodograph spines. Pythagorean Hodograph curves PH curves, introduced in 1990 in [2] are characterized by the specia property that the derivative of their arc ength is a poynomia or rationa function of the curve parameter. This property enabes to compute exacty usefu quantities such as curvature, bending energy or arc ength. For simpe poynomia curve, we can not, in genera, cacuate these quantities and so, we appy some quadrature rues. By a good choice of the quadrature rue, we can reach a satisfying precision but it may be highy time consuming. It is not difficut to design a geometric support which interpoates data points using PH curves. We can appy for exampe L 1 C 1 method from [1] which gives some derivative vaues at each data point and reconstruct a soution as PH quintic spine by C 1 Hermite interpoation See [3]. But C 1 continuity is not sufficient for path panning of mechanica system if trajectory foowing performances are of interest. Indeed, curvature discontinuities, which are not physicay foowed, excite the structure and can significanty deteriorate the accuracy. Thus, we use the G 2 interpoation scheme introduced by Jakič et a. in [7]. By this method, we can ony treat, for the moment, convex data but it uses ow degree spine cubic. The existence of soution is guaranteed and the iterative method to construct the soution converges rapidy. Some artices See for exampe [4], [5], [11] dea with the use of PH curves in CNC. In these artices, authors ony treat the case of a singe PH curve and designed an interpoator based on arc ength computation in order to foow the parametric geometric support. We wi compare the performances in running time and accuracy of such an interpoator for both PH and cassica poynomia curve cases. Then, We wi extend this procedure on a PH spine. This paper is organized as foows. Section II gives some basic information about Pythagorean Hodograph cubic curves. In section III, we present briefy the interpoation scheme that we wi use in our experiment. Section IV is dedicated to the PH interpoator, its appication and a comparison with other interpoators. Some concusions wi be drawn in the ast section. II. PYTHAGOREAN HODOGRAPH CUBIC CURVES A poynomia Pythagorean hodograph curve rτ = xτ, yτ, τ [0, 1] is defined by the property that its hodograph r τ = x τ, y τ satisfies the Pythagorean condition : x 2 τ + y 2 τ = σ 2 τ, 1 for some poynomia στ. This remarkabe property enabes to compute exacty some usefu quantities ike arc ength, curvature or bending energy. To define curves satisfying 1, we use two poynomias, uτ and vτ such that : x τ = u 2 τ v 2 τ, y τ = 2uτvτ, στ = u 2 τ + v 2 τ. For PH cubics, we sha choose uτ and vτ inear functions expressed in the Bernstein basis : uτ = u 0 B 1 0τ + u 1 B 1 1τ, vτ = v 0 B 1 0τ + v 1 B 1 1τ. In order to define reguar curves, we may assume u 1 u v 1 v and u 0 v 1 u 1 v 0 0. By integrating the hodograph, we obtain the formuation of Bézier contro points of a PH cubic curve. Q 1 = Q u2 0 v 2 0, 2u 0 v 0, Q 2 = Q u 0u 1 v 0 v 1, u 0 v 1 + u 1 v 0, Q 3 = Q u2 1 v 2 1, 2u 1 v
3 We have a sufficient and necessary geometric condition on the contro poygon, iustrated in Fig. 1, for a parametric cubic curve to be a PH cubic. The foowing theorem wi be very usefu in the spine interpoation scheme we wi use in section III. Theorem 1: Let rτ be a parametric panar cubic curve with Bézier contro points Q 0, Q 1, Q 2, Q 3. Let L i = Q i Q i+1, i {0, 1, 2} be the engths of the contro poygon edges and θ 1, θ 2 be the contro poygon anges at the interior vertices Q 1 and Q 2. Then rτ is a PH cubic if and ony if : L 1 = L 0 L 2, and, θ 1 = θ 2. 4 The proof of this theorem can be found in [6] section d 0, P 0, P 1, d 1 a convex G 1 Hermite data. P 0 and P 1 are two points in the pane, d 0 and d 1 the associated tangent vectors. We are ooking for λ 0, λ 1 > 0 such that : Q 0 = P 0, Q 1 = P 0 + λ 0 d 0, Q 2 = P 1 λ 1 d 1, Q 3 = P 1, are the contro points of a PH cubic. The authors determine a unique admissibe soution : λ 0 = Λ 0 d 0, P 0, d 1, λ 1 = Λ 1 d 0, P 0, d 1, under the assumption that the anges ϕ 0 = d 0, P 0 and ϕ 1 = d 1, P 0 satisfy the condition : 5 Fig. 1. A PH cubic curve and its contro poygon. ϕ 0 + ϕ 1 < 4 3 π. The goba agorithm of G 2 PH cubic spine determination can be written as foows : Inputs : Data points P 0, P 1,..., P n, end conditions d 0, d n, a toerance ε. Outputs : Tangent vectors at each data points and engths of the first and the third edges of each contro poygon. III. PH G 2 SPLINE INTERPOLATION The probem is to interpoate P 0, P 1,..., P n in R 2 by a G 2 PH cubic spine with end conditions d 0 and d n, respectivey the tangent vectors at P 0 and P n. We wi use L 1 C 1 interpoation, which have good shape-preserving properties, to compute d 0 and d n. We give the main resut of the artice of Jakič et a. [7]. Theorem 2: Let P 0, P 1,..., P n a set of convex data and d 0, d n the tangent vectors at P 0 and P n. The interpoation probem of finding a G 2 PH cubic spine passing through the P i with respect of d 0 and d n has an admissibe without oops or cusps soution if and ony if the anges ϕ 0 = d 1, P 1, ϕ = P 1, P, = 1, 2,..., n 1, ϕ n = P n 1, d n, satisfy ϕ i +ϕ i 1 < M with M = 4π/3 for i = 1, 2,..., n 1. If M = Kπ with, K = π arccos < 4/3, 3 then the soution is unique. The system of non-inear equations, directy obtained by expicitation of 4 on each segment is soved by an iterative method. To understand the agorithm, we may consider the G 1 Hermite interpoation probem with PH cubic. Let 1 λ = 1, 1,..., 1 R 2m, r = 0 2 Compute d [0] = P +1 P 1 P +1 P 1, = 1, 2,..., m 1. 3 Compute λ 2 = Λ 0 d [r], P, d [r] λ 2+1 = Λ 1 d [r] = 0, 1,..., m 1. 4 Compute d [r+1] = 1 λ 2 2, P, d [r] +1, +1, P λ 2+1 d [r] λ 2 1 P 1 λ 2 2 d [r] 1, = 1, 2,..., m 1. 5 D = d [r+1], = 0, 1,..., m. 6 If λ od λ new < ε and D od D new < ε, STOP ese r = r + 1 and go back to 3. We have appied this agorithm to interpoate the ogarithmic spira and a circe. We note in Fig. 2 that the resuting curves preserve we the shape of the data. A. Principe IV. PH INTERPOLATOR In this part, we wish that a mechanica system foows a geometric motion support with respect of a known feedrate. The geometric support is designed by interpoation of some checkpoints. We choose here a set of convex data that we interpoate by the previous procedure. Knowing the feedrate enabes us, by integration, to define a position aw See
4 a b Fig. 2. Circe. Interpoation by G 2 PH cubic spine. a Logarithmic spira b Fig. 3. A simpe feedrate top and its associated position aw bottom. Fig. 3 which means that we know, at time t 0, the distance the system must have run through. The main probem is to determine the geometric parameter which corresponds to the time parameter. In practice, we discretize the time interva [0, T ] with a constant 1 samping period t. We obtain the ist 0 = t 0 < t 1 < < t m = T. For each i {0, 1,..., m}, we have to find the associated geometric parameter of the curve τ i by soving the foowing equation : t i = sτ i, where and s are respectivey the arc ength functions of the position aw and the spine. B. Arc ength computation for a PH curve Let rτ be a PH cubic defined by the coefficients u 0, u 1, v 0, v 1 as in 2. Arc ength at parameter τ [0, 1] may be expressed in the form : where and : sτ = 3 n s k 1 τ 3 k τ k, 6 k k=0 k 1 s 0 = 0, s k = σ j, k = 1,..., 3 j=0 σ 0 = u v 2 0, σ 1 = u 0 u 1 + v 0 v 1, σ 2 = u v We can choose a variabe stepsize. It wi not change the performances of the method. We have in particuar : s1 = σ 0 + σ 1 + σ So contrary to ordinary poynomia curve, we have a cosed form for the arc ength at parameter τ. C. Agorithm We now describe the agorithm we wi use further both in case of a simpe poynomia spine and a PH spine. Inputs : Data points and associated parameters, the interpoating spine defined by some coefficients, the position aw, a time ist. Output : A ist of geometric parameters. 1 Compute k = t k, k = 0, 1,... m. 2 Compute S j, j = 0, 1,..., n the arc engths of the PH cubic spine at each ν j. 3 τ 0 = ν 1 ; r 0 = 1. 4 Compute r k, k = 1, 2,..., n the number of the spine piece where τ k must be. 5 If r k = r k 1, sove : τk ν rk k k 1 = s rk τ k 1 τ k ν rk +1. ν rk +1 ν rk s rk τk 1 ν rk ν rk +1 ν rk,
5 6 If r k r k 1, sove : τk ν rk k S rk = s rk, ν rk +1 ν rk ν rk τ k ν rk +1. Here, s j. is the arc ength function of the j th piece of the spine. The equations to sove in case of PH cubic spine are aso cubics. So, we can sove them exacty by Cardano method. In practice, for a given precision ε = 10 13, we use the Newton-Raphson agorithm which eads to the soution with machine precision in two or three iterations. For poynomia spine, using this agorithm eads to use a quadrature formua. We wi see that in appications, it is more time consuming. V. NUMERICAL EXPERIMENTS A. Traveing a singe PH curve Let us give a PH curve and a feedrate as shown in Fig. 4. This feedrate has been generated with the imited jerk strategy described in [10]. We wi appy the agorithm presented in section IV-C for this curve using the Pythagorean property 1 and in a second experiment, seeing it as a cassica poynomia curve and using a quadrature rue. In both case, we choose a samping period of 10 3 s. Fig. 5. Semiogarithmic errors for position and speed in the singe PH case. Fig. 4. A singe PH curve to trave Top, eft. The desired feedrate Top, right and jerk Bottom. An error anaysis for the PH case is shown in Fig. 5. Fina position error is about m and speed error gobay osciates between and m.s 1. It occurs very sma osciations on the feedrate. We now consider that we do not have a PH curve but ony a genera poynomia one. We obtain simiar accuracy as shown in Fig. 6. The major difference is the running time. For this exampe, treating about 1500 points of discretization takes near to 85 seconds in the poynomia case and ony 6 in the PH case. In order to prove this observation on running time, we aunch severa experiments with different samping periods. The resuts are summarized in Tabe I and Fig. 7. We notice a quasi-inearity property with approximated constant Fig. 6. Semiogarithmic errors for position and speed in the singe poynomia case. in PH case and in poynomia case. So, for the given precision ε, it takes approximativey 14.4 more time to run the interpoator with poynomia spines than PH spines. In poynomia case, running time wi be too high when we wi have to treat an entire spine and not just a singe part. That seems to be the reason that such an interpoator is not used and
6 that interpoators based on Tayor expansions, ess accurate but ess time consuming, are prefered. Samping Number CPU Time, CPU Time, period of points PH poynomia 10 1 s s s s s s s s s s s s s s s 10 2 s s s s s s s s s s s s s s s 10 3 s s s s s s s s s s s s s s s 10 4 s s s TABLE I CPU TIMES OF THE INTERPOLATOR IN SINGLE PH AND SINGLE POLYNOMIAL CASES FOR SOME SAMPLING PERIODS. a Fig. 8. b a PH cubic spine to trave. b The desired feedrate. Once again, we present error graphs in Fig. 9. We ceary distinguish the pieces of the spine in these graphs. When we change of spine piece, the previous error is corrected by starting on the new piece by its extremity and not by the ast computed point. This produces imited jumps for the speed error which are not important. a b Fig. 7. CPU time in function of the number of discretization points, a in the PH case, b in the poynomia case. B. Traveing a PH spine We are now interested in foowing a rea PH G 2 spine with a desired feedrate as shown in Fig. 8. Again, this feedrate has been generated with the strategy described in [10]. Fig. 9. Semiogarithmic errors for position and speed in the PH spine case.
7 We compare the efficiency of the agorithm with an interpoator based on a first order Tayor expansion. The parameters τ i, i = 0, 1,..., m with m = 7959 in our exampe are given by the approximation : τ i+1 = τ i + V i dγτ dτ T s, 9 τ=τi where V i is the desired speed at time t i, Γ the spine obtained by the L 2 1C 2 method see [1] and T s the samping period. With this interpoator, we evidence in Fig. 10 an average position error of 0.3 mm and a fina position error of 5 microns. Error graphs presented in Fig. 11 show the accuracy defaut of the method. Fig. 11. Semiogarithmic errors for position and speed with the first order Tayor interpoator. Fig. 10. Goba resut and fina error position for the first order Tayor interpoator. CONCLUSION In this artice, we have justified the use of Pythagorean hodograph spines in path panning. The associated PH interpoator based on exact arc ength cacuation shows simiar accuracy than a simiar one for simpe poynomia spines but it is amost fifteen times faster. Moreover, it is much more precise than a Tayor interpoator. Thus, PH spines seems to be a good soution for rea time and accuracy issues. REFERENCES [1] Phiippe Auquiert, Oivier Gibaru, Éric Nyiri, Fast L k 1 Ck poynomia spine interpoation agorithm with shape-preserving properties. Computer Aided Geometric Design [2] Rida T. Farouki, T. Sakkais, Pythagorean Hodographs. IBM Journa of Research and Deveopment [3] Rida T. Farouki, Andrew Neff, Hermite interpoation by Pythagorean hodograph quintics. Mathematics of computation 64, n [4] Rida T. Farouki, Sagar Shah, Rea time CNC interpoators for Pythagorean-hodograph curves. Computer Aided Geometric Design [5] Rida T. Farouki, Yi-Feng Tsai, Exact Tayor series coefficient for variabe-feedrate CNC curve interpoators. Computer Aided Geometric Design [6] Rida T. Farouki, Pythagorean-hodograph curves - Agebra and geometry inseparabe. Springer-Verag [7] Gašper Jakič, Jernej Kozak, Marjeta Krajnc, Vito Vitrih & Emi Žagar, On interpoation by panar cubic G 2 Pythagorean-Hodograph spine curves. Mathematics of computation 79, n [8] D.S. Meek, D.J. Waton, G 2 curve design with a pair of Pythagorean Hodograph quintic spira segments. Computer Aided Geometric Design [9] Ade Oabi, Richard Béarée, Éric Nyiri, Oivier Gibaru Enhanced Trajectory Panning For Machining With Industria Six-Axis Robots.IEEE Internationa Conference on Industria Technoogy, ICIT [10] Ade Oabi, Richard Béarée, Oivier Gibaru, Mohamed Damak, Feedrate panning for machining with industria six-axis robots. Contro Engineering Practice, Vo. 18 Issue [11] Yi-Feng Tsai, Rida T. Farouki, Bryan Fedman, Performance anaysis of CNC interpoators for time-dependent feedrates aong PH curves. Computer Aided Geometric Design Future work wi concentrate on appication of this work for rea path and on designing a G 2 PH interpoation scheme for non convex data.
List edge and list total colorings of planar graphs without non-induced 7-cycles
List edge and ist tota coorings of panar graphs without non-induced 7-cyces Aijun Dong, Guizhen Liu, Guojun Li To cite this version: Aijun Dong, Guizhen Liu, Guojun Li. List edge and ist tota coorings
More informationTHINKING IN PYRAMIDS
ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationPolar Snakes: a fast and robust parametric active contour model
Poar Snakes: a fast and robust parametric active contour mode Christophe Coewet To cite this version: Christophe Coewet. Poar Snakes: a fast and robust parametric active contour mode. IEEE Int. Conf. on
More informationTHE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE
THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on
More informationMultiway Regularized Generalized Canonical Correlation Analysis
Mutiway Reguarized Generaized Canonica Correation Anaysis Arthur Tenenhaus, Laurent Le Brusquet, Gisea Lechuga To cite this version: Arthur Tenenhaus, Laurent Le Brusquet, Gisea Lechuga. Mutiway Reguarized
More informationEstablishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
Georgia Southern University Digita Commons@Georgia Southern Mathematica Sciences Facuty Pubications Department of Mathematica Sciences 2009 Estabishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationGeometric Lagrange Interpolation by Planar Cubic Pythagorean-hodograph Curves
Geometric Lagrange Interpolation by Planar Cubic Pythagorean-hodograph Curves Gašper Jaklič a,c, Jernej Kozak a,b, Marjeta Krajnc b, Vito Vitrih c, Emil Žagar a,b, a FMF, University of Ljubljana, Jadranska
More informationCase report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122,
Case report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122, 244902 Juan Olives, Zoubida Hammadi, Roger Morin, Laurent Lapena To cite this version: Juan Olives,
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationA Solution to the 4-bit Parity Problem with a Single Quaternary Neuron
Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria
More informationThe Group Structure on a Smooth Tropical Cubic
The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,
More informationarxiv: v1 [math.ca] 6 Mar 2017
Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationHigh-order approximations to the Mie series for electromagnetic scattering in three dimensions
Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions
More informationThe EM Algorithm applied to determining new limit points of Mahler measures
Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationCollective organization in an adaptative mixture of experts
Coective organization in an adaptative mixture of experts Vincent Vigneron, Christine Fuchen, Jean-Marc Martinez To cite this version: Vincent Vigneron, Christine Fuchen, Jean-Marc Martinez. Coective organization
More informationMultilayer Kerceptron
Mutiayer Kerceptron Zotán Szabó, András Lőrincz Department of Information Systems, Facuty of Informatics Eötvös Loránd University Pázmány Péter sétány 1/C H-1117, Budapest, Hungary e-mai: szzoi@csetehu,
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationTHE DIRECT KINEMATICS OF REDUNDANT PARALLEL ROBOT FOR PREDICTIVE CONTROL
HE DIREC KINEMAICS OF REDUNDAN PARALLEL ROBO FOR PREDICIVE CONROL BELDA KVĚOSLAV, BÖHM JOSEF, VALÁŠEK MICHAEL Department of Adaptive Systems, Institute of Information heory and Automation, Academy of Sciences
More informationHook lengths and shifted parts of partitions
Hook lengths and shifted parts of partitions Guo-Niu Han To cite this version: Guo-Niu Han Hook lengths and shifted parts of partitions The Ramanujan Journal, 009, 9 p HAL Id: hal-00395690
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationTheory of Generalized k-difference Operator and Its Application in Number Theory
Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication
More informationA proposed nonparametric mixture density estimation using B-spline functions
A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),
More informationData Search Algorithms based on Quantum Walk
Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but
More informationChapter 4. Moving Observer Method. 4.1 Overview. 4.2 Theory
Chapter 4 Moving Observer Method 4.1 Overview For a compete description of traffic stream modeing, one woud reuire fow, speed, and density. Obtaining these parameters simutaneousy is a difficut task if
More informationOn the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems
Internationa Journa of Bifurcation and Chaos Vo. 25 No. 10 2015 1550131 10 pages c Word Scientific Pubishing Company DOI: 10.112/S02181271550131X On the Number of Limit Cyces for Discontinuous Generaized
More informationMethylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma.
Methylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma. Loïc De Pontual, Delphine Trochet, Franck Bourdeaut, Sophie Thomas, Heather Etchevers, Agnes Chompret, Véronique Minard,
More informationUI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE
UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,
More informationConvergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems
Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,
More informationIndirect Optimal Control of Dynamical Systems
Computationa Mathematics and Mathematica Physics, Vo. 44, No. 3, 24, pp. 48 439. Transated from Zhurna Vychisite noi Matematiki i Matematicheskoi Fiziki, Vo. 44, No. 3, 24, pp. 444 466. Origina Russian
More informationPhysics 505 Fall Homework Assignment #4 Solutions
Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationOn Newton-Raphson iteration for multiplicative inverses modulo prime powers
On Newton-Raphson iteration for multiplicative inverses modulo prime powers Jean-Guillaume Dumas To cite this version: Jean-Guillaume Dumas. On Newton-Raphson iteration for multiplicative inverses modulo
More informationDisturbance decoupling by measurement feedback
Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,
More informationA simple kinetic equation of swarm formation: blow up and global existence
A simple kinetic equation of swarm formation: blow up and global existence Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot To cite this version: Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot.
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationFull-order observers for linear systems with unknown inputs
Full-order observers for linear systems with unknown inputs Mohamed Darouach, Michel Zasadzinski, Shi Jie Xu To cite this version: Mohamed Darouach, Michel Zasadzinski, Shi Jie Xu. Full-order observers
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationSome Applications on Generalized Hypergeometric and Confluent Hypergeometric Functions
Internationa Journa of Mathematica Anaysis and Appications 0; 5(): 4-34 http://www.aascit.org/journa/ijmaa ISSN: 375-397 Some Appications on Generaized Hypergeometric and Confuent Hypergeometric Functions
More informationMelodic contour estimation with B-spline models using a MDL criterion
Meodic contour estimation with B-spine modes using a MDL criterion Damien Loive, Ney Barbot, Oivier Boeffard IRISA / University of Rennes 1 - ENSSAT 6 rue de Kerampont, B.P. 80518, F-305 Lannion Cedex
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationFOURIER SERIES ON ANY INTERVAL
FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series
More informationImproving the Reliability of a Series-Parallel System Using Modified Weibull Distribution
Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationBP neural network-based sports performance prediction model applied research
Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationResearch of Data Fusion Method of Multi-Sensor Based on Correlation Coefficient of Confidence Distance
Send Orders for Reprints to reprints@benthamscience.ae 340 The Open Cybernetics & Systemics Journa, 015, 9, 340-344 Open Access Research of Data Fusion Method of Muti-Sensor Based on Correation Coefficient
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationA Novel Learning Method for Elman Neural Network Using Local Search
Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama
More informationThomas Lugand. To cite this version: HAL Id: tel
Contribution à la Modélisation et à l Optimisation de la Machine Asynchrone Double Alimentation pour des Applications Hydrauliques de Pompage Turbinage Thomas Lugand To cite this version: Thomas Lugand.
More informationSTA 216 Project: Spline Approach to Discrete Survival Analysis
: Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing
More informationMaintenance activities planning and grouping for complex structure systems
Maintenance activities panning and grouping for compex structure systems Hai Canh u, Phuc Do an, Anne Barros, Christophe Bérenguer To cite this version: Hai Canh u, Phuc Do an, Anne Barros, Christophe
More informationStatistical Learning Theory: A Primer
Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationMathematical Scheme Comparing of. the Three-Level Economical Systems
Appied Mathematica Sciences, Vo. 11, 2017, no. 15, 703-709 IKAI td, www.m-hikari.com https://doi.org/10.12988/ams.2017.7252 Mathematica Scheme Comparing of the Three-eve Economica Systems S.M. Brykaov
More informationReliability: Theory & Applications No.3, September 2006
REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationAvailable online at ScienceDirect. IFAC PapersOnLine 50-1 (2017)
Avaiabe onine at www.sciencedirect.com ScienceDirect IFAC PapersOnLine 50-1 (2017 3412 3417 Stabiization of discrete-time switched inear systems: Lyapunov-Metzer inequaities versus S-procedure characterizations
More informationConsistent linguistic fuzzy preference relation with multi-granular uncertain linguistic information for solving decision making problems
Consistent inguistic fuzzy preference reation with muti-granuar uncertain inguistic information for soving decision making probems Siti mnah Binti Mohd Ridzuan, and Daud Mohamad Citation: IP Conference
More informationSmart Bolometer: Toward Monolithic Bolometer with Smart Functions
Smart Bolometer: Toward Monolithic Bolometer with Smart Functions Matthieu Denoual, Gilles Allègre, Patrick Attia, Olivier De Sagazan To cite this version: Matthieu Denoual, Gilles Allègre, Patrick Attia,
More informationThe FLRW cosmological model revisited: relation of the local time with th e local curvature and consequences on the Heisenberg uncertainty principle
The FLRW cosmological model revisited: relation of the local time with th e local curvature and consequences on the Heisenberg uncertainty principle Nathalie Olivi-Tran, Paul M Gauthier To cite this version:
More informationAnalysis of Boyer and Moore s MJRTY algorithm
Analysis of Boyer and Moore s MJRTY algorithm Laurent Alonso, Edward M. Reingold To cite this version: Laurent Alonso, Edward M. Reingold. Analysis of Boyer and Moore s MJRTY algorithm. Information Processing
More informationAppendix A: MATLAB commands for neural networks
Appendix A: MATLAB commands for neura networks 132 Appendix A: MATLAB commands for neura networks p=importdata('pn.xs'); t=importdata('tn.xs'); [pn,meanp,stdp,tn,meant,stdt]=prestd(p,t); for m=1:10 net=newff(minmax(pn),[m,1],{'tansig','purein'},'trainm');
More informationTraffic data collection
Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationLegendre Polynomials - Lecture 8
Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha
More informationEaster bracelets for years
Easter bracelets for 5700000 years Denis Roegel To cite this version: Denis Roegel. Easter bracelets for 5700000 years. [Research Report] 2014. HAL Id: hal-01009457 https://hal.inria.fr/hal-01009457
More informationTorsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg
Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationLecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling
Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest
More informationRelated Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage
Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is
More informationA new simple recursive algorithm for finding prime numbers using Rosser s theorem
A new simple recursive algorithm for finding prime numbers using Rosser s theorem Rédoane Daoudi To cite this version: Rédoane Daoudi. A new simple recursive algorithm for finding prime numbers using Rosser
More informationDiscrete Bernoulli s Formula and its Applications Arising from Generalized Difference Operator
Int. Journa of Math. Anaysis, Vo. 7, 2013, no. 5, 229-240 Discrete Bernoui s Formua and its Appications Arising from Generaized Difference Operator G. Britto Antony Xavier 1 Department of Mathematics,
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More information