Abstract

This paper presents the design, simulation, and experimental verification of the fractional-order multiscroll Lü chaotic system. We base them on op-amp-based approximations of fractional-order integrators and saturated series of nonlinear functions. The integrators are first-order active realizations tuned to reduce the inaccuracy of the frequency response. By an exponential curve fitting, we got a convenient design equation for realizing fractional-order integrators of orders from 0.1 to 0.95. The results include simulations in SPICE of the mathematical description and the electronic implementation and experimental measurements that confirm them. Monte Carlo and sensitivity tests revealed a robust realization. Contrary to its passive counterparts, the suggested realizations significantly reduce design and implementation efforts by favoring resistors and capacitors with commercial values and reducing hardware requirements.

1. Introduction

Fractional calculus has been known since 1695 when L’Hôpital and Leibniz interchange letters about the noninteger order of the derivative. The fractional calculus is considered as a generalization of the integer-order version. Fractional-order derivatives are characterized by their memory properties providing an excellent approach for describing the physical phenomena elegantly with improved accuracy [1, 2]. As a result, the researchers have proposed novel mathematical models with applications in various fields, such as biology, economics, chaos theory, botany, hidden dynamics, digital circuits, cryptography, control, image processing, wind turbines, viscoelastic studies, and ferroelectric materials [315].

Regarding chaos, there has been recently an increased number of works reporting fractional-order chaotic systems because fractional order provides an extra degree of freedom, which can be useful for generating diverse dynamics such as self-excited attractors, hidden attractors, multistability, and extreme stability [1620]. In particular, the multiscroll chaotic attractors present plenty of complex topological structures contrary to single- or double-scroll attractors. One of the most typical applications is secure communications, e.g., in random number generators, cryptosystems in wireless networks, and image encryption [2123]. Therefore, it is a subject of increasing interest [2430].

The physical realization of chaotic systems is a well-accepted technique to observe the chaotic behaviors and therefore validates the theoretical findings. From an engineering point of view, the challenge consists of implementing the fractional-order multiscroll chaotic systems physically since there is not a fabricated fractional-order capacitor yet commercially available [31, 32]. Hence, in the literature, we can find two approaches for implementing fractional-order differential equations that describe the chaotic systems. On the one side, the embedded methodology uses FPGAs or DSPs to describe the fractional-order dynamical system by implementing a numerical method to solve computationally the system rather than an electronic circuit [33]. Nevertheless, this approach has been demonstrated that induces degradation effects corrupting the chaos generation. On the other side, the analog method employs essential circuit elements for obtaining the fractional-order integrators. The basic idea consists of representing the fractional order with an integer-order transfer function of high degree. As a result, the electronic realization of the fractional integrator depends on a tree-like passive network formed with resistors and capacitors. However, the resulting passive emulator circuits nor contain commercial values neither a suitable design equation between elements (ratio of values), provoking a complexity for a proper tuning of the fractional order for observing chaos behavior. Therefore, the works focused on the physical implementation of fractional-order chaotic systems generating multiscrolls are limited and only found with implementations using passive networks [2428].

The question arising is how to obtain a straightforward approach for implementing fractional-order multiscroll chaotic systems with analog topologies? This paper will give a positive answer to this question. A first response is to use fractional-order integrators derived from an active circuit topology, which can be designed with circuit elements of commercial values or using an appropriate ratio of values. In this regard, some first results were proposed in refs. [31, 32], which reported active realization using integrated circuit technology. However, that approach has not been proved in either single-scroll or multiscroll chaotic systems. Therefore, the main contributions of this paper are as follows: (i) an optimized design methodology is proposed based on continued fraction expansion to obtain an active fractional integrator which can be related in a quick way to the fractional order of the chaotic systems; (ii) to the best knowledge of authors, the first experimental validation of a fractional-order multiscroll chaotic system is proposed based on active circuit topologies; and (iii) additionally, we perform an exhaustive study to demonstrate that the proposed design methodology is robust to statistical variations of element values and insensible to changes in the circuit elements.

The electrical circuits to realize fractional-order integrators are based on “fractances,” electrical elements with fractional-order impedance. These elements produce a magnitude frequency response with roll-off  dB/dec and constant-phase response of degrees [34]. Unfortunately, the nonexistence of circuit elements with fractional response challenges carries out fractional-order transfer functions. The alternative is the approaching of fractances in a desired bandwidth with rational functions by means of methods of approximation such as Muir, Matsuda, Newton, Carlson, Oustaloup, power series (PSE), and continuous fraction expansion (CFE) [3436]. The larger the order of the approximation, the larger the bandwidth where it is valid. Once got the rational function, the synthesis is possible using ladder networks of Cauer or Foster [36], tree structures, or transmission lines [35, 36]. The disadvantage of these realizations is the difficulty to approximate the required elements with commercially available resistances, capacitors, inductors, and even negative impedance converters, which results in bulky realizations [3639]. For instance, Figure 1 shows a third-order realization using the Charef approach for a fractional-order passive fractance with , an error of 1 dB, and a bandwidth from 0.1 rad/s to 1000 rad/s [40]. The circuit requires three resistors and three capacitors of disperse and not available values. To avoid this drawback in this paper, we present two proposals of implementation of with active elements based on algebraic manipulation of a first-order CFE approximation. Despite its inaccuracy, in this work, we explore the use of a compact first-order approximation with obtained by the CFE approach, but tuned to reduce error in the magnitude’s slope response. It has the form [34]

This article is organized as follows. Section 2 provides a brief survey of fractional calculus, approximation of fractances, and op-amp-based building blocks. Section 3 depicts the circuit realization of the fractional-order integrator and presents Monte Carlo and sensitivity tests. Section 4 offers simulation and experimental results of the Lü chaotic oscillator and discusses the advantages of the proposed realization. We draw the conclusion in Section 5.

2. Theoretical Background

This section describes theoretical background such as fractional-order calculus, the Laplace transform of a fractional derivative, and op-amp-based building blocks for implementing fractional-order integrators and chaotic oscillators.

2.1. Fractional Derivative and Fractional Laplacian Operator

The Riemann–Liouville definition for calculation of fractional derivatives and integrals of a function establishes [41, 42]where is the fractional order, , , and is the gamma function. For , , and , one gets the fractional derivative, fractional integral, and identity function. The Laplace transform of (2) with zero initial conditions becomes [41]where stands for Laplace transform of and is the fractional Laplacian operator given by

2.2. Op-Amp-Based Building Blocks

The circuits of Figure 2 are the analog electronic blocks used in this paper to realize fractional-order integrators and chaotic oscillators. The selection of , , , , and as degrees of freedom will allow reducing the number of passive elements with disperse or not commercially available values.

2.2.1. Weighted Differential Amplifier (WDA)

The op-amp-based amplifier of Figure 2(a) amplifies the weighted difference between two voltages and with pondering factors depending on and . By nodal analysis, the output voltage can be calculated as

2.2.2. Weighted Adder Amplifier (WAA)

The amplifier of Figure 2(b) provides the weighted sum of two voltages and given bywhere and control the pondering factors. Note that since this circuit is a noninverter adder based on a voltage series feedback amplifier, it presents a negative feedback loop which adds up stability to the weighted coefficients of and , which in turn are the pondering factors and .

2.2.3. Inverting Integrator (IInv)

The circuit of Figure 2(c) uses capacitive feedback to integrate the voltage through

A capacitor can replace to translate the frequency response by a factor , where

In the time domain, expression (7) is rearranged to be written as

2.2.4. Inverting Amplifier (Inv)

The circuit of Figure 2(d) provides a negative voltage gain given by

This negative gain is indicative of an inverting amplifier with a phase shift between and . The gain can be greater than, less than, or equal to 1 depending on the values of and .

2.2.5. Inverting Weighted Adder Amplifier (IWAA)

The amplifier of Figure 2(e) provides a weighted negative sum of four input voltages , , , and with pondering factors , , , and . By nodal analysis, the output voltage can be calculated as

2.2.6. PWL Function Emulator (PWL)

Multiscroll chaotic oscillator circuits use saturated nonlinear functions (SNLFs) approximated by piecewise linear (PWL) functions [43, 44]. Figure 2(f) shows the SNLF for generating two and four scrolls. The number of scrolls achieved is the number of saturated levels, in this case with amplitudes , , , and . The breakpoints of the SNLF are , , and . Figure 2(g) represents the corresponding electronic realization. To generate two scrolls, the circuit requires two op-amps with two saturated levels at and and a shifted voltage (breakpoint) at . To make four scrolls, it requires four op-amps with four saturated levels at , , , and and three shifted voltages at , , and .

2.2.7. Lead/Lag Phase Network (LL)

The circuit of Figure 2(h) is a particular case of a lead/lag phase network with transfer functionwhere we can use (8) to translate the frequencies of the zero and the pole by a factor . It is important to remark that the circuit of Figure 2(h) is an op-amp-based lead compensator, and for the element values as they are labeled in the network, the transfer function becomes (12), which is a scaled version of the approximation in (1).

3. Fractional-Order Integrators

This section presents two realizations of fractional integrators using the blocks of Figure 2. We will carry out Monte Carlo and sensitivity analyses in them to study how the uncertainty in the output voltage is allocated to the tolerances of the circuit elements. In addition, we will tune some circuit elements to diminish the slope error in the magnitude frequency response caused by the first-order CFE approximation.

3.1. Circuit Implementation

In order to achieve a circuit implementation of the first-order CFE approximation of the fractional integrator given by (1), consider a transfer function in frequency-domain, , expressed in the same form, i.e.,where and are input and output voltages. An algebraic manipulation on (13) leads toand dividing both sides of (14) by and collecting similar terms, we get

Figures 3(a) and 3(b) illustrate the corresponding block diagram and electronic realization using circuits WDA, IInv, and WAA. A comparison of (15) with (5) and (6) results in

To determine and , we apply (7) and (8) after designating the degrees of freedom and . Later, once chosen the degree of freedom , for the order required and using (1), we calculate (circuit elements are denoted in blue in Figure 3(b)). The rest of resistances have a value of . As an example, after adopting the degrees of freedom , , and , the resistances shown in the table of Figure 3(b) result for each order .

Other alternative of implementation is with the circuit of Figure 2(h) which corresponds to the transfer function (12). To determine , we apply (8) after designating the degrees of freedom and . Later, once chosen the degrees of freedom and , for the order required and under , we calculate , , and (circuit elements are denoted in blue in Figure 3(c)).

3.2. Experimental Validation and Slope Adjustment

We simulated the fractional integrator of Figure 3(b) in HSPICE and breadboarded it with op-amps TL081 (gain bandwidth product = 100 kHz, dc gain = 200 V/V, and supply voltage = ±15 V). The input was a sinusoidal voltage with 0.1 V in amplitude and frequency ranged from 1 kHz to 100 kHz. An NI ELVIS II + laboratory platform provided this signal and bias potentials of ±15 V for the op-amps. Figure 4(a) shows the measured frequency responses for . These plots presented a satisfactory agreement with the results in HSPICE and with numerical simulations of (1) in Matlab. In consequence, the first-order CFE approximation of the fractional integrator causes the slope errors of the frequency responses reported in Table 1, with negligible influence of the circuit elements and no idealities of the op-amps. In agreement with Table 1, these errors are considerable (except for , , and ), extending from 0.09 dB to 4.73 dB. This performance is prohibitive to carry out chaotic oscillators, which are susceptible to their parameters and initial conditions. To reduce these errors, for each value of , we did an experimental tuning of (see the table of Figure 3(b)). This is equivalent to adjusting through , getting the frequency responses of Figure 4(b). From this figure and according to Table 1, the adjustment reduces the slope errors to values in the range from −0.09 dB to 0.6 dB, which produces values of closer to the theoretical ones. The corresponding curve fitting of Figure 5 brings us a design equation to calculate respect to , resulting in

For the circuit of Figure 3(c) and using the adjusted values of of Table 1, for , , and , the resistances of values , , and shown in the table of Figure 3(c) are obtained. Here, and . Although there are other alternatives for implementing fractional-order integrators, as discussed in Section 1, the first-order CFE implementation used in this work is very compact and does not require circuit elements with noncommercial values. Therefore, to the best of the authors’ knowledge, this is the first experimental validation of a multiscroll chaotic system based on active fractional-order circuits, and thus, there are not any other works where either simulated or experimental results are reported, and consequently, those cannot be compared with our proposal. In fact, other fractional-order multiscroll chaotic systems based on passive ladder networks as reported in [43, 44] require more hardware, including resistors and capacitors with noncommercial values.

3.3. Monte Carlo Analysis

To verify that the circuit of Figure 3(b) is robust to tolerances, we carried out a Monte Carlo analysis in SPICE. We dealt with tolerances of ±5% and ±20% in the resistors and capacitors using a normal distribution with N = 100. The variable to examine is the slope of the magnitude response for each value . To simplify, Figure 6 depicts only the histograms corresponding to the adjusted circuit with and . We observe distributions similar to a normal distribution. Table 2 outlines the obtained results for each value of and for both cases, with and without adjustment. The resulting standard deviations are small, in the range from 0.027 dB/dec to 0.23 dB/dec, and comparable in both cases. Figure 7(a) illustrates the mean slope obtained for each order. Without tuning, these values move away from the line that describes the theoretical slopes. Figure 7(b) depicts the largest and smallest values that the slope gets without tuning. The maximum error is in the range from 0.5 dB to 4.5 dB. No one order (except ) encloses the curve of the ideal slope. From Figure 7(c), we observe that the greatest and smallest values of the slope with tuning are closer to the ideal line, enclosing it. The maximum error is in the range of 0.6 dB. This way, the error in the first-order CFE approximation without adjustment influences the slope more than the tolerances.

3.4. Sensitivity Analysis

In assessing the sensitivity of the output voltage of the fractional-order integrator of Figure 8 regarding its resistive elements, by performing a normalized sensitivity simulation, the results reported in Table 3 are achieved. Note that only four of the thirteen resistors of the circuit are in Table 3 since these are those who exhibit a “major” impact on the integrator’s sensitivity. As appreciated, different scenarios may occur depending on the value of , but in general terms sensitivity is a little less when the adjustment has been applied. Without loss of generality, we can say that the larger the value of is, the larger the sensitivity becomes. Still, sensitivity is highly acceptable. Let us take the worst case when and a variation of the nominal value of produces a deviation of in the output voltage of the circuit. Even though this value is low, a lesser variation can be attained if a high-precision resistor (with a tolerance ) is used for . We can apply the same criteria for the rest of the resistors in terms of tolerances; however, there is a trade-off between sensitivity and realization cost, and high-precision resistors are also more expensive. It is important to emphasize that the sensitivity of the output voltage at a frequency of 0 Hz , i.e., at DC regime, and is equivalent to determine the sensitivity of the transfer function, , when , i.e., for . And thus, we are evaluating variation as a function of the resistive elements.

To evaluate the sensitivity of the fractional-order integrator regarding capacitor , Figure 9 depicts the behavior of the normalized sensitivity of for both with and without adjustment of parameter . We define the normalized sensitivity aswith nominal values and , and with given by (13). As illustrated, the tendency is rising towards the of variation as changes of its nominal value along with the increment.

According to the results reported above, we can conclude that the proposed fractional-order integrator exhibits a low sensitivity regarding variations of the values of its passive elements.

4. Fractional-Order Multiscroll Lü Chaotic System

In this section, we design electronic circuits based on the building blocks of Figure 2 and the fractional integrators of Figure 3 to realize commensurate fractional-order Lü chaotic system, i.e., system with equal fractional order for integrators. A commensurate realization is acceptable to test the functionality and benefits of using the active fractional-order integrators of Figure 3, but incommensurate fractional-order oscillators, as those reported in [4547] or [48] can also be realized with the proposed approach without increasing hardware requirements. In fact, the proposed approach is an attractive alternative for the design of incommensurate fractional-order oscillators because it is necessary to calculate the parameter using (17) instead of designing ladder networks.

The fractional-order Lü chaotic system is a realization that represents the transition from the Lorenz to the Chen attractor [49, 50]. The Lü system displays the following system of equations:where , , and are state variables, , , , and are positive real constants, is the fractional order, and is a SNLF similar to that represented in Figure 2(f). In [43, 49], the authors determine, for a Lü chaotic oscillator, the suitable incommensurate fractional order to preserve the chaotic behavior with , getting , , and and Lyapunov exponents , , and . Taking as reference that design, we will design our 2-scroll and 4-scroll commensurate realizations with the same parameters , , , and , and with and . The stability can be verified by using the method reported in [51].

Figure 10 depicts the block diagram and design details of the commensurate fractional-order system in (20). Figure 11 shows the corresponding experimental setup with fractional integrators as those depicted in Figure 3(b) and all the circuits realized using TL081 op-amps. The setup comprises a Keysight DSOX3054A oscilloscope, two DC power supplies RIGOL DP832 and BK precision 1550, and the PCB realized with the software Altium Designer 19. The design details of the PWL block for a 4-scroll chaotic attractor are presented in this figure for a saturation voltage . and are the bias voltages of the op-amps TL081, while and are the voltages required by the block . The 2-scroll chaotic attractor can be obtained with this circuit, as is illustrated in Figure 2(g). The four-input IWAA is designed with the degree of freedom and the desired values for , , , and . We incorporated a gain to the first and second approximations of the fractional integrator to add tuning capacity. The parameter was calculated using (17) and the degrees of freedom for the first approximation were , , , , and . For the second approximation, we used , , , and . Values of and were used for and to approximate the frequency response of the integrators to the response of equivalent passive realizations. We can observe that, with the first approximation, the chaotic oscillator requires for implementation 54 resistors, 3 capacitors, and 14 op-amps. With the second approximation, the circuit comprises 27 resistors, 6 capacitors, and 11 op-amps. Figures 12(a) and 12(b) show the SPICE simulation and experimental measurement of the 2-scroll chaotic attractor projected onto xy-plane (phase space portrait) for . Figures 12(c) and 12(d) show the SPICE simulation and experimental measurement of the same projection for . We can observe an agreement between simulation, theoretical, and experimental results. Figures 13(a) and 13(b) show the SPICE simulation and experimental measurement of the projection of the 4-scroll chaotic attractor onto -plane for . Figures 13(c) and 13(d) show the SPICE simulation and experimental measurement of the same projection for . Again we observe agreement between simulation and experimental results.

Table 4 shows a comparison of the hardware requirements of the two proposed realization and the fractional-order oscillators reported in [43, 44]. Despite the active realization of the fractional integrators, we observe that the two proposed implementations reduce the number of required op-amps (reducing power consumption), do not require noncommercial capacitors, and have a comparable number of noncommercial resistors. Unlike the resistors used in [43, 44], the values of the noncommercial resistors of the proposed realizations are easy to approximate with combinations of commercially available resistors, avoiding bulky realizations difficult to redesign. Another advantage is the easy design of the fractional-order integrators using (17) instead of calculating the elements of ladder networks.

5. Conclusion

This paper presented the design and electronic implementation of a commensurate fractional-order multiscroll chaotic oscillator with orders and . The fractional-order integrators realized with tuned first-order active realizations are robust to tolerances in the passive elements and favor the use of resistors and capacitors with commercial values, reducing hardware requirements and design and implementation efforts. The easy design of the fractional-order integrators with improved accuracy using (17) instead of calculating ladder networks could favor the practical realization of other types of systems, such as incommensurate chaotic oscillators, fractional-order filters, fractional-order controllers, and fractional-order memristors.

Finally, the physical realization of fractional-order chaotic systems with hidden dynamics (attractors without equilibrium points, multistability, extreme multistability, etc.) along with their synchronization schemes using the proposed approach will be conducted in future works to explore secure communications based on fractional chaos.

Data Availability

The necessary data are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This work was supported by projects CONACYT (258880), Plan de Trabajo CA (BUAP-CA-276) by SEP-PRODEP, and VIEP-BUAP.