Hybrid Recurrent Laguerre-Orthogonal-Polynomial NN Control System Applied in V-Belt Continuously Variable Transmission System Using Particle Swarm Optimization

Lin, Chih-Hong

doi:https://doi.org/10.1155/2015/106707

Mathematical Problems in Engineering

On this page

Abstract Introduction Experimental Results Conclusions References Copyright Related Articles

Special Issue

Recent Advances on Modeling, Control, and Optimization for Complex Engineering Systems

View this Special Issue

Research Article | Open Access

Volume 2015 | Article ID 106707 | https://doi.org/10.1155/2015/106707

Hybrid Recurrent Laguerre-Orthogonal-Polynomial NN Control System Applied in V-Belt Continuously Variable Transmission System Using Particle Swarm Optimization

Chih-Hong Lin¹

Academic Editor: Kang Li

Received18 May 2014

Accepted01 Sept 2014

Published26 Oct 2015

Abstract

Because the V-belt continuously variable transmission (CVT) system driven by permanent magnet synchronous motor (PMSM) has much unknown nonlinear and time-varying characteristics, the better control performance design for the linear control design is a time consuming procedure. In order to overcome difficulties for design of the linear controllers, the hybrid recurrent Laguerre-orthogonal-polynomial neural network (NN) control system which has online learning ability to respond to the system’s nonlinear and time-varying behaviors is proposed to control PMSM servo-driven V-belt CVT system under the occurrence of the lumped nonlinear load disturbances. The hybrid recurrent Laguerre-orthogonal-polynomial NN control system consists of an inspector control, a recurrent Laguerre-orthogonal-polynomial NN control with adaptive law, and a recouped control with estimated law. Moreover, the adaptive law of online parameters in the recurrent Laguerre-orthogonal-polynomial NN is derived using the Lyapunov stability theorem. Furthermore, the optimal learning rate of the parameters by means of modified particle swarm optimization (PSO) is proposed to achieve fast convergence. Finally, to show the effectiveness of the proposed control scheme, comparative studies are demonstrated by experimental results.

1. Introduction

A V-belt continuously variable transmission (CVT) [1–8] is typically composed of two hydraulically, or spring, actuated variable radii pulleys and a chain, or metal pushing, belt. To launch a vehicle from rest, the input pulley radius will be smaller than the output pulley radius, resulting in a speed reduction and torque multiplication transmitted to the drive shaft. For increased output shaft speed, the pulley radii are inversely manipulated simultaneously (i.e., input pulley radius increases as the output pulley radius decreases) to maintain constant belt length. A CVT may operate at a specific speed while changing the pulleys’ radii to achieve torque multiplication, acceleration, and speed as per the vehicle’s velocity, load requirements, engine power, and gear ratios. This operating profile provides the research motivation for CVT dynamics and nonlinear control algorithms. CVT-based vehicles have been traditionally regulated using a standard proportional integral derivative- (PID-) based controller with measurements of the gear ratio [3]. It has also been demonstrated that this control strategy provides satisfactory performance using gain-scheduling with a large set of points. In addition, numerous fuzzy logic controllers [4] have also been proposed. However, V-belt continuously variable transmission (CVT) system driven by alternating current (AC) motor is yet not shown in any commercial reports so that it provides the research motivation in this study.

The AC motor has several types such as permanent magnet synchronous motor (PMSM), switched reluctance motor (SRM), and induction motor (IM). In order to select the appropriate AC motor for driven V-belt CVT system, high efficiency is one of the most important factories to be selected. The PMSM provides higher efficiency, higher power density, and lower power loss for their size compared to SRM and IM. In addition, field-oriented control is one the most popular control techniques for the PMSM servo-driven system. As a result, torque ripple of the PMSM is lower than the SRM and IM. On the other hand, the PMSM controlled by field-oriented control, which can be achieved fast by four-quadrant operation, is much less sensitive to the parameters variation of the motor [9–11]. Therefore, the PMSM has been widely used in many industrial applications such as robotics, electric power steering, and other mechatronics [9–11].

Artificial neural networks (ANNs) have emerged as a powerful learning technique to perform complex tasks such as highly nonlinear approximations and the control of dynamical systems [12–16]. Some of the prime advantages of using NN are their ability to learn based on optimization of an appropriate error function and their excellent performance for approximation of nonlinear functions. There are different paradigms of NNs proposed by different researchers for the tasks of system identifications and controls [13–16]. One of the major drawbacks of the NN is that it is computationally intensive and needs large number of iterations for its training. In order to reduce the computational complexity, a functional-link NN, which has shown that it is capable of producing similar performance as that of NN but with much less computational cost, was reported in [17–19]. These functional-link NNs [17–19] with faster convergence and lesser computational complexity were executed in the identification and controls of nonlinear dynamic system with satisfactory results. Recently, the Laguerre-functional-expansions combined with NN, which was applied in highly nonlinear approximations and the control of dynamical systems, have been proposed [20–24]. Aadaleesan et al. [20] proposed the Laguerre filter combined with the wavelet network in order to approximate the memoryless nonlinearity. Approximation the linear and nonlinear parts of a Wiener structure by means of the Laguerre filter and the general feed-forward NN was reported in [21]. The Laguerre-functional-expansions feed-forward NN, which employed Laguerre-orthogonal-polynomials in the activation functions of the hidden neurons in order to identify models of the chaotic time series, was proposed by Zou and Xiao [22]. Patra et al. [23] proposed a computationally efficient Laguerre NN, which is based on Laguerre-functional-expansions to autocompensate for the associated nonlinearity and environmental dependence for intelligent sensors, and provide linearized sensor readout even when the motes are operated in harsh environments. Patra et al. [24] present an intelligent technique by means of novel computationally efficient Laguerre NN to compensate for the inherent sensor nonlinearity and the environmental influences. Since the Laguerre NN is a single-layer NN, its computational complexity is found to be much lower than a multilayer perception (MLP). However, these Laguerre-functional-expansions feed-forward NNs without a feedback loops can be used for static function approximation, but they cannot adequately approximate dynamic behaviors found in PMSM servo-driven V-belt CVT system with nonlinear and time-varying characteristics.

The recurrent NN has received increasing attention due to its structural advantage in the modelling of the nonlinear system and dynamic control of the system [25–29]. These networks are capable of effective identification and control of complex process dynamics, but with the expense of large computational complexity. Hence, if each neuron in the recurrent NN is considered as a state in the nonlinear dynamic systems, the self-connection feedback type is able to approximate the dynamic systems efficiently [25–29]. In order to improve the ability of identifying high order systems and reduce computational complexity, the recurrent Laguerre-orthogonal-polynomial NN, which has more advantages than the Laguerre-orthogonal-polynomial NN including better performance, higher accuracy, and dynamic robustness, has been proposed to control the PMSM servo-driven V-belt CVT system with nonlinear and time-varying characteristics in this paper.

Particular swarm optimization (PSO) is a population-based, self-adaptive search optimization technique first introduced by Kennedy and Eberhart [30]. Similar to genetic algorithms [31], an evolutionary algorithm approach, the PSO is an evolutionary optimization tool of swarm intelligence field based on a swarm (population), where each member is seen as a particle, and each particle is a potential solution to the problem under analysis. The motivation for the development of this method was based on the simulation of simplified animal social behaviors such as fish schooling, bird flocking, and so forth. However, unlike in other evolutionary optimization methods, in PSO there is no direct recombination of genetic material between individuals during the search. The PSO algorithm works on the social behavior of particles in the swarm. Therefore, it finds the global best solution by simply adjusting the trajectory of each individual toward its own best location and toward the best particle of the entire swarm at each time step (generation) [30, 32]. Clerc and Kennedy [32] introduced the concept of inertia weight to the original version of PSO, in order to balance the local and global search during the optimization process. Thus, PSO has been widely applied in mathematical modeling, dynamic programming, and system control [33–36] due to simple structure, simple parameter setting, and fast convergence speed. How to improve the convergence speed and how to guarantee the convergence of PSO are the main problems of PSO improvement [37] and are gradually turning into a hot topic in this field. In order to weigh the relationship between local search and global search, Clerc and Kennedy [32] and Eberhart and Shi [38, 39] proposed improved particle swarm optimization with inertial weight to control the exploitation and exploration. Meanwhile, some researchers [40–42] have proposed the topical improved particle swarm algorithm with inertia factor, which is called topical particle swarm optimization. However, the PSO existed in premature convergence problem and the modified PSO is proposed to prevent premature convergence and to acquire optimal learning rate with better convergence in this paper.

In this study the hybrid recurrent Laguerre-orthogonal-polynomial NN control system is developed to control the V-belt CVT system with many nonlinear dynamics [1–8, 43–46], which is driven by PMSM. The hybrid recurrent Laguerre-orthogonal-polynomial NN control system has fast learning property and good generalization capability. The control method, which is not dependent upon the predetermined characteristics of the motor, can adapt to any change in the motor characteristics. The hybrid recurrent Laguerre-orthogonal-polynomial NN control system, which is composed of the inspector control, the recurrent Laguerre-orthogonal-polynomial NN control with adaptive law, and the recouped control, is applied to the V-belt CVT system driven by PMSM. The adaptive law of the online parameter in the recurrent Laguerre-orthogonal-polynomial NN can be derived according to the Lyapunov stability theorem and the gradient descent method. The recurrent Laguerre-orthogonal-polynomial NN has the online learning ability to respond to the system’s nonlinear and time-varying behaviors under the occurrence of the lumped nonlinear external disturbances with parameters variation. Furthermore, two optimal learning rates of the parameters by means of modified PSO are proposed to achieve fast convergence. Finally, the control performances of the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system are verified by experimental results.

The paper is structured as follows: Section 2 provides the configuration of the V-belt CVT system driven by PMSM. Section 3 develops the proposed novel hybrid recurrent Laguerre-orthogonal-polynomial NN control system for controlling the V-belt CVT system driven by PMSM. Section 4 presents the experimental results for comparisons between the proposed control method and PI control method at three cases. Section 5 provides the conclusions.

2. Configuration of System

Since the electric scooter system has much unknown nonlinear uncertainties and parameter variations, such as load torque, rolling resistance, wind resistance, and braking force, the V-belt CVT and clutch in the scooter model can be categorized as functioning in one of two operating modes depending on the speed of the V-belt CVT output axis: disengaged or completely coupled. At the start of the PMSM drive cycles, the scooter is in an idle state. The clutch is initially disengaged, and subsequent transition between modes is controlled by the clutch axis rotational speed. Except for the mechanical losses, the PMSM power is transmitted through the V-belt CVT and clutch to the wheel in the electric scooter.

2.1. Structure of the V-Belt CVT System Driven by PMSM

The development of the V-belt CVT began with rubber V-belts [5]. Despite the fact that rubber V-belt CVTs are not well suited for automotive applications because of their limited torque capacity, there are some interesting concepts on the market. The V-belt CVT consists of a segmented rubber V-belt and two shafts with conical pulleys. The V-belt is clamped between two pairs of conical sheaves. In the V-belt CVT, the transmission ratio is determined by simultaneous adjustment of the running radii of the belt on the pulleys. On each shaft, there is one fixed and one axially moveable sheave. Axial movement of the moveable sheave adjusts the gap between the sheaves and thereby the belt running radius. The input shaft of the V-belt CVT is called the primary shaft which mounted the PMSM, and the output shaft is the secondary shaft which mounted the wheel. The structure of the V-belt CVT is shown in Figure 1. The appearance of the primary pulley side and the secondary pulley side in the V-belt CVT is shown in Figure 1(a), and the cross-section view of the secondary shaft in the V-belt CVT is shown in Figure 1(b). The wheel of electric scooter is connected to the output shaft of the final reduction using a torsion spring, which models the combined stiffness of both the drive shafts. The electric scooter inertia is connected to the wheel using linear damper, which models the tire force. The rolling resistance on the electric scooter is modeled as a load torque. In order to reduce system complexity, the torque dynamic equations in the primary drive shaft and the secondary drive shaft of the V-belt CVT shown in Figure 2 can be simplified as [1–8] in which [1–8] is the lumped nonlinear external disturbances of the secondary drive side on the wheel; is the drive torque of the primary pulley shaft; is the drive torque of the secondary pulley shaft; is the conversion ratio with respect to secondary pulley shaft transferred to primary pulley shaft of V-belt arc length; is the rolling resistance; is the wind resistance; is a braking force; is the total wind velocity; represents the total frictional coefficient of ground surface; and represent the viscous frictional coefficients of the PMSM and the wheel, respectively; and are the inertias of the PMSM and the wheel, respectively; and are the speeds of the PMSM and the wheel, respectively. Then using speed ratio and sliding ratio [1–8], the torque equation can be transformed from the secondary pulley side to the primary pulley side. Therefore, the resultant dynamic equation of the PMSM driven V-belt CVT system from (1) can be simplified as [1–8, 41–45]in which [1–8] is the resultant lumped nonlinear external disturbances with parameter variations; is the fixed load torque; is the resultant parameter variation; is the resultant rolling resistance; is the resultant wind resistance; is the resultant braking force; is the resultant unknown nonlinear load torque; is the resultant viscous frictional coefficient; is the resultant moment of inertia.

(a)

(b)

2.2. Structure of the PMSM Driven System

For convenient design, the stator voltage equations of the PMSM driven V-belt CVT system in the synchronously rotating reference frame can be described as follows [9–11, 42–45]: in which and are the -axis and -axis stator voltages; and are the -axis and -axis stator currents; and are the -axis and -axis stator inductances; is the -axis permanent magnet flux linkage; is the stator resistance; is the electrical angular speed.

The electromagnetic torque of the PMSM driven V-belt CVT system can be described asin which is the number of poles. Due to for a surface-mounted PMSM, the second term of (4) is zero. Moreover, is a constant for surface-mounted PMSM. The rotor flux is produced in the -axis only, while the current vector is generated in the -axis for the field-oriented control. When the -axis rotor flux is a constant and torque angle is [9–11], the maximum torque per ampere can be reached for the field-oriented control. The electromagnetic torque is linearly proportional to the -axis current , which is determined by closed-loop control. The control principle of the PMSM driven system is based on field orientation. The PMSM driven V-belt CVT system with the implementation of field-oriented control can be reduced as , in which is the torque constant. The block diagram of the V-belt CVT system driven by PMSM is shown in Figure 3.

The whole system of the PMSM driven V-belt CVT system can be indicated as follows: a field-oriented institution, a current PI control loop, a sinusoidal PWM control circuit, an interlock circuit and an isolated circuit, an IGBT power module inverter, and a speed control loop. The PI current loop controller is the current loop tracking controller. In order to attain good dynamic response, all gains for well-known PI current loop controller are listed as follows: and through some heuristic knowledge [47–50] on the tuning of the PI controller. The field-oriented institution consists of the coordinate transformation, generation, and lookup table generation. The TMS320C32 DSP control system manufactured by Spinel Technology Corporation is used to implement field-oriented control and speed control. The V-belt CVT system driven by PMSM is manipulated under the lumped external disturbances with nonlinear uncertainties.

3. Design of Hybrid Recurrent Laguerre-Orthogonal-Polynomial NN Control System

Due tononlinear uncertainties of the V-belt CVT system driven by PMSM such as nonlinear friction force of the transmission belt and clutch, rolling resistance, wind resistance, and braking force, these will lead to degeneration of tracking responses in command current and speed of the V-belt CVT system driven by PMSM. These nonlinear uncertainties cause variation of rotor inertia and friction for the V-belt CVT system driven by PMSM. For convenient hybrid recurrent Laguerre-orthogonal-polynomial NN control system design, the dynamic equation of the V-belt CVT system driven by PMSM from (2) can be rewritten asin which is the control effort, that is, the command torque of the PMSM. , , and are three known constants. When the uncertainties including variation of system parameters and external force disturbance occur, the parameters are assumed to be bounded, that is, , , and , where is a known continuous function; and are two known constants. Then, the tracking error can be defined as where represents the desired command rotor speed; is the tracking error between the desired rotor speed and actual rotor speed. If all parameters of the V-belt CVT system driven by PMSM including the lumped nonlinear external disturbances and parameter variation are well known, the ideal control law can be designed asin which is a positive constant. Substituting (6) into (4), the error dynamic equation can be obtained as follows:

The system state can track the desired trajectory gradually if as in (8). However, the hybrid recurrent Laguerre-orthogonal-polynomial NN control system is proposed to control the CVT driven by PMSM under uncertain perturbations. The configuration of the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system is described in Figure 4.

The hybrid recurrent Laguerre-orthogonal-polynomial NN control system is composed of an inspector control system, a recurrent Laguerre-orthogonal-polynomial NN controller, and a recouped controller. The control law is designed as where is the proposed inspected control that is able to stabilize around a predetermined bound area in the states of the controlled system; is the recurrent Laguerre-orthogonal-polynomial NN control which is as the major tracking controller. It is used to imitate an ideal control law. The recouped control is designed to recoup the difference between the ideal control law and the recurrent Laguerre-orthogonal-polynomial NN control. Since the inspected control caused the overdone and chattering effort, the recurrent Laguerre-orthogonal-polynomial NN control and the recouped control are proposed to reduce and smooth the control effort when the system states are inside the predetermined bound area. When the recurrent Laguerre-orthogonal-polynomial NN approximation properties cannot be ensured, the inspected control is able to act in this case.

For the condition of divergence of states, the design of hybrid recurrent Laguerre-orthogonal-polynomial NN control system is essential to stretch the divergent states back to the predestinated bound area. The hybrid recurrent Laguerre-orthogonal-polynomial NN control system can uniformly approximate the ideal control law inside the bound area. Then stability of the hybrid recurrent Laguerre-orthogonal-polynomial NN control system can be warranted. An error dynamic equation from (5) to (9) can be acquired as

Firstly, the inspected control can be designed asin which is a sign function. When the recurrent Laguerre-orthogonal-polynomial NN approximation properties cannot be ensured, the inspected control law is able to act in this case; that is, . Due to the inadequate bound values, for example, , , , and sign function, the inspected control can produce in overdone and chattering effort. Therefore, the recurrent Laguerre-orthogonal-polynomial NN control and the recouped control can be devised to conquer the mentioned blemish. The recurrent Laguerre-orthogonal-polynomial NN control raised to imitate the ideal control . Then the recouped control posed to recoup the difference between the ideal control and the recurrent Laguerre-orthogonal-polynomial NN control .

Secondly, the architecture of the proposed three-layer recurrent Laguerre-orthogonal-polynomial NN is depicted in Figure 5. It is composed of an input layer, a hidden layer, and an output layer. The activation functions and signal actions of nodes in each layer of the recurrent Laguerre-orthogonal-polynomial NN can be described as follows.

First Layer: Input Layer. Each node in this layer is indicated by using , which multiplies by each other between each other for input signals. Then outputs signals are the results of product. The input and the output for each node in this layer are expressed asThe is the tracking error between the desired speed and the rotor speed . The is the tracking error change. The is the recurrent weight between output layer and input layer. The denotes the number of iterations. The is the output value of the output layer in the recurrent Laguerre-orthogonal-polynomial NN.

Second Layer: Hidden Layer. The single node th in this layer is labeled with . The net input and the net output for node th of the hidden layer are expressed asLaguerre-orthogonal polynomials [20–24] are selected for activation function of the hidden layer. The Laguerre-orthogonal polynomials are denoted by , where is the order of expansion and is the argument of the polynomial, is the number of nodes, and is the self-connecting feedback gain of the hidden layer which is selected between 0 and 1. The zero, the first, and the second order Laguerre-orthogonal polynomials are given by , , and , respectively. The higher order polynomials are given by and . The higher order Laguerre-orthogonal polynomials may be generated by the recursive formula given by .

Third Layer: Output Layer. The single node th in this layer is labeled with . It computes the overall output as the summation of all input signals. The net input and the net output for node th in this layer are expressed aswhere is the connective weight between the hidden layer and the output layer; is the activation function, which is selected as a linear function; represents the th input to the node of output layer. The output value of the recurrent Laguerre-orthogonal-polynomial NN can be represented as . Then the output value of the recurrent Laguerre-orthogonal-polynomial NN, , can be denoted asin which is the adjustable weight parameters vector between the hidden layer and the output layer of the recurrent Laguerre-orthogonal-polynomial NN. is the inputs vector in the output layer of the recurrent Laguerre-orthogonal-polynomial NN, in which is determined by the selected Laguerre-orthogonal-polynomials.

Thirdly, in order to evolve the recouped control , a minimum approximation error is defined asin which is an ideal weight vector to reach minimum approximation error. It is assumed that absolute value of is less than a small positive value ; that is, . Then, the error dynamic equation from (10) can be rewritten asThen, the Lyapunov function is selected asin which is a learning rate, is an adaptation gain, and is the bound estimated error. Differentiating the Lyapunov function with respect to and using (16), then (17) can be rewritten asIn order to obtain , the adaptive law and the recouped controller with estimation law can be designed as follows:Substituting (11), (20), and (21) into (19) and using (11) with , then (19) can be represented asSubstituting (22) into (23), then (23) can be obtained asFrom (24), the is a negative semidefinite; that is, . It implies that and are bounded. In addition, the function is defined asIntegrating (25) with respect to , then Due to the fact that is bounded and is nonincreasing and bounded, thenDifferentiating (25) with respect to givesDue to the fact that all the variables in the right side of (17) are bounded. It implies that is also bounded. Then, is a uniformly continuous function [51, 52]. It is denoted that by using Barbalat’s lemma [51, 52]. Therefore, as . Furthermore, In order to avoid chattering phenomenon of recouped controller , the sign function can be replaced by the equation , where in which and are positive constants.

According to Lyapunov stability theorem and the gradient descent method, an online parameter training methodology of the recurrent Laguerre-orthogonal-polynomial NN can be derived and trained effectively. Then the parameter of adaptive law can be computed by the gradient descent method. The adaptive law shown in (20) calls for a proper choice of the learning rate. For a small value of learning rate, the convergence of controller parameter can be guaranteed but the convergent speed is very slow. In order to achieve the fast convergence of output tracking error, the modified PSO is used to obtain optimal learning rate for training parameters of the recurrent Laguerre-orthogonal-polynomial NN. First, the parameter of adaptive law shown in (20) can be rewritten asThe central part of the training algorithm for the recurrent Laguerre-orthogonal-polynomial NN is concerned with how to obtain recursively a gradient vector in which each element in the training algorithm is defined as the derivative of an energy function with respect to a parameter of the network. This is done by means of the chain rule, because the gradient vector is calculated in the direction opposite to the flow of the output of each node. In order to describe the online training algorithm of the recurrent Laguerre-orthogonal-polynomial NN, a cost function is defined asAccording to the gradient descent method, the adaptive law of the weight also can be represented asComparing (30) with (32) yields . Then, the convergence analysis in the following theorem is to derive specific learning rate to assure convergence of the output tracking error.

The adaptation law of recurrent weight using the gradient descent method can be updated asin which is a learning rate. Then, the convergence analysis in the following theorem is to derive specific learning rate to assure convergence of the output tracking error. In order to obtain better learning rate, therefore, the modified PSO is proposed to search for a better learning rate of the weights of the recurrent Laguerre-orthogonal-polynomial NN.

The PSO [30–42], which has three parameters as two acceleration coefficients , and inertia weight , has a significant impact on performance of the algorithm, especially the impact of inertia weight. The impact is different on different conditions and is also different at different times under the same condition. The PSO with larger inertia weight has a faster convergence speed and works well in global search, while the PSO with smaller inertia weight can reach a more accurate optimum value but only works well in local search. Now dynamic inertia weight modification is used to train the appropriate value of , in order to coordinate between search accuracy and search speed. The modified PSO algorithm depends in its implementation in the following two relations:where indicates the iterations; is the current position of particle in hyperspace and with respect to learning rate , ; is the current speed of particle ; is the inertia weight with respective to ; and are acceleration positive constants; , , and are random numbers obtained from a uniform random distribution function in the interval ; is the initial positive constant in the interval ; and are the initial positive constants in the interval ; denotes the number of generations; and represent the best previous position of the th particle and the position of the best particle among all particles in the population, respectively. Usually, is modified gradually in descending within the limit of 0.8 < < 1.4 [40], so that search space can be changed steadily from the global to the local. Decreasing inertia weight particle swarm optimization is a topical algorithm, of which inertia weight decreases linearly from 0.9 to 0.4 [38, 39]. Some scholars propose the increasing inertia weight particle swarm optimization of which inertia weight increases linearly from 0.4 to 0.9 [40, 41]. is the constriction factor introduced by Eberhart and Shi [38, 39] to avoid the swarm from premature convergence and to ensure stability of the system.

4. Experimental Results

The whole system of the DSP-based control system for the PMSM driven V-belt CVT system is shown in Figure 1. The control algorithm was executed by a TMS320C32 DSP control system including four channels of D/A, eight channels of programmable PWM, and one encoder interface circuit. The IGBT power module voltage source inverter is executed by current-controlled SPWM with a switching frequency of 15 kHz. The specification of PMSM is a three-phase 48 V, 750 W, 16.5 A, and 3600 rpm. The parameters of the PMSM are given as follows as , mH, Nms, Nms/rad, and Nm/A by using open circuit test, short test, rotor block test, and loading test. Owing to inherent uncertainty in V-belt CVT system (e.g., the lumped nonlinear external disturbances and parameter variations) and output current limitation of battery capacity, the PMSM can only operate at 376.8 rad/s (3600 rpm) to avoid burning IGBT modules for V-belt CVT system at high speed perturbation.

To show the control performance of the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system, two cases are provided in the experimentation. One is the 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation and the other is the 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations . In order to achieve good transient and steady-state control performance all gains of the well-known PI controller are and through some heuristic knowledge [47–50] on the tuning of the PI controller at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation for the speed tracking. The experimental results of the well-known PI controller for the V-belt CVT system driven by PMSM at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation and 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations are shown in Figures 6, 7, and 8, where tracking responses of the command rotor speed , the desired command rotor speed , and the measured rotor speed are shown in Figures 6(a) and 7(a); tracking responses of the speed error are shown in Figures 6(b) and 7(b); tracking responses of the speed error amplification are shown in Figures 6(c) and 7(c). Since the low speed operation is the same as the nominal case due to smaller disturbance, the response of speed shown in Figure 6(a) has better tracking performance. Moreover, the degenerate tracking responses of speed shown in Figure 7(a) is obvious due to bigger nonlinear disturbance (e.g., rolling resistance, wind resistance, and parameter variation) at high speed operation. From the experimental results, sluggish tracking responses of speed and current are obtained for the V-belt CVT system driven by PMSM using the well-known PI controller. The linear controller has the weak robustness under bigger nonlinear disturbance because of no appropriate gains tuning or no degenerate nonlinear effect. In addition, responses of the command electromagnetic torque at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation and 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations are shown in Figures 8(a) and 8(b), respectively. The dynamic response of command electromagnetic torque results in great torque ripple due to V-belt CVT system with nonlinear disturbance, such as V-belt shaking friction and action frictions between primary pulley and second pulley.

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

The control gains of the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system are , . All control gains of the hybrid recurrent Laguerre-orthogonal-polynomial NN control system are chosen to achieve the best transient control performance in experimentation considering the requirement of stability. The parameter adjustment process remains continually active for the duration of the experimentation. The structure of the recurrent Laguerre-orthogonal-polynomial NN controller has 2 nodes, 3 nodes, and 1 node in the input layer, the hidden layer, and the output layer, respectively. The experimental results of the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system for the V-belt CVT system driven by PMSM at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation and 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations are shown in Figures 9, 10, and 11, where tracking responses of the command rotor speed , the desired command rotor speed , and the measured rotor speed are shown in Figures 9(a) and 10(a); tracking responses of the speed error are shown in Figures 9(b) and 10(b); tracking responses of the speed error amplification are shown in Figures 9(c) and 10(c). Since the low speed operation is the same as the nominal case due to smaller disturbance, the speed tracking response shown in Figure 9(a) has better tracking performance. The better speed tracking response is shown in Figure 10(a) under the occurrence of bigger lumped nonlinear external disturbances and parameter variation. From the experimental results, the accurate tracking performance is obtained for the V-belt CVT system driven by PMSM using the hybrid recurrent Laguerre-orthogonal-polynomial NN control system owing to the online adaptive mechanism of the recurrent Laguerre-orthogonal-polynomial NN and action of the recouped controller. In addition, responses of the command electromagnetic torque at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation and 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations are shown in Figures 11(a) and 11(b), respectively. In addition, in order to decrease inertia weight particle swarm optimization, the initial coefficient of inertia weight is set as 0.4. Two initial coefficients and of the constriction factor are set as 0.3 to avoid the swarm from premature convergence and to ensure stability of the system. Thus, the convergence responses of two learning rates and of the recurrent Laguerre-orthogonal-polynomial NN using modified PSO at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation are shown in Figures 12(a) and 12(b), respectively. The convergence responses of two learning rates and of the recurrent Laguerre-orthogonal-polynomial NN using modified PSO at 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations are shown in Figures 13(a) and 13(b), respectively. The convergence response of two learning rates and of the recurrent Laguerre-orthogonal-polynomial NN using modified PSO have faster convergence speed than the fixed learning rate of the recurrent Laguerre-orthogonal-polynomial NN. The dynamic response of command electromagnetic torque results lower torque ripple by on-line adjusted of the recurrent Laguerre-orthogonal-polynomial NN to cope with the high-frequency unmodeled dynamic of the V-belt CVT system with nonlinear disturbance, such as V-belt shaking friction and action frictions between primary pulley and second pulley. Therefore, these results show that the hybrid recurrent Laguerre-orthogonal-polynomial NN control system has better control performance than the well-known PI controller under high speed perturbation for the V-belt CVT system driven by PMSM. Moreover, the numerical results of the iteration performance for two learning rates and by means of modified PSO method from 0 s to 0.6 s with 240 iterations (runs) at 188.4 rad/s (1800 rpm) case under the lumped nonlinear external disturbances with parameter variation and 376.8 rad/s (3600 rpm) case under the lumped nonlinear external disturbances with twice parameter variations are shown in Figures 14 and 15, respectively. In all cases, the numerical results of the iteration performance show that modified PSO can reach the global optimum effectively and avoid premature convergence in less than 100 iterations. From above, it can be easily observed that the proposed modified PSO can provide a more accurate optimal solution and converges to the criteria with a greater probability.

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

The measured rotor speed responses due to step disturbance torque are given finally. The condition under load torque disturbance with adding load is tested by using the PI controller and the hybrid recurrent Laguerre-orthogonal-polynomial NN control system. The experimental result of load adjustment using well-known PI controller under load torque disturbance with adding load at command rotor speed 376.8 rad/s (3600 rpm) is shown in Figure 16. The experimental results of the measured rotor speed response using well-known PI controller under load disturbance with adding load at 376.8 rad/s (3600 rpm) is shown in Figure 16(a). The experimental results of the measured current in phase using well-known PI controller under load disturbance with adding load at 376.8 rad/s (3600 rpm) is shown in Figure 16(b). The experimental result of load adjustment using the hybrid recurrent Laguerre-orthogonal-polynomial NN control system under load torque disturbance with adding load at command rotor speed 376.8 rad/s (3600 rpm) is shown in Figure 17. The experimental results of the measured rotor speed response using the hybrid recurrent Laguerre-orthogonal-polynomial NN control system under load disturbance with adding load at 376.8 rad/s (3600 rpm) is shown in Figure 17(a). The experimental results of the measured current in phase using the hybrid recurrent Laguerre-orthogonal-polynomial NN control system under load disturbance with adding load at 376.8 rad/s (3600 rpm) is shown in Figure 17(b). From the experimental results, the degenerated responses under load disturbance are much improved using the hybrid recurrent Laguerre-orthogonal-polynomial NN control system. From experimental results, transient response of the hybrid recurrent Laguerre-orthogonal-polynomial NN control system has faster convergence and better load regulation than the well-known PI controller at load regulation.

(a)

(b)

(a)

(b)

In addition, the control performance comparisons of the well-known PI controller and the hybrid recurrent Laguerre-orthogonal-polynomial NN control system are summarized in Table 1 for experimental results at three test cases. In Table 1, the hybrid recurrent Laguerre-orthogonal-polynomial NN control system results in smaller tracking errors with respect to the well-known PI controller. According to the tabulated measurements, the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system indeed yields the superior control performance compared to the well-known PI controller.

5. Conclusions

The hybrid recurrent Laguerre-orthogonal-polynomial NN control system has been successfully developed to control the V-belt CVT system driven by PMSM with robust control characteristics in this paper. The hybrid recurrent Laguerre-orthogonal-polynomial NN control system with inspected control based on the uncertainty bounds of the controlled system was designed to stabilize the system states around a predetermined bound area. To drop the excessive and chattering resulted by control efforts, the hybrid recurrent Laguerre-orthogonal-polynomial NN control system, which is composed of the inspected control, the recurrent Laguerre-orthogonal-polynomial NN control, and the recouped control, was proposed to reduce and smooth the control effort when the system states were inside the predetermined bound area.

The main contributions of this paper are as follows: the dynamic models for the V-belt CVT driven by PMSM with unknown nonlinear and time-varying characteristics were successfully derived; the hybrid recurrent Laguerre-orthogonal-polynomial NN control system for the V-belt CVT system driven by PMSM under the occurrence of the lumped nonlinear load disturbances was successfully applied to enhance robustness; the adaptive law of online parameters tuning in the recurrent Laguerre-orthogonal-polynomial NN and the estimation law of the recouped controller by using the Lyapunov stability theorem was successfully derived; two optimal learning rates of connective weights and recurrent weights in the recurrent Laguerre-orthogonal-polynomial NN by means of modified PSO are successfully applied to speed up convergence; the hybrid recurrent Laguerre-orthogonal-polynomial NN control system, which has better online learning capability in order to fast capture the system’s nonlinear and time-varying behaviors, is successfully developed; the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system has lower torque ripple than the well-known PI controller.

Finally, the control performance of the proposed hybrid recurrent Laguerre-orthogonal-polynomial NN control system is more suitable than the well-known PI controller for the V-belt CVT system driven by PMSM.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The author would like to acknowledge the financial support of the Ministry of Science and Technology in Taiwan, through its Grant MOST 103-2221-E-239-016.

References

C. Y. Tseng, L. W. Chen, Y. T. Lin, and J. Y. Li, “A hybrid dynamic simulation model for urban scooters with a mechanical-type CVT,” in Proceedings of the IEEE International Conference on Automation and Logistics (ICAL '08), pp. 515–519, Qingdao, China, September 2008.
View at: Publisher Site | Google Scholar
C. Y. Tseng, Y. F. Lue, Y. T. Lin, J. C. Siao, C. H. Tsai, and L. M. Fu, “Dynamic simulation model for hybrid electric scooters,” in Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '09), pp. 1464–1469, Seoul, Republic of Korea, July 2009.
View at: Publisher Site | Google Scholar
L. Guzzella and A. M. Schmid, “Feedback linearization of spark-ignition engines with continuously variable transmissions,” IEEE Transactions on Control Systems Technology, vol. 3, no. 1, pp. 54–58, 1995.
View at: Publisher Site | Google Scholar
W. Kim and G. Vachtsevanos, “Fuzzy logic ratio control for a CVT hydraulic module,” in Proceedings of the IEEE International Symposium on Intelligent Control, pp. 151–156, Rio, Greece.
View at: Google Scholar
G. Carbone, L. Mangialardi, B. Bonsen, C. Tursi, and P. A. Veenhuizen, “CVT dynamics: theory and experiments,” Mechanism and Machine Theory, vol. 42, no. 4, pp. 409–428, 2007.
View at: Publisher Site | Google Scholar
H. Sattler, “Efficiency of metal chain and V-belt CVT,” in Proceedings of the International Conference on Continuously Variable Power Transmissions, pp. 99–104, Eindhoven, The Netherlands, 1999.
View at: Google Scholar
G. Carbone, L. Mangialardi, and G. Mantriota, “The influence of pulley deformations on the shifting mechanism of metal belt CVT,” Journal of Mechanical Design, vol. 127, no. 1, pp. 103–113, 2005.
View at: Publisher Site | Google Scholar
N. Srivastava and I. Haque, “A review on belt and chain continuously variable transmissions (CVT): dynamics and control,” Mechanism and Machine Theory, vol. 44, no. 1, pp. 19–41, 2009.
View at: Publisher Site | Google Scholar
D. W. Novotny and T. A. Lipo, Vector Control and Dynamics of AC Drives, Oxford University Press, New York, NY, USA, 1996.
R. Krishnan, Electric Motor Drives: Modeling, Analysis, and Control, Prentice Hall, Upper Saddle River, NJ, USA, 2001.
F. J. Lin, “Real-time IP position controller design with torque feedforward control for pm synchronous motor,” IEEE Transactions on Industrial Electronics, vol. 44, no. 3, pp. 398–407, 1997.
View at: Publisher Site | Google Scholar
K. S. Narendra and K. Parthasarathy, “Identification and control of dynamical systems using neural networks,” IEEE Transactions on Neural Networks, vol. 1, no. 1, pp. 4–27, 1990.
View at: Publisher Site | Google Scholar
P. S. Sastry, G. Santharam, and K. P. Unnikrishnan, “Memory neuron networks for identification and control of dynamical systems,” IEEE Transactions on Neural Networks, vol. 5, no. 2, pp. 306–319, 1994.
View at: Publisher Site | Google Scholar
R. Griñó, G. Cembrano, and C. Torras, “Nonlinear system identification using additive dynamic neural networks—two on-line approaches,” IEEE Transactions on Circuits and Systems I, vol. 47, no. 2, pp. 150–165, 2000.
View at: Publisher Site | Google Scholar
Y. Li, “Stability of almost periodic solution for a general class of discontinuous neural networks with mixed time-varying delays,” Mathematical Problems in Engineering, vol. 2013, Article ID 236189, 12 pages, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
W. Dong, X. Liu, and Y. Li, “Analysis of stiffened penstock external pressure stability based on immune algorithm and neural network,” Mathematical Problems in Engineering, vol. 2014, Article ID 823653, 11 pages, 2014.
View at: Publisher Site | Google Scholar
Y. Pao, Adaptive Pattern Recognition and Neural Networks, Addison-Wesley, Boston, Mass, USA, 1989.
View at: MathSciNet
Y. H. Pao and S. M. Phillips, “The functional link net and learning optimal control,” Neurocomputing, vol. 9, no. 2, pp. 149–164, 1995.
View at: Publisher Site | Google Scholar | Zentralblatt MATH
J. C. Patra, R. N. Pal, B. N. Chatterji, and G. Panda, “Identification of nonlinear dynamic systems using functional link artificial neural networks,” IEEE Transactions on Systems, Man, and Cybernetics B: Cybernetics, vol. 29, no. 2, pp. 254–262, 1999.
View at: Publisher Site | Google Scholar
P. Aadaleesan, N. Miglan, R. Sharma, and P. Saha, “Nonlinear system identification using Wiener type Laguerre-Wavelet network model,” Chemical Engineering Science, vol. 63, no. 15, pp. 3932–3941, 2008.
View at: Publisher Site | Google Scholar
S. Mahmoodi, J. Poshtan, M. R. Jahed-Motlagh, and A. Montazeri, “Nonlinear model predictive control of a pH neutralization process based on Wiener-Laguerre model,” Chemical Engineering Journal, vol. 146, no. 3, pp. 328–337, 2009.
View at: Publisher Site | Google Scholar
A. Zou and X. Xiao, “An asynchronous encryption arithmetic based on Laguerre chaotic neural networks,” in Proceedings of the WRI Global Congress on Intelligent Systems (GCIS '09), pp. 36–39, Xiamen, China, May 2009.
View at: Publisher Site | Google Scholar
J. C. Patra, C. Bornand, and P. K. Meher, “Laguerre neural network-based smart sensors for wireless sensor networks,” in Proceeding of the IEEE Intrumentation and Measurement Technology Conference (I2MTC '09), pp. 832–837, Singapore, May 2009.
View at: Publisher Site | Google Scholar
J. C. Patra, P. K. Meher, and G. Chakraborty, “Development of laguerre neural-network-based intelligent sensors for wireless sensor networks,” IEEE Transactions on Instrumentation and Measurement, vol. 60, no. 3, pp. 725–734, 2011.
View at: Publisher Site | Google Scholar
T. W. S. Chow and Y. Fang, “A recurrent neural-network-based real-time learning control strategy applying to nonlinear systems with unknown dynamics,” IEEE Transactions on Industrial Electronics, vol. 45, no. 1, pp. 151–161, 1998.
View at: Publisher Site | Google Scholar
M. A. Brdyś and G. J. Kulawski, “Dynamic neural controllers for induction motor,” IEEE Transactions on Neural Networks, vol. 10, no. 2, pp. 340–355, 1999.
View at: Publisher Site | Google Scholar
X. D. Li, J. K. L. Ho, and T. W. S. Chow, “Approximation of dynamical time-variant systems by continuous-time recurrent neural networks,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 52, no. 10, pp. 656–660, 2005.
View at: Publisher Site | Google Scholar
C. H. Lu and C. C. Tsai, “Adaptive predictive control with recurrent neural network for industrial processes: an application to temperature control of a variable-frequency oil-cooling machine,” IEEE Transactions on Industrial Electronics, vol. 55, no. 3, pp. 1366–1375, 2008.
View at: Publisher Site | Google Scholar
C. F. Hsu, “Intelligent position tracking control for LCM drive using stable online self-constructing recurrent neural network controller with bound architecture,” Control Engineering Practice, vol. 17, no. 6, pp. 714–722, 2009.
View at: Publisher Site | Google Scholar
J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proceedings of the IEEE International Conference on Neural Networks, vol. 4, pp. 1942–1948, Perth, Australia, December 1995.
View at: Google Scholar
D. E. Goldberg, The Design of Innovation: Lessons from and for Competent Genetic Algorithms, Norwell, Kluwer Academic Publishers, Boston, Mass, USA, 2002.
View at: Publisher Site | MathSciNet
M. Clerc and J. Kennedy, “The particle swarm-explosion, stability, and convergence in a multidimensional complex space,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 1, pp. 58–73, 2002.
View at: Publisher Site | Google Scholar
A. B. de Carvalho, A. Pozo, and S. R. Vergilio, “A symbolic fault-prediction model based on multiobjective particle swarm optimization,” Journal of Systems and Software, vol. 83, no. 5, pp. 868–882, 2010.
View at: Publisher Site | Google Scholar
Q. Li, W. Chen, Y. Wang, S. Liu, and J. Jia, “Parameter identification for PEM fuel-cell mechanism model based on effective informed adaptive particle swarm optimization,” IEEE Transactions on Industrial Electronics, vol. 58, no. 6, pp. 2410–2419, 2011.
View at: Publisher Site | Google Scholar
C. T. Cheng, S. L. Liao, Z. T. Tang, and M. Y. Zhao, “Comparison of particle swarm optimization and dynamic programming for large scale hydro unit load dispatch,” Energy Conversion and Management, vol. 50, no. 12, pp. 3007–3014, 2009.
View at: Publisher Site | Google Scholar
Y. X. Liao, J. H. She, and M. Wu, “Integrated hybrid-PSO and fuzzy-NN decoupling control for temperature of reheating furnace,” IEEE Transactions on Industrial Electronics, vol. 56, no. 7, pp. 2704–2714, 2009.
View at: Publisher Site | Google Scholar
E. Zitzler, K. Deb, and L. Thiele, “Comparison of multiobjective evolutionary algorithms: empirical results,” Evolutionary Computation, vol. 8, no. 2, pp. 173–195, 2000.
View at: Publisher Site | Google Scholar
R. C. Eberhart and Y. Shi, “Comparison between genetic algorithms and particle swarm optimization,” in Proceedings of the 7th International Conference Evolutionary Programming, pp. 611–616, San Diego, Calif, USA, 1998.
View at: Google Scholar
R. C. Eberhart and Y. Shi, “Comparing inertia weights and constriction factors in particle swarm optimization,” in Proceedings of the Congress on Evolutionary Computation (CEC '00), pp. 84–88, La Jolla, Calif, USA, July 2000.
View at: Google Scholar
H. Gao and W. Xu, “A new particle swarm algorithm and its globally convergent modifications,” IEEE Transactions on Systems, Man, and Cybernetics B: Cybernetics, vol. 41, no. 5, pp. 1334–1351, 2011.
View at: Publisher Site | Google Scholar
T. Y. Sun, C. C. Liu, S. J. Tsai, S. T. Hsieh, and K. Y. Li, “Cluster guide particle swarm optimization (CGPSO) for underdetermined blind source separation with advanced conditions,” IEEE Transactions on Evolutionary Computation, vol. 15, no. 6, pp. 798–811, 2011.
View at: Publisher Site | Google Scholar
Y. Zhang, X. Xiong, and Q. Zhang, “An improved self-adaptive PSO algorithm with detection function for multimodal function optimization problems,” Mathematical Problems in Engineering, vol. 2013, Article ID 716952, 8 pages, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
C. H. Lin, P. H. Chiang, C. S. Tseng, Y. L. Liu, and M. Y. Lee, “Hybrid recurrent fuzzy neural network control for permanent magnet synchronous motor applied in electric scooter,” in Proceeding of the International Power Electronics Conference, pp. 1371–1376, Sapporo, Japan, June 2010.
View at: Publisher Site | Google Scholar
C. H. Lin and C. P. Lin, “The hybrid RFNN control for a PMSM drive electric scooter using rotor flux estimator,” International Journal of Electrical Power and Energy Systems, vol. 51, pp. 213–223, 2013.
View at: Publisher Site | Google Scholar
C. H. Lin, “Hybrid recurrent wavelet neural network control of PMSM servo-drive system for electric scooter,” International Journal of Control, Automation and Systems, vol. 12, no. 1, pp. 177–187, 2014.
View at: Publisher Site | Google Scholar
C. H. Lin, “Novel adaptive recurrent Legendre neural network control for PMSM servo-drive electric scooter,” ASME-Journal of Dynamic Systems, Measurement and Control, vol. 137, no. 1, 12 pages, 2015.
View at: Publisher Site | Google Scholar
J. G. Ziegler and N. B. Nichols, “Optimum settings for automatic controllers,” Transactions of the ASME, vol. 64, pp. 759–768, 1942.
View at: Google Scholar
K. J. Astrom and T. Hagglund, PID Controller: Theory, Design, and Tuning, Instrument Society of America, Research Triangle Park, Raleigh, NC, USA, 1995.
T. Hagglund and K. J. Astrom, “Revisiting the Ziegler-Nichols tuning rules for PI control,” Asian Journal of Control, vol. 4, no. 4, pp. 364–380, 2002.
View at: Google Scholar
T. Hagglund and K. J. Astrom, “Revisiting the Ziegler-Nichols tuning rules for PI control. Part II. The frequency response method,” Asian Journal of Control, vol. 6, no. 4, pp. 469–482, 2004.
View at: Google Scholar
J. J. E. Slotine and W. Li, Applied Nonlinear Control, Prentice Hall, Englewood Cliffs, NJ, USA, 1991.
K. J. Astrom and B. Wittenmark, Adaptive Control, Addison-Wesley, New York, NY, USA, 1995.

Copyright

Copyright © 2015 Chih-Hong Lin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1343

Downloads

1634

Citations