Abstract
In this paper, a novel prototype laboratory is presented for engineering education, in which experiments are based on the fractional calculus. The prototypes of analog and digital fractional-order proportional-integral-derivative (PID) controllers are built in the laboratory. These fractional-order PID controllers are applied to linear and nonlinear plants to demonstrate the effectiveness of fractional-order calculus in real time. These experiments are designed, developed, and implemented on the analog and digital platforms. These controllers are integrated to control the DC motor, brushless DC motor, and magnetic levitation modules through hardware-in-loop as well as stand-alone systems. The analog type of fractional-order PID implementation is carried out by using passive components (i.e. resistances and capacitances) with an operational amplifier. However, real-time digital implementation is carried out using field-programmable gate array and digital signal processor. This paper describes how the experiments on fractional calculus can be tailored for graduate, undergraduate students’ education and extended for research in this emerging area.
Keywords
Introduction
The concept of fractional calculus and fractional-order circuits is emerging slowly in the engineering domain. 1 Fractional calculus acts as a generalization to the classical calculus. 2 It shows many beneficial applications in several fields. It includes control systems, industrial drives, magnetic levitation (MagLev), inverted pendulum, electric motors, robotics, and biomedical applications as presented in the literature.3–14 These concepts are somewhat new to circuit theory and have exhibited applications in fractional-order hold circuit, 4 power electronics, 15 integrator and differentiator circuits,16,17 and filters18,19 with potentially many other engineering applications. Chen et al. 20 offered a modeling and analysis method for the fractional-order DC–DC converters. Dimeas et al. 21 implemented fractional-order controller in analog form for a DC motor. Debnath 22 has precisely reviewed the applications of fractional calculus in diverse domain of engineering and science. Fractional-order models of supercapacitors, batteries, and fuel cells have been surveyed in Freeborn et al. 23
Fractional differential equations characterize the fractional systems. 6 In order to solve the fractional differential equations, a combination of Laplace transform and fractional calculus methods can be used. 3 It is a known fact that the fractional systems involve infinite memory; however, the integer-order systems are local and point property. 5 It is significant to express the fractional-order system using a finite difference equation. For that, integer-order approximations are mainly used in continuous-time domain.3–6 Some methods are available that can be used to compute the integer-order approximation of fractional operator.3,24–28 The approach to realize analog fractional systems is available in the literature.25–27 Adhikary et al. 29 reviewed some significant research work, along with their success and limitations, in the field of fractional circuit realization.
The design of digital fractional system includes the discretization of fractional-order operators.
28
Direct discretization and indirect discretization are the two generally practiced methods. Generally, the ideal frequency response is approximated using either an infinite-length impulse response (IIR) filter or a finite-length impulse response (FIR) filter.
3
FIR filters are often selected due to their effective implementation aspects. As compared to IIR filters, FIR filters have no constraint on maximum sampling rate and they are always stable.
30
Experimentation plays an important role in learning and research.
31
Fractional systems can be implemented using a microcontroller, digital signal processor (DSP), field-programmable gate array (FPGA), etc. Jiang et al.
32
proposed that the FPGA is an effective platform to realize high throughput, high quality approximations of fractional systems. In Muresan et al.,
33
an FPGA-based fractional controller is developed for a class of DC motors. The major objectives of this paper include:
A common algorithm of constant phase shaping is proposed to select the pole-zero (p-z) pairs for analog and digital implementation. Analog fractional integrator circuits are implemented using commercially available passive components. The analog fractional-order PID (FO-PID) controller is tested with DC motor emulator as a plant. Digital fractional integrator is implemented on FPGA. The digital FO-PID controllers are applied to DC motor, brushless DC (BLDC) motor, and MagLev system.
To address the all above objectives, the “Fractional Calculus Engineering Laboratory” is developed at Visvesvaraya National Institute of Technology, Nagpur, India. This work is organized as follows: the preliminaries of a fractional calculus are given in the following section. The laboratory experiments are discussed in “Examples of laboratory experiment” section. The paper is concluded in “Conclusion” section with the summary of lab outcomes.
Preliminaries
Fractional calculus investigates the feasibility of using the real/complex number powers of the integration operator. The basics of fractional calculus are available in the literature.3–6 The real-order operator is inferred as follows in equation (1)
Method of approximation
Every transfer function can be represented by its pole-zero pairs.3–6 The Bode’s phase plot of fractional transfer function is constant at

Flowchart of proposed method.
The basic steps of the algorithm are as follows:
Step I: Consider the values of α and the desired frequency band (ωL, ωH). Step II: Calculate the pole-zero pairs using equation (4). Step III: Plot asymptotic and exact Bode plot. Check the average value of phase angle. If the average value of the angle over the desired frequency band is in required tolerance band then go for either analog or digital implementation. Step IV: For analog implementation, find out the values of resistances and capacitances. Step V: For digital implementation, convert the continuous time transfer function into discrete time transfer function using “Tustin” as discretization method and choosing sample time.
Examples of laboratory experiment
Various experiments are designed based on the fractional calculus and presented next.
Lab 1. Design of analog fractional integrator
The aim of the experiment is to realize the fractional integrator using analog components. The algorithm, as shown in Figure 1, gives p–z pairs for α = −0.2, −0.5, and −0.8 as listed in Table 1. The corresponding phase plot, with frequency band ωL = 100 rad/s and ωH = 10000 rad/s, is presented in Figure 2. The rationally approximated transfer function is realized by the cascading of basic two-port networks, where one-port network is having first-order impedance function with a zero. Partial fraction expansion of Znew gives

Combined logarithmic phase plot for integrator orders 0.2, 0.5, and 0.8.

Proposed circuit implementation of one p–z pair.
p–z pairs obtained using pole-zero approximation method.

Fractional-order integrator circuit with six p–z pair.
R–C values used for implementation of s−0.5.
The realized analog circuit for an integrator is shown in Figure 5. Experimental setup to test the fractional integrator circuits is shown in Figure 6. It consists of function generator, digital storage oscilloscope, and ±15 V power supply. The input signals are applied from the function generator to confirm the integrator behavior.

Photograph of analog fractional integrator.

Experimental setup of analog fractional integrator. DSO: digital storage oscilloscope.
The circuit is tested for the square wave input. Figure 7 shows time response captured for fractional integrator circuits. The gains are set to get readable plots.

Time response (square) of analog fractional integrator.
Learning: This experiment explains about the design of analog fractional integrator to the students. However, for analog implementation of each fractional integrator, different circuits are required. Component’s derating has occurred in the case of analog implementation which affects the results. Hence, digital implementation is suggested and practiced.
Lab 2. Design of digital fractional integrator
The discretization is a key concept in realization of a digital fractional integrator. For indirect discretization, two steps are needed, i.e. frequency-domain fitting in continuous time and then discretization of the fit continuous time transfer function. 28 The p–z pairs, obtained by proposed algorithm as shown in Table 1, are discretized using first-order hold, impulse invariant, zero-order hold, matched pole-zero pair, and “Tustin” approximation. 34 Bode’s plot for s−0.5 digital fractional integrator is presented in Figure 8. It interprets that digital integrator with “Tustin” approximation matches with an integrator in continuous time and also it retains the stability.

Bode’s plot for s−0.5 digital fractional integrator. FOH: first-order hold; ZOH: zero-order hold.
The constant phase p–z interlacing algorithm for digital fractional integrator of α = −0.2, −0.5, −0.8 within the desired band of frequency ωL = 100 rad/s and ωH = 10000 rad/s gives discrete p–z pairs as in Table 3 with gain values of 0.11666, 0.0089841, and 0.00036712 respectively.
Discrete p–z pairs obtained using approximation.
For the digital implementation of a fractional integrator, finite impulse response (FIR) filter structure is proposed. FIR filter gives good results while implementing on digital platform. 30 FIR approximations can be derived by choosing the order of filter “N.” The more terms an approximation uses, it will represent the original operator more accurately. Here, the order of filter “N” is chosen as 64, 128, and 1024.
Design example:
Consider the value of α; say −0.5. Approximate it using the proposed method. The continuous time transfer function will come as
Discretize equation (5) using “Tustin” approximation method and choosing sampling time Ts = 0.0001 s
Take the values of discrete p–z pairs and convert it into FIR filter form choosing order of filter N. By considering N = 64, the transfer function will come as shown in equation (7)
where x[n] is the the input signal, y[n] is the output signal, and bi is the value of the impulse response at the ith instant for
The digital integrator is verified using FPGA-in-the-loop (FIL) wizard of MATLAB/Simulink. Experimental setup to test the digital circuits is shown in Figure 9. It consists of a DE2-115 FPGA development board and Intel Core i7 Processor. In Figure 10, the block diagram of FIL implementation is shown. Here, the digital integrator is realized using Altera’s DE2-115 FPGA development and education board. The digital integrator is successfully tested and the obtained results are presented. The digitally implemented fractional integrator with N = 64 gives the response as shown in Figure 11. N = 64 is chosen to reduce the computation burden, while the circuit is tested up to 1024. The similar response, as observed with analog circuit, is seen with a digital implementation which shows the potential of proposed algorithm and its implementation.

Experimental setup of digital fractional integrator. FIL: FPGA-in-the-loop; FPGA: field-programmable gate array.

FIL block diagram. D/A: digital-to-analog; FIL: FPGA-in-the-loop.

Time response (square) of digital fractional integrator.
Learning: This experiment explains about the design of digital fractional integrator to the students. From the results obtained, it can be stated that the value of α indicates the strength of integration operation. For higher value of α, the output approaches to the output of integrator of order one, whereas, for small α, the output is nearer to the input signal.
Advantages of digital implementation of fractional integrator: An important advantage offered by digital realization is in the flexibility of its modifying controller characteristics. The ability to “redesign” the fractional operator/controller by changing software (rather than hardware) is an important feature of digital realization. Implementation of fractional operator for wide range of frequency is possible which was earlier constrained by the limitation of analog realization. However, the digital realization with digital processors, the range can be extended to the desired limit with virtually limitless computing power.
Limitations of digital implementation of fractional integrator: The main issues associated with the digital realization are related to the effects of sampling and quantization. The digital realization involves a cascade of analog-to-digital (A/D) conversion system followed by a digital processor which is, in turn, followed by digital-to-analog (D/A) conversion system. The A/D conversion process involves the sampling of analog signal at discrete instants of time separated by sampling time Ts. The D/A conversion process involves reconstructing continuous time signals from the samples given by the digital processor. The ill effects of sampling can be reduced, if not completely eliminated, by sampling at a very high rate.
With the availability of low-cost, high-performance digital processors and interfacing hardware, the implementation problems in digital realization do not pose a serious threat to its usefulness. The advantages of digital realization overweigh its implementation drawbacks for most of the applications.
Lab 3. Analog FO-PID controller for DC motor emulator
To verify the analog FO-PID controller, emulator of DC motor (Plant) is developed. More details are available in Aware et al. 26 Plant consists of DC motor and load with specification as speed = 2000 r/min, Ra = 2 W, L = 3 mH, J = 1.78e−4 kg m2, K = 1.02, V = 24 V. The experimental setup is shown in Figure 12. The response of DC motor model with analog FO-PID controller is shown in Figure 13. The FO-PID performance indicates that the FO-PID with α = ±0.5 could make it possible to maintain the desired control on the output.

Experimental setup of analog FO-PID with emulator. DSO: digital storage oscilloscope; FO-PID: fractional-order PID.

Result of DC motor emulator with FO-PID.
Learning: This experiment explains about the design of analog FO-PID controller to the students. In this experiment, DC motor emulator is used in place of actual DC motor to put more emphasis on controller implementation part. The results can be compared with integer-order PID controller also. The isodamping property 5 of FO-PID controller can be observed by varying DC forward gain.
Lab 4. Digital FO-PID controller for hardware-in-the-loop (HIL) MagLev system
MagLev system is basically an electromagnetic system which levitates ferromagnetic objects in space by magnetic force induced due to the electric current flowing through coils around a solenoid. The continuous time approximation method is discretized and a novel digital FO-PID controller is developed. More details about the scheme are discussed in Chopade et al. 7 This digital FO-PID controller is used to control highly nonlinear MagLev system in HIL mode as given in Figure 14. The controller is implemented in MATLAB/Simulink.

Hardware-in-loop MATLAB/Simulink interface of MagLev. MagLev: magnetic levitation.

Position control of MagLev with PID controller 7 .
The performance analysis for digital PID and FO-PID controllers has been carried out. The desired and measured ball position control with PID controller is shown in Figure 15. The result in Figure 16 shows that the better control over position accuracy with lesser control efforts “u(t)” can be obtained with FO-PID controller.

Position control of MagLev with FO-PID controller.
Learning: This experiment explains about the design of digital FO-PID controller using MATLAB/Simulink. Control effort “u(t)” minimization is verified with this experiment. Control efforts are more in case of PID controller, hence, enveloped response is seen. The PID controller takes previous values in integration and works with instantaneous values at a point for differentiation. Whereas, the fractional-order integrations and differentiation take previous values into consideration with decreasing weights and therefore it has a memory effect in it. This memory effect gives ideal filtering action and thus, the smooth control signal is achieved using the FO-PID controller as it is filtered. Hence, FO-PID controller shows a clean signal. Fractional calculus-based system does the control action with a lesser effort than the conventional classical calculus-based controllers, therefore they are more efficient. Digital FO-PID controller is found to improve the position accuracy and control efforts of MagLev.
Lab 5. FIL digital FO-PID controller for BLDC motor
The speed control scheme for FO-PID controller-fed BLDC motor drive is shown in Figure 17.

Block diagram of BLDC motor speed control scheme. PMBLDC: Permanent Magnet Brushless DC Motor.
A digital FO-PID controller-fed BLDC motor drive with four-switch inverter is modeled in this experiment. The FO-PID controller is used as a speed controller. As the speed reference is given, the FO-PID controller provides the control signal in accordance with the generated speed error. The generated reference current is given as input to the current controller. Further, the current controller provides switching pulses, which are obtained by comparing the reference input signal and the different modes provided by the motor module to operate the switches of an inverter. Here, the scheme is having FIL. The digital FO-PID controller is implemented on FPGA DE2-115 board. More details are available in Khubalkar et al. 9 The FIL block is shown in Figure 18. The speed response of PID and FO-PID controlled scheme is displayed in Figure 19. The scheme with FO-PID controller is showing improved results than its integer counterpart.

FIL diagram of BLDC motor speed control scheme. A/D: analog-to-digital; D/A: digital-to-analog; FIL: FPGA-in-the-loop; FO-PID: fractional-order PID; PMBLDC: Permanent Magnet Brushless DC Motor.

Speed response BLDC motor at transient state. FO-PID: fractional-order PID; PID: Proportional-Integral-Derivative.
Learning: This experiment explains about the design of digital FO-PID controller using FPGA to the students. In this experiment, conventional six-switch inverter topology can be used. As the design is digital, various forms of fractional-order controllers can also be verified.
Lab 6. Digital FO-PID controller for DC motor
Speed control scheme of buck converter-fed DC motor drive is shown in Figure 20. In Figure 21, experimental setup of digital FO-PID-fed DC motor is displayed.

Block diagram of DC motor speed control scheme. FO-PID: fractional-order PID; PWM: Pulse width modulation.

Experimental setup of DC motor drive with FO-PID.
Here, a reference speed of motor is selected as set point for the controller. The controller gives the control signal calculating the error generated in the control scheme. The PWM pulses are generated at particular duty cycle corresponding to the controller’s output. This signal operates the switch of DC–DC converter and controls the armature voltage of the DC motor. The controllers are implemented on TMS320f28377s DSP platform and used to control the speed of DC motor. More details about the scheme are available in Khubalkar et al.8,35 Speed control of DC motor with reference speed of 1000 r/min is shown in Figure 22. The response improves considerably with FO-PID as shown in the results.

Speed response DC motor at transient state. FO-PID: fractional-order PID; PID: Proportional-Integral-Derivative.
Learning: This experiment explains about the design of stand-alone digital FO-PID controller using DSP to the students. Fixed point and floating point DSPs can be used to implement the controller. By changing the coding of DSP, different forms of FO-PID controllers can be tested on the setup. The FO-PID controllers are robust against plant parameter variation. This property can be studied while performing this experiment on two different DC motors of same rating but having different plant parameters. 35
Utilization in research
All these experiments are effective in learning and research in the domain of fractional calculus. The various properties viz. iso-damping property, robustness against plant parameter variation, control effort minimization, etc. of fractional-order controllers are verified using the laboratory experiments. Undergraduates, postgraduates, and research scholars are the beneficiaries of this unique laboratory.
Conclusion
The research contribution of the laboratory is to provide efficient and practical way to implement fractional-order systems. This laboratory provides setup to overcome the barriers between simulations of fractional systems and real-time solutions. The pole-zero interlacing algorithm is realized and effectively used for developing the analog and digital fractional integrators. In the similar way, fractional differentiators are also developed and included in the laboratory. The analog circuit realization with commercial R–C components is achieved. FPGA and DSPs are very effective platforms for implementing fractional systems. The designed analog and digital FO-PID controllers are tested on various test benches. Analog controller is tested with DC motor emulator. Digital FO-PID is tested with MagLev as HIL. Then, the controller is tested with BLDC motor as FIL mode. Finally, stand-alone digital FO-PID is tested with 1.5 kW DC motor drive. All these experiments give a real-time exposure to fractional calculus and designed in the way that learner can understand the concepts effectively.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by Board of Research in Nuclear Sciences, project sanction no. 2012/36/69.
