arxiv: v1 [quant-ph] 19 Jan 2014

Size: px

Start display at page:

Download "arxiv: v1 [quant-ph] 19 Jan 2014"

Alfred Ford
6 years ago
Views:

1 Analytical Solution of Mathieu Equation Ditri Yerchuck (a), Alla Dovlatova (b), Yauhen Yerchak(c), Felix Borovik (a) (a) - Heat-Mass Transfer Institute of National Acadey of Sciences of RB, Brovka Str., 15, Minsk, 007, RB (b) - M.V.Loonosov Moscow State University, Moscow, (c) Belarusian State University, Nezavisiosti Avenue 4, Minsk, 0030, RB (Dated: January, 014) The general solution of the hoogeneous daped Mathieu equation in the analytical for, allowing its practical using in any applications, including superconductivity studies, without nuerical calculations has been found. arxiv: v1 [quant-ph] 19 Jan 014 PACS nubers: Keywords: Mathieu equation, superconductivity I. INTRODUCTION AND BACKGROUND It is known, that a nuber of physical phenoena can be described atheatically by Mathieu equation. For instance, the equation of otion for the flux lattice in the supeconductor theory is the following [1] d y(t) +η dy(t) +( +kcosωt)y(t) = cos Ωt, c (1) where B is the agnetic induction in the saple, is the total ass (per unit length) of the flux lattice, η is the viscosity coefficient, ω is the frequency of the odulating agnetic field, c is the light velocity, J 0 and Ω are the aplitude and the frequency of the icrowave current, respectively, and k are aplitudes of constant coponent and an alternating coponent in the function K(t) in the relationship F(t) = K(t)y(t) () between the force F(t) and sall displaceent y(t) of the flux lattice fro its equilibriu position, that is K(t) = +kcosωt. (3) The equation (1) is the inhoogeneous daped Mathieu equation. Mathieu equation is well know in the theory of differential equations, see for exaple [], [3], [4], [5]. At that,the so-called general Mathieu equation is the following equation where P(t) is a periodical function with the period, equal to π, µ is so-called characteristic index, depending on the values of h and θ. The function y(t) = exp( µt)p( t), (6) represents itself the second solution. The solutions (5) and(6) are linearly independent, they produce the fundaental syste of the solutions, excluding the case, when iµ Z, that is, to the set of whole nubers [4]. Further, the solution (5) is written forally in the for of the following infinite series y(t) = n= c n exp(µ+ni)t, (7) in which for the set of coefficients {c n } the following recurrent relations θc n 1 +[h+(µ+ni) ]c n θc n+1 = 0,n Z (8) were obtained. It is the only algorith for the nuerical solution. In fact, the analytical solution of the general Mathieu equation (4) was not found. The existence the only algorithn for the nuerical solution is inconveniently for the practical usage of Mathieu equation, especially in the cases of physical applications when analytical dependences are required to understand the physical processes. For instance, the authors of the work [1] have preferred instead of trying to solve the equation (1) nuerically to use the linearized equation d y(t) +[h θcos(t)]y(t) = 0, (4) where h and θ are real or coplex constants. The known solution of the general Mathieu equation (4) is built in the for y(t) = exp(µt)p(t), (5) d δy +ηdδy+δy = kcosωty 0, (9) in which y 0 is the particular solution of (1) for the case of k = 0. Let us reark, that an (in)hoogeneous daped Mathieu equation (4) can be reduced to an (in)hoogeneous general Mathieu equation by the following transforation of the function y(t): y(t) = w(t)exp( 1 ηt).

2 The authors of the work [1] have found, to the first order in k, η, and ω, that the value of an electric field E(t) induced is E(t) = B J 0 Ω [1 ǫcos(ωt φ)]sin(ωt α), (10) c where ǫ = k, tanφ = ηω, and tanα = ηω. So, the equation (10) indicates that the induced E(t) field is aplitude odulated. In other words, the concrete physical echanis and its analytical description were established, alhough for rather restricted ranges of the paraeters used. Really, the equation (10) is valid for k, ω Ω, and Ω. At the sae tie, the authors of [1] reark, that when static agnetic field H becoes too large [larger than 800 G], the flux structure becoes probably too coplex for the siple odel proposed to reain valid. It is clear on the given exaple, that an analytical solution of Mathieu equation reains to be very actual for its applications in physical sciences and in engineering. On the other hand, an analytical solution of Mathieu equation has also the atheatical theoretical aspect. It is deterined by the fact that the solution of a nuber of differential equations is reduced to the solution of Mathieu equation. They, for exaple, are 1) (1 t ) d y(t) +t dy(t) +(at +b)y(t) = 0, (11) The transforation of variable t = cosz leads to the Mathieu equation ) d y(z) dz +(a+b+acosz)y(z) = 0, (1) t(t 1) d y(t) +(t 1) dy(t) +(at+b)y(t) = 0, (13) The transforation of variable t = cos z leads to the Mathieu equation 3) d y(z) dz (a+b+acosz)y(z) = 0, (14) d y(t) +(asinλt+b)y(t) = 0, (15) The transforation of variable λt = z + π leads to the Mathieu equation 4) The equations d y(z) dz ( 4b λ + 4a λ cosz)y(z) = 0, (16) d y(t) +(asin t+b)y(t) = 0, d y(t) +(acos t+b)y(t) = 0 (17) are transfored to Mathieu equations by using of trigonoetric forulae sin t = 1 cost, cos t = 1+cost (18) correspondingly. The ai of the given work is to find the analytical solution of the general Mathieu equation. II. RESULTS Theore. The general solution of the general Mathieu equation can be represented analytically to be the superposition of Bessel functions of the first and the second kinds. Proof. It is evident, that the equation (1) can be represented in the for d + 1 ηdy + 1 y + 1 kexp(t)y c cosωt+ [ d + 1 ηdy + 1 y + 1 kexp( t)y c cosωt] = (ˆL 1 + ˆL )y = 0, where ˆL 1 and ˆL are differential operators (19) ˆL 1 = 1 d + 1 η d kexp(t) c cosωt ˆL = [ 1 d + 1 η d kexp( t) c cosωt]. (0) It is also evident that the partial solution of the starting equation will correspond to the intersection of sets of the solutions satysfying siultaneously to the equations

3 3 ˆL = 0, ˆL y = 0. So, we have to solve the equations d + 1 ηdy + 1 y + 1 kexp(t)y = c cosωt, d + 1 ηdy + 1 y + 1 kexp( t)y = c cosωt or in the equivalent for The case d y + η dy + y + k exp(t)y = c cosωt, d y + η c cosωt, dy + y + k exp( t)y = (1) () (3) that is the siple differential equation of the first order, the solution of which is evident. The hoogeneous equations, corresponding to inhoogeneous equations (1) that is the equations d y + η dy + y + k exp(t)y = 0 d y + η dy + y + k exp( t)y = 0 (7) can be solved strictly. Let us designate η = a K0 = b, k = c, = λ, = λ. Then, according to [6] the solutions are y(t) [1] = exp( η t)[c[1] Y ν( exp(t ))] 1 J ν( for the first equation in (7) and C [] Y ν ( exp( t ))], exp(t ))+ (8) y(t) [] = exp( η t)[c[] 1 J ν ( exp( t ))+ (9) d + 1 ηdy + 1 y + 1 kexp(t)y c cosωt = [ d + 1 ηdy + 1 y + 1 kexp( t)y c cosωt] (4) can be also studied. There are a nuber of variants to siplify the solution of (4) using the syetry of y(t). In the case, if y(t) is even, dy is uneven, d y is even fro (4) we obtain by t t d y +y + 1 kexp(t )y + 1 kexp( t )y = cosωt. c (5) Given case is equivalent to the task above forulated by the relations (19) - (1). Further, if y(t) is even, dy is uneven, d y is even fro (4) we obtain by t t η dy cosωt = 0, (6) c for the second equation in (7). Here J ν ( Y ν ( exp(t )),J ν ( exp( t )), exp(t )),Y ν ( exp( t )) (30) are Bessel functions of the first kind and of the second kind of index ν = ( η ) 4 k for the equation (30) and of index ν = ( η ) 4 k for the equation (31). Then the set of functions, satisfying to the relation (31) (3) y(t) [1] y(t) [] (33)

4 4 will be the solution of the hoogeneous equation, corresponding to the starting inhoogeneous equation (1). Given conclusion is correct, if y(t) [1] y(t) []. In particular, it is satisfied for the case 1 = 0, C [] 1 = 0 and by Y ν( C [] that is, in the case J ν ( C [] K0 Y ν( exp(t )) = exp( t )), exp( t ))cosπν J ν( K0 J ν ( K0 exp( t ))cosπν +J ν( exp( t K0 (34) )) = exp( t )), (35) Hence we obtain the conditions, by which the relation exp( η t)y ν( exp(t )) (36) will be the solution of the hoogeneous daped Mathieu equation. They are that is, by cosπν C[] cosπν C[1] cos πν 1 0. = 0, = 0, = 0 (37) (38) Consequently, ν Z, that is ν has to belong to the set of whole nubers. Therefore, the relation (36) is the solution of hoogeneous daped Mathieu equation by 0 C and by ν Z. It eans, that there are restrictions on possible values of paraeters in the starting hoogeneous daped Mathieu equation. They are the following It is evident that by η = 0 we obtain the solution of the hoogeneous undaped Mathieu equation, that is the solution of the general Mathieu equation. To obtain the fundaental syste of the solutions of the hoogeneous daped Mathieu equation, we have to find the second linearly independent solution. It can be easily done, if to take into account, that J ν ( ( 1) ν J ν ( exp( t )) = J ν( exp(t )) = exp(t )) (41) for ν Z. Since the replaceent is equivalent to J ν ( exp( t )) J ν( exp(t )), (4) we obtain, that by = 0, C [] = 0 the expression (33) that is 1 exp( η t)j ν( exp(t )) (43) will be the second linearly independent solution of the hoogeneous daped Mathieu equation 1 0 C and by ν Z at 1 = C [] 1 for even ν Z and at 1 = C [] 1 for uneven ν Z, that is at 1 = C [] 1 sgn( 1)ν. Its linear independence fro the first solution follows fro linear independence of Bessel functions on the first and the second kinds. At fixed 1 and the relations (36) and (43) represent theselves the fundaental syste of the solution of hoogeneous daped Mathieu equation. So, the general solution of hoogeneous daped Mathieu equation is exp( η t)[c[1] Y ν( 1 J ν( exp(t ))] exp(t ))+ (44) that is ( η ν = ( ) 4 k ) Z, 4 k ν = ( ( η ) ) Z ω (39) (40) The theore is proved. Thus, we have found the general solution of the hoogeneous daped Mathieu equation in the essentially ore siple for, allowing its practical using in any applications without nuerical calculations, since Bessel functions are well known.

5 5 [1] Ki P.W., de Graaf A M, Chen J. T., Friedan E.J, and Ki S. H.,Phase Reversal and Modulated Flux Motion in Superconducting Thin Fils, Phys.Rev.B, 6, N 3 (197) [] Kotowski G, Z.Angew.Math.Mech., 3 (1943) 13 [3] Meixner J and Schäfke F W, Mathieusche Funktionen und Sphäroidfunktionen, Springer, Berlin,1954 [4] Batean H,Erdelyi A, Higher Transcendental Functions, V.3, New York, Toronto, London, MC Graw-Hill Book Copany, Inc., 1955, M., Nauka, 1967, 300 pp [5] Mac Lachlan, Theory and Applications of Mathieu Functions, M., Izd.In.Lit., 1953, 476 pp [6] Zaitsev V F, Polyanin A D, Hand-Book on Ordinary Differential Equations, M, Fizatlit, 001, 576 pp

SIMPLE HARMONIC MOTION: NEWTON S LAW

SIMPLE HARMONIC MOTION: NEWTON S LAW siple not siple PRIOR READING: Main 1.1, 2.1 Taylor 5.1, 5.2 http://www.yoops.org/twocw/it/nr/rdonlyres/physics/8-012fall-2005/7cce46ac-405d-4652-a724-64f831e70388/0/chp_physi_pndul.jpg