Multidimensional Taylor Network Optimal Control of MIMO Nonlinear Systems without Models for Tracking by Output Feedback.
In most input-output control systems, an alteration of one input signal may trigger the change of multiple output signals [1, 2]. However, in the actual control process, the system's stability is highly demanded [3-5], along with its desirable dynamic performance and multiple tracking with less or even zero errors . Because of the complex internal structure of the multi-input and multioutput (MIMO) system, satisfactory results are hardly gained as a result .
Support vector machine stemming from statistical theory is good at handling classification and regression problems. Its great progress has been observed in the fields of nonlinear control , fault diagnosis , and process modeling . However, at present, it is only applicable to the single variable system due to its complexity.
Recently, great advances for adaptive output feedback control and robust control of nonlinear systems have been achieved [11-13]. Good approximation property of the neural networks provides a new access to adaptive output feedback control [14-16]. For general MIMO nonlinear systems, especially the systems without models, the control methods are always based on neural networks [17-19]. The neuron of neural networks includes the exponential function, which contributes to the complexity of calculation, making the neural network control unable to meet the real-time requirements. Though the larger the number of neural network nodes, the smaller the approximation error, neural network with a large hidden node number tends to complicate the control and the computation of the control system. So, it is difficult to realize the system's real-time control.
Polynomial control method is an intuitive controller design approach. Its design parameters are of clear physical sense . However, it requires an accurate model of the controlled object for solving MIMO nonlinear control problems.
The multidimensional Taylor network (MTN, whose idea was proposed by Hong-Sen Yan in 2010 and realization was completed by Bo Zhou who is Yan's Ph.D. student) is commonly applied to the analysis of time series prediction [21-28]. MTN, a simple function of the state and input, is good at analyzing and solving the problems in point due to its polynomials. Besides, MTN only involves multiplication and addition, so its simple computation makes desirable real-time control possible. In fact, its computation complexity is nearly equal to that of Taylor expansion of a single neuron of neural network. The idea of MTN optimal control was proposed by Yan . The optimal adjustment controller based on MTN was then developed . However, the parameters of the controller are fixed. The graduate students supervised by Yan have used the simple MTN to solve some practical control engineering problems successfully, such as the plane flight , the tank firing control in high speed motion , the axisymmetric cruise missile flight for attacking static target [33, 34], and ship roll stabilization , all involving MIMO nonlinear systems with strong disturbance. Simulation results show that the simple MTN promises faster response speed, stronger anti-interference capability, and better external stability than PID [31-35], PID neural network [31,33], neural network , sliding mode control [31, 34, 35], and active disturbance rejection controller .
As has already been pointed out, the existing control methods are known to have some shortcomings, such as needing an accurate model of the controlled object, mainly aiming at the SISO nonlinear system. Some MTN control strategies have also been developed. However, they are just based on the simple MTN (i.e., PID plus the sum of their second-order monomials plus PID, each item of which is multiplied by its corresponding parameter), and the PID parameters are chosen as the initial parameters of MTN controller via the minimum principle, mainly targeting the SISO nonlinear systems. Proposed in the present study is the MIMO feedback control, which, unlike other approaches, does not require the specific format of MTN. Each input of the controlled object corresponds to a subset of internal weights of Multidimensional Taylor Network Controller (MTNC). Even without knowing the internal characteristic parameters of the controlled object, and only by adjusting the internal weights of MTNC, the outputs of the closed-loop system can be made to track the desired signals effectively. Superior to the BP neural network self-adaption reconstitution controller, the MTNC tracks the expectation output curves more satisfactorily as well as suppressing the disturbance more efficiently, promising a better real-time dynamic performance.
The main contributions of this paper are listed as follows:
(1) For the first time, the MTN approach is used to solve the tracking control problem of MIMO nonlinear systems. The computation burden is released due to the simple structure of MTN.
(2) Even with modeling errors, the MTNC proposed in this paper can still be able to satisfy certain performance indexes with excellent dynamic performance.
(3) Compared with the existing control approaches, the proposed control method is simple and accurate.
2. Problem Statement
Consider the following MIMO nonlinear system:
x (k + 1) = f (x (k), u (k)) y (k) = h (x (k), u (k)). (1)
After differential homeomorphic transformation, (1) can be converted into (2) . Without losing generality, consider the following MIMO nonlinear system:
[mathematical expression not reproducible], (2)
where [[r.sub.1], [r.sub.2], [r.sub.m]] is the relative order of the system, and [r.sub.1] + [r.sub.2] +... + [r.sub.m] = n; y = [[y.sub.1], [y.sub.2]..., [y.sub.m]].sup.T] [member of] [R.sup.m] is the output vector of the system; [mathematical expression not reproducible] is the state vector of the system; u = [[u.sub.1], [u.sub.2], ..., [u.sub.m]].sup.T] [member of] [R.sup.m] is the input vector of the system; [f.sub.i] (x(k)), i = 1, 2, ..., m is the unknown bounded nonlinear mapping; [g.sub.ij] (x(k)), i, j = 1, 2, ..., m is the known bounded nonlinear mapping.
Equation (1) can be rewritten as
[mathematical expression not reproducible], (3)
where [mathematical expression not reproducible] is the control gain matrix; d = [[[d.sub.1], [d.sub.2], ..., [d.sub.m]].sup.T] is the unknown and bounded disturbance.
From the concept of the dominant input, the m x m MIMO system can be seen as m multi-input single output system, and the ith subsystem is
[mathematical expression not reproducible]. (4)
Select a dominant input in the m input, and set [u.sub.i](k). The rest is treated as interference; then we have
[mathematical expression not reproducible], (5)
[mathematical expression not reproducible], (6)
and (4) is equivalent to
[mathematical expression not reproducible]. (7)
Assuming only output can be measured, and
[g.sub.ii] (x (k)) [greater than or equal to] [g.sub.0] > 0,
[absolute value of dsi] [less than or equal to] D, where D is an unknown constant,
the expected output of the system is
[mathematical expression not reproducible], (9)
the desired state vector of the ith subsystem is
[mathematical expression not reproducible], (10)
and the tracking error is
[e.sub.i] (k) = [y.sub.di], (k) - [y.sub.i] (k). (11)
To prove that the tracking error converges to 0, set the tracking error vector
[e.sub.i] (k) = [x.sub.di] (k) - [x.sub.i] (k). (12)
Control target: for the affine nonlinear MIMO systems, without state observation, design the adaptive MTN optimal output feedback controller to make the nonlinear system (2) stable, enabling its output to track the desired signal and its parameters and tracking error to be uniformly bounded.
Define the filtering error as
[s.sub.i] (k) = [LAMBDA] * [e.sub.i] (k) = [[tau], [tau], ...,[tau]] * [e.sub.i] (k),
i = 1, 2, ... m, (13)
where [LAMBDA] = [[tau], [tau], ..., [tau]], [tau] > 0 is the design constant.
Theorem 1. Assuming system (5) with [f.sub.i] (x(k)) and [g.sub.ii](x(k)) is known and [d.sub.si] = 0, the tracking error then converges to 0 if the optimal control input is set as
[u.sup.*.sub.i] (k) = Ts (k) * [e.sub.i] (k) + [u.sup.*.sub.in] (k), (14)
where [u.sup.*.sub.in] (k) = ([y.sub.di](k + [r.sub.i] - [f.sub.i](x(k)) - [e.sub.i](k))/[g.sub.ii](x(k)), Ts(k) is the design constant, and Ts(k) * [e.sub.i] (k) * [g.sub.ii](x(k)) > 0.
Proof. Taking the Lyapunov function candidate as
[s.sub.i] (k) = [LAMBDA] * [e.sub.i] (k) = [[tau], [tau], ..., [tau]]* [e.sub.i] (k), (15)
[mathematical expression not reproducible]. (16)
Substituting [u.sup.*.sub.i] (k) and [d.sub.si] = 0 into (7) yields
[DELTA][s.sub.i] (k) = -Ts (k) * [e.sub.i] (k) * [g.sub.ii] (x (k)) < 0. (17)
According to the Lyapunov theorem, this result means [lim.sub.k[right arrow][infinity]] [s.sub.i](k) = 0.
That completes the proof.
However, [d.sub.si] [not equal to] 0. So other methods are needed to suppress the interference. As [f.sub.i] (x(k)) is unknown, [u.sup.*.sub.in] (k) is unobtainable. Thus, the MTN is used to approximate [u.sup.*.sub.in] (k), and [u.sub.si](k) is added to interference suppression.
3. Multidimensional Taylor Network
MTN can approximate any nonlinear functions with the finite point of discontinuity. Compared with the existing methods based on neural networks, MTN has the following merits: (1) being neatly structured; (2) being good at representing or approximating to nonlinear dynamical systems; (3) guaranteeing real-time control by only addition or multiplication operations allowable. In addition, as with the neural networks, only the internal weights of MTNC are required to be adjusted to make the outputs of the closed-loop system track the desired signals effectively.
By the function approximation of MTN, we obtain
[mathematical expression not reproducible], (18)
where approximation error [epsilon] satisfies the following conditions:
[absolute value of [epsilon]] [less than or equal to] [[epsilon].sub.0], (19)
where [[epsilon].sub.0] is the given normal number, and
= [[[e.sub.1] (k), [e.sub.i] (k-1), ..., [e.sub.i] (k), ..., [e.sub.m] (k),...].sup.T]. (20)
For the convenience of writing, denote the number of elements of [eta](k) as [n.sub.z], and we have
[mathematical expression not reproducible]. (21)
The basic structure of MTN is shown in Figure 1.
The output of MTN [u.sub.in](k) can be expressed as
[mathematical expression not reproducible], (22)
where N([n.sub.z], t) is the total number of the expansions, [w.sub.ij] (k) is the weight of the jth product term, [[lambda].sub.(s,j)] is the power of [z.sub.s] (k) in the jth product term, and [[summation].sup.n.sub.s=1] [[lambda].sub.s,j] [less than or equal to] t.
Setting [mathematical expression not reproducible] gives
[u.sup.*.sub.in](k) = [W.sup.*.sub.i] (k) * [THETA] ([eta]/(k)) + [epsilon]. (23)
The optimal parameter [W.sup.*.sub.i] (k) is defined as
[mathematical expression not reproducible], (24)
where [mathematical expression not reproducible] is the design constant.
Then the MTN output feedback controller can be written as
[mathematical expression not reproducible], (25)
where [mathematical expression not reproducible] is the output of the MTN, [[??].sub.i] (k) is the estimated value of the optimal weight [W.sup.*.sub.i] (k), and [u.sub.si] is the robust control term:
[u.sub.si] (k) = D + [g.sub.0] x [[epsilon].sub.0]/[g.sup.0] x sgn ([s.sub.i](k)). (26)
The control structure is shown in Figure 2.
4. Parameter Adjustment
To ensure that the output [y.sub.i](k) of the controlled object follows the desired input [r.sub.i](k), the parameters [W.sub.i](k) of MTN controller are adjusted by back-propagation learning algorithm.
Specific adjustments go as follows:
[mathematical expression not reproducible]. (27)
And we get
[w.sub.ij] (k+1) = [w.sub.ij] (k) - [[alpha].sub.j] (k) [partial derivative]E/[w.sub.ij](k) (28)
[mathematical expression not reproducible]. (29)
The system structure is unknown. However, the influence of the real system by output [u.sub.i](k) of controller can be obtained, so [partial derivative][y.sub.i](k)/[partial derivative][u.sub.i](k) can be replaced by the sign function
sgn [[y.sub.i] (k+1) - [y.sub.i] (k)/[u.sub.i]](k) - [u.sub.i](k-1)]. (30)
This alternative is feasible, as the numerical size of [partial derivative][y.sub.i] (k)/ [partial derivative][y.sub.i] (k), as an intermediate product term, can be adjusted by step [[alpha].sub.i](k). But it does not change the direction of the adjustment. Simultaneously, it can prevent the numerical size from going too large by the two-step output of the MTN convergence.
Step [[alpha].sub.i](k) uses the following adaptive adjustment strategy:
[[alpha].sub.j] (0) = [[alpha].sub.j0], (31)
where [[alpha].sub.j0] is a given constant as the initial value,
[mathematical expression not reproducible], (32)
where [beta] > 1, 0 < [gamma] < 1, and [parallel][partial derivative]E/[w.sub.ij](k)[parallel] is the norm of [partial derivative]E/[w.sub.ij](k).
5. Simulation Results
For demonstration of the effectiveness of the proposed method, consider the following nonlinear system. No specific parameters of the controlled object are needed in the control process.
[mathematical expression not reproducible], (33)
where Ts = 0.1.
The initial parameters of the system are
[mathematical expression not reproducible]. (34)
5.1. Step Response. For the proposed MTNC, when the dimension n of the controller is equal to 2, the unit step response curve is shown in Figure 3.
The unit step response curve of the neural network self-adaption reconstitution controller is presented in Figure 3, which reveals that the BP neural network self-adaption reconstitution controller responds slowly, while the MTNC has a much faster response speed. Figure 4 illustrates the absolute tracking error curve.
5.2. Antijamming Capability. For demonstration of the antijamming capability of the controller, when k is greater than 60, a sinusoidal signal with an amplitude of 0.12 is superimposed on [x.sub.2] and [x.sub.4], respectively, and the response curve is shown in Figure 5.
From Figure 5, it can be seen that the MTNC is of higher antijamming capability.
The absolute error curve is displayed in Figure 6.
5.3. Tracking Ability. When k > 10, with the expected output curves [y.sub.d1] and [y.sub.d2] overlaying the cosine signal with an amplitude of 0.1 to verify the follow-up response performance of the controller, the response curves are shown in Figure 7.
As indicated by Figure 7, the proposed design approach tracks the desired signal more quickly.
5.4. Robustness. The amplitude of disturbance is expanded 5 times to verify the robustness of the controller, and the response curve is given in Figure 8.
Experimental results demonstrate that the MTNC based on the multidimensional Taylor network is of better robustness for disturbance.
5.5. Summary. In short, by realizing the optimal closed-loop tracking control of the MIMO nonlinear systems through output feedback with no model involved, the MTNC guarantees its faster response speed, robustness, and antidisturbance capability.
Besides, the MTNC, with good real-time capability and low computational complexity, does not contain the exponential function, so it is easy to achieve real-time control. Taking ten nodes as an example and supposing that the highest order of the exponential function expansion of the neural network is three, the controller of the neural network needs 98 additions and 160 multiplications, while the MTNC needs only 18 additions and 30 multiplications. Therefore, the reduction of calculation complexity greatly improves the real-time performance of the proposed controller.
The MIMO nonlinear system control method based on the multidimensional Taylor network has been proposed here. Compared with the traditional methods, MTN, as a kind of dynamic model, can be seen as a system model. In the control process, the closed-loop tracking control of nonlinear system was realized by output feedback without state observation. Due to the Lyapunov function, the tracking error turned out to be nearly zero. The MTNC parameter adjustment algorithm without the specific parameters of the controlled object was developed as well.
During the experiments, the MTNC tracked the expected output curves satisfactorily, capable of suppressing the disturbance efficiently when a sinusoidal disturbance was superimposed on [x.sub.2] and [x.sub.4]; the real-time dynamic performance of the controller has been significantly improved due to the reduction of calculation complexity.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Qi-Ming Sun and Hong-Sen Yan contributed equally to this work.
This work was supported in part by the National Natural Science Foundation of China under Grants 61673112 and 60934008, the Fundamental Research Funds for the Key Universities of China under Grants 2242017K10003 and 2242014K10031, and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). The authors thank Professor Li Lu for valuable comments and suggestions.
 L. Mo and M. M. Bayoumi, "Adaptive pole placement control of MIMO systems," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 38, no. 6, pp. 967-970, 1993.
 S. S. Xiong and Q. X. Zhu, "Decentralized risk-sensitive design for large-scale stochastic interconnected systems with time-varying delays," Journal of The Franklin Institute, vol. 353, no. 7, pp. 1527-1552, 2016.
 G. Pin and T. Parisini, "Networked predictive control of uncertain constrained nonlinear systems: recursive feasibility and input-to-state stability analysis," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 56, no. 1, pp. 72-87, 2011.
 B. Wang and Q. X. Zhu, "Stability analysis of Markov switched stochastic differential equations with both stable and unstable subsystems," Systems & Control Letters, vol. 105, pp. 55-61,2017
 Q. X. Zhu and Q. Y. Zhang, "Pth moment exponential stabilisation of hybrid stochastic differential equations by feedback controls based on discrete-time state observations with a time delay," IET Control Theory Applications, vol. 11, no. 12, pp. 1992-2003, 2017.
 A. Mannava, S. N. Balakrishnan, L. Tang, and R. G. Landers, "Optimal tracking control of motion systems," IEEE Transactions on Control Systems Technology, vol. 20, no. 6, pp. 1548-1558, 2012.
 Q. X. Zhu, "Razumikhin-type theorem for stochastic functional differential equations with Levy noise and Markov switching," International Journal of Control, vol. 90, no. 8, pp. 1703-1712, 2017.
 C. J. Ong, D. Sui, and E. G. Gilbert, "Enlarging the terminal region of nonlinear model predictive control using the support vector machine method," Automatica, vol. 42, no. 6, pp. 1011-1016, 2006.
 H. Yi, X. F. Song, B. Jiang, and D. C. Wang, "Support vector machine based on nodes refined decision directed acyclic graph and its application to fault diagnosis," Acta Automatica Sinica, vol. 36, no. 3, pp. 427-432, 2010.
 I. Goethals, K. Pelckmans, J. A. K. Suykens, and B. D. Moor, "Identification of MIMO Hammerstein models using least squares support vector machines," Automatica, vol. 41, no. 7, pp. 1263-1272, 2005.
 J. X. Li, C. Y. Gao, Y. M. Fang, C. Guo, and W. B. Zhang, "Adaptive regulation of a class of uncertain nonlinear systems with nonstrict feedback form," Mathematical Problems in Engineering, Article ID 9748489, pp. 1-7, 2016.
 N. Hovakimyan, B. J. Yang, and A. J. Calise, "Adaptive output feedback control methodology applicable to non-minimum phase nonlinear systems," Automatica, vol. 42, no. 4, pp. 513-522, 2006.
 H. Wang and Q. X. Zhu, "Finite-time stabilization of high-order stochastic nonlinear systems in strict-feedback form," Automatica, vol. 54, pp. 284-291, 2015.
 J. Vance and S. Jagannathan, "Discrete-time neural network output feedback control of nonlinear discrete-time systems in non-strict form," Automatica, vol. 44, no. 4, pp. 1020-1027, 2008.
 Z. Wu, J. Lu, J. Shi, Y. Liu, and Q. Zhou, "Robust adaptive neural control of morphing aircraft with prescribed performance," Mathematical Problems in Engineering, pp. 1-16, 2017
 Q. X. Zhu and J. D. Cao, "Exponential stability of stochastic neural networks with both Markovian jump parameters and mixed time delays," IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 41, no. 2, pp. 341-353, 2011.
 Q. Yang, Z. Yang, and Y. Sun, "Universal neural network control of MIMO uncertain nonlinear systems," IEEE Transactions on Neural Networks and Learning Systems, vol. 23, no. 7, pp. 1163-1169, 2012.
 Q. X. Zhu and J. D. Cao, "Stability analysis of markovian jump stochastic BAM neural networks with impulse control and mixed time delays," IEEE Transactions on Neural Networks and Learning Systems, vol. 23, no. 3, pp. 467-479, 2012.
 Q. X. Zhu, J. D. Cao, and R. Rakkiyappan, "Exponential input-to-state stability of stochastic Cohen-Grossberg neural networks with mixed delays," Nonlinear Dynamics, vol. 79, no. 2, pp. 1085-1098, 2015.
 C. Samia and S. Lassaad, "Polynomial approach to nonlinear-system control via the exact state-input linearization method," in Proceedings of the 16th International Conference on Sciences and Techniques of Automatic Control and Computer Engineering, STA 2015, pp. 905-911, tun, December 2015.
 Y. Lin, H. S. Yan, and B. Zhou, "Nonlinear time series prediction method based on multi-dimensional Taylor network and its applications," Control and Decision, vol. 29, no. 5, pp. 795-801, 2014, (in Chinese).
 Y. Lin, H. S. Yan, and B. Zhou, "A novel modeling method based on multi-dimensional Taylor network and its application in time series prediction," Advanced Materials Research, vol. 940, pp. 480-484, 2014.
 Y. Lin, H. S. Yan, and B. Zhou, "A novel nonlinear time series prediction method and its application in structural vibration response prediction," Applied Mechanics and Materials, vol. 599601, pp. 1918-1921, 2014.
 B. Zhou and H. S. Yan, "Financial time series forecasting based on wavelet and multi-dimensional Taylor network dynamics model," Systems Engineering--Theory and Practice, vol. 33, no. 10, pp. 2654-2662, 2013, (in Chinese).
 B. Zhou and H. S. Yan, "Time series forecasting based on the empirical mode decomposition multi-dimensional Taylor network model," in Proceedings of the 2013 9th International Conference on Natural Computation, ICNC 2013, pp. 1194-1198, chn, July 2013.
 B. Zhou and H. S. Yan, "Nonlinear system identification and prediction based on dynamics cluster multi-dimensional Taylor network model," Control and Decision, vol. 29, no. 1, pp. 33-38, 2014, (in Chinese).
 B. Zhou and H. S. Yan, "A dynamics model based on intermittent feedback multi-dimensional Taylor network," ACTA Automatica Sinica, vol. 40, no. 7, pp. 1517-1521, 2014, (in Chinese).
 H. S. Yan and A. M. Kang, "Asymptotic tracking and dynamic regulation of SISO nonlinear system based on discrete multidimensional Taylor network," in IET Control Theory and Applications, pp. 10-1049, 2017
 H. S. Yan, "MTN (multi-dimensional Taylor network) optimal control," Tech. Rep., Research Institute of Control and Optimization of Manufacturing Systems, School of Automation, Southeast University, Nanjing, China, 2010, (in Chinese).
 Q. M. Sun and H. S. Yan, "Optimal adjustment control of SISO nonlinear systems based on multi-dimensional taylor network only by output feedback," Advanced Materials Research, vol. 1049-1050, pp. 1389-1391, 2014.
 Y. Zhou, Research and simulation of multi-dimensional Taylor network optimal control for the flight of the plane [M.S. thesis], Southeast University, Dhaka, Bangladesh, 2016, (in Chinese).
 T. Jin, Research and simulation of multi-dimensional Taylor network optimal control for the tank firing control in high speed motion [M.S. thesis], Southeast University, Dhaka, Bangladesh, 2016, (in Chinese).
 L. C. Xia, Research into multi-dimensional Taylor network optimal control of the axisymmetric cruise missile flight for attacking static target [M.S. thesis], Southeast University, Dhaka, Bangladesh, 2016, (in Chinese).
 J. X. Zhang, Research and simulation of the cruise missile flight trajectory control [M.S. thesis], Southeast University, Dhaka, Bangladesh, 2015, (in Chinese).
 Y. Yang, Study on multi-dimensional Taylor network optimal control for ship roll stabilization [M.S. thesis], Southeast University, Dhaka, Bangladesh, 2016, (in Chinese).
 G. Bartolini and E. Punta, "Reduced-order observer in the sliding-mode control of nonlinear nonaffine systems," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 55, no. 10, pp. 2368-2373, 2010.
Qi-Ming Sun (1,2) and Hong-Sen Yan (1,2)
(1) School of Automation, Southeast University, Nanjing 210096, China
(2) Key Laboratory of Measurement and Control of Complex Systems of Engineering, Ministry of Education, Nanjing 210096, China
Correspondence should be addressed to Hong-Sen Yan; email@example.com
Received 15 June 2017; Accepted 14 September 2017; Published 30 October 2017
Academic Editor: Quanxin Zhu
Caption: Figure 1: The basic structure of MTN.
Caption: Figure 2: The control structure.
Caption: Figure 3: Output comparison.
Caption: Figure 4: Absolute error curve.
Caption: Figure 5: Output comparison.
Caption: Figure 6: Absolute error curve.
Caption: Figure 7: Output comparison.
Caption: Figure 8: Output comparison.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Sun, Qi-Ming; Yan, Hong-Sen|
|Publication:||Mathematical Problems in Engineering|
|Date:||Jan 1, 2017|
|Previous Article:||Analysis of China's Primary Energy Structure and Emissions Reduction Targets by 2030 Based on Multiobjective Programming.|
|Next Article:||Continuous Analog of Accelerated OS-EM Algorithm for Computed Tomography.|