The performance of laser beams in tracking Lissajous scan trajectories is severely limited by beam jitter. To enhance the performance of fast steering mirror (FSM) control in tracking Lissajous scan trajectories, this paper proposed a fractional order active disturbance rejection controller (FOADRC) and verified its effectiveness in improving system scanning tracking accuracy. A dynamic mathematical model of fast steering mirror was studied and the design of parameters for the control mode of the closed-loop system was determined. A reduced order linear active disturbance rejection controller suitable for FSM systems was designed and the corresponding fractional order proportional differentiation (FOPD) controller were determined according to the mathematical model. The use of the designed controller enabled high-performance tracking of high-frequency Lissajous scanning curves (X-axis 500 Hz, Y-axis 350 Hz) and met the need for high-frequency repetitive scanning. The controller has the characteristics of simple implementation and low computational complexity and is suitable for closed-loop control applications in engineering.

The Stepped Frequency Waveform (SFW) is commonly employed in radar technology to synthesize wideband signals through the aggregation of narrow-band pulses, thereby achieving high-resolution range profiles without the need to increase the radar's instantaneous bandwidth. However, the inherent large time-bandwidth product of SFW introduces substantial ranging errors and energy dispersion, which significantly hinders its efficacy in detecting high-speed objects. This paper introduces a pioneering velocity estimation technique utilizing the fractional Fourier transform (FrFT) to address these limitations. Leveraging the characteristic of the Doppler signal from a moving target, which manifests as a chirp signal with a rate proportional to the target's velocity, the FrFT is utilized for precise velocity estimation. Subsequent to this, velocity compensation is applied using the deduced metrics, followed by the application of the inverse fast Fourier transform (iFFT) to pinpoint the target's exact location. To optimize the computational efficiency of determining the FrFT's optimal order, we propose an iterative algorithm founded on the golden section search method. The effectiveness of the proposed approach is verified by simulation data, and the results demonstrate that the proposed approach can accurately estimate the velocity and the range of the high-speed targets with a relatively low computational complexity.

The article outlines a creative approach by fusing ideas from fuzzy contractions, graph mappings, and Kannan mappings to create a brand-new notion known as "Kannan-graph-fuzzy contraction." This combination of concepts makes it possible to create a unique framework for resolving real-world issues controlled by nonlinear models. It is anticipated that using Kannan-graph-fuzzy contractions will offer a fresh approach to discussing solution existence in situations where these ideas overlap. This strategy might have a significant impact on a lot of real-world applications, proving how practical multidisciplinary techniques are for tackling issues in the real world originating from technology and engineering. The article includes several illustrated examples based on computer simulations, making the results more explicit.

The key objective of the present paper is to propose a numerical scheme based on the homotopy analysis transform technique to analyze the time-fractional nonlinear predator-prey population model. The population model is coupled fractional order nonlinear partial differential equations often employed to narrate the dynamics of biological systems in which two species interact, first is a predator and the second is a prey. The proposed scheme provides the series solution with great freedom and flexibility by choosing appropriate parameters. The convergence of results is free from small or large parameters. Three examples are discussed to demonstrate the correctness and efficiency of the used computational approach.

Analyzing the impacts of urban landscape patterns on the thermal environment has become one of the key research areas in addressing urban heat islands (UHI) and improving the living environment. A case study was carried out in Fuzhou, Fujian Province of China, to quantify the spatial patterns of urban landscape and land surface temperature (LST) in different development areas. The contribution of the thermal environment and the factors influencing the heterogeneity of LST were analyzed and discussed. In this study, bi-temporal Landsat imagery was selected to calculate LST, percent impervious surface area (ISA), and fractional vegetation cover (FVC). The urban area was further divided into three concentric urban zones, ranging from the city center to the urban periphery, based on urban development densities. The spatial pattern of LST and its variance were analyzed and compared between different zones and different dates. The thermal environment contribution of different zones was also quantified to indicate the change in urban landscape pattern resulting from urban expansion in different zones. Furthermore, Geodetector was used to explore the single factors and interaction factors controlling the spatial patterns of LST in each zone. The results showed that (ⅰ) Urban expansion primarily increased in zone 2 and zone 3, and the areal proportion of high and sub-high LST areas increased from 56.11% and 21.08% to 62.03% and 32.49% in zone 2 and zone 3, respectively, from 2004 to 2021; (ⅱ) The heat effect contribution of zones 2 and 3 reached from 75.16% in 2004 to 89.40% in 2021, indicating that the increase of ISA with >LSTmean was more pronounced in zone 3 and zone 2 during the period; (ⅲ) The driving factors of LST spatial distribution were regionally different because of the different landscape patterns, and the explanatory power for the heterogeneity of LST in zone 1 was weaker than that in zone 2 and zone 3 in the study area; (ⅳ) The interaction of different factors had a higher explanatory power on the spatial distribution of LST than a single factor in each zone because the distributions of land cover types are heterogeneous in urban area. The results of the study can be used to improve urban planning for urban ecology and UHI mitigation.

In this study, we first discover some new concept coordinated UD-convex mappings with fuzzy-number values. After that, we look into Hermite-Hadamard type inequalities via fuzzy-number-valued coordinated UD-convex fuzzy-number-valued mapping (coordinated UD-convex FNVM). In the case of coordinated UD-convex FNVM, novel conclusions are derived by making particular decisions in recently proven inequalities. Additionally, it is demonstrated that the recently discovered inequalities are expansions of comparable findings in the literature. It is important to note that the main outcomes are validated by nontrivial examples.

In this paper, existence and stability results of non-local integro-differential equation with Hilfer derivative of order ω∈(1,2) have been investigated. In this derivative is applied in RLC circuit equation of various fractional order. To proposed problem the existence solutions are derived using Schaefer’s fixed point theorem and the Banach contraction principle is using for uniqueness results. To demonstrate the effectiveness and applicability of our theoretical conclusions and two-step Lagrange polynomial interpolation were used to solve four numerical results.

The multivariate Gaussian random fields with matrix-based scaling laws are widely used for inference in statistics and many applied areas. In such contexts interests are often in symmetry and in the rates of change of spatial surfaces in any given direction. This article analyzes the almost sure sample function behavior for operator fractional Brownian motion, including multivariate fractional Brownian motion. We obtain the estimations of small ball probability and the strongly locally nondeterministic for operator fractional Brownian motion in any given direction. Applying these estimates we obtain Chung's laws of the iterated logarithm for spatial surfaces of operator fractional Brownian motion. Our results show that the precise rates of change of spatial surfaces are completely determined by the self-similarity exponent.

Four discrete models using the exact spectral derivative discretization finite difference (ESDDFD) method are proposed for a chaotic five-dimensional, conformable fractional derivative financial system incorporating ethics and market confidence. Since the system considered was recently studied using the conformable Euler finite difference (CEFD) method and found to be hyperchaotic, and the CEFD method was recently shown to be valid only at fractional index , the source of the hyperchaos is in question. Through numerical experiments, illustration is presented that the hyperchaos previously detected is in part an artifact of the CEFD method as it is absent from the ESDDFD models.

The fractional telegraph partial differential equation with fractional Atangana-Baleanu-Caputo (ABC) derivative is studied. Laplace method is used to find the exact solution of this equation. Stability inequalities are proved for the exact solution. Difference schemes for the implicit finite method are constructed. The implicit finite method is used to deal with modelling the fractional telegraph differential equation defined by Caputo fractional of Atangana-Baleanu (AB) derivative for different interval. Stability of difference schemes for this problem is proved by the matrix method. Numerical results with respect to the exact solution confirm the accuracy and effectiveness of the proposed method.

The cell concentration measurement on a culture assay using bioimpedance is a very useful tool, but it is complex to translate impedance to cell concentration values. The purpose of this work is to find a method to obtain in real time the cell concentration values for a given cell-culture assay using an oscillator as the measurement circuit. From a basic cell-electrode model, enhanced models of the cell culture immersed in a saline solution (culture medium) can be derived. These models can be used in a fitting routine to real time estimation of the cell concentration in a cell culture, using the oscillation frequency and amplitude delivered by measurement circuits proposed by the authors. Based on real experimental data (frequency and amplitude of oscillations) obtained by connecting the cell culture to an oscillator as a load, the fitting routine is simulated, and real time data of cell concentration is achieved. These results are compared with concentration data found by traditional optical methods for counting. In addition, the error obtained is divided and analyzed in two parts: in the first part of the experiment (when the few cells are adapting to the culture medium) and the second part of the experiment (when the cells grow exponentially until they completely cover the well).Low error values are obtained in the growth phase of the cell culture (the relevant phase), therefore the results obtained are considered to be promising, proving that the fitting routine is valid, and that the cell concentration can be measured in real time using an oscillator.

The evolution of integrable classical systems leads to conserved quantities and vanishing Poisson brackets. In contrast, such invariants do not exist in the dynamics of non-integrable systems, which include (but are not limited to) deterministic models with long-term chaotic behavior. The object of this review is to briefly survey the mathematical background of nonintegrability and its role in the physics unfolding well-above the Standard Model (SM) scale.

Suitable measures of grazing impacts on ground cover, that enable separation of the effects of climatic variations, are needed to inform land managers and policy makers across the arid rangelands of the Northern Territory of Australia. This work developed and tested a time-series, change-point detection method for application to time series of vegetation fractional cover derived from Landsat data to identify irregular and episodic ground-cover growth cycles. These cycles were classified to distinguish grazing impacts from that of climate variability. A measure of grazing impact was developed using a multivariate technique to quantify the rate and degree of ground cover change. The method was successful in detecting both long term (> 3 years) and short term (< 3 years) growth cycles. Growth cycle detection was assessed against rainfall surplus measures indicating a relationship with high rainfall periods. Ground cover change associated with grazing impacts was also assessed against field measurements of ground cover indicating a relationship between both field and remotely sensed ground cover. Cause and effects between grazing practices and ground cover resilience can now be explored in isolation to climatic drivers. This is important to the long term balance between ground cover utilisation and overall landscape function and resilience.

Androgenetic alopecia, the most common cause of hair loss affecting both men and women, is typically treated using pharmaceutical options such as minoxidil and finasteride. While these medications work for many individuals, they are not suitable options for all. To date, the only non-pharmaceutical option that the United States Food and Drug Administration has cleared as a treatment for androgenetic alopecia is low-level laser therapy (LLLT). Various clinical trials utilizing LLLT devices of various types are available. However, a myriad of other physical treatments for this form of hair loss are reported in the literature. This review evaluated the effectiveness of microneedling, pulsed-electromagnetic field (PEMF) therapy, low-level laser therapy (LLLT), fractional laser therapy, and non-ablative laser therapy for the treatment of androgenetic alopecia (AGA). It also explores the potential of multimodal treatments combining these physical therapies. LLLT has the majority of evidence in the literature as a physical therapy for androgenetic alopecia. However, other physical treatments, such as non-ablative laser treatments, and multimodal approaches, such as PEMF-LLLT, seem to have the potential to be equally or more promising and merit further exploration.

In this paper, the time-resolved model of photoacoustic signal for samples with a complex inner structure is derived including local non-equilibrium of structural elements with multiple degrees of freedom, i.e. structural entropy of the system. The local non-equilibrium is taken into account through the fractional operator. By analyzing the model for two types of time-dependent excitation, a very short pulse and a very long pulse, it is shown that the rates of non-equilibrium relaxations in complex samples can be measured by applying the derived model and time-domain measurements. Limitations of the model and further directions of its development are discussed.

In the present work we study the existence of solutions for hybrid nonlinear fractional integrodifferential equations. We developed an algorithm by using the operator theoretical techniques in order to obtain the approximate solutions. The main results depend on the Dhage iteration method that were incorporated with the modern hybrid fixed point theorems. The approximate solutions were obtained by using Lipschitz conditions and weaker form of mixed partial continuity. Further, we provide some examples to explain the hypotheses and the related results.

To obtain drop-in fuel properties from non-feed biomass, we herein report the catalytic hydrotreatment of microalgae biocrude, produced from hydrothermal liquefaction (HTL) of Spirulina. Our contribution focuses on the effect of temperature, initial H2 pressure, and residence time on the removal of heteroatoms (O and N). In contrast to common hydrotreating experimental protocols at batch scale, we devised a set of two-level factorial experiments and studied the most influential parameters affecting the removal of heteroatoms. It was found that up to 350 °C, the degree of deoxygenation (de-O) is mainly driven by temperature, whereas the degree of denitrogenation (de-N) also relies on initial H2 pressure and temperature-pressure interaction.Based on this, complete deoxygenation was obtained at mild operating conditions (350 °C), reaching a concurrent 47 % denitrogenation. Moreover, three optimized experiments are reported with 100 % removal of oxygen. In addition, the analysis by GC-MS and Sim-Dis gives insight to the fuel quality. The distribution of heteroatom N in lower (<340 °C) and higher (>340 °C) fractional cuts is studied by a fractional distillation unit following ASTM D-1160. Final results show that 63-68 % of nitrogen is concentrated in higher fractional cuts.

Antimicrobial resistance (AMR) has been a serious threat to human health, and combination therapy is proved to be an economic and effective strategy to fight the resistance. However, the abuse of drug combinations would conversely accelerate the spread of AMR. In our previous work, it had been concluded that the mutant selection indexes (SIs) of one agent against a specific bacterial strain are closely related to the proportions of two agents in a drug combination. To discover probable correlations, predictors and laws for further proposing feasible principles and schemes guiding the AMR-preventing practice, here three aspects were further explored. First, the power function (y=ax^b, a > 0) correlation between the SI (y) of one agent and the ratio value (x) of two agents in a drug combination was further established based on the mathematical and statistical analyses for those experimental data, and two rules a₁ × MIC₁ = a₂ × MIC₂ and b₁ + b₂ = -1 were discovered from both equations of y=a₁x^b1 and y=a₂x^b2 respectively for two agents in drug combinations. Simultaneously, it was found that one agent with larger MPC alone for drug combinations show greater potency for narrowing itself MSW and preventing the resistance. Second, a new concept as mutation-preventing selection index (MPSI) was proposed and used for evaluating the mutation-preventing potency difference of two agents in drug combinations, and the positive correlation between the MPSI and the mutant prevention concentration (MPC) or minimal inhibitory concentration (MIC) was subsequently established. Inspired by this, the significantly positive correlation, contrary to previous reports, between the MIC and the corresponding MPC of antimicrobial agents against pathogenic bacteria was established using one hundred and eighty-one of data pairs reported. These results together of above three aspects indicate that the MPCs in alone and combination are very important indexes for drug combinations to predict the mutation-preventing effects and the trajectories of collateral sensitivity, and while the MPC of an agent can be roughly calculated from its corresponding MIC. Subsequently, the former conclusion was further verified and improved by the antibiotic exposure to forty-three groups designed as different drug concentrations and various proportions. The results further proposed that the C/MPC for the agent with larger proportion in drug combinations can be considered as a predictor and is the key to judge whether the resistance and the collateral sensitivity occur to two agents. Based on these above correlations, laws, and their verification experiments, some principles were proposed, and a diagram of the mutation-preventing effects and the resistant trajectories for drug combinations with different concentrations and ratios of two agents was presented. Simultaneously, the reciprocal of MPC alone (1/MPC), proposed as the stress factors of two agents in drug combinations, together with their SI in combination, is the key to predict the mutation-preventing potency and control the trajectories of collateral sensitivity. Finally, a preliminary scheme for antimicrobial combinations preventing the AMR was further proposed for subsequent improvement research and clinic popularization, based on the above analyses and discussion. Moreover, some similar conclusions were speculated for triple or multiple drug combinations.

Background: Coronary artery disease (CAD) remains a significant global health concern, being one of the most important contributors to morbidity and mortality worldwide. Accurate assessment of coronary lesions is crucial for guiding appropriate interventions and optimizing patient outcomes. Patients with acute myocardial infarction represent a population which is rarely selected in studies of intracoronary hemodynamics. Fractional Flow Reserve (FFR) has emerged as a valuable physiological index for evaluating the functional significance of coronary lesions, although the data of intravascular hemodynamic evaluation in patients with acute myocardial infarction is scarce [1-3]. In recent years, advancements in intravascular imaging techniques, such as Optical Coherence Tomography (OCT), have provided detailed insights into coronary anatomy both in acute or elective settings. [4-7]. Methods: This article aims to explore the correlation between different non-culprit vessel OCT-parameters and FFR value in the patients with recent acute myocardial infarction. The study included 114 patients with acute myocardial infarction (ST elevation myocardial infarction and non ST elevation myocardial infarction)for whom primary PCI (percutaneous coronary intervention) was performed for the culprit lesion. The non-culprit vessels with intermediate coronary stenoses (157 lesions) were evaluated in a staged procedure (during the index hospitalization) with OCT, FFR and revascularization was performed accordingly. Minimal luminal area (MLA) and percent luminal area stenosis (PSA) were analyzed in relation with the FFR value. Results: There is good correlation between OCT-derived parameters and FFR values in patients with recent acute myocardial infarction. Minimal luminal area (MLA) and percent luminal stenosis area (PSA) were independent predictors of FFR< 0,80 in multivariate logistic regression analysis. A value of 2,08 mm2 for the MLA and 76,85 % for the PSA have proved to be the best cut-offs for predicting a FFR< 0,80. Conclusions: OCT parameters are valuable additional tools in current angiographic practice, including patients with acute myocardial infarction, with the ability to bring numerous direct imagistic and indirect hemodynamic information.

In this paper, new types of M-fractional wave solutions of mathematical physics model named as truncated M-fractional (1+1)-dimensional non-linear simplified Modified Camassa-Holm model are achieved by applying the modified simplest equation (MSE), Sardar sub-equation and generalized Kudryashov techniques. The gained solutions containing dark, bright, periodic and mixed wave solitons. Effect of fractional order derivative is also discussed. Achieved wave solitons are verified by Mathematica tool. Few of the gained wave solitons are also described through 2-dimensional, 3-dimensional and contour graphs through Mathematica tool. The gained solutions are helpful for the further development of concerned model. Finally, these techniques are simple, fruitful and effective to deal with non-linear FPDEs.

This work is a top-level summary of several contributions published in the last three decades. It makes the case that complex dynamics of nonlinear systems lies at the heart of foundational physics.

The main objective of this paper is to introduce k-fractional derivative operator by using the definition of k-beta function. We establish some results related to the newly defined fractional operator such as Mellin transform and relations to k-hypergeometric and k-Appell's functions. Also, we investigate the k-fractional derivative of k-Mittag-Leffler and Wright hypergeometric functions.

Complex modeling has received significant attention in recent years and is increasingly used to explain the statistical phenomenon with increasing and decreasing fluctuations such as the similarity or difference of spike protein charge patterns of coronaviruses. Different from the existing covariance or correlation coefficient methods in traditional integer dimension construction, this study proposes a simplified novel fractional dimension derivation with the exact Excel tool algorithm. It involves the fractional center moment extension to covariance, which ends up a complex covariance coefficient that is better than the Pearson correlation coefficient, in the sense that the nonlinearity relationship can be further depicted. The spike protein sequences of coronaviruses were obtained from the GenBank and GISAID database, including the coronaviruses from pangolin, bat, canine, swine (three variants), feline, tiger, SARS-CoV-1, MERS, and SARS-CoV-2 (including the strains of Wuhan, Beijing, New York, German, and UK variant B.1.1.7) were used as the representative examples in this study. By examining the values above and below the average/mean based on the positive and negative charge patterns of the amino acid residues of the spike proteins from coronaviruses, the proposed algorithm provides deep insights into the nonlinear evolving trends of spike proteins for understanding the viral evolution and identifying the protein characteristics associated with viral fatality. The calculation results demonstrate that the complex covariance coefficient analyzed by this algorithm is capable of distinguishing the subtle nonlinear differences in the spike protein charge patterns with reference to Wuhan strain SARS-CoV-2 for which the Pearson correlation coefficient may overlook. Our analysis reveals the unique convergent (positive correlative) to divergent (negative correlative) domain center positions of each virus. The convergent or conserved region may be critical to the viral stability or viability; while the divergent region is highly variable between coronaviruses suggesting high frequency of mutations in this region. The analyses show that the conserved center region of SARS-CoV-1 spike protein is located at amino acid residues 900, but shifted to the amino acid residues 700 in MERS spike protein, and then to amino acid residues 600 in SARS-COV-2 spike protein, indicating the evolvement of the coronaviruses. Interestingly, the conserved center region of the spike protein in SARS-COV-2 variant B.1.1.7 shifted back to amino acid residues 700, suggesting this variant is more virulent than the original SARS-COV-2 strain. Another important characteristic our study reveals is that the distance between the divergent mean and the maximal divergent point in each of the viruses (MERS>SARS-CoV-1>SARS-CoV-2) is proportional to viral fatality rate. This algorithm may help to understand and analyze the evolving trends and critical characteristics of SARS-COV-2 variants, other coronaviral proteins and viruses.

The Landau damping effect was observed in collisionless plasma, as a microscopic resonant mechanism between electromagnetic radiation and the collective modes. In this paper we demonstrate the occurrence of the Landau damping at macroscopic scale in the interaction between water waves and anharmonic lattice of magnetic buoys. By coupling the Navier-Stokes equations for incompressible fluid with the nonlinear dynamics of an anharmonic magnetic lattice we obtain a resonant transfer of momentum and energy between the two systems. The velocity of the flow is obtained in the Stokes approximation with Basset type of drag force. The dynamics of the buoys is calculated in the surfactant approximation for a specific frequency, then we use Fourier analysis to obtain the general time variable interaction. After involving an integral Dirichlet transform we obtain the time dependent expression of the drag force, the interaction waves-lattice with a new term in the form of a Caputo fractional derivative. We compare the results of the model with experiments performed in a wave tank with free floating magnetic buoys under the action of small amplitude gravitational waves. This configuration can be applied in studies for the attenuation with resonant damping of rogue waves, storms or tsunamis.

In this paper, some new properties of the fundamental solution to the multi-dimensional space- and time-fractional diffusion-wave equation are deduced. We start with the Mellin-Barnes representation of the fundamental solution that was derived in the previous publications of the author. The Mellin-Barnes integral is used to get two new representations of the fundamental solution in form of the Mellin convolution of the special functions of the Wright type. Moreover, some new closed form formulas for particular cases of the fundamental solution are derived. In particular, we solve an open problem of representation of the fundamental solution to the two-dimensional neutral-fractional diffusion-wave equation in terms of the known special functions.

In this paper, the new optical wave solutions of truncated M-fractional perturbed Kundu-Eckhaus model with full non-linearity are obtained by utilizing the expa function technique and modified extended tanh expansion function technique. The solutions are in the form of dark soliton, bright soliton, singular solitons and other form of solutions. The gained solutions are helpful for the further development of concerned model. The obtained results have also been presented graphically in both two-dimensional and three-dimensional formats to discuss the dynamical features as well as the parametric dependence of the constructed solutions. The study offers a highly spectacular and acceptable techniques to combine various intriguing wave demonstrations for more sophisticated models of the modern day.

In this paper, using the two integers that describe the stationary 2-dimensional wave and the charge quantization along with the balance between the Lorentz force and electrical force, we succeed in deriving the fractional quantum Hall effect and the integer quantum Hall effect; we find that the latter exists as a special case of the former. Moreover, using the derived expression describing the fractional quantum Hall effect, a relationship between the plateau in the resistivity of the sample and the applied magnetic field is obtained. The findings of this model agree well with experimental measurements. Because the two integers that describe the stationary 2-dimensional wave and the charge quantization along with the force balance have concrete physical meanings in this work, we could provide a clear picture of the origin of both the integer quantum Hall effect and the fractional quantum Hall effect.

Various factors, such as uncertainty of component parameters and uncertainty of sensor meas-urement values, contribute to the difficulty of fault detection in the landing gear retrac-tion/extension hydraulic system. In this paper, we introduce linear fractional transformation technology and uncertainty analysis theory for the construction of the diagnostic bond graph of the landing gear retraction/extension hydraulic system. In this way, interval analytical redundancy relations and fault signature matrix can be derived. Using the fault signature matrix, existing faults of the system can be preliminarily detected and isolated. Additionally, interval analytical re-dundancy relations can be used to detect system faults in detail, and cases analysis can be carried out to determine if the actuator is externally or internally leaky, and if the landing gear selector valve is reversing stuck. Compared to the traditional analytical redundancy relations, this method takes into account the negative factors of uncertainty; and compared to the traditional absolute diagnostic threshold, the interval diagnostic threshold is more accurate and sensitive.

In the face of the upcoming 30th anniversary of econophysics, we review our contributions and other related works on the modeling of the long–range memory phenomenon in physical, economic, and other social complex systems. Our group has shown that the long–range memory phenomenon can be reproduced using various Markov processes, such as point processes, stochastic differential equations and agent–based models. Reproduced well enough to match other statistical properties of the financial markets, such as return and trading activity distributions and first–passage time distributions. Research has lead us to question whether the observed long–range memory is a result of actual long–range memory process or just a consequence of non–linearity of Markov processes. As our most recent result we discuss the long–range memory of the order flow data in the financial markets and other social systems from the perspective of the fractional Lèvy stable motion. We test widely used long-range memory estimators on discrete fractional Lèvy stable motion represented by the ARFIMA sample series. Our newly obtained results seem indicate that new estimators of self–similarity and long–range memory for analyzing systems with non–Gaussian distributions have to be developed.

Creating self-learning algorithms, developing deep neural networks and improving other methods that "learn" for various areas of human activity is the main goal of the theory of machine learning. It helps to replace the human with a machine, aiming to increase the quality of production. The theory of artificial neural networks, which already have replaced the humans in problems of detection of moving objects, recognition of images or sounds, time series prediction, big data analysis and numerical methods remains the most dispersed branch of the theory of machine learning. Certainly, for each area of human activity it is necessary to select appropriate neural network architectures, methods of data processing and some novel tools from applied mathematics. But the universal problem for all these neural networks with specific data is the achieving the highest accuracy in short time. Such problem can be resolved by increasing sizes of architectures and improving data preprocessing, where the accuracy rises with the training time. But there is a possibility to increase the accuracy without time growing, applying optimization methods. In this survey we demonstrate existing optimization algorithms of all types, which can be used in neural networks. There are presented modifications of basic optimization algorithms, such as stochastic gradient descent, adaptive moment estimation, Newton and quasi-Newton optimization methods. But the most recent optimization algorithms are related to information geometry, for Fisher-Rao and Bregman metrics. This approach in optimization extended the theory of classic neural networks to quantum and complex-valued neural networks, due to geometric and probabilistic tools. There are provided applications of all introduced optimization algorithms, what delighted many kinds of neural networks, which can be improved by including any advanced approaches in minimization of the loss function. Afterwards, we demonstrated ways of developing optimization algorithms in further researches, engaging neural networks with progressive architectures. Classical gradient based optimizers can be replaced by fractional order, bilevel and, even, gradient free optimization methods. There is a possibility to add such analogues in graph, spiking, complex-valued, quantum and wavelet neural networks. Besides the usual problems of image recognition, time series prediction, object detection, there are many are other tasks for modern theory of machine learning, such as solving problem of quantum computations, partial differential and integro-differential equations, stochastic processes and Brownian motion, making decisions and computer algebra.