Open access peer-reviewed chapter

Collision Avoidance Law Using Information Amount

By Seiya Ueno and Takehiro Higuchi

Submitted: December 10th 2010Reviewed: April 28th 2011Published: September 9th 2011

DOI: 10.5772/24668

Downloaded: 1742

1. Introduction

The collision avoidance control has been one of the key technology for future transportation. Recently, many unmanned systems are developed in shapes of robots, cars, ships, aircraft, etc. In these environments, proper navigation and control systems including collision avoidance is needed. This paper is on collision avoidance control law for air vehicles under uncertain information. The control law uses information amount as one of the physical parameter for control system.

In the field of guidance, navigation, and control, collision avoidance of automated transportation system has been one of main interest of researchers. Many researches started from collision avoidance of ships (Ciletti et al., 1997) where collision avoidance has been one of the problems due to the increasing demand for the naval transportations. Wide varieties of studies on collision avoidance are treated in fields of robots (Fukuda & Kubota, 1999), cars (Hiraoka et al., 2009a, 2009b) and satellites. Some of these researches treat avoidance problems with the formation control which requires the cooperative information control (Slater et al., 2006; Stipanovic et al., 2007).

In the field of aeronautics, the Traffic alert and Collision Avoidance System (TCAS) has been one of the references for the collision avoidance. TCAS exchanges the information of aircrafts and advises the aircraft to avoid in vertical direction. For the conflicts in collision avoidance control, Frazzoli et al. (2001) have shown feasible strategy to treat the conflict problem. Gates (2009) has proposed rule-based collision avoidance control strategies for real-time online collision avoidance. Miele et al. (2010) has proposed collision avoidance control for case of abort landing with low computational load which can be calculated by on-board computer.

Conventional avoidance problems assume that all information about avoidance (intruders and environments) is certain. Therefore, control law is designed based on certain information. However in real cases, all information may not be correct and most of it is uncertain. These uncertainty of information differes by the relative position of the evader and intruder or the absolute position of intruder. There has been no research on control law to deal with uncertain information. This paper proposes control law that treats uncertain information. New parameters quantifying information amount are defined for this purpose. The proposed control law provides new performance by enabling the aircraft to obtain information and to check the certainty of the information.

Two different cases of numerical simulations are used to investigate the usage of the information amount. The first case defines the problem as the uncertainty of the information changes by the relative position of the evader and the target. The problem treats the case where the amount of infomation changes by relative position, for example, the flight in fog or smoke. The information is clearer as the evader gets closer to the fog. These uncertainties are quantified and used as parameters for collision avoidance control law. The second case defines the problem as the uncertainty of the infomation is given as absolute position. The infomation can or cannot be obtained by the position itself, for example, the flight around urban buildings or moutains. In both cases, the information amount is obtained from focused area assigned by the user. Using the information amount, the control law is designed for safer flight of the air vehicles.

2. Information amount

In this study, amounts of information are treated as parameters for the control law. First, the focused area : SE is treated as the region of the area that the user focuses. This area can be large if the vehicle is moving fast or very small if the vehicle is in urban area moving very slowly. The cleared area : SC is the area where the information are certain. In the cleared area, all of the infomation is available, meaning if there is an intruder in that region, the evader can obtain all the infomation of the intruder. In the other hand, the blurred area : SB is the area where infomation is uncertain.

From these parameters of the areas, the infomation amount is derived qunatatively as physical value to be used. One of the important factor used in this paper is information localization : IL which is calculated from the amount of cleared area out of the focused area. The IL is,

ILSC/SEE1

and the schematic image of this areas are shown in Fig.1.

Another important factor for the information amount is information acquisition requirement : IR which is requirement of the IL for safe flight. The evader selects the proper amount for the vechicle to obtain. The higher the value is, more infomation have to be obtained. In the other hand, evader have to move more when the value is high. This amount can be changed by the requirement of the user.

The infomation amounts can be changed according to the users request and experience. If the evader is moving fast, the SE might be large but the IR could be small. If the evader is trying to turn the corner, SE might be small but the IR could be very large. These amounts are similar to that of human sence of avoiding the dangers, which makes this amount unique and useful for colision avoidance.

Figure 1.

Schematic image of information amounts.

3. Collision avoidance law

The total system of collision avoidance law in this paper consists of three types of control laws. They are actual collision avoidance, information gathering, and cource keeping. They are switched by the risk of collision and amount of infomation obtained. All of the simulation are in 2-dimensions and either acceralation or angular velocity of the vehicle is used as input variables. Fig.2 shows the basic definition of variables and constants used in this paper.

Figure 2.

Definitions of variables and constants.

The risk of collision is described numerically for collision avoidance control law. Two values are introduced in this paper. One is Range to Closest Point to Approach: RCPA, which is the minimum range between two vehicles when their velocities and directional angles are kept at present value. The closest point is shown in Fig.3. This value indicates the future risk of collision. The other is Time to Closest Point to Approach: TCPA, which is time to the range between two vehicles is RCPA. Even if RCPA is small, it is not necessary to avoid quickly

Figure 3.

Image of closest point to approach.

when TCPA is large, in the othere hand, when the TCPA is small, the evader have to start the motion very quickly. The collision avoidance control law is derived from the combination of these two risk functions.

RCPA and TCPA are driven using the relative position and velocity in the body fixed coordinate shown in Fig.3. The TCPA is derived as,

TCPA=(xRuR+yRvR)/(uR2+vR2)E2

Then, the relative position and range at closest point of approach is given as,

(xCyC)=(xRyR)+TCPA(uRvR)E3
RCPA=xC2+yC2E4

The risk function is defined as the following equation.

φR02/(R02+RCPA2)E5

where R0 represents the safety range. Collision will occur when φ=1. When the minimum range between two vehicle is equal to R0, φ =0.5. The collision avoidance control law is designed to reduce φ less than 0.5 in a period of TCPA. These parameters that show the risk of collisions are used for collision avoidance. When the risk is high, the direct collision avoidance control is activated to avoid the collision. One of the examples of collision avoidance control law is described in the following.

The collision avoidance control law satisfies the following equation.

φ˙=φ/TCE6

TC is time constant that is derived from the following requirement.

TC=TCPA/ln(2φ0)E7

where φ0 is the initial value of risk function. The left hand side of Eq. (6) can be derived from the derivative of Eq. (5) and eliminating the effect of the angular velocity of the evader. The effect of the angular velocity is momentary, where they return to their original values after the avoidance. Let the absolute velocity and angle of direction of the evader be, V0 and y0, and the intruder be, V1 and y1, the relative velocity of the vehicles can be expressed as,

(uRvR)=(yRxR)ψ˙0+(cos(ψ1ψ0)1sin(ψ1ψ0)0)(V1V0)E8

The first factor of the right hand side is the angular velocity of the evader, so by eliminating this factor, the following can be derived by taking the derivative of the relative velocity.

ddt(uRvR)=(vRuR)ψ˙0E9

As total, the angular velocity for collision avoidance can be derived as,

(ω0)A=1pCvRqCuRR022φTAE10

where,

pC=xCpx+yCpyqC=xCqx+yCqypx=vR(uRpovRqi)/(uR2+vR2)qx=uR(2vRpiyR(uR2+vR2))/(uR2+vR2)py=vR(2uRpixR(uR2+vR2))/(uR2+vR2)qy=uR(vRpouRpi)/(uR2+vR2)pi=xRuR+yRvRpo=xRvRyRuRE11

This collision avoidance control law activates when the risk is high. In other words, other 2 control laws, infomation gathering and course keeping laws are used when the risk is low. In the following 2 sections, the different types of infomation gathering control laws are introduced depending on the difference of the uncertainties up ahead. The course keeping control law is used to keep the original course, which is not important in this paper, so will not be explained in details.

4. Uncertainty depending on relative position

The information gathering control with uncertainty depending on relative position is introduced in this section. The uncertainty depending on relative position stands for the cases where the infomation that can be obtained are defined as function of relative distance to an uncertainty. This is applicable for the flights in the fog or smokes where the uncertainty differs by the distance, closer you are, clear infomation you can obtain. The control target is fixed wing aircraft and the control input is angular velocity. The control law is designed from fuzzy logic to realize the fuzziness of the infomation. First, the additional parameters of the uncertain infomation is explained in this section. Followed by the control law and control results.

4.1. Uncertain parameters

The control law uses additional parameters for infomation in this section. The basic parameters were IL and IR which was explained in section 2. The following parameters quantifying uncertainty are added to design the control law dealing with uncertain information. Uncertainty of information depends on the target existence and location. Therefore, first, uncertainty parameters are defined separately. Then, the uncertainty coefficient for the control law is obtained from the uncertainty parameters for existence and location.

4.1.1. Information probability - IP

Infomation probability is a parameter describing the probability, possibility or likelihood of the target existence. IP takes a fixed value from 0 to 1 and is assigned by the user before avoidance. When IP = 0, there is no probability of existence. On the other hand, the target existence is certain when IP = 1.

4.1.2. Information clarity - IC

Infomation clarity is a parameter describing the clarity of target existence. IC varies with the quantity of information. When IC = 0, there is no information about existence. On the other hand, information on target existence becomes clear when IC = 1. When the information has uncertainly, IC varies from 0 to 1 with relative distance between the target and evader. For example, IC varies as shown in Fig. 4 when visibility is obscured by fog. Information becomes clear (IC = 1) at a certain distance and worsens gradually with further distance. No information is provided (IC = 0) beyond a certain distance. Figure 1 shows an environment when visibility is just barely secured at 5000m. The target is found definitely when the relative distance is less than 4000m. As the relative distance increases, visibility worsens gradually and the target cannot be found when the relative distance exceeds 5000m.

Figure 4.

Information clarity : The figure shows the case when the visibility is 4000 – 5000[m]. The information is clear when the distance is lower than 4000[m] and cannot be obtained when the distance is over 5000[m]. The information clarity changes linearly between two regions.

4.1.3. Information truth - IT

Infomation truth is a parameter describing truth of the target existence. The value determines whether the target exists or not. It takes a value of either 0 or 1. When IT = 0, there is no target. On the other hand, when IT = 1, there is a target.

4.1.4. Information location accuracy - IA

Infomation location accuracy is a measure of the area in which the target exists. For example, in a 2D model, IA is a circle with a radius of arbitrary length. As shown in Fig. 5, it is assumed that IA depends on the radius of the zone containing the target. IA takes a value between 0 (existence zone is vast) to 1 (existence zone is very small). In summary, the existence zone is the domain where the target may exist.

4.2. Application to control law

The control law using the uncertain infomation is introduced. From the viewpoint of complexity and difficulty, it is wrong to design a whole new control law adopting

Figure 5.

Information location accuracy: The figure shows that the unknown target is within 5000[m] of radius with highest possibility at the center.

uncertainty. Therefore, the control law to deal with uncertainty simply by introducing technique to the conventional control law with only minor modification is proposed.

4.2.1. Uncertainty coefficient

As a first step in designing a control law to deal with uncertainty, the uncertainty coefficient k based on the uncertainty parameters of existence and location (IP, IC and IT) is brought in. k is a function of these parameters expressed as,

k=f(IP,IC,IT)E12

As shown in Fig. 6, k is a coefficient introduced for the following reasons. Under the conventional control law, target information is certain and the evader flies a course for either target existence or target absence. On the other hand, when the information is uncertain, the evader flies somewhere between target existence and target absence. It is assumed that k has three components as follows, depending on relative distance: part based on original estimate (corresponding to IC = 0), part based on information that gradually becomes clear (corresponding to 0 < IC < 1), and part based on clear information (corresponding to IC = 1).

Figure 6.

The trajectory of uncertain condition and input.

In the situation in Fig. 4, when each uncertain parameter is given, the value of k becomes as shown in Fig. 7. k takes a constant value depending on IP when Rr > 5000, because the amount of information quantity does not vary in this condition. Information gradually becomes clear when 4000 < Rr < 5000. Variation of IC reflects variation of k. k takes either 0 or 1 depending on only IT when Rr < 4000, because the information becomes certain. IT is finally decided whether the target exists or not. However, if it is assumed that IT approaches a true value gradually as the information becomes clear, k is determined in real time. In such a circumstance, control input (angular velocity) ωuncert takes a value between target existence and target absence. Therefore ωuncert is expressed by the following equation using k,

wuncert=kwavo+(1k)wnotavoE13

where ωavo and ωnotavo are the control input to the evader for target existence and target absence, respectively. Avoidance control depends on the area where the target may exist, thus ωavo is a function of IA.

Figure 7.

Image of uncertainty coefficient with uncertain information.

4.2.2. Information acquisition requirement

The IR explained in section 2 is now brought in for the control. As explained in section 2, the IR is the requirement of IL by the user. IR represents the degree of need to obtain information on the focused area. IR is determined by the user and takes a value from 0 to 1. When IR = 0, there is no requirement. On the other hand, all information is required when IR = 1. As shown in Fig. 8, the clear region in the focused area becomes small as IR approaches 0. In contrast, as IR approaches 1, the cleared part becomes large. If the conventional control law is modified by introduction of the uncertainty coefficient, k, and the information acquisition requirement, IR, design of a new control law to deal with uncertainty is comparatively easy.

Figure 8.

Image of information acquisition requirement. Left figure shows the case where IR is small (vehicle passes in far region from the focused area), right figure shows the case where IR is large (vehicle passes in closer region from focused area).

4.3. Simulation result of relative position - in-fog problem

An example of avoidance problem is uncertainty of information defined in the relative coordinate (body fixed) system is shown. The problem is assumed to be in-fog problem, where there is area where the information is uncertain upahead.

4.3.1. Statement of problem

The evader cannot see the target beyond a certain distance because visibility is obscured by an obstacle like fog. The problem is defined as two-dimensional in the horizontal plane. The evader flies on a straight course with constant velocity towards a target that may exist in existence zone as shown in Fig. 9. Visibility is defined as a function of relative distance from the evader. When the relative distance is smaller than a certain distance, for example 4000m, the evader can see the target clearly. However, the evader cannot see the target when the relative distance is larger than a certain distance, for example 5000m. Visibility changes gradually between these two areas. The information clearness IC is defined depending on the relative distance to the target as for visibility as explained in Fig.4. The target is close to evader’s course, but information about existence and position are uncertain. Therefore, the target existence is given as the information probability, IP, and the position is given as the target existence zone (circlular region with radius of rE).

Figure 9.

Image of the problem, the target is assumed to be inside the fog where, the information obtained from evader is uncertain.

4.3.2. Initial conditions and requirement

The initial position of the evader and target existence zone are shown in Fig. 9. Other constants are shown in Table 1. Figure 4 is used for IC. Also the required separation between the evader and the target is set from 3000m to 4000m.

4.3.3. Initial conditions and requirement

Figures 10 and 11 show the avoidance trajectory and angular velocity, respectively. Solid lines represent the results for the proposed control law; dashed lines represent the results for the conventional control law. The figure shows two cases for the conventional control law: avoidance with correct information; and avoidance with incorrect information where target appears suddenly without information.

The avoidance trajectories in Fig. 10 show that avoidance using the conventional control law with incorrect information causes significant delay because the evader does not avoid until the target is found. On the other hand, the avoidance trajectories produced by the proposed

ParametersSymbolsValuesUncertain parametersValues
VelocityvA250[m/s]IP0.25
Initial position(0, 0)ICShown in Fig. 4
Radius of existence zonerE1000[m]IT1
Center of existence zone(25000, 2000)IAShown in Fig.5
True position of target(25500, 2500)IR0.2, 0.9

Table 1.

Parameters and constants for the simulation.

control law depend on the value of IR, because IR indicates the degree of necessary information. When IR is large, the evader must fly closer to the target and does not take early avoidance. Therefore, two stage avoidance occurs when IR = 0.9. The first stage is infomation gathering based on current infomation and information acquisition; the second stage is avoidance after finding the target. Figure 11 shows the angular velocities for avoidance. Information uncertainty reduces the sudden and severe avoidance that occurs using the conventional control law with incorrect information.

Both of the figures show that the proposed control law was able to increase the safetiness and reliability of the flight in uncertain information defined in relative position from the evader.

Figure 10.

Avoidance trajectory of the cases with/without proposed control law using uncertain parameters (IP = 0.25).

Figure 11.

Angular velocity for avoidance of the cases with/without proposed control law using uncertain parameters (IP = 0.25).

5. Uncertainty depending on absolute position

In this section, the uncertainty depending on absolute position is treated. Different from the uncertainty that differs by relative position as explained in section 4, the information does not change due to the environment. For example, when aircraft is going around a mountain or a helicopter going around the buildings, the information does not change due to relative position.

The information parameters IL and IR which was explained in section 2 is used in this section to see the effect of the information amount. The collision avoidance control explained in section 3 is brought in when the vehicle enters critical condition, otherwise the following either infomation gathering control or course keeping control takes place.

5.1. Design of information gathering control law

Design of information gathering control law is derived using a model in Fig.12 and Fig.13. Figure 12 is the vehicle in the ground fixed coordinate and Fig.13 shows the vehicle in body fixed coodinate. Focused area depends on the speed and direction of the vehicle, so the distance xP can be set by the users descision. The angle θP also depends on the level of safety which can be chosen. The area of focused area is given as follows.

SE=xP2tanθPE14

The shadow area cannot be seen from the vehicle. So, the area of shadow area is assumed as the following equation using the focused area.

SB=12(xP2xA2)(tanθPtanθA)E15

Figure 12.

Vehicle in ground fixed coordinate system.

Figure 13.

Vehicle in body fixed coordinate system.

where,

tanθA=yA/xAE16

xA and yA are relative coordinate in the body fixed frame. The time derivatives of xA and yA are given as follows.

ddt(xAyA)=(yAxA)ω0(10)V0E17

where ω0 and V0 are angular velocity and velocity of the vehicle, respectively. To derive the dynamical property of information amount, time derivative of the IL is derived.

dILdt=SBS˙ES˙BSESE2E18

SB and SE are given in Eqs.(15) and (16). It is necessary to define the xP and θP of the focused area in order to derive the time derivative. The distance of focused area xP is proportional to the velocity of helicopter V0.

xP=V0ΔtE19

where Δt is a constant with the unit of time. This shows that the vehicle moves the distance of xP in the period of Δt. The angle θP is set as constant, thus the vehicle focuses wide area in case of high speed.

The time derivatives of the SB and SE are derived.

ddt(SBSE)=(ABD0)(V˙0V0ω0)+(C0)=(ABD0)(axay)+(C0)E20

where,

A=V0(Δt)2(tanθPtanθA)B=[xAyA(tanθPtanθA)+12(xP2xA2)(1+tan2θA)]/V0C=[xA2(tanθPtanθA)12(xP2xA2)tanθA]V0/xAD=2V0(Δt)2tanθPE21

ax and ay are horizontal acceleration in the body fixed frame. Substituting Eq. (21) into Eq.(19), the time derivative of IL is given.

dILdt=Eax+Fay+GE22

where,

E=(SBDSEA)/SE2F=B/SEG=C/SEE23

Therefore, the information amount of safety can be changed by ax and ay. This shows that the direction of acceleration changes the information amount. The information gathering control law is required to keep the information amount higher than the specified value IR. Thus the shortage of information amount is defined in the following equation.

IE={IRIL(IR>IL)0(IRIL)E24

In the case of IL is less than IR, the controller is required to increase IL. On the other hand, in the case of IL is greater than IR, the high decreasing rate of IL is not desired. Therefore the following control law satisfies the both cases.

I˙L(I˙L)minE25

where,

(I˙L)min=kE(IRIL)E26

kE is a feedback gain. The left hand side of Eq. (26) is given in Eq. (23). Equation (23) has two input variables, ax and ay. Thus the minimum norm of input vector is chosen.

(ax)I=EuN,(ay)I=FuNE27

where,

uN=(kE(IRIL)G)/(E2+F2)E28

This control law uses the same feedback gain in the cases of IR>IL and IR < IL. The vehicle is required to keep the desired velocity and direction when the information amount satisfies the requirement. In this case, the following feedback law is used.

(ax)D=kV(VRV0),(ay)D=kψ(ψRψ0)V0E29

where kV and kψ are feedback gains. Finally, the information gathering control law is given as follows.

E(ax)D+F(ay)D+G(I˙L)minE30
ax=(ax)D,ay=(ay)DE31
E(ax)D+F(ay)D+G<(I˙L)minE32
ax=(ax)I,ay=(ay)IE33

The İL in the inequality conditions are given in Eq. (26). The control law is called Information Amount FeedBack (IAFB) for it feeds back the information amount as one of the parameters for information gathering.

5.2. Simulation result of absolute position

The simulation result using the collision avoidance control with IAFB is introduced. Two different cases of similation will be shown in this scetion. The first case is the case with the helicopters. The velocity of the vehicle can be changed directly by the control law. The second case is the case with the fixed wing aircraft. The input is given as the angular velocity and the velocity itself is kept as constant.

5.2.1. Simulation result of helicopters

Figure 14 shows the initial condition of the evader and intruder. The intruder cannot be seen from evader at beginning of the control. The intruder is incoming from behind the obstacle with velocity of 10[m/s] and 20[m] away from the obstacle. The evader starts from 150[m] away from the obstacle with various position defined by y(0). Figure 15 shows the simulation result of the avoidance for different initial conditions. Figure 16 shows the case without IAFB for comparison.

Figure 14.

Initial condition of evader and intruder.

Figure 15.

Trajectory of evader with IAFB starting from different initial positions.

Figure 16.

Trajectory of evader without IAFC starting from different initial positions.

In the first half of the control, the evader starts to obtain the information behind the obstacle. After the intruder is found approaching, the evader decides to evade either in front or back of the intruder depending on the estimated trajectory of the intruder passing in the way. In this case, the results were split into exactly two groups where evader decelerates and passes behind the intruder or accelerates and passes in front of intruder. In Fig.16, the case without IAFB, the first half of the information gathering does not occur, so the helicopter avoids the intruder after they find the incoming vehicle. Figure 17 is comparison of the minimum distance when the two vehicles pass each other. The result with IAFB shows higher level of avoidance due to the earlier motion of gathering information which leads to easier avoidance and faster recognition of the intruder.

Figure 17.

Trajectory of evader without IAFC starting from different initial positions.

5.2.2. Simulation result of fixed wing aircraft

For fixed wing type aircraft, it is not easy and efficient to change the velocity so often. The control input for these types are changed to angular velocity input. Basic input is same as the one described in section 5.1. Most of the conditions are same as that of the case of the helicopters except that the cruising speed of evader and intruder is 100[m/s] and the results are compared with different IR and different course of the intruder.

Figure 18.

Trajectory of evader with different IR.

Figure 18 shows the trajectory of evader with different IR, in the case of IR =0 this is same as the case without IAFB, we can see that the trajectory changes by the amount of information required.

Figure 19 shows the time history of the relative distance of the vehicles. The results show that the relative distance decreases very quickly in the case without IAFB and the minimum distance between the vehicles are shorter than the others. This clearly shows the effectiveness of the IAFB.

Figure 19.

Time history of relative distance between intruder and evader with different IR.

Figure 20 shows the trajectory of evader when the intruder starts from different positions. The IR is set as 0.9 for this simulation. In the cases when the intruder is far from the obstacle, the trajectory is smoother because the intruder is found quicker. The other two cases makes sharp turns due to the slower finding of the intruder.

Figure 20.

Trajectory of evader when the intruder approaches from different course.

6. Conclusion

Two cases of collision avoidance control is simulated to see the effect of the information amount as parameter for control. One was that uncertainty of the information changes by the relative position of the evader and the target and the other was that uncertainty of the infomation is given as absolute position. Both cases have shown smoother and safer trajectories than the conventional control laws. The simulation results have shown that the control laws using information amounts does not rely on the coodinates. The motion of the aircraft show similar trajectories to that of humans to obtain safe margin to gain information when they do not have enough information.

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Seiya Ueno and Takehiro Higuchi (September 9th 2011). Collision Avoidance Law Using Information Amount, Numerical Analysis - Theory and Application, Jan Awrejcewicz, IntechOpen, DOI: 10.5772/24668. Available from:

chapter statistics

1742total chapter downloads

1Crossref citations

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

Finite Element and Finite Difference Methods for Elliptic and Parabolic Differential Equations

By Aklilu T. G. Giorges

Related Book

First chapter

Application of the Lyapunov Exponents and Wavelets to Study and Control of Plates and Shells

By J. Awrejcewicz, V.А. Krysko, I.V. Papkova, Т.V. Yakovleva, N.A. Zagniboroda, М.V. Zhigalov, A.V. Krysko, V. Dobriyan, E.Yu. Krylovа and S.A. Mitskevich

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