Improved control of integrating cascade processes with time delays using fractional-order internal model controller with the Smith predictor

This article proposes an enhanced control of integrating cascade processes. The importance of the scheme shows by controlling a double integrating process with time delay under significant parametric uncertainties and load disturbances. A novel design is proposed based on the Smith predictor principle and uses a fractional-order internal model controller in the outer loop. The required fractional-order tuning parameter follows the desired gain and phase margins. At the same time, another tuning parameter, such as the fractional filter’s time constant, is calculated from the desired performance constraint. Numerical analysis and comparison are performed to showcase the feature of the hybrid structure. The simulation results show that the fractional-order internal model controller–based scheme provides enhanced tracking and faster regulation capabilities and works well under nominal and parameter uncertainty conditions.


Introduction
This article makes original contributions to two areas of controller design.They are second-order integrating plus time delay (IPTD) plants and fractional-order theory into internal model control (IMC) design.Therefore, the background concerning each of the two domains is first reviewed.The series-cascade control structure (SCCS) is a prime focus of this work.
The second-order or double integrating plant with time delay is ramping and relatively slows due to its considerable time constant.This plant primarily possesses two poles at origin, has a non-self-regulating nature, and is disturbed easily with a minor load disturbance.It is due to poor performance rejecting load disturbance by a conventional feedback system and deviating the controlled variable from a setpoint.Therefore, the controlling method requires a structural change or additional controllers to handle such plants.In this area, a cascade control structure has been adopted widely in plant industries, such as level, temperature, pressure, and flow control loops.Franks and Worley 1 had introduced an SCCS.The closed-loop performance was improved using the inner loop by rejecting disturbances quickly.Such cascade control structure consists of two control loops.The inner-loop controller is called the secondary or slave controller, whereas the outer-loop controller is called the primary or master controller.In a typical cascade structure, faster dynamics of the inner loop help in speedier disturbance attenuation, thus minimizing the effect of disturbance before they affect the primary output.Another objective of cascade control is to minimize the sensitivity of the primary plant gain variations.However, the control strategy in cascade is more complex than a single-loop control.The structure also measures an additional output and controller to be tuned.
In industries, two different cascade controllers can be seen.The control and disturbance signals affect the primary output response through the secondary or inner loop in series-cascade control.In parallel-cascade control, the control and disturbance signals simultaneously affect primary and secondary outputs. 2 Researchers in the last decade discussed the design and analysis of cascade control techniques for improved tuning and performances [3][4][5][6][7][8] and references within.It was found in the literature also a plant with input time delay or dead time could add more complexity to cascade control.The input delay increases the phase lag, leading to a decrease in gain margin and phase margin, finally leading to instability.The work with this kind of plant is a difficult task.
Our focus in this article is the cascade control method for integrating time delay plants.Therefore, the research on series-cascade control design with integrating plants is reviewed in the following.In Kaya and Atherton, 9 a cascade integrating time delay plant was controlled using a proportional-integral (PI)-proportional-derivative (PD) scheme with the Smith predictor (SP) in the outer loop and IMC in the inner loop.An SP-based cascade control structure was seen in Uma et al., 10 in which a secondary loop was again designed using IMC.Nevertheless, a primary loop was tuned using the direct synthesis technique for better setpoint and load disturbance rejection.A modified cascade control system was proposed in Padhan and Majhi 11 using SP and IMC rules and involved two controllers and a setpoint filter.Another scheme presented in Nandong and Zang 12 for multi-scale cascade control consisting of four controllers was applicable to self-regulating, double integrating, and unstable plants.An online tuning method in cascade control was developed by Jeng 13 for higher-order plants.An SP strategy was continued for cascade control to enhance the performance.In C xakırog˘lu et al., 14 higher-order models were decomposed into first-order models, and depending on the order of the model, a separate loop was established.A modified cascade controller is based on the method of moments and the Routh stability criterion developed in Raja and Ali. 15 Recently, a frequency domain method with model matching was discussed for cascade controllers 16 in a nonlinear continuous stirred tank reactor.
8][19] The methods with FOC provide an extra degree of freedom in tuning parameters.The fractional-order proportional-integral-derivative (FOPID) presented by Podlubny 20 has two additional parameters: fractional differential and fractional integral orders.Since then, many tuning approaches based on FOPID structures have been proposed in the literature.It has also seen some variations with fractional-order derivatives in the SP and IMC schemes.The fractional-order PI plus D controller was developed for integrating plants using the explicit formulae of the complex root boundary. 21This method and some references within have shown better tracking performance than the classical approaches.It has been found recently a type of controller like FOI 1Àl D l using a fractional integrator and derivative actions.Oustaloup recursive approximation (ORA) method is used to implement this controller. 22ractional calculus was also applied in the IMC theory.Arya and Chakrabarty 18 discussed fractionalorder internal model controller (FOIMC) scheme using the gain and phase margins' specifications.Similarly, Ranjan and Mehta 23 presented a superiority using a fractional-order tilt double derivative (FOTDD) in the IMC framework.Also, a new double-loop control approach was proposed for a time-delayed unstable system that uses a PD/P controller in its inner loop and FOIMC in its outer loop. 24To continue, Shweta et al. 25 developed a novel dual-loop hybrid control method for second-order unbounded plants with dead time and zeros.In this method, an external-loop controller was designed using the FOIMC, whereas the internal-loop controller used a simple proportional-integral-derivative (PID).Kumar et al. 26 suggested FOIMC for achieving satisfactory servo performance and PD to reject the disturbance rejection.It is observed the performance was improved with dual-loop and fractional-order than before. 27A compensator-based series-cascade control with FOIMC was proposed for unstable processes in Mukherjee et al. 28 In Chandran et al. 7 and Pashaei and Bagheri, 29 the secondary controller was designed from FOIMC.In general, it is noticed that the fractional-order filter in IMC can provide a more flexibility to modify the output performances.However, there are some limitations to the available methods and those are presented below.
The literature review suggested that the available methods often consume much time in computing parameters when the cascade control is used.Also, a few methods are available for integrating plants; specifically, double integrating plants are significantly less.If the integrating plant is affected by load disturbance, it will further deteriorate the relationship between input and output.Moreover, many real-time plant models contain time delays.Even considering a large time delay in plants, the perturbation may lead to poor performance.A considerable number of works have been proposed on this subject.However, they may require more tuning parameters or controllers in a loop or complex design approach.Despite some modifications suggested in the structure, the results showed poor disturbance rejection and robustness.These investigations have assisted in understanding the necessary and developing a new control scheme for cascade controlling of integrating plants.Briefly, the contributions of this article comprise as follows: The improved structure is suggested with only one fractional-order parameter, and a hybrid scheme contains the merits of both IMC and SP principles.This method can handle highly sensitive processes under parametric uncertainties and load disturbances.The fractional-order parameter in the controller could be selected based on gain margin (g m ) and phase margin (f m ) specifications.The desirable range is provided for ease of selection tasks.The fractionalorder IMC filter's time constant is tuned by the balanced, objective function under constraints.The robustness study is performed to demonstrate the advantages of the proposed scheme, which includes parameter perturbation, external and internal load disturbances, and noisy output signal.The quantitative analysis also shows superior performances compared with recently developed works.
The article is organized as follows.Section ''Preliminary'' represents the general IMC scheme and tilt-integral-derivative (TID) controller.In sections ''Plant models, controller, and tuning procedures'' and ''The proposed controller design in real-time,'' the proposed controlling scheme is discussed in detail.Then, the selection of tuning parameters is suggested in section ''Tuning of l and b for the proposed controller structure.''Section ''Robustness analysis and stability'' presents the robustness and stability analysis of the proposed controller.The numerical simulation is performed in section ''Simulation results,'' compared with other methods, and finally, the conclusion is given in section ''Conclusion.''

Preliminary
This section presents the basic IMC structure and equations related to it and also some details about the seriescascade control system.

The general IMC control system
The well-known IMC scheme is shown in Figure 1, where G p (s), G m (s), Q(s), and G c (s) representing the plant, plant model, IMC controller, and feedback controller, respectively.The design procedure using the IMC principle can be seen briefly below.
Step 1: factorize the plant model, G m into two parts.G À m as invertible and minimum phase.G + m as noninvertible and non-minimum phase parts (delays and right half-plane zeros) Step 2: take inversion for invertible part G À m .
Step 3: add a low-pass filter f(s) in order to make IMC controller proper.Then, the IMC controller becomes Step 4: the simplified form of the feedback controller is obtained as The fractional-order tilt-integral-derivative controller The well-known FOPID involves the tuning of five parameters: commonly gains K p , K i , K d and fractional orders a and m.It has been demonstrated in the literature that it can offer additional flexibility to meet design control goals.The TID controller also belongs to the same fractional-order family.Despite having only four parameters, the structure is almost identical to the widely used PID.Only the proportional gain is replaced by the tilted gain (K t ) having a transfer function 1=s 1 n , where n is a non-zero real number. 30The standard transfer function of TID controller can be represented as where K t is the tilted gain, K i is the integral gain, and K d is the derivative gain.The tuning parameter, n should be between the numbers 2 and 3. A TID controller exhibits more accurate tuning, quick disturbance rejection, and minimizes parameter uncertainty when compared to conventional controllers.The study of this topic has recently attracted a lot of interest.A robust 3-degree-of-freedom (3DOF)-TID controller 31 as well as integral derivative-tilted (ID-T) controller 32 developed for load frequency control in interconnected power systems was some of the works developed recently.The transfer function of the controller can be written below where a and m are the positive real orders for integral and derivative, respectively.Here, we want to apply the IMC principle designing the fractional-order tilt-integral-derivative (FOTID) of form in equation (5).
Aiming to make a simple tuning approach, even though the controller seems to be complex with three fractional orders.The proposed design is discussed in the following sections for integrating plant models.We attempt to obtain one common tuning parameter to calculate other parameters and use the plant information.

Series-cascade control system
Figure 2 shows the typical structure of a series-cascade control system, which consists of two control loops, namely, the inner loop with plant G p2 and the outer loop with plant G p1 .Let us assume the load disturbances entering inner and outer loops with values D 1 , D 2 , and D 3 as seen in Figure 2. The primary plantcontrolled variable Y 1 (s) with setpoint R 1 (s) is used by the primary controller G c1 to build a setpoint for the secondary controller G c2 .The general guideline behind this configuration is that the inner loop should be faster than the outer loop, that is, the disturbance affecting the inner loop is compensated before it affects the primary plant.The outer loop can be designed for setpoint tracking, but the problem arises when the outer plant will have a large time delay.In such cases, the SP can be used for compensating the delay.

Plant models, controller, and tuning procedures
The series plant is mainly combined with the stable inner plant, whereas the primary plant may be stable, unstable, or integrating by nature.In our work, we have considered the following plant models to study.Let the IPTD model as The double integrating plus time delay (DIPTD) model is The integrating second-order plus time delay (ISOPTD) model is The proposed SCCS scheme is given in the block diagram (Figure 3).As seen in the figure, it has two main controllers and a filter compensator.First, the inner-loop controller is designed using the IMC principle.Following the inner loop, the outer primary controller is developed as an FOC.Aiming to design a new fractional controller for integrating plants to make the overall system more stable and robust with load disturbance inputs.Let us express the primary and secondary controllers as The G c1 in equation ( 9) is called a fractional-order integrator and derivative (FOID) controller and a, m are the fractional orders of the integral and differential.Also, one can use the gains K c1 =Ti1 = K i1 and K c1 T d1 = K d1 for simplifying design steps.The G c2 in equation ( 10) is the classical PI controller.

The inner-loop controller design
The model of a secondary process G m2 (s) having a small time delay can be controlled effectively using the wellknown simple internal model control (SIMC) tuning. 33et us consider the stable first-order plus time delay (FOPTD) model as Then, following the IMC technique from section ''The general IMC control system,'' one can obtain the IMC controller for the inner loop as where l 2 is a single tuning parameter and defined also as a secondary closed-loop time constant.From equation (3), the secondary controller can be obtained as below after approximating e Àu 2 s to (1 À u 2 s) using the Taylor series expansion The objective of the inner controller is to reject the load disturbance before affecting the outer plant response.According to the SIMC tuning rule, a fast response with good robustness is obtained when Comparing with the classical PI in equation (10), we obtain the parameter values as below Now, the inner-loop transfer function is given by Substitution G p2 and G c2 in equation ( 16) gives After replacing the small value time delay e Àu 2 s = 1 À u 2 s, the approximated inner-loop becomes Since it is straightforward to simplify the inner-loop, the primary controller can be designed from the following section.

The proposed outer-loop controller design
To design the main loop controller, the overall plant transfer function can be considered after the modified inner loop as below where u m represents the overall time delay, and Ĝp1 represents the total plant transfer function.
The overall primary closed-loop transfer function for servo response can be written as By considering the model exactly represents the process ( Ĝp1 (s) = G m e Àu m s ), equation ( 20) can be represented as Case 1: IPTD plant.Let us take G p1 as equation ( 6), where K 1 is the static gain and u 1 is a dead time.Using equation ( 19), the overall transfer function for the outer loop becomes where u m = 2u 2 + u 1 .Aiming to design a fractionalorder IMC filter, we have chosen the form of IMC filter with fractional-order differentiator as below where l represents a fractional filter time constant (l .0) and b is a positive real number, that is, the fractional order b, lies in the range 0 ł b ł 2. The value of b varies with different plants, and the choice of b as per plant type is further discussed in section ''Choice of b.''Now, according to the IMC rule, one can obtain the outer-loop controller G c1 using equation (3).Following the basic rule of the IMC, the FOC becomes The above equation can be seen as a fractional-order D controller as below where K d1 = 1=K 1 l.The overall closed-loop transfer function after substituting in equation ( 21) can be obtained as Case 2: DIPTD plant.The DIPTD plant in equation ( 7) can obtain the outer controller as below after the inner-loop transformation Finally, the outer controller is obtained as below Case 3: ISOPTD.The ISOPTD plant in equation ( 8) is formed as the outer controller as below The above equation can be rearranged as below This can be denoted as an FOID controller and its parameters become As per the transfer function models of the outer process, the controller settings are tabulated in Table 1.

The proposed controller design in realtime
Let us take first Ĝp1 (s) from equation (22), where K 1 is the static gain and u m is a dead time.This IPTD model can be approximated into the FOPTD model as below where g is considered as a constant with a high value, say g = 100.As per the new filter equation ( 23), the resulting IMC controller becomes IPTD: integrating plus time delay; DIPTD: double integrating plus time delay; ISOPTD: integrating second-order plus time delay.
It can be realized as the feedback controller from equation (3) below The above equation can be rearranged with respect to TD controller in equation (5).In this case, one can see the expression as Comparing equations ( 37) and ( 5), the TD controller values are obtained as below The DIPTD plant in equation (7) in the same way can be approximated as In the same way, one can obtain the feedback controller as shown below The above equation can be arranged with respect to TD controller as Comparing equations (41) and ( 5), the TD controller values are obtained as below It can be seen for IPTD and DIPTD process models that the integral part is not required.Thus, the parameter K i = 0.
The ISOPTD plant in equation ( 8) can be represented as The feedback controller can be obtained as The above equation can be arranged with respect to TID controller as Comparing equations (45) and ( 5), the FOTID parameters become The proposed controller settings as per practical cases can be noted as shown in Table 2.

The disturbance filter for an inner loop
The literature suggested some restrictions with IMC and SIMC methods.They perform poorly with disturbance inputs and have less flexibility to improve performance for robustness. 3,15It is essential to improve load rejections, especially the plants are in cascade.In our work, a disturbance filter is required in the inner loop.The tuning of a filter should be simple and based on model parameters.For that purpose, we have included a disturbance filter from the inverse model and converted it into a proper second-order filter.A general form 1=(cs + 1) 2 can be included as below into the filter, F d (s) as

Plant
Primary controller IPTD: integrating plus time delay; DIPTD: double integrating plus time delay; ISOPTD: integrating second-order plus time delay.
Since the inner plant model is a first-order stable, the above equation becomes The time delay can be approximated with a firstorder Taylor series expansion, and then F d is expressed by A smaller value of c results in a faster recovery, which is recommended for the cascade plants.Considering the ease of design step, c is chosen as u 2 =2 in this work.

The robustness filter for an outer loop
The model uncertainties are critical for any controller scheme.As seen in Figure 3, the filter, namely, F f , is included to act on errors between the actual and predicted outputs.It is also proven that a filter helps to deal with a time delay issue. 34Let us take F f (s) = 1=t f s + 1.It is noted that a small value of t f constant yields a fast regulatory response and small settling time.Whereas a high value can result in more settling time and sluggish response.For our plant type, it is aimed to focus on robustness.After a numerical study on all integrating plant types, t f is directly related to plant time delay, and it can be taken simply as 2u m .
Tuning of l and b for the proposed controller structure A plant industry is well practiced with a classical controller design from the maximum sensitivity and margins specifications. 35These specifications can handle various dynamics including small and large time delay plants.The basic concept is to convert an open-loop transfer function by choosing controller's zeros equal to poles of the model.In the proposed work, we have considered a fractional-order pole with order b and the loop equation becomes Comparing equations ( 50) and (51), one can easily calculate, l = 1=K.It is to remind that a parameter l appears from the fractional-order filter equation (23).Now, the basic definitions for the gain and phase margins can be written as below where the gain margin and phase margin as g m and f m and their crossover frequencies as v gc and v pc , respectively.Substituting equation (51) into equations ( 52)-( 55), following set of expressions is obtained Using equations ( 57) and (58), one can write After multiplying u m both sides of equation ( 60) and then further modification with equations ( 56) and (59), the new phase and gain margin relation can be obtained as below Differently from a classical relation, the above expression provides an extra parameter to tune for desired g m and f m .This is well related to the new fractional-order ID controller.Same way, equations ( 56) and (57) can be modified for l as below Now the possible values of b are observed using relation in equation (61).As recommended values for g m are between 2 and 5 and for f m are between 308 and 608.The interesting result can be seen in Figure 4.

Choice of b
According to numerical values, the typical choice of b shall be between 0:8 and 1:5.So, the fractional-order IMC filter's parameters are bounded by recommended ranges for the gain and phase margins.As per this requirement and model type, we have found for IPTD plant selection from 0:8 \ b \ 1:0, for DIPTD plant 0:8 \ b \ 1:5 and for ISOPTD plant 1:0 \ b \ 1:5.This enables the estimation task to be faster and overcomes the problem of selection gain and phase margins.

Final tuning based on a performance index
The performance criterion is defined by the following cost function equation (63) to fulfill dual requirements for given b where the first term calculates the control signal variations and the second term is the integral of the timesquared error.U j and Y 1 ( j) are the outer plant input and output at time t = t j , with M data points.It can be seen from equation (63) that the index can tradeoffs between control signal variations and transient setpoint errors by choosing any suitable value of n 2 (0, 1).In order to guarantee robustness in the cascade control system, the parameters are calculated from the minimum of the cost function.

Summary of tuning procedure
The two non-linear expressions from equations ( 61) and ( 62) are important relation in this study.The robustness and stability of any closed-loop control system are related with g m and f m values.The proposed method's required parameters are initially chosen from the first choice of g m and f m .The major steps for the complete control scheme are summarized in Table 3.

Robustness analysis and stability
It is necessary to analyze the robustness and stability analysis of any closed-loop system in the presence of model uncertainties.The model uncertainties are widespread in the plant industries.The robust stability of the closed-loop system can be analyzed with the small gain theorem.We have considered uncertainties in plant gain and dead time.Now, the condition to satisfy for the closed-loop system to be stable is 36 jjD m ( jw)T( jw where T( jw) is the complementary sensitivity function, and D m ( jw) is the bound on the plant multiplicative uncertainty.The uncertainty in the plant can be expressed as where G m ( jw) is the nominal model that is used for controller design, and G p1 ( jw) is the actual plant.If uncertainty is present in both the plant gain and dead time, then the tuning parameters have to be selected such that For any integrating model discussed here, the complementary sensitivity function for the proposed scheme is given by Table 3. Tuning algorithm.
Step 2: choose f m and g m and compute b to satisfy equations ( 61) and (62).
Step 3: calculate J min from equation (63) for suggested range of l.
Step 4: choose l for the lowest cost function value.
Step 5: calculate G c1 parameters as per the plant type.
Also, to ensure robust closed-loop performance, the constraints to be followed by the sensitivity and complementary sensitivity functions are where w m ( jw) is the uncertainty bound on the sensitivity function, which is given as S( jw) = 1 À T( jw).Therefore, tuning parameters must be selected so that the resulting controller satisfies the robust performance and stability limitations.

Simulation results
Recent works presented by Siddiqui et al., 16,37 Raja and Ali, 15 C xakırog˘lu et al., 14 Jeng, 13 Padhan and Majhi, 11 and Uma et al. 10 were considered for comparative analysis.The response parameters such as rise time (t r ), settling time (t s , s), and integral of time-weighted squared error (ITSE) are calculated to demonstrate performances.A t r is defined as the time the output takes from 10% to 90% of its steady-state value.The settling time is measured from the output to reach and steady within a given tolerance band.We have considered the 5% band throughout all examples.For evaluating the manipulated input usage, we compute the total variation (TV) in the control signal, which is the sum of all its moves up and down.It is a good measure of the signal's smoothness and should be as small as possible.If required, a setpoint filter should be added to remove the undesirable overshoot.A low-pass filter of the form 1=(hs + 1) is simple, and parameter h is selected appropriately by the designer. 38The responses and numerical analysis with and without a setpoint filter are shown in all examples.Also, the white noise having a variance of 0.001 is added as a measurement noise to verify the control scheme's robustness against the noisy output.An extended investigation of the proposed scheme in realtime can be seen from the third example following details from section ''The proposed controller design in real-time.''Example 1.First example is extensively studied in the literature.The outer loop G p1 (s) = 2e À2s =s and inner loop, G p2 (s) = 4e Às =s + 1 are considered as a series plant.Recent method by Siddiqui et al. 16 presented two types of controllers for the same plant.In a method 1, the primary G c1 (s) = 0:06 + 0:0027=s + 0:02s and secondary G c2 (s) = 0:255 + 0:212=s were suggested.Then, in a method 2, the controllers were taken as G c1 (s) = 0:06 + 0:002s and G c2 (s) = 0:29 + 0:17=s.Uma et al. 10 also have studied the same process models with control settings: the primary, G c1 (s) = (1:8 + 0:5=s + 1:93s)½0:5s + 1=0:6561s 2 + 2:960s + 1, secondary, G c2 (s) = s + 1=2s + 4, and disturbance rejection controller, G cd (s) = (0:0436 + :1206s)½0:75s 2 + 1:5s + 1=1:8314s 2 + 0:7686s + 1.
The cost function values with b are shown in Figure 5.The corresponding responses are plotted in Figure 6.The results show that one set of parameters is chosen among possible values of b to provide the most optimal responses with balanced control performance.Considering this observation, b has selected as 0:95.The corresponding g m and f m are 2:3 and 55:178, respectively.The final controller settings are K c2 = 0:125, K i2 = 0:125, and K d1 = 0:806.The disturbance rejection filter obtained as F d = s 2 + 2s + 1=s 2 + 4s + 4 and that of robustness filter is F f (s) = 1=8s + 1.The setpoint filter of the form f s = 1=(3s + 1) is also added to improve the setpoint smoothness.A positive load disturbance of 0:05 is added at 100, 300, and 500 s, respectively.The closed-loop nominal responses and control responses along with other approaches were shown in Figures 7  and 8, respectively.It can be seen the proposed method results a lower overshoot and improved disturbance rejection performances.Table 4 provided the same agreement with the responses by minimum ITSE and TV values.Furthermore, 10% uncertainty is applied in plant parameters gain and delay, in both primary and secondary models.Figure 9 again depicted the robustness from parameter uncertainties.Figure 10 shows the magnitude of the complementary sensitivity function for different perturbations.It is evident that the robust  stability margin lowers as perturbation increases, and the controller cannot meet the robust stability requirement stated in equation (66), even with a 20% uncertainty in both gain and time delay.When we tested the methods under the measurement noise with variance = 0.001, Siddiqui et al. 16 and Uma et al.'s 10 method obtained large variations in a control signal.The noise test also proved the proposed scheme's robustness as seen in Figures 11 and 12.
Example 2. Controlling a double integrating cascading plant is more challenging.A series plant, G p1 (s) = e À0:7s =s 2 and G p2 (s) = e À0:3s =(s + 1) was studied by C xakırog˘lu et al. 14 and Siddiqui et al. 16 In C xakırog˘lu et al., 14 the output decomposition method is presented with two outer-loop controllers, G c1a (s) = 0:833 and G c1b (s) = 0:333, and inner controller as G c2 (s) = 3:33s + 1=s.Whereas in Siddiqui et al., 16 G c2 = 3:74 + 7:95=s and G c1 = 0:16 + 0:02=s + 0:56s were suggested.After adopting the proposed controller, l = 1:28 and b = 1:02 are obtained from the chosen g m = 2:5 and f m = 52:288.As seen in Figure 13, this choice provided a minimum cost value.Then, the controller parameters are K c2 = 1:67, K i2 = 1:67, and K d1 = 0:781.The setpoint filter f s = 1=(1:2s + 1) is added to improve the smoothness of the response.The disturbance rejection filter F d = 0:3s 2 + 1:3s + 1=0:0225s 2 + 0:3s + 1 is used in the secondary loop for enhancing disturbance rejection and that of robustness filter is F f (s) = 1=2:6s + 1.The step disturbances of magnitude 20.5 are applied at 100 and 200 s, while 20.1 is applied at 300 s.From Figure 14, it is observed that the proposed method shows a similar setpoint response with C xakırog˘lu et al.However, the disturbance rejection capability and its control efforts are superior to both methods, seen in Figure 15.Table 5 shows the quantitative analysis    Furthermore, the perturbations test with 10% changes in gain and time delay for both primary and secondary plants is shown in Figure 16.One can see     From Figure 17, it can be inferred that the maximum permitted perturbation is 70% because any additional rise in perturbation causes output from the system to become excessively oscillatory.When we tested the methods under the measurement noise with variance = 0.001, Siddiqui et al.'s 16 method obtained large variations in a control signal.The noise test also proved the proposed scheme's robustness as per Figures 18 and 19.Example 3. A large time delay cascade integrating plant, with G p1 (s) = e À6:5672s =s(3:4945s + 1) and G p2 (s) = 2e À2s =s + 1 studied by Raja and Ali. 15 They have applied a series-cascade control with inner and outer loops as classical PIs.The parameter settings were K c1 = 0:039, T i1 = 845:34, K c2 = 0:2717, and T i2 = 1:656.Also, a stabilizing proportional controller, K p = 0:0012 was added in the structure.For the same plant, Padhan and Majhi 11 designed two controllers and a setpoint filter.They have given, a setpoint filter 6:989s 3 + 5:4945s 2 + s + 1=628:806s    9:2060s + 1=45:3194s 3 + 29:0943s 2 + 7:8726s + 1).As it could be seen that their method suggested a higherorder function.
For the same plant, our method suggested a fractional controller with f m = 53:168 and g m = 5.Following the design steps, we have obtained parameters l = 3:38 and b = 1:2.Figure 20 presents the index value with chosen l and b values.Then, the inner-loop and outer-loop parameters become K c2 = 0:125, K i2 = 0:125, K i1 = 0:145, and K d1 = 0:506.From Table 2, the TID controller values can be obtained based on minimum cost index value l = 1, as K t1 = 0:01, K i1 = 1:03, and K d1 = 3:494.The disturbance rejection filter F d = 2s 2 + 3s + 1=2s 2 + 4s + 2 and that of robustness filter is F f (s) = 1=21:1344s + 1.The disturbance inputs are considered with a magnitude of 0:2 at t = 200 s, at t = 400 s and 0:1 at t = 600 s.The servo and regulatory responses and control efforts are shown in Figures 21  and 22, respectively.The setpoint filter of the form f s = 1=(10s + 1) is added to improve smoothness of response.To illustrate the robustness test, perturbation of 10% in the time delay and plant gain of both       6, it is observed that the output obtains the desired setpoint value quickly with minimum TV and optimal ITSE values.In general, the closed-loop responses from the proposed method are superior to the improved load disturbance rejections.The outputs and control variations with measurement noise in Figures 24 and 25 show that Padhan and Majhi's method shows more control variations compared to others.In addition to performance analysis, robustness stability analysis for the proposed scheme has also been done and is shown in Figure 26.
Example 4. A higher-order integrating plant with G p1 (s) = e À3s =s(10s + 1)(2s + 1) and G p2 (s) = e À0:5s =s + 1 studied by Siddiqui et al. 37 and Jeng. 13Jeng suggested a master controller with parameters as K c1 = 0:0862, T i1 = 56:34, and T d1 = 8:50 and slave controller as K c2 = 1:336 and T i2 = 0:833.Siddiqui et al. suggested secondary controller settings as G c2 (s) = 2:68 + 2:81=s + 0:45s and primary controller settings as G c2 (s) = 0:1 + 0:001=s + 0:65s. 37ow, the proposed method designs the fractional IMC filter using f m = 59:078 and g m = 3.Following the procedure, we have obtained l = 3:1 and b = 1:01 as per values from Figure 27.The controller settings are K p2 = 1:0, K i2 = 1:0, K i1 = 0:323, and K d1 = 3:226.The setpoint filter f s = 1=(8s + 1), the robustness filter F f (s) = 1=12s + 1, and the disturbance filter F d = 0:5s 2 + 1:5s + 1=0:0625s 2 + 0:5s + 1 are included as per the structure.The closed-loop responses under the presence of disturbances 20.5 at 300 and 500 s and 20.05 at 700 s are shown in Figure 28.The    corresponding control signals are plotted in Figure 29.For encapsulating the system robustness, a perturbation of 10% is introduced in gain and dead time of G p2 and G p1 .The test results are given in Figure 30 with other methods' performance.The numerical values are measured and listed in Table 7. From these results, it can be     proven that the new fractional IMC with SP structure provides better tracking and robustness in comparison to the recent literature method.When we test the methods with measurement noise as shown in Figures 31 and  32, the Siddiqui et al.'s 37 method again shows worst performance and large variations in the control signal.It can be noticed from Figure 33 that maximum allowable perturbation is allowable up to 50%.

Conclusion
This article developed a hybrid structure with IMC and SP concepts for integrating plants with time delay.The suggested method includes the fractional-order IMC filter, having an additional degree for improving performance.This design is simple, with two tuning parameters, which follow good robustness under significant disturbance inputs and parameter perturbation.Moreover, a selection task becomes manageable using the balanced cost function and margins (g m , f m ) relationship.The system's performance can be implemented in all classes of integrating plants, including thirdorder plants in the outer loop.Quantitative analysis is carried out using ITSE and TV, and the proposed method shows promising results compared to previous works with the same or more tuning parameters.As    one can see, there is no or minimal overshoot even without a setpoint filter, which is the key advantage compared to other approaches.The designer can choose the setpoint value so that desired response concerning settling time and rise time can be accomplished.Using the proposed controllers, simulation results show robust stability, better control performance, and disturbance rejection.Now, research is underway to extend the proposed method to fractional-order integrating plant models and other types of plants with unstable and positive/ negative zeroes.Research can be done to develop a fractional controller for plants with delay, even though the SP makes the system complex.Future research can also be conducted for series-cascade unstable models and the related extension of the proposed methodology.The authors believe that the method will provide a considerable effort toward fractional-order control in industrial use.

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.

Figure 3 .
Figure 3. Proposed series-cascade control structure for integrating plants.

Ĝp1
(s)G c1 (s) = Ke Àu m s =s bð50Þwhere K represents loop transfer function gain and can be determined based on gain and phase margin specifications.Now, as per plant models and their controller structures in section ''The proposed outer-loop controller design,'' the open-loop transfer function can be written as Ĝp1 (s)G c1 (s) = e Àu m s ls b ð51Þ

Figure 4 .
Figure 4. Region of b as per g m and f m specifications.

Figure 5 .
Figure 5. Calculated cost function values with selected b.Figure 6. Responses for selected b.

Figure 6 .
Figure 5. Calculated cost function values with selected b.Figure 6. Responses for selected b.

Figure 10 .
Figure 10.Magnitude plot of the complementary sensitivity function for Example 1.

Figure 13 .
Figure 13.Variation of l and J min to obtain b for Example 2.

Figure 17 .
Figure 17.Magnitude plot of the complementary sensitivity function for Example 2.

Figure 18 .
Figure 18.Servo and regulatory responses with measurement noise for Example 2: (a) proposed, (b) proposed with f s , (c) C xakırog˘lu et al., 14 and (d) Siddiqui et al.-2.16

Figure 20 .
Figure 20.Variation of l and J min to obtain b for Example 3.

Figure 22 .
Figure 22.Control signal for Example 3: (a) proposed, (b) proposed with f s , (c) proposed (real-time), (d) Raja and Ali, 15 and (e) Padhan and Majhi. 11 primary and secondary plants have been considered at once and shown in Figure23.From the numerical analysis shown in Table

Figure 26 .
Figure 26.Magnitude plot of the complementary sensitivity function for Example 3.

Figure 27 .
Figure 27.Variation of l and J min to obtain b for Example 4.

Figure 28 .
Figure 28.Servo as well as regulatory responses for Example 4: (a) proposed, (b) proposed with f s , (c) Siddiqui et al., 37 and (d) Jeng et al.13

Figure 31 .
Figure 31.Servo and regulatory responses with measurement noise for Example 4: (a) proposed, (b) proposed with f s , (c) Siddiqui et al., 37 and (d) Jeng et al.13

Figure 33 .
Figure 33.Magnitude plot of the complementary sensitivity function for Example 4.

Table 1 .
The proposed controller settings as per process model.

Table 4 .
Performance comparisons from Example 1.
ITSE: integral of time-weighted squared error; TV: total variation.

Table 5 .
Performance comparison for Example 2.

Table 6 .
Performance comparison for Example 3.

Table 7 .
Performance comparison for Example 4.