Theory of radar imaging using time-frequency distribution

The concept of radar was developed for the estimation of the distance (range) and velocity of a target from a receiver. The distance measurement is obtained by measuring the time taken for the transmitted signal to propagate to the target and return to the receiver. The target's velocity is determined by measuring the Doppler induced frequency shift of the returned signal caused by the rate of change of the time- delay from the target. As researchers further developed conventional radar systems it become apparent that additional information was contained in the backscattered signal and that this information could in fact be used to describe the shape of the target itself. It is due to the fact that a target can be considered to be a collection of individual point scatterers, each of which has its own velocity and time- delay. DelayDoppler parameter estimation of each of these point scatterers thus corresponds to a mapping of the target's range and cross range, thus producing an image of the target. Much research has been done in this area since the early radar imaging work of the 1960s. At present there are two main categories into which radar imaging falls. The first of these is related to the case where the backscattered signal is considered to be deterministic. The second is related to the case where the backscattered signal is of a stochastic nature. In both cases the information which describes the target's scattering function is extracted by the use of the ambiguity function, a function which correlates the backscattered signal in time and frequency with the transmitted signal. In practical situations, it is often necessary to have the transmitter and the receiver of the radar system sited at different locations. The problem in these situations is 'that a reference signal must then be present in order to calculate the ambiguity function. This causes an additional problem in that detailed phase information about the transmitted signal is then required at the receiver. It is this latter problem which has led to the investigation of radar imaging using time- frequency distributions. As will be shown in this thesis, the phase information about the transmitted signal can be extracted from the backscattered signal using time- frequency distributions. The principle aim of this thesis was in the development, and subsequent discussion into the theory of radar imaging, using time- frequency distributions. Consideration is first given to the case where the target is diffuse, ie. where the backscattered signal has temporal stationarity and a spatially white power spectral density. The complementary situation is also investigated, ie. where the target is no longer diffuse, but some degree of correlation exists between the time- frequency points. Computer simulations are presented to demonstrate the concepts and theories developed in the thesis. For the proposed radar system to be practically realisable, both the time- frequency distributions and the associated algorithms developed must be able to be implemented in a timely manner. For this reason an optical architecture is proposed. This architecture is specifically designed to obtain the required time and frequency resolution when using laser radar imaging. The complex light amplitude distributions produced by this architecture have been computer simulated using an optical compiler.

Nonlinear Dynamics of Beams with Stochastic Parameter Variations

The aim of this paper is to investigate the steady state response of beams under the action of random support motions. The study is of relevance in the context of earthquake response of extended land based structures such as pipelines and long span bridges, and, secondary systems such as piping networks in nuclear power plant installations. The following complicating features are accounted for in the response analysis: (a) differential support motions: this is characterized in terms of cross power spectral density functions associated with distinct support motions, (b) nonlinear support conditions, and (c) stochastically inhomogeneous stiffness and mass variations of the beam structure; questions on non-Gaussian models for these variations are considered. The method of stochastic finite elements is combined with equivalent linearization technique and Monte Carlo simulations to obtain response moments.

The relationship between the radio and gamma-ray emission of blazars

Blazars are active galaxies with a jet closely oriented to our line of sight. They are powerful, variable emitters from radio to gamma-ray wavelengths. Although the general picture of synchrotron emission at low energies and inverse Compton at high energies is well established, important aspects of blazars are not well understood. In particular, the location of the gamma-ray emission region is not clearly established, with some theories favoring a location close to the central engine, while others place it at parsec scales in the radio jet.

We developed a program to locate the gamma-ray emission site in blazars, through the study of correlated variations between their gamma-ray and radio-wave emission. Correlated variations are expected when there is a relation between emission processes at both bands, while delays tell us about the relative location of their energy generation zones. Monitoring at 15 GHz using the Owens Valley Radio Observatory 40 meter telescope started in mid-2007. The program monitors 1593 blazars twice per week, including all blazars detected by the Fermi Gamma-ray Space Telescope (Fermi) north of -20 degrees declination. This program complements the continuous monitoring of gamma-rays by Fermi.

Three year long gamma-ray light curves for bright Fermi blazars are cross-correlated with four years of radio monitoring. The significance of cross-correlation peaks is investigated using simulations that account for the uneven sampling and noise properties of the light curves, which are modeled as red-noise processes with a simple power-law power spectral density. We found that out of 86 sources with high quality data, only three show significant correlations (AO 0235+164, B2 2308+34 and PKS 1502+106). Additionally, we find a significant correlation for Mrk 421 when including the strong gamma-ray/radio flare of late 2012. In all four cases radio variations lag gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. For PKS 1502+106 we locate the gamma-ray emission site parsecs away from the central engine, thus disfavoring the model of Blandford and Levinson (1995), while other cases are inconclusive. These findings show that continuous monitoring over long time periods is required to understand the cross-correlation between gamma-ray and radio-wave variability in most blazars.

A mathematical and experimental investigation of microcycle spectral estimates of semiconductor flicker noise

The experimental portion of this thesis tries to estimate the density of the power spectrum of very low frequency semiconductor noise, from 10^-6.3 cps to 1. cps with a greater accuracy than that achieved in previous similar attempts: it is concluded that the spectrum is 1/f^α with α approximately 1.3 over most of the frequency range, but appearing to have a value of about 1 in the lowest decade. The noise sources are, among others, the first stage circuits of a grounded input silicon epitaxial operational amplifier. This thesis also investigates a peculiar form of stationarity which seems to distinguish flicker noise from other semiconductor noise.

In order to decrease by an order of magnitude the pernicious effects of temperature drifts, semiconductor "aging", and possible mechanical failures associated with prolonged periods of data taking, 10 independent noise sources were time-multiplexed and their spectral estimates were subsequently averaged. If the sources have similar spectra, it is demonstrated that this reduces the necessary data-taking time by a factor of 10 for a given accuracy.

In view of the measured high temperature sensitivity of the noise sources, it was necessary to combine the passive attenuation of a special-material container with active control. The noise sources were placed in a copper-epoxy container of high heat capacity and medium heat conductivity, and that container was immersed in a temperature controlled circulating ethylene-glycol bath.

Other spectra of interest, estimated from data taken concurrently with the semiconductor noise data were the spectra of the bath's controlled temperature, the semiconductor surface temperature, and the power supply voltage amplitude fluctuations. A brief description of the equipment constructed to obtain the aforementioned data is included.

The analytical portion of this work is concerned with the following questions: what is the best final spectral density estimate given 10 statistically independent ones of varying quality and magnitude? How can the Blackman and Tukey algorithm which is used for spectral estimation in this work be improved upon? How can non-equidistant sampling reduce data processing cost? Should one try to remove common trands shared by supposedly statistically independent noise sources and, if so, what are the mathematical difficulties involved? What is a physically plausible mathematical model that can account for flicker noise and what are the mathematical implications on its statistical properties? Finally, the variance of the spectral estimate obtained through the Blackman/Tukey algorithm is analyzed in greater detail; the variance is shown to diverge for α ≥ 1 in an assumed power spectrum of k/|f|^α, unless the assumed spectrum is "truncated".

Synthesis and characterization of photo-cross-linked hydrogels based on biodegradable polyphosphoesters and poly(ethylene glycol) copolymers

Novel biodegradable hydrogels by photo-cross-linking macromers based on polyphosphoesters and poly(ethylene glycol) (PEG) are reported. Photo-cross-linkable macromers were synthesized by ring-opening polymerization of the cyclic phosphoester monomer 2-(2-oxo-1,3,2-dioxaphospholoyloxy) ethyl methacrylate (OPEMA) using PEG as the initiator and stannous octoate as the catalyst. The macrorners were characterized by H-1 NMR, Fourier transform infrared spectroscopy, and gel permeation chromatography measurements. The content of polyphosphoester in the macromer was controlled by varying the feed ratio of OPEMA to PEG. Hydrogels were fabricated by exposing aqueous solutions of macromers with 0.05% (w/w) photoinitiator to UV light irradiation, and their swelling kinetics as well as degradation behaviors were evaluated. The results demonstrated that cross-linking density and pH values strongly affected the degradation rates. The macromers was compatible to osteoblast cells, not exhibiting significant cytotoxicity up to 0.5 mg/mL. "Live/dead" cell staining assay also demonstrated that a large majority of the osteoblast cells remained viable after encapsulation into the hydrogel constructs, showing their potential as tissue engineering scaffolds.

Effect of cross channel salinity imbalance due to freshwater influence in a narrow estuary

Senior thesis written for Oceanography 445

Adaptive modeling and high quality spectral estimation for speech enhancement

In this work an adaptive modeling and spectral estimation scheme based on a dual Discrete Kalman Filtering (DKF) is proposed for speech enhancement. Both speech and noise signals are modeled by an autoregressive structure which provides an underlying time frame dependency and improves time-frequency resolution. The model parameters are arranged to obtain a combined state-space model and are also used to calculate instantaneous power spectral density estimates. The speech enhancement is performed by a dual discrete Kalman filter that simultaneously gives estimates for the models and the signals. This approach is particularly useful as a pre-processing module for parametric based speech recognition systems that rely on spectral time dependent models. The system performance has been evaluated by a set of human listeners and by spectral distances. In both cases the use of this pre-processing module has led to improved results.

GLS Detrending, Efficient Unit Root Tests and Structural Change

We extend the class of M-tests for a unit root analyzed by Perron and Ng (1996) and Ng and Perron (1997) to the case where a change in the trend function is allowed to occur at an unknown time. These tests M(GLS) adopt the GLS detrending approach of Dufour and King (1991) and Elliott, Rothenberg and Stock (1996) (ERS). Following Perron (1989), we consider two models : one allowing for a change in slope and the other for both a change in intercept and slope. We derive the asymptotic distribution of the tests as well as that of the feasible point optimal tests PT(GLS) suggested by ERS. The asymptotic critical values of the tests are tabulated. Also, we compute the non-centrality parameter used for the local GLS detrending that permits the tests to have 50% asymptotic power at that value. We show that the M(GLS) and PT(GLS) tests have an asymptotic power function close to the power envelope. An extensive simulation study analyzes the size and power in finite samples under various methods to select the truncation lag for the autoregressive spectral density estimator. An empirical application is also provided.

NMR study of slow motions of HDA hydrocarbons chains inside lamellar structures

Measurements of H-1 and C-13 Nuclear Magnetic Resonance (NMR) for the nano-composite materials formed by the intercalation of hexadecylamine (HDA) in metal oxides (TiO2, V2O5 and MoO3), are reported. The H-1 NMR spin-lattice relaxation in the rotating frame was described by using the spectral density due to Davidson and Cole, which incorporates a distribution of correlation times characterized by a width parameter epsilon. The fitting of the data was obtained for epsilon = 0.74, indicating that the correlation times are distributed over a narrow range in this system. High-resolution C-13 NMR techniques were used to resolve the NMR lines of middle-chain methylene groups in the spectra and variable contact time cross-polarization {H-1-}C-13 experiments were employed to analyze the reorientation dynamics of the CH3 and CH2 groups in the HDA chains.

Spectral properties of temporally aggregated long memory processes

This paper derives the spectral density function of aggregated long memory processes in light of the aliasing effect. The results are different from previous analyses in the literature and a small simulation exercise provides evidence in our favour. The main result point to that flow aggregates from long memory processes shall be less biased than stock ones, although both retain the degree of long memory. This result is illustrated with the daily US Dollar/ French Franc exchange rate series.

Periodic breathing during ascent to extreme altitude quantified by spectral analysis of the respiratory volume signal

High altitude periodic breathing (PB) shares some common pathophysiologic aspects with sleep apnea, Cheyne-Stokes respiration and PB in heart failure patients. Methods that allow quantifying instabilities of respiratory control provide valuable insights in physiologic mechanisms and help to identify therapeutic targets. Under the hypothesis that high altitude PB appears even during physical activity and can be identified in comparison to visual analysis in conditions of low SNR, this study aims to identify PB by characterizing the respiratory pattern through the respiratory volume signal. A number of spectral parameters are extracted from the power spectral density (PSD) of the volume signal, derived from respiratory inductive plethysmography and evaluated through a linear discriminant analysis. A dataset of 34 healthy mountaineers ascending to Mt. Muztagh Ata, China (7,546 m) visually labeled as PB and non periodic breathing (nPB) is analyzed. All climbing periods within all the ascents are considered (total climbing periods: 371 nPB and 40 PB). The best crossvalidated result classifying PB and nPB is obtained with Pm (power of the modulation frequency band) and R (ratio between modulation and respiration power) with an accuracy of 80.3% and area under the receiver operating characteristic curve of 84.5%. Comparing the subjects from 1(st) and 2(nd) ascents (at the same altitudes but the latter more acclimatized) the effect of acclimatization is evaluated. SaO(2) and periodic breathing cycles significantly increased with acclimatization (p-value < 0.05). Higher Pm and higher respiratory frequencies are observed at lower SaO(2), through a significant negative correlation (p-value < 0.01). Higher Pm is observed at climbing periods visually labeled as PB with > 5 periodic breathing cycles through a significant positive correlation (p-value < 0.01). Our data demonstrate that quantification of the respiratory volume signal using spectral analysis is suitable to identify effects of hypobaric hypoxia on control of breathing.

Alkenone analyses of sediment cores from the western Arabian Sea

Records of total organic carbon (TOC) and C37 alkenones were used as indicators for past primary productivity in the western and eastern Arabian Sea. Data from GeoB 3005, an open ocean site in the western Arabian Sea upwelling area, are compared with similar records of GeoB 3007 from the Owen Ridge, Ocean Drilling Program (ODP) Site 723 from the continental margin off Oman and MD 900963 from the eastern Arabian Sea. TOC/C37 alkenone records together with other proxies used to reconstruct upwelling intensity, indicate periods of high productivity in tune with precessional forcing all over the Arabian Sea. Based on their phase-relationship to variations in boreal summer insolation they can be divided into three groups. In the western Arabian Sea the precession-related phasing is different between productivity proxies and those for summer monsoon wind strength and upwelling intensity. TOC and C37 alkenone records from the western Arabian Sea lag the other monsoonal indicators by about 5 kyr, but lead productivity indicators from the eastern Arabian Sea by 3 kyr. Based on the differences in phase relationships associated with the precessional cycling between productivity and monsoonal proxies in the western Arabian Sea it is proposed that the TOC/C37 alkenone signal in the western Arabian Sea document a combined signal of moderate SW monsoon winds and of strengthened and prolonged NE monsoon winds. In the eastern Arabian Sea the phasing hints to coincidence between maximum productivity and stronger NE monsoon winds associated with precession-related maxima in ice volume. In contrast, variations in paleoproductivity at site GeoB 3007 from the Owen Ridge indicate productivity maxima during glacial substages 8.2, 6.2 and 2.2, whereas precessionrelated changes are of only minor importance at this location. The results of frequency analyses confirm that productivity at site GeoB 3007 responds predominantly to glacialinterglacial climate changes, while site GeoB 3005 from the open ocean upwelling region near the Gulf of Aden is dominated by precessional insolation. A possible explanation for the pattern revealed at the Owen Ridge is the periodic NW-SE displacement of the Findlater Jet axis, which separates the region of open ocean upwelling to the northwest from downwelling to the southeast ofthe jet. The carbon isotopes of planktic foraminifera reflect nutrient related d13C variations of dissolved inorganic carbon. The difference between the planktic foraminifera Globigerinoides ruber (w), living in the upper 50 m of the water column, and the deeper Iiving Neogloboquadrina dutertrei (Delta d13Cr-d) of core GeoB 3005 displays nutrient variations in the upwelling area near the Gulf of Aden. The results of cross-spectral analyses between Deltad13Cr-d of GeoB 3005 and proxies for SW monsoon intensity indicate, too, a dissociation of productivity from monsoonal upwelling intensity. Instead, productivity depends mainly on the availability of nutrients, while upwelling intensity of sub-surface water masses seems to be of only secondary importance. Additionally, sea surface temperatures (SSTs) were reconstructed using the unsaturation ratio of C37 alkenones. Alkenone SSTs reflect annual mean temperatures rather than explicitly the season of upwelling. This is evident from alkenone SSTs in a transect of surface sediments extending from the inner Gulf of Aden into the western Arabian Sea. The alkenone-derived SST records of GeoB 3005 from the open ocean upwelling region near the Gulf of Aden and GeoB 3007 from the Owen Ridge reveal similar variations with high SSTs during interglacial and low SSTs during glacial periods. The glacial marine oxygen isotope stage (MIS) 6 remains relatively warm and was not as cold as MIS 3 to 4 and 8 according to the alkenone SST. Similar variation-patterns were reconstructed in the coastal upwelling area off Oman for ODP Site 723 as weIl as in the eastern Arabian Sea for MD 900963, where upwelling is not as pronounced as in the western Arabian Sea. Spectral-analyses indicate that SST changes are in good agreement with the modulation of low-latitude precessional insolation changes by eccentricity. However, they do not show the pronounced cydicity in the precessional frequency band, which is characteristic for variations in paleoproductivity. Although the overall variation pattern is very similar, a dose comparison between the western (GeoB 3005) and the eastern Arabian Sea (MD 900963) shows larger differences between both sites during cold intervals than during periods of warm SSTs. This is attributed to a more effective cooling of surface waters in the western Arabian Sea by prolonged NE monsoon winds during times of expanded Northern Hemisphere ice-sheets, thereby lowering the annual mean SSTs stronger than in the eastern Arabian Sea.

Vortex-induced vibrations on bridges : the interaction mechanism in complex sections and a new proposed semi-empirical model

El presente trabajo de investigación se ocupa del estudio de las vibraciones verticales inducidas por vórtices (VIV) en aquellos puentes que, por sus características geométricas y propiedades dinámicas, muestran cierta sensibilidad este tipo de fenómeno aeroelástico. El objeto principal es el análisis del mecanismo de interacción viento-estructura sobre secciones no fuseladas de geometría simple, con objeto de realizar una adecuada caracterización del problema y poder abordar posteriormente el análisis de otras secciones de geometría más compleja, representativas de los principales elementos estructurales de los puentes, como arcos, tableros, torres y pilas. Este aspecto es fundamental durante la fase de diseño del puente, donde deberán tenerse en cuenta también una serie de detalles que pueden influir significativamente su sensibilidad ante problemas aerodinámicos, como la morfología y dimensiones principales de la sección transversal del tablero, la disposición de barreras de seguridad y barreras cortaviento, o las riostras que unen diferentes elementos estructurales. La configuración de dos elementos en tándem o la construcción de un puente en las inmediaciones de otro existente son otros aspectos a considerar respecto a la sensibilidad frente a efectos aeroelásticos. El estudio se ha llevado a cabo principalmente mediante la implementación de simulaciones numéricas que reproducen la interacción entre la corriente de aire y secciones representativas de modelos estructurales, a partir de un código CFD basado en el método de las partículas de vórtices (VPM), siguiendo por tanto un esquema Lagrangiano. Los resultados han sido validados con datos experimentales existentes, valores procedentes de ensayos en túnel de viento y registros reales a partir de diferentes casos de estudio: Alconétar (2006), Niterói (1980), Trans- Tokyo Bay (1995) y Volgogrado (2010). Finalmente, se propone un modelo semi-empírico para la estimación del rango de velocidades críticas y amplitudes de oscilación basado en la utilización de las derivadas de flameo de Scanlan, y la densidad espectral de las fuerzas aerodinámicas en el dominio de la frecuencia. The present research work concerns the study of vertical vortex-induced vibrations (VIV) in bridges which show certain sensitivity to this type of aeroelastic phenomenon. It focuses on the analysis of the wind-structure interaction mechanism on bluff sections, with the objective of making a good characterisation of the problem and subsequently addressing the analysis of sections with a complex geometry, which are representative of the bridge structural elements, such as arches, decks, towers and piers. This issue is of relative importance during the bridge design phase, since minor details of the aforementioned elements can significantly influence its sensitivity to aerodynamic problems. The shape and main dimensions of the deck cross section, the addition of safety barriers and windshields, the presence of braces to enhance the structure mechanical properties, the utilisation of cross sections in tandem arrangement, or the erection of a new bridge in the vicinity of another existing one are some of the aspects to be considered regarding the sensitivity to the aeroelastic effects. The study has been carried out mainly through the implementation of numerical simulations that reproduces the interaction between the airflow and the representative cross section of a structural bridge model, by the use of a CFD code based on the vortex particle method (VPM), thus following a Lagrangian scheme. The results have been validated with existing experimental data, values from wind tunnel tests and full scale observations from the different case studies: Alconétar (2006), Niterói (1980), Trans-Tokyo Bay (1995) and Volgograd (2010). Finally, a new semi-empirical model is proposed for the estimation of the critical wind velocity ranges and oscillation amplitudes based on the use of the Scanlan’s flutter derivatives and the power spectral density of aerodynamic force time history in the frequency domain.

Integration of autonomic and local mechanisms in regulating cardiovascular responses to heating and cooling in a reptile (Crocodylus porosus)

Reptiles change heart rate and blood flow patterns in response to heating and cooling, thereby decreasing the behavioural cost of thermoregulation. We tested the hypothesis that locally produced vasoactive substances, nitric oxide and prostaglandins, mediate the cardiovascular response of reptiles to heat. Heart rate and blood pressure were measured in eight crocodiles (Crocodylus porosus) during heating and cooling and while sequentially inhibiting nitric-oxide synthase and cyclooxygenase enzymes. Heart rate and blood pressure were significantly higher during heating than during cooling in all treatments. Power spectral density of heart rate and blood pressure increased significantly during heating and cooling compared to the preceding period of thermal equilibrium. Spectral density of heart rate in the high frequency band (0.19-0.70 Hz) was significantly greater during cooling in the saline treatment compared to when nitric-oxide synthase and cyclooxygenase enzymes were inhibited. Cross spectral analysis showed that changes in blood pressure preceded heart rate changes at low frequencies (

Innovative measurement methodology in the study of weathering of automotive coatings

An initial aim of this project was to evaluate the conventional techniques used in the analysis of newly prepared environmentally friendly water-borne automotive coatings and compare them with solvent-borne coatings having comparable formulations. The investigation was carried out on microtuned layers as well as on complete automotive multi-layer paint systems. Methods used included the very traditional methods of gloss and hardness and the commonly used photo-oxidation index (from FTIR spectral analysis). All methods enabled the durability to weathering of the automotive coatings to be initially investigated. However, a primary aim of this work was to develop methods for analysing the early stages of chemical and property changes in both the solvent-borne and water-borne coating systems that take place during outdoor natural weathering exposures and under accelerated artificial exposures. This was achieved by using dynamic mechanical analysis (DMA), in both tension mode on the microtomed films (on all depths of the coating systems from the uppermost clear-coat right down to the electron-coat) and bending mode of the full (unmicrotomed) systems, as well as MALDI-Tof analysis on the movement of the stabilisers in the full systems. Changes in glass transition temperature and relative cross-link density were determined after weathering and these were related to changes in the chemistries of the binder systems of the coatings after weathering. Concentration profiles of the UV-stabilisers (UVA and HALS) in the coating systems were analysed as a consequence of migration in the coating systems in separate microtomed layers of the paint samples (depth profiling) after weathering and diffusion co-efficient and solubility parameters were determined for the UV stabilisers in the coating systems. The methods developed were used to determine the various physical and chemical changes that take place during weathering of the different (water-borne and solvent-borne) systems (photoxidation). The solvent-borne formulations showed less changes after weathering (both natural and accelerated) than the corresponding water-borne formulations due to the lower level of cross-links in the binders of the water-borne systems. The silver systems examined were more durable than the blue systems due to the reflecting power of the aluminium and the lower temperature of the silver coatings.