Subspace system identification framework for the analysis of multimoded propagation of THz-transient signals

We provide a system identification framework for the analysis of THz-transient data. The subspace identification algorithm for both deterministic and stochastic systems is used to model the time-domain responses of structures under broadband excitation. Structures with additional time delays can be modelled within the state-space framework using additional state variables. We compare the numerical stability of the commonly used least-squares ARX models to that of the subspace N4SID algorithm by using examples of fourth-order and eighth-order systems under pulse and chirp excitation conditions. These models correspond to structures having two and four modes simultaneously propagating respectively. We show that chirp excitation combined with the subspace identification algorithm can provide a better identification of the underlying mode dynamics than the ARX model does as the complexity of the system increases. The use of an identified state-space model for mode demixing, upon transformation to a decoupled realization form is illustrated. Applications of state-space models and the N4SID algorithm to THz transient spectroscopy as well as to optical systems are highlighted.

Optimization of apodization functions in terahertz transient spectrometry

A quadratic programming optimization procedure for designing asymmetric apodization windows tailored to the shape of time-domain sample waveforms recorded using a terahertz transient spectrometer is proposed. By artificially degrading the waveforms, the performance of the designed window in both the time and the frequency domains is compared with that of conventional rectangular, triangular (Mertz), and Hamming windows. Examples of window optimization assuming Gaussian functions as the building elements of the apodization window are provided. The formulation is sufficiently general to accommodate other basis functions. (C) 2007 Optical Society of America

MIMO Wiener model identification for large scale fading of wireless mobile communications links

The large scale fading of wireless mobile communications links is modelled assuming the mobile receiver motion is described by a dynamic linear system in state-space. The geometric relations involved in the attenuation and multi-path propagation of the electric field are described by a static non-linear mapping. A Wiener system subspace identification algorithm in conjunction with polynomial regression is used to identify a model from time-domain estimates of the field intensity assuming a multitude of emitters and an antenna array at the receiver end.

Comparison of subspace and ARX models of a waveguide's terahertz transient response after optimal wavelet filtering

A quasi-optical deembedding technique for characterizing waveguides is demonstrated using wide-band time-resolved terahertz spectroscopy. A transfer function representation is adopted for the description of the signal in the input and output port of the waveguides. The time-domain responses were discretized and the waveguide transfer function was obtained through a parametric approach in the z-domain after describing the system with an AutoRegressive with eXogenous input (ARX), as well as with a state-space model. Prior to the identification procedure, filtering was performed in the wavelet domain to minimize both signal distortion, as well as the noise propagating in the ARX and subspace models. The optimal filtering procedure used in the wavelet domain for the recorded time-domain signatures is described in detail. The effect of filtering prior to the identification procedures is elucidated with the aid of pole-zero diagrams. Models derived from measurements of terahertz transients in a precision WR-8 waveguide adjustable short are presented.

Measurement of propagation constant in waveguides with wideband coherent terahertz spectroscopy

A quasi-optical technique for characterizing micromachined waveguides is demonstrated with wideband time-resolved terahertz spectroscopy. A transfer-function representation is adopted for the description of the relation between the signals in the input and output port of the waveguides. The time-domain responses were discretized, and the waveguide transfer function was obtained through a parametric approach in the z domain after describing the system with an autoregressive with exogenous input model. The a priori assumption of the number of modes propagating in the structure was inferred from comparisons of the theoretical with the measured characteristic impedance as well as with parsimony arguments. Measurements for a precision WR-8 waveguide-adjustable short as well as for G-band reduced-height micromachined waveguides are presented. (C) 2003 Optical Society of America.

Heterosite effects in novel heteronuclear clusters [OS2Ru(CO)(11)(PPh3)] and Os2Ru(CO)(10)(2-acetylpyridine-N-isopropylimine)]

A new synthetic route towards the mixed-metal cluster [OS2Ru(CO)(12)] is described together with the syntheses of its PPh3 and iPr-AcPy (iPr-AcPy = 2-acetylpyridine-N-isopropylimine) derivatives. The molecular structures of the novel clusters [Os2Ru(CO)(11)(PPh3)] and [Os2Ru(CO)(10)(iPr-AcPy)] were determined on the basis of crystalline solid solutions of the Os2Ru and corresponding Os-3 species. The structures reveal that coordination of the Lewis bases occurs exclusively at the ruthenium site of [Os2Ru(CO)(12)], which is in agreement with density functional theory (DFT) calculations on several structural isomers of these compounds. According to the time-dependent DFT results, the lowest optically accessible excited state of [Os2Ru(CO)(10)(iPr-AcPy)] has a prevailing sigma(Ru-Os-2)pi*(iPr-AcPy) character, with a partial sigma sigma*(Ru-Os-2) contribution. In weakly coordinating 2-chlorobutane, the excited state has a lifetime tau = 10.4 +/- 1.2 ps and produces biradicals considerably faster than observed for [Os-3(CO)10(iPr-AcPy) (tau = 25.3 +/- 0.7ps)]. In coordinating acetonitrile, the excited state of [Os2Ru(CO)(10)(iPr-AcPy)] decays mono-exponentially with a lifetime tau = 2.1 +/- 0.2 ps. In contrast to [Os-3(CO)(10)(iPr-AcPy)] that forms biradicals as the main primary photoproduct even in strongly coordinating solvents, zwitterion formation from the solvated lowest excited state is observed for the heterometallic cluster. This is concluded from time-resolved absorption studies in the microsecond time domain. Due to the lower tendency of the coordinatively unsaturated Ru+(CO)(2)(iPr-AcPy-/0) moiety to bind a Lewis base, the heteronuclear biradical and zwitterionic photoproducts live significantly shorter than their triosmium counterparts. The influence of the weaker Os-2-Ru(iPr-AcPy) bond on the redox reactivity is clearly reflected in very reactive radical anions formed upon electrochemical reduction of [Os2Ru(CO)(10)(iPr-AcPy)]. The dimer [-OS(CO)(4)-Os(CO)(4)-Ru(CO)(2)(iPr-AcPy)](2)(2-) is the only IR-detectable intermediate reduction product. The dinuclear complex [Os-2(CO)(8)](2-) and insoluble [Ru(CO)(2)(iPr-AcPy)](n), are the ultimate reduction products, proving fragmentation of the OS2Ru core.

Distributed water ingress and water potential measurements using fibre optics

We report on a distributed moisture detection scheme which uses a cable design based on waterswellable hydrogel polymers. The cable modulates the loss characteristic of light guided within a multi-mode optical fibre in response to relative water potentials in the surrounding environment. Interrogation of the cable using conventional optical time-domain reflectometry (OTDR) instruments allows water ingress points to be identified and located with a spatial resolution of 50 cm. The system has been tested in a simulated tendon duct grouting experiment as a means of mapping the extent of fill along the duct during the grouting process. Voided regions were detected and identified to within 50 cm. A series of salt solutions has been used to determine the sensor behaviour over a range of water potentials. These experiments predict that measurements of soil moisture content can be made over the range 0 to – 1500 kPa. Preliminary data on soil measurements have shown that the sensor can detect water pressure changes with a resolution of 45 kPa. Applications for the sensor include quality assurance of grouting procedures, verification of waterproofing barriers and soil moisture content determination (for load-bearing calculations).

Multiband OFDM modulation and demodulation for ultra wideband communications

This chapter considers the Multiband Orthogonal Frequency Division Multiplexing (MB- OFDM) modulation and demodulation with the intention to optimize the Ultra-Wideband (UWB) system performance. OFDM is a type of multicarrier modulation and becomes the most important aspect for the MB-OFDM system performance. It is also a low cost digital signal component efficiently using Fast Fourier Transform (FFT) algorithm to implement the multicarrier orthogonality. Within the MB-OFDM approach, the OFDM modulation is employed in each 528 MHz wide band to transmit the data across the different bands while also using the frequency hopping technique across different bands. Each parallel bit stream can be mapped onto one of the OFDM subcarriers. Quadrature Phase Shift Keying (QPSK) and Dual Carrier Modulation (DCM) are currently used as the modulation schemes for MB-OFDM in the ECMA-368 defined UWB radio platform. A dual QPSK soft-demapper is suitable for ECMA-368 that exploits the inherent Time-Domain Spreading (TDS) and guard symbol subcarrier diversity to improve the receiver performance, yet merges decoding operations together to minimize hardware and power requirements. There are several methods to demap the DCM, which are soft bit demapping, Maximum Likelihood (ML) soft bit demapping, and Log Likelihood Ratio (LLR) demapping. The Channel State Information (CSI) aided scheme coupled with the band hopping information is used as a further technique to improve the DCM demapping performance. ECMA-368 offers up to 480 Mb/s instantaneous bit rate to the Medium Access Control (MAC) layer, but depending on radio channel conditions dropped packets unfortunately result in a lower throughput. An alternative high data rate modulation scheme termed Dual Circular 32-QAM that fits within the configuration of the current standard increasing system throughput thus maintaining the high rate throughput even with a moderate level of dropped packets.

Multimodel estimate of the global terrestrial water balance: setup and first results

Six land surface models and five global hydrological models participate in a model intercomparison project (WaterMIP), which for the first time compares simulation results of these different classes of models in a consistent way. In this paper the simulation setup is described and aspects of the multi-model global terrestrial water balance are presented. All models were run at 0.5 degree spatial resolution for the global land areas for a 15-year period (1985-1999) using a newly-developed global meteorological dataset. Simulated global terrestrial evapotranspiration, excluding Greenland and Antarctica, ranges from 415 to 586 mm year-1 (60,000 to 85,000 km3 year-1) and simulated runoff ranges from 290 to 457 mm year-1 (42,000 to 66,000 km3 year-1). Both the mean and median runoff fractions for the land surface models are lower than those of the global hydrological models, although the range is wider. Significant simulation differences between land surface and global hydrological models are found to be caused by the snow scheme employed. The physically-based energy balance approach used by land surface models generally results in lower snow water equivalent values than the conceptual degree-day approach used by global hydrological models. Some differences in simulated runoff and evapotranspiration are explained by model parameterizations, although the processes included and parameterizations used are not distinct to either land surface models or global hydrological models. The results show that differences between model are major sources of uncertainty. Climate change impact studies thus need to use not only multiple climate models, but also some other measure of uncertainty, (e.g. multiple impact models).

Well-posed two-point initial-boundary value problems with arbitrary boundary conditions

We study initial-boundary value problems for linear evolution equations of arbitrary spatial order, subject to arbitrary linear boundary conditions and posed on a rectangular 1-space, 1-time domain. We give a new characterisation of the boundary conditions that specify well-posed problems using Fokas' transform method. We also give a sufficient condition guaranteeing that the solution can be represented using a series. The relevant condition, the analyticity at infinity of certain meromorphic functions within particular sectors, is significantly more concrete and easier to test than the previous criterion, based on the existence of admissible functions.

The impact of resolution on the adjustment and decadal variability of the Atlantic Meridional Overturning Circulation in a coupled climate model

Variations in the Atlantic Meridional Overturning Circulation (MOC) exert an important influence on climate, particularly on decadal time scales. Simulation of the MOC in coupled climate models is compromised, to a degree that is unknown, by their lack of fidelity in resolving some of the key processes involved. There is an overarching need to increase the resolution and fidelity of climate models, but also to assess how increases in resolution influence the simulation of key phenomena such as the MOC. In this study we investigate the impact of significantly increasing the (ocean and atmosphere) resolution of a coupled climate model on the simulation of MOC variability by comparing high and low resolution versions of the same model. In both versions, decadal variability of the MOC is closely linked to density anomalies that propagate from the Labrador Sea southward along the deep western boundary. We demonstrate that the MOC adjustment proceeds more rapidly in the higher resolution model due the increased speed of western boundary waves. However, the response of the Atlantic Sea Surface Temperatures (SSTs) to MOC variations is relatively robust - in pattern if not in magnitude - across the two resolutions. The MOC also excites a coupled ocean-atmosphere response in the tropical Atlantic in both model versions. In the higher resolution model, but not the lower resolution model, there is evidence of a significant response in the extratropical atmosphere over the North Atlantic 6 years after a maximum in the MOC. In both models there is evidence of a weak negative feedback on deep density anomalies in the Labrador Sea, and hence on the MOC (with a time scale of approximately ten years). Our results highlight the need for further work to understand the decadal variability of the MOC and its simulation in climate models.

A model of UK domestic appliance-use using Artificial Neural Networks

The development of an Artificial Neural Network model of UK domestic appliance energy consumption is presented. The model uses diary-style appliance use data and a survey questionnaire collected from 51 households during the summer of 2010. It also incorporates measured energy data and is sensitive to socioeconomic, physical dwelling and temperature variables. A prototype model is constructed in MATLAB using a two layer feed forward network with backpropagation training and has a12:10:24architecture.Model outputs include appliance load profiles which can be applied to the fields of energy planning (micro renewables and smart grids), building simulation tools and energy policy.

Decision making for HVAC&R system selection for a typical office building in the UK

Demands for thermal comfort, better indoor air quality together with lower environmental impacts have had ascending trends in the last decade. In many circumstances, these demands could not be fully covered through the soft approach of bioclimatic design like optimisation of the building orientation and internal layout. This is mostly because of the dense urban environment and building internal energy loads. In such cases, heating, ventilation, air-conditioning and refrigeration (HVAC&R) systems make a key role to fulfill the requirements of indoor environment. Therefore, it is required to select the most proper HVAC&R system. In this study, a robust decision making approach for HVAC&R system selection is proposed. Technical performance, economic aspect and environmental impacts of 36 permutations of primary and secondary systems are taken into account to choose the most proper HVAC&R system for a case study office building. The building is a representative for the dominant form of office buildings in the UK. Dynamic performance evaluation of HVAC&R alternatives using TRNSYS package together with life cycle energy cost analysis provides a reliable basis for decision making. Six scenarios broadly cover the decision makers' attitudes on HVAC&R system selection which are analysed through Analytical Hierarchy Process (AHP). One of the significant outcomes reveals that, despite both the higher energy demand and more investment requirements associated with compound heating, cooling and power system (CCHP); this system is one of the top ranked alternatives due to the lower energy cost and C02 emissions. The sensitivity analysis reveals that in all six scenarios, the first five top ranked alternatives are not changed. Finally, the proposed approach and the results could be used by researchers and designers especially in the early stages of a design process in which all involved bodies face the lack of time, information and tools for evaluation of a variety of systems.

A combined engineering and statistical model of UK domestic appliance electrical load profiles

The development of a combined engineering and statistical Artificial Neural Network model of UK domestic appliance load profiles is presented. The model uses diary-style appliance use data and a survey questionnaire collected from 51 suburban households and 46 rural households during the summer of 2010 and2011 respectively. It also incorporates measured energy data and is sensitive to socioeconomic, physical dwelling and temperature variables. A prototype model is constructed in MATLAB using a two layer feed forward network with back propagation training which has a 12:10:24 architecture. Model outputs include appliance load profiles which can be applied to the fields of energy planning (microrenewables and smart grids), building simulation tools and energy policy.

Reduction of vestigial sideband effects due to multipath transmission in terrestrial broadcast television

In terrestrial television transmission multiple paths of various lengths can occur between the transmitter and the receiver. Such paths occur because of reflections from objects outside the direct transmission path. The multipath signals arriving at the receiver are all detected along with the intended signal causing time displaced replicas called 'ghosts' to appear on the television picture. With an increasing number of people living within built up areas, ghosting is becoming commonplace and therefore deghosting is becoming increasingly important. This thesis uses a deterministic time domain approach to deghosting, resulting in a simple solution to the problem of removing ghosts. A new video detector is presented which reduces the synchronous detector local oscillator phase error, caused by any practical size of ghost, to a lower level than has ever previously been achieved. From the new detector, dispersion of the video signal is minimised and a known closed-form time domain description of the individual ghost components within the detected video is subsequently obtained. Developed from mathematical descriptions of the detected video, a new specific deghoster filter structure is presented which is capable of removing both inphase (I) and also the phase quadrature (Q) induced ghost signals derived from the VSB operation. The new deghoster filter requires much less hardware than any previous deghoster which is capable of removing both I and Q ghost components. A new channel identification algorithm was also required and written which is based upon simple correlation techniques to find the delay and complex amplitude characteristics of individual ghosts. The result of the channel identification is then passed to the new I and Q deghoster filter for ghost cancellation. Generated from the research work performed for this thesis, five papers have been published. D