Optimal feedback control for a class of fed-batch fermentation processes using switched dynamical system approach

: This paper considers an optimal feedback control problem for a class of fed-batch fermentation processes.


Introduction
Over the past decades, the use of biochemical reactors and correlation techniques has increased greatly because of their fruitful application in converting biomass or cells into pharmaceutical or chemical products, such as vaccines [1], antibiotics [2], beverages [3], and industrial solvents [4]. Among various classes or operation regions of bioreactors, the fed-batch modes have extensively used in the biotechnological industry due to its considerable economic profits [5][6][7]. The main objective of these reactors is to achieve a given or maximum concentration of production at the end of the operation, which can be implemented by using some suitable feed rates [8][9][10]. Thus, in order to ensure economic benefit and product quality of the fed-batch processes, the process control of this units is an very important topic for the engineers [11][12][13].
Switched dynamical systems provide a flexible modeling method for a variety of different types of engineering systems, such as financial system [14], train control system [15], hybrid electric vehicle [16], chemical process system [17], and biological system [18][19][20][21]. Generally speaking, switched dynamical systems are formed by some continuous-time or discrete-time subsystems and a switching rule [22]. There usually exist four types of switching rules as follows: time-dependent switching [23], state-dependent switching [24], average dwell time switching [25], and minimum dwell time switching [26]. Recently, switched dynamical system optimal control problems are becoming increasingly attractive due to their significance in theory and industry production [27][28][29][30]. Because of the discrete nature of switching rules, it is very challenging that switched dynamical system optimal control problems are solved by directly using the classical optimal control approaches such as the maximum principle and the dynamic programming method [31][32][33][34]. In additions, analytical methods also can not be applied to obtain an solution for switched dynamical system optimal control problems due to their nonlinear nature [35][36][37]. Thus, in recent work, two kinds of well-known numerical optimization algorithms are developed for switched dynamical system optimal control problems to obtain numerical solutions. One is the bi-level algorithm [38,39]. The other is the embedding algorithm [40,41]. Besides above two kinds of well-known numerical optimization algorithms, many other available numerical optimization algorithms are also developed for obtaining the solution of switched dynamical system optimal control problems [42]. Unfortunately, most of these numerical optimization algorithms depend on the following assumption: the time-dependent switching strategy is used to design the switching rules, which implies that the system dynamic must be continuously differentiable with respect to the system state [43][44][45]. However, this assumption is not reasonable, since some small perturbations of the system state may lead to the dynamic equations being changed discontinuously. Thus, the solution obtained is usually not optimal. In additions, although these approaches have demonstrated to be effective by solving many practical problems, they only obtaining an open loop control [46][47][48][49][50][51][52][53]. Unfortunately, such open loop controls are not usually robust in practice. Thus, an optimal feedback controller is more and more popular.
In this paper, we consider an optimal feedback control problem for a class of fed-batch fermentation processes by using switched dynamical system approach. Our main contributions are as follows. Firstly, a dynamic optimization problem for a class of fed-batch fermentation processes is modeled as a switched dynamical system optimal control problem, and a general state-feedback controller is designed for this dynamic optimization problem. Unlike the existing works, the state-dependent switching method is applied to design the switching rule, and the structure of this state-feedback controller is not restricted to a particular form. In generally, the traditional methods for obtaining an optimal feedback control require solving the well-known Hamilton-Jacobi-Bellman partial differential equation, which is a very difficult issue even for unconstrained optimal control problems. Then, in order to overcome this difficulty, this problem is transformed into a mixed-integer optimal control problem by introducing a discrete-valued function. Furthermore, each of these discrete variables is represented by using a set of 0-1 variables. Then, by using a quadratic constraint, these 0-1 variables are relaxed such that they are continuous on the closed interval [0, 1]. Accordingly, the original mixed-integer optimal control problem is transformed into a nonlinear parameter optimization problem, which can be solved by using any gradient-based numerical optimization algorithm. Unlike the existing works, the constraint introduced for these 0-1 variables are at most quadratic. Thus, it does not increase the number of locally optimal solutions of the original problem. During the past decades, many iterative approaches have been proposed for solving the nonlinear parameter optimization problem by using the information of the objective function. The idea of these iterative approaches is usually that a iterative sequence is generated such that the corresponding objective function value sequence is monotonically decreasing. However, the existing algorithms have the following disadvantage: if an iteration is trapped to a curved narrow valley bottom of the objective function, then the iterative methods will lose their efficiency due to the target with objective function value monotonically decreasing may leading to very short iterative steps. Next, in order to overcome this challenge, an improved gradient-based algorithm is developed based on a novel search approach. In this novel search approach, it is not required that the objective function value sequence is always monotonically decreasing. And a large number of numerical experiments shows that this novel search approach can effectively improve the convergence speed of this algorithm, when an iteration is trapped to a curved narrow valley bottom of the objective function. Finally, an optimal feedback control problem of 1, 3-propanediol fermentation processes is provided to illustrate the effectiveness of this method developed by this paper. Numerical simulation results show that this method developed by this paper is low time-consuming, has faster convergence speed, and obtains a better result than the existing approaches.
The rest of this paper is organized as follows. Section 2 presents the optimal feedback control problem for a class of fed-batch fermentation processes.
In Section 3, by introducing a discrete-valued function and using a relaxation technique, this problem is transformed into a nonlinear parameter optimization problem, which can be solved by using any gradient-based numerical optimization algorithm. An improved gradient-based numerical optimization algorithm are developed in Section 4. In Section 5, the convergence results of this numerical optimization algorithm are established. In Section 6, an optimal feedback control problem of 1, 3-propanediol fermentation processes is provided to illustrate the effectiveness of this algorithm developed by this paper.

Problem formulation
In this section, a general state-feedback controller is proposed for a class of fed-batch fermentation process dynamic optimization problems, which will be modeled as an optimal control problem of switched dynamical systems under state-dependent switching.
Let α 1 = α 11 , · · · , α 1r 1 T ∈ R r 1 and α 2 = α 21 , · · · , α 2r 2 T ∈ R r 2 be two parameter vectors satisfying and b j α 2r j b j , j = 1, · · · , r 2 , (2.2) respectively, where a i ,ā i , i = 1, · · · , r 1 ; b j ,b j , j = 1, · · · , r 2 present given constants. Suppose that t f > 0 presents a given terminal time. Then, a class of fed-batch fermentation process dynamic optimization problems can be described as choose two parameter vectors α 1 ∈ R r 1 , α 2 ∈ R r 2 , and a general statefeedback controller to minimize the objective function subject to the switched dynamical system under state-dependent switching with the initial condition where x (t) ∈ R n presents the system state; x 0 presents a given initial system state; u (t) ∈ R m presents the control input; ϑ = ϑ 1 , · · · , ϑ r 1 T ∈ R r 3 presents a state-feedback parameter vector satisfying c k ϑ k c k , k = 1, · · · , r 3 , (2.7) c k andc k , k = 1, · · · , r present given constants. Υ : R n × R r → R m ; φ : R n → R, f 1 : R n × 0, t f → R n , f 2 : R n × R m × 0, t f → R n , g 1 : R n × R r 1 × 0, t f → R n , g 2 : R n × R r 2 × 0, t f → R n present five continuously differentiable functions. For convenience, this problem is called as Problem 1. Remark 1. In the switched dynamical system (2.5), Subsystem 1 presents the batch mode, during which there exists no input feed (i.e., control input) u (t), and Subsystem 2 presents the feeding mode, during which there exists input feed (i.e., control input) u (t). This fed-batch fermentation process will oscillate between Subsystem 1 (the batch mode) and Subsystem 2 (the feeding mode), and g 1 (x (t) , α 1 , t) = 0 and g 2 (x (t) , α 2 , t) = 0 present the active conditions of Subsystems 1 and 2, respectively. Remark 2. Note that an integral term, which is used to measure the system running cost, can be easily incorporated into the objective function (2.4) by augmenting the switched dynamical system (2.5) with an additional system state variable (see Chapter 8 of this work [54] ). Thus, it is not a serious restriction that the integral term does not appear in the objective function (2.4). Remark 3. The structure for this general state-feedback controller (2.3) can be governed by the given continuously differentiable function Υ, and the state-feedback parameter vector ϑ is decision variable vector, which will be chosen optimally. For example, the linear state-feedback controller described by u (t) = K x (t) is a very common state-feedback controller, where K ∈ R m×n presents a state-feedback gain matrix to be found optimally.

Problem transformation
In Problem 1, the state-dependent switching strategy is adopted to design the switching rule, which is unlike the existing switched dynamical system optimal control problem. Then, the solution of Problem 1 can not be obtained by directly using the existing numerical computation approaches for switched dynamical systems optimal control problem, in which the switching rule is designed by using time-dependent strategy. In order to overcome this difficulty, by introducing a discrete-valued function, the problem will be transformed into a equivalent nonlinear dynamical system optimal control problem with discrete and continuous variables in this subsection.
Firstly, by substituting the general state-feedback controller (2.3) into the switched dynamical system (2.5), Problem 1 can be equivalently written as the following problem: subject to the switched dynamical system under state-dependent switching and the three bound constraints (2.1), (2.2) and (2.7), wheref 2 Next, note that the solution of Problem 1 can not be obtained by directly using the existing numerical computation approaches for switched dynamical systems optimal control problem, in which the switching rule is designed by using time-dependent strategy and not state-dependent strategy. In order to overcome this difficulty, a novel discrete-valued function y (t) is introduced as follows: Then, Problem 2 can be transformed into the following equivalent optimization problem with discrete and continuous variables: subject to the nonlinear dynamical system the equality constraint and the three bound constraints (2.1), (2.2), and (2.7).

Problem relaxation
Note that standard nonlinear numerical optimization algorithms are usually developed for nonlinear optimization problems only with continuous variables, for example the sequential quadratic programming algorithm, the interior-point method, and so on. Thus, the solution of Problem 3, which has discrete and continuous variables, can not be obtained by directly using these existing standard algorithms. In order to overcome this difficulty, this subsection will introduce a relaxation problem, which has only continuous variables. Define where Then, a theorem can be established as follows. Theorem 1. If the nonnegative functions σ 1 (t) and σ 2 (t) satisfy the following equality: then two results can be obtained as follows: (1) For any t ∈ 0, t f , the function P (σ (t)) is nonnegative; (1) By using the equality (3.8) and the Cauchy-Schwarz inequality, we have Note that the functions σ 1 (t) and σ 2 (t) are nonnegative. Then, squaring both sides of the inequality (3.9) yields which implies that for any t ∈ 0, t f , the function P (σ (t)) is nonnegative.
(2) The correctness of the second part for Theorem 1 only need to prove the following result: for any t ∈ 0, t f , P (σ (t)) = 0 has solutions σ i * (t) = 1 for one i * ∈ {1, 2} and σ i (t) = 0 for the other i ∈ {1, 2} and i i * . Define Then, the inequality (3.9) can be equivalently transformed into as follows: where · and · present the vector dot product and the Euclidean norm, respectively. Note that the holds if and only if there exists a constant β ∈ R such that v 1 (t) = βv 2 (t) . (3.12) By using the equality (3.8), one obtain v 1 (t) 0 and v 2 (t) 0, where 0 presents the zero vector. Then, β is a nonzero constant and the equality (3.12) implies Furthermore, the constant β can be set equal to one integer i * ∈ {1, 2}, and for the other From the two equalities (3.8) and (3.15), we obtain σ i * (t) = 1. This completes the proof of Theorem 1. Now, Problem 3 can be rewritten as a relaxation problem as follows: subject to the nonlinear dynamical system the two equality constraints By using Theorem 1, one can derive that Problems 3 and 4 are equivalent.

A nonlinear parameter optimization problem
Note that the bound constraint (3.20) is essentially some continuous-time inequality constraints. Thus, the solution of Problem 4 can not also be obtained by directly using the existing standard algorithms. In order to obtain the solution of Problem 4, this subsection will introduce a nonlinear parameter optimization problem, which has some continuous-time equality constraints and several bound constraints.
Suppose that τ i presents the ith switching time. Then, one have where M 1 presents a given fixed integer. It is important to note that the switching times are not independent optimization variables, whose values can be obtained indirectly by using the state trajectory of the switched dynamical system (2.5). Then, Problem 4 can be transformed into an equivalent optimization problem as follows: subject to the nonlinear dynamical system 28) and the three bound constraints (2.1), (2.2), and (2.7), where ξ i 1 and ξ i 2 present, respectively, the values of σ 1 (t) and and χ I (t) is given by which is a function defined on the subinterval I ⊂ 0, t f . Due to the switching times being unknown, it is very challenging to acquire the gradient of the objective function (3.23). In order to overcome this challenge, the following time-scaling transformation is developed to transform variable switching times into fixed times: Suppose that the function t (s) : [0, M] → R is continuously differentiable and is governed by the following equation: with the boundary condition t (0) = 0, (3.31) where θ i is the subsystem dwell time on the ith subinterval [i − 1, i) ⊂ 0, t f . In general, the transformation (3.30)-(3.31) is referred to as a time-scaling transformation.
Then, by using the time-scaling transform (3.30) and (3.31), we can rewrite Problem 5 as the following equivalent nonlinear parameter optimization problem, which has fixed switching times .
subject to the nonlinear dynamical system

An improved gradient-based numerical optimization algorithm
In this section, an improved gradient-based numerical optimization algorithm will be proposed for obtaining the solution of Problem 1.

A penalty problem
In order to handle the continuous-time equality constraints (3.35) and (3.36), by adopting the idea of l 1 penalty function [55], Problem 6 will be written as a nonlinear parameter optimization problem with a linear equality constraint and several simple bounded constraints in this subsection. Problem 7. Choose (α 1 , α 2 , ϑ, ξ, θ) ∈ R r 1 × R r 2 × R r 3 × R 2M × R M to minimize the objective function where γ > 0 presents the penalty parameter. The idea of l 1 penalty function [47] indicates that any solution of Problem 7 is also a solution of Problem 6. In additions, it is straightforward to acquire the gradient of the linear function in the equality constraint (3.27), and the gradient of the objective function (4.1) will be presented in Section 4.2. Thus, the solution of Problem 7 can be achieved by applying any gradient-based numerical computation method.

Remark 4.
During the past decades, many iterative approaches have been proposed for solving the nonlinear parameter optimization problem by using the information of the objective function [57]. The idea of these iterative approaches is usually that a iterative sequence is generated such that the corresponding objective function value sequence is monotonically decreasing. However, the existing algorithms have the following disadvantage: if an iteration is trapped to a curved narrow valley bottom of the objective function, then the iterative methods will lose their efficiency due to the target with objective function value monotonically decreasing may leading to very short iterative steps. Then, in order to overcome this challenge, an improved gradient-based algorithm is developed based on a novel search approach in Algorithm 1. In this novel search approach, it is not required that the objective function value sequence is always monotonically decreasing. In additions, an improved adaptive strategy for the memory element N k described by (4.12), which is used in (4.13), is proposed in iterative processes in Algorithm 1. The corresponding explanation on the equality (4.12) is as follows. If the 1st condition described by (4.12) holds, then it implies that the iteration is trapped to a curved narrow valley bottom of the objective function. Thus, in order to avoid creeping along the bottom of this narrow curved valley, the value of the memory element N k should be increased. If the 2nd condition described by (4.12) holds, then the value of the memory element N k is better to remain unchanged. If the 3rd condition described by (4.12) holds, then it implies that the iteration is in a flat region. Thus, in order to decrease the objective function value, the value of the memory element N k will be decreased. Above discussions imply that the novel search approach described in Algorithm 1 is also an adaptive method. Remark 5. The sufficient descent condition is extremely important for the convergence of any gradientbased numerical optimization algorithm. Thus, the goal of lines 12-16 described in Algorithm 1 is avoiding uphill directions and keeping {ρ k } uniformly bounded. As a matter of fact, for any k, ρ min ρ k ρ max and d k = −ρ k g k ensure that there are two constants l 1 > 0 and l 2 > 0 such that d k satisfies the following two conditions:

Convergence analysis
This section will establish the convergence results of Algorithm 1 developed by Section 4. In order to establish the convergence results of this algorithm, we suppose that the following two conditions hold: Assumption 1. J γ is a continuous differentiable function and bounded below on R r 1 × R r 2 × R r 3 × R 2M × R M . Assumption 2. For any η 1 ∈ Ω and η 2 ∈ Ω, there is a constant l 3 such that where Ω presents a open set and g (η) presents the gradient of J γ (η). Theorem 3. Suppose that Assumptions 1 and 2 hold. Let {η k } be a sequence obtained by using Algorithm 1. Then, there is a constant ς > 0 such that the following inequality holds: If ς 0 > 0, then by using the inequalities (4.14) and (4.15), one can obtain Let ς be defined by ς = l 1 ς 0 . Then, the proof of Theorem 1 is complete for ς 0 > 0. If ς 0 = 0, then there is a subset Λ ⊆ {0, 1, 2, · · ·} such that the following equality holds: lim k∈Λ, k→∞ ω k = 0, (5.4) which indicates that there exists ak such that the following inequality holds: for any k >k and k ∈ Λ. Let ω = ω k . Then, the inequality (4.14) does not hold. That is, one can obtain Applying the mean value theorem to the left-hand side of the inequality (5.7) yields where 0 ζ k 1. From the inequality (5.8), one obtain (g (η k + ζ k ωd k )) T d k > µ (g k ) T d k .
Next, by using 0 q (Ap) A, the inequality (5.2), and the inequality (5.16), one can derive 17) which implies that the inequality (5.14) is also true for j + 1. Then, the inequality (5.14) is true for 1 j A by using mathematical induction. Thus, the inequality (5.12) holds.
In additions, Assumption 1 shows J γ being a continuous differentiable function and bounded below on R r 1 × R r 2 × R r 3 × R 2M × R M , which indicates that Then, summing the inequality (5.12) over p yields Thus, the inequality (5.13) holds. This completes the proof of Lemma 1. Theorem 4. Suppose that these conditions of Theorem 3 are true. Then, the following equality holds: where g (η k ) presents the gradient of the objective function J γ described by (4.1) at the point η k . Proof. Firstly, the following result will be proved: there is a constant l 4 such that By using Assumptions 1 and 2, one can obtain

Numerical results
In this section, an optimal feedback control problem of 1, 3-propanediol fermentation processes is provided to illustrate the effectiveness of the approach developed by Sections 2-5, and the numerical simulations are all implemented on a personal computer with Intel Pentium Skylake dual core processor i5-6200U CPU(2.3GHz).
The 1, 3-propanediol fermentation process can be described by switching between two subsystem: batch subsystem and feeding subsystem. There exists no input feed during the batch subsystem, while alkali and glycerol will be added to the fermentor during the feeding subsystem. In generally, the subsystem switching will happen, if the glycerol concentration reaches the given upper and lower thresholds. By using the result of the work [58], the 1, 3-propanediol fermentation process can be modeled as the following switched dynamical system under state-dependent switching: where t f denotes the given terminal time; the system states x 1 (t), x 2 (t), x 3 (t), x 4 (t) denote the volume of fluid (L), the concentration of biomass (gL −1 ), the concentration of glycerol (mmolL −1 ), the concentration of 1,3-propanediol (mmolL −1 ), respectively; the control input u (t) denotes the feeding rate (Lh −1 ); ] T denotes the system state vector; Subsystem 1 and Subsystem 2 denote the batch subsystem and the feeding subsystem, respectively; α 1 and α 2 (two parameters that need to be optimized) denote the upper and lower of the glycerol concentration, respectively; and the functions f 1 (x (t) , t), f 2 (x (t) , u (t) , t) are given by Subsystem 1 is essentially a natural fermentation process due to no input feed. The functions ϕ, ∆ 1 , and ∆ 2 are defined by which denote the growth rate of cell, the consumption rate of substrate, and the formation rate of 1,3-propanediol, respectively. In the equality (6.4), the parameters x * 3 and x * 4 denote the critical concentrations of glycerol and 1,3-propanediol, respectively; h 1 , h 2 , h 3 , Y 1 , Y 2 , Y 3 , Z 1 , Z 2 , l 7 , and l 8 are given parameters. Note that the feeding subsystem doesn't only consist of the natural fermentation process. Thus, the function f 2 (x (t) , u (t) , t) is provided to describe the process dynamics because of the control input feed in Subsystem 2. In the equality (6.3), the given parameters l 5 and l 6 denote the proportion and concentration of glycerol in the control input feed, respectively.
In generally, as the increase of the biomass, the consumption of glycerol also increases. Then, during Subsystem 1 (batch subsystem), the concentration of glycerol will eventually become too low due to no new glycerol being added. Thus, Subsystem 1 will switch to Subsystem 2 (feeding subsystem), when the equality x 3 (t) − α 2 = 0 (the active condition of Subsystem 2) satisfies. On the other hand, during Subsystem 2 (feeding subsystem), the concentration of glycerol will eventually become too high due to new glycerol being added. This will inhibit the growth of cell. Thus, Subsystem 2 will switch to Subsystem 1 (batch subsystem), when the equality x 3 (t) − α 1 = 0 (the active condition of Subsystem 1) satisfies.
Suppose that the feeding rate u (t), the upper of the glycerol concentration α 1 , and the lower of the glycerol concentration α 2 satisfy the following bound constraints: Suppose that the control input u (t) takes the piecewise state-feedback controller Our main objective is to maximize the concentration of 1,3-propanediol at the terminal time t f . Thus, the optimal feedback control problem of 1, 3-propanediol fermentation processes can be presented as follows: choose a control input u (t) to minimize the objective function J (u (t)) = −x 4 (t) subject to the switched dynamical system described by (6.1) with with the initial condition x (0) = x 0 and the bound constraints (6.7-6.9). Then, the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3) is adopted to solve the optimal feedback control problem of 1, 3-propanediol fermentation processes by using Matlab 2010a. The optimal objective function value is J * = −x 4 t f = −1265.5597 and the optimal values of the parameters α 1 and α 2 are 584.3908 and 246.5423, respectively. The optimal feedback gain matrixes K * i , i = 1, · · · , 9 are presented by Note that Problem 6 is an optimal control problem of nonlinear dynamical systems with state constraints. Thus, the finite difference approximation approach developed by Nikoobin and Moradi [59] can also be applied for solving this dynamic optimization problem of 1, 3-propanediol fermentation processes. In order to compare with the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3), the finite difference approximation approach developed by Nikoobin and Moradi [59] is also adopted for solving this dynamic optimization problem of 1, 3-propanediol fermentation process with the same model parameters under the same condition, and the numerical comparison results are presented by Figure 5 and Table 1. The improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4. 3) The finite difference approximation approach developed by Nikoobin and Moradi [59] Figure 5. Convergence rates for the finite difference approximation approach developed by Nikoobin and Moradi [59] and the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3).  Figure 5 shows that the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3) takes only 67 iterations to obtain the satisfactory result x 4 (t f ) = 1265.5597, while the finite difference approximation approach developed by Nikoobin and Moradi [59] takes 139 iterations to achieve the satisfactory result x 4 (t f ) = 1052.9140 . That is, the iterations of the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3) is reduced by 51.7986%. In additions, Table 1 also shows that the result x 4 (t f ) = 1052.9140 obtained by using the finite difference approximation approach developed by Nikoobin and Moradi [59] is not superior to the result (x 4 (t f ) = 1265.5597) obtained by using the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3) with saving 60.4695% computation time.
In conclusion, the above numerical simulation results show that the improved gradient-based numerical optimization algorithm (Algorithm 1 described by Section 4.3) is low time-consuming, has faster convergence speed, and can obtain a better numerical optimization than the finite difference approximation approach developed by Nikoobin and Moradi [59]. That is, an effective numerical optimization algorithm is presented for solving the dynamic optimization problem of 1, 3-propanediol fermentation process.

Conclusions
In this paper, the dynamic optimization problem for a class of fed-batch fermentation processes is modeled as an optimal control problem of switched dynamical systems under state-dependent switching, and a general state-feedback controller is designed for this dynamic optimization problem. Then, by introducing a discrete-valued function and using a relaxation technique, this problem is transformed into a nonlinear parameter optimization problem. Next, an improved gradient-based algorithm is developed based on a novel search approach, and a large number of numerical experiments show that this novel search approach can effectively improve the convergence speed of this algorithm, when an iteration is trapped to a curved narrow valley bottom of the objective function. Finally, an optimal feedback control problem of 1, 3-propanediol fermentation processes is provided to illustrate the effectiveness of this method developed by this paper, and the numerical simulation results show that this method developed by this paper is low time-consuming, has faster convergence speed, and obtains a better result than the existing approaches. In the future, we will continue to study the dynamic optimization problem for a class of fed-batch fermentation processes with uncertainty constraints.