Transceiver designs and analysis for LTI, LTV and broadcast channels - new matrix decompositions and majorization theory

Signal processing techniques play important roles in the design of digital communication systems. These include information manipulation, transmitter signal processing, channel estimation, channel equalization and receiver signal processing. By interacting with communication theory and system implementing technologies, signal processing specialists develop efficient schemes for various communication problems by wisely exploiting various mathematical tools such as analysis, probability theory, matrix theory, optimization theory, and many others. In recent years, researchers realized that multiple-input multiple-output (MIMO) channel models are applicable to a wide range of different physical communications channels. Using the elegant matrix-vector notations, many MIMO transceiver (including the precoder and equalizer) design problems can be solved by matrix and optimization theory. Furthermore, the researchers showed that the majorization theory and matrix decompositions, such as singular value decomposition (SVD), geometric mean decomposition (GMD) and generalized triangular decomposition (GTD), provide unified frameworks for solving many of the point-to-point MIMO transceiver design problems.

In this thesis, we consider the transceiver design problems for linear time invariant (LTI) flat MIMO channels, linear time-varying narrowband MIMO channels, flat MIMO broadcast channels, and doubly selective scalar channels. Additionally, the channel estimation problem is also considered. The main contributions of this dissertation are the development of new matrix decompositions, and the uses of the matrix decompositions and majorization theory toward the practical transmit-receive scheme designs for transceiver optimization problems. Elegant solutions are obtained, novel transceiver structures are developed, ingenious algorithms are proposed, and performance analyses are derived.

The first part of the thesis focuses on transceiver design with LTI flat MIMO channels. We propose a novel matrix decomposition which decomposes a complex matrix as a product of several sets of semi-unitary matrices and upper triangular matrices in an iterative manner. The complexity of the new decomposition, generalized geometric mean decomposition (GGMD), is always less than or equal to that of geometric mean decomposition (GMD). The optimal GGMD parameters which yield the minimal complexity are derived. Based on the channel state information (CSI) at both the transmitter (CSIT) and receiver (CSIR), GGMD is used to design a butterfly structured decision feedback equalizer (DFE) MIMO transceiver which achieves the minimum average mean square error (MSE) under the total transmit power constraint. A novel iterative receiving detection algorithm for the specific receiver is also proposed. For the application to cyclic prefix (CP) systems in which the SVD of the equivalent channel matrix can be easily computed, the proposed GGMD transceiver has K/log_2(K) times complexity advantage over the GMD transceiver, where K is the number of data symbols per data block and is a power of 2. The performance analysis shows that the GGMD DFE transceiver can convert a MIMO channel into a set of parallel subchannels with the same bias and signal to interference plus noise ratios (SINRs). Hence, the average bit rate error (BER) is automatically minimized without the need for bit allocation. Moreover, the proposed transceiver can achieve the channel capacity simply by applying independent scalar Gaussian codes of the same rate at subchannels.

In the second part of the thesis, we focus on MIMO transceiver design for slowly time-varying MIMO channels with zero-forcing or MMSE criterion. Even though the GGMD/GMD DFE transceivers work for slowly time-varying MIMO channels by exploiting the instantaneous CSI at both ends, their performance is by no means optimal since the temporal diversity of the time-varying channels is not exploited. Based on the GTD, we develop space-time GTD (ST-GTD) for the decomposition of linear time-varying flat MIMO channels. Under the assumption that CSIT, CSIR and channel prediction are available, by using the proposed ST-GTD, we develop space-time geometric mean decomposition (ST-GMD) DFE transceivers under the zero-forcing or MMSE criterion. Under perfect channel prediction, the new system minimizes both the average MSE at the detector in each space-time (ST) block (which consists of several coherence blocks), and the average per ST-block BER in the moderate high SNR region. Moreover, the ST-GMD DFE transceiver designed under an MMSE criterion maximizes Gaussian mutual information over the equivalent channel seen by each ST-block. In general, the newly proposed transceivers perform better than the GGMD-based systems since the super-imposed temporal precoder is able to exploit the temporal diversity of time-varying channels. For practical applications, a novel ST-GTD based system which does not require channel prediction but shares the same asymptotic BER performance with the ST-GMD DFE transceiver is also proposed.

The third part of the thesis considers two quality of service (QoS) transceiver design problems for flat MIMO broadcast channels. The first one is the power minimization problem (min-power) with a total bitrate constraint and per-stream BER constraints. The second problem is the rate maximization problem (max-rate) with a total transmit power constraint and per-stream BER constraints. Exploiting a particular class of joint triangularization (JT), we are able to jointly optimize the bit allocation and the broadcast DFE transceiver for the min-power and max-rate problems. The resulting optimal designs are called the minimum power JT broadcast DFE transceiver (MPJT) and maximum rate JT broadcast DFE transceiver (MRJT), respectively. In addition to the optimal designs, two suboptimal designs based on QR decomposition are proposed. They are realizable for arbitrary number of users.

Finally, we investigate the design of a discrete Fourier transform (DFT) modulated filterbank transceiver (DFT-FBT) with LTV scalar channels. For both cases with known LTV channels and unknown wide sense stationary uncorrelated scattering (WSSUS) statistical channels, we show how to optimize the transmitting and receiving prototypes of a DFT-FBT such that the SINR at the receiver is maximized. Also, a novel pilot-aided subspace channel estimation algorithm is proposed for the orthogonal frequency division multiplexing (OFDM) systems with quasi-stationary multi-path Rayleigh fading channels. Using the concept of a difference co-array, the new technique can construct M^2 co-pilots from M physical pilot tones with alternating pilot placement. Subspace methods, such as MUSIC and ESPRIT, can be used to estimate the multipath delays and the number of identifiable paths is up to O(M^2), theoretically. With the delay information, a MMSE estimator for frequency response is derived. It is shown through simulations that the proposed method outperforms the conventional subspace channel estimator when the number of multipaths is greater than or equal to the number of physical pilots minus one.

Nonstationary response of structures and its application to earthquake engineering

This thesis presents a simplified state-variable method to solve for the nonstationary response of linear MDOF systems subjected to a modulated stationary excitation in both time and frequency domains. The resulting covariance matrix and evolutionary spectral density matrix of the response may be expressed as a product of a constant system matrix and a time-dependent matrix, the latter can be explicitly evaluated for most envelopes currently prevailing in engineering. The stationary correlation matrix of the response may be found by taking the limit of the covariance response when a unit step envelope is used. The reliability analysis can then be performed based on the first two moments of the response obtained.

The method presented facilitates obtaining explicit solutions for general linear MDOF systems and is flexible enough to be applied to different stochastic models of excitation such as the stationary models, modulated stationary models, filtered stationary models, and filtered modulated stationary models and their stochastic equivalents including the random pulse train model, filtered shot noise, and some ARMA models in earthquake engineering. This approach may also be readily incorporated into finite element codes for random vibration analysis of linear structures.

A set of explicit solutions for the response of simple linear structures subjected to modulated white noise earthquake models with four different envelopes are presented as illustration. In addition, the method has been applied to three selected topics of interest in earthquake engineering, namely, nonstationary analysis of primary-secondary systems with classical or nonclassical dampings, soil layer response and related structural reliability analysis, and the effect of the vertical components on seismic performance of structures. For all the three cases, explicit solutions are obtained, dynamic characteristics of structures are investigated, and some suggestions are given for aseismic design of structures.

Stationary random response of multidegree-of-freedom systems

An approximate approach is presented for determining the stationary random response of a general multidegree-of-freedom nonlinear system under stationary Gaussian excitation. This approach relies on defining an equivalent linear system for the nonlinear system. Two particular systems which possess exact solutions have been solved by this approach, and it is concluded that this approach can generate reasonable solutions even for systems with fairly large nonlinearities. The approximate approach has also been applied to two examples for which no exact or approximate solutions were previously available.

Also presented is a matrix algebra approach for determining the stationary random response of a general multidegree-of-freedom linear system. Its derivation involves only matrix algebra and some properties of the instantaneous correlation matricies of a stationary process. It is therefore very direct and straightforward. The application of this matrix algebra approach is in general simpler than that of commonly used approaches.

Coupled foundation-superstructure analysis and influence of building stiffness on foundation response

This paper proposes a simple method to include superstructure stiffness in foundation analyses. The method involves extracting a small "condensed structural matrix" from finite element models of the superstructure, which can then be incorporated into pile group or piled raft analyses using common approaches such as elastic continuum or load transfer methods. The matrix condensation method directly couples structural and geotechnical analyses, and eliminates the need for iterative analyses between structural and geotechnical engineers. Effectiveness of the approach is illustrated through analyses of several buildings designed with a typical floor plan but with varying heights. The parametric study illustrates that superstructure stiffness can have a significant influence on foundation settlement estimates, and the stiffening effects are dominated by the lower stories of the superstructure. The proposed method aims to bridge the gap between structural and geotechnical analyses. Also, being a computationally simple and accurate approach, it is applicable to parametric or optimization studies that would otherwise involve large amounts of analyses. © 2010 ASCE.

The use of spatial impulse responses to characterise flexible forming processes with mobile tools

A novel test method for the characterisation of flexible forming processes is proposed and applied to four flexible forming processes: Incremental Sheet Forming (ISF), conventional spinning, the English wheel and power hammer. The proposed method is developed in analogy with time-domain control engineering, where a system is characterised by its impulse response. The spatial impulse response is used to characterise the change in workpiece deformation created by a process, but has also been applied with a strain spectrogram, as a novel way to characterise a process and the physical effect it has on the workpiece. Physical and numerical trials to study the effects of process and material parameters on spatial impulse response lead to three main conclusions. Incremental sheet forming is particularly sensitive to process parameters. The English wheel and power hammer are strongly similar and largely insensitive to both process and material parameters. Spinning develops in two stages and is sensitive to most process parameters, but insensitive to prior deformation. Finally, the proposed method could be applied to modelling, classification of existing and novel processes, product-process matching and closed-loop control of flexible forming processes. © 2012 Elsevier B.V.

Seismic response of stone masonry spires: Analytical modeling

Stone masonry spires are vulnerable to seismic loading. Computational methods are often used to predict the dynamic linear elastic response of masonry towers and spires, but this approach is only applicable until the first masonry joint begins to open, limiting the ability to predict collapse. In this paper, analytical modeling is used to investigate the uplift, rocking and collapse of stone spires. General equations for static equilibrium of the spire under lateral acceleration are first presented, and provide a reasonable lower bound for predicting collapse. The dynamic response is then considered through elastic modal analysis and rigid body rocking. Together, these methods are used to provide uplift curves and single impulse overturning collapse curves for a complete range of possible spire geometries. Results are used to evaluate the historic collapse of two specific stone spires. © 2012 Elsevier Ltd.

Seismic response of stone masonry spires: Computational and experimental modeling

Analytical methods provide a global context from which to understand the dynamics of stone spires, but computational and experimental methods are useful to predict more specific behavior of multiple block structures. In this paper, the spire of St. Mary Magdalene church in Waltham-on-the-Wolds, UK, which was damaged in the 2008 Lincolnshire Earthquake, is used as a case study. Both a physical model and a discrete element computational model of the spire were created and used to investigate collapse under constant horizontal acceleration, impulse base motion, and earthquake ground motion. Results indicate that the global behavior compares well with analytical modeling, but local block displacements evident in DEM and experimental results also reduce the stability of the structure. In this context, the observed damage to St. Mary Magdalene church is evaluated and discussed. © 2012 Elsevier Ltd.

The influence of matrix integrity on stress-fiber remodeling in 3D.

Matrix anisotropy is important for long term in vivo functionality. However, it is not fully understood how to guide matrix anisotropy in vitro. Experiments suggest actin-mediated cell traction contributes. Although F-actin in 2D displays a stretch-avoidance response, 3D data are lacking. We questioned how cyclic stretch influences F-actin and collagen orientation in 3D. Small-scale cell-populated fibrous tissues were statically constrained and/or cyclically stretched with or without biochemical agents. A rectangular array of silicone posts attached to flexible membranes constrained a mixture of cells, collagen I and matrigel. F-actin orientation was quantified using fiber-tracking software, fitted using a bi-model distribution function. F-actin was biaxially distributed with static constraint. Surprisingly, uniaxial cyclic stretch, only induced a strong stretch-avoidance response (alignment perpendicular to stretching) at tissue surfaces and not in the core. Surface alignment was absent when a ROCK-inhibitor was added, but also when tissues were only statically constrained. Stretch-avoidance was also observed in the tissue core upon MMP1-induced matrix perturbation. Further, a strong stretch-avoidance response was obtained for F-actin and collagen, for immediate cyclic stretching, i.e. stretching before polymerization of the collagen. Results suggest that F-actin stress-fibers avoid cyclic stretch in 3D, unless collagen contact guidance dictates otherwise.

The influence of matrix integrity on stress-fiber remodeling in 3D

Matrix anisotropy is important for long term in vivo functionality. However, it is not fully understood how to guide matrix anisotropy in vitro. Experiments suggest actin-mediated cell traction contributes. Although F-actin in 2D displays a stretch-avoidance response, 3D data are lacking. We questioned how cyclic stretch influences F-actin and collagen orientation in 3D. Small-scale cell-populated fibrous tissues were statically constrained and/or cyclically stretched with or without biochemical agents. A rectangular array of silicone posts attached to flexible membranes constrained a mixture of cells, collagen I and matrigel. F-actin orientation was quantified using fiber-tracking software, fitted using a bi-model distribution function. F-actin was biaxially distributed with static constraint. Surprisingly, uniaxial cyclic stretch, only induced a strong stretch-avoidance response (alignment perpendicular to stretching) at tissue surfaces and not in the core. Surface alignment was absent when a ROCK-inhibitor was added, but also when tissues were only statically constrained. Stretch-avoidance was also observed in the tissue core upon MMP1-induced matrix perturbation. Further, a strong stretch-avoidance response was obtained for F-actin and collagen, for immediate cyclic stretching, i.e. stretching before polymerization of the collagen. Results suggest that F-actin stress-fibers avoid cyclic stretch in 3D, unless collagen contact guidance dictates otherwise. © 2012 Elsevier Ltd.

In vitro osteoblast response to ferritic stainless steel fiber networks for magneto-active layers on implants

The use of a porous coating on prosthetic components to encourage bone ingrowth is an important way of improving uncemented implant fixation. Enhanced fixation may be achieved by the use of porous magneto-active layers on the surface of prosthetic implants, which would deform elastically on application of a magnetic field, generating internal stresses within the in-growing bone. This approach requires a ferromagnetic material able to support osteoblast attachment, proliferation, differentiation, and mineralization. In this study, the human osteoblast responses to ferromagnetic 444 stainless steel networks were considered alongside those to nonmagnetic 316L (medical grade) stainless steel networks. While both networks had similar porosities, 444 networks were made from coarser fibers, resulting in larger inter-fiber spaces. The networks were analyzed for cell morphology, distribution, proliferation, and differentiation, extracellular matrix production and the formation of mineralized nodules. Cell culture was performed in both the presence of osteogenic supplements, to encourage cell differentiation, and in their absence. It was found that fiber size affected osteoblast morphology, cytoskeleton organization and proliferation at the early stages of culture. The larger inter-fiber spaces in the 444 networks resulted in better spatial distribution of the extracellular matrix. The addition of osteogenic supplements enhanced cell differentiation and reduced cell proliferation thereby preventing the differences in proliferation observed in the absence of osteogenic supplements. The results demonstrated that 444 networks elicited favorable responses from human osteoblasts, and thus show potential for use as magnetically active porous coatings for advanced bone implant applications. © 2012 Wiley Periodicals, Inc.

The high strain rate response of Ultra High Molecular-weight Polyethylene: From fibre to laminate

The effect of strain rate upon the uniaxial response of Ultra High Molecular-weight Polyethylene (UHMWPE) fibres, yarns and laminates of lay-up [0/90]48 has been measured in both the 0/90 and ±45 configurations. The tensile strength of the matrix-dominated ±45 laminate is two orders of magnitude less than that of the fibre-dominated 0/90 laminate, and is more sensitive to strain rate. A piezoelectric force sensor device was developed to obtain the high strain rate data, and this achieved a rise time of less than 1 μs. It is found that the failure strength (and failure strain) of the yarn is almost insensitive to strain rate within the range (10 -1-103 s-1). At low strain rates (below 10 -1 s-1), creep of the yarn dominates and the failure strain increases with diminishing strain rate. The tensile strength of the dry yarn exceeds that of the laminate by about 20%. Tests on single fibres exceed the strength of the yarn by 20%. © 2013 Elsevier Ltd. All rights reserved.

The effect of shear strength on the ballistic response of laminated composite plates

The ballistic performance of clamped circular carbon fibre reinforced polymer (CFRP) and Ultra High Molecular Weight Polyethylene (UHMWPE) fibre composite plates of equal areal mass and 0/90 lay-up were measured and compared with that of monolithic 304 stainless steel plates. The effect of matrix shear strength upon the dynamic response was explored by testing: (i) CFRP plates with both a cured and uncured matrix and (ii) UHMWPE laminates with identical fibres but with two matrices of different shear strength. The response of these plates when subjected to mid-span, normal impact by a steel ball was measured via a dynamic high speed shadow moiré technique. Travelling hinges emanate from the impact location and travel towards the supports. The anisotropic nature of the composite plate results in the hinges travelling fastest along the fibre directions and this results in square-shaped moiré fringes in the 0/90 plates. Projectile penetration of the UHMWPE and the uncured CFRP plates occurs in a progressive manner, such that the number of failed plies increases with increasing velocity. The cured CFRP plate, of high matrix shear strength, fails by cone-crack formation at low velocities, and at higher velocities by a combination of cone-crack formation and communition of plies beneath the projectile. On an equal areal mass basis, the low shear strength UHMWPE plate has the highest ballistic limit followed by the high matrix shear strength UHMWPE plate, the uncured CFRP, the steel plate and finally the cured CFRP plate. We demonstrate that the high shear strength UHMWPE plate exhibits Cunniff-type ballistic limit scaling. However, the observed Cunniff velocity is significantly lower than that estimated from the laminate properties. The data presented here reveals that the Cunniff velocity is limited in its ability to characterise the ballistic performance of fibre composite plates as this velocity is independent of the shear properties of the composites: the ballistic limit of fibre composite plates increases with decreasing matrix shear strength for both CFRP and UHMWPE plates. © 2013 Elsevier Masson SAS. All rights reserved.

Transient response of structures with uncertain properties to nonlinear shock loading

A method is presented to predict the transient response of a structure at the driving point following an impact or a shock loading. The displacement and the contact force are calculated solving the discrete convolution between the impulse response and the contact force itself, expressed in terms of a nonlinear Hertzian contact stiffness. Application of random point process theory allows the calculation of the impulse response function from knowledge of the modal density and the geometric characteristics of the structure only. The theory is applied to a wide range of structures and results are experimentally verified for the case of a rigid object hitting a beam, a plate, a thin and a thick cylinder and for the impact between two cylinders. The modal density of the flexural modes for a thick slender cylinder is derived analytically. Good agreement is found between experimental, simulated and published results, showing the reliability of the method for a wide range of situations including impacts and pyroshock applications. © 2013 Elsevier Ltd. All rights reserved.

Predicting the response of structures to transient shock loading.

This work concerns the prediction of the response of an uncertain structure to a load of short duration. Assuming an ensemble of structures with small random variations about a nominal form, a mean impulse response can be found using only the modal density of the structure. The mean impulse response turns out to be the same as the response of an infinite structure: the response is calculated by taking into account the direct field only, without reflections. Considering the short duration of an impulsive loading, the approach is reasonable before the effect of the reverberant field becomes important. The convolution between the mean impulse response and the shock loading is solved in discrete time to calculate the response at the driving point and at remote points. Experimental and numerical examples are presented to validate the theory presented for simple structures such as beams, plates, and cylinders.

Dynamics of filopodium-like protrusion and endothelial cellular motility on onedimensional extracellular matrix fibrils

Endothelial filopodia play key roles in guiding the tubular sprouting during angiogenesis. However, their dynamic morphological characteristics, with the associated implications in cell motility, have been subjected to limited investigations. In this work, the interaction between endothelial cells and extracellular matrix fibrils was recapitulated in vitro, where a specific focus was paid to derive the key morphological parameters to define the dynamics of filopodium-like protrusion during cell motility. Based on one-dimensional gelatin fibrils patterned by near-field electrospinning (NFES), we study the response of endothelial cells (EA.hy926) under normal culture or ROCK inhibition. It is shown that the behaviour of temporal protrusion length versus cell motility can be divided into distinct modes. Persistent migration was found to be one of the modes which permitted cell displacement for over 300 μm at a speed of approximately 1 μm min-1. ROCK inhibition resulted in abnormally long protrusions and diminished the persistent migration, but dramatically increased the speeds of protrusion extension and retraction. Finally, we also report the breakage of protrusion during cell motility, and examine its phenotypic behaviours. © 2014 The Author(s) Published by the Royal Society. All rights reserved.