\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Levenberg -Marquardt's Algorithm used for PID Controller Parameters Optimization},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{2015 2015-05-15 (revised: 27 Year 2015 15 May 2015)}
\def\TheID{\makeatother }
\def\TheDate{2015 2015-05-15}
\title{Levenberg -Marquardt's Algorithm used for PID Controller Parameters Optimization}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \author[1]{Ahmed Said Abd  El-hamid}

             \affil[1]{  National Research Center}

\renewcommand\Authands{ and }

\date{\small \em Received: 11 April 2015 Accepted: 3 May 2015 Published: 15 May 2015}

\maketitle


\begin{abstract}
        


The determination of parameters of controllers is an important problem in automatic control systems. In this paper, the Levenberg Marquardt (LM) Algorithm is used to effectively solve this problem with reasonable computational effort. The Levenberg Marquardt (LM) Algorithm for optimization of three term (PID) controller parameters with dynamic model of pH neutralization process is presented. The main goal is to show the merits of Levenberg Marquardt algorithm optimization and to determine its suitability in the area of control systems. Lastly, the application of this approach to the calculation of the parameters of PID controller shows that the Levenberg Marquardt (LM) algorithm has a better dynamic performance of pH neutralization process.

\end{abstract}


\keywords{parameters optimization, PID controller, pH neutralization process, levenberg marquardt (LM) algorithm.}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\let\tabcellsep& 	 	 		 
\section[{Introduction}]{Introduction}\par
n mathematics and computing, the Levenberg-Marquardt Algorithm (LMA), also known as the damped least-squares (DLS) method, is used to solve non-linear least squares problems. These minimization problems arise especially in least squares curve fitting. The LMA interpolates between the Gauss-Newton Algorithm (GNA) and the method of gradient descent.\par
The LMA is more robust than the GNA, which means that in many cases it finds a solution even if it starts very far off the final minimum. For well-behaved functions and reasonable starting parameters, the LMA tends to be a bit slower than the GNA. LMA can also be viewed as Gauss-Newton using a trust region approach. The LMA is a very popular curve-fitting algorithm used in many software applications for solving generic curve-fitting problems. However, as for many fitting algorithms, the LMA finds only a local minimum, which is not necessarily the global minimum \hyperref[b0]{[1]}. The primary application of the Levenberg-Marquardt algorithm is in the least squares curve fitting problem: given a set of ?? empirical datum pairs of independent and dependent variables, (????, ????) optimize the parameters ?? of the model curve ð??"ð??"(??, ??) so that the sum of the squares of the deviations becomes minimal.??(??) = ?[?? ?? ? ð??"ð??"(?? ?? , ??)] 2 ?? ??=1\par
The Levenberg-Marquardt (LM) algorithm is an iterative technique that locates the minimum of a multivariate function that is expressed as the sum of squares of non-linear real-valued functions \hyperref[b3]{[4,}\hyperref[b5]{6]}.\par
In the other hand, a common problem in control system design is establishing the appropriate value of controller gains. In general a low value of gain produces a slow system response, while high gain values can cause an excessively-oscillatory response with the possibility of instability. Somewhere between these extremes is a value of gain that produces the best system response. The essential function of a feedback control system is to reduce the error, ??(??) between any variable and its demanded value to zero as quickly as possible. Therefore, any criterion used to measure the quality of system response must take into account the variation of the error over the whole range of time. Four basic criteria are in common use:Integral of absolute error (IAE) = ? |??(??)|???? ? 0 Integral of squared error (ISE) = ? \{??(??)\} 2 ? 0 ???? Integral of time multiplied by absolute error (ITAE) = ? ?? ? 0 |??(??)|???? Integral of time multiplied by squared error (ISE) = ? ?? ? 0 \{??(??)\} 2 ????\par
For any of the possible criteria, the best response corresponds to the minimum value of the chosen criterion. Note that in all cases it is either the absolute error or the squared error which is involved, straightforward integration of the error would produce zero result, even if the system response was a constant amplitude oscillation. IAE is often used where digital simulation of a system is being employed, but it is inapplicable for analytical work, because the absolute value of an error function is not generally analytic in form. This problem is overcome by the ISE criterion. The ITAE and ITSE have an additional time multiplier of the error function, which emphasizes long-duration errors, and therefore these criteria are most often applied in systems requiring a fast settling time \hyperref[b2]{[3,}\hyperref[b4]{5]}.  
\section[{C}]{C}\par
In this paper, the shape of the complete closed loop response, from time ?? = 0, until steady state has been reached, could be used for the formulation of a dynamic performance criterion. The simple criteria of this category are based on the entire response of the process and the integral of the Square Error (ISE) criterion used here, where?????? = ? ?? 2 (??) ? 0 ???? = ? ?Y sp (t)-Y(t)? 2 ? 0 ???? (1)\par
Where ??(??) = ?? ???? (??) -??(??) is the deviation (error) of the response from the desired set point.\par
The ideal continuous time domain PID controller for a SISO process is expressed in the Laplace domain as follows:??(??) = ?? ?? + ?? ?? ?? + ?? ?? ??\textbf{(2)}\par
With K p = proportional gain, K i = integral time constant and K d = derivative time constant. 
\section[{II.}]{II.} 
\section[{Mathematical Model of ph Neutralization Process}]{Mathematical Model of ph Neutralization Process}\par
Consider a pH neutralization process as shown in Fig.  {\ref 1}. The flow rates of acid, buffer, base and effluent streams are denoted by ?? 1 , ?? 2 , ?? 3 , and ?? 4 , respectively. Output of the process is the pH value of the effluent stream, and the flow rate of base stream, ?? 3 is the control input. A dynamic model is derived using the conservation laws and reactions equilibrium. The modeling assumptions include perfect mixing, constant volume of the neutralization tank ( V ), and complete solubility of the ions involved. The chemical reactions in the system are as follows \hyperref[b7]{[8]}: The equilibrium constants for these reactions are:???? 1 = [?????? 3 ? ][?? + ] [?? 2 ???? 3 ] ???? 2 = [???? 3 2? ][?? + ] [?????? 3 ? ] ???? = [?? + ] + [???? ? ] + [???? 3 2? ]\par
The chemical equilibrium equations are modeled using the reaction invariant concept. For this system, concentrations of reaction invariants are defined as:?? 1 = [???? 3 ? ] ?? 2 = [???? + ] ?? 3 = [?? 2 ???? 3 ] + [?????? 3 ? ] + [???? 3 2? ]\par
Denoting, ?? = ???? the ions neutrality balance in the tank results the following static equation:?(??, ??) = ??? 1 + ?? 2 + ?? 3 ?? ??3 + 10 ??? ? 10 ??????? ?? = 0 (3) ?? ??3 = 2 + 10 ???? 2 ??? 1 + 10 ???? 2 ??? + 10 ???? 1 +???? 2 ?2?? ???? 1 = ? log 10 ???? 1 ???? 2 = ? log 10 ???? 2\par
The dynamic equations are given by:?????1 ???? = ?? 1 ?? (?? 11 ? ?? 1 ) + ?? 2 ?? (?? 21 ? ?? 1 ) + ?? 3 ?? (? 1 ? ?? 1 ) (4) ?????2 ???? = ?? 1 ?? (?? 12 ? ?? 2 ) + ?? 2 ?? (?? 22 ? ?? 2 ) + ?? 3 ?? (? 2 ? ?? 2 ) (5) ?????3 ???? = ?? 1 ?? (?? 13 ? ?? 3 ) + ?? 2 ?? (?? 23 ? ?? 3 ) + ?? 3 ?? (? 3 ? ?? 3 ) (6)\par
Where: © 
\section[{Simulation Results}]{Simulation Results}\par
The closed loop control system was solved using Levenberg-Marquardt's optimization approach with sampling time of 0.001 s. The simulation method combines SIMULINK module for pH neutralization model and M-file for LMA approach. A list of M-file programs used in the paper is provided in Appendix 1 and 2. Figure  {\ref (1 and 2)} show the responses of the pH neutralization obtained with change in pH set point. The optimal gains of PID controller are calculated to minimize the error function which described in equation \hyperref[b1]{(2)}. Also the values of gains of PID controller are plotted. The response of pH value tends to set point value with minimum steady state error and the values of gains tend to minimum values once the error reaches to zero value. Figure \hyperref[b2]{(3)} show the change in set point from pH = 6 to 7, noticed that the spike value of manipulated variable (q 3 ) and the values of gains (Kp and Ki) tend to increasing while the gain (Kd) tends to decreasing. In other hand, figure \hyperref[b3]{(4)} show the change in set point from 7 to 6, the value of manipulated variable (q 3 ) and the value of gains (Kp, Ki and Kd) tend to decreasing.    
\section[{Conclusion}]{Conclusion}\par
The paper presents an application of the Levenberg-Marquardt Algorithm (LMA) to optimization of parameters (Kp, Ki, Kd) of the PID controller structure according to minimum of integral square error. The simulated results were obtained of parameters by means of computer program implemented in Matlab software. As an example, the optimization of parameters of PID controllers with reference to a ph Neutralization process was presented. \begin{figure}[htbp]
\noindent\textbf{12}\includegraphics[]{image-2.png}
\caption{\label{fig_1}Fig. 1 : 2 ?}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-3.png}
\caption{\label{fig_2}}\end{figure}
    		 		\backmatter  			 \par
Appendix (1)   \%Determination of PID controller prameters using "Levenberg-Marquardt \%Algorithm" \%Step 1: Write an M-file tracklsq.m. function [Kp,Ki,Kd] = phpidL2015k ph2015 \% Load the simulink model   {\ref plot(t,y1,'r',t,y2,'-.b','LineWidth',2}); title('levenberg-marquardt') legend(\{'pH','pH setpoint'\},'Location','SE','FontSize',8') xlabel('Time (min)','FontSize',12') ylabel('pH','FontSize',12) grid on subplot(2,1,2),plot(t,y3,'r','LineWidth',2); xlabel('Time (min)','FontSize',12') ylabel('q3 (ml/s)','FontSize',12') Appendix (2) \% Read Mat. Files of PID gains \% ki(First number, second number)first number points to raw number \% second number points to col. number \% load ki.mat; \%\# assume this contains a matrix called ki load kp.mat load kd.mat for i=1:1:51 \%\# numbers of columns x(i) = ki(1,i); \%\# numbers from all rows, column 1, into X x1 			 			  				\begin{bibitemlist}{1}
\bibitem[Meng and Song ()]{b7}\label{b7} 	 		‘Fast Genetic Algorithms Used for PID Parameter Optimization’.  		 			X Meng 		,  		 			Song 		.  	 	 		\textit{Proceedings of IEEE International Conference on Automatic and Logistic},  				 (IEEE International Conference on Automatic and Logistic)  		2007. p. .  	 
\bibitem[Fuzzy-neural predictive control using Levenberg-Marquardt optimization approach]{b0}\label{b0} 	 		\textit{Fuzzy-neural predictive control using Levenberg-Marquardt optimization approach},  		 	 
\bibitem[Mendis et al. ()]{b6}\label{b6} 	 		‘Learning Generalized Weighted Relevance Aggregation Operators Using Levenberg-Marquardt’.  		 			B S U Mendis 		,  		 			. T D Gedeon 		,  		 			. L Kóczy 		.  	 	 		\textit{Proc. of the Sixth IEEE International Conference on Hybrid Intelligent Systems},  				 (of the Sixth IEEE International Conference on Hybrid Intelligent Systems)  		2006.  	 
\bibitem[Cominos and Munro (2002)]{b2}\label{b2} 	 		‘PID Controllers: Recent Tuning Methods and Design to Specification’.  		 			P Cominos 		,  		 			N Munro 		.  	 	 		\textit{IEE Pro-ceedings Control Theory \& Applications}  		January 2002. 149  (I) .  	 
\bibitem[Turkoglu et al. ()]{b4}\label{b4} 	 		‘PID Parameter Optimization of an UAV Longitudinal Flight Control System’.  		 			Kamran Turkoglu 		,  		 			Ugur Ozdemir 		,  		 			Melike Nikbay 		,  		 			Elbrous M Jafarov 		.  	 	 		\textit{World Academy of Science, Engineering and Technology}  		2008. 21.  	 
\bibitem[Ayman and Aly ()]{b1}\label{b1} 	 		‘PID Parameters Optimization Using Genetic Algorithm Technique for Electrohydraulic Servo Control System’.  		 			A Ayman 		,  		 			Aly 		.  	 	 		\textit{Intelligent Control and Automation}  		2011. 2 p. .  	 
\bibitem[Kada and Ghazzawi (2011)]{b5}\label{b5} 	 		‘Robust PID Controller Design for an UAV Flight Control System’.  		 			B Kada 		,  		 			Y Ghazzawi 		.  	 	 		\textit{WCECS}  		2011. October 19-21, 2011.  	 
\bibitem[Wang et al. ()]{b3}\label{b3} 	 		‘Survey of Advanced PID Parameter Tuning Methods’.  		 			W Wang 		,  		 			J T Zhang 		,  		 			T Y A Chai 		.  	 	 		\textit{Acta Auto-matica}  		2000. 26  (3)  p. .  	 
\end{bibitemlist}
 			 		 	 
\end{document}