A Method of Setting the Penalization Constants in the Suboptimal Linear Quadratic Tracking Method

Size: px
Start display at page:

Download "A Method of Setting the Penalization Constants in the Suboptimal Linear Quadratic Tracking Method"

Transcription

1 XXVI. ASR '21 Seminar, Instruments and Control, Ostrava, Aril 26-27, 21 Paer 57 A Method of Setting the Penalization Constants in the Subotimal Linear Quadratic Tracking Method PERŮTKA, Karel Ing., Deartment of Control Theory, Institute of Information Technologies, Tomas Bata University in Zlín, nám. T.G. Masaryka 275, Zlín, , Czech Reublic kerutka@ft.utb.cz, htt://.utb.cz Abstract: In this aer, a simle solution of setting the enalization constants in the quadratic functional of the subotimal linear quadratic tracking method is resented. Subotimal linear quadratic method gives very good results in the adative control systems in comarison ith the dynamics inversion method and the ole lacement method and allos us to reduce the order of the controlled system. The enalization constants in this method had to be chosen manually. In this aer, one of the ossibilities, ho to set the constants of the quadratic functional, is given. It is based on a standard regression model. This roosed method has been verified by simulation in Matlab. The simulation results demonstrate the usefulness of the roosed method. Keyords: dynamics inversion method, enalization constants, ole lacement method, subotimal linear quadratic tracking 1 Introduction SISO (single inut, single outut) system is one of the most frequently used systems for control, esecially hen comuters drive the control rocess. It has an otimal behavior and a lot of methods for control have been introduced and verified at the SISO systems (Balátě 1996). When the standard controllers ith fixed, though adjustable, arameters are unsatisfactory, a self-tuning controller (STC) is an alternative for control (Bobál et al. 1999, Bobál et al. 2). The main idea on hich an STC is founded is a combination of a recursive identification rocedure and a choice of controller synthesis. The advantages of SISO system together ith the adative controller can be alied to the decentralized control. Perůtka (2) comared three techniques for the controller synthesis, the dynamics inversion method (Bobál et al. 2), the ole lacement method (Bobál et al. 1999) and the subotimal linear quadratic (LQ) tracking method (Dostál and Bobál 1999). All of them give good behavior, but only the subotimal LQ tracking method allos us good control of reduced order system (Perůtka 2).This method as introduced by Dostál and Bobál (1999). A comutation of all oles in every counting eriod and stabilizing controller for unsteady systems are the big advantages of this method. Hoever, the enalization constants in the quadratic functional have to be set manually. In this aer, one of the ossibilities, ho to adjust these constants is given. This ne method is based on a standard regression model and reflect the influence of overshooting of the outut signal. Firstly, the subotimal LQ tracking method is briefly given. This is folloed by the descrition of the roosed method. Finally, simulation result of the roosed method is shon

2 2 Subotimal LQ tracking method LQ tracking means the linear quadratic deterministic control and e reduce the roblem into the subotimal tracking. There are to ays of solving this roblem (Dostál and Bobál 1999). The first one is based on an estimation of state-sace vectors and has to solve the Riccati equations. Another ay uses the inut-outut descrition, the external linear model (ELM) of the controlled system. Sectral factorization and other comlex oerations ith the multinominals are alied in this method. It is requested a feedback controller Q as the controlled system ould be stable in the otimal tracking roblem. Thus, e ant to minimize the quadratic functional J 2 { () + 2 µ e t ϕu () t } dt (1) here µ, ϕ > are the enalization constants. This functional can be ritten in agreement ith the Parseval theorem as j * { e () s µ e() s + u () s ϕu() s } 1 * J 2πj j ds (2) The controller synthesis ith the otimal LQ tracking method gives a controller, hich is ractically nonroductive because of the deendence of the arameters on the initial values. Therefore, the controller synthesis method giving the subotimal tracking of the reference signal for a grou of reference signals ill be derived. Let us suose the Lalace form of the reference signal funtion as here h () s sf ( h ) deg( ) f () s () s deg, () f (3) (4) Figure 1 - Block scheme of the 1 DOF control Comensator Q(s) is an integrator in Figure 1. Lalace form of the signals in the controlled system is b o y y u +, u fu ou a a, q u e (5) here o y are initial conditions for y and o u are initial conditions for u. We otain equations - 2 -

3 e u d q d ( af ) bo u fo y ( af ) bo u fo y (6) (7) After subtituting s instead f and (3) into (6) and (7) e obtain e m f d (8) here u q m f d as + bq d m ah bf o sf o u y (9) (1) (11) We obtain a conroller via subotimal tracking method, as is folloed. We ill count stable multinominals g and n as a result of the setral factorizations * * * ( as) ϕas + b µ b g g n * n a * a (12) (13) Multinominals q, in the controller transfer function are given ith solving these diohantic equations g * q - v * as b * µn (14) g * v * b (as) * ϕn (15) here deg (v) < deg (g), suosing that f solving only the one olynomial equation does not divide a, e can reduce the roblem on as + bq gn (16) 3 A method of setting the enalization constants ϕ, µ of the functional J for the subotimal LQ tracking method Let us suose f(ϕ,µ) (17) here is a maximum value of the overshooting of the outut signal in the ste resonse of the hole system, ϕ, µ are the enalization constants of the functional J. Furthemore, let us suose the global minimun of the equation (17). This function deendence for the system given in the simulation results section is shon in the table and figures belo

4 Table 1 - Deendence of the overshooting on the enalization constants µ µ,4,2 4,2,5,3 5,3,6,4 6,4,7,6 7,6,8,7 8,7 1,9 1,9 2,17 2,17 3,22 3,22 4,23 4,23 1,34 1,34 1,47 1,47 1,61 1,61 ϕ 1 ϕ 1,4,3,2,1,12Ln(µ) +,99 r, µ Figure 2 - Deendence of the overshooting on the enalization constant µ for ϕ 1,4,3,2,1,12Ln(µ) -,137 r, µ Figure 3 - Deendence of the overshooting on the enalization constant µ for ϕ 1-4 -

5 As e can see from Table 1 and Figure 2 and 3, e can choose one enalization constant, e.g. ϕ 1 and the second one change. Intersection of the logarithm regression curve ith the x axis, hich is used for roortional systems, gives us a first estimation of the otimal constant. If the overshooting still remains e ill do the correction ith the equation () k,8µ ( k 1) µ, k 1,2 K here k is the correction ste. (18) 4 Simulation examle Simulation as realized on the roortional system ith transfer function G () s s 2 6,2 + 12,2s ,7364 (,11s + 1)(,4544s + 1) (t), y(t),.4 y(t) (t) time (s) Figure 4 - Subotimal LQ tracking ith the enalization constants set manually (t), y(t),.4 y(t) (t) time (s) Figure 5 - Subotimal LQ tracking ith the enalization constants set by the roosed method - 5 -

6 .4 1 (t), 2 (t) y 2 (t) y 1 (t) čas (s) Figure 6 - A comarison of the both methods 5 Discussion In this aer, a method for setting the enalization constants of the quadratic functional of the subotimal LQ tracking method is shon. It is clear from Table 1 and Figure 2 and 3, that is suficient to count only one enalization constant, hile the other one can be chosen, the regression equations in Figure 2 and 3 are the same. The main advantage of the roosed method is, that this method seed u the control (see Figure 6). Red curve in Figure 6, hich is the detail of the Figure 2 and 3, reresents the outut signal ith the roosed algorithm and the green one the standard subotimal LQ tracking algorithm. 6 Conclusions The subotimal LQ tracking method as introduced by Dostál and Bobál (1999) and gives very good behavior. The method roosed in this aer seed u the control of this method. Sectral factotization is an imortant asect for stabilization of the system. Another advantage is its accurancy hich as verified by simulation in Matlab Simulink on a SISO system. These characteristic demostrate the usefulness of the roosed method. References BALÁTĚ, J Selected arts of automatic control. Brno: VUT Brno, 1996 (in Czech). BOBÁL, V., J. BÖHM, J. PROKOP and J. FESSL Practical Asects of Self-Tuning Controllers: Algorithms and Imlementation. VUTIUM Press, Brno: University of Technology, 1999 (in Czech). BOBÁL, V., DOSTÁL, P., MACHÁČEK, J., VÍTEČKOVÁ, M. 2. Self-Tuning PID Controllers Based on Dynamics Inversion Method. In: Proc. IFAC Worksho on Digital Control, University of Catalonia, Sain, 2 DOSTÁL, P., BOBÁL, V The subotimal tracking roblem in linear systems. In: 7 th Conference on Control and Automation, Haifa, Israel, 1999, PERŮTKA, K. 2. The Decentralized adative continuous-time control. Diloma ork. Zlín: VUT Brno, FT in Zlín, 2 (in Czech)

Robust Performance Design of PID Controllers with Inverse Multiplicative Uncertainty

Robust 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 information

Multivariable Generalized Predictive Scheme for Gas Turbine Control in Combined Cycle Power Plant

Multivariable Generalized Predictive Scheme for Gas Turbine Control in Combined Cycle Power Plant Multivariable Generalized Predictive Scheme for Gas urbine Control in Combined Cycle Power Plant L.X.Niu and X.J.Liu Deartment of Automation North China Electric Power University Beiing, China, 006 e-mail

More information

DEVELOPMENT AND VALIDATION OF A VERSATILE METHOD FOR THE CALCULATION OF HEAT TRANSFER IN WATER-BASED RADIANT SYSTEMS METHODS

DEVELOPMENT AND VALIDATION OF A VERSATILE METHOD FOR THE CALCULATION OF HEAT TRANSFER IN WATER-BASED RADIANT SYSTEMS METHODS Eleventh International IBPSA Conference Glasgo, Scotland July 27-30, 2009 DEVELOPMENT AND VALIDATION OF A VERSATILE METHOD FOR THE CALCULATION OF HEAT TRANSFER IN WATER-BASED RADIANT SYSTEMS Massimiliano

More information

Distributed Rule-Based Inference in the Presence of Redundant Information

Distributed 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 information

Feedback-error control

Feedback-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 information

On Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.**

On 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 information

ME 375 System Modeling and Analysis. Homework 11 Solution. Out: 18 November 2011 Due: 30 November 2011 = + +

ME 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 information

Solved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.

Solved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points. Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the

More information

Robust Beamforming via Matrix Completion

Robust Beamforming via Matrix Completion Robust Beamforming via Matrix Comletion Shunqiao Sun and Athina P. Petroulu Deartment of Electrical and Comuter Engineering Rutgers, the State University of Ne Jersey Piscataay, Ne Jersey 8854-858 Email:

More information

Roots Blower with Gradually Expanding Outlet Gap: Mathematical Modelling and Performance Simulation Yingjie Cai 1, a, Ligang Yao 2, b

Roots Blower with Gradually Expanding Outlet Gap: Mathematical Modelling and Performance Simulation Yingjie Cai 1, a, Ligang Yao 2, b 2nd International Conference on Advances in Mechanical Engineering and Industrial Informatics (AMEII 216) Roots Bloer ith Gradually Exanding Outlet Ga: Mathematical Modelling and Performance Simulation

More information

Finding Shortest Hamiltonian Path is in P. Abstract

Finding Shortest Hamiltonian Path is in P. Abstract Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs

More information

MODELLING, SIMULATION AND ROBUST ANALYSIS OF THE TEMPERATURE PROCESS CONTROL

MODELLING, 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 information

ADAPTIVE PID CONTROLLER WITH ON LINE IDENTIFICATION

ADAPTIVE PID CONTROLLER WITH ON LINE IDENTIFICATION Journal of ELECTRICAL ENGINEERING, VOL. 53, NO. 9-10, 00, 33 40 ADAPTIVE PID CONTROLLER WITH ON LINE IDENTIFICATION Jiří Macháče Vladimír Bobál A digital adaptive PID controller algorithm which contains

More information

THE 3-DOF helicopter system is a benchmark laboratory

THE 3-DOF helicopter system is a benchmark laboratory Vol:8, No:8, 14 LQR Based PID Controller Design for 3-DOF Helicoter System Santosh Kr. Choudhary International Science Index, Electrical and Information Engineering Vol:8, No:8, 14 waset.org/publication/9999411

More information

Research on Evaluation Method of Organization s Performance Based on Comparative Advantage Characteristics

Research on Evaluation Method of Organization s Performance Based on Comparative Advantage Characteristics Vol.1, No.10, Ar 01,.67-7 Research on Evaluation Method of Organization s Performance Based on Comarative Advantage Characteristics WEN Xin 1, JIA Jianfeng and ZHAO Xi nan 3 Abstract It as under the guidance

More information

PREDICTIVE CONTROL OF TIME-DELAY PROCESSES

PREDICTIVE CONTROL OF TIME-DELAY PROCESSES PREDICIVE CONROL OF IME-DELAY PROCESSES Marek Kubalčík, Vladimír Bobál omas Bata Universit in Zlín Facult of Alied Informatics Nad Stráněmí 45 76 5 Zlín Czech Reublic E-mail: kubalcikl@fai.utb.cz KEYWORDS

More information

y p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2

y p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2 ROBUSTNESS OF FLEXIBLE SYSTEMS WITH COMPONENT-LEVEL UNCERTAINTIES Peiman G. Maghami Λ NASA Langley Research Center, Hamton, VA 368 Robustness of flexible systems in the resence of model uncertainties at

More information

Frequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer

Frequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer Frequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer C.P. an + F. Crusca # M. Aldeen * + School of Engineering, Monash University Malaysia, 2 Jalan Kolej, Bandar Sunway, 4650 Petaling,

More information

Recursive Estimation of the Preisach Density function for a Smart Actuator

Recursive Estimation of the Preisach Density function for a Smart Actuator Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator

More information

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

MODELING 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 information

Robust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System

Robust 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 information

Algebraic Parameter Estimation of Damped Exponentials

Algebraic Parameter Estimation of Damped Exponentials Algebraic Parameter Estimation of Damed Exonentials Aline Neves, Maria Miranda, Mamadou Mbou To cite this version: Aline Neves, Maria Miranda, Mamadou Mbou Algebraic Parameter Estimation of Damed Exonentials

More information

Approximating min-max k-clustering

Approximating 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 information

Damage Identification from Power Spectrum Density Transmissibility

Damage Identification from Power Spectrum Density Transmissibility 6th Euroean Worksho on Structural Health Monitoring - h.3.d.3 More info about this article: htt://www.ndt.net/?id=14083 Damage Identification from Power Sectrum Density ransmissibility Y. ZHOU, R. PERERA

More information

Controllability and Resiliency Analysis in Heat Exchanger Networks

Controllability 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 information

Keywords: Vocal Tract; Lattice model; Reflection coefficients; Linear Prediction; Levinson algorithm.

Keywords: Vocal Tract; Lattice model; Reflection coefficients; Linear Prediction; Levinson algorithm. Volume 3, Issue 6, June 213 ISSN: 2277 128X International Journal of Advanced Research in Comuter Science and Software Engineering Research Paer Available online at: www.ijarcsse.com Lattice Filter Model

More information

On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm

On 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 information

Multivariable PID Control Design For Wastewater Systems

Multivariable 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 information

Feedback-Based Iterative Learning Control for MIMO LTI Systems

Feedback-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 information

Optimal Choice of Weighting Factors in Adaptive Linear Quadratic Control

Optimal Choice of Weighting Factors in Adaptive Linear Quadratic Control International Journal of Automation and Computing 11(3), June 2014, 241-248 DOI: 10.1007/s11633-014-0786-5 Optimal Choice of Weighting Factors in Adaptive Linear Quadratic Control Jiri Vojtesek Petr Dostal

More information

Passive Identification is Non Stationary Objects With Closed Loop Control

Passive Identification is Non Stationary Objects With Closed Loop Control IOP Conerence Series: Materials Science and Engineering PAPER OPEN ACCESS Passive Identiication is Non Stationary Obects With Closed Loo Control To cite this article: Valeriy F Dyadik et al 2016 IOP Con.

More information

Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier

Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier Australian Journal of Basic and Alied Sciences, 5(12): 2010-2020, 2011 ISSN 1991-8178 Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doed Fiber Amlifier

More information

Figure : An 8 bridge design grid. (a) Run this model using LOQO. What is the otimal comliance? What is the running time?

Figure : An 8 bridge design grid. (a) Run this model using LOQO. What is the otimal comliance? What is the running time? 5.094/SMA53 Systems Otimization: Models and Comutation Assignment 5 (00 o i n ts) Due Aril 7, 004 Some Convex Analysis (0 o i n ts) (a) Given ositive scalars L and E, consider the following set in three-dimensional

More information

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model

Shadow 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

State Estimation with ARMarkov Models

State 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 information

Computer arithmetic. Intensive Computation. Annalisa Massini 2017/2018

Computer arithmetic. Intensive Computation. Annalisa Massini 2017/2018 Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3

More information

Oil Temperature Control System PID Controller Algorithm Analysis Research on Sliding Gear Reducer

Oil 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 information

Design of NARMA L-2 Control of Nonlinear Inverted Pendulum

Design 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 information

A Simple Fuzzy PI Control of Dual-Motor Driving Servo System

A Simple Fuzzy PI Control of Dual-Motor Driving Servo System MATEC Web of Conferences 04, 0006 (07) DOI: 0.05/ matecconf/07040006 IC4M & ICDES 07 A Simle Fuzzy PI Control of Dual-Motor Driving Servo System Haibo Zhao,,a, Chengguang Wang 3 Engineering Technology

More information

STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS

STABILITY 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 information

Indirect Rotor Field Orientation Vector Control for Induction Motor Drives in the Absence of Current Sensors

Indirect Rotor Field Orientation Vector Control for Induction Motor Drives in the Absence of Current Sensors Indirect Rotor Field Orientation Vector Control for Induction Motor Drives in the Absence of Current Sensors Z. S. WANG *, S. L. HO ** * College of Electrical Engineering, Zhejiang University, Hangzhou

More information

4 Scheduling. Outline of the chapter. 4.1 Preliminaries

4 Scheduling. Outline of the chapter. 4.1 Preliminaries 4 Scheduling In this section, e consider so-called Scheduling roblems I.e., if there are altogether M machines or resources for each machine, a roduction sequence of all N jobs has to be found as ell as

More information

Linear diophantine equations for discrete tomography

Linear diophantine equations for discrete tomography Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,

More information

ALTERNATIVE SOLUTION TO THE QUARTIC EQUATION by Farid A. Chouery 1, P.E. 2006, All rights reserved

ALTERNATIVE SOLUTION TO THE QUARTIC EQUATION by Farid A. Chouery 1, P.E. 2006, All rights reserved ALTERNATIVE SOLUTION TO THE QUARTIC EQUATION b Farid A. Chouer, P.E. 006, All rights reserved Abstract A new method to obtain a closed form solution of the fourth order olnomial equation is roosed in this

More information

LQ CONTROL OF HEAT EXCHANGER DESIGN AND SIMULATION

LQ CONTROL OF HEAT EXCHANGER DESIGN AND SIMULATION LQ CONTROL OF HEAT EXCHANGER DESIGN AND SIMULATION Vladimír Bobál,, Petr Dostál,, Marek Kubalčík and Stanislav Talaš Tomas Bata University in Zlín Centre of Polymer Systems, University Institute Department

More information

A Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split

A Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com

More information

Modeling and Model Predictive Control of Nonlinear Hydraulic System

Modeling and Model Predictive Control of Nonlinear Hydraulic System Modeling and Model Predictive Control of Nonlinear Hydraulic System Petr Chalupa, Jakub Novák Department of Process Control, Faculty of Applied Informatics, Tomas Bata University in Zlin, nám. T. G. Masaryka

More information

LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL

LINEAR 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 information

Improving AOR Method for a Class of Two-by-Two Linear Systems

Improving AOR Method for a Class of Two-by-Two Linear Systems Alied Mathematics 2 2 236-24 doi:4236/am22226 Published Online February 2 (htt://scirporg/journal/am) Imroving AOR Method for a Class of To-by-To Linear Systems Abstract Cuixia Li Shiliang Wu 2 College

More information

4. Score normalization technical details We now discuss the technical details of the score normalization method.

4. 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 information

Spectral Analysis by Stationary Time Series Modeling

Spectral Analysis by Stationary Time Series Modeling Chater 6 Sectral Analysis by Stationary Time Series Modeling Choosing a arametric model among all the existing models is by itself a difficult roblem. Generally, this is a riori information about the signal

More information

MODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS

MODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A

More information

Simulation of the PIR detector active function

Simulation of the PIR detector active function 04036 (016) DOI: 10.1051/ matecconf/0167604036 CSCC 016 Simulation of the PIR detector active function Rudolf Drga 1,a, Dagmar Janacova and Hana Charvatova 3 1 Tomas Bata University in Zlín, Faculty of

More information

Observer/Kalman Filter Time Varying System Identification

Observer/Kalman Filter Time Varying System Identification Observer/Kalman Filter Time Varying System Identification Manoranjan Majji Texas A&M University, College Station, Texas, USA Jer-Nan Juang 2 National Cheng Kung University, Tainan, Taiwan and John L. Junins

More information

Research of PMU Optimal Placement in Power Systems

Research of PMU Optimal Placement in Power Systems Proceedings of the 5th WSEAS/IASME Int. Conf. on SYSTEMS THEORY and SCIENTIFIC COMPUTATION, Malta, Setember 15-17, 2005 (38-43) Research of PMU Otimal Placement in Power Systems TIAN-TIAN CAI, QIAN AI

More information

The Graph Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule

The Graph Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule The Grah Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule STEFAN D. BRUDA Deartment of Comuter Science Bisho s University Lennoxville, Quebec J1M 1Z7 CANADA bruda@cs.ubishos.ca

More information

Linear Matrix Inequalities Approach to Input Covariance Constraint Control With Application to Electronic Throttle

Linear Matrix Inequalities Approach to Input Covariance Constraint Control With Application to Electronic Throttle Ali Khudhair Al-Jiboory East Lansing, MI 88 Andrew White East Lansing, MI 88 Shueng Zhang East Lansing, MI 88 Guoming Zhu Deartment of Electrical and Comuter Engineering, East Lansing, MI 88 Jongeun Choi

More information

Dynamic System Eigenvalue Extraction using a Linear Echo State Network for Small-Signal Stability Analysis a Novel Application

Dynamic System Eigenvalue Extraction using a Linear Echo State Network for Small-Signal Stability Analysis a Novel Application Dynamic System Eigenvalue Extraction using a Linear Echo State Network for Small-Signal Stability Analysis a Novel Alication Jiaqi Liang, Jing Dai, Ganesh K. Venayagamoorthy, and Ronald G. Harley Abstract

More information

Impact Damage Detection in Composites using Nonlinear Vibro-Acoustic Wave Modulations and Cointegration Analysis

Impact Damage Detection in Composites using Nonlinear Vibro-Acoustic Wave Modulations and Cointegration Analysis 11th Euroean Conference on Non-Destructive esting (ECND 214), October 6-1, 214, Prague, Czech Reublic More Info at Oen Access Database www.ndt.net/?id=16448 Imact Damage Detection in Comosites using Nonlinear

More information

Simulation of 3-Phase 2- Stator Induction Motor Using MATLAB Platform

Simulation of 3-Phase 2- Stator Induction Motor Using MATLAB Platform International Journal of Alied Engineering Research ISSN 0973-456 Volume 3, Number (08). 9437-944 Simulation of 3-Phase - Stator Induction Motor Using MATLAB Platform Pallavi R.Burande Deartment of Electrical

More information

Using a Computational Intelligence Hybrid Approach to Recognize the Faults of Variance Shifts for a Manufacturing Process

Using a Computational Intelligence Hybrid Approach to Recognize the Faults of Variance Shifts for a Manufacturing Process Journal of Industrial and Intelligent Information Vol. 4, No. 2, March 26 Using a Comutational Intelligence Hybrid Aroach to Recognize the Faults of Variance hifts for a Manufacturing Process Yuehjen E.

More information

Symbolic Dynamical Model of Average Queue Size of Random Early Detection Algorithm

Symbolic Dynamical Model of Average Queue Size of Random Early Detection Algorithm Symbolic Dynamical odel of Average Queue Size of Rom Early Detection Algorithm Charlotte Yuk-Fan Ho Telehone: 44 () 788 5555 ext. 4333 Fax: 44 () 788 7997 Email: c.ho@mul.ac.uk School of athematical Sciences

More information

Determination of the Best Apodization Function and Grating Length of Linearly Chirped Fiber Bragg Grating for Dispersion Compensation

Determination of the Best Apodization Function and Grating Length of Linearly Chirped Fiber Bragg Grating for Dispersion Compensation 84 JOURNAL OF COMMUNICATIONS, VOL. 7, NO., NOVEMBER Determination of the Best Aodization Function and Grating Length of Linearly Chired Fiber Bragg Grating for Disersion Comensation Sher Shermin A. Khan

More information

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations

Preconditioning 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 information

MODELING AND KINETICS OF PROTEIN ENZYMATIC HYDROLYSIS FROM AMARANTH FLOUR AND WHOLE GRAIN

MODELING AND KINETICS OF PROTEIN ENZYMATIC HYDROLYSIS FROM AMARANTH FLOUR AND WHOLE GRAIN 1 th International Research/Exert Conference Trends in the Develoment of Machinery and Associated Technology TMT 008, Istanbul, Turkey, 6-30 August, 008 MODELING AND KINETICS OF PROTEIN ENZYMATIC HYDROLYSIS

More information

Statistical Methods for River Runoff Prediction

Statistical Methods for River Runoff Prediction Water Resources, Vol. 32, No. 2, 25,. 115 126. Translated from Vodnye Resursy, Vol. 32, No. 2, 25,. 133 145. Original Russian Text Coyright 25 by Pisarenko, Lyubushin, Bolgov, Rukavishnikova, Kanyu, Kanevskii,

More information

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test)

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test) Chater 225 Tests for Two Proortions in a Stratified Design (Cochran/Mantel-Haenszel Test) Introduction In a stratified design, the subects are selected from two or more strata which are formed from imortant

More information

Detection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform

Detection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform Detection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform Chaoquan Chen and Guoing Qiu School of Comuter Science and IT Jubilee Camus, University of Nottingham Nottingham

More information

Genetic Algorithm Based PID Optimization in Batch Process Control

Genetic Algorithm Based PID Optimization in Batch Process Control International Conference on Comuter Alications and Industrial Electronics (ICCAIE ) Genetic Algorithm Based PID Otimization in Batch Process Control.K. Tan Y.K. Chin H.J. Tham K.T.K. Teo odelling, Simulation

More information

A New Method of DDB Logical Structure Synthesis Using Distributed Tabu Search

A New Method of DDB Logical Structure Synthesis Using Distributed Tabu Search A New Method of DDB Logical Structure Synthesis Using Distributed Tabu Search Eduard Babkin and Margarita Karunina 2, National Research University Higher School of Economics Det of nformation Systems and

More information

COMPARISON OF VARIOUS OPTIMIZATION TECHNIQUES FOR DESIGN FIR DIGITAL FILTERS

COMPARISON OF VARIOUS OPTIMIZATION TECHNIQUES FOR DESIGN FIR DIGITAL FILTERS NCCI 1 -National Conference on Comutational Instrumentation CSIO Chandigarh, INDIA, 19- March 1 COMPARISON OF VARIOUS OPIMIZAION ECHNIQUES FOR DESIGN FIR DIGIAL FILERS Amanjeet Panghal 1, Nitin Mittal,Devender

More information

CONSTRUCTIVE NEURAL NETWORKS IN FORECASTING WEEKLY RIVER FLOW

CONSTRUCTIVE NEURAL NETWORKS IN FORECASTING WEEKLY RIVER FLOW CONSTRUCTIVE NEURAL NETWORKS IN FORECASTING WEEKLY RIVER FLOW MÊUSER VALENÇA Universidade Salgado de Oliveira - UNIVERSO /Comanhia Hidroelétrica do São Francisco - Chesf Rua Gasar Peres, 47/4 CDU, Recife,

More information

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO) Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment

More information

Analysis of Fractional order PID controller for Ceramic Infrared Heater

Analysis 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 information

Position Control of Induction Motors by Exact Feedback Linearization *

Position 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 information

Factor Analysis of Convective Heat Transfer for a Horizontal Tube in the Turbulent Flow Region Using Artificial Neural Network

Factor Analysis of Convective Heat Transfer for a Horizontal Tube in the Turbulent Flow Region Using Artificial Neural Network COMPUTATIONAL METHODS IN ENGINEERING AND SCIENCE EPMESC X, Aug. -3, 6, Sanya, Hainan,China 6 Tsinghua University ess & Sringer-Verlag Factor Analysis of Convective Heat Transfer for a Horizontal Tube in

More information

Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation

Uniformly 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 information

F. Augustin, P. Rentrop Technische Universität München, Centre for Mathematical Sciences

F. Augustin, P. Rentrop Technische Universität München, Centre for Mathematical Sciences NUMERICS OF THE VAN-DER-POL EQUATION WITH RANDOM PARAMETER F. Augustin, P. Rentro Technische Universität München, Centre for Mathematical Sciences Abstract. In this article, the roblem of long-term integration

More information

Journal of Chemical and Pharmaceutical Research, 2014, 6(5): Research Article

Journal of Chemical and Pharmaceutical Research, 2014, 6(5): Research Article Available online www.jocr.com Journal of Chemical and harmaceutical Research, 04, 6(5):904-909 Research Article ISSN : 0975-7384 CODEN(USA) : JCRC5 Robot soccer match location rediction and the alied research

More information

Estimation of the large covariance matrix with two-step monotone missing data

Estimation 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 information

On Optimization of Power Coefficient of HAWT

On Optimization of Power Coefficient of HAWT Journal of Power and Energy Engineering, 14,, 198- Published Online Aril 14 in Scies htt://wwwscirorg/journal/jee htt://dxdoiorg/1436/jee1448 On Otimization of Power Coefficient of HAWT Marat Z Dosaev

More information

The Recursive Fitting of Multivariate. Complex Subset ARX Models

The Recursive Fitting of Multivariate. Complex Subset ARX Models lied Mathematical Sciences, Vol. 1, 2007, no. 23, 1129-1143 The Recursive Fitting of Multivariate Comlex Subset RX Models Jack Penm School of Finance and lied Statistics NU College of Business & conomics

More information

arxiv: v1 [quant-ph] 20 Jun 2017

arxiv: 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 information

Open Access Design and Optimization of Feedforward-PID Integrated Controller for Pre-Injection Numerical Control Air-Powered Engine

Open Access Design and Optimization of Feedforward-PID Integrated Controller for Pre-Injection Numerical Control Air-Powered Engine Send Orders for Rerints to rerints@benthamscience.ae 56 The Oen Automation and Control Systems Journal, 5, 7, 56-6 Oen Access Design and Otimization of Feedforard-PID Integrated Controller for Pre-Injection

More information

Research of power plant parameter based on the Principal Component Analysis method

Research of power plant parameter based on the Principal Component Analysis method Research of ower lant arameter based on the Princial Comonent Analysis method Yang Yang *a, Di Zhang b a b School of Engineering, Bohai University, Liaoning Jinzhou, 3; Liaoning Datang international Jinzhou

More information

The Noise Power Ratio - Theory and ADC Testing

The Noise Power Ratio - Theory and ADC Testing The Noise Power Ratio - Theory and ADC Testing FH Irons, KJ Riley, and DM Hummels Abstract This aer develos theory behind the noise ower ratio (NPR) testing of ADCs. A mid-riser formulation is used for

More information

1. Introduction. 2. Background of elliptic curve group. Identity-based Digital Signature Scheme Without Bilinear Pairings

1. Introduction. 2. Background of elliptic curve group. Identity-based Digital Signature Scheme Without Bilinear Pairings Identity-based Digital Signature Scheme Without Bilinear Pairings He Debiao, Chen Jianhua, Hu Jin School of Mathematics Statistics, Wuhan niversity, Wuhan, Hubei, China, 43007 Abstract: Many identity-based

More information

Generalized Coiflets: A New Family of Orthonormal Wavelets

Generalized Coiflets: A New Family of Orthonormal Wavelets Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University

More information

Uncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning

Uncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment

More information

Wolfgang POESSNECKER and Ulrich GROSS*

Wolfgang 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 information

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MATHEMATICAL 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 information

Modeling of Hydraulic Control Valves

Modeling of Hydraulic Control Valves Modeling of Hydraulic Control Valves PETR CHALUPA, JAKUB NOVAK, VLADIMIR BOBAL Department of Process Control, Faculty of Applied Informatics Tomas Bata University in Zlin nam. T.G. Masaryka 5555, 76 1

More information

2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.

2-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 information

Pulse Propagation in Optical Fibers using the Moment Method

Pulse Propagation in Optical Fibers using the Moment Method Pulse Proagation in Otical Fibers using the Moment Method Bruno Miguel Viçoso Gonçalves das Mercês, Instituto Suerior Técnico Abstract The scoe of this aer is to use the semianalytic technique of the Moment

More information

Cybernetic Interpretation of the Riemann Zeta Function

Cybernetic Interpretation of the Riemann Zeta Function Cybernetic Interretation of the Riemann Zeta Function Petr Klán, Det. of System Analysis, University of Economics in Prague, Czech Reublic, etr.klan@vse.cz arxiv:602.05507v [cs.sy] 2 Feb 206 Abstract:

More information

EE451/551: Digital Control. Relationship Between s and z Planes. The Relationship Between s and z Planes 11/10/2011

EE451/551: Digital Control. Relationship Between s and z Planes. The Relationship Between s and z Planes 11/10/2011 /0/0 EE45/55: Digital Control Chater 6: Digital Control System Design he Relationshi Between s and Planes As noted reviously: s j e e e e r s j where r e and If an analog system has oles at: s n jn a jd

More information

Generation of Linear Models using Simulation Results

Generation of Linear Models using Simulation Results 4. IMACS-Symosium MATHMOD, Wien, 5..003,. 436-443 Generation of Linear Models using Simulation Results Georg Otte, Sven Reitz, Joachim Haase Fraunhofer Institute for Integrated Circuits, Branch Lab Design

More information

Adaptive Control of Fluid Inside CSTR Using Continuous-Time and Discrete-Time Identification Model and Different Control Configurations

Adaptive Control of Fluid Inside CSTR Using Continuous-Time and Discrete-Time Identification Model and Different Control Configurations Adaptive Control of Fluid Inside CSTR Using Continuous-Time and Discrete-Time Identification Model and Different Control Configurations JIRI VOJTESEK, PETR DOSTAL Tomas Bata University in Zlin Faculty

More information

Mathematical Efficiency Modeling of Static Power Converters

Mathematical Efficiency Modeling of Static Power Converters Fabrício Hoff Duont Regional Integrated University of Uer Uruguai and Missions (URI Av. Assis Brasil, 9, 980 000 Frederico Westhalen, RS, Brazil Contact: fhd@ieee.org Mathematical Efficiency Modeling of

More information

A 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 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 information