One Class of Splitting Iterative Schemes
|
|
- Adelia Robinson
- 5 years ago
- Views:
Transcription
1 One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi of a new cla of iterative method for elliptic problem. Thee cheme are baed on a general plitting method, which decompoe a multidimenional parabolic problem into a ytem of one dimenional implicit problem. We ue a pectral tability analyi and invetigate the convergence order of two iterative cheme. Finally, ome reult of numerical experiment are preented. 1 Introduction Iterative method for olving dicrete elliptic problem can be viewed a finite difference cheme for non tationary parabolic problem. The mot important difference i that for elliptic problem we elect the time tep parameter according the requirement of the convergence to the tationary olution and can ignore the approximation error. For plitting method for elliptic problem we refer to [4, 6]. In particular, [5, 6] involve an alternating direction method, [4] preent factorization cheme. The convergence rate can be increaed if optimal nontationary parameter are ued for the definition of each iteration, ee e.g. [6]. Recently the multicomponent verion of alternating direction method were propoed in [1]. Thee cheme are alo ued for olving multidimenional elliptic problem [2]. While in the previou paper the tability of thee plitting cheme wa invetigated by the energy method, we will ue the pectral tability analyi. For ymmetric problem it give neceary and ufficient convergence condition and enable u to find optimal value of iterative parameter. Such analyi wa ued alo in [3] The content of thi paper i organized a follow. In Section 2 we formulate the multicomponent iterative cheme. The convergence of 2D cheme i invetigated in Section 3. In ection 4 we invetigate the tability of 3D iterative cheme, the analyi i done uing numerical experiment. Finally, in Section 5 we tudy the convergence of the p-dimenional multicomponent iterative cheme and prove the energy tability etimate. In Section 6, the Seidel type cheme i formulated and invetigated. P.M.A. Sloot et al. Ed.: ICCS 2002, LNCS 2330, pp , Springer-Verlag Berlin Heidelberg 2002
2 One Cla of Splitting Iterative Scheme Multicomponent iterative cheme Let the trictly elliptic problem be given by i=1 2 u x 2 i = fx, x Q, ux =0 x Q, where Q =0, 1 p. We introduce a uniform grid in Q and approximate the elliptic problem by the following finite difference cheme A α y = f, 1 here A α denote the approximation of a differential operator uing tandard central difference formula. Let, and denote the inner product and the L 2 norm of dicrete function, repectively. Generally, we introduce the following aumption: A α = A α, 0 <m α y A α y, y M α y. It i well known that for the Laplace operator thee aumption are atified with 8 m α, M α 4 h 2, 2 where h i the patial tep ize. In order to olve the ytem of linear equation 1 we define p unknown function y α, α =1, 2,...,p. Then the multicomponent alternating direction MAD cheme i given a ee, [1]: y α ỹ + pa α y α y α + ỹ = 1 p y α, A β yβ = f,, 2,...,p, 3 where y α i the -th iteration of y α. We note that all equation 3 can be olved in parallel. We will invetigate the convergence of the MAD cheme by uing the following error norm: e = ỹ y, r = A αyα f.
3 v. 424 R. Ciegi and V. Pakalnyte 3 Spectral tability analyi of 2D iterative cheme In thi ection we conider a two-dimenional iterative cheme 1 y α ỹ +2A α y α y α + 2 A β yβ = f,, 2. 4 Let denote the error function e α = y α y, α =1, 2,...,p, ẽ = ỹ y. Then the error function atify the following MAD cheme: e α ẽ +2A α e α e α + ẽ = e α, 2 A β eβ =0, α =1, 2, 5 To apply the dicrete von Neumann tability criteria to problem 5, we let e α = N 1 j=1 N 1 k=1 dα,jk injπx 1 inkπx 2, α =1, 2, 6 where N i the number of grid point in one-dimenional grid. It i well known, that injπx α are eigenvector of the operator A α, i.e.: A α injπx α =λ j injπx α, 8 λ 1 λ 2 λ N 1 4 h 2. 7 If we replace e α and e α in 5 by expreion of the form given by equation 6 and ue 7, we get the following matrix equation where d jk i the column vector of pectral coefficient d1,jk d jk =, d 2,jk and Q 2 i the tability matrix of MAD cheme 0.5+λ j 1+2λ j Q 2 = 0.5 λ j 1+2λ k d jk= Q 2 djk, λ k 1+2λ j 0.5+λ k 1+2λ k.
4 One Cla of Splitting Iterative Scheme 425 Now we conider the neceary condition for the tability of iterative MAD cheme 4. Since Q 2 i not ymmetric matrix, the dicrete von Neumann tability criteria can not prove that they are alo ufficient for tability of MAD cheme. Theorem 1. All eigenvalue of tability matrix Q 2 atify inequalitie q jk < 1, 1 j, k N 1 unconditionally for any value of parameter and h. Proof. Uing imple computation we get that eigenvalue q of the amplification matrix Q 2 atify the quadratic equation q 2 q + Then the eigenvalue of Q 2 are q 1,2 = 1 1 ± 2 λ j + λ k 1 + 2λ j 1 + 2λ k =0. 1 2λ j 1 2λ k 1 + 2λ j 1 + 2λ k and they are obviouly both le than 1. The theorem i proved. The convergence order of MDAiterative cheme depend on the parameter. We will ue quai - optimal parameter 0 which olve the minimization problem w 0 =min max 1 2λ λ 1 λ λ N 1 1+2λ. Since the function ga = 1 a 1+a i trictly decreaing, we find 0 from the equation After imple computation we get 1 2λ 1 = 2λ N λ 1 1+2λ N 1 0 = 1 2 h λ 1 λ N 1 4π and all eigenvalue of the amplification matrix Q 2 atify the inequality 1 q jk π h. 9 λ 1 /λ N 1 2 Thee etimate are imilar to convergence etimate obtained for the tandard iterative cheme of alternating direction.
5 v. 426 R. Ciegi and V. Pakalnyte 4 Spectral tability analyi of 3D iterative cheme A it wa tated in ection 3, the von Neumann tability criteria give only ufficient tability condition of MAD iterative cheme. In thi ection we will ue the pectral tability analyi for 3D iterative cheme 3. But intead of finding eigenvalue of an amplification matrix Q 3 we ue direct numerical experiment and find the optimal value of the parameter. Such methodology give u a poibility to invetigate the convergence order of MAD cheme in the uual l 2 norm. Let conider the model problem 3 A α y =λ j + λ k + λ l injπx 1 inkπx 2 inlπx 3, 10 which ha the exact olution y = injπx 1 inkπx 2 inlπx 3. Then the olution of 3D MAD cheme can be repreented a y α = dα injπx 1 inkπx 2 inlπx 3, α =1, 2, 3, where d α can be computed explicitly dβ 1 3 d α= 1+3λ α 3 λ β dβ 1 +3λ α dα, α =1, 2, 3. Then the error of the th iteration d i etimated by the following formula d1 + d2 + d3 e = 1 3. Let S,λ j,λ k,λ l be the number of iteration required to achieve the accuracy e ε for given eigenvalue λ α. Then we invetigate the whole interval [m, M], which characterize the tiffne of the problem and compute S = max m λ α M S,λ j,λ k,λ l. Thi problem i olved approximately by computing S,λ j,λ k,λ l for K +1 3 combination of eigenvalue λ α = m + im m/k. Firt, we invetigated the dependence of the number of iteration S on the parameter. It wa proved that there exit the optimal value 0, uch that S 0 S
6 One Cla of Splitting Iterative Scheme 427 Table 1. The number of iteration a a function of ε S and thi value atifie the following condition max m λ α M S 0,λ j,λ k,λ l =S 0,m,m,m. 11 In Table 1 we preent number of iteration S for different value of. Thee experiment were done with m =10,M = The optimal value of the parameter depend lightly on ε. In Table 2 we preent optimal number of iteration S 0 for different value of ε. Table 2. The optimal value of a a function of ε ε 0 S Finally, we invetigated the dependence of the convergence rate of the MAD iterative cheme on m and M. In Table 3 we preent optimal value of and number of iteration S 0 for different pectral interval. We ued ε =10 4 in thee experiment. It follow from reult preented in Table 3 that 0 = c mm, S 0 =O M. m The above concluion agree well with reult of ection 3.
7 v. 428 R. Ciegi and V. Pakalnyte Table 3. The optimal value of a a function of m and M m M 0 S Error etimate for p-dimenional MAD cheme In thi ection we conider p-dimenional iterative cheme 3. Let u introduce the following notation: y y α y α αt = v 2 3 = α,, α>β, ỹ ỹ ỹ t =, v α,β = y α y β, v α,β 2, Q p y= r p 2 2 v 2 3. In the following theorem we etimate the convergence rate of MAD iterative cheme. Theorem 2. Iterative cheme 3 produce a equence converging unconditionally to the olution of problem 2 and the convergence rate i etimated a Q p y 1 q Q p y, 1 q =min 1+mp, 1+, 12 2M p where m and M are the pectral etimate of the operator A: m = min 1 α p m α, M = max 1 α p M α. Proof. Multiplying both ide of 3 by y αt and adding all equalitie we get y α ỹ, y αt + p A αyαt, y αt +p A βyβ f, ỹ t =0. 13
8 One Cla of Splitting Iterative Scheme 429 The firt term of 13 can be rewritten a I 1 = 1 p = p ỹ t p y α y β, y αt + 1 p α,, α>β v α,β, v t α,β = p ỹ t p v t p 2 v 2 3 v 2 3 By adding equation 3 we get that ỹ t = A α y α f, y β y β, y αt. 14 hence uing the third term of 13 we can prove that p ỹ t 2 + p A βyβ f, ỹ t p r 2 r The econd term of 13 i etimated imilarly. Thu the convergence rate etimate 12 follow trivially. The theorem i proved. A a corollary of Theorem 2 we can find the optimal value of parameter 0 =1/p 2mM, which i obtained from the equation 1+pm =1+ 6 Seidel-type iterative cheme 1 2M p. In thi ection we invetigate the convergence rate of Seidel-type iterative cheme: y α y α + α A β y β + β=α+1 y 1 = y 1, A β yβ = f,, 2,...,p, 16 yα y α =0.5 + y α Spectral tability analyi of 2D cheme To apply the dicrete von Neumann tability criteria to problem 16, we write the global error a a erie: e α = N 1 j=1 N 1 k=1 dα,jk injπx 1 inkπx 2, α =1, 2,
9 v. 430 R. Ciegi and V. Pakalnyte Subtituting thi expanion into 16, we obtain the equation for coefficient d jk= Q 2 djk. 17 where d jk i the column vector of pectral coefficient and Q 2 i the tability matrix of cheme 16 1 λ k 1+λ j 1+λ j Q 2 = 0.51 λ j 0.5+λ j + 2 λ j λ k. 1 + λ j 1 + λ k 1 + λ j 1 + λ k Now we conider the neceary condition for the tability of cheme 16. The eigenvalue of the amplification matrix Q 2 atify the quadratic equation q λ j 1+ q = λ j 1 + λ k 1+λ j Theorem 3. All eigenvalue of tability matrix Q 2 atify inequalitie q jk < 1, 1 j, k N 1 unconditionally for any value of parameter and h. Proof. Application of the Hurwitz criterion give that q jk 1i atified if and only if 0.5 < 1, 1+λ j λ j λ j 1 + λ k < λ j Simple computation prove that both inequalitie are atified unconditionally. The theorem i proved. 6.2 Spectral tability analyi of 3D iterative cheme Let conider the model problem 10. The olution of 3D cheme 16 can be repreented a y α = dα injπx 1 inkπx 2 inlπx 3, α =1, 2, 3, where d α,, 2, 3, are computed explicitly 1 d α= dα + α 1 dα 1 λ β d β 1+λ α 2 e = max dα 1 α 3 3 β=α+1 λ β dβ + We etimate the error of the th iteration dα by the following formula. 3 λ β.
10 One Cla of Splitting Iterative Scheme 431 Table 4. The optimal value of a a function of m and M m M 0 S We invetigated numerically the dependence of the convergence rate of the iterative cheme 16 on m and M. In Table 4 we preent optimal value of and number of iteration S 0 for different pectral interval. We ued ε =10 4 in thee experiment. It follow from reult preented in Table 4 that 0 = c mm, S 0 =O M m The convergence rate of the Seidel type cheme 16 i the ame a of cheme 3.. Reference 1. Abrahin, V.N.: On the tability of multicomponent additive direction method. Differentialnyje uravnenyja in Ruian. 2. Abrahin, V.N., Zhadaeva, N.G.: On the convergence rate of economical iterative method for tationary problem of mathematical phyic. Differentialnyje uravnenyja in Ruian. 3. Aleinikova, T., Čiegi, R.: Invetigation of new iterative method for olving multidimenional elliptic equation. Differentialnyje uravnenyja in Ruian. 4. Marchuk, G.I. Splitting method. Nauka, Mocow, in Ruian. 5. Peaceman, D., Rachford, H.: The numerical olution of parabolic and elliptic differential equation. SIAM, Samarkii, A.A., Nikolajev, E.S. Method for olving difference equation. Nauka, Mocow, 1978.in Ruian.
SOME RESULTS ON INFINITE POWER TOWERS
NNTDM 16 2010) 3, 18-24 SOME RESULTS ON INFINITE POWER TOWERS Mladen Vailev - Miana 5, V. Hugo Str., Sofia 1124, Bulgaria E-mail:miana@abv.bg Abtract To my friend Kratyu Gumnerov In the paper the infinite
More informationMulti-dimensional Fuzzy Euler Approximation
Mathematica Aeterna, Vol 7, 2017, no 2, 163-176 Multi-dimenional Fuzzy Euler Approximation Yangyang Hao College of Mathematic and Information Science Hebei Univerity, Baoding 071002, China hdhyywa@163com
More informationManprit Kaur and Arun Kumar
CUBIC X-SPLINE INTERPOLATORY FUNCTIONS Manprit Kaur and Arun Kumar manpreet2410@gmail.com, arun04@rediffmail.com Department of Mathematic and Computer Science, R. D. Univerity, Jabalpur, INDIA. Abtract:
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationA Simplified Methodology for the Synthesis of Adaptive Flight Control Systems
A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More informationON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS. Volker Ziegler Technische Universität Graz, Austria
GLASNIK MATEMATIČKI Vol. 1(61)(006), 9 30 ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS Volker Ziegler Techniche Univerität Graz, Autria Abtract. We conider the parameterized Thue
More informationON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS
ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS VOLKER ZIEGLER Abtract We conider the parameterized Thue equation X X 3 Y (ab + (a + bx Y abxy 3 + a b Y = ±1, where a, b 1 Z uch that
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More informationTail estimates for sums of variables sampled by a random walk
Tail etimate for um of variable ampled by a random walk arxiv:math/0608740v mathpr] 11 Oct 006 Roy Wagner April 1, 008 Abtract We prove tail etimate for variable i f(x i), where (X i ) i i the trajectory
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationA Full-Newton Step Primal-Dual Interior Point Algorithm for Linear Complementarity Problems *
ISSN 76-7659, England, UK Journal of Information and Computing Science Vol 5, No,, pp 35-33 A Full-Newton Step Primal-Dual Interior Point Algorithm for Linear Complementarity Problem * Lipu Zhang and Yinghong
More informationTheoretical Computer Science. Optimal algorithms for online scheduling with bounded rearrangement at the end
Theoretical Computer Science 4 (0) 669 678 Content lit available at SciVere ScienceDirect Theoretical Computer Science journal homepage: www.elevier.com/locate/tc Optimal algorithm for online cheduling
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationComputers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order
Computer and Mathematic with Application 64 (2012) 2262 2274 Content lit available at SciVere ScienceDirect Computer and Mathematic with Application journal homepage: wwweleviercom/locate/camwa Sharp algebraic
More informationSome Sets of GCF ϵ Expansions Whose Parameter ϵ Fetch the Marginal Value
Journal of Mathematical Reearch with Application May, 205, Vol 35, No 3, pp 256 262 DOI:03770/jin:2095-26520503002 Http://jmredluteducn Some Set of GCF ϵ Expanion Whoe Parameter ϵ Fetch the Marginal Value
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationConvergence criteria and optimization techniques for beam moments
Pure Appl. Opt. 7 (1998) 1221 1230. Printed in the UK PII: S0963-9659(98)90684-5 Convergence criteria and optimization technique for beam moment G Gbur and P S Carney Department of Phyic and Atronomy and
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationDISCRETE ROUGH PATHS AND LIMIT THEOREMS
DISCRETE ROUGH PATHS AND LIMIT THEOREMS YANGHUI LIU AND SAMY TINDEL Abtract. In thi article, we conider it theorem for ome weighted type random um (or dicrete rough integral). We introduce a general tranfer
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationNORMAL MODE ANALYSIS OF SECOND-ORDER PROJECTION METHODS FOR INCOMPRESSIBLE FLOWS. Jae-Hong Pyo and Jie Shen. (Communicated by Z.
DISCRETE AND CONTINUOUS Webite: http://aimcience.org DYNAMICAL SYSTEMS SERIES B Volume 5 Number 3 Augut 005 pp. 817 840 NORMAL MODE ANALYSIS OF SECOND-ORDER PROJECTION METHODS FOR INCOMPRESSIBLE FLOWS
More informationAsymptotic behavior of solutions of mixed problem for linear thermo-elastic systems with microtemperatures
Mathematica Aeterna, Vol. 8, 18, no. 4, 7-38 Aymptotic behavior of olution of mixed problem for linear thermo-elatic ytem with microtemperature Gulhan Kh. Shafiyeva Baku State Univerity Intitute of Mathematic
More informationFOURIER SERIES AND PERIODIC SOLUTIONS OF DIFFERENTIAL EQUATIONS
FOURIER SERIES AND PERIODIC SOLUTIONS OF DIFFERENTIAL EQUATIONS Nguyen Thanh Lan Department of Mathematic Wetern Kentucky Univerity Email: lan.nguyen@wku.edu ABSTRACT: We ue Fourier erie to find a neceary
More informationChapter 7. Root Locus Analysis
Chapter 7 Root Locu Analyi jw + KGH ( ) GH ( ) - K 0 z O 4 p 2 p 3 p Root Locu Analyi The root of the cloed-loop characteritic equation define the ytem characteritic repone. Their location in the complex
More informationEfficient Methods of Doppler Processing for Coexisting Land and Weather Clutter
Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr
More informationUnbounded solutions of second order discrete BVPs on infinite intervals
Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b
More informationROBUST CONTROLLER DESIGN WITH HARD INPUT CONSTRAINTS. TIME DOMAIN APPROACH. Received August 2015; revised November 2015
International Journal of Innovative Computing, Information and Control ICIC International c 2016 ISSN 1349-4198 Volume 12, Number 1, February 2016 pp. 161 170 ROBUST CONTROLLER DESIGN WITH HARD INPUT CONSTRAINTS.
More informationFebruary 5, :53 WSPC/INSTRUCTION FILE Mild solution for quasilinear pde
February 5, 14 1:53 WSPC/INSTRUCTION FILE Mild olution for quailinear pde Infinite Dimenional Analyi, Quantum Probability and Related Topic c World Scientific Publihing Company STOCHASTIC QUASI-LINEAR
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationEstimation of Current Population Variance in Two Successive Occasions
ISSN 684-8403 Journal of Statitic Volume 7, 00, pp. 54-65 Etimation of Current Population Variance in Two Succeive Occaion Abtract Muhammad Azam, Qamruz Zaman, Salahuddin 3 and Javed Shabbir 4 The problem
More informationExplicit formulae for J pectral factor for well-poed linear ytem Ruth F. Curtain Amol J. Saane Department of Mathematic Department of Mathematic Univerity of Groningen Univerity of Twente P.O. Box 800
More informationRobustness analysis for the boundary control of the string equation
Routne analyi for the oundary control of the tring equation Martin GUGAT Mario SIGALOTTI and Mariu TUCSNAK I INTRODUCTION AND MAIN RESULTS In thi paper we conider the infinite dimenional ytem determined
More informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
More informationThe combined Laplace-homotopy analysis method for partial differential equations
Available online at wwwir-publicationcom/jmc J Math Computer Sci 6 (26), 88 2 Reearch Article The combined Laplace-homotopy analyi method for partial differential equation Javad Vahidi Department of Mathematic,
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More informationStudy of a Freely Falling Ellipse with a Variety of Aspect Ratios and Initial Angles
Study of a Freely Falling Ellipe with a Variety of Apect Ratio and Initial Angle Dedy Zulhidayat Noor*, Ming-Jyh Chern*, Tzyy-Leng Horng** *Department of Mechanical Engineering, National Taiwan Univerity
More informationChapter 13. Root Locus Introduction
Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will
More informationONLINE APPENDIX: TESTABLE IMPLICATIONS OF TRANSLATION INVARIANCE AND HOMOTHETICITY: VARIATIONAL, MAXMIN, CARA AND CRRA PREFERENCES
ONLINE APPENDIX: TESTABLE IMPLICATIONS OF TRANSLATION INVARIANCE AND HOMOTHETICITY: VARIATIONAL, MAXMIN, CARA AND CRRA PREFERENCES CHRISTOPHER P. CHAMBERS, FEDERICO ECHENIQUE, AND KOTA SAITO In thi online
More informationBy Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago
Submitted to the Annal of Applied Statitic SUPPLEMENTARY APPENDIX TO BAYESIAN METHODS FOR GENETIC ASSOCIATION ANALYSIS WITH HETEROGENEOUS SUBGROUPS: FROM META-ANALYSES TO GENE-ENVIRONMENT INTERACTIONS
More informationFactor Analysis with Poisson Output
Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationLateral vibration of footbridges under crowd-loading: Continuous crowd modeling approach
ateral vibration of footbridge under crowd-loading: Continuou crowd modeling approach Joanna Bodgi, a, Silvano Erlicher,b and Pierre Argoul,c Intitut NAVIER, ENPC, 6 et 8 av. B. Pacal, Cité Decarte, Champ
More informationApplication of Laplace Adomian Decomposition Method on Linear and Nonlinear System of PDEs
Applied Mathematical Science, Vol. 5, 2011, no. 27, 1307-1315 Application of Laplace Adomian Decompoition Method on Linear and Nonlinear Sytem of PDE Jaem Fadaei Mathematic Department, Shahid Bahonar Univerity
More informationUNIQUE CONTINUATION FOR A QUASILINEAR ELLIPTIC EQUATION IN THE PLANE
UNIQUE CONTINUATION FOR A QUASILINEAR ELLIPTIC EQUATION IN THE PLANE SEPPO GRANLUND AND NIKO MAROLA Abtract. We conider planar olution to certain quailinear elliptic equation ubject to the Dirichlet boundary
More informationSobolev-Il in Inequality for a Class of Generalized Shift Subadditive Operators
Nonlinear Analyi and Differential Equation, Vol. 5, 217, no. 2, 75-88 HIKAI Ltd, www.m-hikari.com http://doi.org/1.12988/nade.217.61299 Sobolev-Il in Inequality for a Cla of Generalized Shift Subadditive
More informationOVERFLOW PROBABILITY IN AN ATM QUEUE WITH SELF-SIMILAR INPUT TRAFFIC
Copyright by IEEE OVERFLOW PROBABILITY IN AN ATM QUEUE WITH SELF-SIMILAR INPUT TRAFFIC Bori Tybakov Intitute for Problem in Information Tranmiion Ruian Academy of Science Mocow, Ruia e-mail: bt@ippi ac
More informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationTo appear in International Journal of Numerical Methods in Fluids in Stability analysis of numerical interface conditions in uid-structure therm
To appear in International Journal of Numerical Method in Fluid in 997. Stability analyi of numerical interface condition in uid-tructure thermal analyi M. B. Gile Oxford Univerity Computing Laboratory
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hoted at the Radboud Repoitory of the Radboud Univerity Nijmegen The following full text i an author' verion which may differ from the publiher' verion. For additional information about thi publication
More informationTechnical Appendix: Auxiliary Results and Proofs
A Technical Appendix: Auxiliary Reult and Proof Lemma A. The following propertie hold for q (j) = F r [c + ( ( )) ] de- ned in Lemma. (i) q (j) >, 8 (; ]; (ii) R q (j)d = ( ) q (j) + R q (j)d ; (iii) R
More informationAdvanced methods for ODEs and DAEs
Lecture : Implicit Runge Kutta method Bojana Roić, 9. April 7 What you need to know before thi lecture numerical integration: Lecture from ODE iterative olver: Lecture 5-8 from ODE 9. April 7 Bojana Roić
More informationarxiv: v1 [math.nt] 18 Jan 2016
A functional relation for L-function of graph equivalent to the Riemann Hypothei for Dirichlet L-function arxiv:60.04573v [math.nt] 8 Jan 206 Fabien Friedli September 9, 208 Abtract In thi note we define
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationDomain Optimization Analysis in Linear Elastic Problems * (Approach Using Traction Method)
Domain Optimization Analyi in Linear Elatic Problem * (Approach Uing Traction Method) Hideyuki AZEGAMI * and Zhi Chang WU *2 We preent a numerical analyi and reult uing the traction method for optimizing
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationPacific Journal of Mathematics
Pacific Journal of Mathematic OSCILLAION AND NONOSCILLAION OF FORCED SECOND ORDER DYNAMIC EQUAIONS MARIN BOHNER AND CHRISOPHER C. ISDELL Volume 230 No. March 2007 PACIFIC JOURNAL OF MAHEMAICS Vol. 230,
More informationMinimal state space realization of MIMO systems in the max algebra
KULeuven Department of Electrical Engineering (ESAT) SISTA Technical report 94-54 Minimal tate pace realization of MIMO ytem in the max algebra B De Schutter and B De Moor If you want to cite thi report
More informationSOLVING THE KONDO PROBLEM FOR COMPLEX MESOSCOPIC SYSTEMS
SOLVING THE KONDO POBLEM FO COMPLEX MESOSCOPIC SYSTEMS V. DINU and M. ÞOLEA National Intitute of Material Phyic, Bucharet-Magurele P.O. Box MG-7, omania eceived February 21, 2005 Firt we preent the calculation
More informationin a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim
Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College
More informationMulticolor Sunflowers
Multicolor Sunflower Dhruv Mubayi Lujia Wang October 19, 2017 Abtract A unflower i a collection of ditinct et uch that the interection of any two of them i the ame a the common interection C of all of
More informationLinear Quadratic Stochastic Differential Games under Asymmetric Value of Information
Preprint of the 2th World Congre The International Federation of Automatic Control Touloue, France, July 9-4, 27 Linear Quadratic Stochatic Differential Game under Aymmetric Value of Information Dipankar
More informationOn the Relationship Between Continuous- and Discrete-Time Control Systems
On the Relationhip Between Continuou- and Dicrete-Time Control Sytem V.M. Veliov Reearch Report 21-13 November, 21 Operation Reearch and Control Sytem Intitute of Mathematical Method in Economic Vienna
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationEquivalent Strain in Simple Shear Deformations
Equivalent Strain in Simple Shear Deformation Yan Beygelzimer Donetk Intitute of Phyic and Engineering The National Academy of Science of Ukraine Abtract We how that imple hear and pure hear form two group
More informationImproving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France,
More informationA FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT
A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: Zenon Medina-Cetina International Centre for Geohazard / Norwegian Geotechnical Intitute Roger
More informationThe type 3 nonuniform FFT and its applications
Journal of Computational Phyic 206 (2005) 1 5 Short Note The type 3 nonuniform FFT and it application June-Yub Lee a, *, Lelie Greengard b a Department of Mathematic, Ewha Woman Univerity, 11-1 Daehyundong,
More informationQuasi-optimal convergence rates for adaptive boundary element methods with data approximation, Part I: Weakly-singular integral equation
ASC Report No. 24/2013 Quai-optimal convergence rate for adaptive boundary element method with data approximation, Part I: Weakly-ingular integral equation M. Feichl, T. Führer, M. Karkulik, J.M. Melenk,
More information3.2. The Ramsey model - linearization, backward integration, relaxation (complete)
3.2. The Ramey model - linearization, backward integration, relaxation (complete) (Lecture Note, Thoma Steger, Univerity of Leipzig, ummer term ) Ramey Model: linearized verion. Model etup, firt-order
More informationPOINCARE INEQUALITY AND CAMPANATO ESTIMATES FOR WEAK SOLUTIONS OF PARABOLIC EQUATIONS
Electronic Journal of Differential Equation, Vol. 206 (206), No. 204, pp. 8. ISSN: 072-669. URL: http://ejde.math.txtate.edu or http://ejde.math.unt.edu POINCARE INEQUALITY AND CAMPANATO ESTIMATES FOR
More informationDigital Control System
Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)
More informationTHE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi
International Journal of Innovative Computing, Information Control ICIC International c 206 ISSN 349-498 Volume 2, Number 2, April 206 pp. 357 370 THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY
More informationResearch Article Least-Mean-Square Receding Horizon Estimation
Mathematical Problem in Engineering Volume 212, Article ID 631759, 19 page doi:1.1155/212/631759 Reearch Article Leat-Mean-Square Receding Horizon Etimation Bokyu Kwon 1 and Soohee Han 2 1 Department of
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationNonlocal initial value problems for implicit differential equations with Hilfer Hadamard fractional derivative
Nonlinear Analyi: Modelling and Control, Vol. 23, No. 3, 34 360 ISSN 392-53 http://doi.org/0.5388/na.208.3.4 Nonlocal initial value problem for implicit differential equation with Hilfer Hadamard fractional
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationarxiv: v3 [math.na] 20 Dec 2017
On a plitting method for the Zakharov ytem Ludwig Gauckler 1 Verion of 15 November 217 arxiv:167.7556v3 [math.na] 2 Dec 217 Abtract An error analyi of a plitting method applied to the Zakharov ytem i given.
More informationSecond order conditions for k means clustering: partitions vs. centroids
Second order condition for k mean clutering: partition v. centroid Ziqiu Su Jacob Kogan Abtract Familie of centroid baed clutering method are capable of handling high dimenional pare data ariing in Text
More informationFinding the location of switched capacitor banks in distribution systems based on wavelet transform
UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh
More informationLinear System Fundamentals
Linear Sytem Fundamental MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Content Sytem Repreentation Stability Concept
More informationOn Uniform Exponential Trichotomy of Evolution Operators in Banach Spaces
On Uniform Exponential Trichotomy of Evolution Operator in Banach Space Mihail Megan, Codruta Stoica To cite thi verion: Mihail Megan, Codruta Stoica. On Uniform Exponential Trichotomy of Evolution Operator
More informationOn The Approximate Solution of Linear Fuzzy Volterra-Integro Differential Equations of the Second Kind
AUSTRALIAN JOURNAL OF BASIC AND APPLIED SCIENCES ISSN:99-878 EISSN: 39-8 Journal ome page: www.ajbaweb.com On Te Approimate Solution of Linear Fuzzy Volterra-Integro Differential Equation of te Second
More informationA PROOF OF TWO CONJECTURES RELATED TO THE ERDÖS-DEBRUNNER INEQUALITY
Volume 8 2007, Iue 3, Article 68, 3 pp. A PROOF OF TWO CONJECTURES RELATED TO THE ERDÖS-DEBRUNNER INEQUALITY C. L. FRENZEN, E. J. IONASCU, AND P. STĂNICĂ DEPARTMENT OF APPLIED MATHEMATICS NAVAL POSTGRADUATE
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More information