Open access peer-reviewed chapter

Use of Fuzzy Logic for Design and Control of Nonlinear MIMO Systems

By Pavol Fedor and Daniela Perduková

Submitted: October 28th 2016Reviewed: February 22nd 2017Published: August 30th 2017

DOI: 10.5772/68050

Downloaded: 351


Standard analytical methods are often ineffective or even useless for design of nonlinear control systems with imprecisely known parameters. The use of fuzzy logic principles presents one possible way to control such systems which can be used both for modeling and design of the control. The advantage of using this method consists in its simplicity and easy way of developing the algorithm, which in the phase of designing the controllers and also for modeling the features of the designed structures, allows the use of computer technology. Simplicity of the proposed structure (usually with the PI controllers) and determination of their parameters without any need for complex mathematical description present another considerable advantage of the used method. This chapter presents two typical examples of designing the control of nonlinear multi‐input multi‐output (MIMO) systems from the field of mechatronic systems based on fuzzy logic principles.


  • fuzzy PI torque controller
  • asynchronous motor
  • continuous line
  • black‐box modeling
  • fuzzy model‐based control

1. Introduction

There are many processes in technological practice the analytical description of which is rather complicated. This can be due to their complexity, nonlinearity, transfer lags, complicated measurement of important parameters, etc. However, information on the performance of these processes can often be obtained experimentally (by suitably chosen measurements or by monitoring their responses to the control activities of the operator). In these situations, fuzzy systems can always be considered as an alternative for system modeling and control.

Applying fuzzy logic in fuzzy controller design is very often a suitable possibility of solving problem issues in control in various fields of industry because these controllers are an effective tool for achieving high‐quality properties of the controlled systems [19]. The disadvantage in this case is the unsystematic approach to their synthesis and a relatively demanding analysis of their stability. A fuzzy controller design is primarily based on the fuzzification of its range of inputs and the setting up of rules of its behavior within this range. The behavior of classic fuzzy controllers was designed on basis of linguistic rules obtained from experts. However, this knowledge is not always easily obtained, especially in cases of higher‐order nonlinear systems [10, 11]. For this reason, special attention has been focused in recent years on the design of fuzzy control systems that are not based on the search for expert linguistic rules [1217]. Control methods based on the controlled system fuzzy model have many modifications that depend on the particular application [6, 1317], while the quality of the fuzzy model of the controlled system is also of significant importance.

It has been proved that fuzzy modeling can be recognized as one of the nonlinear black‐box modeling techniques [11, 12, 1820]. When designing a black‐box fuzzy system, it is necessary to identify its qualitative properties only on the basis of experimentally measured data, while neither its structure nor its parameters are known. That often results in problems with inconsistency of the database, problems with covering the entire space of possible inputs, etc. [2022], which makes the fuzzy model unusable in practical applications. In the design of a black‐box fuzzy model of a dynamic system, a suitable method for the selection of qualitative properties from the collected database always needs to be applied. The functional dependencies between inputs and outputs can then be used for developing a suitable nonparametric fuzzy model of the process that can be applied in the design of their control [2327].

Two typical examples for designing the control of nonlinear MIMO systems using fuzzy approach are presented in this chapter:

  • The design of fuzzy PI torque controller of the PI type for a drive with an induction motor, whose parameters and rules are obtained by searching control input of such a vector which is optimal in terms of the selected criterion of optimality.

  • The design of control for middle part of a continuous line for material processing by tension, where the continuous line presents a nonlinear MIMO system. Its control requires to ensure decoupled control of individual subsystems, because the output quality of the processed material depends directly on quality of the control. The controllers of the subsystems ensuring such decoupling usually are of complex structures and, when designing them by analytic methods, they are often unrealizable. When the continuous line is presented by a fuzzy model, it is possible to design simple controllers of the PI type ensuring high‐quality dynamical properties of the controlled system.

2. Design of fuzzy torque controller for asynchronous motor drive

An asynchronous motor represents a strongly nonlinear fifth‐order system, whose good quality vector torque control is solved by relatively complex mathematical transformations and leads to a complicated control structure [28, 29]. Therefore, a fuzzy system for design of torque controller for asynchronous motor drive has been used. The fuzzy controller design is based on the concept of identifying the time sequence of the input signal into the controlled system that will provide the control target in terms of the selected optimality criterion. Fuzzy controller design method is characterized by simplicity, and quality of control is appropriate to the considered drive.

API‐type discrete controller is generally described by the equation:


where uk and uk−1 are values of controller output in the relevant sampling steps, ek and ek−1 are values of the control error, and q0 and q1 are parameters of controller [28]. From these follows, it is possible structure, shown in Figure 1.

Figure 1.

Block diagram of fuzzy controller.

A discrete fuzzy PI controller can be described, for example, by the following rules:


where quantities ek, ek−1, and duk are fuzzy variables that describe the relevant workspace of the fuzzy PI controller. The fuzzy controller design procedure consists of the following three steps:

  1. Step 1. Finding the optimal sequence of input values

Fuzzy rules and fuzzification of the fuzzy PI controller workspace can be identified by means of relations expressed by triplets [ekek−1duk] that have been obtained for its optimal behavior. This behavior is represented by the time sequence of input vector duopt, at which the value of the optimality criterion is minimal. It is suitable to choose this criterion in the next form of the integral of the quadratic deviation of the system output y from the desired value w:


When obtaining the sequence of values of vector duopt, we apply to the drive input various sequences of input vector and evaluate the criterion value Eq. (3). This is a standard optimization task that can be solved, for example, by suitable geometric division of the drive workspace (which leads to rather high computing demands in the design process), or by applying the genetic algorithm method, which significantly speeds up the whole process of identification.

  1. Step 2. Finding the database of optimal data

Having found the optimal input sequence duopt for the controlled drive, we then set up the database of triplets [ekek−1duk], which describes relations between the inputs and the outputs of the optimal fuzzy PI controller.

  1. Step 3. Designing the fuzzy controller from the optimal data database

From the obtained triplets [ekek−1duk] of optimal data, it is possible to design a concrete fuzzy controller of various types using standard procedures of clustering the data into significant clusters and describing them by means of rules.

In the concrete application of the said procedure in a drive with asynchronous motor, we will use its analytical model (see Refs. [28, 29]). If we consider a rotating system which rotates with the frequency of the motor’s stator field (usually marked by coordinates x, y), the mathematical model of the device is described by the following equations:


Used symbols:

isx, isy components of stator current space vector is

irx, iry components of rotor current space vector ir

usx, usy components of stator voltage space vector us

ωm the motor mechanical angular speed

ω1 angular frequency of the stator voltage

p number of pole pairs (p = 2)

ω2 slip angular speed ω2 = ω1ωm

Rs, Rr stator and rotor phase resistance

Ψs, Ψr stator and rotor magnetic flux

Me electrical motor moment

AC drive parameters are given in the Appendix.

As a standard, asynchronous motors are supplied from static voltage frequency converters in which the stator frequency and voltage rate are U1/ω1.

When connected directly to the power supply network, the motor shows a large increase in torque and also in current (Figure 2).

Figure 2.

Direct connection of drive with AM.

Let the aim of the torque controller design be to adjust the slip ω2 (i.e., the difference of ω1ωm) according to the desired torque value. The structure of the controlled system will be as shown in Figure 3.

Figure 3.

Structure of the controlled drive.

For finding the optimal input sequence duopt, we shall use the diagram shown in Figure 4.

Figure 4.

Block diagram for finding vector duopt.

We search for the input signal sequence duopt through such changes of drive input voltage and frequency that will lead to minimal value of the criterion according to Eq. (3). This goal can be achieved, for example, by application of the genetic algorithm method, which efficiently enables finding the extreme of the selected function in a given space of mutations. Having selected input signal sampling time 50 ms and desired value of motor torque Mz = 30, we obtained the optimal input vector [0.165 0.145 0.155 0.12 0.105 0.095 0 0 0 0 0 0], as illustrated in Figure 5.

Figure 5.

Optimal start‐up of drive with AM.

The said optimal input signal sequence was used for generating the database of triplets [ωmkωmk−1 (ω1ωm)k] for the design of the P‐type fuzzy controller. Using the Anfisedit tool in the Matlab program, a static Sugeno‐type fuzzy system was designed, based on the acquired database of measured data ([ωmkωmk−1 (ω1ωm)k]). This system has three zones for input data fuzzification, and its internal structure can be seen in Figure 6.

Figure 6.

Structure of AM fuzzy controller.

The resulting torque control structure of AM with the designed fuzzy controller is shown in Figure 7.

Figure 7.

Structure of torque control in AM with fuzzy controller.

The start‐up of the drive with desired torque Mz = 30 Nm using the designed fuzzy controller is shown in Figure 8.

Figure 8.

Optimal start‐up of drive with AM with fuzzy controller.

The comparison of Figures 5 and 8 shows that the application of the fuzzy controller resulted in the achievement of the desired torque responses of the AM at start‐up.

2.1. Discussion

The fuzzy controller designed in this chapter is a PI controller, and it provides optimal dynamics in terms of the selected criterion Eq. (3). The design procedure consists of three steps:

In the first step, we search for such input sequence of input vector u of the asynchronous motor that would provide minimal value of function Eq. (3). This can be achieved by modeling various vector u input sequences in the system with asynchronous motor, and by a simulation of its behavior for a given sequence. The selection of sampling time for a change of the input vector value and by that also of the number of its samples has to be defined according to the Shannon‐Kotelnikov theorem. Obviously, with an increase in the number of samples also the time required for simulation and for optimal vector uopt identification grows. As the whole process of optimal vector identification can be fully automated and the performance of computing means nowadays is sufficient, there is no essential problem to identify the optimal control input signal for a concrete drive with asynchronous motor.

In the second step, using the database obtained by application of the input vector uopt to the asynchronous motor, the input/output relations of the optimal PI controller for the given controlled system are defined. Various input and output signals (e.g., depending on their actual measurability …) can be chosen for setting up the suitable database. This makes it possible to adjust the controller design to a concrete drive with concrete technical capabilities.

The relations obtained in the second step of the design procedure, written down in the form of a suitable input and output signals database (mostly in a table), are described in the third step by means of fuzzy logic principles. This way a fuzzy PI controller similar to an optimal continuous PI controller is constructed. Standard computing means for working with fuzzy systems are employed in this step, such as the fuzzy toolbox in Matlab.

The whole procedure of fuzzy PI controller design was verified by simulation of its properties in the concrete control of a drive with asynchronous motor. The results of simulation experiments show that the controller, in spite of its simplicity and the uncomplicated computer oriented design procedure applied, enables considerable improvement in the control circuit dynamic properties also in case of strongly nonlinear higher order controlled systems.

3. Multi‐motor drive optimal control using a fuzzy model‐based approach

Typical representative of multi‐motor drive is the middle part of the continuous line, where the individual working machines are coupled with each other through the material. It can be lines for processing continuous flows of material (e.g., sheet metal strips, tubes, processing lines in paper mills, and printing works) by material traction in the field of elastic or plastic deformation, which influences the material’s mechanical properties. It means that the multiple motor drives are complex and coupled MIMO nonlinear systems. Therefore, due to the complexity of their mathematical models, which parameters are difficult to identify, the development of effective control systems is quite complicated task. This chapter presents the design of optimal control of continuous production line using a fuzzy model‐based approach.

The structure of the middle part of the continuous line (further referred to as CL) is shown in Figure 9. The structure includes DC motors powered through static transistor converters TC. The working machines of the line are driven by the motors through gearbox j; v1 and v2 are machine rolls circumferential velocities, and F12 is the tension in the web of material between the two machines. The main line disturbances are tensions before and after the middle part of the considered line which are affecting the first and second drive (F01 and F23). Kv is circumferential velocity sensors, KF is tension sensor, uv1 and uv2 are outputs from velocity sensors, and uF12 is output of tension sensor. The controlling voltages u1, u2 of converters present the input variables of the system. The tension in the web of material F12 and the web of material velocity v2 is the output variables (let us consider y1 = uF12 and y2 = uv2).

Figure 9.

Structure of middle section of continuous line.

The described system with the mechanical coupling of two machines presents a third‐order nonlinear MIMO system with two inputs and two outputs (Figure 10), the parameters of the system change depending on the mechanical properties of the material and on the speed of its motion. Defining precise parameters of this nonlinear system analytically presents a rather demanding task, and therefore, it is suitable to use for its description a fuzzy system (model) built only on basis of its measured input/output data.

Figure 10.

Middle section of continuous line as MIMO system.

Various fuzzy system structures consisting of static fuzzy subsystems and their dynamic parts can be found in the literature. In setting up the structure of the fuzzy model of a continuous line, we used its state description, where the given state of the system and the given input allow us to define the subsequent state, which can be expressed mathematically by the following equation:


where u is the model’s input quantities vector, x is the state quantities vector, f is the searched for static vector function of the controlled system, and k is representing the sampling step.

Construction of the CL fuzzy model consists in determining the fuzzy approximation of this function on basis of the obtained CL inputs and outputs database. Considering the choice of CL input, state and output quantities presented in Figure 10, the structure of the proposed CL fuzzy model is shown in Figure 11.

Figure 11.

Structure of the discrete CL fuzzy model.

The whole design of CL optimal control consists of two steps:

  1. Step 1. The design of the fuzzy model for the middle section of the continuous line.

The first step in the design of the fuzzy model for the middle section of the continuous line is the establishment of a consistent database from measured inputs and their corresponding outputs, which covers its entire assumed work space and describes the behavior of the modeled system. For establishing a consistent database, we can use, for example, the method of dividing the input range into n‐levels and generating n(n − 1) transient trajectories between them [12, 31, 32], or the method of exciting the system by a random input signal [10, 30] in case it is not possible (e.g., for operational reasons) to apply a pre‐defined input signal at the system’s input. Knowledge of the structure or of the parameters of the modeled system is not required in either of these methods. To define suitable sampling time T (according to the Shannon‐Kotelnikov theorem) and approximate times for transitions for database measurement, we performed identification measurements on the physical model of the CL with input signals u1 and u2. Their value and performance are shown in Figure 12. Responses of the CL physical model output quantities to the input signals are illustrated in Figure 13.

Figure 12.

Input signals u1 and u2 for identification measurements.

Figure 13.

Identification response of CL middle section outputs.

The responses of the CL to inputs u1 and u2 (Figure 13) show that this system includes a fast tension subsystem and a slow speed subsystem. The range of input values for input u1 (web tension) is assumed within the interval [−1, 1] and the range of input values for input u2 (line speed) within the interval [−4, 4]. The database for CL fuzzy model set up will be generated so that line speed (input u2) will increase in steps each 12 s, and each one second the faster (oscillating) part of the system will be excited by input u1. The plot of input signals for generating the CL fuzzy model database is shown in Figure 14; Figure 15 shows the output quantities corresponding to these inputs.

Figure 14.

Identification of transitions for CL fuzzy model database generation.

Figure 15.

CL output variables corresponding identification transitions of inputs.

The database for CL fuzzy model was generated as demonstrated in Figure 16. With sampling time T = 0.1 s, we obtained a database with 1000 samples.

Figure 16.

Generating database for CL fuzzy model.

This measured database can be used to search for two FIS structures of the given nonlinear system which best describe the measured relations between [u1k−1, u2k−1, x1k−1, x2k−1, x3k−1] → dy1k, and [u1k−1, u2k−1, x1k−1, x2k−1, x3k−1] → dy2k.

Using the measured database, the particular fuzzy model can be designed by standardly known procedures of cluster analysis and adaptive approaches to improve the quality of modeling and reduce development time. The fundamental features of cluster analysis are reduction of the number of fuzzy rules and provision of good initial rule parameters. For our purpose from the large number of methods for adaptive fuzzy networks development [3336], we chose the adaptive neuro‐fuzzy inference system (ANFIS) with subtractive clustering [14], which is a fast and robust data analysis method, having the following parameters: range of influence = 0.4, squash factor = 1.25, accept ratio = 0.4, reject ratio = 0.01. Subtractive clustering determines the optimal clusters [34] in a multi‐dimensional input/output space that accurately represent the data [34, 37] and CL behavior. The ANFIS approach uses Gaussian functions for fuzzy sets, linear functions for the rule outputs, and Sugeno’s inference mechanism [15]. The results were two static Sugeno‐type fuzzy systems with two rules for each output quantity as is shown in Figure 17.

Figure 17.

CL fuzzy model—SUGENO type with 2 rules.

The thus obtained fuzzy systems were implemented into the final continuous line fuzzy model structure, as illustrated in Figure 11.

To verify the correctness of the CL fuzzy model, randomly generated signals u1 and u2 were applied to its input, as demonstrated in Figure 18.

Figure 18.

Performance of randomly generated signals u1 and u2.

The comparison of the fuzzy model outputs and CL physical model outputs for these inputs is shown in Figure 19.

Figure 19.

Performance of CL fuzzy model outputs y1 and y2 for randomly generated inputs u1 and u2.

The obtained results confirm that the designed fuzzy model very well approximates the performance of the continuous line also for randomly generated inputs and can be further used for the design of CL control.

  1. Step 2. Design of optimal controller for middle section of continuous line.

The principal aim of CL control consists in achieving good dynamic control of tension in the material, with the speed of material movement being in accord with the pre‐set CL speed. As it has been said above, this is in fact a nonlinear MIMO system an important feature of which is mutual influencing of the individual input and state quantities that can result in bad quality or even in the destruction of the material being processed. This fact makes the controller design methods and their subsequent resulting structures often very complex and presents an obstacle to their wider practical application in industry. Therefore, our aim was to design a simple CL controller that would ensure the desired dynamics in terms of the selected criterion for systems that are only described by input/output relationships, that is, on basis of their fuzzy model.

For control of middle section of CL (for which fuzzy model was designed), we chose the simplest control structure consisting of two standard PI controllers (one for tension control F12 and one for output velocity control v2), as illustrated in Figure 20.

Figure 20.

Continuous line PI controller structure.

Processing of material in a CL is usually carried out in operation cycles during which a required amount of prepared material is processed (e.g., a roll of paper, a sheet metal coil.) An operation cycle includes three stages—line start‐up, line running at constant processing speed, and line delayed shut‐off.

The objective of the optimization is to find such vector K[KPF, KIF, KPV, KIV] of the CL controller parameters for which the selected optimization criterion for a given CL operation cycle would be minimum. Most often this criterion is selected in quadratic form according to following equation


where e1 is the deviation between the desired and real tensile force of CL, e2 is the deviation between the desired and real output velocity of the CL. Coefficients C1 and C2 determine the importance placed on the control errors of the particular outputs. In this case, we chose C1 = 5 and C2 = 1, which in terms of physics can be interpreted as larger emphasis put on the quality of regulation of error e1 (tension in the strip of material which primarily determines its final quality). The optimal value of parameters of vector K is identified in the space of real values of gain of proportional and integrating elements of the particular controllers.

Let us note that what we are looking for is the extreme of the function of various variables, where the value of the criterial function for the individual vector K is determined by simulation on basis of the experimentally constructed CL fuzzy model according to Figure 21.

Figure 21.

Simulation diagram for computing the value of criterial function J(K).

Several procedures can be applied for the purpose of optimization (e.g., genetic algorithm methods, and network charts). Thanks to today’s availability of high‐performance computing means, we chose the method of even geometrical division of the parameter space into equal intervals and of systematic searching within the whole range of the space. The advantage of this approach is that we can always identify the global minimum of the function Eq. (6), the disadvantage may be the time and computing demands in case vector K has several parameters, and the division of the space is more dense.

At the start of the optimization process, we determined the initial values of controller parameters KPF0 = 1, KIF0 = 1, KPV0 = 1, KIV0 = 10, and we divided the parameter space of individual gains of controllers by increments ΔKPF = 1, ΔKIF = 5, ΔKPV = 1, ΔKIV = 5. The maximum values of individual parameters were defined as KPFmax = 10, KIFmax = 100, KPVmax = 10, KIVmax = 100—this was based on physically standard values of proportional and integrating elements of PI controllers.

The value of the criterial function for initial values of vector K = K0 and for the chosen CL operation cycle is equal to J(K0) = 37.19. The time responses of CL physical model output quantities corresponding to this initial PI controller parameter setting are illustrated in Figure 22.

Figure 22.

Time responses of CL outputs during operation cycle for K0.

For finding optimal controller parameters, an m‐file was created in Matlab program environment. At the end of the search process for optimal vector of CL controller parameters, the value of the optimization criterion was J(Kopt) = 0.3906, which corresponds with optimal CL controller parameter values Kopt [9, 20, 7, 80]. Time responses of CL physical model output quantities for optimally set parameters of its controller are shown in Figure 23.

Figure 23.

Time responses of CL outputs during operation cycle for Kopt.

3.1. Discussion

The proposed controller has been verified by experimental measurements on a real system which presented the physical model of the continuous line (the parameters of the CL physical model are specified in the Appendix). Figures 22 and 23 illustrate experimental results of the control of the continuous line middle section for selected operational cycle. In industrial practice, the required tension in the strip of material is set the first and then the line starts up to the desired operational speed.

Figure 22 shows the selected CL operation cycle in which first the desired value of tension in the strip of material is set to 0.8 N and at time 4 s the line starts up to reach the operational speed. At time 20 s, failure F01 = 0.7*FN occurs at entry of the line middle section (caused e.g., by a change in material thickness—material weld). The Figure shows that at initial values of parameters of CL speed and tension PI controllers, the deviation in tension from the desired value is up to 10%, the speed deviation is up to 25%, and at certain moments, the line also runs in the opposite direction. Autonomy and invariance in terms of failures are poor in dynamic states. On the contrary, when we set optimal values of parameters of CL speed and tension PI controllers, we can see—as illustrated in Figure 23—that tension in the line is maintained also in dynamic states within the range of 2% (which ensures high material processing quality during the whole operation cycle), and line speed only briefly falls outside the desired value by approx. 8% at a moment of influence by an external step disturbance. Optimal setting of the CL PI controller parameters therefore ensures good quality dynamics, autonomy and invariance of the controlled system against failures.

For the control of continuous line tension and velocity, a very simple control structure with two PI controllers was designed. We looked for four optimal parameters in the structure, such that would best satisfy the chosen quadratic optimality criterion for the given operation cycle of the line.

The quality of the designed controllers depends on the quality of the constructed fuzzy model which very well approximates the performance of the modeled system and can be further employed in the design of various CL control structures and also in the identification of non‐measurable additive disturbances influencing the system, principally in real time.

The quality of the proposed controller depends on a large extent on a good quality of the nonlinear system fuzzy model which is constructed in the first step of the design procedure The model is constructed only on basis of suitably measured relations between the system’s inputs and outputs, without the necessity of preliminary knowledge of its internal structure and parameters. The fuzzy model design is based on the basic idea of dynamic system description in state space.

The quality of the proposed controller depends on a large extent on a good quality of the nonlinear system fuzzy model which is constructed in the first step of the design procedure (see Section 4). The model is constructed only on basis of suitably measured relations between the system’s inputs and outputs, without the necessity of preliminary knowledge of its internal structure and parameters. The fuzzy model design is based on the basic idea of dynamic system description in state space.

With this method, no principal limitations for the investigated system’s nonlinearities are defined, and therefore, there is good reason to assume that the presented method will find wide use in multi‐motor drives in steel industry, paper‐making, printing and textile industries, in the production of synthetic fibers and foils in the chemical industry and in other industries.


The authors wish to thank the project VEGA 1/0464/15 for its support.

AC drive parameters:

PN = 15 kWU1N = 220 VI1N = 29.5 AnN = 1450 rev/min
MN = 98.8 NmJ = 0.11 kgm2M = 0,064 HR1 = 0.178 Ω
R2 = 0.36 ΩK11 = 277.08 H–1K12 = –269 H–1K12 = –269 H–1

UN = 24 VnN = 3650 rpm‐sRa = 0.7 ΩIN = 8.5 APN = 140 WFN = 250 N
La = 0.1 MhMN = 0.39 NmJ = 0.002 kgm2j = 24 = 0.043 VsImax = 20 A

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Pavol Fedor and Daniela Perduková (August 30th 2017). Use of Fuzzy Logic for Design and Control of Nonlinear MIMO Systems, Modern Fuzzy Control Systems and Its Applications, S. Ramakrishnan, IntechOpen, DOI: 10.5772/68050. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

351total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Vibration Suppression Controller of Multi-Mass Resonance System Using Fuzzy Controller

By Hidehiro Ikeda

Related Book

First chapter

A Real-Time Speech Enhancement Front-End for Multi-Talker Reverberated Scenarios

By Rudy Rotili, Emanuele Principi, Stefano Squartini and Francesco Piazza

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us