**1. Introduction**

In order to reduce the waste of fossil energy and CO2 emissions, many countries have released different policies. China has released a policy document to fulfill its target of reaching peak carbon emissions by 2030. The United States of America released policy to cut carbon emissions in half by 2030. A lot of renewable power technologies were also proposed. However, most applied energy still comes from the combustion of fossil fuels. In addition, the control performance has a close relationship with the utilization of fossil fuels [1]. In this paper, the fossil fuels in ships are focused. Many of the ships are equipped with internal combustion engine, or gas turbine. However, for large scale ships, the power systems based on boiler-turbine still occupy a large proportion [2–4]. For example, many aircraft carriers are powered with boiler-turbine system. Hence, a lot of academics and companies are doing research to improve the control performce of the boiler-turbine system [5–9], of which there are three manipulated variables and three controlled variables with complicated dynamics. For the marine steam power plant, the disturbance and the energy required changes more frequently compared than that on land. However, there is not so much research about the control for the marine boiler-turbine system.

In the boiler-turbine system, the interactions of rvariables and constraints are the mainly reasons which make it difficult to obtain a satisfied control performance. The input variables for the system are the flow rates of fuel, steam to the turbine, and feedwater to the boiler, while the output variables are the steam pressure in the drum of the boiler, power required of the turbine, and the water level of the drum. The constraints are the limitations for the actuator, including the upper and lower bound, and rate limiter. In addition, power

**Citation:** Zha, S.; Wang, S.; Cajo, R.; Ren, W.; Li, B. Power Tracking Control of Marine Boiler-Turbine System Based on Fractional Order Model Predictive Control Algorithm. *J. Mar. Sci. Eng.* **2022**, *10*, 1307. https://doi.org/10.3390/jmse 10091307

Academic Editors: Bowen Xing, Marco Altosole, Maria Acanfora and Flavio Balsamo

Received: 7 August 2022 Accepted: 12 September 2022 Published: 15 September 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**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/).

requirement changes a lot, which is usually treated as disturbance. In order to compensate unknown disturbances in the boiler-turbine system, a high order sliding mode observer was designed for a baseline exponentially stable feedback controller [9]. To improve the economy of the boiler-turbine system, the economic index was utilized directly in the cost function, and a global economic optimum routine was obtained for the system [10]. In the literature [11], the model of the boiler-turbine system was linearized and decoupled with an adaptive feedback linearization method, and a second order sliding mode controller was designed to deal with the disturbances and uncertainties. Ref. [12] applied an online policy iteration integral reinforcement learning method to the boiler-turbine system, and optimal tracking control performance was obtained.

The model predictive control( MPC) has the advantages in dealing with the nonlinear dynamics, interactions and constraints problems [13], hence, MPC is a preferred choice for researchers. The MPC is studied in many fields such as building energy management [14], landscape office lighting regulatory system [15,16], wind turbines [17], tank-system [18], pressure oscillation adsorption process [19,20], permanent-magnet synchronous motor [21], autonomous underwater vehicle [22], and so on. For the boiler-turbine system, there are also some applications of MPC. A nonlinear model predictive control method was designed for boiler-turbine system with a data driven model [23], and the optimal problem was solved by immune genetic algorithm. Ref. [24] presented a zone economic model predictive cotroller to fulfill the economic target of the boiler-turbine system. Fuzzy model predictive control was designed to realize load tracking and economy of the boiler-turbine system, and a fuzzy model was used to approximate the nonlinear dynamics of the system [25]. Other MPC applications on boiler-turbine system can be found in [26–28].

For the MPC method, the weighting factors have a significant different effect on the control performance, and the number of the weighting factors is large, which makes it diffcult to obtain a good choice for these parameters. For example, if the *Nip* denotes the prediction horizon for the *i*th output, and *Njc* for the *j*th control horizon, the number of weighting factor will be ∑*<sup>n</sup> <sup>i</sup>*=<sup>1</sup> *Nip* + <sup>∑</sup>*<sup>m</sup> <sup>j</sup>*=<sup>1</sup> *Njc*, where *m* and *n* are the numbers of output and input. The commonly used method to tune the weighting factors is trial and error or choose them empirically. Another alternative way to choose the control parameters is optimization method [29,30]. However, due to the high dimension of the weighting factors, it is difficult to optimize them. So, in this paper, the fractional order model predictive control(FOMPC) for the boiler-turbine system is proposed. By two fractional order papermeters (one for the tracking error, and another for the control effort), the weight factors matrices can be obtained, which reduces the difficulties in weighting factors configuration. The boilerturbine system is a nonlinear multiple inputs and multiple outputs system, so the nonlinear distributed structure of MPC is studied for the system.

The rest of the paper is structured as follows. The boiler-turbine system is formulated in Section 2; Section 3 presents the nonlinear distributed MPC with EPSAC framework. Fractional order MPC is designed for the boiler-turbine system in the Section 4. The simulation experiments are presented in Section 5. The last section gives the conclusions.

#### **2. Boiler-Turbine System**

The boiler-turbine system is a core part in the power plant, and Figure 1 shows the structure of the system [9]. The details elements are indicated in the figure, and they are listed as follows: 1—drum; 2—superheater; 3—water spray desuperheater; 4—valve for the steam to turbine; 5—turbine high-pressure cylinders; 6—foward control valve; 7—backward control valve; 8—turbine middle- and low- pressure cylinders; 9—shafting; 10—condenser; 11—replenish water; 12—condensate pump; 13—low-pressure heater; 14—deaerator; 15—feed water pump; 16—high-pressure heater; 17—feed water valve; 18—economizer; 19—downcomers; 20—water-cooled walls; 21—furnace; 22—heat flow control valve; 23—nozzle; 24—blower; 25—preheater for air; 26—air conditioner; 27—flue gas baffle; 28—induced draft fan; 29—flue; 30—gearbox; 31—turbine.

**Figure 1.** The structure of a boiler-turbine unit.

Due to the lack of data for marine boiler-turbine sytem, it is difficult to obtain its model. The structure is similar to that on land. Hence, a model of on load power plant is chosen for study, and the nonlinear model of the boiler-turbine system is shown as follows according to the literature [31]:

$$
\dot{x}(t) = F(x(t)) + G(x(t))u(t), \tag{1}
$$

and the *F*(*x*(*t*)), *G*(*x*(*t*)) are defined as follows:

$$F(\mathbf{x}(t)) = \begin{bmatrix} 0 \\ -0.1\mathbf{x}\_2 - 0.016\mathbf{x}\_1^{9/8} \\ 0.0022\mathbf{x}\_1 \end{bmatrix} \tag{2}$$

$$G(x(t)) = \begin{bmatrix} 0.9 & -0.0018x\_1^{9/8} & -0.15\\ 0 & 0.073x\_1^{9/8} & 0\\ 0 & -0.0129x\_1 & 1.6588 \end{bmatrix} \tag{3}$$

where the inputs *u* = [*u*1, *u*2, *u*3] *<sup>T</sup>* for the system are the valve opening of fuel, steam to the turbine and feedwater to the drum. The states are drum steam pressure, power required for the turbine and steam water density denoted by *x* = [*x*1, *x*2, *x*3] *<sup>T</sup>*. The outputs are drum steam pressure, power required for the turbine and water level in the drum. The level of the drum can be calculated as:

$$L = 0.05(0.13073x\_3 + 100a\_s + q\_c/9 - 67.975) \tag{4}$$

and *qe* = (0.854*u*<sup>2</sup> <sup>−</sup> 0.147)*x*<sup>1</sup> <sup>+</sup> 45.59*u*<sup>1</sup> <sup>−</sup> 2.51*u*<sup>3</sup> <sup>−</sup> 2.096; *<sup>α</sup><sup>s</sup>* <sup>=</sup> (1−0.001538*x*3)(0.8*x*1−25.6) *<sup>x</sup>*3(1.0394−0.0012304*x*1) . The *qe* and *α<sup>s</sup>* denote the evaporation rate and steam quality, respectively.

The rates and amplitudes limitaion for the inputs are listed as follows:

$$\begin{cases} -0.007 \le \frac{du\_1}{dt} \le 0.007 & 0 \le u\_1 \le 1\\ -2.0 \le \frac{du\_2}{dt} \le 0.02 & 0 \le u\_2 \le 1\\ -0.05 \le \frac{du\_3}{dt} \le 0.05 & 0 \le u\_3 \le 1 \end{cases} \tag{5}$$

For the boiler-turbine, there are different operating points. To evaluate the effectiveness of the proposed method, experiments around the following operating points shown in Table 1 are carried out.


**Table 1.** Operating points for the boiler-turbine system.

The drum water level should always be kept as zero meters.

#### **3. Nonlinear Distributed MPC for the Boiler-Turbine System**

According to the introduction of the boiler-turbine, it can be found that this system is a nonlinear multiple inputs multiple outputs system. Hence, the nonlinear MPC is designed for the system with distributed structure.

#### *3.1. The Basic of the EPSAC*

This part presents the basic of EPSAC. For more details about the EPSAC, it can be found in the refs. [32–34]. For a discrete system, the system output can be expressed as:

$$y(t) = x(t) + w(t)\tag{6}$$

where *y*(*t*) is the system output; *x*(*t*) is the model output and the *w*(*t*) is the disturbances. *x*(*t*) can be calculated according to the model of the system as follow:

$$\mathbf{x}(t) = f[\mathbf{x}(t-1), \mathbf{x}(t-2), \dots, \mathbf{u}(t-1), \mathbf{u}(t-2), \dots] \tag{7}$$

In Equation (7), the *f*(*x*, *u*) denotes the model of the system, *x*(*t* − *i*) and *u*(*t* − *i*) *i* = 1, 2, . . . indicate the past model outputs and inputs.

In the EPSAC, the input scenario for the future is composed with two parts:

$$
\mu(t+k|t) = \mu\_{base}(t+k|t) + \delta\mu(t+k|t) \tag{8}
$$

where the *ubase*(*t* + *k*|*t*) and *δu*(*t* + *k*|*t*) are the basic and optimized future control actions. Then the future system output can be predicted as:

$$y(t+k|t) = y\_{base}(t+k|t) + y\_{opt}(t+k|t) \tag{9}$$

where *ybase*(*t* + *k*|*t*) is the result of the basic future control action; *ubase*(*t* + *k*|*t*) and *yopt*(*t* + *k*|*t*) can be calculated according to the optimized future control action *δu*(*t* + *k*|*t*).

The *yopt*(*t* + *k*|*t*) can be obtained with:

$$y\_{opt}(t+k|t) = h\_k \delta u(t|t) + h\_{k-1} \delta u(t+1|t) + \dots + g\_{k-N\_\varepsilon+1} \delta u(t+N\_\varepsilon - 1|t) \tag{10}$$

In Equation (10), the *hi* and *gi* are the impulse response and step response coefficients of the system, respectively; *Nc* and *Np* are the control horizon and the prediction horizon, respectively. The system output can be re-written in matrix form:

$$\mathbf{Y} = \mathbf{\hat{Y}} + \mathbf{G}\mathbf{U} \tag{11}$$

where **<sup>Y</sup>** = [*y*(*<sup>t</sup>* <sup>+</sup> *<sup>N</sup>*1|*t*)... *<sup>y</sup>*(*<sup>t</sup>* <sup>+</sup> *Np*|*t*)]*T*, **<sup>U</sup>** = [*δu*(*t*|*t*)... *<sup>δ</sup>u*(*<sup>t</sup>* <sup>+</sup> *Nc* <sup>−</sup> <sup>1</sup>|*t*)]*T*, **<sup>Y</sup>**¯ = [*ybase*(*<sup>t</sup>* <sup>+</sup> *<sup>N</sup>*1|*t*)... *ybase*(*<sup>t</sup>* + *NP*|*t*)]*T*; *<sup>N</sup>*<sup>1</sup> is the time delay of the system, and

$$\mathbf{G} = \begin{bmatrix} h\_{N\_1} & h\_{N\_1-1} & \dots & \mathcal{S}\_{N\_1-N\_c+1} \\ h\_{N\_1+1} & h\_{N\_1} & \dots & \dots \\ \dots & \dots & \dots & \dots \\ h\_{N\_P} & h\_{N\_P-1} & \dots & \mathcal{S}\_{N\_P-N\_c+1} \end{bmatrix}$$

The disturbance term *w*(*t*) in (6) includes all the effects on the system output. It can be modeled by a colored noise process as:

$$w(t+k|t) = \frac{\mathcal{C}(q^{-1})}{D(q^{-1})} w\_f(t+k|t) \tag{12}$$

where *q*−<sup>1</sup> is the backward shift operator.

In this work, the *C*(*q*−1) and *D*(*q*−1) are designed as follows:

$$\frac{C(q^{-1})}{D(q^{-1})} = \frac{1}{(1 - q^{-1})(1 - a e^{+j\alpha} q^{-1})(1 - a e^{-j\alpha} q^{-1})}\tag{13}$$

with *α* = 2*π f*0*Ts* and *a* ≈ 1; *Ts* is the sampling time and *a* ≤ 1 for stability. The cost function for the boiler-turbine system can be defined as:

$$J\_{\rm MPC} = \sum\_{k=N\_1}^{N\_2} p\_k [r(t+k|t) - y(t+k|t)]^2 + \sum\_{k=1}^{N\_u} q\_k \triangle u(t+k)^2 \tag{14}$$

The *pk* and *qk* are nonnegative weighting factors, and they are usually kept as constants. The matrix form of Equation (14) can be written as:

$$J\_{MPC} = (\mathbf{R} - \mathbf{Y})^T \mathbf{P} (\mathbf{R} - \mathbf{Y}) + \mathbf{U}^T \mathbf{Q} \mathbf{U} = (\mathbf{R} - \mathbf{\hat{Y}} - \mathbf{G} \mathbf{U})^T \mathbf{P} (\mathbf{R} - \mathbf{\hat{Y}} - \mathbf{G} \mathbf{U}) + \mathbf{U}^T \mathbf{Q} \mathbf{U} \tag{15}$$

with *<sup>P</sup>* = *diag*(*p*1, *<sup>p</sup>*2,..., *<sup>p</sup>*(*N*2−*N*1+1)) and *<sup>Q</sup>* = *diag*(*q*1, *<sup>q</sup>*2,..., *qNu* ).

For systems with constraint, the optimization problem can be solved with quadratic programming. Otherwise, the results of the optimal input part, which are indicated by *δu*(*t* + *k*|*t*), can be obtained as:

$$\mathbf{U}\_{\rm MPC}^{\*} = (\mathbf{G}^{T}\mathbf{P}\mathbf{G} + \mathbf{Q})^{-1}\mathbf{G}^{T}\mathbf{P}(\mathbf{R} - \mathbf{Y}) \tag{16}$$

#### *3.2. The Fractional Order MPC*

For the fractional order MPC, the cost function is designed as:

$$J\_{\rm FOMPC} = \gamma^{\gamma} \, I\_{\rm N\_1}^{\rm N\_2} p\_k [r(t+k|t) - y(t+k|t)]^2 + {}^{\lambda} \, I\_1^{\rm N\_c} q\_k \triangle u(t+k)^2 \tag{17}$$

where *<sup>γ</sup> I N*<sup>2</sup> *<sup>N</sup>*<sup>1</sup> and *<sup>λ</sup> <sup>I</sup> Nc* <sup>1</sup> indicate fractional order integral with fraction order of *γ* and *λ*; [*N*1, *N*2] and [1, *Nc*] are the integration intervals.

According to [35], the Equation (17) can be written by:

$$\begin{split} \mathbf{J}\_{\text{F}\text{OMPC}} &= (\mathbf{R} - \mathbf{Y})^{T} \mathbf{P} \mathbf{T} (T\_{5\prime} \gamma) (\mathbf{R} - \mathbf{Y}) + \mathbf{U}^{T} \mathbf{Q} \mathbf{A} (T\_{5\prime} \lambda) \mathbf{U} \\ &= (\mathbf{R} - \mathbf{\tilde{Y}} - \mathbf{G} \mathbf{U})^{T} \mathbf{P} \mathbf{T} (T\_{5\prime} \gamma) (\mathbf{R} - \mathbf{\tilde{Y}} - \mathbf{G} \mathbf{U}) + \mathbf{U}^{T} \mathbf{Q} \mathbf{A} (T\_{5\prime} \lambda) \mathbf{U} \end{split} \tag{18}$$

$$\Gamma(T\_s, \gamma) = T\_s^\gamma \operatorname{diag}(m\_{N\_2 - N\_1 \prime} m\_{N\_2 - N\_1 - 1 \prime}, \dots, m\_{1 \prime} m\_0) \tag{19}$$

$$\mathbf{A}(T\_{\mathbb{S}\prime}\lambda) = T\_s^{\lambda} \operatorname{diag}(m\_{\mathbb{N}\_{\mathbb{C}\prime}} m\_{\mathbb{N}\_{\mathbb{C}}-1\prime}, \dots, m\_1, m\_0) \tag{20}$$

The *mi* in Equations (19) and (20) with fractional order *α* can be calculated as:

$$m\_{\dot{j}} = \omega\_{\dot{j}}^{(-\mathfrak{a})} - \omega\_{\dot{j}-\mathfrak{a}}^{(-\mathfrak{a})} \tag{21}$$

where *n* is the number of the *mi*, and *ω* can be calculated with:

$$
\omega\_j^{-a} = \begin{cases}
(1 - (1 - a)/j)\omega\_{j-1}^{(-a)} & j > 0; \\
1 & j = 0; \\
0 & j < 0.
\end{cases}
\tag{22}
$$

According to Equation (19) to Equation (22), the weight matrix in the cost function *JFOMPC* can be easily tuned with fractional order *γ* and *λ*. For system with constraints, the optimization for input sequence can be solved with quadratic programming. For the case without constraints, the results for the optimal input *δu*(*t* + *k*|*t*) can be calculated as:

$$\mathbf{U}\_{\mathbf{FOMPC}}^{\*} = (\mathbf{G}^{T}\mathbf{P}(\Gamma + \Gamma^{T})\mathbf{G} + \mathbf{Q}(\Lambda + \Lambda^{T}))^{-1}\mathbf{G}^{T}\mathbf{P}(\Gamma + \Gamma^{T})(\mathbf{R} - \tilde{\mathbf{Y}}) \tag{23}$$

*3.3. Application of the Fractional Order EPSAC to the Nonlinear MIMO System with Distributed Structure*

The fractional order MPC introduced above is for the linear case. In order to apply the FOMPC to the boiler-turbine system, the nonlinear FOMPC is studied. According to the Equation (8) and (9), the principle of superposition is applied for linear system. To get over the superposition, the optimal future input *δu*(*t* + *k*|*t*) should be removed iteratively smaller tends to zero [36]. The procedure for the nonlinear MPC is summarized as follows:


The flow chart of nonlinear MPC is shown in Figure 2.

**Figure 2.** The flow chart of nonlinear MPC.

The boiler-turbine is a MIMO system, and there are strong interactions between variables. In order to calculate the optimal input sequence, the effect from coupling variables should be considered and the communication network should be established. In this work, the distributed structure is applied. The pseudocode is provided in Algorithm 1.

The boiler-turbine system is a nonlinear MIMO system, hence, the nonlinear MPC with distributed scheme will be applied. According to the procedure of nonlinear MPC and algorithm for the distributed MPC, the optimal future input sequence should fulfill the follow conditions:

$$\begin{cases} ||\delta \mathbf{U}\_i^{iter+1} - \delta \mathbf{U}\_i^{iter}|| \leqslant \varepsilon\_i \\ \delta \mathbf{U}\_i^{iter+1} \approx 0. \end{cases} \tag{24}$$

#### **Algorithm 1** Algorithm for the distributed MPC


#### **4. Simulation of the Fractional Order MPC on Boiler-Turbine System**

This section shows the simulation results of the fractional order MPC. Firstly, different fractional order terms are applied to different loops. Then, the best fractional order terms are applied to the drum steam pressure loop, power for turbine loop and water level loop, and the results are compared with the integer order MPC. Finally, the results are discussed.

The parameters configuration are listed in Table 2.

**Table 2.** Parameters for the MPC.


In Table 2, *Nci* and *Npi* (*i* = 1, 2, 3) are control horizon and prediction horizon, respectively; *Ns* is the number of simulation steps. The termination conditions for the nonlinear iteration are set as: *δ***U***iter <sup>i</sup>* - 0.05; or the iteration times *iternmpc* > 5. The termination conditions for the distributed MPC are set as: ||*δ***U***iter*+<sup>1</sup> *<sup>i</sup>* − *<sup>δ</sup>***U***iter <sup>i</sup>* || -0.005; or *iter* > 5.

#### *4.1. The Influence of Fractional Order Terms to the Different Loops*

In order to test the effect of different fractional order on the control performance, different fractional order terms are introduced to the cost function for each loop. The details for fractional order terms are listed in Table 3. In this work, the *γ* for the reference tracking and *λ* for the control effort are chosen the same for simplification.

**Table 3.** Fractional order terms for each loop.


According to the results shown in Figures 3–5, it can be seen that the effectiveness of different fractional order terms varies a lot. For the drum steam pressure control, the best fractional order term is 1.5; for the required power control, it is 1; and for the drum water level control, the fractional order term of 1.7 is the best.

**Figure 3.** The drum steam pressure with different fractional orders.

**Figure 4.** The required power for turbine with different fractional orders.

**Figure 5.** The drum water level with different fractional orders.

#### *4.2. The Influence of Fractional Order Terms to the Different Loops*

This section shows the comparison experiment between FOMPC and traditional MPC. The fractional order terms applied in the boiler-turbine system are chosen as 1.5, 1, 1.7 for the three loops, respectively.

In order to evaluate the effectiveness of the proposed method, the following performance indexes are compared, including Integrated Absolute Relative Error (*IARE*), Integral Secondary control output (*ISU*), Ratio of Integrated Absolute Relative Error (*RIARE*), Ratio of Integral Secondary control output (*RISU*) and combined index (*J*).

$$IARE\_i = \sum\_{k=0}^{N\_b - 1} |r\_i(k) - y\_i(k)| / r\_i(k) \quad (i = 1, 2, 3) \tag{25}$$

$$ISIL\_i = \sum\_{k=0}^{N\_i - 1} \left( u\_i(k) - u\_{ssi}(k) \right)^2 \quad (i = 1, 2, 3) \tag{26}$$

$$RIARE\_i(C\_2, C\_1) = \frac{IARE\_i(C\_2)}{IARE\_i(C\_1)} \quad (i = 1, 2, 3) \tag{27}$$

$$RISLI\_i(\mathbb{C}\_2, \mathbb{C}\_1) = \frac{ISL\_i(\mathbb{C}\_2)}{ISL\_i(\mathbb{C}\_1)} \quad (i = 1, 2, 3) \tag{28}$$

$$f(\mathbb{C}\_2, \mathbb{C}\_1) = \frac{1}{3} \sum\_{i=1}^3 \frac{w\_1 RIARE\_i(\mathbb{C}\_2, \mathbb{C}\_1) + w\_2 RISL\_i(\mathbb{C}\_2, \mathbb{C}\_1)}{w\_1 + w\_2} \tag{29}$$

where *ussi* is the steady state value of *i*th input; *C*1,*C*<sup>2</sup> are the two compared controllers; the weighting factors *w*<sup>1</sup> and *w*<sup>2</sup> in equation (29) are chosen as *w*<sup>1</sup> = *w*<sup>2</sup> = 0.5.

The combined index *J* for the reference tracking case is 0.6564. From Figure 6, Tables 4 and 5, it shows that the IARE and ISU of fractional order MPC are both smaller than that of traditional MPC, which means the fractional order MPC can obtain better performance with less control effort changes. By choosing suitable fractional order terms, the fractional order MPC shows superiority compared with integer order MPC.

**Figure 6.** Outputs and inputs of the boiler-turbine system with FOMPC and MPC in the reference tracking experiment (the outputs are listed on the left hand, and the inputs are listed on the right hand). (**a**) control loop for drum steam pressure, (**b**) control loop for required power, (**c**) control loop for drum water level.


**Table 4.** Performance indexes for *IARE* and *ISU* in the power reference tracking experiment.

**Table 5.** Performance indexes for *RIARE* and *RISU* in the power reference tracking experiment (MPC is the *C*<sup>2</sup> and FOMPC the *C*<sup>1</sup> according to Equations (27) and (28)).


#### **5. Conclusions**

This paper proposed a fractional order model predictive controller for the boilerturbine system. Due to the nonlinearity and multiple variables of the boiler-turbine, the nonlinear MPC with distributed scheme is designed, and the termination conditions are given. The integer order cost function is replaced with the fractional order cost function, which simplified the configration of the weighting factor matrices in the cost function. The number of weighting factors required to be tuned decreases from *Np* + *Nc* to two. According to the simulation for power tracking, it is proved that the fractional order MPC improves the control performance compared with the traditional MPC method. In this work, better control performance is obtained with fractional order MPC; however, how the fractional order effects the control performance is not clear, which can be researched further in the future.

**Author Contributions:** Methodology, S.Z.; software, S.W.; validation, S.W., R.C. and W.R.; formal analysis, S.W., R.C. and B.L.; writing—original draft preparation, S.Z.; writing—review and editing, S.Z., S.W., R.C., W.R. and B.L.; supervision, S.Z.; project administration, S.Z. and B.L.; funding acquisition, S.Z. and B.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Natural Science Foundation of Heilongjiang under grant KY10400210217, Foundation of Fundamental Strengthening Program for Technical Field under grant 2021-JCJQ-JJ-0026, Fundamental Research Funds for the Central Universities undergrant 3072020CFT1501 and Foundation of High-level scientific research guidance project of Harbin Engineering University under grant 3072022QBZ0406.

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

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

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

#### **References**

