Open access peer-reviewed chapter

# Optimal Trajectory Synthesis and Tracking Control for Spacecraft Large Attitude Manoeuvers

By Ranjan Vepa

Submitted: January 28th 2019Reviewed: April 24th 2019Published: January 15th 2020

DOI: 10.5772/intechopen.86498

## Abstract

The classical approach to the problem of synthesizing an optimal attitude manoeuver trajectory, involves the use of the calculus of variations and the use Lagrange multipliers or co-states. The nonlinear large attitude manoeuver trajectory is controlled by a set of nonlinear evolving co-states. In this paper, following a review of the methodologies available for trajectory synthesis followed by tracking control, the optimal trajectory for a typical optimal attitude manoeuver is synthesized by solving for the states and co-states defined by a two point boundary value problem. Gravity gradient torques are included as a matter of course. Following the synthesis of the optimal attitude-rate trajectory, tracking control laws are synthesized by re-formulating the optimal control as a feedback law. The approximate linear tracking feedback controls are evaluated by relating the optimal state and co-state vector by a linear relation. The control laws are synthesized numerically. The problem of optimal attitude orientation trajectory synthesis is also addressed. The methodologies are applied to typical sample problems and results are presented. Quantitative comparisons of the results of the methods are made to the results obtained by the application of other linear and nonlinear methods, to illustrate the key features of the methodologies.

### Keywords

• attitude manoeuvers
• optimal manoeuver trajectory
• trajectory optimization
• trajectory tracking
• feedback control laws

## 1. Introduction

The need for designing fast attitude and angular rate acquisition manoeuvers for a spacecraft with restricted or low actuation torqueing capacity arises in many space recent applications. Spacecraft are usually equipped with an attitude control system (ACS), which operates in one of two modes; the first mode involves maneuvering for attitude or angular rate acquisition while the second is to ensure stability. In the first mode, the ACS is responsible for acquiring and tracking an attitude or an angular rate state trajectory which could include a steady state. The requirements for this mode are set by the need to remotely capture an orbiting body, de-tumble a spacecraft, synchronize with another orbiting body or point at a specific direction in space. Although a large number of publications have appeared before the end of the last millennium on the subject of attitude stabilization and feedback control, a few recent papers have focused on the construction of optimal maneuvering trajectories synthesis for attitude or angular rate acquisition. There have been some publications related to the synthesis of optimal maneuvering trajectories for attitude or angular rate acquisition during the last two decades of the preceding millennium. Yet some significant advances have been made in the early part of this century. This includes papers by Lee et al. [1], Yoshida et al. [2], Aghili [3, 4], Yang and Wu [5], Liu et al. [6], and Zhang et al. [7], who have considered the maneuvering for attitude or angular rate acquisition problems using classical methodologies. Sharma and Tewari [8] have addressed the issue of nonlinear tracking of spacecraft attitude manoeuvers while Hegrenas et al. [9] have considered the maneuvering for attitude or angular rate acquisition problem by means of explicit model predictive control, via a nonlinear programming approach.

In dealing with the optimal attitude trajectory synthesis, several real world effects such as gravity gradient torques are generally neglected. Neglecting the gravity gradient torques can have a serious effect on the trajectory synthesis problem as (i) the gravity gradient torques can influence the stability of the spacecraft and (ii) they tend to couple the attitude rate dynamics with the quaternion kinematics. For this reason, it is not always advisable to ignore these torques on grounds of “smallness” as even the smallest of these perturbations can not only trigger instability but also induce the bifurcation of the orbit. It is the gravity that is primarily responsible for the orbital motion and the attitude stability of a spacecraft. The importance of gravity gradient torques has also been underscored by Lobo et al. [10].

The classical methodologies for trajectory synthesis compare well with other nonlinear and deterministic artificial intelligence approaches such as those developed by Sands et al. [11], Nakatani and Sands [12] and Baker et al. [13].

In this paper, the optimal trajectory for a typical attitude rate manoeuver is synthesized by solving for the states and co-states defined by a two point boundary value problem. Gravity gradient torques are included as a matter of course. Following the synthesis of the optimal attitude-rate trajectory, tracking control laws are synthesized by re-formulating the optimal control as a feedback law. The approximate linear tracking feedback controller gains are evaluated by relating the optimal state and co-state vector by a linear relation. The feedback control laws are synthesized numerically. The problem of optimal attitude orientation trajectory synthesis is also addressed. The optimization methodologies are applied to typical sample problems and results are presented. Quantitative comparisons of the results of the optimization method are made to the results obtained by the application of other linear and nonlinear control methods, to illustrate the key features of the methodologies.

## 2. Spacecraft attitude dynamics and quaternion kinematics

In matrix form, when the inertia matrix is not diagonal the equations of attitude motion of chaser spacecraft are

Iω̇+ΩIω=M+Mgg+Md,E1

where Iis the moment of inertia matrix which is assumed to be

I=I11I12I13I12I22I23I13I23I33,ωω1ω2ω3,Ω=0ω3ω2ω30ω1ω2ω10E2

Mdare the disturbance torques and Mggare the gravity gradient torques.

It is important to emphasize that the targets dynamics are irrelevant to us as there is little or no chance of acquiring the target’s inertia properties. However the target’s angular velocity vector is assumed to be given by ωd, its attitude quaternion relative to the chaser’s body frame is assumed to be qdor its relative attitude quaternion Δq, relative to the chaser’s body frame, can in principle be measured from within the chaser spacecraft.

Expressions for the gravity gradient moment are obtained assuming that zaxis of the spacecraft body is nominally pointing to the Earth. The direction vector the center of gravity of the spacecraft pointing to the Earth is given by the last column of TBR, the transformation from the Earth orbiting frame to the body fixed frame of the spacecraft as

c=c1c2c3T.E3

The corresponding cross product operator c×is defined as

c×=0c3c2c30c1c2c10.E4

Hence the gravity gradient moments acting on the spacecraft and manipulator body are:

Mgg=3n2c×Ic3n2c×Ic=LggMggNggT.E5

Thus,

Mgg=LggMggNgg=3n2c2c1I31+c22I32+c2c3I33I22c3c1I21c32I23c3c1I11I33+c2c3I12+c32I13c12I31c1c2I32c12I21+c1c2I22I11+c1c3I23c22I12c2c3I13.E6

If we express the transformation from the orbiting frame to the body coordinates in terms of an attitude quaternion of the chaser spacecraft with components ε1, ε2, ε3and ηas

TBRq=η2+ε12ε22ε322ε1ε2+ε3η2ε1ε3ε2η2ε1ε2ε3ηη2ε12+ε22ε322ε2ε3+ε1η2ε1ε3+ε2η2ε2ε3ε1ηη2ε12ε22+ε32,E7

then from the last column, the Earth pointing direction vector is:

c=c1c2c3=2ε1ε3ε2η2ε2ε3+ε1ηη2ε12ε22+ε32.E8

The quaternion kinematics satisfies

dqdt=12Aωωq,qTq=1E9

where the quaternion q=ε1ε2ε3ηT, consists of a vector part, ε=ε1ε2ε3Tand the scalar ηso,

q=εηandAω=ΩωωωT0,Ωω=0ω3ω2ω30ω1ω2ω10.E10

The quaternion kinematics may also be compactly expressed as

dqdt=ddtεη=12Γqω,E11
Γq=ηI3×3+SεεT,Sε=0ε3ε2ε30ε1ε2ε10,E12

where I3×3is the 3×3unit matrix. These relations may be inverted as

ω=2ηI3×3+STεεε̇η̇T=2Γ1qε̇η̇T.E13

The desired attitude quaternion relative to the chaser’s body frame which is assumed to be qdand the relative attitude quaternion Δq, relative to the chaser’s body frame are related to the chasers attitude by

qd=Δqq.E14

Given two quaternions, Δq=Δq1Δq2Δq3Δq4T, q=ε1ε2ε3ηT, the quaternion product qd=Δqqis defined as

qd=q1dq2dq3dq4d=ηε3ε2ε1ε3ηε1ε2ε2ε1ηε3ε1ε2ε3ηΔq1Δq2Δq3Δq4=Δq4Δq3Δq2Δq1Δq3Δq4Δq1Δq2Δq2Δq1Δq4Δq3Δq1Δq2Δq3Δq4ε1ε2ε3η.E15

Hence it is expressed in matrix form as

qd=C0Δq,C0=ηε3ε2ε1ε3ηε1ε2ε2ε1ηε3ε1ε2ε3η,C0=ηI3×3+STεq1:3q1:3η.E16

Similarly,

qd=q1dq2dq3dq4d=Δq4Δq3Δq2Δq1Δq3Δq4Δq1Δq2Δq2Δq1Δq4Δq3Δq1Δq2Δq3Δq4ε1ε2ε3η=Δq4I3×3+SΔq1:3Δq1:3Δq1:3Δq4q.E17

The traditional method of defining the attitude of a spacecraft is by the use of Euler angle sequences. The conversion of Euler angles defined as Euler angle sequences, may be converted to an equivalent quaternion set, using well-known conversion formulae, such as, those given by Smeresky et al. [14].

## 3. Formulation of the optimal angular rate trajectory synthesis problem

The first task is to formulate the optimal control problem, so it can be solved numerically. This is briefly reviewed. The attitude equations of the spacecraft may be expressed in state space form as

dqdt=12Aωωq,qTq=1,E18
Iω̇+ΩIω=Mc+Mgg+Md,E19

where Mcis the control torque vector acting on the spacecraft. For our purposes it will be assumed that it can be expressed as, Mc=ITWu, where Tis the scalar magnitude of the specific torque or torque per unit inertia, Wis a symmetric, positive-definite, torque transformation weighting matrix and udefines the direction of the torque vector. It is the attitude steering control input to the spacecraft. Thus Eq. (19) is

ω̇=I1MggqI1ΩIω+TWu,E20

with

u=sinαcosβcosαcosβsinβT.E21

When one is interested in the problem of finding the steering control

u=ut,t0ttf,E22

the torque direction time history is sought, such that it minimizes the cost functional:

J=0.5ωtωdTQfωtωdt=tf=Φωtt=tf,E23

subject to, Eqs. (18), (20) and (21).

Introducing the single state vector, x=qTωTT, so the Eqs. (18), (20) and (22) are expressed as

ddtxT=ddtqTωT=fT,qTq1g=0E24

To solve the optimization problem, seven Lagrangian multipliers or co-states are introduced given by the two vectors λqtand λωt, and a scalar λcdenoted by a single column vector, λt. Following Bryson and Ho [15], a Hamiltonian function is defined as

H=λTtfTtgT.E25

Hence,

H=λqT12Aωωq+λωTI1MggqI1ΩIω+λωTTWu+λcg.E26

The necessary conditions (Bryson and Ho [15], Conway [16]) for the first variation of Jto be zero include the co-state differential equations

ddtλTt=Hx.E27

Explicitly the co-state equations are

λ̇qit=Hqi=λqT12Aωωqiq+λωTI1qiMggq,E28
λ̇ωit=Hωi=λqT12ωiAωωqλωTI1ωiΩIω,E29
λ̇c=0.E30

The optimality conditions are

Hα=λωTTWcosαcosβsinαcosβ0T=0,E31

and

Hβ=λωTTWsinαsinβcosαsinβcosβT=0.E32

Hence,

Wλω=WλωsinαcosβcosαcosβsinβT=Wλωu.E33

Thus the two-parameter control vector u, can be expressed as

u=Wλω/Wλω.E34

The choice of the sign in Eq. (34) will depend on the direction of the desired torque, forward or reverse torque. Thus the closed-loop equations of motion are

ddtx=f,qTq1g=0,u=Wλω/Wλω.E35

To complete the definition of the optimal solution, the boundary conditions at t=tffor the co-state system are found by applying the transversality conditions.

The transversality conditions ensure that the initial and final states are selected optimally within the feasible regions of the states. For the transversality conditions, one may write

λqtf=∂Φωtqt=tf=0,E36
λωtf=∂Φωtωt=tf=Qfωtfωd.E37
λctf=0.E38

Thus, λct=0is a feasible solution. If it can be ensured that the constraint on the quaternion, qTq1g=0is satisfied at each and every integration time step, λctcould be set to zero for all time.

The solution to the state and co-state equations, subject to the optimal control defined by Eq. (34) and the boundary conditions defined by Eqs. (36) and (37), may be found by solving a two point boundary value problem (TPBVP). This can be done using MATLAB’s function, bvp4c.m.

Re-considering Eq. (34), the control input vector may be expressed as

u=Wλω/Wλω=r1BTλ,E39

where BTis a projection matrix relating λto Wλωand r=Wλω. It is important to recognize that the weighting matrix Walso needs to be chosen. On one hand it provides a set of free parameters so one can construct an optimal solution, while on the other it makes the solution much harder to obtain. Its choice is discussed further in Section 8.

## 4. Feedback implementation of optimal co-states

Of interest at this stage is to be able to implement the controller, obtained in the last section, as a feedback control law. Thus, inspired by linear optimal control, it can be assumed that, locally, the co-state vector λis proportional to the state vector x. Hence one may express

λ=Px.E40

Furthermore it is assumed that the matrix P is slowly varying and hence does not change appreciably as the time tchanges from the current time tto t±Δtand to t±nΔt, n=2,3,4. Thus the matrix P may be obtained by evaluating λand xin Eq. (40) at the times t±nΔt, n=0,1,,4, where Δtis a reasonably small time step, over which both λand xchange appreciably. The solution for the matrix Pis obtained by differencing the data and solving an over determined system of linear equations by a least squares approach over a moving time window. The matrix Pis also constrained to be a symmetric non-negative definite matrix. Moreover xis expressed as

x=xxref+xref,E41

where xrefis the optimal transfer trajectory for xobtained by solving the system of equations for the states, parameters and co-states defined in Section 3, which together constitute a TPBVP. By solving the system of equations for the states, parameters and co-states defined in Section 3, one also seeks the trajectory coordinates of the reference trajectory and the total transfer time. It is important that the spacecraft’s attitude controller is able to track the reference trajectory. Consequently the control input is expressed as

u=r1BTPxxref+r1BTλref,λref=Pxref.E42

The second of Eq. (42) is used to obtain the matrix, P. Like in Eq. (34), the choice of the sign in the first of Eq. (42) will depend on the direction of the desired torque, forward or reverse torque.

## 5. Simplified formulation of the optimal angular rate trajectory synthesis problem

Some authors (for example, see Aghili [4]) have formulated the attitude rate acquisition problem without including the gravity gradient torques. Thus the quaternion kinematics could be ignored. Thus in this case one could set λqt0, which results in considerable simplification of the trajectory synthesis problem. The downside of the approach is that the quaternion kinematics is ignored and consequently the quaternion attitude could be quite arbitrary. Quite often after a de-tumbling manoeuver, a precise orientation must be acquired. The required attitude could be acquired in an independent manoeuver and the methodology for this is developed in the next section. The associated tracking problem, which involves tracking the complete state vector set point, must then be separately addressed. Typically this is done by using a barrier Lyapunov function as illustrated by Vepa [17].

## 6. Optimal attitude orientation acquisition trajectory synthesis

To begin with the quaternion kinematics is given by Eqs. (9) and (11), and can be expressed in one of two alternate forms as.

dqdt=12Aωωq=12Γqω.E43

In Eq. (43), the angular velocity vector is treated as a control variable and expressed as

ω=ωmaxu,E44

where the direction vector uis parametrized by an equation similar to Eq. (21). Thus,

u=sinαcosβcosαcosβsinβT.E45

When one is interested in the problem of finding the directional control

u=ut,t0ttf,E46

the angular velocity direction time history is sought, such that it minimizes the cost functional:

J=0.5qtqdTQfqtqdt=tf=Φqtt=tf,E47

subject to, Eqs. (43), (44) and (45). The corresponding Hamiltonian function is

H=λqT12Aωωq=λqT12Γqω=ωmax2λqTΓqu.E48

The corresponding co-state differential equations are

ddtλqTt=Hq=12λqTAωω.E49

By using an argument similar to the one used in developing Eqs. (31)(34), the optimal control is given by

u=ΓTqλq/ΓTqλq=Γqλq/Γqλq.E50

For the co-state boundary conditions one has

λqtf=∂Φqtqt=tf=Qfqtfqd.E51

Once the control is found from Eqs. (50) and (44) is used to define the angular velocity vector and Eqs. (19) and (20) to define the optimal input control torque.

## 7. Shape based optimal trajectory synthesis

An alternative approach to the optimization based on the integration of co-states is to use a shape based approach as outlined by Caubet and Biggs [18, 19]. For purposes of comparison the shape based approach serves as a useful alternative. In a shape based approach, each of the quaternion components are expressed as a summation of polynomials in terms of a time variable, multiplied by coefficients which may be determined by applying the relevant boundary conditions at the initial and final values of the time variable over a finite time frame. Thus, for example, the quaternion components are expressed as

qi=qi0+q̇i0tfttf+q̇i0q̇it=tftfttf21ttf+qifqi0q̇i0tft2tf232ttf+j=5N+4e¯ij4j4ttf2j3ttf3+ttfj1,E52

where the coefficients e¯ij4are yet to be determined. They are determined by minimizing the cost function

J=01q12+q22+q32+q421dttf2.E53

Once all the coefficients of the quaternion components qiare determined, the angular velocity vector is defined by the inverse of the relation given by Eq. (43) which is

ω=q4q3q2q1q3q4q1q2q2q1q4q3q̇1q̇2q̇3q̇4.E54

The angular velocity vector ωis evaluated for a range of non-dimensional time values between 0 and 1. From the ratio of the maximum of this set, defining the angular velocity time history in terms of the non-dimensional time variable, and the maximum allowable angular velocity magnitude, the length of the time frame tfover which the control torques must be applied may be found. From the angular velocity vector ωthe torques that must be applied to the satellite including the gravity gradient torques may be found. From a range of choices for N(say 1N6) in the Eq. (52) defining the quaternion, the one that gives the lowest value for tfis selected.

## 8. Typical simulation examples

The first example considered the attitude dynamics is defined by Eqs. (18) and (19). Thus the gravity gradient torques acting on the spacecraft are included in the dynamic model and they are responsible for coupling the attitude quaternion kinematics and the angular velocity dynamics. The objective is to spin the spacecraft so the final angular velocity vector is given by ω=111Trads/s. The initial angular velocity vector is ω=000T. The spacecraft is fitted with magnetic torque actuators and the maximum three axis torques are limited to Tc=0.6211TNm. The diagonal non-zero elements of the weighting matrix Win Eq. (34) and the principal moment of inertia of the spacecraft are respectively given by

Wdiag=0.240.91,I=2.272933.273310.3989kgm2.E55

In all examples the solution of the TPBVP is done using MATLAB’s function, bvp4c.m. Whenever there was a need to solve an initial value problem, the equations were integrated using MATLAB’s ode45.m.

In all cases, the time variable was made non-dimensional so it raged from [0–1]. The integration time step was chosen to be relatively small initially (Δt=0.0001), and automatically and iteratively reduced linearly as the final time was approached. The iterations were terminated when no further improvement in the accuracy of the predicted final time was feasible.

Figure 1 shows the time history of the reference quaternion components and Figure 2 shows the corresponding angular velocities (p, q, r). Figures 3 and 4 show the corresponding, actual, quaternion components and the corresponding angular velocity components, where an approximate optimal linear feedback law based on Eq. (42) is used to track the reference trajectory.

In the next example, the simplified attitude dynamics is used with the gravity gradient torques neglected. This decouples the angular rate dynamics from the attitude quaternion dynamics, which need to be considered for synthesizing the reference trajectories. The reference angular velocities are then integrated to obtain the spacecraft’s quaternion attitude time history. Figure 5 shows the reference angular velocity components.

Figure 6 shows the errors in the actual angular velocity components when compared with corresponding reference values and Figure 7 the corresponding quaternion components. Figure 8 shows the attitude in terms of the Euler axis and the Euler principal angle components.

In the final example it is desired to alter the attitude quaternion of the spacecraft, so as to point the spacecraft in a desired direction. In this case on the quaternion kinematics defined by Eq. (44) are used. The maximum angular velocity of the spacecraft is assumed to be limited to 0.001 rad/s. The desired pointing direction is specified as a yaw, roll pitch Euler angle sequence given by ψθϕ=22°25°30°corresponding to the components of the quaternion qd=0.20680.25180.16820.9303.

In Figure 9 are shown the reference optimal quaternion components and in Figure 10 are shown the corresponding angular velocity components. The optimum torque components required to affect the attitude change are shown in Figure 11. These include compensation for the gravity gradient torques. They show that they could be easily achieved by low thrust electric actuators such as electro-spray thrusters.

When small reaction wheels are used Eq. (19) may be modified to include the momentum of the wheels and the control inputs to the wheels could also be estimated. If reaction wheels are used much larger torques are possible and the time over which they are used could be shortened. In Figure 12 are shown the attitude time history components in terms of the Euler axis and Euler principal angle components.

This example is also solved using the shape based approach discussed briefly in Section 7. In Figure 13 are shown the required applied torque components obtained by the shape based approach with N=1in Eq. (52). The time frame over which the control must be applied is tf=733.3s, which is the lowest for all Nconsidered and is about the same as the time required by the approach based on the integration of co-states. However the reference torque components shown in Figure 13, are much larger than those plotted in Figure 11, and for this reason, they are not referred to as ‘optimum torques.’ The corresponding gravity gradient torques, acting on a satellite orbiting the Earth at the geostationary orbit radius, are also shown in Figure 14. Although the gravity gradient torques are of the same orders of magnitude as the reference optimum torque components in Figure 11, they are much smaller than the corresponding torque components obtained by the shape based approach and shown in Figure 13. It must be recognized that the gravity gradient torques become much larger as the spacecraft orbits the Earth at a much closer orbit radius.

## 9. Discussion and conclusions

A close examination of the results in Figures 9 and 12 shows that while the quaternion component time histories are not linearly varying, the Euler axis and Euler principal angle components are almost linear. This allows for linear extrapolation of the trajectories if when desired. It is also observed that the acquisition of the Euler axis is relative fast in comparison with growth rate of the Euler angle which is relatively slower. The kinematics of the Euler axis seems to represent a fast sub-system while the growth of the Euler angle represents the slow sub-system. This observation, facilitates the construction of approximate sub-optimal trajectories where in the Euler axis is acquired instantly and the Euler angle increases or decreases linearly with time. Once a sub-optimal solution can be defined in terms of interpolating polynomial it can also be further optimized by using Eq. (13) and shape based methods as those proposed by Caubet and Biggs [18, 19], quite rapidly and if need be, by the pseudo spectral method, or other direct collocation methods. The advantage of further optimization using shape based methods is that the precise shape of the desired output could be achieved avoiding overshoot. However depending on the choice of the output shape function, the control could be restrictive and so the magnitudes of the torques required could be much larger in comparison with the co-states approach.

For the preceding example, where a set angular velocity components were desired, the velocity components and the corresponding quaternion components are shown in Figures 5, 6, 7 as they vary with time. Figure 8 also shows that the Euler axis and Euler principal angle components are varying as quadratic functions of time. A similar conclusion cannot be drawn as far as the angular velocity components and the components of the quaternion. This again is extremely useful in applying low order polynomials for developing formulae for extrapolating the optimal trajectories, by converting the quaternion components to the domain of the Euler axis and Euler principal angle components. It also facilitates the integration of various optimal segments into a single trajectory over an extended time frame.

From the first example, comparing Figures 1 and 2 with Figures 3 and 4, it is seen that the optimum tracking feedback control law obtained by linearly approximating the relationship between the states and co-states by Eq. (40), performs well. The errors between these two sets of trajectories, the reference trajectory and the actual tracked trajectory, are always within 5% of the corresponding reference value, over the time frame of the plots.

In this paper, either the simplest form of the attitude dynamics or the basic kinematic equations alone are used to construct the optimal trajectories. The required control torques are obtained from the inverse dynamic relations. The usefulness of transforming the attitude representation to the Euler axis and Euler principal angle components, as it facilitates the application of low order polynomials for the construction of approximate sub-optimal trajectories, is demonstrated. Furthermore it is shown how optimal feedback control laws may be constructed from the solution for the optimal trajectories, for tracking the reference trajectories.

## Conflict of interest

There are no conflicts of interest to declare.

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Ranjan Vepa (January 15th 2020). Optimal Trajectory Synthesis and Tracking Control for Spacecraft Large Attitude Manoeuvers, Advances in Spacecraft Attitude Control, Timothy Sands, IntechOpen, DOI: 10.5772/intechopen.86498. Available from:

### Related Content

#### Advances in Spacecraft Attitude Control

Edited by Timothy Sands

Next chapter

#### Parameter Optimization for Spacecraft Attitude Stabilization Using Magnetorquers

By Renato Bruni and Fabio Celani

#### Deterministic Artificial Intelligence

Edited by Timothy Sands

First chapter

#### Stochastic Artificial Intelligence: Review Article

By T.D. Raheni and P. Thirumoorthi

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.