A new algorithm to accelerate harmonic analysis of time series

The Lomb periodogram has been traditionally a tool that allows us to elucidate if a frequency turns out to be important for explaining the behaviour of a given time series. Many linear and nonlinear reiterative harmonic processes that are used for studying the spectral content of a time series take into account this periodogram in order to avoid including spurious frequencies in their models due to the leakage problem of energy from one frequency to others. However, the estimation of the periodogram requires long computation time that makes the harmonic analysis slower when we deal with certain time series. Here we propose an algorithm that accelerates the extraction of the most remarkable frequencies from the periodogram, avoiding its whole estimation of the harmonic process at each iteration. This algorithm allows the user to perform a specific analysis of a given scalar time series. As a result, we obtain a functional model made of (1) a trend component, (2) a linear combination of Fourier terms, and (3) the so-called mixed secular terms by reducing the computation time of the estimation of the periodogram.

Spectral analysis of seabed elevation from the Danish Wadden Sea

Bedforms both reflect and influence shallow water hydrodynamics and sediment dynamics. A correct characterization of their spatial distribution and dimensions is required for the understanding, assessment and prediction of numerous coastal processes. A method to parameterize geometrical characteristics using two-dimensional (2D) spectral analysis is presented and tested on seabed elevation data from the Knudedyb tidal inlet in the Danish Wadden Sea, where large compound bedforms are found. The bathymetric data were divided into 20x20 m areas on which a 2D spectral analysis was applied. The most energetic peak of the 2D spectrum was found and its energy, frequency and direction were calculated. A power-law was fitted to the average of slices taken through the 2D spectrum; its slope and y-intercept were calculated. Using these results the test area was morphologically classified into 4 distinct morphological regions. The most energetic peak and the slope and intercept of the power-law showed high values above the crest of the primary bedforms and scour holes, low values in areas without bedforms, and intermediate values in areas with secondary bedforms. The secondary bedform dimensions and orientations were calculated. An area of 700x700 m was used to determine the characteristics of the primary bedforms. However, they were less distinctively characterized compared to the secondary bedforms due to relatively large variations in their orientations and wavelengths. The method is thus appropriate for morphological classification of the seabed and for bedform characterization, being most efficient in areas characterized by bedforms with regular dimensions and directions.

Statistical properties of fast Fourier transform coefficients computed from real-valued, covariance-stationary, periodic random sequences /

"July 1976."

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Character recognition systems using Fourier transformation in incoherent light

Many attempts have been made to overcome problems involved in character recognition which have resulted in the manufacture of character reading machines. An investigation into a new approach to character recognition is described. Features for recognition are Fourier coefficients. These are generated optically by convolving characters with periodic gratings. The development of hardware to enable automatic measurement of contrast and position of periodic shadows produced by the convolution is described. Fourier coefficients of character sets were measured, many of which are tabulated. Their analysis revealed that a few low frequency sampling points could be selected to recognise sets of numerals. Limited treatment is given to show the effect of type face variations on the values of coefficients which culminated in the location of six sampling frequencies used as features to recognise numerals in two type fonts. Finally, the construction of two character recognition machines is compared and contrasted. The first is a pilot plant based on a test bed optical Fourier analyser, while the second is a more streamlined machine d(3signed for high speed reading. Reasons to indicate that the latter machine would be the most suitable to adapt for industrial and commercial applications are discussed.

Quantification of hydroxycinnamic acids and lignin in perennial forage and energy grasses by Fourier-transform infrared spectroscopy and partial least squares regression

Levels of lignin and hydroxycinnamic acid wall components in three genera of forage grasses (Lolium,Festuca and Dactylis) have been accurately predicted by Fourier-transform infrared spectroscopy using partial least squares models correlated to analytical measurements. Different models were derived that predicted the concentrations of acid detergent lignin, total hydroxycinnamic acids, total ferulate monomers plus dimers, p-coumarate and ferulate dimers in independent spectral test data from methanol extracted samples of perennial forage grass with accuracies of 92.8%, 86.5%, 86.1%, 59.7% and 84.7% respectively, and analysis of model projection scores showed that the models relied generally on spectral features that are known absorptions of these compounds. Acid detergent lignin was predicted in samples of two species of energy grass, (Phalaris arundinacea and Pancium virgatum) with an accuracy of 84.5%.

The analysis of organic matter in oil-shales

The organic matter in five oil shales (three from the Kimmeridge Clay sequence, one from the Oxford Clay sequence and one from the Julia Creek deposits in Australia) has been isolated by acid demineralisation, separated into kerogens and bitumens by solvent extraction and then characterised in some detail by chromatographic, spectroscopic and degradative techniques. Kerogens cannot be characterised as easily as bitumens because of their insolubility, and hence before any detailed molecular information can be obtained from them they must be degraded into lower molecular weight, more soluble components. Unfortunately, the determination of kerogen structures has all too often involved degradations that were far too harsh and which lead to destruction of much of the structural information. For this reason a number of milder more selective degradative procedures have been tested and used to probe the structure of kerogens. These are: 1. Lithium aluminium hydride reduction. - This procedure is commonly used to remove pyrite from kerogens and it may also increase their solubility by reduction of labile functional groups. Although reduction of the kerogens was confirmed, increases in solubility were correlated with pyrite content and not kerogen reduction. 2. O-methylation in the presence of a phase transfer catalyst. - By the removal of hydrogen bond interactions via O-methylation, it was possible to determine the contribution of such secondary interactions to the insolubility of the kerogens. Problems were encountered with the use of the phase transfer catalyst. 3. Stepwise alkaline potassium permanganate oxidation. - Significant kerogen dissolution was achieved using this procedure but uncontrolled oxidation of initial oxidation products proved to be a problem. A comparison with the peroxytrifluoroaceticacid oxidation of these kerogens was made. 4. Peroxytrifluoroacetic acid oxidation. - This was used because it preferentially degrades aromatic rings whilst leaving any benzylic positions intact. Considerable conversion of the kerogens into soluble products was achieved with this procedure. At all stages of degradation the products were fully characterised where possible using a variety of techniques including elemental analysis, solution state 1H and 13C nuclear magnetic resonance, solid state 13C nuclear magnetic resonance, gel-permeationchromatography, gas chromatography-mass spectroscopy, fourier transform infra-red spectroscopy and some ultra violet-visible spectroscopy.

Orbit and altimetric corrections for the ERS satellites through analysis of single and dual satellite crossovers

Due to the failure of PRARE the orbital accuracy of ERS-1 is typically 10-15 cm radially as compared to 3-4cm for TOPEX/Poseidon. To gain the most from these simultaneous datasets it is necessary to improve the orbital accuracy of ERS-1 so that it is commensurate with that of TOPEX/Poseidon. For the integration of these two datasets it is also necessary to determine the altimeter and sea state biases for each of the satellites. Several models for the sea state bias of ERS-1 are considered by analysis of the ERS-1 single satellite crossovers. The model adopted consists of the sea state bias as a percentage of the significant wave height, namely 5.95%. The removal of ERS-1 orbit error and recovery of an ERS-1 - TOPEX/Poseidon relative bias are both achieved by analysis of dual crossover residuals. The gravitational field based radial orbit error is modelled by a finite Fourier expansion series with the dominant frequencies determined by analysis of the JGM-2 co-variance matrix. Periodic and secular terms to model the errors due to atmospheric density, solar radiation pressure and initial state vector mis-modelling are also solved for. Validation of the dataset unification consists of comparing the mean sea surface topographies and annual variabilities derived from both the corrected and uncorrected ERS-1 orbits with those derived from TOPEX/Poseidon. The global and regional geographically fixed/variable orbit errors are also analysed pre and post correction, and a significant reduction is noted. Finally the use of dual/single satellite crossovers and repeat pass data, for the calibration of ERS-2 with respect to ERS-1 and TOPEX/Poseidon is shown by calculating the ERS-1/2 sea state and relative biases.

Numerical and statistical analysis of long-haul undersea optical communication systems

Firstly, we numerically model a practical 20 Gb/s undersea configuration employing the Return-to-Zero Differential Phase Shift Keying data format. The modelling is completed using the Split-Step Fourier Method to solve the Generalised Nonlinear Schrdinger Equation. We optimise the dispersion map and per-channel launch power of these channels and investigate how the choice of pre/post compensation can influence the performance. After obtaining these optimal configurations, we investigate the Bit Error Rate estimation of these systems and we see that estimation based on Gaussian electrical current systems is appropriate for systems of this type, indicating quasi-linear behaviour. The introduction of narrower pulses due to the deployment of quasi-linear transmission decreases the tolerance to chromatic dispersion and intra-channel nonlinearity. We used tools from Mathematical Statistics to study the behaviour of these channels in order to develop new methods to estimate Bit Error Rate. In the final section, we consider the estimation of Eye Closure Penalty, a popular measure of signal distortion. Using a numerical example and assuming the symmetry of eye closure, we see that we can simply estimate Eye Closure Penalty using Gaussian statistics. We also see that the statistics of the logical ones dominates the statistics of the logical ones dominates the statistics of signal distortion in the case of Return-to-Zero On-Off Keying configurations.

Proteomic analysis of a noninvasive human model of acute inflammation and its resolution:the twenty-one day gingivitis model

The 21-day experimental gingivitis model, an established noninvasive model of inflammation in response to increasing bacterial accumulation in humans, is designed to enable the study of both the induction and resolution of inflammation. Here, we have analyzed gingival crevicular fluid, an oral fluid comprising a serum transudate and tissue exudates, by LC-MS/MS using Fourier transform ion cyclotron resonance mass spectrometry and iTRAQ isobaric mass tags, to establish meta-proteomic profiles of inflammation-induced changes in proteins in healthy young volunteers. Across the course of experimentally induced gingivitis, we identified 16 bacterial and 186 human proteins. Although abundances of the bacterial proteins identified did not vary temporally, Fusobacterium outer membrane proteins were detected. Fusobacterium species have previously been associated with periodontal health or disease. The human proteins identified spanned a wide range of compartments (both extracellular and intracellular) and functions, including serum proteins, proteins displaying antibacterial properties, and proteins with functions associated with cellular transcription, DNA binding, the cytoskeleton, cell adhesion, and cilia. PolySNAP3 clustering software was used in a multilayered analytical approach. Clusters of proteins that associated with changes to the clinical parameters included neuronal and synapse associated proteins.

A Generalized Convolution with a Weight Function for the Fourier Cosine and Sine Transforms

A generalized convolution with a weight function for the Fourier cosine and sine transforms is introduced. Its properties and applications to solving a system of integral equations are considered.

On the Range of the Fourier Transform Associated with the Spherical Mean Operator

We characterize the range of some spaces of functions by the Fourier transform associated with the spherical mean operator R and we give a new description of the Schwartz spaces. Next, we prove a Paley-Wiener and a Paley-Wiener-Schawrtz theorems.

Dynamic analysis of submerged microscale plates:the effects of acoustic radiation and viscous dissipation

The aim of this paper is to study the dynamic characteristics of micromechanical rectangular plates used as sensing elements in a viscous compressible fluid. A novel modelling procedure for the plate- fluid interaction problem is developed on the basis of linearized Navier-Stokes equations and noslip conditions. Analytical expression for the fluidloading impedance is obtained using a double Fourier transform approach. This modelling work provides us an analytical means to study the effects of inertial loading, acoustic radiation and viscous dissipation of the fluid acting on the vibration of microplates. The numerical simulation is conducted on microplates with different boundary conditions and fluids with different viscosities. The simulation results reveal that the acoustic radiation dominates the damping mechanism of the submerged microplates. It is also proved that microplates offer better sensitivities (Q-factors) than the conventional beam type microcantilevers beingmass sensing platforms in a viscous fluid environment. The frequency response features of microplates under highly viscous fluid loading are studied using the present model. The dynamics of the microplates with all edges clamped are less influenced by the highly viscous dissipation of the fluid than the microplates with other types of boundary conditions.

Dynamically mechanical and nano-impact (fatigue) analysis of touch screen thin films deposited on polyethylene terephthalate substrate

Nano-scale touch screen thin film have not been thoroughly investigated in terms of dynamic impact analysis under various strain rates. This research is focused on two different thin films, Zinc Oxide (ZnO) film and Indium Tin Oxide (ITO) film, deposited on Polyethylene Terephthalate (PET) substrate for the standard touch screen panels. Dynamic Mechanical Analysis (DMA) was performed on the ZnO film coated PET substrates. Nano-impact (fatigue) testing was performed on ITO film coated PET substrates. Other analysis includes hardness and the elastic modulus measurements, atomic force microscopy (AFM), Fourier Transform Infrared Spectroscopy (FTIR) and the Scanning Electron Microscopy (SEM) of the film surface.
Ten delta of DMA is described as the ratio of loss modulus (viscous properties) and storage modulus (elastic properties) of the material and its peak against time identifies the glass transition temperature (Tg). Thus, in essence the Tg recognizes changes from glassy to rubber state of the material and for our sample ZnO film, Tg was found as 388.3 K. The DMA results also showed that the Ten delta curve for Tg increases monotonically in the viscoelastic state (before Tg) and decreases sharply in the rubber state (after Tg) until recrystallization of ZnO takes place. This led to an interpretation that enhanced ductility can be achieved by negating the strength of the material.
For the nano-impact testing using the ITO coated PET, the damage started with the crack initiation and propagation. The interpretation of the nano-impact results depended on the characteristics of the loading history. Under the nano-impact loading, the surface structure of ITO film suffered from several forms of failure damages that range from deformation to catastrophic failures. It is concluded that in such type of application, the films should have low residual stress to prevent deformation, good adhesive strength, durable and good resistance to wear.