**Contents**


### **About the Editors**

### **Jorge M. Cruz-Duarte**

Jorge M. Cruz-Duarte (PhD) is a research professor (since 2021) in the Research Group on Advanced Artificial Intelligence at the Tecnologico de Monterrey (TEC) and a member of the ´ Mexican National System of Researchers (SNI-CONACyT), IEEE, the Mexican Mathematical Society (SMM), and the Mexican Academy in Computer Sciences (AMEXCOMP). Jorge received his B.Sc. and M.Sc. degrees in Electronic Engineering from the Universidad Industrial de Santander (UIS) and his Ph.D. in Electrical Engineering from the Universidad de Guanajuato (UGto). He was on a Short Research Stay at the University POLITEHNICA of Bucharest (UPB). From 2019 to 2021, he was on a Postdoctoral Stay at TEC in collaboration with the Chinese Academy of Sciences (CAS). His research interests include automatic algorithm design, heuristic methods, data science, optimization, mathematical methods, thermodynamics, digital signal processing, electronic thermal management, and fractional calculus.

### **Porfirio Toledo-Hern´andez**

Porfirio Toledo-Hernandez (Ph.D.) is a professor in the Faculty of Mathematics at the University ´ of Veracruz. He is a member of the Mexican National System of Researchers (SNI-CONACyT) and the Mexican Mathematical Society (SMM). Porfirio obtained his Ph.D. in basic mathematics at the Center for Research in Mathematics (CIMAT) in 2011. He has developed research work in mathematical modeling and optimization in topics such as game theory, differential equations, dynamical systems, biological models, and behavior analysis.

### **Preface to "The 5th Mexican Workshop on Fractional Calculus"**

The Mexican Workshop on Fractional Calculus (MWFC), which originated in the Irapuato-Salamanca Campus Engineering Division of the University of Guanajuato in 2012, has an international scope. Its main objective is to present the latest theoretical advances and applications of fractional calculus to different areas of Science and Engineering. It is also an opportunity for free scientific discussion among young people and experts. Indisputable personalities in the area of fractional calculus have participated in these workshops and have stimulated young scientists. The workshops are aimed at high-level and postgraduate students, as well as young researchers and consolidated researchers. These workshops are held every two years at different universities around the country.

> **Prefaced by Prof. J. Juan Rosales-Garc´ıa Department of Electrical Engineering, DICIS-UG**

### *Editorial* **Fractional Calculus in Mexico: The 5th Mexican Workshop on Fractional Calculus (MWFC) †**

**Jorge M. Cruz-Duarte 1,\* and Porfirio Toledo-Hernández 2**


The Mexican Workshop on Fractional Calculus (MWFC) is a bi-annual international workshop and the largest Latin American technical event in the field of fractional calculus in Mexico. The focus of MWFC is to motivate critical discussions of ideas about abstract and non-integer operators, chaotic systems, complex behavior, and their applications in several areas. MWFC offers an ideal forum for intellectuals and curious people to present their latest research insights. In its 5th version, MWFC was held in Monterrey, one of the most influential cities in México and home of the Tecnológico de Monterrey (TEC), where MWFC was hosted in a hybrid format. During the 5th MWFC, we attended six plenary talks from distinguished international scientists, as well as two tutorials, a poster session, and four panels of accepted papers.

In submitting conference proceedings to *Computer Sciences and Mathematics Forum*, we certify to the publisher that all papers published in this volume have been subjected to peer review administered by the volume editors. Expert referees conducted reviews to the professional and scientific standards expected of a proceedings journal.


**Conflicts of Interest:** The authors declare no conflict of interest.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

**Citation:** Cruz-Duarte, J.M.;

Toledo-Hernández, P. Fractional Calculus in Mexico: The 5th Mexican Workshop on Fractional Calculus (MWFC). *Comput. Sci. Math. Forum* **2022**, *4*, 4007. https://doi.org/10.3390/ cmsf2022004007

 Published: 3 February 2023

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

•

1

### *Proceeding Paper* **Further Remarks on Irrational Systems and Their Applications †**

**Adrián-Josué Guel-Cortez 1,\*, César-Fernando Méndez-Barrios 2, Diego Torres-García 2 and Liliana Félix 2**


**Abstract:** Irrational Systems (ISs) are transfer functions that include terms with irrational exponents. Since such systems are ubiquitous and can be seen when solving partial differential equations, fractional-order differential equations, or non-linear differential equations; their nature seems to be strongly linked with a low-order description of distributed parameter systems. This makes ISs an appealing option for model-reduction applications and controls. In this work, we review some of the fundamental concepts behind a set of ISs that are of core importance in their stability analysis and control design. Specifically, we introduce the notion of multivalued functions, branch points, time response, and stability regions, as well as some practical applications where these systems can be encountered. The theory is accompanied by some numerical examples or simulations.

**Keywords:** irrational systems; fractional-order control; model-reduction methods

### **1. Introduction**

Irrational Systems (ISs) can be found when solving partial differential equations, fractional-order differential equations, or non-linear differential equations [1]. ISs have also been called implicit operators since they often come from solving a second-order polynomial whose solution describes the total impedance of infinite linear lumped-element networks (for further details, see [2]). In addition, fractional behavior or non-exponential decay can also be associated with ISs due to their time-response, often related to special functions such as Bessel or error functions [3]. Hence, ISs' nature is strongly linked to a low-order description of distributed parameter systems, making ISs an appealing option for control and model-reduction applications.

In the literature, ISs have been applied as mathematical models in different scenarios. For instance, ref. [4] uses fractances in a lumped model of the cardiovascular system, leading to the description of different types of heart anomalous behaviors. In [3,5,6], ISs are introduced to describe robotic formations in the form of tree-like networks and ladder-like networks. Ref. [7] models pipeline-infinite networks that converge to an IS. Similarly, ref. [8] shows an application of ISs to model electrical line transmissions.

When talking about control applications, in our previous work [9], we discuss the stability and control for a type of IS driven by fractional-order controls. Then, in [10], we formalize our analysis and apply the fractional-order controls of the type proportionalintegral (PI) and proportional-derivative (PD), including their fragility analysis.

As the conceptualisation and application of ISs to engineering are still under development, in this work, we summarize the fundamental concepts, assumptions, and limits of ISs. The summary aims to briefly explain IS stability analysis and control based on our previous and undergoing investigations. The work is organized as follows: Section 2 defines the fundamental concepts and details the origins of ISs. It also discusses the intrinsic connection between fractional calculus and ISs. Furthermore, it states the fundamental hypothesis when applying ISs in system modeling. Section 3 describes the conditions of

**Citation:** Guel-Cortez, A.-J.; Méndez-Barrios, C.-F.; Torres-García, D.; Félix, L. Further Remarks on Irrational Systems and Their Applications. *Comput. Sci. Math. Forum* **2022**, *4*, 5. https://doi.org/ 10.3390/cmsf2022004005

Academic Editors: Jorge M. Cruz-Duarte and Porfirio Toledo-Hernández

Published: 22 December 2022

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

stability for ISs. Section 4 briefly explains how to design low-order controllers of the type PD-*μ* (i.e., fractional-order proportional-derivative controls) for the kind of ISs that the paper focuses on. Section 5 gives some examples of ISs' control. Finally, Section 6 provides concluding remarks and outlines future work.

### **2. Preliminaries**

Before formally describing the set of ISs that this work will discuss, we need to introduce the following concepts:

**Definition 1** (Multivalued function [11,12])**.** *A multivalued function is any complex function F* : C → C *satisfying*

$$F[z(r, \theta + 2\pi)] \neq F[z(r, \theta)],\tag{1}$$

*where r and θ are the magnitude and argument of the complex variable z* ∈ C*.*

**Definition 2** (Branch points and branch cuts [13])**.** *The Branch point (BP) or point of accumulation is defined as the point with the smallest magnitude for which a function is multivalued. Another definition would be: a branch point is a point such that the function is discontinuous when going around an arbitrarily small circuit around this point.*

It is now easier to understand the following definition of an IS:

**Definition 3** (Irrational system)**.** *An irrational system is a multi-valued transfer function G*(*s*) *with one or more terms raised to the power α* ∈ Q*.*

In this work, we consider ISs described by

$$G(\mathbf{s}) = \frac{N(\mathbf{s}) + \sqrt{P(\mathbf{s})}}{D(\mathbf{s}) + \sqrt{Q(\mathbf{s})}'} \tag{2}$$

where *N*(*s*) = ∑*mk*=<sup>0</sup> *bks<sup>k</sup>*, *D*(*s*) = ∑*nk*=<sup>0</sup> *aksk*, *ai*, *bi*, *an* = 0 are arbitrary real numbers, and *n* ≥ *m*. Moreover, *P*(*s*) and *Q*(*s*) are second-order polynomials with positive real coefficients defined as *P*(*s*) = ∑*Nk*=<sup>0</sup> *βks<sup>N</sup>* and *Q*(*s*) = ∑*Nk*=<sup>0</sup> *<sup>γ</sup>ksN*, respectively.

Clearly, according to Definition 3, the open-loop system (2) and IS with *α* = 12 , and branch points are located at the roots of *P*(*s*) and *Q*(*s*).

### *Origins and Connection with Fractional Calculus*

Irrational transfer functions such as (2) are proven to appear when modeling infinite networks of lumped elements (for a complete description, see [1,7]). The proof of such a statement, omitted here for brevity, requires the following assumptions:



**Figure 1.** Tree-like network of *N* layers that can be described by an ISs transfer function. In the network, it is necessary to have L∞ and L∈ to be linear operators. Note that all end-points *xout* are in the same position. The movement is in one-dimension.

In spite of such necessary conditions, ISs can still be applied in many scenarios, epecially when the following hypotheses are considered:

**Hypothesis 1.** *The accuracy of 0-dimensional models of complex dynamical systems can increase by adding "networks" as lumped elements in the model.*

**Hypothesis 2.** *ISs can be used as model reductions for large-scale dynamical systems.*

The validity of Hypothesis 1 has been proved in scenarios such as the cardiovascular system [4,14] or muscle/joint modeling [15].

On the other hand, Hypothesis 2 has been explored in robot formations and transmission lines (for further details, see [6,16–18]). Furthermore, in some cases the ISs description leads to basic fractional-order transfer functions (for instance, see [18]).

For a better understanding of Hypotheses 1 and 2, consider the circuit-like description of the cardiovascular system shown in Figure 2a and the graphical description of a formation of mobile robots in one dimension in Figure 2b. In the case of Figure 2a, we can assume that one of the elements is itself an infinite network of linear elements with an impedance equal to 1 *<sup>C</sup><sup>α</sup>Fs<sup>α</sup>* (i.e., of fractional-order); then, its mathematical model is given by

$${}\_{0}\mathcal{D}\_{t}^{a}P\_{a} = \frac{Q\_{a}}{\mathcal{C}\_{F}^{a}} - \frac{P\_{a}}{R\mathcal{C}\_{F}^{a}}.\tag{3}$$

Furthermore, by using the Caputo definition of the fractional derivative operator <sup>0</sup>D*<sup>α</sup>t* of order 0 < *α* < 1, the time response of system (3) is given by

$$P\_a(t) = P\_a(0)t^{a-1} \mathbf{E}\_{a,a}(-\frac{1}{RC\_F^a}t^a) + \frac{1}{C\_F^a} \int\_0^t Q\_a(t-\tau)\tau^{a-1} \mathbf{E}\_{a,a}(-\frac{1}{RC\_F^a}\tau^a)d\tau,\tag{4}$$

where **<sup>E</sup>***<sup>α</sup>*,*<sup>α</sup>*(*z*) is the Mittag–Leffler function of the complex value *z* [19]. Equation (4) describes the arterial pressure as a power law equation with a diffusive term. Even though this kind of representation implies some physical challenges that have been recently discussed through various works (for instance, see [20,21] and the references therein), it

is clear that we can substitute any IS's impedance instead of a simple lumped element to expand the capabilities of 0-dimensional models, as discussed in Hypothesis 1.

 **Figure 2.** Examples of the application of Hypothesis 1 and 2 in realistic scenarios. (**a**) Model reduction of the cardiovascular system by an electrical system using a fractance. (**b**) Ladder network description of mobile robots described by mechanical elements and driven by PID controls [6].

Regarding Hypothesis 2, Figure 2b shows a PID-driven robotic formation in onedimension (for example, see [22]) where all robots are described as simple mechanical elements. If we regard the network as infinite, the transfer function relating the leader's position and the last robot in the formation would be an IS [6]. Is this representation a good approximation of the actual transfer function when the number of robots in the network is finite? The answer is still not conclusive and requires further investigation (a preliminary analysis is given in [3]).

### **3. Stability Analysis**

One of the grea<sup>t</sup> advantages of ISs mathematical models is that we can easily study their stability if the following Theorem is considered:

**Theorem 1** ([8])**.** *A given multivalued transfer function is stable if and only if it has no pole in* C+ *and no branch points in* C−*. Here,* C+ *and* C− *stand for the closed right half plane (RHP) and the open RHP of the first Riemann sheet in the complex plane, respectively.*

Briefly, Theorem 1 states that the BPs should not be located in the right-hand side of the complex plane to achieve the IS's stability. To grasp this conclusion, let us consider the following simple examples:

**Example 1.** *Consider the multivalued transfer function*

$$G(s) = \frac{1}{\sqrt{s+k}'} \tag{5}$$

*where k* ∈ R*. The impulse response of* (5) *is given by (proof given in Appendix A)*

$$y(t) = \mathcal{E}^{\epsilon - 1} \left[ \frac{1}{\sqrt{s + k}} \right] = \begin{cases} \frac{\epsilon^{-kt}}{\sqrt{\pi t}} & k > 0 \\ \frac{\epsilon^{kt}}{\sqrt{\pi t}} & k < 0 \\ \frac{1}{\sqrt{\pi t}} & k = 0 \end{cases} \tag{6}$$

In Example 1, the IS (5) is stable iff *k* > 0. Note that as mentioned earlier, the BP is located at the root of the radicand *s* + *k*.

**Example 2.** *Consider the multivalued transfer function*

$$G(s) = \frac{1}{\sqrt{s^2 + k}},\tag{7}$$

*where k* ∈ R*. The impulse response of* (7) *is given by (see Appendix B)*

$$y(t) = \mathcal{X}^{-1}\left[\frac{1}{\sqrt{s^2 + k}}\right] = \begin{cases} l\_0(\sqrt{k}t) & k > 0\\ l\_0(t\sqrt{k}t) & k < 0\\ 1 & k = 0 \end{cases} \tag{8}$$

Equation (7) is an example where the open-loop IS is stable as the BPs are located in the imaginary-axis and not in the right-hand side of the complex plane.

### **4. Control Design**

*PDμ Control*

Once the stability conditions for ISs are established, it is possible to design stabilising low-order controllers by following the *D*-composition method [23]. The following summarizes the procedure, but a complete guide can be found at references [10,24,25].

First, consider a low-order controller, for instance, the fractional-order PD control whose transfer function is well-known to be [9]

$$
\mathcal{L}(s) = k\_{\mathcal{P}} + k\_d s^\mu. \tag{9}
$$

Then, compute the characteristic polynomial of (2). In this case defined as

$$
\Delta(\mathbf{s}) := D(\mathbf{s}) + \sqrt{Q(\mathbf{s})} + \left(N(\mathbf{s}) + \sqrt{P(\mathbf{s})}\right) \left(k\_p + k\_d \mathbf{s}^\mu\right). \tag{10}
$$

Now, substitute the stability boundary of the complex plane *s* = *iω* in the characteristic polynomial and solve for the control gains.

**Remark 1.** *Note that the substitution of s* = *iω requires one to consider the cases ω* = 0*, ω* → ∞ *and ω* = 0 *separately [25]. In this manner, we create various sets of stability boundaries in the controller's parameter plane that permit us to enclose a stability region (See Figure 3).*

**Figure 3.** Example of the D-composition method. The method maps the complex plane stability region to the controller parameters' plane. In this case, the plane has not stability boundary at *s* → ∞.

### **5. Applications**

### *5.1. Control of IS*

We now apply the D-composition method to some specific ISs. Note that, in all examples, the time response of the ISs closed-loop is obtained numerically using the inverse numerical Laplace transform proposed in [26]. Furthermore, fixed parameters have been selected randomly as the example's purpose is to prove the results' validity.

### *5.2. Bessel*

Consider the Laplace transform of the Bessel function of order zero described as

$$\frac{1}{\sqrt{s^2+1}}.\tag{11}$$

Therefore, the characteristic function is given by

$$
\Delta(\mathbf{s}) = \sqrt{s^2 + 1} + k\_p + k\_d \mathbf{s}^\mu. \tag{12}
$$

Setting *μ* = 0.3 and solving (12) for *kp* and *kd*, we obtain two stability boundaries shown in red and blue in Figure 4a. The red curve corresponds to the case where *ω* = 0 before solving Δ(*iω*) = 0, while the blue corresponds to the case where *ω* = 0. Figure 4b shows the closed-loop response of a set of controller gains inside and outside of what is found to be the IS's stability region.

**Figure 4.** Stability analysis of system (11). (**a**) Stability region (gray) of the closed-loop system with *μ*= 0.3. (**b**) Time response for control gains inside different regions on the parameter's plane.

### *5.3. First Order IS*

Another example of an IS that can be controlled by a fractional-order PD control is

$$\frac{\sqrt{3s+1}}{s+\sqrt{2s+1}}.\tag{13}$$

In this case, the characteristic equation is

$$
\Delta(\mathbf{s}) = \mathbf{s} + \sqrt{2\mathbf{s} + 1} + \sqrt{3\mathbf{s} + 1}. \tag{14}
$$

Likewise, in the previous application we set *μ* = 0.4 and solved (14) for *kp* and *kd* to obtain the stability boundaries shown in Figure 5a. Again, the red curve corresponds to the case where *ω* = 0 before solving Δ(*iω*) = 0, while the blue corresponds to the case where

*ω* = 0. Figure 5b shows the closed-loop response of a set of controller gains inside and outside of what is found to be the IS's stability region.

### **6. Conclusions**

In this work, we have briefly described the concepts, hypotheses, and assumptions behind the use of ISs. In addition, in our applications, we present the stability analysis and control of ISs. As can be seen by the reader, future work may take various paths, including the design of different low-order controls for ISs, the application of ISs to model other complex phenomena, or the creation of irrational controls that could stabilize ISs whose BPs are in the right-hand side of the complex plane.

**Author Contributions:** Conceptualization, A.-J.G.-C. and C.-F.M.-B.; methodology, A.-J.G.-C.; software, A.-J.G.-C. and D.T.-G.; validation, A.-J.G.-C., C.-F.M.-B., and D.T.-G.; formal analysis, A.-J.G.-C. and C.-F.M.-B.; investigation, A.-J.G.-C.; resources, A.-J.G.-C. and L.F.; data curation, A.-J.G.-C., D.T.-G., and L.F.; writing—original draft preparation, A.-J.G.-C. and D.T.-G.; writing—review and editing, C.-F.M.-B. and L.F.; visualization, A.-J.G.-C.; supervision, A.-J.G.-C. and C.-F.M.-B.; project administration, A.-J.G.-C. and C.-F.M.-B.; and funding acquisition, L.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Faculty of Engineering, Autonomous University of San Luis Potosi gran<sup>t</sup> number POA-AME-IMT2022.

**Institutional Review Board Statement:** Not applicable.

**InformedConsentStatement:**Notapplicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.
