Multi-rotor with suspended load: System dynamics and control toolbox

There is an increasing demand for Unmanned Aerial Systems (UAS) to carry suspended loads as this can provide significant benefits to several applications in agriculture, law enforcement and construction. The load impact on the underlying system dynamics should not be neglected as significant feedback forces may be induced on the vehicle during certain flight manoeuvres. The constant variation in operating point induced by the slung load also causes conventional controllers to demand increased control effort. Much research has focused on standard multi-rotor position and attitude control with and without a slung load. However, predictive control schemes, such as Nonlinear Model Predictive Control (NMPC), have not yet been fully explored. To this end, we present a novel controller for safe and precise operation of multi-rotors with heavy slung load in three dimensions. The paper describes a System Dynamics and Control Simulation Toolbox for use with MATLAB/SIMULINK which includes a detailed simulation of the multi-rotor and slung load as well as a predictive controller to manage the nonlinear dynamics whilst accounting for system constraints. It is demonstrated that the controller simultaneously tracks specified waypoints and actively damps large slung load oscillations. A linear-quadratic regulator (LQR) is derived and control performance is compared. Results show the improved performance of the predictive controller for a larger flight envelope, including aggressive manoeuvres and large slung load displacements. The computational cost remains relatively small, amenable to practical implementations.

Improving PLDA speaker verification using WMFD and linear-weighted approaches in limited microphone data conditions

This paper proposes the addition of a weighted median Fisher discriminator (WMFD) projection prior to length-normalised Gaussian probabilistic linear discriminant analysis (GPLDA) modelling in order to compensate the additional session variation. In limited microphone data conditions, a linear-weighted approach is introduced to increase the influence of microphone speech dataset. The linear-weighted WMFD-projected GPLDA system shows improvements in EER and DCF values over the pooled LDA- and WMFD-projected GPLDA systems in inter-view-interview condition as WMFD projection extracts more speaker discriminant information with limited number of sessions/ speaker data, and linear-weighted GPLDA approach estimates reliable model parameters with limited microphone data.

Quantification of stability in an agility drill using linear and nonlinear measures of variability

This study implemented linear and nonlinear methods of measuring variability to determine differences in stability of two groups of skilled (n = 10) and unskilled (n = 10) participants performing 3m forward/backward shuttle agility drill. We also determined whether stability measures differed between the forward and backward segments of the drill. Finally, we sought to investigate whether local dynamic stability, measured using largest finite-time Lyapunov exponents, changed from distal to proximal lower extremity segments. Three-dimensional coordinates of five lower extremity markers data were recorded. Results revealed that the Lyapunov exponents were lower (P < 0.05) for skilled participants at all joint markers indicative of higher levels of local dynamic stability. Additionally, stability of motion did not differ between forward and backward segments of the drill (P > 0.05), signifying that almost the same control strategy was used in forward and backward directions by all participants, regardless of skill level. Furthermore, local dynamic stability increased from distal to proximal joints (P < 0.05) indicating that stability of proximal segments are prioritized by the neuromuscular control system. Finally, skilled participants displayed greater foot placement standard deviation values (P < 0.05), indicative of adaptation to task constraints. The results of this study provide new methods for sport scientists, coaches to characterize stability in agility drill performance.

Study of robust H(infinity) filtering application in loosely coupled INS/GPS system

Since a celebrate linear minimum mean square (MMS) Kalman filter in integration GPS/INS system cannot guarantee the robustness performance, a H(infinity) filtering with respect to polytopic uncertainty is designed. The purpose of this paper is to give an illustration of this application and a contrast with traditional Kalman filter. A game theory H(infinity) filter is first reviewed; next we utilize linear matrix inequalities (LMI) approach to design the robust H(infinity) filter. For the special INS/GPS model, unstable model case is considered. We give an explanation for Kalman filter divergence under uncertain dynamic system and simultaneously investigate the relationship between H(infinity) filter and Kalman filter. A loosely coupled INS/GPS simulation system is given here to verify this application. Result shows that the robust H(infinity) filter has a better performance when system suffers uncertainty; also it is more robust compared to the conventional Kalman filter.

Spin-state equilibria in the LCoO3 system from 59Co NMR

Spin-state equilibria in the whole set of LCoO3 (where L stands for a rare-earth metal or Y) have been investigated with the use of 59Co NMR as a probe for the polycrystalline samples (except Ce) in the temperature interval 110-550 K and frequency range 3- 11.6 MHz. Besides confirming the coexistence of the high-spin—low-spin state in this temperature range, a quadrupolar interaction of ∼0.1 -0.5 MHz has been detected for the first time from 59Co NMR. The NMR line shape is found to depend strongly on the relative magnitude of the magnetic and quadrupolar interactions present. Analysis of the powder pattern reveals two basically different types of transferred hyperfine interaction between the lighter and heavier members of the rare-earth series. The first three members of the lighter rare-earth metals La, Pr (rhombohedral), and Nd (tetragonal), exhibit second-order quadrupolar interaction with a zero-asymmetry parameter at lower temperatures. Above a critical temperature TS (dependent on the size of the rare-earth ion), the quadrupolar interaction becomes temperature dependent and eventually gives rise to a first-order interaction thus indicating a possible second-order phase change. Sm and Eu (orthorhombic) exhibit also a second-order quadrupolar interaction with a nonzero asymmetry parameter ((η∼0.47)) at 300 K, while the orthorhombic second-half members (Dy,..., Lu and Y) exhibit first-order quadrupolar interaction at all temperatures. Normal paramagnetic behavior, i.e., a linear variation of Kiso with T-1, has been observed in the heavier rare-earth cobaltites (Er,..., Lu and Y), whereas an anomalous variation has been observed in (La,..., Nd)CoO3. Thus, Kiso increases with increasing temperature in PrCoO3 and NdCoO3. These observations corroborate the model of the spin-state equilibria in LCoO3 originally proposed by Raccah and Goodenough. A high-spin—low-spin ratio, r=1, can be stabilized in the perovskite structure by a cooperative displacement of the oxygen atoms from the high-spin towards the low-spin cation. Where this ordering into high- and low-spin sublattices occurs at r=1, one can anticipate equivalent displacement of all near-neighbor oxygen atoms towards a low-spin cobalt ion. Thus the heavier LCoO3 exhibits a small temperature-independent first-order quadrupolar interaction. Where r<1, the high- and low-spin states are disordered, giving rise to a temperature-dependent second-order quadrupolar interaction with an anomalous Kiso for the lighter LCoO3.

On the direct parallel solution of systems of linear equations: New algorithms and systolic structures

The paper presents two new algorithms for the direct parallel solution of systems of linear equations. The algorithms employ a novel recursive doubling technique to obtain solutions to an nth-order system in n steps with no more than 2n(n −1) processors. Comparing their performance with the Gaussian elimination algorithm (GE), we show that they are almost 100% faster than the latter. This speedup is achieved by dispensing with all the computation involved in the back-substitution phase of GE. It is also shown that the new algorithms exhibit error characteristics which are superior to GE. An n(n + 1) systolic array structure is proposed for the implementation of the new algorithms. We show that complete solutions can be obtained, through these single-phase solution methods, in 5n−log2n−4 computational steps, without the need for intermediate I/O operations.

Asymptotic behavior of a hierarchical system of learning automata

Learning automata arranged in a two-level hierarchy are considered. The automata operate in a stationary random environment and update their action probabilities according to the linear-reward- -penalty algorithm at each level. Unlike some hierarchical systems previously proposed, no information transfer exists from one level to another, and yet the hierarchy possesses good convergence properties. Using weak-convergence concepts it is shown that for large time and small values of parameters in the algorithm, the evolution of the optimal path probability can be represented by a diffusion whose parameters can be computed explicitly.

Analysis of an internally mountable accelerometer balance system for use with non-isotropic models in shock tunnels

Knowledge of drag force is an important design parameter in aerodynamics. Measurement of aerodynamic forces at hypersonic speed is a challenge and usually ground test facilities like shock tunnels are used to carry out such tests. Accelerometer based force balances are commonly employed for measuring aerodynamic drag around bodies in hypersonic shock tunnels. In this study, we present an analysis of the effect of model material on the performance of an accelerometer balance used for measurement of drag in impulse facilities. From the experimental studies performed on models constructed out of Bakelite HYLEM and Aluminum, it is clear that the rigid body assumption does not hold good during the short testing duration available in shock tunnels. This is notwithstanding the fact that the rubber bush used for supporting the model allows unconstrained motion of the model during the short testing time available in the shock tunnel. The vibrations induced in the model on impact loading in the shock tunnel are damped out in metallic model, resulting in a smooth acceleration signal, while the signal become noisy and non-linear when we use non-isotropic materials like Bakelite HYLEM. This also implies that careful analysis and proper data reduction methodologies are necessary for measuring aerodynamic drag for non-metallic models in shock tunnels. The results from the drag measurements carried out using a 60 degrees half angle blunt cone is given in the present analysis.

Soil properties and crop yields in a dryland Vertisol sown with cotton-based crop rotations

Information on the effects of growing cotton (Gossypium hirsutum L.)-based crop rotations on soil quality of dryland Vertisols is sparse. The objective of this study was to quantify the effects of growing cereal and leguminous crops in rotation with dryland cotton on physical and chemical properties of a grey Vertisol near Warra, SE Queensland, Australia. The experimental treatments, selected after consultations with local cotton growers, were continuous cotton (T1), cotton-sorghum (Sorghum bicolor (L.) Moench.) (T2), cotton-wheat (Triticum aestivum L.) double cropped (T3), cotton-chickpea (Cicer arietinum L.) double cropped followed by wheat (T4) and cotton-wheat (T5). From 1993 to 1996 land preparation was by chisel ploughing to about 0.2 m followed by two to four cultivations with a Gyral tyne cultivator. Thereafter all crops were sown with zero tillage except for cultivation with a chisel plough to about 0.07-0.1 m after cotton picking to control heliothis moth pupae. Soil was sampled from 1996 to 2004 and physical (air-filled porosity of oven-dried soil, an indicator of soil compaction; plastic limit; linear shrinkage; dispersion index) and chemical (pH in 0.01 M CaCl2, organic carbon, exchangeable Ca, Mg, K and Na contents) properties measured. Crop rotation affected soil properties only with respect to exchangeable Na content and air-filled porosity. In the surface 0.15 m during 2000 and 2001 lowest air-filled porosity occurred with T1 (average of 34.6 m3/100 m3) and the highest with T3 (average of 38.9 m3/100 m3). Air-filled porosity decreased in the same depth between 1997 and 1998 from 45.0 to 36.1 m3/100 m3, presumably due to smearing and compaction caused by shallow cultivation in wet soil. In the subsoil, T1 and T2 frequently had lower air-filled porosity values in comparison with T3, T4 and T5, particularly during the early stages of the experiment, although values under T1 increased subsequently. In general, compaction was less under rotations which included a wheat crop (T3, T4, T5). For example, average air-filled porosity (in m3/100 m3) in the 0.15-0.30 m depth from 1996 to 1999 was 19.8 with both T1 and T2, and 21.2 with T3, 21.1 with T4 and 21.5 with T5. From 2000 to 2004, average air-filled porosity (in m3/100 m3) in the same depth was 21.3 with T1, 19.0 with T2, 19.8 with T3, 20.0 with T4 and 20.5 with T5. The rotation which included chickpea (T4) resulted in the lowest exchangeable Na content, although differences among rotations were small. Where only a cereal crop with a fibrous root system was sown in rotation with cotton (T2, T3, T5) linear shrinkage in the 0.45-0.60 m depth was lower than in rotations, which included tap-rooted crops such as chickpea (T4) or continuous cotton (T1). Dispersion index and organic carbon decreased, and plastic limit increased with time. Soil organic carbon stocks decreased at a rate of 1.2 Mg/ha/year. Lowest average cotton lint yield occurred with T2 (0.54 Mg/ha) and highest wheat yield with T3 (2.8 Mg/ha). Rotations which include a wheat crop are more likely to result in better soil structure and cotton lint yield than cotton-sorghum or continuous cotton.

Evaluation of a remote drafting system for regulating sheep access to supplement

Remote drafting technology now available for sheep makes possible targeted supplementation of individuals within a grazing flock. This system was evaluated by using 68 Merino wethers grazing dry-season, native Mitchell grass pasture (predominantly Astrebla spp.) as a group and receiving access to lupin grain through a remote drafter 0, 1, 2, 4 or 7 days/week for 8 weeks. The sole paddock watering point was separately fenced and access was via a one-way flow gate. Sheep exited the watering point through a remote drafter operated by solar power and were drafted by radio frequency identification (RFID) tag, according to treatment, either back into the paddock or into a common supplement yard where lupins were provided ad libitum in a self-feeder. Sheep were drafted into the supplement yard on only their first time through the drafter during the prescribed 24-h period and exited the supplement yard via one-way flow gates in their own time. The remote drafter operated with a high accuracy, with only 2.1% incorrect drafts recorded during the experimental period out of a total of 7027 sheep passes through the remote drafter. The actual number of accesses to supplement for each treatment group, in order, were generally less than that intended, i.e. 0.02, 0.69, 1.98, 3.35 and 6.04 days/week. Deviations from the intended number of accesses to supplement were mainly due to sheep not coming through to water on their allocated day of treatment access, although some instances were due to incorrect drafts. There was a non-linear response in growth rate to increased frequency of access to lupins with the growth rate response plateauing at similar to 3 actual accesses per week, corresponding to a growth rate of 72.5 g/head. day. This experiment has demonstrated the application of the remote drafting supplementation system for the first time under grazing conditions and with the drafter operated completely from solar power. The experiment demonstrates a growth response to increasing frequency of access to supplement and provides a starting point with which to begin to develop feeding strategies to achieve sheep weight-change targets.

Exact solutions of coupled multispecies linear reaction–diffusion equations on a uniformly growing domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.

A conditionally linearized Monte Carlo filter in non-linear structural dynamics

State and parameter estimations of non-linear dynamical systems, based on incomplete and noisy measurements, are considered using Monte Carlo simulations. Given the measurements. the proposed method obtains the marginalized posterior distribution of an appropriately chosen (ideally small) subset of the state vector using a particle filter. Samples (particles) of the marginalized states are then used to construct a family of conditionally linearized system of equations and thus obtain the posterior distribution of the states using a bank of Kalman filters. Discrete process equations for the marginalized states are derived through truncated Ito-Taylor expansions. Increased analyticity and reduced dispersion of weights computed over a smaller sample space of marginalized states are the key features of the filter that help achieve smaller sample variance of the estimates. Numerical illustrations are provided for state/parameter estimations of a Duffing oscillator and a 3-DOF non-linear oscillator. Performance of the filter in parameter estimation is also assessed using measurements obtained through experiments on simple models in the laboratory. Despite an added computational cost, the results verify that the proposed filter generally produces estimates with lower sample variance over the standard sequential importance sampling (SIS) filter.

Radiation induced depinning of a charge density wave system

The frequency-dependent response of a pinned charge density wave is considered in terms of forced vibration of an oscillator held in an anharmonic well. It is shown that the effective pinning-frequency can be reduced by applying a d.c. field. If a strong a.c. field, superposed on a d.c. field is applied on such a system “jumps” can be observed in the frequency dependent response of the system. The conditions at which these “jumps” occur are investigated with reference to NbSe3. The possibility of observing such phenomena in other systems like superionic conductors, non-linear dielectrics like ferroelectrics is pointed out. The characteristics are expressed in terms of some “scaled variables” — in terms of which the characteristics show a universal behaviour.

Approximate analysis of coupled non-linear non-conservative systems subjected to step-function excitation

The paper deals with the approximate analysis of non-linear non-conservative systems oftwo degrees of freedom subjected to step-function excitation. The method of averaging of Krylov and Bogoliubov is used to arrive at the approximate equations for amplitude and phase. An example of a spring-mass-damper system is presented to illustrate the method and a comparison with numerical results brings out the validity of the approach.

Paraxial Maxwell beams: transformation by general linear optical systems

Extending the work of earlier papers on the relativistic-front description of paraxial optics and the formulation of Fourier optics for vector waves consistent with the Maxwell equations, we generalize the Jones calculus of axial plane waves to describe the action of the most general linear optical system on paraxial Maxwell fields. Several examples are worked out, and in each case it is shown that the formalism leads to physically correct results. The importance of retaining the small components of the field vectors along the axis of the system for a consistent description is emphasized.