<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ICA</journal-id><journal-title-group><journal-title>Intelligent Control and Automation</journal-title></journal-title-group><issn pub-type="epub">2153-0653</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ica.2014.53012</article-id><article-id pub-id-type="publisher-id">ICA-48455</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>COMPUTER SCIENCE &amp; COMMUNICATIONS</subject><subject>ENGINEERING</subject><subject>PHYSICS &amp; MATHEMATICS</subject></subj-group></article-categories><title-group><article-title>A Study on Autotuning Controller for Servo System</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nguyen</surname><given-names>Hoang Giap</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jin-Ho</surname><given-names>Shin</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Won-Ho</surname><given-names>Kim</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Intelligent System Engineering, Graduate School of Dong-eui University, Busan, Korea</addr-line></aff><aff id="aff2"><addr-line>Department of Mechatronics Engineering, Dong-eui University, Busan, Korea</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>hoanggiap1005@gmail.com(NHG)</email>;<email>jhshin7@deu.ac.kr(JS)</email>;<email>kwh@deu.ac.kr(WK)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>03</day><month>07</month><year>2014</year></pub-date><volume>05</volume><issue>03</issue><fpage>102</fpage><lpage>110</lpage><history><date date-type="received"><day>16</day>	<month>May</month>	<year>2014</year></date><date date-type="rev-recd"><day>20</day>	<month>June</month>	<year>2014</year>	</date><date date-type="accepted"><day>3</day>	<month>July</month>	<year>2014</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
	This paper introduces a PID Autotuning controller using intelligent
neural network control based on relay feedback approach. The proposed
controller takes advantage of offline learning, in which the initial knowledge
of control system is recognized by the relay feedback approach, and the online
learning capability of neural network controller helps the control system
respond quickly to the dynamics changes. The robustness and motion tracking
performance are validated through simulation results.
</p></abstract><kwd-group><kwd>PID</kwd><kwd> Neural Network</kwd><kwd> Relay Feedback</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Several control design approaches have been developed as independent sub-fields of motion control domain, such as: adaptive control, robust control, model predictive control, lead-lag compensation. However, controller of PID type is the most commonly used in motion control [<xref ref-type="bibr" rid="scirp.48455-ref1">1</xref>] . Numbers of tuning methods for PID controller have been introduced since Ziegler and Nichols published their PID tuning rule in 1942. And until now it is still one of the most popular schemes. In this method, the control rules are defined by the ultimate gain and period, which are obtained by observing the oscillation of proportional closed-loop controller. This proposed scheme is difficult to carry out in motion control application since the amplitude of the oscillation is hard to keep under control.</p><p>A method for automatic determination of ultimate gain and period called relay feedback is used to overcome such drawbacks [<xref ref-type="bibr" rid="scirp.48455-ref2">2</xref>] -[<xref ref-type="bibr" rid="scirp.48455-ref5">5</xref>] . The ultimate gain and period of the system are determined by observing the ultimate frequency where the phase angle is -PI under a closed-loop relay feedback control. Although it is straightforward to find the PID control parameters from the calculated ultimate gain and period based on desired amplitude margin, the result is not optimum.</p><p>Conventional PID controller is well-suited for the linear control system which behaves in a predictable way. However, the parameters of PID controller are difficult to adapt to a wide range of uncertainty. Therefore, the robustness and performance of PID deteriorate in some industrial control applications, where the systems are normally subjected to rapid and unpredictable ways. Consequently, it is mandatory to incorporate the intelligent control function to the PID controller to adapt with the dynamics change of the control system. There are several researches and commercial products which apply such advance functions [<xref ref-type="bibr" rid="scirp.48455-ref6">6</xref>] -[<xref ref-type="bibr" rid="scirp.48455-ref12">12</xref>] .</p><p>In [<xref ref-type="bibr" rid="scirp.48455-ref13">13</xref>] , a phase locked loop identifier is designed to solve different open-loop frequency domain non-para- metric identification problems by giving exact estimation of the gain crossover point. In [<xref ref-type="bibr" rid="scirp.48455-ref14">14</xref>] , the author presents a new strategy of PID autotuning based on the combination of relay feedback methodology and operational optimization. For the operational optimization, data base of the control system is generated to develop the radial basis function network model. This model is applied to a nonlinear program to obtain the optimized control output. The relay feedback is applied to obtain the ultimate gain and frequency of control system. Then, the operating regions are separated, and the process data generated in the previous step are used to identify the frequency response model for each operating region.</p><p>In this research, we introduce an autotuning of PID controller based on radial basis function neural network (RBFNN) and relay feedback approach. The controller is designed as follows: firstly, the relay feedback control is carried out to analyze the system dynamic and calculate the ultimate gain and ultimate frequency of the control system by measuring the output amplitude and period; then a RBFNN with gradient descent algorithm is designed to approximate the system dynamic. Finally, the ultimate gain and ultimate frequency are incorporated in the neural network to adjust the coefficients related to the PID control parameters.</p></sec><sec id="s2"><title>2. Design of PID Autotuning Controller</title><sec id="s2_1"><title>2.1. Architecture of PID Autotuning Controller</title><p>The proposed PID AutoTuning controller using radial basis neural network based on relay feedback approach (ATNNRF) control structure is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. At the first phase of control, the relay feedback control mode is activated. The ultimate frequency and ultimate gain are obtained by the peak detection and calculation of the times between zero crossings of the control input level.</p><p>At the second phase of control, a neural network with back propagation algorithm is used to adjust the PID gains continuously.</p></sec><sec id="s2_2"><title>2.2. Relay Feedback Control</title><p>The ultimate gain and period of the system are determined by observing the ultimate frequency where the phase angle is -PI under a closed-loop relay feedback control.</p><p>Relay feedback control signal are shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><fig id="fig1"><label>Figure 1</label><caption><p> ATNNRF control system</p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\07325525-eeb0-46b6-b96b-b50ebf681439.png"/></fig><fig id="fig2"><label>Figure 2</label><caption><p> Relay feedback control signal</p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\125241b4-52e4-4db0-a769-54499ef1664d.png"/></fig><p>The relay feedback control is carried out as following:</p><disp-formula id="scirp.48455-formula2031"><label>(1)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ee28f4f1-2f8e-468a-a80b-4a9e4feffd63.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ac9bd673-f318-41b9-99f8-369f002ec6d7.png" xlink:type="simple"/></inline-formula> is the control input of the system at current time; <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\9239d752-85f6-490f-ac5b-4cd0f931f38e.png" xlink:type="simple"/></inline-formula>is the hysterical level at which the control input switch to the opposite position; <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\6d310c53-1e32-499d-824c-7ac19f8964b3.png" xlink:type="simple"/></inline-formula>is the position error of the control system at current time.</p><p>Applying relay feedback control (1) to the motion system, we obtain a periodic signal with amplitude <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\7dc84612-8f79-453d-801f-eaf6a71d6e35.png" xlink:type="simple"/></inline-formula> and period<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\e03bc6ea-e3e8-4c7a-a5c8-8c9c1ace0c99.png" xlink:type="simple"/></inline-formula>. In order to identify the ultimate gain from the relay feedback control, the describing function analysis is used.</p><p>Because the amplitude <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\34fe91a4-0ac9-43a9-83fd-925d860a054f.png" xlink:type="simple"/></inline-formula> and period <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\9bc120a7-4f4e-4af7-8cce-3f37e6d8ce28.png" xlink:type="simple"/></inline-formula> of the oscillation are measured from the relay feedback control, the ultimate gain <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\39ba8749-a12e-41df-9985-329ca0dad146.png" xlink:type="simple"/></inline-formula> and ultimate frequency <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\d0b94372-b0bd-44c0-8a76-81c0b597fdf0.png" xlink:type="simple"/></inline-formula> can be obtained as following:</p><disp-formula id="scirp.48455-formula2032"><label>(2)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\9843c751-d1a8-4053-a390-b6b8e8634437.png"/></disp-formula><disp-formula id="scirp.48455-formula2033"><label>(3)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\6caaebe1-7c28-464e-baa0-c2345645b3b6.png"/></disp-formula></sec><sec id="s2_3"><title>2.3. Radial Basis Function Neural Network</title><p>Since neural networks can approximate any continuous function over a compact set with high accuracy and fast learning capability, they have been recognized as a powerful tool for the applications of control system. It is proven that the standard 2-layer neural network consists of two layers of weights, threshold, a hidden layer and an output layer has sufficient generality for closed-loop control purpose.</p><p>In particular, the radial basis function neural network (RBFNN) is a feed-forward 2-layer neural network and demonstrates good control performance in the presence of unmodeled dynamics. The RBFNN can be considered as a local approximation model of control system. Moreover, the capability of online learning and rapid convergence of RBFNN makes it feasible in the application of control system.</p><p>The general RBFNN shown in <xref ref-type="fig" rid="fig3">Figure 3</xref> can be presented in matrix form:</p><disp-formula id="scirp.48455-formula2034"><label>(4)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ab35b1cc-347b-463a-b7f3-bcddc9965048.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\f387168e-8e03-43ed-a8bb-3fa4de4dbd07.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\482c1b75-28f0-4a34-8c83-b9f77618ed2d.png" xlink:type="simple"/></inline-formula> are the weight vectors, <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ba34e742-e1d3-4bc0-b51e-9f7b0580ce84.png" xlink:type="simple"/></inline-formula>is the input vector,</p><p><inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\569d6ccd-4c6f-45a7-a625-5c95aead01c6.png" xlink:type="simple"/></inline-formula>is the output vector, and <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\02d5574d-4dbd-4540-9cf7-a7e59682e75e.png" xlink:type="simple"/></inline-formula> denotes the activation function vector.</p><p>We know that the main disadvantage of the multi-layer neural network is highly nonlinear in parameter. Hence, we consider hereafter the neural network of fixed<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\0cd10ea8-9657-45ba-adb8-f60729443dff.png" xlink:type="simple"/></inline-formula>, which makes the neural network linearly parameterized. Defining<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\c0a6c676-a995-40d8-bd9e-d9ffe6b9efcc.png" xlink:type="simple"/></inline-formula>, we have:</p><fig id="fig3"><label>Figure 3</label><caption><p> Radial basis function neural network</p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\d74cfc8c-f65f-4c7c-9c95-a3c8b9bc5d08.png"/></fig><disp-formula id="scirp.48455-formula2035"><label>(5)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\b81f5bf3-a526-436a-85f9-87c254193a21.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\d162e999-ed2d-46ec-847e-49b6c4f4938c.png" xlink:type="simple"/></inline-formula> is the generalized basis function vector. In this paper, we use the radial basis function with Gaussian form as following:</p><disp-formula id="scirp.48455-formula2036"><label>(6)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\edb60d84-b323-464d-9611-299939cdf63c.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\8c57b472-ba7c-45ac-ade9-17d1705bcf01.png" xlink:type="simple"/></inline-formula> is the node center of the i-th radial basis function, and <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\6cb85471-bdc6-4bf1-bdb5-199fd428570f.png" xlink:type="simple"/></inline-formula> is the width coefficient.</p></sec><sec id="s2_4"><title>2.4. The Incorporation of Relay Feedback to the Autotuning Controller Using Neural Network</title><p>The proposed PID autotuning controller based on neural network and relay feedback (ATNNRF) is designed by the following procedures:</p><p>-Firstly, the relay feedback control is carried out to analyze the system dynamic and calculate the ultimate gain and ultimate frequency of the control system by measuring the output amplitude and period.</p><p>-Secondly, a RBFNN with gradient descent algorithm is designed to approximate the system dynamic. In this step, the Jacobian matrix, denotes the sensitivity of control system output to the control input is obtained.</p><p>-Lastly, the ultimate gain and ultimate frequency calculated in the first step are incorporated with the Jacobian matrix to adjust the coefficients related to the PID control parameters.</p><p>Consequently, the PID gains are tuned automatically by the self-learning capability of neural network. It is proved from the simulations results that the proposed controller has the strong robustness, high adaptibility and high performance compared with other current researches and commercial products.</p><p>The design of proposed ATNNRF is described as below.</p><p>The conventional PID controller can be represented with the discrete form as followings:</p><disp-formula id="scirp.48455-formula2037"><label>(7)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\fb982036-184b-426a-a7f2-dccd87d885ce.png"/></disp-formula><p>It is assumed that the closed-loop characteristics are invariant to the time constant.</p><p>Since the PID tuning parameters are uniform functions of ultimate gain and ultimate period as described in (2) and (3), they can be scaled for control system that has different time constant as followings:</p><disp-formula id="scirp.48455-formula2038"><label>(8)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\59d2bd3f-48e4-485c-b644-c310a27007d9.png"/></disp-formula><disp-formula id="scirp.48455-formula2039"><label>(9)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\b2eee344-3034-4c20-a1ac-80fd0b9f94a8.png"/></disp-formula><disp-formula id="scirp.48455-formula2040"><label>(10)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\287fd00d-4840-412e-a786-a4a827b0481c.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ef2a8cea-98c4-4256-b130-aa96f9ae5c2a.png" xlink:type="simple"/></inline-formula> are the scale factors that define the tuning rule.</p><p>The output of PID controller in (7) can be rewritten as the following discrete form:</p><disp-formula id="scirp.48455-formula2041"><label>(11)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\3072a752-7cab-4fe0-9f13-f16cc547b82e.png"/></disp-formula><disp-formula id="scirp.48455-formula2042"><label>(12)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\5231f3b4-c1fc-4684-ade9-498b9018d951.png"/></disp-formula><disp-formula id="scirp.48455-formula2043"><label>(13)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\dd1e61e4-f80c-4947-8e7b-b04046b0c651.png"/></disp-formula><disp-formula id="scirp.48455-formula2044"><label>(14)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\d3738955-fda2-4aec-a6f0-55eb0cc8fc98.png"/></disp-formula><disp-formula id="scirp.48455-formula2045"><label>(15)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\513cb6ec-2404-40cc-b972-62625e365980.png"/></disp-formula><disp-formula id="scirp.48455-formula2046"><label>(16)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\1906106b-763a-4fef-a3fb-9359a1de45b2.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ab017699-24d0-4d71-97fb-e318423aa0e6.png" xlink:type="simple"/></inline-formula> is the error of closed-loop control system.</p><p>Applying the tuning rule from (8), (9), and (10), we obtain the update control output as follow:</p><disp-formula id="scirp.48455-formula2047"><label>(17)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\aa962bb0-441f-4e87-82f7-93a26d3aaed7.png"/></disp-formula><p>Consider a RBFNN as in <xref ref-type="fig" rid="fig3">Figure 3</xref>. The input vector is presented as:</p><disp-formula id="scirp.48455-formula2048"><label>(18)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\edd7b1e2-9207-4314-98bd-358fc562ce7a.png"/></disp-formula><p>And the output of RBFNN is given by:</p><disp-formula id="scirp.48455-formula2049"><label>(19)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\8ddd0269-9251-4931-99c2-aaced3440eec.png"/></disp-formula><p>Define the performance cost function of RBFNN as:</p><disp-formula id="scirp.48455-formula2050"><label>(20)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\56d8642a-8647-4057-bd7c-0eb2bc925760.png"/></disp-formula><disp-formula id="scirp.48455-formula2051"><label>(21)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\0c4ba2c2-7196-41c6-81ef-eb8c309c10a5.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\e5da059f-ee1b-4814-bd54-1b2892548ad3.png" xlink:type="simple"/></inline-formula> is denoted as the approximation error.</p><p>Apply the gradient descent method, we can obtain the updating parameters of weight vector, node centers and width coefficients of RBFNN as followings:</p><disp-formula id="scirp.48455-formula2052"><label>(22)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\2d296e16-0265-4ada-b44f-1ff17497c3b0.png"/></disp-formula><disp-formula id="scirp.48455-formula2053"><label>(23)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\42bd2946-0db6-4c4b-912a-37ae667e91dc.png"/></disp-formula><disp-formula id="scirp.48455-formula2054"><label>(24)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\3a93dda5-dc7b-40b4-b509-721ec2be11f8.png"/></disp-formula><disp-formula id="scirp.48455-formula2055"><label>(25)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\1ba7c997-b0ca-4cf1-b54e-39d55d0b29e5.png"/></disp-formula><disp-formula id="scirp.48455-formula2056"><label>(26)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\9d8d7cb2-09a9-4a8c-9e4d-b3d9b4f7c7b6.png"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\898e149c-3e6f-48a5-ba46-e40ceb15be35.png" xlink:type="simple"/></inline-formula> is the learning rate of RBFNN, and <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\7b55d250-11d5-426b-97ad-6f3b9f2513a9.png" xlink:type="simple"/></inline-formula> is the momentum gene.</p><p>The Jacobian matrix for the model identification is obtained as follows:</p><disp-formula id="scirp.48455-formula2057"><label>(27)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\3621dc42-81cf-43be-86c9-573942720bcb.png"/></disp-formula><p>Gradient descent is an iterative method that is given an initial point, and follows the negative of the gradient in order to move the point toward a desired local minimum.</p><p>In this paper, we assume that the approximation error in (21) converges to zero by gradient descent algorithm, then (27) can be rewrite as follows:</p><disp-formula id="scirp.48455-formula2058"><label>(28)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\14b5e549-f5d8-4216-9e48-8621187a5eda.png"/></disp-formula><p>The PID control parameters can be obtained through the calculation of scale factor as described as followings:</p><p>Define the cost function of system as following:</p><disp-formula id="scirp.48455-formula2059"><label>(29)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\22cad965-062d-4fd3-ac22-516d372a7292.png"/></disp-formula><p>According to the gradient descent method, the scale factors can be updated by the following rules:</p><disp-formula id="scirp.48455-formula2060"><label>(30)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\19107aa1-7337-4d0b-b0a9-3267dd70bdef.png"/></disp-formula><disp-formula id="scirp.48455-formula2061"><label>(31)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\68959f44-66c8-4db4-9717-51c8e20e4c05.png"/></disp-formula><disp-formula id="scirp.48455-formula2062"><label>(32)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\e5e283d7-2584-42f9-bc3b-ff707743074d.png"/></disp-formula><p>From (17) and (27), we can obtain:</p><disp-formula id="scirp.48455-formula2063"><label>(33)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\aeea429c-aa21-4cae-9ca9-9157c580055f.png"/></disp-formula><disp-formula id="scirp.48455-formula2064"><label>(34)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\08d65573-5327-41b3-81b8-cd0ed9008ddb.png"/></disp-formula><disp-formula id="scirp.48455-formula2065"><label>(35)</label><inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\469a1b87-2218-4ce2-a751-d2dbf4a6ecab.png"/></disp-formula><p>Consequently, the PID control parameters can be obtained and the system is automatically tuned.</p></sec></sec><sec id="s3"><title>3. Simulation Results</title><p>The simulation is carried out by Matlab/Simulink with required inputs to achieve desired level of control and results.</p><p>The simulation of relay feedback control is shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. The results of ultimate gain and ultimate frequency are figured out by measuring the output amplitude and period as follows:<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\e8903f7f-5cb8-4929-bee0-03e1916a4a25.png" xlink:type="simple"/></inline-formula>;<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\051205c6-8f42-469c-8cec-09f7e8f42d42.png" xlink:type="simple"/></inline-formula>;<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\ed6737c5-d131-4427-8f68-4049cef6114a.png" xlink:type="simple"/></inline-formula>.</p><p>Therefore, the ultimate gain and ultimate frequency are:</p><p><inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\34d36800-ada0-41a4-991a-b1fc2d4280dd.png" xlink:type="simple"/></inline-formula>;<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\2b0bb075-5e09-45f6-aaa1-668f7b1aee21.png" xlink:type="simple"/></inline-formula>.</p><p>The learning rate and momentum gene of neural network are selected as followings:<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\a7419eb9-5cd7-46d7-9692-0f0fdfaecee1.png" xlink:type="simple"/></inline-formula>;<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\af401026-6e03-4d18-a396-189d7c5816d4.png" xlink:type="simple"/></inline-formula>.</p><p>Number of hidden neural is set as<inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\425dfd05-2683-413d-a57a-ac51aa28b1b1.png" xlink:type="simple"/></inline-formula>.</p><p>Sine wave with frequency of 0.25 Hz is used as reference profile for motion control system.</p><p>The disturbance is a random noise signal.</p><p>The simulation results including position tracking, position tracking error, and control input of proposed ATNNRF are shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>. The position tracking of in <xref ref-type="fig" rid="fig5">Figure 5</xref>(a) shows that the actual position of the control system has tracked the desired sine wave position, and the position tracking error in <xref ref-type="fig" rid="fig5">Figure 5</xref>(b) proves the robustness and high tracking ability of the proposed controller.</p><p>The simulated PID autotuning parameters and cost function of ATNNRF at 0.25 Hz are shown in <xref ref-type="fig" rid="fig6">Figure 6</xref>.</p><fig id="fig4"><label>Figure 4</label><caption><p> Simulation of Relay feedback control</p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\fe9420d2-7128-4747-853a-7837236cc07d.png"/></fig><fig-group id="fig5"><caption><title>Figure 5</title><p> Simulated response of ATNNRF control system at 0.25 Hz. (a) Position tracking; (b) Error; (c) Control input</p></caption><fig id ="fig5_1"><label>(a)</label><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\6682921f-e3e4-49f2-9266-6a50cadc6421.png"/></fig></fig-group><fig id="fig6"><label>Figure 6</label><caption><p> Simulated PID autotuning parameters of ATNNRF control system at 0.25 Hz. (a) Kp Gain; (b) Ki Gain; (c) Kd Gain; (d) Cost function of RBFNN</p></caption><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\34330ef2-d247-4390-9f28-f698c4d8b054.png"/></fig><p>The fast convergence of the simulated cost function <inline-formula><inline-graphic xlink:href="http://file.scirp.org/Html/htmlimages\3-7900323x\5cce1951-c01e-4480-8f63-cf6a14688ce3.png" xlink:type="simple"/></inline-formula> of ATNNRF has proven the effectiveness of the proposed controller scheme.</p></sec><sec id="s4"><title>4. Conclusion</title><p>This paper successfully demonstrates the application of an autotuning of PID controller based on neural network and relay feedback approach. First, the principle of relay feedback control is introduced to find the ultimate gain and ultimate frequency for the initial conditions of the neural network controller. Then, the results of relay feedback control are incorporated with the neural network to adjust the PID gains online through the gradient descent algorithm. An embedded motion control board with high performance DSP has been developed to implement the proposed control algorithms. The simulation results have proven the robustness and high tracking ability of the proposed controller.</p></sec><sec id="s5"><title>Acknowledgements</title><p>This work was supported by Dong-Eui University Foundation Grant (grant number 2014AA441).</p></sec></body><back><ref-list><title>References</title><ref id="scirp.48455-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">DESBOROUGH, L. AND MILLER, R. (2002) INCREASING CUSTOMER VALUE OF INDUSTRIAL CONTROL PERFORMANCE MONITORING— HONEYWELL EXPERIENCE. IN: RAWLINGS, J.B., OGUNNAIKE, B.A. AND EATON, J.W., EDS., 6TH INTERNATIONAL CONFERENCE ON CHEMICAL PROCESS CONTROL, AICHE SYMP., SERIES 326, AICHE, NEW YORK, 172-192.</mixed-citation></ref><ref id="scirp.48455-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">LEVY, S., KOROTKIN, S., HADAD, K., ELLENBOGEN, A., ARAD, M. AND KADMON, Y. (2012) PID AUTOTUNING USING RELAY FEEDBACK. IEEE 27TH CONVENTION OF ELECTRICAL &amp; ELECTRONICS ENGINEERS IN ISRAEL (IEEEI), EILAT, 14-17 NOVEMBER 2012, 1-4.</mixed-citation></ref><ref id="scirp.48455-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">TANG, W., WANG, Z.F., ZHANG, J.G. AND WANG, M.X. (2010) A RELAY FEEDBACK BASED AUTOTUNING PID CONTROLLER AND ITS APPLICATION TO PULP CONSISTENCY CONTROL. 29TH CHINESE CONTROL CONFERENCE (CCC), BEIJING, 29-31 JULY 2010, 3406-3410.</mixed-citation></ref><ref id="scirp.48455-ref4"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>TANG</surname><given-names> W.</given-names></name>,<name name-style="western"><surname> SHI</surname><given-names> S.J. </given-names></name>,<name name-style="western"><surname> WANG</surname><given-names> M.X. </given-names></name>,<etal>et al</etal>. (<year>2003</year>)<article-title>TANG, W., SHI, S.J. AND WANG, M.X.  RELAY FEEDBACK BASED AUTOTUNING PID CONTROL OF PAPER BASIS WEIGHT IN PAPER MAKING PROCESS</article-title><source> JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS</source><volume> 14</volume>,<fpage> 63</fpage>-<lpage>69</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.48455-ref5"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>MAJHI</surname><given-names> S. </given-names></name>,<name name-style="western"><surname> ATHERTON</surname><given-names> D.P. </given-names></name>,<etal>et al</etal>. (<year>1999</year>)<article-title>MAJHI, S. AND ATHERTON, D.P.  AUTOTUNING AND CONTROLLER DESIGN FOR PROCESSES WITH SMALL TIME DELAYS</article-title><source> IEE PROCEEDINGS OF CONTROL THEORY AND APPLICATIONS</source><volume> 146</volume>,<fpage> 415</fpage>-<lpage>425</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.48455-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">ZHU, A.-D., XING, J.-G. AND ZHAO, W.-M. (2010) AN AUTOTUNING PID CONTROL ALGORITHM FOR SENSORLESS DC MOTORS. INTERNATIONAL CONFERENCE ON ELECTRICAL AND CONTROL ENGINEERING (ICECE), WUHAN, 25-27 JUNE 2010, 369-372.</mixed-citation></ref><ref id="scirp.48455-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">VRANA, S. AND SULC, B. (2011) PID CONTROLLER AUTOTUNING BASED ON NONLINEAR TUNING RULES. 12TH INTERNATIONAL CARPATHIAN CONTROL CONFERENCE (ICCC), VELKE KARLOVICE, 25-28 MAY 2011, 443-446.</mixed-citation></ref><ref id="scirp.48455-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">ALAVI, S.M.S., AKBARZADEH, A. AND FARUGHIAN, A. (2011) AUTO-TUNING SMITH-PREDICTIVE CONTROL OF THREE-TANKS SYSTEM BASED ON MODEL REFERENCE ADAPTIVE SYSTEM. 2ND INTERNATIONAL CONFERENCE ON CONTROL, INSTRUMENTATION AND AUTOMATION (ICCIA), SHIRAZ, 27-29 DECEMBER 2011, 711-714.</mixed-citation></ref><ref id="scirp.48455-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">TOLEDO, C.F.M., LIMA, J.M.G. AND DA SILVA ARANTES, M. (2012) A MULTI-POPULATION GENETIC ALGORITHM APPROACH FOR PID CONTROLLER AUTO-TUNING. IEEE 17TH CONFERENCE ON EMERGING TECHNOLOGIES &amp; FACTORY AUTOMATION (ETFA), KRAKOW, 17-21 SEPTEMBER 2012, 1-8.</mixed-citation></ref><ref id="scirp.48455-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">ALAVIJE, H.S. AND AKHBARI, M. (2013) AUTO-TUNING PID CONTROLLER FOR LOW COST FAULT TOLERANT MOTOR DRIVE OF ELECTRIC VEHICLE. 4TH POWER ELECTRONICS, DRIVE SYSTEMS AND TECHNOLOGIES CONFERENCE (PEDSTC), TEHRAN, 13-14 FEBRUARY 2013, 38-43.</mixed-citation></ref><ref id="scirp.48455-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">WU, Y.-J., JUNG, Y.-G. AND LIM, Y.-C. (2012) AUTO-TUNING FUZZY PD CONTROL SCHEME FOR OUTPUT VOLTAGE CONTROL OF THREE-PHASE Z-SOURCE INVERTER. IEEE INTERNATIONAL SYMPOSIUM ON INDUSTRIAL ELECTRONICS (ISIE), HANGZHOU, 28-31 MAY 2012, 222-227.</mixed-citation></ref><ref id="scirp.48455-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">BLONDIN, M.-J. AND SICARD, P. (2013) PID CONTROLLERS AND ANTI-WINDUP SYSTEMS TUNING USING ANT COLONY OPTIMIZATION. 15TH EUROPEAN CONFERENCE ON POWER ELECTRONICS AND APPLICATIONS (EPE), LILLE, 2-6 SEPTEMBER 2013, 1-10.</mixed-citation></ref><ref id="scirp.48455-ref13"><label>13</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>CROWE</surname><given-names> J. </given-names></name>,<name name-style="western"><surname> JOHNSON</surname><given-names> M.A. </given-names></name>,<etal>et al</etal>. (<year>1999</year>)<article-title>CROWE, J. AND JOHNSON, M.A.  A NEW NONPARAMETRIC IDENTIFICATION PROCEDURE FOR ONLINE CONTROLLER TUNING</article-title><source> PROCEEDINGS OF THE 1999 AMERICAN CONTROL CONFERENCE</source><volume> 5</volume>,<fpage> 3337</fpage>-<lpage>3341</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.48455-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">CHEON, Y. LEE, D., LEE, I.-B. AND SUNG, S.W. (2013) A NEW PID AUTO-TUNING STRATEGY WITH OPERATIONAL OPTIMIZATION FOR MCFC SYSTEMS. 9TH ASIAN CONTROL CONFERENCE (ASCC), ISTANBUL, 23-26 JUNE 2013, 1-6.</mixed-citation></ref></ref-list></back></article>