HIGH PERFORMANCE, LOW COST SUBSPACE DECOMPOSITION AND POLYNOMIAL ROOTING FOR REAL TIME DIRECTION OF ARRIVAL ESTIMATION: ANALYSIS AND IMPLEMENTATION

This thesis develops high performance real-time signal processing modules for direction of arrival (DOA) estimation for localization systems. It proposes highly parallel algorithms for performing subspace decomposition and polynomial rooting, which are otherwise traditionally implemented using sequential algorithms. The proposed algorithms address the emerging need for real-time localization for a wide range of applications. As the antenna array size increases, the complexity of signal processing algorithms increases, making it increasingly difficult to satisfy the real-time constraints. This thesis addresses real-time implementation by proposing parallel algorithms, that maintain considerable improvement over traditional algorithms, especially for systems with larger number of antenna array elements. Singular value decomposition (SVD) and polynomial rooting are two computationally complex steps and act as the bottleneck to achieving real-time performance. The proposed algorithms are suitable for implementation on field programmable gated arrays (FPGAs), single instruction multiple data (SIMD) hardware or application specific integrated chips (ASICs), which offer large number of processing elements that can be exploited for parallel processing. The designs proposed in this thesis are modular, easily expandable and easy to implement. Firstly, this thesis proposes a fast converging SVD algorithm. The proposed method reduces the number of iterations it takes to converge to correct singular values, thus achieving closer to real-time performance. A general algorithm and a modular system design are provided making it easy for designers to replicate and extend the design to larger matrix sizes. Moreover, the method is highly parallel, which can be exploited in various hardware platforms mentioned earlier. A fixed point implementation of proposed SVD algorithm is presented. The FPGA design is pipelined to the maximum extent to increase the maximum achievable frequency of operation. The system was developed with the objective of achieving high throughput. Various modern cores available in FPGAs were used to maximize the performance and details of these modules are presented in detail. Finally, a parallel polynomial rooting technique based on Newton’s method applicable exclusively to root-MUSIC polynomials is proposed. Unique characteristics of root-MUSIC polynomial’s complex dynamics were exploited to derive this polynomial rooting method. The technique exhibits parallelism and converges to the desired root within fixed number of iterations, making this suitable for polynomial rooting of large degree polynomials. We believe this is the first time that complex dynamics of root-MUSIC polynomial were analyzed to propose an algorithm. In all, the thesis addresses two major bottlenecks in a direction of arrival estimation system, by providing simple, high throughput, parallel algorithms.

Fading pdf of free-space optical communication system with pointing error

A free-space optical (FSO) laser communication system with perfect fast-tracking experiences random power fading due to atmospheric turbulence. For a FSO communication system without fast-tracking or with imperfect fast-tracking, the fading probability density function (pdf) is also affected by the pointing error. In this thesis, the overall fading pdfs of FSO communication system with pointing errors are calculated using an analytical method based on the fast-tracked on-axis and off-axis fading pdfs and the fast-tracked beam profile of a turbulence channel. The overall fading pdf is firstly studied for the FSO communication system with collimated laser beam. Large-scale numerical wave-optics simulations are performed to verify the analytically calculated fading pdf with collimated beam under various turbulence channels and pointing errors. The calculated overall fading pdfs are almost identical to the directly simulated fading pdfs. The calculated overall fading pdfs are also compared with the gamma-gamma (GG) and the log-normal (LN) fading pdf models. They fit better than both the GG and LN fading pdf models under different receiver aperture sizes in all the studied cases. Further, the analytical method is expanded to the FSO communication system with beam diverging angle case. It is shown that the gamma pdf model is still valid for the fast-tracked on-axis and off-axis fading pdfs with point-like receiver aperture when the laser beam is propagated with beam diverging angle. Large-scale numerical wave-optics simulations prove that the analytically calculated fading pdfs perfectly fit the overall fading pdfs for both focused and diverged beam cases. The influence of the fast-tracked on-axis and off-axis fading pdfs, the fast-tracked beam profile, and the pointing error on the overall fading pdf is also discussed. At last, the analytical method is compared with the previous heuristic fading pdf models proposed since 1970s. Although some of previously proposed fading pdf models provide close fit to the experiment and simulation data, these close fits only exist under particular conditions. Only analytical method shows accurate fit to the directly simulated fading pdfs under different turbulence strength, propagation distances, receiver aperture sizes and pointing errors.

Statistics of the received power for free space optical channels

Free space optical (FSO) communication links can experience extreme signal degradation due to atmospheric turbulence induced spatial and temporal irradiance fuctuations (scintillation) in the laser wavefront. In addition, turbulence can cause the laser beam centroid to wander resulting in power fading, and sometimes complete loss of the signal. Spreading of the laser beam and jitter are also artifacts of atmospheric turbulence. To accurately predict the signal fading that occurs in a laser communication system and to get a true picture of how this affects crucial performance parameters like bit error rate (BER) it is important to analyze the probability density function (PDF) of the integrated irradiance fuctuations at the receiver. In addition, it is desirable to find a theoretical distribution that accurately models these ?uctuations under all propagation conditions. The PDF of integrated irradiance fuctuations is calculated from numerical wave-optic simulations of a laser after propagating through atmospheric turbulence to investigate the evolution of the distribution as the aperture diameter is increased. The simulation data distribution is compared to theoretical gamma-gamma and lognormal PDF models under a variety of scintillation regimes from weak to very strong. Our results show that the gamma-gamma PDF provides a good fit to the simulated data distribution for all aperture sizes studied from weak through moderate scintillation. In strong scintillation, the gamma-gamma PDF is a better fit to the distribution for point-like apertures and the lognormal PDF is a better fit for apertures the size of the atmospheric spatial coherence radius ρ0 or larger. In addition, the PDF of received power from a Gaussian laser beam, which has been adaptively compensated at the transmitter before propagation to the receiver of a FSO link in the moderate scintillation regime is investigated. The complexity of the adaptive optics (AO) system is increased in order to investigate the changes in the distribution of the received power and how this affects the BER. For the 10 km link, due to the non-reciprocal nature of the propagation path the optimal beam to transmit is unknown. These results show that a low-order level of complexity in the AO provides a better estimate for the optimal beam to transmit than a higher order for non-reciprocal paths. For the 20 km link distance it was found that, although minimal, all AO complexity levels provided an equivalent improvement in BER and that no AO complexity provided the correction needed for the optimal beam to transmit. Finally, the temporal power spectral density of received power from a FSO communication link is investigated. Simulated and experimental results for the coherence time calculated from the temporal correlation function are presented. Results for both simulation and experimental data show that the coherence time increases as the receiving aperture diameter increases. For finite apertures the coherence time increases as the communication link distance is increased. We conjecture that this is due to the increasing speckle size within the pupil plane of the receiving aperture for an increasing link distance.

Finite element study on fracture patterns in human scaphoid wrist bone due to free fall

Scaphoid is one of the 8 carpal bones found adjacent to the thumb supported proximally by Radius bone. During the free fall, on outstretched hand, the impact load gets transferred to the scaphoid at its free anterior end. Unique arrangement of other carpal bones in the palm is also one of the reasons for the load to get transferred to scaphoid. About half of the total load acting upon carpal bone gets transferred to scaphoid at its distal pole. There are about 10 to 12 clinically observed fracture pattern in the scaphoid due to free fall. The aim of the study is to determine the orientation of the load, magnitude of the load and the corresponding fracture pattern. This study includes both static and dynamic finite element models validated by experiments. The scaphoid model has been prepared from CT scans of a 27 year old person. The 2D slices of the CT scans have been converted to 3D model by using MIMICS software. There are four cases of loading studied which are considered to occur clinically more frequently. In case (i) the load is applied at the posterior end at distal pole whereas in case (ii), (iii) and (iv), the load is applied at anterior end at different directions. The model is given a fixed boundary condition at the region which is supported by Radius bone during the impact. Same loading and boundary conditions have been used in both static and dynamic explicit finite element analysis. The site of fracture initiation and path of fracture propagation have been identified by using max principal stress / gradient and max principal strain / gradient criterion respectively in static and dynamic explicit finite element analysis. Static and dynamic impact experiments were performed on the polyurethane foam specimens to validate the finite element results. Experimental results such as load at fracture, site of fracture initiation and path of fracture propagation have been compared with the results of finite element analysis. Four different types of fracture patterns observed in clinical studies have been identified in this study.

AN EXPERIMENTAL STUDY OF PARTIAL LOW RISK DEPLOYMENT IN A FIXED VEHICLE ENVIRONMENT

This study will look at the passenger air bag (PAB) performance in a fix vehicle environment using Partial Low Risk Deployment (PLRD) as a strategy. This development will follow test methods against actual baseline vehicle data and Federal Motor Vehicle Safety Standards 208 (FMVSS 208). FMVSS 208 states that PAB compliance in vehicle crash testing can be met using one of three deployment methods. The primary method suppresses PAB deployment, with the use of a seat weight sensor or occupant classification sensor (OCS), for three-year old and six-year old occupants including the presence of a child seat. A second method, PLRD allows deployment on all size occupants suppressing only for the presents of a child seat. A third method is Low Risk Deployment (LRD) which allows PAB deployment in all conditions, all statures including any/all child seats. This study outlines a PLRD development solution for achieving FMVSS 208 performance. The results of this study should provide an option for system implementation including opportunities for system efficiency and other considerations. The objective is to achieve performance levels similar too or incrementally better than the baseline vehicles National Crash Assessment Program (NCAP) Star rating. In addition, to define systemic flexibility where restraint features can be added or removed while improving occupant performance consistency to the baseline. A certified vehicles’ air bag system will typically remain in production until the vehicle platform is redesigned. The strategy to enable the PLRD hypothesis will be to first match the baseline out of position occupant performance (OOP) for the three and six-year old requirements. Second, improve the 35mph belted 5th percentile female NCAP star rating over the baseline vehicle. Third establish an equivalent FMVSS 208 certification for the 25mph unbelted 50th percentile male. FMVSS 208 high-speed requirement defines the federal minimum crash performance required for meeting frontal vehicle crash-test compliance. The intent of NCAP 5-Star rating is to provide the consumer with information about crash protection, beyond what is required by federal law. In this study, two vehicles segments were used for testing to compare and contrast to their baseline vehicles performance. Case Study 1 (CS1) used a cross over vehicle platform and Case Study 2 (CS2) used a small vehicle segment platform as their baselines. In each case study, the restraints systems were from different restraint supplier manufactures and each case contained that suppliers approach to PLRD. CS1 incorporated a downsized twins shaped bag, a carryover inflator, standard vents, and a strategic positioned bag diffuser to help disperse the flow of gas to improve OOP. The twin shaped bag with two segregated sections (lobes) to enabled high-speed baseline performance correlation on the HYGE Sled. CS2 used an A-Symmetric (square shape) PAB with standard size vents, including a passive vent, to obtain OOP similar to the baseline. The A-Symmetric shape bag also helped to enabled high-speed baseline performance improvements in HYGE Sled testing in CS2. The anticipated CS1 baseline vehicle-pulse-index (VPI) target was in the range of 65-67. However, actual dynamic vehicle (barrier) testing was overshadowed with the highest crash pulse from the previous tested vehicles with a VPI of 71. The result from the 35mph NCAP Barrier test was a solid 4-Star (4.7 Star) respectfully. In CS2, the vehicle HYGE Sled development VPI range, from the baseline was 61-62 respectively. Actual NCAP test produced a chest deflection result of 26mm versus the anticipated baseline target of 12mm. The initial assessment of this condition was thought to be due to the vehicles significant VPI increase to 67. A subsequent root cause investigation confirmed a data integrity issue due to the instrumentation. In an effort to establish a true vehicle test data point a second NCAP test was performed but faced similar instrumentation issues. As a result, the chest deflect hit the target of 12.1mm; however a femur load spike, similar to the baseline, now skewed the results. With noted level of performance improvement in chest deflection, the NCAP star was assessed as directional for 5-Star capable performance. With an actual rating of 3-Star due to instrumentation, using data extrapolation raised the ratings to 5-Star. In both cases, no structural changes were made to the surrogate vehicle and the results in each case matched their perspective baseline vehicle platforms. These results proved the PLRD is viable for further development and production implementation.

Free Troposphere Ozone and Carbon Monoxide over the North Atlantic for 2001-2011

Carbon Monoxide (CO) and Ozone (O3) are considered to be one of the most important atmospheric pollutants in the troposphere with both having significant effects on human health. Both are included in the U.S. E.P.A list of criteria pollutants. CO is primarily emitted in the source region whereas O3 can be formed near the source, during transport of the pollution plumes containing O3 precursors or in a receptor region as the plumes subside. The long chemical lifetimes of both CO and O3 enable them to be transported over long distances. This transport is important on continental scales as well, commonly referred to as inter-continental transport and affects the concentrations of both CO and O3 in downwind receptor regions, thereby having significant implications for their air quality standards. Over the period 2001-2011, there have been decreases in the anthropogenic emissions of CO and NOx in North America and Europe whereas the emissions over Asia have increased. How these emission trends have affected concentrations at remote sites located downwind of these continents is an important question. The PICO-NARE observatory located on the Pico Mountain in Azores, Portugal is frequently impacted by North American pollution outflow (both anthropogenic and biomass burning) and is a unique site to investigate long range transport from North America. This study uses in-situ observations of CO and O3 for the period 2001-2011 at PICO-NARE coupled with output from the full chemistry (with normal and fixed anthropogenic emissions) and tagged CO simulations in GEOS-Chem, a global 3-D chemical transport model of atmospheric composition driven by meteorological input from the Goddard Earth Observing System (GEOS) of the NASA Global Modeling and Assimilation Office, to determine and interpret the trends in CO and O3 concentrations over the past decade. These trends would be useful in ascertaining the impacts emission reductions in the United States have had over Pico and in general over the North Atlantic. A regression model with sinusoidal functions and a linear trend term was fit to the in-situ observations and the GEOS-Chem output for CO and O3 at Pico respectively. The regression model yielded decreasing trends for CO and O3 with the observations (-0.314 ppbv/year & -0.208 ppbv/year respectively) and the full chemistry simulation with normal emissions (-0.343 ppbv/year & -0.526 ppbv/year respectively). Based on analysis of the results from the full chemistry simulation with fixed anthropogenic emissions and the tagged CO simulation it was concluded that the decreasing trends in CO were a consequence of the anthropogenic emission changes in regions such as USA and Asia. The emission reductions in USA are countered by Asian increases but the former have a greater impact resulting in decreasing trends for CO at PICO-NARE. For O3 however, it is the increase in water vapor content (which increases O3 destruction) along the pathways of transport from North America to PICO-NARE as well as around the site that has resulted in decreasing trends over this period. This decrease is offset by increase in O3 concentrations due to anthropogenic influence which could be due to increasing Asian emissions of O3 precursors as these emissions have decreased over the US. However, the anthropogenic influence does not change the final direction of the trend. It can thus be concluded that CO and O3 concentrations at PICO-NARE have decreased over 2001-2011.