Том 30, № 2 (2022)

Purpose of this work is to study the initial-boundary value problem for a parabolic functional-differential equation in an annular region, which describes the dynamics of phase modulation of a light wave passing through a thin layer of a nonlinear Kerr-type medium in an optical system with a feedback loop, with a rotation transformation (corresponds the involution operator) and the Neumann conditions on the boundary in the class of periodic functions. A more detailed study is made of spatially inhomogeneous stationary solutions bifurcating from a spatially homogeneous stationary solution as a result of a bifurcation of the “fork” type and time-periodic solutions of the “traveling wave” type. Methods. To represent the original equation in the form of nonlinear integral equations, the Green’s function is used. The method of central manifolds is used to prove the theorem on the existence of solutions of the indicated equation in a neighborhood of the bifurcation parameter and to study their asymptotic form. Numerical modeling of spatially inhomogeneous solutions and traveling waves was carried out using the Galerkin method. Results. Integral representations of the considered problem are obtained depending on the form of the linearized operator. Using the method of central manifolds, a theorem on the existence and asymptotic form of solutions of the initial-boundary value problem for a functional-differential equation of parabolic type with an involution operator on an annulus is proved. As a result of numerical modeling based on Galerkin approximations, in the problem under consideration, approximate spatially inhomogeneous stationary solutions and time-periodic solutions of the traveling wave type are constructed. Conclusion. The proposed scheme is applicable not only to involutive rotation operators and Neumann conditions on the boundary of the ring, but also to other boundary conditions and circular domains. The representation of the initial-boundary value problem in the form of nonlinear integral equations of the second kind allows one to more simply find the coefficients of asymptotic expansions, prove existence and uniqueness theorems, and also use a different number of expansion coefficients of the nonlinear component in the right-hand side of the original equation in the neighborhood of the selected solution (for example, stationary). Visualization of the numerical solution confirms the theoretical calculations and shows the possibility of forming complex phase structures.

The objective of this study is to show the possibility of using the smoothing cardiointervalograms (CIG) method which is solely time domain analysis of CIG to separate and display the influence of various mechanisms of human physiological regulation systems on his heart rate. Methods.This paper shows the possibility of using the method of smoothing the cardiointervalogram by means of a moving average for its subsequent decomposition into slow and fast components. Decomposition results are visualized by line graphs and pseudo-phase portraits. Visualization settings allow us to isolate unique transients and calculate its timing. The method is applied to data obtained under different subject functional states and differing in the level of adaptation risks, the presence or absence of stress. For analysis were selected stress episodes detected using the information and telecommunication technology of event-related cardiac telemetry (ITT ERCT) presented by the Internet resource “StressMonitor”. Results.For the numerical series of RR-intervals, a clear division into fast and slow components is obtained. An algorithm for identifying the frequency content of heart rate variability has been formulated and tested. A visualization method is proposed that is convenient for comparing data obtained for different patients. A pseudo-phase portrait pattern corresponding to the moment of stress onset is found. The proposed method reduced the discreteness of identifying the stress onset moment from 10 seconds to single heart beats. Conclusion. The correspondence of the results to the verified ITT ERCT method and the Baevsky–Chernikova concept of adaptive risk has been demonstrated. This confirms the possibility of using the time cardiointervalograms smoothing method for the analysis of heart rate variability.

Purpose is to study the mechanisms leading to genetic divergence (stable genetic differences between two adjacent populations). We considered the following classical model situation. Populations are panmictic with Mendelian rules of inheritance. The action of natural selection (differences in fitness) on each of population is the same and is determined by the genotypes of only one diallel locus. We assume that adjacent generations do not overlap and genetic transformations can be described by a discrete time model. This model describes the change in the concentration of one of the alleles in each population and the ratio (weight) of first population to the total size. Methods. We used the analogue of saddle charts to construct parametric portraits showing the domains of qualitatively different dynamic modes. The study is supplemented with phase portraits, basins of attraction and bifurcation diagrams. Results. We found that the model dynamic regimes qualitatively coincide with the regimes of a similar model with continuous time, but only for a weak migration. With a strong coupling, fluctuations of the phase variables are possible. We showed that the genetic divergence is possible only with reduced fitness of heterozygotes and is the result of a series of bifurcations: pitchfork bifurcation, period doubling, or saddle-node bifurcation. After these qualitative changes, the dynamics become bi- or quadstable. In the first case, the solutions corresponding to the genetic divergence are unstable and are just a part of the transient process to monomorphic state. In the second case, the divergence is stable and appears as 2-cycle for a strong migration coupling. Conclusion. In neighboring populations, movement towards an asymptotic genetic structure (monomorphism, polymorphism or divergence) can be strictly monotonous or in the form of damped unstable or undamped stable fluctuations with a period of 2 for biologically significant parameters. For insignificant parameters, we found a complex dynamics (chaos) that consist of divergent fluctuations around fixed points and quasi-random transitions between them.

Purpose of this study introduces a multielectrode array (MEA) registration system in order to generate electric field images of the episodic discharges generated by weakly electric fish. A multielectrode registration system has several important features: the design of the multielectrode lattice, the amplifier circuit, the choice of reference points for differential measurements, the recovery of the absolute values of the electric field potentials, and the application of principal components analysis. Methods. There are several advantages of our MEA registration as compared with the traditional twoelectrode registration: (a) the signal-to-noise ratio is significantly increased, (b) it is possible to construct the spatial distribution of the electric field for a single electric discharge, (c) the signals’ sources can be easily separated and identified, and (d) quantitative data on the electrical potential distribution can be obtained throughout the entire experimental tank. Results. The results illustrate an example of applied MEA registration. Electric discharges were recorded from a weakly electric catfish, Clarias gariepinus, using an array of 8 x 8 electrodes at a sampling rate of 20 kHz. Data show oscillograms and two-dimensional plots of the spatial distribution of the electrical field.