Least-Squares Spectral Analysis Theory Summary
|
|
- Gabriel Horton
- 5 years ago
- Views:
Transcription
1 Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation, Department o Geodesy and Geomatics Engineering, Technical Report No. 282, University o New Brunswick, Fredericton, New Brunswick, Canada, 208 pp. 1. LEAST-SQUARES SPECTRAL ANALYSIS Least Squares Spectral Analysis is a powerul sotware developed at the University o New Brunswick, Fredericton. LSSA was irst developed by Vaníček (1968, 1971) as an alternative to the classical Fourier methods (Pagiatakis, 1998). It has been revised by Wells, Vaníček and Pagiatakis in 1985 in order to eliminate certain bugs and to be more versatile (Wells, 1985). In this study version LSSA v is used developed by Spiros Pagiatakis. LSSA has several advantages over the commonly known Fourier transorm. Table 1 summarizes these and also the limitations o the Fourier transorm based on Pagiatakis [1998, 2008]. Table 1. Advantages o LSSA in comparison with the classical Fourier transorm Time series with unequally spaced values can be analyzed without preprocessing Time series does not have to be continuous (allowance o data gaps) Time series with an associated covariance matrix can be analyzed Time series can have datum shits and trends LSSA FOURIER TRANSFORM
2 No limitations or the length o the time series Higher accuracy achieved with longer time series Statistical testing on the signiicance o spectral peaks can be perormed Perhaps the most essential advantage o LSSA is the allowance o data gaps in the observed time series, what is common in GPS pseudorange and carrier-phase observations. The original time series can be directly used without interpolating or adding artiicial data to ill the data gaps, hence not degrading the results o the spectral analysis. LSSA also permits to add bigger weights to statistically more signiicant observations and smaller weight to the other ones. This can be useul in geodetic applications when analyzing pseudorange and carrier-phase observations, assigning more weight to the more precise carrier-phase and less weight to pseudorange observations. Let us discuss the LSSA and the other great properties which this sotware oers in the next section. 4.1 Theory o signal and noise The observed time series can be considered to be composed o signal, what we are interested to study, and noise which obscure the original signal. The noise can be random or periodic. Idealizing, we can think about the random noise as a white noise what is uncorrelated, has a constant spectral density, and it may or may not have a Gaussian distribution. Usually in practice we deal with a non-white random noise. Systematic noise is a noise what can be described by a certain unctional orm. It can be periodic (colored) or non-periodic. Non-periodic noise can be datum-shits (osets) and trends (linear, exponential, quadratic, etc.), and it causes the statistical properties o the series to be a unction o time (Pagiatakis, 1998)
3 4.2 Mathematics o the least-squares spectrum Let us consider a discrete time series (ti) in Hilbert space, (H) where ti is a vector o observation times and i 12,,, m, m is the number o data points in the time series. We assume that a ully populated covariance matrix C is available or the time series. One o the main obectives o LSSA is to detect unknown periodic signals in, especially when contains systematic variations o unknown magnitude whose unctional orms are known. Ater Vaníček (1971) we can call these known constituents systematic noise since their presence is nuisance rom the point o view o detecting the unknown periodic signal (Vaníček, 1981). The idea is to approximate with another unction g using the least squares principle, so that the dierences between and g (the vector o residuals r) are minimum in the least-squares sense. Thus the time series can be modeled by g as ollows (Pagiatakis, 1998): g Φ x, (9) where x is the vector o unknowns and is the design matrix (matrix o normal equations) which consist o several column vectors: 12,,m] called based unctions, each o which is a know unction o the same dimension as (Wells, 1985). Matrix speciies the unctional orm o both signal and (systematic) noise. Using the standard least-squares notation we can write or the residual series (Pagiatakis, 1998)
4 T 1 C 1 Φ Φ T C 1 ˆ r gˆ Φ Φ. (10) In equation (10), ĝ is the orthogonal proection o onto the subspace S (rom Hilbert space) generated by the column vector o (see Figure 1 or illustration). It ollows rom the proection theorem that ˆr gˆ. This means that has been decomposed into a signal ĝ and noise rˆ. Figure 1. Orthogonal proection o vectors ĝ and rˆ In Figure 1 the red arrow represents the orthogonal proection o ĝ into. Then the inner product (dot product) o these two vectors is:, ˆg ˆg cos (11) and the orthogonal proection is then:,gˆ ˆg cos. (12) - 4 -
5 From equation (11) one can conclude that i the angle between ĝ and is 0 degrees, then ĝ ully represents. Thus the ratio o the length o the orthogonal proection multiplied by 100% and the length o tells how much the least-squares approximation o g represents in terms o percent. The ratio is smaller than unity and can be expressed as ollows: s, gˆ, gˆ 2 T C 1 gˆ T C 1 (13) where s is the least-squares spectrum in percentage ( 0 s 1). The larger the s (closer to 1) the better is the least-squares it to the data (Pagiatakis, 2008). In spectral analysis it is usual to search or periodic signals which can be expressed in terms o sine and cosine base unctions. Speciically or LSSA we know (t) and we use: Φ cos t i,sin t i (14) where is a spectral requency or which spectral values spectral requencies can be written as: orthogonal proection o onto S will be dierent or each independently rom the rest (Pagiatakis, 1998). For each compute: s are desired. The vector o where 12,,, k (Wells, 1985). The and each is tried or which we want s we - 5 -
6 g( 1 2 ) xˆ cos t xˆ sin t, (15) where xˆ 1 ˆx is determined rom (Wells, 1985): xˆ 2 T 1 1 T 1 Φ C Φ Φ C ˆ x. (16) Then the LSS is deined by T 1 C gˆ, 12,,, k. (17) T C s 1 When calculating the LSS there is a simultaneous least-squares solution or the parameters o the process. This is a rigorous approach to the hidden periodicities: the parameters o the assumed linear system driven by noise are determined simultaneously with the amplitudes and phases o the periodic components, and with other parameters that describe systematic noise (Pagiatakis, 1998). For example the approximating unction can be a sine wave with known requency: g acos( 2 ) (18) where a is the amplitude and is the phase or a speciic requency which gives the best approximation (least-squares it) to. For dierent requencies the LSS s varies (s is a unction o - 6 -
7 requency) and so it becomes the representation o how well the sine wave its the data. The time series is best approximated when s equals to 1 or a certain predeined requency (Pagiatakis, 2008). 4.3 Least-squares spectrum evaluation Table 1 also lists that LSSA allows perorming statistical testing on the signiicance o spectral peaks. This is very important since we need to know which peak is statistically signiicant and which one can be suppressed. Vaníček [1971] derived the expected (mean) spectral value o white noise in the LSS or series consisting o statistically independent random values. He pointed out the possibility o deriving magnitudes (threshold values) above which spectral peaks are statistically signiicant (Pagiatakis, 1998). Postulating that series has been derived rom a population o random variables ollowing the multidimensional normal distribution n(0, C). This null hypothesis (H0) implies that i: v 2Fv,, Qn Qs 2 (19) then the series will contain statistically insigniicant noise. In equation (18) Q is a random variable with subscripts n and s reerring to signal and noise, respectively; v is the degree o reedom and F stands or F-distribution and is the signiicance level. The alternative hypothesis H1 is: v 2Fv,, Qn Qs 2. (20) - 7 -
8 It makes sense to use only the lower tail end o F, since large values o the ratio Q n Q s (upper tail o F) imply the signal is signiicantly larger than the noise and hence can not be detected. Then the statistically signiicant spectral peaks will satisy the ollowing inequality (Pagiatakis, 1998): 1 v s 1 F v, 2, 2. (21) Once a statistically signiicant signal is detected it should be related to the physics o the observed system. Once it is understood the detected signal becomes systematic noise and then urther hidden periodicities can be searched (Vaníček, 1981). 4.4 Input and output parameters or LSSA The input parameters to compute the spectrum consist o the time series, the limits and density o the spectral band to be produced and the known constituents base unctions (representing the systematic noise unctional orms) (Vaníček, 1985). The input ile, which contains the userdeined parameters, is called lssa.in. There are eight blocks o parameters that need to be speciied or the analysis o the series beore each run. Let us shortly describe the parameters to be deined in lssa.in (Vaníček, 1985). First, the user must enter the name o the ile to be analyzed, or example data_ser.dat, the number o data points, the units o time o the series and the units o the values o the series. The - 8 -
9 irst line in the data_ser.dat is a text line containing the name o the proect, ollowed by three columns. The irst column stands or the time o the series (years, days, seconds, etc.), the second column is the time series itsel and the third column represents the standard deviation o the time series i known. I the standard deviations are unknown the user can replace them with values he/she may think are reasonable. The units o the standard deviations are the same as the units o the time series. Next, the user can identiy the known constituents in the time series according to his/her best knowledge. There are our optional types o known constituents build in the algorithm. The sotware makes it possible or the user to speciy his/her own unctions as well. The known constituents build in LSSA are the ollows: Random constant (Datum Shits). The user can speciy the number o datum biases in the time series and the epochs when the shits are suspected. The minimum number o shits is one which represents the beginning o the time series. Linear trend. The user can decide whether to calculate linear trend or not. Periodic signals. Certain periodic constituents can be orced (itted) to the series. The result is a sine wave given by its amplitude and phase, along with a statistical test on its signiicance. The user can enter the number o sine waves to be orced, the periods and names o each sine wave. User speciied I the user believes that some nonlinear trend (exponential, polynomial) exists he/she can speciy it here. In general, a polynomial it o maximum order 5 can be achieved
10 The other parameters to be deined in lssa.in are: Processes Certain random processes can be tried on the series. I an autoregressive (AR) process up to the order 5 is suspected then the user should switch to 1. In this case the coeicient o the AR process along with a statistical test is calculated. The latter indicates on output whether it is signiicant or not. In addition a statistical test is perormed on the calculated coeicient to speciy whether it is dierent rom unity (random walk). I the coeicient is statistically equal to unity then it is a random walk. Characteristics o series This block gives inormation about the covariance o the series. Only a diagonal covariance matrix is accepted. I the a-priori variance actor o the input series is known the user should switch to 1. I so, the values o the series are considered equally weighted and the weight is set to unity. Otherwise the weight o each value o the series is calculated as the inverse o the variance. As mentioned above the standard deviations o the time series are given in the 3 rd column o the input series. Spectrum characteristics The user can deine the output spectrum, as well as the number o spectral values in each band. The user should keep in mind that the spectrum is band limited between the undamental requency and the Nyquist requency. Statistics The user can deine the critical level or statistical testing. Usually values =0.05 or=0.01 are used or critical levels to detect signiicant peaks in spectrum
11 Ater executing LSSA using the input ile lssa.in with properly deined parameters, several output iles are created. These are: lssa.out, residual.dat, spectrum.dat and hist.dat. Lssa.out contains the summary o results together with the statistical tests. The spectrum is described in three dierent orms: percentage variance (least-squares spectrum), power spectral density in db (what is identical with the Fourier transorm when the series is equally spaced), and power spectral density in unit 2 /requency (where unit is the unit o the time series e.g. cycles/day etc.). The resulting spectrum is given in six columns which list the period o spectrum in units o the time series, requency in cycles/unit time, idelity (in units o time). I a signiicant peak is present in the spectrum it shows up in this column. The 4 th column stands or the least-squares spectrum in percentage variance, next ollows with the power spectral density in db and last, the power spectral density in unit 2 /requency is given. The spectrum is also printed with asterisks or quick identiication o peaks Residual.dat lists the input and residual series. There are our columns which contain the time, the input time series, normalized residual series (ater removing trends, periods, processes, etc.) and the standard deviation o the residuals
12 Spectrum.dat contains the aorementioned 6 columns o the output spectrum as described in lssa.out. This ile can be easily exported into Matlab or other mathematical package used or plotting or any other urther analysis. Hist.dat contains the histogram o the normalized residual series
The achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More informationAH 2700A. Attenuator Pair Ratio for C vs Frequency. Option-E 50 Hz-20 khz Ultra-precision Capacitance/Loss Bridge
0 E ttenuator Pair Ratio or vs requency NEEN-ERLN 700 Option-E 0-0 k Ultra-precision apacitance/loss ridge ttenuator Ratio Pair Uncertainty o in ppm or ll Usable Pairs o Taps 0 0 0. 0. 0. 07/08/0 E E E
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationIn many diverse fields physical data is collected or analysed as Fourier components.
1. Fourier Methods In many diverse ields physical data is collected or analysed as Fourier components. In this section we briely discuss the mathematics o Fourier series and Fourier transorms. 1. Fourier
More informationObjectives. By the time the student is finished with this section of the workbook, he/she should be able
FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationLab on Taylor Polynomials. This Lab is accompanied by an Answer Sheet that you are to complete and turn in to your instructor.
Lab on Taylor Polynomials This Lab is accompanied by an Answer Sheet that you are to complete and turn in to your instructor. In this Lab we will approimate complicated unctions by simple unctions. The
More informationAnalog Computing Technique
Analog Computing Technique by obert Paz Chapter Programming Principles and Techniques. Analog Computers and Simulation An analog computer can be used to solve various types o problems. It solves them in
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More informationAdditional exercises in Stationary Stochastic Processes
Mathematical Statistics, Centre or Mathematical Sciences Lund University Additional exercises 8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationCOMP 408/508. Computer Vision Fall 2017 PCA for Recognition
COMP 408/508 Computer Vision Fall 07 PCA or Recognition Recall: Color Gradient by PCA v λ ( G G, ) x x x R R v, v : eigenvectors o D D with v ^v (, ) x x λ, λ : eigenvalues o D D with λ >λ v λ ( B B, )
More informationThe Deutsch-Jozsa Problem: De-quantization and entanglement
The Deutsch-Jozsa Problem: De-quantization and entanglement Alastair A. Abbott Department o Computer Science University o Auckland, New Zealand May 31, 009 Abstract The Deustch-Jozsa problem is one o the
More informationTLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.
TLT-5/56 COMMUNICATION THEORY, Exercise 3, Fall Problem. The "random walk" was modelled as a random sequence [ n] where W[i] are binary i.i.d. random variables with P[W[i] = s] = p (orward step with probability
More informationIntroduction to Analog And Digital Communications
Introduction to Analog And Digital Communications Second Edition Simon Haykin, Michael Moher Chapter Fourier Representation o Signals and Systems.1 The Fourier Transorm. Properties o the Fourier Transorm.3
More informationCHAPTER 8 ANALYSIS OF AVERAGE SQUARED DIFFERENCE SURFACES
CAPTER 8 ANALYSS O AVERAGE SQUARED DERENCE SURACES n Chapters 5, 6, and 7, the Spectral it algorithm was used to estimate both scatterer size and total attenuation rom the backscattered waveorms by minimizing
More information10. Joint Moments and Joint Characteristic Functions
10. Joint Moments and Joint Characteristic Functions Following section 6, in this section we shall introduce various parameters to compactly represent the inormation contained in the joint p.d. o two r.vs.
More informationRATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions
RATIONAL FUNCTIONS Finding Asymptotes..347 The Domain....350 Finding Intercepts.....35 Graphing Rational Functions... 35 345 Objectives The ollowing is a list o objectives or this section o the workbook.
More information9.1 The Square Root Function
Section 9.1 The Square Root Function 869 9.1 The Square Root Function In this section we turn our attention to the square root unction, the unction deined b the equation () =. (1) We begin the section
More informationAsymptote. 2 Problems 2 Methods
Asymptote Problems Methods Problems Assume we have the ollowing transer unction which has a zero at =, a pole at = and a pole at =. We are going to look at two problems: problem is where >> and problem
More informationRoberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points
Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More informationAn Ensemble Kalman Smoother for Nonlinear Dynamics
1852 MONTHLY WEATHER REVIEW VOLUME 128 An Ensemble Kalman Smoother or Nonlinear Dynamics GEIR EVENSEN Nansen Environmental and Remote Sensing Center, Bergen, Norway PETER JAN VAN LEEUWEN Institute or Marine
More information2 Frequency-Domain Analysis
2 requency-domain Analysis Electrical engineers live in the two worlds, so to speak, o time and requency. requency-domain analysis is an extremely valuable tool to the communications engineer, more so
More informationProbabilistic Model of Error in Fixed-Point Arithmetic Gaussian Pyramid
Probabilistic Model o Error in Fixed-Point Arithmetic Gaussian Pyramid Antoine Méler John A. Ruiz-Hernandez James L. Crowley INRIA Grenoble - Rhône-Alpes 655 avenue de l Europe 38 334 Saint Ismier Cedex
More informationPower Spectral Analysis of Elementary Cellular Automata
Power Spectral Analysis o Elementary Cellular Automata Shigeru Ninagawa Division o Inormation and Computer Science, Kanazawa Institute o Technology, 7- Ohgigaoka, Nonoichi, Ishikawa 92-850, Japan Spectral
More informationReceived: 30 July 2017; Accepted: 29 September 2017; Published: 8 October 2017
mathematics Article Least-Squares Solution o Linear Dierential Equations Daniele Mortari ID Aerospace Engineering, Texas A&M University, College Station, TX 77843, USA; mortari@tamu.edu; Tel.: +1-979-845-734
More informationOn GNSS Ambiguity Acceptance Tests
International Global Navigation Satellite Systems Society IGNSS Symposium 27 The University o New South Wales, Sydney, Australia 4 6 December, 27 On GNSS Ambiguity Acceptance Tests P.J.G. Teunissen Delt
More informationShaker Rigs revisited
Shaker Rigs revisited One o the irst articles I ever wrote or Racecar Engineering was on 4 post and 7 post shaker rigs. That article was over 5 years a go and a lot o water has lowed under the bridge and
More informationAcoustic forcing of flexural waves and acoustic fields for a thin plate in a fluid
Acoustic orcing o leural waves and acoustic ields or a thin plate in a luid Darryl MCMAHON Maritime Division, Deence Science and Technology Organisation, HMAS Stirling, WA Australia ABSTACT Consistency
More informationIntrinsic Small-Signal Equivalent Circuit of GaAs MESFET s
Intrinsic Small-Signal Equivalent Circuit o GaAs MESFET s M KAMECHE *, M FEHAM M MELIANI, N BENAHMED, S DALI * National Centre o Space Techniques, Algeria Telecom Laboratory, University o Tlemcen, Algeria
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationBiplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 6 Offprint
Biplots in Practice MICHAEL GREENACRE Proessor o Statistics at the Pompeu Fabra University Chapter 6 Oprint Principal Component Analysis Biplots First published: September 010 ISBN: 978-84-93846-8-6 Supporting
More informationA Comparative Study of Non-separable Wavelet and Tensor-product. Wavelet; Image Compression
Copyright c 007 Tech Science Press CMES, vol., no., pp.91-96, 007 A Comparative Study o Non-separable Wavelet and Tensor-product Wavelet in Image Compression Jun Zhang 1 Abstract: The most commonly used
More informationAerodynamic Admittance Function of Tall Buildings
Aerodynamic Admittance Function o Tall Buildings in hou a Ahsan Kareem b a alou Engineering Int l, Inc., 75 W. Campbell Rd, Richardson, T, USA b Nataz odeling Laboratory, Uniersity o Notre Dame, Notre
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS
ROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS I. Thirunavukkarasu 1, V. I. George 1, G. Saravana Kumar 1 and A. Ramakalyan 2 1 Department o Instrumentation
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationOPTIMAL PLACEMENT AND UTILIZATION OF PHASOR MEASUREMENTS FOR STATE ESTIMATION
OPTIMAL PLACEMENT AND UTILIZATION OF PHASOR MEASUREMENTS FOR STATE ESTIMATION Xu Bei, Yeo Jun Yoon and Ali Abur Teas A&M University College Station, Teas, U.S.A. abur@ee.tamu.edu Abstract This paper presents
More informationApplication of Wavelet Transform Modulus Maxima in Raman Distributed Temperature Sensors
PHOTONIC SENSORS / Vol. 4, No. 2, 2014: 142 146 Application o Wavelet Transorm Modulus Maxima in Raman Distributed Temperature Sensors Zongliang WANG, Jun CHANG *, Sasa ZHANG, Sha LUO, Cuanwu JIA, Boning
More informationBayesian Technique for Reducing Uncertainty in Fatigue Failure Model
9IDM- Bayesian Technique or Reducing Uncertainty in Fatigue Failure Model Sriram Pattabhiraman and Nam H. Kim University o Florida, Gainesville, FL, 36 Copyright 8 SAE International ABSTRACT In this paper,
More informationFigure 3.1 Effect on frequency spectrum of increasing period T 0. Consider the amplitude spectrum of a periodic waveform as shown in Figure 3.2.
3. Fourier ransorm From Fourier Series to Fourier ransorm [, 2] In communication systems, we oten deal with non-periodic signals. An extension o the time-requency relationship to a non-periodic signal
More informationarxiv:quant-ph/ v2 12 Jan 2006
Quantum Inormation and Computation, Vol., No. (25) c Rinton Press A low-map model or analyzing pseudothresholds in ault-tolerant quantum computing arxiv:quant-ph/58176v2 12 Jan 26 Krysta M. Svore Columbia
More informationGaussian Process Regression Models for Predicting Stock Trends
Gaussian Process Regression Models or Predicting Stock Trends M. Todd Farrell Andrew Correa December 5, 7 Introduction Historical stock price data is a massive amount o time-series data with little-to-no
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationLecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters
1 Lecture Outline Basics o Spatial Filtering Smoothing Spatial Filters Averaging ilters Order-Statistics ilters Sharpening Spatial Filters Laplacian ilters High-boost ilters Gradient Masks Combining Spatial
More informationEstimation and detection of a periodic signal
Estimation and detection o a periodic signal Daniel Aronsson, Erik Björnemo, Mathias Johansson Signals and Systems Group, Uppsala University, Sweden, e-mail: Daniel.Aronsson,Erik.Bjornemo,Mathias.Johansson}@Angstrom.uu.se
More informationFigure 1: Bayesian network for problem 1. P (A = t) = 0.3 (1) P (C = t) = 0.6 (2) Table 1: CPTs for problem 1. (c) P (E B) C P (D = t) f 0.9 t 0.
Probabilistic Artiicial Intelligence Problem Set 3 Oct 27, 2017 1. Variable elimination In this exercise you will use variable elimination to perorm inerence on a bayesian network. Consider the network
More informationA Systematic Approach to Frequency Compensation of the Voltage Loop in Boost PFC Pre- regulators.
A Systematic Approach to Frequency Compensation o the Voltage Loop in oost PFC Pre- regulators. Claudio Adragna, STMicroelectronics, Italy Abstract Venable s -actor method is a systematic procedure that
More informationReferences Ideal Nyquist Channel and Raised Cosine Spectrum Chapter 4.5, 4.11, S. Haykin, Communication Systems, Wiley.
Baseand Data Transmission III Reerences Ideal yquist Channel and Raised Cosine Spectrum Chapter 4.5, 4., S. Haykin, Communication Systems, iley. Equalization Chapter 9., F. G. Stremler, Communication Systems,
More informationOBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL
OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell
More informationCORRESPONDENCE ANALYSIS
CORRESPONDENCE ANALYSIS INTUITIVE THEORETICAL PRESENTATION BASIC RATIONALE DATA PREPARATION INITIAL TRANSFORAMATION OF THE INPUT MATRIX INTO PROFILES DEFINITION OF GEOMETRIC CONCEPTS (MASS, DISTANCE AND
More informationUNCERTAINTY EVALUATION OF SINUSOIDAL FORCE MEASUREMENT
XXI IMEKO World Congress Measurement in Research and Industry August 30 eptember 4, 05, Prague, Czech Republic UNCERTAINTY EVALUATION OF INUOIDAL FORCE MEAUREMENT Christian chlegel, Gabriela Kiekenap,Rol
More informationNew Functions from Old Functions
.3 New Functions rom Old Functions In this section we start with the basic unctions we discussed in Section. and obtain new unctions b shiting, stretching, and relecting their graphs. We also show how
More informationSHORT COMMUNICATION. Diversity partitions in 3-way sorting: functions, Venn diagram mappings, typical additive series, and examples
COMMUNITY ECOLOGY 7(): 53-57, 006 585-8553/$0.00 Akadémiai Kiadó, Budapest DOI: 0.556/ComEc.7.006.. SHORT COMMUNICATION Diversity partitions in 3-way sorting: unctions, Venn diagram mappings, typical additive
More informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Algorithm. partitioning. procedure
Development o a Mathematical Algorithm Transormation Process Mathematical Model partitioning procedure Mathematical Algorithm The mathematical algorithm is a plan (or blueprint) that identiies the independent
More informationLight, Quantum Mechanics and the Atom
Light, Quantum Mechanics and the Atom Light Light is something that is amiliar to most us. It is the way in which we are able to see the world around us. Light can be thought as a wave and, as a consequence,
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationA UNIFIED FRAMEWORK FOR MULTICHANNEL FAST QRD-LS ADAPTIVE FILTERS BASED ON BACKWARD PREDICTION ERRORS
A UNIFIED FRAMEWORK FOR MULTICHANNEL FAST QRD-LS ADAPTIVE FILTERS BASED ON BACKWARD PREDICTION ERRORS César A Medina S,José A Apolinário Jr y, and Marcio G Siqueira IME Department o Electrical Engineering
More informationDigital Image Processing. Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009
Digital Image Processing Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 009 Outline Image Enhancement in Spatial Domain Spatial Filtering Smoothing Filters Median Filter
More informationNON-EQUILIBRIUM REACTION RATES IN THE MACROSCOPIC CHEMISTRY METHOD FOR DSMC CALCULATIONS. March 19, 2007
NON-EQUILIBRIUM REACTION RATES IN THE MACROSCOPIC CHEMISTRY METHOD FOR DSMC CALCULATIONS M.J. GOLDSWORTHY, M.N. MACROSSAN & M.M ABDEL-JAWAD March 9, 7 The Direct Simulation Monte Carlo (DSMC) method is
More informationSTAT 801: Mathematical Statistics. Hypothesis Testing
STAT 801: Mathematical Statistics Hypothesis Testing Hypothesis testing: a statistical problem where you must choose, on the basis o data X, between two alternatives. We ormalize this as the problem o
More informationSPOC: An Innovative Beamforming Method
SPOC: An Innovative Beamorming Method Benjamin Shapo General Dynamics Ann Arbor, MI ben.shapo@gd-ais.com Roy Bethel The MITRE Corporation McLean, VA rbethel@mitre.org ABSTRACT The purpose o a radar or
More informationTFY4102 Exam Fall 2015
FY40 Eam Fall 05 Short answer (4 points each) ) Bernoulli's equation relating luid low and pressure is based on a) conservation o momentum b) conservation o energy c) conservation o mass along the low
More informationSliding Mode Control and Feedback Linearization for Non-regular Systems
Sliding Mode Control and Feedback Linearization or Non-regular Systems Fu Zhang Benito R. Fernández Pieter J. Mosterman Timothy Josserand The MathWorks, Inc. Natick, MA, 01760 University o Texas at Austin,
More informationEnhancement Using Local Histogram
Enhancement Using Local Histogram Used to enhance details over small portions o the image. Deine a square or rectangular neighborhood hose center moves rom piel to piel. Compute local histogram based on
More informationDiscrete-Time Fourier Transform (DTFT)
Connexions module: m047 Discrete-Time Fourier Transorm DTFT) Don Johnson This work is produced by The Connexions Project and licensed under the Creative Commons Attribution License Abstract Discussion
More informationThe second-order diffraction-radiation for a piston-like arrangement
The second-order diraction-radiation or a piston-like arrangement I.K. Chatjigeorgiou & S.A. Mavrakos Laboratory o Floating Structures and Mooring Systems, School o Naval Architecture and Marine Engineering,
More informationSymbolic-Numeric Methods for Improving Structural Analysis of DAEs
Symbolic-Numeric Methods or Improving Structural Analysis o DAEs Guangning Tan, Nedialko S. Nedialkov, and John D. Pryce Abstract Systems o dierential-algebraic equations (DAEs) are generated routinely
More informationRandom Error Analysis of Inertial Sensors output Based on Allan Variance Shaochen Li1, a, Xiaojing Du2,b and Junyi Zhai3,c
International Conerence on Civil, Transportation and Environment (ICCTE 06) Random Error Analysis o Inertial Sensors output Based on Allan Variance Shaochen Li, a, Xiaojing Du, and Junyi Zhai3,c School
More informationA Fourier Transform Model in Excel #1
A Fourier Transorm Model in Ecel # -This is a tutorial about the implementation o a Fourier transorm in Ecel. This irst part goes over adjustments in the general Fourier transorm ormula to be applicable
More informationMEASUREMENT UNCERTAINTIES
MEASUREMENT UNCERTAINTIES What distinguishes science rom philosophy is that it is grounded in experimental observations. These observations are most striking when they take the orm o a quantitative measurement.
More information3. Several Random Variables
. Several Random Variables. Two Random Variables. Conditional Probabilit--Revisited. Statistical Independence.4 Correlation between Random Variables. Densit unction o the Sum o Two Random Variables. Probabilit
More informationEstimation of the precision by the tidal analysis programs ETERNA and VAV
11189 Estimation o the precision by the tidal analysis programs and B. Ducarme*, L. Vandercoilden*, A.P. Venedikov** * Royal Observatory o Belgium, ducarme@oma.be ** Geophysical Institute & Central Laboratory
More informationProduct Matrix MSR Codes with Bandwidth Adaptive Exact Repair
1 Product Matrix MSR Codes with Bandwidth Adaptive Exact Repair Kaveh Mahdaviani, Soheil Mohajer, and Ashish Khisti ECE Dept, University o Toronto, Toronto, ON M5S3G4, Canada ECE Dept, University o Minnesota,
More information9.3 Graphing Functions by Plotting Points, The Domain and Range of Functions
9. Graphing Functions by Plotting Points, The Domain and Range o Functions Now that we have a basic idea o what unctions are and how to deal with them, we would like to start talking about the graph o
More informationATP-EMTP evaluation of relaying algorithms for series-compensated line
AP-EMP evaluation o relaying algorithms or series-compensated line M. M. Saha, E. osolowski, J. Izykowski, P. Pierz Abstract A group o the algorithms aimed at measuring a ault-loop impedance or digital
More informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationSimulation of earthquake ground motion for generation of artificial accelerograms
Earthquake Resistant Engineering Structures V 83 Simulation o earthquake ground motion or generation o artiicial accelerograms J. A. Abdalla & Y. M. ag-elhassan Department o Civil Engineering, American
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 19 Modeling Topics plan: Modeling (linear/non- linear least squares) Bayesian inference Bayesian approaches to spectral esbmabon;
More informationStanding Waves If the same type of waves move through a common region and their frequencies, f, are the same then so are their wavelengths, λ.
Standing Waves I the same type o waves move through a common region and their requencies,, are the same then so are their wavelengths,. This ollows rom: v=. Since the waves move through a common region,
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationChapter 4 Imaging. Lecture 21. d (110) Chem 793, Fall 2011, L. Ma
Chapter 4 Imaging Lecture 21 d (110) Imaging Imaging in the TEM Diraction Contrast in TEM Image HRTEM (High Resolution Transmission Electron Microscopy) Imaging or phase contrast imaging STEM imaging a
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationA Brief Survey on Semi-supervised Learning with Graph Regularization
000 00 002 003 004 005 006 007 008 009 00 0 02 03 04 05 06 07 08 09 020 02 022 023 024 025 026 027 028 029 030 03 032 033 034 035 036 037 038 039 040 04 042 043 044 045 046 047 048 049 050 05 052 053 A
More informationTHE use of radio frequency channels assigned to primary. Traffic-Aware Channel Sensing Order in Dynamic Spectrum Access Networks
EEE JOURNAL ON SELECTED AREAS N COMMUNCATONS, VOL. X, NO. X, X 01X 1 Traic-Aware Channel Sensing Order in Dynamic Spectrum Access Networks Chun-Hao Liu, Jason A. Tran, Student Member, EEE, Przemysław Pawełczak,
More informationProducts and Convolutions of Gaussian Probability Density Functions
Tina Memo No. 003-003 Internal Report Products and Convolutions o Gaussian Probability Density Functions P.A. Bromiley Last updated / 9 / 03 Imaging Science and Biomedical Engineering Division, Medical
More informationDetecting Doubly Compressed JPEG Images by Factor Histogram
APSIPA ASC 2 Xi an Detecting Doubly Compressed JPEG Images by Factor Histogram Jianquan Yang a, Guopu Zhu a, and Jiwu Huang b a Shenzhen Institutes o Advanced Technology, Chinese Academy o Sciences, Shenzhen,
More informationMISS DISTANCE GENERALIZED VARIANCE NON-CENTRAL CHI DISTRIBUTION. Ken Chan ABSTRACT
MISS DISTANCE GENERALIZED VARIANCE NON-CENTRAL CI DISTRIBUTION en Chan E-Mail: ChanAerospace@verizon.net ABSTRACT In many current practical applications, the probability o collision is oten considered
More informationMath Review and Lessons in Calculus
Math Review and Lessons in Calculus Agenda Rules o Eponents Functions Inverses Limits Calculus Rules o Eponents 0 Zero Eponent Rule a * b ab Product Rule * 3 5 a / b a-b Quotient Rule 5 / 3 -a / a Negative
More informationDesign criteria for Fiber Reinforced Rubber Bearings
Design criteria or Fiber Reinorced Rubber Bearings J. M. Kelly Earthquake Engineering Research Center University o Caliornia, Berkeley A. Calabrese & G. Serino Department o Structural Engineering University
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationMath-Essentials Unit 3 Review. Equations and Transformations of the Linear, Quadratic, Absolute Value, Square Root, and Cube Functions
Math-Essentials Unit Review Equations and Transormations o the Linear, Quadratic, Absolute Value, Square Root, and Cube Functions Vocabulary Relation: A mapping or pairing o input values to output values.
More informationL8. LINEAR PANEL DATA MODELS UNDER STRICT AND WEAK EXOGENEITY. 1. A Structural Linear Panel Model
Victor Chernozhukov and Ivan Fernandez-Val. 14.382 Econometrics. Spring 2017. Massachusetts Institute o Technology: MIT OpenCourseWare, https://ocw.mit.edu. License: Creative Commons BY-NC-SA. 14.382 L8.
More informationINPUT GROUND MOTION SELECTION FOR XIAOWAN HIGH ARCH DAM
3 th World Conerence on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 2633 INPUT GROUND MOTION LECTION FOR XIAOWAN HIGH ARCH DAM CHEN HOUQUN, LI MIN 2, ZHANG BAIYAN 3 SUMMARY In
More informationAnalytical expressions for field astigmatism in decentered two mirror telescopes and application to the collimation of the ESO VLT
ASTRONOMY & ASTROPHYSICS MAY II 000, PAGE 57 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 44, 57 67 000) Analytical expressions or ield astigmatism in decentered two mirror telescopes and application
More informationMicro-canonical ensemble model of particles obeying Bose-Einstein and Fermi-Dirac statistics
Indian Journal o Pure & Applied Physics Vol. 4, October 004, pp. 749-757 Micro-canonical ensemble model o particles obeying Bose-Einstein and Fermi-Dirac statistics Y K Ayodo, K M Khanna & T W Sakwa Department
More informationNONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION
NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION Steven Ball Science Applications International Corporation Columbia, MD email: sball@nmtedu Steve Schaer Department o Mathematics
More information