Local Dynamics of Logistic Equation with Delay and Diffusion

: The behavior of all the solutions of the logistic equation with delay and diffusion in a sufﬁciently small positive neighborhood of the equilibrium state is studied. It is assumed that the Andronov–Hopf bifurcation conditions are met for the coefﬁcients of the problem. Small perturbations of all coefﬁcients are considered, including the delay coefﬁcient and the coefﬁcients of the boundary conditions. The conditions are studied when these perturbations depend on the spatial variable and when they are time-periodic functions. Equations on the central manifold are con-structed as the main results. Their nonlocal dynamics determines the behavior of all the solutions of the original boundary value problem in a sufﬁciently small neighborhood of the equilibrium state. The ability to control the dynamics of the original problem using the phase change in the perturbing force is set. The numerical and analytical results regarding the dynamics of the system with parametric perturbation are obtained. The asymptotic formulas for the solutions of the original boundary value problem are given.


Introduction
The logistic equation with delaẏ is the natural generalization of the classical logistic equation. For example, it describes well the dynamics of biological population level changes [1][2][3][4]. All parameters in (1) are positive, and only non-negative solutions are considered, i.e., u(t) ≥ 0. The parameter r is called the Malthusian coefficient, T is the delay time (for example, the age of individuals puberty). The parameter a characterizes the resistance of the external environment. The literature on the study of the Equation (1) is extensive (see, for example, Refs. [5][6][7][8][9]). We note that the zero solution in (1) is unstable, and the positive equilibrium state u 0 = a −1 is asymptotically stable under the condition 0 < rT ≤ π 2 . It was shown in [5] that the equilibrium state is globally stable for 0 < rT < 3 2 , i.e., each solution with a positive initial function tends to u 0 as t → ∞. In the same paper [5], it was hypothesized that the global stability takes place in the wider range of parameters 0 < rT ≤ 37 24 . This hypothesis is proved in [7,8]. The equilibrium state is unstable for rT > π 2 , and (1) contains the stable cycle u 0 (t). The asymptotic behavior of the cycle has the form u 0 (t) = a −1 for 0 < rT − π 2 1. The cycle u 0 (t) is relaxation for rT 1. Its asymptotic behavior is given in [9]. We also note that the number of unstable cycles in (1) grows indefinitely as rT → ∞.
We consider the logistic equation with delay and diffusion with the boundary conditions We fix the space C [−T,0] × W 2 2[0,1] as the space of initial conditions that are the nonnegative functions φ(t, x) for t ∈ [−T, 0], x ∈ [0, 1] satisfying the condition (4). The boundary value problem (3) and (4) naturally generalizes the Equation (1) for the problems of spatial interaction within one population. The boundary value problems of this type have been studied by many authors (see, for example, Refs. [10][11][12]). A fairly complete overview is given in [10].
First, we assume that all the coefficients are independent of t and x. The boundary value problem (3) and (4) linearized on u 0 has the form The location of the roots of the characteristic quasi-polinomial for (5) is of great importance for the stability of the equilibrium state u 0 of the boundary value problem (3) and (4). All roots of (6) have negative real parts under the condition 0 < rT < π 2 , and (6) has the root with a positive real part as rT > π 2 . Thus, in the first case, the solution u 0 is asymptotically stable, and is unstable in the second case. We consider the critical case when the equality holds for d > 0, r = r 0 and T = T 0 . In this case, the Equation (6) has a pair of pure imaginary roots ±iω 0 for k = 0 where ω 0 = 2π/T 0 , and all its other roots have negative real parts. We present the observation of the boundary value problem (3) and (4) dynamics in the neighborhood of the equilibrium state u 0 for the coefficients d and a, which are close to some values d 0 and a 0 , and the parameters r and T close to r 0 and T 0 for which the equality (7) holds. Let where is the small positive parameter We note that the positive equilibrium state u 0 = u 0 ( ) has the form u 0 = a − 0 1 1 − a 1 a −1 0 + . . . . Under the above conditions, the two-dimensional stable local invariant integral manifold M( ) exists in some sufficiently small and -independent neighborhood of the equilibrium state u 0 ( ) [13,14]. On it, the boundary value problem (3) and (4) can be represented in the form up to terms of the order O( 2 ) where τ = t, and the equalities are met for the coeffiicients λ 1 and σ.
In this paper, the Andronov-Hopf bifurcation is studied in more complex situations. In Sections 2-4, we suppose that the coefficients in (3) and (4) are close to the constants and depend on the spatial variable.
where all the functions are sufficiently smooth. The Andronov-Hopf bifurcation under the conditions (13) is considered in Section 2. In Section 3, the perturbations also concern the boundary conditions. The Dirichlet boundary conditions are considered in Section 4. Section 5 considers non-autonomous perturbations. In particular, the problem of parametric resonance is considered.

Andronov-Hopf Bifurcation under Neumann Type Boundary Conditions
First, we find the asymptotic behavior of the positive equilibrium state K(x, ) as → 0. The equation holds for it. Since the value K ≡ 1 is a simple root of the nonlinear function in (14) for = 0, it is natural to look for K(x, ) in the form We substitute (15) into (3) and (4) and collect the coefficients at the same powers of . We successively find K 1 (x), K 2 (x), . . . from the resulting relations. So, collecting the coefficients at the first power of , we obtain From here, The following result holds. Theorem 1. The boundary value problem (3) and (4) has the positive equilibrium state K(x, ) for which the asymptotic equality (15) holds.

Linear Analysis
We linearize (3) and (4) at K(x, ). As a result, we obtain the boundary value problem ∂v ∂x x=0 Setting v = exp(λt)w(x), we obtain the characteristic equation where For = 0, we have the equation with the constant coefficients where The equalities µ 0 n = d 0 (πn) 2 , w n (x) = cos πnx, n = 0, 1, 2, . . . hold for the eigenvalues µ 0 n and the eigenfunctions w n (x) of the boundary value problem (20). Taking into consideration the relation (21), we obtain the following statement for the roots of the equations We note that the eigenfunction w 0 (x) ≡ 1 corresponds to the roots λ ± 0 . We find the asymptotic behavior of those two complex conjugated roots λ ± ( ) that tend to ±iπ(2T 0 ) −1 as → 0. Let in (19).
Collecting the coefficients at the first power of , we obtain the equation where The boundary value problem (25) is solvable if and only if the equality holds. Taking into consideration (26), we come to the conclusion that Collecting the coefficients at 2 (in (19) with the equalities (23) and (24)), we obtain λ 2 , . . . The formula (27) can be simplified. From (16), it follows that 1 0 (a 1 (x) + K 1 (x))dx = 0;

Boundary Value Problem (3) and (4) Dynamics in Neighborhood of Equilibrium K(x, ) under Conditions (13)
The characteristic Equation (19) has two roots close to the imaginary axis under these conditions, and all of its remaining roots lie to the left of the imaginary axis and are separated from this axis as → 0. Then, in the small neighborhood of K(x, ), the local stable two-dimensional integral manifold exists (see, for example, Refs. [13,14]) where the boundary value problem (3) and (4) can be written as one complex first-order equation to within O( 2 ). Here, τ = t is a "slow" time, and the "amplitude" ξ(τ) is related to the solutions of (3) and (4) by the asymptotic equality u(t, x, ) = K(x, ) where u 3 is periodic with respect to t. We substitute (29) in (3). Then, collecting the coefficients at the same powers of , we obtain the correct equality for 1/2 . At the second step, equating the coefficients of 1 on the left and right sides of the formal equality, we obtain the expressions for u 20 and u 21 : At the next step, we arrive at the equation in u 3 (t, τ, x). From the condition of its solvability in the class of π/2T 0 -periodic with respect to t functions, we obtain the equality (28) for ξ(τ), where β = λ 1 , and γ = σ as in the Equation (10).
Let us summarize what has been said.
Then, for all sufficiently small , all the solutions from someindependent neighborhood of the equilibrium K(x, ) tend to it as t → ∞. If λ 1 > 0, then there is the stable cycle u 0 (t, x, ) in the neighborhood of K(x, ) for which We note that the coefficients d 1 (x) and a 1 (x) do not affect the values of the coefficients in the Equation (28). The coefficients r 1 (x) and T 1 (x) determine the value of β = λ 1 only through their mean values on the segment x ∈ [0, 1].

Bifurcations under Disturbances in Boundary Conditions
Instead of the Neumann boundary conditions (4), we consider the "perturbed" boundary conditions for the Equation (3). Here, α 1 , α 2 , γ 1 , γ 2 are some fixed constants. For simplicity, all the coefficients in (3) are assumed constant for simplicity: d = d 0 , a = 1, r = r 0 , T = T 0 and r 0 T 0 = π 2 . First, we find the asymptotic behavior of the positive equilibrium K(x, ) of the boundary value problem (3) and (31). For this, we put K(x, ) = 1 + K 1 (x) + . . .. We substitute this expression in (3) and (31). Then, we obtain the boundary value problem . We find from here that The characteristic equation of the boundary value problem linearized on K(x, ) has a pair of pure imaginary roots ±iπ(2T 0 ) −1 as = 0. All its other roots have negative real parts.
We find the asymptotic behavior of those two complex conjugate roots λ ± ( ) of this characteristic equation0 which tend to λ ± 0 as → 0. Let 0 We look for the eigenfunction corresponding to the root λ + ( ) in the form We substitute (33) and (34) into Equation (19) with the boundary conditions Collecting the coefficients at the first power of in the resulting formal identity, we obtain the boundary value problem (19) From this, we find that At the final stage, we use the line of reasoning from the previous section. We substitute the formal series (29) in (3) and (31) and collect the coefficients at the same powers of . For 1/2 , we obtain the correct equality. At the next step, we find the expressions for u 20 and u 21 which are also defined by the equalities (30). Then, we obtain the equation for u 3 (t, τ, x). From its solvability condition in the class of 4T 0 -periodic with respect to t functions, we obtain the Equation (28) where It remains to repeat Theorem 2 with λ 1 replaced by λ 10 as a final statement. We note that the parameters of the boundary conditions determine the value of λ 10 , and all bifurcation effects depend only on the relationα 2 −α 1 , i.e., on α 2 − α 1 + γ 1 − γ 2 .

Andronov-Hopf Bifurcation in Case of Dirichlet Boundary Conditions
We consider the bifurcations in the equilibrium neighborhood of the boundary value problem Here, we assume that all the coefficients are close to the constant values The values of d 0 , r 0 and T 0 are determined below from the condition for the existance of the critical case in the zero equilibrium (38) and (39) stability problem. The characteristic equation of the linearized at zero boundary value problem has the form where µ = −λ − r exp(−λT). The smallest value µ 0 of the boundary value problem (40) is equal to π 2 d 0 as = 0, and the corresponding eigenfunction is sin πx.
The bifurcation values d 0 , r 0 and T 0 are found from the condition λ = iω 0 , i.e., from the equation From this, we obtain The value ω 0 is determined from the following equation, and the value r 0 is determined from the equality Below, we fix the positive constants d 0 and T 0 arbitrarily and denote by r 0 , the expression (42).
For small values of , the eigenvalue λ( ), which is close to iω 0 as → 0, and the corresponding eigenfunction w(x, ) are sought in the form We substitute (43) in (38) and (39) and collect the terms at the first power of to find λ 1 and w 1 (x). Then, The boundary value problem (44) is solvable if and only if the equality holds. Given (45) in it, we find the value of λ 1 : After that, we define w 1 (x): We use the above formal description when investigating the solutions of (38) and (39) from a sufficiently small neighborhood of the zero equilibrium. The series is an analogue of the asymptotic series (29). The expression complex conjugate to the previous term is denoted by cc here. We substitute (48) into (38) and (39). Performing standard actions, we obtain the equalities for u 20 (x) and u 21 (x) at the second step: Finally, we obtain the equations for u 30 (x) and u 31 (x) at the next step. From their solvability condition, we obtain the Equation (10) in which the formula (46) holds for λ 1 , and the formula holds for σ.
The main statement repeats the formulation of Theorem 2 in this case too.

Remark 1.
Unlike the results of the previous sections, the role of the coefficient d 1 (x) increases for the boundary conditions (39). Its mean value is included in the formula (46) for the coefficient λ 1 .

Remark 2.
The next statement follows from the results presented above related to the linearized problem for (38) and (39). For any T 0 > 0 and r 0 > π 2T 0 , we can choose such d 0 that the Andronov-Hopf bifurcation condition holds: the characteristic equation has a pair of roots on the imaginary axis for = 0, and all its other roots have negative real parts. In particular, this means that the increase in the diffusion coefficient leads to zero equilibrium stability in (38) and (39).

Bifurcations under Time-Periodic Perturbations
In order to illustrate more clearly the new features that arise in the case of time-periodic perturbations in the logistic equation with delay and diffusion, we investigate the simplest cases. We present the main general results in Section 5.1. The obtained results are applied to solve the problem of parametric resonance in Section 5.2. A more complicated problem of parametric resonance with a two-frequency perturbation is studied in Section 5.3.

Main Results
We consider the equation with the periodic boundary conditions Here, r 0 = π 2 , d 0 > 0 and the functions α j (t) (j = 1, 2) are periodic with the period h. In addition, α j (t) = α j0 + α j1 (t) and For = 0, the critical case occurs in the problem of zero equilibrium stability in (49) and (50): the corresponding characteristic equation has two roots λ ± 0 = ±i π 2 on the imaginary axis, and all its other roots have a negative real part. According to the above technique for studying the solutions from zero neighborhood of the boundary value problem (49) and (50), we introduce into consideration the formal series Here, τ = t is a "slow" time, and the dependence on the argument t is 4-periodic. We substitute (52) in (49) and (50) and collect the coefficients at the same powers of in the resulting formal identity. We obtain the correct equality for 1/2 . For the first power of , we obtain the equation for u 2 (t, τ), and the formula for this function coincides with (30). Collecting the coefficients of 3/2 at the third step, we obtain the equation for u 3 : We use the criterion for the boundary value problem (53)-(55) solvability in the class of 4-periodic with respect to t functions. For this, it is necessary and sufficient that the equality holds. Here, the value σ is the same as above. We note that the Equation (56) takes the simplest form when the functions α 1,2 (t) do not contain harmonics with the doubled frequency of the natural oscillations exp ±2i π 2 t . Here, A = 1 + i π 2 −1 d 0 (α 20 − α 10 ). The solutions of (57) can be obtained in an explicit form. Under the condition all the solutions of (57) tend to zero as t → ∞. Hence, we conclude that for sufficiently small , all the solutions of (49) and (50) from a sufficiently small and independent of neighborhood of the zero equilibrium also tend to zero as t → ∞. If A > 0, then (57) contains the stable cycle This means that the zero equilibrium of the boundary value problem (49) and (50) is unstable and the stable cycle Then, we consider the case when the functions α 1,2 (t) have harmonics with the frequencies ±iπt, i.e., where b j (t) do not contain harmonics with zero and ±iπ frequencies. The Equation (56) takes the form

Problem of Parametric Resonance
In this section, the results obtained in Section 5.1 are used to study the problem of parametric resonance. Let the formulas α j (t) = α j0 + b j cos ωt (59) hold for the periodic functions α 1,2 (t).
We assume here that α 20 − α 10 < 0, i.e., the zero solution in (49) and (50) is asymptotically stable in the absence of nonautonomous components in (51), and for all sufficiently small , all the solutions of (49) and (50) from a sufficiently small and independent neighborhood of zero equilibrium tend to zero as t → ∞. The same conclusion is valid under the conditions when b j = 0 and the inequality holds for the arbitrary parameter ω.
We consider the case of ω = π. The condition of zero solution stability in (49) and (50) is defined by the zero solution stability in (58). The corresponding criterion is formulated in terms of the matrix eigenvalues. By virtue of (60), we have the equality tr C = 2 A = 2d 0 1 + π 2 4 −1 (α 20 − α 10 ) < 0. Therefore, both eigenvalues of the matrix C have negative real parts under the condition Thus, the zero equilibrium in (58) and in (49) and (50) is asymptotically stable. If ∆ < 0, then the zero equilibrium in (58) and in (49) and (50) is unstable (for small ). We assume below that We investigate the problem of parametric resonance, i.e., the problem of the parameter ω entire set selection for which the zero solution of (49) and (50) is unstable. For this, we arbitrarily fix the parameter z and put in (49) and (50).
To ensure that the oscillation frequency does not depend on under the boundary conditions (50), we change the time (1 + z)t = t 1 . (64) As a result, we obtain the boundary value problem To study the solutions of this boundary value problem, we use the above algorithm. Let in (65) and (66) where τ = t and u j (t, τ, x) are 4-periodic with respect to t. Again, we obtain the formula (30) for u 2 (t, τ), and from the equation solvability condition with respect to u 3 , we get the final equation for determining ξ ∂ξ ∂τ where the values A, B and σ are the same as in (58). The stability of the zero solution is determined by the eigenvalues of the matrix We note that tr C 1 = tr C < 0 and det C 1 = π 2 4 z 2 + mz + ∆ where m = d 0 (α 20 − α 10 ) 2 1 + π 2 4 −1 . Let us formulate a simple statement about the roots of the equation det C 1 = 0.

Lemma 2. The inequality
holds under the conditions The inequalities (70) define the set of those values of the frequency ω of (63) for which the zero solution in (65) and (66) (for small ) and in (68) is unstable. It is at these frequencies that parametric resonance is observed. The dynamics of the (68) solutions have been studied by many authors (see, for example, Refs. [18,19]).

Example. Dynamics at Same Frequencies and Different Phases of External Influences
As an example, we consider a slightly different from (59) situation when where φ ∈ [0, 2π] and the equality holds.
Performing the above sequence of actions in order to determine the local dynamics of the boundary value problem (49) and (50), we obtain the normal form, i.e., the equation for the unknown "amplitude" ξ(τ) (τ = t): in which the coefficients α and σ are the same as in (68), and Replacing ξ = v exp i π 2 wτ , we obtain the equation and α 0 = α − 1 2 iπw. The stability of the zero solution in (73) is determined by the sign of the expression ρ(φ) where The zero solution in (73) is asymptotically stable under the condition ρ(φ) > 0, and it is unstable when ρ(φ) < 0. It remains to note that the function ρ(φ) is alternating in some range of parameters variation in (73) as φ ∈ [0, 2π]. Thus, the phase shift in (71) can significantly change the dynamic behavior of the original boundary value problem (49) and (50).

Parametric Resonance for Two-Frequency Perturbation
In contrast to the periodic perturbations in (59), here we assume that the frequencies of each of the function α 1 (t) and α 2 (t) are different: In addition, we assume that the inequality (60) holds and each of the values ω 1 and ω 2 is close to the doubled natural oscillations frequency, i.e., We note that the paper [20] investigated the problem of the parametric resonance for a two-frequency perturbation in the linear equation with delay. In the paper [21], this issue was studied in a nonlinear statement of the problem for the logistic equation with delay where the coefficients r and T were perturbed by two-frequency actions.
The algorithmic approach here repeats the constructions of the previous section. We consider the formal asymptotic series (67) and substitute it in (49) and (50) taking into account the equalities (74) and (75). We collect the coefficients at the same powers of in the resulting formal identity. At the second step, we obtain the formulas (30) here too. From the condition of solvability of the equation with respect to u 3 , we obtain the equation for ξ(τ) determining: )). The parameter σ is the same as in (58). Let, for definiteness, w 2 > w 1 . We perform the following changes in (76): As a result, we obtain the equation with 2π-periodic coefficients We recall that α 0 < 0, σ < 0. We investigate the dynamics of the Equation (77) for various values of the parameters µ, w 1 , α 1 , α 2 . It should be noted that the "crude" periodic solution of (77) corresponds to the two-dimensional torus of the same stability.
For sufficiently large µ, we apply the well-known averaging principle [22]. The stability properties are then determined by the averaged equation Thus, the case of large values of µ is reduced to the case of the parametric resonance with a single-frequency perturbation.
Then, according to the papers [23,24] results, the stability of the zero solution of the system (78) is closely related to the behavior of the function Three options are distinguished in the system (78) study depending on whether the function ρ(τ 1 ) is positive, negative, or alternating. Taking into account the form of the matrix B(τ 1 ) in (78), we have Thus, any of these conditions can be fulfilled by appropriate change of α 0 . Following [21], we describe the dynamics of the system (78). As it turned out, there are three following options for the behavior of this system solutions: 1. Let ρ(τ) > 0 for all τ > 0. Then, there exists µ 0 > 0 such that for 0 < µ ≤ µ 0 , the zero solution of the system (78) is asymptotically stable. 2. Let ρ(τ) < 0 for all τ > 0. Then, there exists µ 0 > 0 such that for 0 < µ ≤ µ 0 , the zero solution of (78) is asymptotically stable (unstable) if 3. Finally, in the case of the alternating function ρ(τ), the problem (78) is a system with turning points. Moreover, if the condition holds then the stability and instability of the system (78) solutions alternate infinitely for µ → 0. (We note that if the inequality is replaced by the opposite in (82), then the zero solution is asymptotically stable for sufficiently small µ). The properties listed above of the system (78) solutions allow us to draw the conclusions about the local stability or instability of the zero solution of the nonlinear Equation (77) for sufficiently large or small enough µ. The phase rearrangements of the Equation (77) are not asymptotically analizable at conditionally "average" values of the parameter µ. In this regard, it is appropriate to apply numerical methods. The acceptability of their application for finding the stable modes of the Equation (77) is due, in particular, to its dissipativity ( σ < 0). The latter property makes it possible to choose the initial conditions only from some neighborhood of the zero point of the phase plane. The detailed numerical study of the Equation (77) is given in [21]. The following conclusions were formulated on its basis.

•
In the case of the parametric resonance with a two-frequency perturbation, the dynamics is fundamentally more complicated than in the case of a single-frequency perturbation. Cycles and regions (in the parameter space) of the irregular behavior of the solutions appear much more complex in their form.

•
If ρ(τ 1 ) < 0 i.e., the asymptotic stability of the zero solution exists, no non-stationary steady-state solutions were found. • The dynamic behavior is relatively simple for the sufficiently large µ. We have the complicated relaxation stable cycles for the sufficiently small µ. They are either symmetric or two stable symmetric cycles coexisting with respect to each other. • Irregular oscillations are specific for the "average" values of the parameter µ. • Irregular oscillations in (77) are realized in the relatively narrow ranges of µ variation. In this case, the phase rearrangements are related with the symmetry-loss bifurcations, with the cascades of period-doubling bifurcations and with the bifurcations of separatrix splitting.

Delayed Perturbation Influence on Boundary Value Problem (49) and (50) Dynamic Behavior
We suppose that the small perturbation containing the delay is added to the right-hand side of the Equation (49) Here, d > 0 and r = π 2 . Thus, the Andronov-Hopf bifurcation occurs. For fixed h, f 1 and f 2 , the problem of local dynamics of (83) and (84) in the neighborhood of the zero equilibrium is solved in the same was as in Sections 2 and 3. We obtain the normal form with the same coefficient σ as in (10). Under the condition that the delay value h (83) is large enough, the problem becomes much more complicated. Let us show this below.
Let the equality be met for some fixed value h 1 > 0. We consider the characteristic equation for the linearized at zero boundary value problem (83) and (84): Under the condition (86), infinitely many of its roots tend to the imaginary axis as → 0.
In the case (86), the algorithm for the local dynamics study is the same as above [25,26]. We consider the formal series (11) for T 0 = 1. Then, as a result of standard actions, we obtain the normal form similar to (10), which is the equation with delay The main result is that according to the asymptotic formula (11), the nonlocal dynamics of the Equation (88) determines the local behavior of the boundary value problem (83) and (84) solution. An important distinction from (10) is not only in the appearance of the containing delay term. The expression (88) includes the parameter exp(iφ), which infinitely runs over all values on the unit circle on the complex plane as → 0. Therefore, we can only assert that the determined for some φ = φ 0 , some steady state in (88) corresponds to the solution of the boundary value problem (83) and (84) with asymptotic behavior (11) for the sufficiently small values of = n where n = πh 1 [2(φ 0 + 2πn)] −1 .
The dynamic behavior of (88) may differ for different values of φ. This means that an infinite process of forward and reverse bifurcations can occur in (83) and (84) as → 0. We note that interesting results were obtained on the Equation (88) dynamics under the condition h 1 1 in [25,27].

Conclusions
The questions regarding the behavior of all the solutions of the logistic equation with delay and diffusion in a small neighborhood of a positive equilibrium are considered. It is assumed that the most common and significant conditions are fulfilled when the critical Andronov-Hopf case is realized in the problem of the equilibrium stability. All coefficients differ from the constant values by some rather small function of the spatial or time variable. In particular, we note that the perturbations of both the delay coefficient and the boundary conditions coefficients are considered.
The normal forms are constructed. They are the special complex ordinary differential first-order equations on the central manifold, the nonlocal dynamics that determines the behavior of all the solutions of the original boudary value problem in a small neighborhood of the equilibrium. The asymptotic formulas that couple the solutions of the parabolic boundary value problem with the normal form solutions are presented.
In the case when the small time-periodic perturbation with the frequency close to the double that of natural oscillations affects the equation coefficients, the equation becomes more complicated on the central manifold. The problem of parametric resonance is considered as an application. In particular, it is shown that under certain conditions, the dynamics of the original boundary value problem can be controlled by means of a phase shift in a small periodic perturbation. The most complicated results concern the problem of parametric resonance with a two-frequency perturbation. The conditions under which the local dynamics of the problem under consideration can be complex and irregular are revealed with the help of numerical and analytical methods.