ASPECTS OF POLE PLACEMENT TECHNIQUE IN SYMMETRICAL OPTIMUM METHOD FOR PID CONTROLLER DESIGN
|
|
- Evan Golden
- 6 years ago
- Views:
Transcription
1 ASES OF OLE LAEMEN EHNIQUE IN SYMMERIAL OIMUM MEHOD FOR ID ONROLLER DESIGN Viorel Nicolau *, onstantin Miholca *, Dorel Aiordachioaie *, Emil eanga ** * Deartment of Electronics and elecommunications, Dunarea de Jos University of Galati, 47 Domneasca Street, Galati, 88, Galati, ROMANIA Viorel.Nicolau@ugal.ro, onstantin.miholca@ugal.ro, Dorel.Aiordachioaie@ugal.ro ** Research enter in Advanced rocess ontrol Systems, Dunarea de Jos University of Galati, 47 Domneasca Street, Galati, 88, Galati, ROMANIA Emil.eanga@ugal.ro Abstract: he aer resents an analytical aroach to the design of ID controllers by combining ole lacement with symmetrical otimum method, for the integration lus first-order lant model. he desired closed-loo transfer function (c.l.t.f.) contains a second-order oscillating system and a lead-delay comensator. It is shown that the ero value of c.l.t.f. deends on the real-ole value of c.l.t.f. and in addition, there is only one ole value, which satisfies the assumtions of symmetrical otimum method. In these conditions, the analytical eressions of the controller arameters can be simlified. he method is alied to design a ID autoilot for heading control of a shi with firstorder Nomoto model. oyright 5 IFA Keywords: ID control, ole lacement, symmetrical otimum method, autoilot. INRODUION Due to widesread industrial use of ID controllers, it is clear that even a small ercentage imrovement in ID design could have a major imact worldwide (Silva, et al., ). uning of ID controllers is a difficult task, as a three-arameter model should be defined and it must be accurate at higher frequencies (Astrom and Hagglund, 995). Although, analytical methods are more convenient than grahical methods based on frequency diagrams, in industry, most controllers are tuned using frequency resonse methods (ang and Ortega, 993). Analytical methods rely on low-order models characteried by a small number of arameters. he most emloyed models are the integration lus firstorder model, which is used for thermal and electromechanical rocesses, and the first-order lus dead-time model, which is used for chemical rocesses (Datta, et al., ). In this aer, the integration lus first-order model tye is used for shi dynamics modelling. Also, analytical design of ID autoilot for heading control is considered. If ole lacement method (M) is used to synthesie the ID controller, the first ste is to secify some erformance conditions of the closedloo system, which lead to the eression of the closed-loo transfer function (c.l.t.f.) (Yu and Salgado, 3). In this aer, the desired c.l.t.f. contains a second-order oscillating system and a ole-ero air, with real and negative values. Alying M, the ero value deends on the ole value, and the controller arameters deend on the arameters of c.l.t.f. he resulting oen-loo transfer function (o.l.t.f.) contains a double-integrator element and a ole-ero air, with real and negative values. Hence, the symmetrical otimum method (SOM) can be used (Kessler, 958).
2 Imosing symmetrical characteristics of the oenloo transfer function, the analytical eressions of the controller arameters can be simlified. he goal of this aer is to find the ole-ero values of c.l.t.f. and the simlified analytical eressions of ID controller arameters, which satisfy two simultaneous conditions: the desired close-loo transfer function and symmetrical characteristics of the oen-loo transfer function. he aer is organied as follows. Section rovides mathematical models used in simulations. In section 3, the analytical eressions of ID controller arameters are obtained using M. In section 4, the eressions of the controller arameters are simlified, imosing symmetrical characteristics of the o.l.t.f. Section 5 describes the simulation results, using a ID autoilot for heading control of a shi. onclusions are resented in section 6.. MAHEMAIAL MODELS onsider the classical structure of the control loo without disturbances, as shown in Fig.. he lant model contains an integrator and the controller is of ID tye. r + - ε ontroller (ID) Fig.. lassical structure of the control loo he erformance conditions of the closed-loo system can be secified imosing the eression of system transfer function. In general, a second order reference model is chosen to aroimate the behaviour of the closed-loo system: H (, () s + ζ s + where > is the natural frequency and ζ > is the daming coefficient. Because the lant model contains an integrator and another one is included into the ID controller, the oen-loo transfer function contains a doubleintegrator, which can not be obtained with c.l.t.f. given in (). herefore, the reference model must be comleted with a lead-delay comensator, which contains a ole-ero air, with real and negative values (eanga, et al., ): ( s + ) H s () ( s + ζ s + ) ( s + ) where > and >. u lant H( (with integrator) y he lant model contains an integrator, and it is characteried by a dominant time constant ( ) and a gain coefficient (k ). he eression of the model deends on the rocess tye. a) If the rocess is fast, then the small time constants can not be neglected and the model contains an equivalent small time constant ( Σ ), corresonding to the sum of arasitic time constants: H ( s ( s ( s +, (3) ) where Σ <. In this case, the ID controller is of the form: H (, (4) where < < Σ. ( s ( s he oen-loo transfer function is: ( ( (5) s ( s ( s Using ole cancellation, the non-ero dominant ole of the lant model is cancelled by choosing: (6) hus, only two controller arameters must be determined: k and. It can be observed that, if the rocess does not have any non-ero dominant ole (model without time constant ), then the controller is of I tye (without time constant ) and the same arameters must be determined (k and ). b) If the rocess is slow, the equivalent small time constant ( Σ ) can be neglected and the lant model is: H ( (7) s ( s he ID controller contains a sulementary degree of freedom and it is of the form: + H ( (, (8) s + where Σ < < <. he oen-loo transfer function is: + H ( ( (9) s + s ( s Again, using ole cancellation, the non-ero dominant ole of the lant model is cancelled, resulting equation (6):. In this case, three controller arameters must be determined: k, and. his is a more general case because the time constant is not imosed by the rocess and it can be chosen. So, in this aer, the lant model given in (7) is used for comutations, but discussions are made also for model given in (3).
3 In both cases, the oen-loo transfer functions, given in (5) and (9), have similar eressions: ( ( s H ( H ( () s With this oen-loo transfer function, the symmetrical otimum method (SOM) can be used. he time constant has different meanings: in the first case, it reresents the equivalent small time constant ( Σ ) imosed by the rocess, while in the second case, it is a controller arameter ( ). 3. ID ONROLLER DESIGN USING OLE LAEMEN MEHOD onsider the control system illustrated in Fig. with the desired closed-loo transfer function given in (). roosition : For every ole value (s -) of the desired closed-loo transfer function given in (), there is only one ero value (s -) for which the oen-loo transfer function has a double-ole in origin, and in addition, the ero frequency is smaller than the ole frequency: <. he roosition demonstration includes the net lemma results. Lemma : he necessary and sufficient condition, for the eistence of a double-ole in origin for the oen-loo transfer function of a control system illustrated in Fig., starting from a desired c.l.t.f. of the form given in (), is:,,, > ζ ζ () roof L: he transfer function of the oen-loo system can be comuted starting from the desired closed-loo transfer function (Nicolau, 4): H ( H ( H ( ζ s ζ s + s From () it can be observed that a double-ole in origin is obtained if the equation below is satisfied: ζ, (3) which is equivalent with: ζ + (4) ( s + ) () From (4) it results the necessary and sufficient condition indicated in () (q.e.d.). Imlicitly, the unique ero value results, whose eression deends on the selected ole value: s,,, ζ > (5) ζ For every frequency () of the ole, the corresonding frequency () of the ero is smaller: > < (6) ζ So, for the lead-delay comensator introduced into the desired c.l.t.f. given in (), the hase-lead effect is dominantly. In this case, the real values of ole-ero air and conjugate comle oles, of the desired c.l.t.f. given in (), are illustrated in Fig.. s o s ζ s s ζ Im Fig.. oles and ero of the desired c.l.t.f. Re Using () and (3) in (), the eression of the oen-loo transfer function is obtained: ( ζ ) s + ζ + H s (7) s s + + ζ [ ] Denote by and, resectively, the ero and ole frequencies of the oen-loo transfer function: ζ + (8) ζ he oen-loo transfer function can be rewritten: (ζ ) ( s + ) H ( (9) s [ s + ] utting into evidence the time constants, the oenloo transfer function can be rewritten, like in (): ζ s + H ( () + s (ζ ) + s ζ + From (), using (5) or (9) corresonding to the lant model indicated in (3) or (7), resectively, it results: + ( s s + s s + ζ s +, () s (ζ + ) s + ζ + where Σ or.
4 Equation () can be reduced to an equality of two olynomials of 3 rd order in s variable: (ζ + ) ζ + ( s ( s s + ζ + + s s s +, () where no ole cancellation was considered. he equality must be true for every frequency, resulting a four equation system: ( ζ + ) (3.) ζ (3.) ζ + + ( + ) ζ + ζ + + ( + ) (3.3) ζ (3.4) ζ + he solutions of the equation system are the ID controller arameters (Nicolau, 4): ζ k (4.) k ζ + ζ (4.) (4.3) (4.4) ζ + It can be observed that the solution (4.3) reresents the ole cancellation condition, considered in (6), and it does not deend on the ole value. If the rocess is slow and the equivalent small time constant ( Σ ) is ignored, the time constant reresents a controller arameter given in (4.4). he time constants must satisfy the inequalities: Σ < < <, (5) which can be transosed into frequency domain: < < < (6) herefore, the reference model must be chosen so that the arameters of c.l.t.f. (, ζ and ) to satisfy the system of inequalities: ζ + < (7.) ζ + > (7.) ζ + > (7.3) If the time constant Σ is imosed by the rocess and the lant model given in (3) is considered, then the solution (4.4) becomes:, (8) ζ + which reresents a sulementary condition for arameters of c.l.t.f. (, ζ and ). In addition, the time constants must satisfy the inequalities: < Σ <, (9) which can be transosed into frequency domain: < < (3) In this case, the reference model must be chosen so that the arameters of c.l.t.f. (, ζ and ) to satisfy the following conditions: ζ + (3.) ζ + > (3.) ζ > (3.3) It can be observed that the first condition in (3) is more restrictive than the corresonding one in (7), while the last two conditions are the same. 4. SYMMERIAL HARAERISIS OF OEN-LOO RANSFER FUNION he ID controller arameters deend on the arameters of c.l.t.f. (, ζ and ). In general, and ζ characterie the desired system behaviour and they have fied values, while the ole value can be chosen. Secific ole values can be imosed by using sulementary conditions. In this aer, the conditions for choosing the ole value refer to the symmetrical otimum method, which simlify the eressions of ID arameters. he goal is to find that ole value of the c.l.t.f., which satisfies the assumtions of symmetrical otimum method around natural frequency, for the transfer function of oen-loo system given in (9). Using this value, the eressions of ID arameters in (4) are simlified. roosition : here is only one admissible value for the ole (s - ) of c.l.t.f. given in (), so that the corresonding o.l.t.f. given in (7) to have symmetrical characteristics around :,, ζ (3), > roof: For the secified oen-loo transfer function, the symmetry of magnitude-frequency characteristic around natural frequency imlies the symmetry of hase-frequency characteristic. herefore, only the symmetry of former characteristic must be imosed.
5 he general form of the symmetrical otimum method imoses two conditions for magnitudefrequency characteristic: a) the central frequency must be equally laced between ero and ole frequencies on the -base logarithmic scale: (33) b) for central frequency, the magnitude-frequency characteristic of o.l.t.f. must have db: H ( j ) (34) Using (8) in (33), the first condition becomes: (ζ + ) (35) ζ From (35), it results:,,, ζ > (36) So, the condition (33) is satisfied if. For the second condition in (34), the magnitude of oen-loo transfer function in frequency is comuted from (9): (ζ ) + H ( j ) (37) + he frequencies and are relaced with their eressions from (8), resulting: + 4ζ + 4ζ H ( j ) (38) + 4ζ + 4ζ Using (38) in (34), the same solution in (36) is obtained: oncluding, there is only one admissible value for the ole of c.l.t.f., so that the corresonding o.l.t.f. to have symmetrical characteristics around (q.e.d.). From (), it results: (39) ζ + he eression of c.l.t.f. becomes: (ζ s + ζ + H s (4) ( s + ζ s + ) ( s + ) Also, from (8), the ero and ole frequencies of o.l.t.f. are obtained:, ( ζ (4) ζ + he oen-loo transfer function can be rewritten: (ζ s + ζ + H s (4) s s + (ζ [ ] he real values of ole-ero air and conjugate comle oles, of the c.l.t.f. given in (4), are illustrated in Fig. 3. s s s ζ o s ζ + Im ζ Re Fig. 3. oles and ero of the c.l.t.f. with he osition of the ero s deend on ζ + the arameter ζ: - if ζ,, then < < ζ and the ero is laced between the two oints: s and s ζ, resectively; - if ζ, then ζ. his is the articular case of the Kessler s symmetrical otimum method; - if ζ >, then > > ζ and the ero is laced to the right of the oint s ζ. his is the case illustrated in Fig. 3. Knowing the ole value of c.l.t.f. ( ), the ID controller arameters result from (4): (43.) ζ + (43.) (43.3) (43.4) (ζ he arameters in (43) corresond to the lant model given in (7) and ID controller given in (8). In this case, the conditions from (7) deend on the arameters and ζ. Hence, the reference model must be chosen so that the arameters of c.l.t.f. ( and ζ) to satisfy the system of inequalities: ζ + < (44.) Σ ζ + > (44.) ζ + > (44.3)
6 roosition 3: In the case of symmetrical characteristics of the o.l.t.f. given in (4) around the natural frequency, the hase margin and the distance between the frequency oints on the -base logarithmic scale deend only on the arameter ζ. roof: he distance between the frequency oints on the -base logarithmic scale can be easily obtained, using (4) in (33): ζ + (45) he hase margin is: ϕm π + arg( j) ) (46) Using the o.l.t.f. given in (4), results: ϕ m arctg(ζ arctg (47) ζ + It can be observed that, for articular case of the Kessler s symmetrical otimum method ( ζ. 5 ), the distance between frequency oints is equal with an octave and the hase margin is ϕ [deg]. 5. SIMULAION RESULS For simulations, the heading control roblem of a shi is considered, using a ID autoilot. he shi model is linear, being identified for a shi seed of knots (Nicolau, 4). It is a first order Nomoto model of the form given in (7): ψ ( H (, (48) δ ( s ( s where ψ( and δ( reresent the Lalace transforms of yaw angle and rudder angle, resectively. he shi model arameters are: k.834 [s - ], [s] (49) he autoilot model is given in (8) and the desired c.l.t.f. is given in (). he arameters and ζ are chosen from erformance conditions (Fossen, 994): ζ.9,. [rad/s] (5) Starting from the desired c.l.t.f. and imosing symmetrical characteristics of the o.l.t.f., the eressions in (4) and (4) are obtained. he ste resonse of the c.l.t.f. is illustrated in Fig. 4. From (43), the autoilot arameters are obtained: m k., 8 [s], [s], [s] Fig. 4. Ste resonse of the closed-loo system onsidering [s], the conditions in (44) are satisfied. he symmetrical characteristics of the o.l.t.f. are illustrated in Fig. 5. It can be observed that the hase margin is ϕ [deg]. Σ m Fig. 5. Symmetrical characteristics of the o.l.t.f. 6. ONLUSIONS here is only one ossible air for the ole-ero values of c.l.t.f. so that the corresonding arameters of ID controller to satisfy two simultaneous conditions: the desired behaviour of close-loo system and symmetrical characteristics of the oenloo transfer function. REFERENES Astrom, K. and. Hagglund (995). ID ontrollers: heory, Design, and uning. Research riangle ark, N: Instrument Society of America. eanga, E,. Nichita, L. rotin and N. A. utululis (). heorie de la ommande des Systemes. Ed. ehnica, Bucharest, Romania. Datta, A., M.. Ho and S.. Bhattacharyya (). Structure and Synthesis of ID ontrollers. Londin, UK: Sringer-Verlag. Fossen,. I. (994). Guidance and ontrol of Ocean Vehicles. John Wiley and Sons Ltd, NY. Kessler,. (958). Das Symmetrische Otimum. Regelungstechnik, 6, and Nicolau, V. (4). ontributions in Advanced Automatic ontrol of Naval Systems. h.d.hesis (in romanian), University of Galati, Romania. Silva, G. J., A. Datta and S.. Bhattacharyya (). New Results on the Synthesis of ID ontrollers. IEEE ransactions on Automatic ontrol, 47, no., ang, Y. and R. Ortega (993). Adative uning to Frequency Resonse Secification. Automatica, 9, Yu, J. I. and M. E. Salgado (3). From lassical to State-Feedback-Based ontrollers. IEEE ontrol Systems Magaine, 3, no. 4,
Research Article On PID Controller Design by Combining Pole Placement Technique with Symmetrical Optimum Criterion
Mathematical Problems in Engineering Volume 23, Article ID 36827, 8 pages http://dx.doi.org/.55/23/36827 Research Article On PID Controller Design by Combining Pole Placement Technique with Symmetrical
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationLINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL
LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL Mohammad Bozorg Deatment of Mechanical Engineering University of Yazd P. O. Box 89195-741 Yazd Iran Fax: +98-351-750110
More informationRobust Performance Design of PID Controllers with Inverse Multiplicative Uncertainty
American Control Conference on O'Farrell Street San Francisco CA USA June 9 - July Robust Performance Design of PID Controllers with Inverse Multilicative Uncertainty Tooran Emami John M Watkins Senior
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationMODELLING, SIMULATION AND ROBUST ANALYSIS OF THE TEMPERATURE PROCESS CONTROL
The 6 th edition of the Interdiscilinarity in Engineering International Conference Petru Maior University of Tîrgu Mureş, Romania, 2012 MODELLING, SIMULATION AND ROBUST ANALYSIS OF THE TEMPERATURE PROCESS
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationFeedback-Based Iterative Learning Control for MIMO LTI Systems
International Journal of Control, Feedback-Based Automation, Iterative and Systems, Learning vol. Control 6, no., for. MIMO 69-77, LTI Systems Aril 8 69 Feedback-Based Iterative Learning Control for MIMO
More informationRobust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System
Vol.7, No.7 (4),.37-38 htt://dx.doi.org/.457/ica.4.7.7.3 Robust Predictive Control of Inut Constraints and Interference Suression for Semi-Trailer System Zhao, Yang Electronic and Information Technology
More informationOn Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.**
On Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.** * echnology and Innovation Centre, Level 4, Deartment of Electronic and Electrical Engineering, University of Strathclyde,
More informationDesign of a Heading Autopilot for Mariner Class Ship with Wave Filtering Based on Passive Observer
Design of a Heading Autopilot for Mariner Class Ship with Wave Filtering Based on Passive Observer 1 Mridul Pande, K K Mangrulkar 1, Aerospace Engg Dept DIAT (DU), Pune Email: 1 mridul_pande000@yahoo.com
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationOil Temperature Control System PID Controller Algorithm Analysis Research on Sliding Gear Reducer
Key Engineering Materials Online: 2014-08-11 SSN: 1662-9795, Vol. 621, 357-364 doi:10.4028/www.scientific.net/kem.621.357 2014 rans ech Publications, Switzerland Oil emerature Control System PD Controller
More informationMinimax Design of Nonnegative Finite Impulse Response Filters
Minimax Design of Nonnegative Finite Imulse Resonse Filters Xiaoing Lai, Anke Xue Institute of Information and Control Hangzhou Dianzi University Hangzhou, 3118 China e-mail: laix@hdu.edu.cn; akxue@hdu.edu.cn
More informationROBUST ADAPTIVE MOTION CONTROL WITH ENVIRONMENTAL DISTURBANCE REJECTION FOR PERTURBED UNDERWATER VEHICLES
Journal of Marine Science and echnology DOI: 0.69/JMS-03-05-3 ROBUS ADAPIVE MOION CONROL WIH ENVIRONMENAL DISURBANCE REJECION FOR PERURBED UNDERWAER VEHICLES Hamid Reza Koofigar his article has been eer
More informationExpansion in n 1 for percolation critical values on the n-cube and Z n : the first three terms
Eansion in n for ercolation critical values on the n-cube and Z n : the first three terms Remco van der Hofstad Gordon Slade December 22, 2003 Abstract Let c (Q n ) and c (Z n ) denote the critical values
More informationL p Norms and the Sinc Function
L Norms and the Sinc Function D. Borwein, J. M. Borwein and I. E. Leonard March 6, 9 Introduction The sinc function is a real valued function defined on the real line R by the following eression: sin if
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationPosition Control of Induction Motors by Exact Feedback Linearization *
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 8 No Sofia 008 Position Control of Induction Motors by Exact Feedback Linearization * Kostadin Kostov Stanislav Enev Farhat
More informationME scope Application Note 16
ME scoe Alication Note 16 Integration & Differentiation of FFs and Mode Shaes NOTE: The stes used in this Alication Note can be dulicated using any Package that includes the VES-36 Advanced Signal Processing
More informationarxiv: v1 [quant-ph] 20 Jun 2017
A Direct Couling Coherent Quantum Observer for an Oscillatory Quantum Plant Ian R Petersen arxiv:76648v quant-h Jun 7 Abstract A direct couling coherent observer is constructed for a linear quantum lant
More informationParticipation Factors. However, it does not give the influence of each state on the mode.
Particiation Factors he mode shae, as indicated by the right eigenvector, gives the relative hase of each state in a articular mode. However, it does not give the influence of each state on the mode. We
More informationShadow Computing: An Energy-Aware Fault Tolerant Computing Model
Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms
More information#A47 INTEGERS 15 (2015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS
#A47 INTEGERS 15 (015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS Mihai Ciu Simion Stoilow Institute of Mathematics of the Romanian Academy, Research Unit No. 5,
More informationControllability and Resiliency Analysis in Heat Exchanger Networks
609 A ublication of CHEMICAL ENGINEERING RANSACIONS VOL. 6, 07 Guest Editors: Petar S Varbanov, Rongxin Su, Hon Loong Lam, Xia Liu, Jiří J Klemeš Coyright 07, AIDIC Servizi S.r.l. ISBN 978-88-95608-5-8;
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationOn Character Sums of Binary Quadratic Forms 1 2. Mei-Chu Chang 3. Abstract. We establish character sum bounds of the form.
On Character Sums of Binary Quadratic Forms 2 Mei-Chu Chang 3 Abstract. We establish character sum bounds of the form χ(x 2 + ky 2 ) < τ H 2, a x a+h b y b+h where χ is a nontrivial character (mod ), 4
More informationImproved Bounds on Bell Numbers and on Moments of Sums of Random Variables
Imroved Bounds on Bell Numbers and on Moments of Sums of Random Variables Daniel Berend Tamir Tassa Abstract We rovide bounds for moments of sums of sequences of indeendent random variables. Concentrating
More informationDesign of NARMA L-2 Control of Nonlinear Inverted Pendulum
International Research Journal of Alied and Basic Sciences 016 Available online at www.irjabs.com ISSN 51-838X / Vol, 10 (6): 679-684 Science Exlorer Publications Design of NARMA L- Control of Nonlinear
More informationThe inverse Goldbach problem
1 The inverse Goldbach roblem by Christian Elsholtz Submission Setember 7, 2000 (this version includes galley corrections). Aeared in Mathematika 2001. Abstract We imrove the uer and lower bounds of the
More informationMODULAR LINEAR TRANSVERSE FLUX RELUCTANCE MOTORS
MODULAR LINEAR TRANSVERSE FLUX RELUCTANCE MOTORS Dan-Cristian POPA, Vasile IANCU, Loránd SZABÓ, Deartment of Electrical Machines, Technical University of Cluj-Naoca RO-400020 Cluj-Naoca, Romania; e-mail:
More informationNotes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle
McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8
More informationGeneralized analysis method of engine suspensions based on bond graph modeling and feedback control theory
Generalized analysis method of ine susensions based on bond grah modeling and feedback ntrol theory P.Y. RICHARD *, C. RAMU-ERMENT **, J. BUION *, X. MOREAU **, M. LE FOL *** *Equie Automatique des ystèmes
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
More informationANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE
ANALYTICAL MODEL FOR THE BYPASS ALE IN A LOOP HEAT PIPE Michel Seetjens & Camilo Rindt Laboratory for Energy Technology Mechanical Engineering Deartment Eindhoven University of Technology The Netherlands
More informationThe Fekete Szegő theorem with splitting conditions: Part I
ACTA ARITHMETICA XCIII.2 (2000) The Fekete Szegő theorem with slitting conditions: Part I by Robert Rumely (Athens, GA) A classical theorem of Fekete and Szegő [4] says that if E is a comact set in the
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationMultivariable PID Control Design For Wastewater Systems
Proceedings of the 5th Mediterranean Conference on Control & Automation, July 7-9, 7, Athens - Greece 4- Multivariable PID Control Design For Wastewater Systems N A Wahab, M R Katebi and J Balderud Industrial
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationRUN-TO-RUN CONTROL AND PERFORMANCE MONITORING OF OVERLAY IN SEMICONDUCTOR MANUFACTURING. 3 Department of Chemical Engineering
Coyright 2002 IFAC 15th Triennial World Congress, Barcelona, Sain RUN-TO-RUN CONTROL AND PERFORMANCE MONITORING OF OVERLAY IN SEMICONDUCTOR MANUFACTURING C.A. Bode 1, B.S. Ko 2, and T.F. Edgar 3 1 Advanced
More informationSome Unitary Space Time Codes From Sphere Packing Theory With Optimal Diversity Product of Code Size
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 5, NO., DECEMBER 4 336 Some Unitary Sace Time Codes From Shere Packing Theory With Otimal Diversity Product of Code Size Haiquan Wang, Genyuan Wang, and Xiang-Gen
More informationE-companion to A risk- and ambiguity-averse extension of the max-min newsvendor order formula
e-comanion to Han Du and Zuluaga: Etension of Scarf s ma-min order formula ec E-comanion to A risk- and ambiguity-averse etension of the ma-min newsvendor order formula Qiaoming Han School of Mathematics
More informationResearch Article Estimation of Observer Parameters for Dynamic Positioning Ships
Mathematical Problems in Engineering Volume 13, Article ID 17363, 7 ages htt://dx.doi.org/1.1155/13/17363 Research Article Estimation of Observer Parameters for Dynamic Positioning Shis Xiaogong Lin, Yehai
More informationMagnetostrictive Dynamic Vibration Absorber (DVA) for Passive and Active Damping
aer : 59 /. Magnetostrictive ynamic Vibration Absorber (VA) for Passive and Active aming C. May a,. uhnen b, P. Pagliarulo b and H. Janocha b a Centre for nnovative Production (ZP), Saarbrücken, ermany,
More informationRadar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D.
Radar Dish ME 304 CONTROL SYSTEMS Mechanical Engineering Deartment, Middle East Technical University Armature controlled dc motor Outside θ D outut Inside θ r inut r θ m Gearbox Control Transmitter θ D
More informationFUZZY CONTROL. Khurshid Ahmad, Professor of Computer Science, Department of Computer Science Trinity College, Dublin-2, IRELAND Nov 15-16th, 2011
: Mamdani & akagi-sugeno Controllers Khurshid Ahmad, Professor of Comuter Science, Deartment of Comuter Science rinit College, Dublin-, IRELAND Nov 5-6th, 0 htts://www.cs.tcd.ie/khurshid.ahmad/eaching/eaching.html
More informationRecovering preferences in the household production framework: The case of averting behavior
Udo Ebert Recovering references in the household roduction framework: The case of averting behavior February 2002 * Address: Deartment of Economics, University of Oldenburg, D-26 Oldenburg, ermany Tel.:
More informationAnytime communication over the Gilbert-Eliot channel with noiseless feedback
Anytime communication over the Gilbert-Eliot channel with noiseless feedback Anant Sahai, Salman Avestimehr, Paolo Minero Deartment of Electrical Engineering and Comuter Sciences University of California
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
More informationMODELING OF THE DYNAMIC RESPONSE OF A FRANCIS TURBINE
MODELING OF THE DYNAMIC RESPONSE OF A FRANCIS TURBINE Paolo PENNACCHI, Steven CHATTERTON *, Andrea VANIA Politecnico di Milano, Deartment of Mechanical Engineering - Via G. La Masa, - 256 - Milano, Italy.
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationOptimal array pattern synthesis with desired magnitude response
Otimal array attern synthesis with desired magnitude resonse A.M. Pasqual a, J.R. Arruda a and P. erzog b a Universidade Estadual de Caminas, Rua Mendeleiev, 00, Cidade Universitária Zeferino Vaz, 13083-970
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationCOMMUNICATION BETWEEN SHAREHOLDERS 1
COMMUNICATION BTWN SHARHOLDRS 1 A B. O A : A D Lemma B.1. U to µ Z r 2 σ2 Z + σ2 X 2r ω 2 an additive constant that does not deend on a or θ, the agents ayoffs can be written as: 2r rθa ω2 + θ µ Y rcov
More informationPROFIT MAXIMIZATION. π = p y Σ n i=1 w i x i (2)
PROFIT MAXIMIZATION DEFINITION OF A NEOCLASSICAL FIRM A neoclassical firm is an organization that controls the transformation of inuts (resources it owns or urchases into oututs or roducts (valued roducts
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationOptimal Recognition Algorithm for Cameras of Lasers Evanescent
Otimal Recognition Algorithm for Cameras of Lasers Evanescent T. Gaudo * Abstract An algorithm based on the Bayesian aroach to detect and recognise off-axis ulse laser beams roagating in the atmoshere
More informationME 375 System Modeling and Analysis. Homework 11 Solution. Out: 18 November 2011 Due: 30 November 2011 = + +
Out: 8 November Due: 3 November Problem : You are given the following system: Gs () =. s + s+ a) Using Lalace and Inverse Lalace, calculate the unit ste resonse of this system (assume zero initial conditions).
More informationOnline Appendix to Accompany AComparisonof Traditional and Open-Access Appointment Scheduling Policies
Online Aendix to Accomany AComarisonof Traditional and Oen-Access Aointment Scheduling Policies Lawrence W. Robinson Johnson Graduate School of Management Cornell University Ithaca, NY 14853-6201 lwr2@cornell.edu
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationOptimal Integrated Control and Scheduling of Systems with Communication Constraints
Otimal Integrated Control and Scheduling of Systems with Communication Constraints Mohamed El Mongi Ben Gaid, Arben Çela and Yskandar Hamam Abstract This aer addresses the roblem of the otimal control
More informationInequalities for the L 1 Deviation of the Empirical Distribution
Inequalities for the L 1 Deviation of the Emirical Distribution Tsachy Weissman, Erik Ordentlich, Gadiel Seroussi, Sergio Verdu, Marcelo J. Weinberger June 13, 2003 Abstract We derive bounds on the robability
More informationAn analytical approximation method for the stabilizing solution of the Hamilton-Jacobi equation based on stable manifold theory
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July -3, 27 ThA8.4 An analytical aroimation method for the stabilizing solution of the Hamilton-Jacobi
More informationAnalysis of Fractional order PID controller for Ceramic Infrared Heater
06 IJEDR Volume 4, Issue ISSN: 3-9939 Analysis of Fractional order PID controller for Ceramic Infrared Heater Vineet Shekher, V.S. Guta, Sumit Saroha EEE, Deartment SRM, University, NCR Camus, Ghaziabad,
More informationA Social Welfare Optimal Sequential Allocation Procedure
A Social Welfare Otimal Sequential Allocation Procedure Thomas Kalinowsi Universität Rostoc, Germany Nina Narodytsa and Toby Walsh NICTA and UNSW, Australia May 2, 201 Abstract We consider a simle sequential
More informationAn Estimate For Heilbronn s Exponential Sum
An Estimate For Heilbronn s Exonential Sum D.R. Heath-Brown Magdalen College, Oxford For Heini Halberstam, on his retirement Let be a rime, and set e(x) = ex(2πix). Heilbronn s exonential sum is defined
More informationHotelling s Two- Sample T 2
Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test
More informationDecoding First-Spike Latency: A Likelihood Approach. Rick L. Jenison
eurocomuting 38-40 (00) 39-48 Decoding First-Sike Latency: A Likelihood Aroach Rick L. Jenison Deartment of Psychology University of Wisconsin Madison WI 53706 enison@wavelet.sych.wisc.edu ABSTRACT First-sike
More informationANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM
ANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM JOHN BINDER Abstract. In this aer, we rove Dirichlet s theorem that, given any air h, k with h, k) =, there are infinitely many rime numbers congruent to
More informationREVIEW PERIOD REORDER POINT PROBABILISTIC INVENTORY SYSTEM FOR DETERIORATING ITEMS WITH THE MIXTURE OF BACKORDERS AND LOST SALES
T. Vasudha WARRIER, PhD Nita H. SHAH, PhD Deartment of Mathematics, Gujarat University Ahmedabad - 380009, Gujarat, INDIA E-mail : nita_sha_h@rediffmail.com REVIEW PERIOD REORDER POINT PROBABILISTIC INVENTORY
More informationSwitching Control of Air-Fuel Ratio in Spark Ignition Engines
American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, FrB.4 Switching Control of Air-Fuel Ratio in Sar Ignition Engines Denis V. Efimov, Member, IEEE, Hosein Javaherian, Vladimir
More informationON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE
MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE
More informationSynoptic Meteorology I: The Geostrophic Approximation. 30 September, 7 October 2014
The Equations of Motion Synotic Meteorology I: The Geostrohic Aroimation 30 Setember, 7 October 2014 In their most general form, and resented without formal derivation, the equations of motion alicable
More informationUniform Sample Generations from Contractive Block Toeplitz Matrices
IEEE TRASACTIOS O AUTOMATIC COTROL, VOL 5, O 9, SEPTEMBER 6 559 Uniform Samle Generations from Contractive Bloc Toelitz Matrices Tong Zhou and Chao Feng Abstract This note deals with generating a series
More informationCENTRALIZED AND DECENTRALIZED ADAPTIVE FAULT-TOLERANT CONTROL APPLIED TO INTERCONNECTED AND NETWORKED CONTROL SYSTEM
CENRALIZED AND DECENRALIZED ADAPIVE FAUL-OLERAN CONROL APPLIED O INERCONNECED AND NEWORKED CONROL SYSEM Amina CHALLOUF (a), Adel ELLILI (b), Christohe AUBRUN (c), Mohamed Naceur ABDELKRIM (d) (a,d) Unité
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More information2.2 BUDGET-CONSTRAINED CHOICE WITH TWO COMMODITIES
Essential Miroeonomis -- 22 BUDGET-CONSTRAINED CHOICE WITH TWO COMMODITIES Continuity of demand 2 Inome effets 6 Quasi-linear, Cobb-Douglas and CES referenes 9 Eenditure funtion 4 Substitution effets and
More informationA Note on Guaranteed Sparse Recovery via l 1 -Minimization
A Note on Guaranteed Sarse Recovery via l -Minimization Simon Foucart, Université Pierre et Marie Curie Abstract It is roved that every s-sarse vector x C N can be recovered from the measurement vector
More informationŽ. Ž. Ž. 2 QUADRATIC AND INVERSE REGRESSIONS FOR WISHART DISTRIBUTIONS 1
The Annals of Statistics 1998, Vol. 6, No., 573595 QUADRATIC AND INVERSE REGRESSIONS FOR WISHART DISTRIBUTIONS 1 BY GERARD LETAC AND HELENE ` MASSAM Universite Paul Sabatier and York University If U and
More informationp-adic Measures and Bernoulli Numbers
-Adic Measures and Bernoulli Numbers Adam Bowers Introduction The constants B k in the Taylor series exansion t e t = t k B k k! k=0 are known as the Bernoulli numbers. The first few are,, 6, 0, 30, 0,
More informationMean Square Stability Analysis of Sampled-Data Supervisory Control Systems
17th IEEE International Conference on Control Alications Part of 28 IEEE Multi-conference on Systems and Control San Antonio, Texas, USA, Setember 3-5, 28 WeA21 Mean Square Stability Analysis of Samled-Data
More informationA PARTICLE SWARM OPTIMIZATION APPROACH FOR TUNING OF SISO PID CONTROL LOOPS NELENDRAN PILLAY
A PARTICLE SWARM OPTIMIZATION APPROACH FOR TUNING OF SISO PID CONTROL LOOPS NELENDRAN PILLAY 2008 A PARTICLE SWARM OPTIMIZATION APPROACH FOR TUNING OF SISO PID CONTROL LOOPS By Nelendran Pillay Student
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationVOLTAGE COLLAPSE AVOIDANCE IN POWER SYSTEMS: A RECEDING HORIZON APPROACH
Intelligent Automation and Soft Comuting, Vol. 12, No. X,. 1-14, 26 Coyright 26, SI Press Printed in the USA. All rights reserved VOLAGE COLLAPSE AVOIDANCE IN POWER SYSEMS: A RECEDING HORIZON APPROACH
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationBEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH
BEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH DORIN BUCUR, ALESSANDRO GIACOMINI, AND PAOLA TREBESCHI Abstract For Ω R N oen bounded and with a Lischitz boundary, and
More informationThe non-stochastic multi-armed bandit problem
Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at
More informationLocation of solutions for quasi-linear elliptic equations with general gradient dependence
Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations
More informationComplex Analysis Homework 1
Comlex Analysis Homework 1 Steve Clanton Sarah Crimi January 27, 2009 Problem Claim. If two integers can be exressed as the sum of two squares, then so can their roduct. Proof. Call the two squares that
More informationLower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data
Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment
More information