Introduction to Least Squares Adjustment for geodetic VLBI
|
|
- August Chandler
- 5 years ago
- Views:
Transcription
1 Introduction to Least Squares Adjustment for geodetic VLBI Matthias Schartner a, David Mayer a a TU Wien, Department of Geodesy and Geoinformation
2 Least Squares Adjustment why? observation is τ (baseline) unknown parameters are: station positions EOP source positions troposphere... How do we link our observations to our unknown parameters? AVN school 1/ 19 (45)
3 A very simple example three sample points data fitting polynome of order 2: y = ax 2 + bx + c observations: (x i, y i ) unknowns: a, b, c AVN school 2/ 19 (45)
4 A very simple example (x 1, y 1 ) (x 2, y 2 ) (x 3, y 3 ) y 1 = ax1 2 + bx 1 + c y 2 = ax2 2 + bx 2 + c y 3 = ax3 2 + bx 3 + c Or in matrix notation: y 1 x y 2 = 1 2 x 1 1 a x2 2 x 2 1 b x3 2 x 3 1 c y AVN school 3/ 19 (45)
5 A very simple example y 1 x y 2 = 1 2 x 1 1 a x2 2 x 2 1 b y 3 x3 2 x 3 1 c solve the equations a x b = 1 2 x y 1 x2 2 x 2 1 y 2 c x3 2 x 3 1 y 3 coefficients for y = ax 2 + bx + c AVN school 4/ 19 (45)
6 A very simple example Warning! Every observation has an error How accurate are our coefficients a, b and c? What can we do? If we have more observations we can average out the error AVN school 5/ 19 (45)
7 A very simple example What happens if I have more observations? y = ax 2 + bx + c three unknowns a, b, c (n unk ) many observations (n obs ) n obs > n unk How do we get our unknown parameters? AVN school 6/ 19 (45)
8 A very simple example simply pick three: (x 1, y 1 ) (x 2, y 2 ) (x 3, y 3 ) y 1 = ax1 2 + bx 1 + c y 2 = ax2 2 + bx 2 + c y 3 = ax3 2 + bx 3 + c Or in matrix notation: y 1 x y 2 = 1 2 x 1 1 a x2 2 x 2 1 b x3 2 x 3 1 c y AVN school 7/ 19 (45)
9 A very simple example y 1 x y 2 = 1 2 x 1 1 a x2 2 x 2 1 b y 3 x3 2 x 3 1 c solve the equations a x b = 1 2 x y 1 x2 2 x 2 1 y 2 c x3 2 x 3 1 y 3 coefficients for y = ax 2 + bx + c AVN school 8/ 19 (45)
10 A very simple example Why not this three? AVN school 9/ 19 (45)
11 A very simple example 50 or this three: Our goal: we want to use every available observations AVN school 10/ 19 (45)
12 A very simple example y 1 x 2 1 x 1 1 y 2 x 2 = 2 x 2 1 a b.... c y n xn 2 x n 1 }{{}}{{}}{{} l A x A R n obs n unk x R n unk l R n obs A x = l (1) AVN school 11/ 19 (45)
13 A very simple example trick and solution A A x = A l (2) ( 1 x = A A) A l (3) residuals v: v i = ax 2 i + bx i + c y i v = A x l (4) minimization v 2 = v v (5) AVN school 12/ 19 (45)
14 A very simple example AVN school 13/ 19 (45)
15 A very simple example conclusion / goal get the most reasonable values for our unknown parameters take all available information more observations then unknown parameters get some error margins How can we do this? Warning! every observation has an error AVN school 14/ 19 (45)
16 Least Squares Methode pros easy while powerful minimize residuals rigorous statistical approach individual weighting (28) accuracy information (32) manipulate geodetic datum (37) combine results of various LSMs (42) cons only for linear functions liniarization via Taylor polynome (22) a priori values iterations (31) we need to model our observations o c outliers requirements for LSM partial derivatives (24) AVN school 15/ 19 (45)
17 What do we want to estimate? station coordinates source coordinates earth orientation parameters tropospheric parameters clock behavior and many more... Many unknown parameters AVN school 16/ 19 (45)
18 What do we want to estimate? Example: Scan with 5 stations = 10 observations Estimates: Troposphere: (15p.) zenith wet delay (5p.) north gradient (5p.) east gradient (5p.) clock (5p.) Parameters Way more unknown parameters then observations Requirement for LSM We need more observations then unknowns! n obs > n unk AVN school 17/ 19 (45)
19 What do we want to estimate? Solution piecewise linear offsets (34) we estimate parameters only every x minutes linear interpolation in between reduces number of estimated parameters! Improvement constraints (35) zwd sta1,03:00 = zwd sta1,04:00 ±1.5 [cm] }{{} constraint AVN school 18/ 19 (45)
20 Conclusion What should you remeber? we calculate our solutions using Least Squares Adjustment we minimize (squared sum of) residuals we need to model our observations! theoretical delay o c we need partial derivatives for our Jacobian matrix A we need a priori values X 0 (usually not a big deal) we estimate additions x to a priori values X 0 we need to eliminate outliers AVN school 19/ 19 (45)
21 Example P2 P3 5 stations 3 with known coordinates 2 with unknown coordinates P5 P4 P AVN school 20/ 19 (45)
22 Example 2 measured distance s ij unknowns: (x4, y4) (x5, y5) n unk = 4 observations: s ij n obs = 7 A R 7 4 x R 4 l R 7 n obs > n unk P2 P3 P5 P4 P AVN school 21/ 19 (45)
23 Some basics about LSM We need to model our system: s ij = (x j x i ) 2 + (y j y i ) 2 linearization not linear Taylor expansion: T,x0 (x) = f (x 0 ) + f (x 0 ) 1! = n=0 f (n) (x 0 ) (x x 0 ) n n! liniarized observation equation What do we need? We need partial derivative We need a priori values x 0 (x x 0 ) + f (x 0 ) (x x 0 ) ! AVN school 22/ 19 (45)
24 Some basics about LSM definitions observations: L (computed) approximated observations: L 0 reduced observation vector (observed minus computed): l l = L L 0 (6) estimated unknowns: ˆX approximated unknowns: X 0 reduced unknowns: x x = ˆX X 0 (7) Designmatrix or Jacobian matrix: A (partials) AVN school 23/ 19 (45)
25 Back to Example 2 unknown: x 4, y 4, x 5, y 5 known: x 1, y 1, x 2, y 2, x 3, y 3 approximate values: x 4,0, y 4,0, x 5,0, y 5,0 s ij = (x j x i ) 2 + (y j y i ) 2 ( ) s14 ( x 4 ) s15 x 4 A = (. ) s45 x 4 ( ) s14 y 4 ( ) s15 y 4 (. ) s45 y 4 ( ) s14 x 5 ( ) s15 x 5 (. ) s45 x 5 ( ) s14 ( y 5 ) s15 y 5 (. ) s45 y 5 s 14 x 4 = x 4,0 x 1 ( x4,0 x1) 2 + ( y 4,0 y 1 ) 2 s 14 y 4 = y 4,0 y 1 ( x4,0 x1) 2 + ( y 4,0 y 1 ) 2 s 14 x 5 = s 14 y 5 = AVN school 24/ 19 (45)
26 Example 2 unknown: x 4, y 4, x 5, y 5 known: x 1, y 1, x 2, y 2, x 3, y 3 approximate values: x 4,0, y 4,0, x 5,0, y 5,0 reduced observation vector (observed minus computed) l = L L 0 s 14 s 14,0 s 15 s 15,0 l =. s 45 s 45,0 s 14,0 = (x 4,0 x 1 ) 2 + (y 4,0 y 1 ) AVN school 25/ 19 (45)
27 Example 2 reduced observation vector (observed minus computed) l = L L 0 Jacobian matrix (partials): A reduced unknowns: x = X X 0 trick Ax = l A Ax = A l ( 1 x = A A) A l (8) What do we estimate? We don t estimate whole unknown parameter ˆX but only additions x to some a priori values X AVN school 26/ 19 (45)
28 Example The coordinates of our unknown points change P5(ˆx 5, ŷ 5 ) P4(ˆx 4, ŷ 4 ) AVN school 27/ 19 (45)
29 Some more basics about LSM usually not every observations is equally accurate: Covariance matrix Q ll R n obs n obs, Weight matrix P: σ 2 l 1 σ l1,l 2... σ l1,l n σ l1,l Q ll = 2 σ 2 l 2... σ l2,l n σ l1,l n σ l2,l n... σl 2 n if observations are independent σ li,l j = 0 σ 2 l σ 2 Q ll = l P = 1. s σl 2 n covariance matrix }{{} usually1 Q 1 ll (9) AVN school 28/ 19 (45)
30 Some more basics about LSM covariance matrix Variance-Covariance matrix Q ll, Weight matrix P R n obs n obs : P = 1 s0 2 Qll 1 = }{{} usually σl σl σln 2 good observation small standard deviation σ li and variance high weight in P matrix σ 2 l i (10) AVN school 29/ 19 (45)
31 Some more basics about LSM Final formulas: A} {{ PA} x = A} {{ Pl } N b equations (11) x = N 1 b (12) Important for LSM Normal equation matrix N = right hand side b = A Pl we minimize v Pv ( ) A PA b R n unk N R n unk n unk ( 1 ˆX = X 0 + A PA) A Pl = X 0 + N 1 b (13) AVN school 30/ 19 (45)
32 Some more basics about LSM iteration Sometimes a iteration is necessary. inaccurate a priori values X 0 numerical reasons You can use ˆX again as new a priori values X 0 and redo everything X 0,n+1 = ˆX n = X 0,n + N 1 n b n (14) AVN school 31/ 19 (45)
33 How good are our estimates? a posteriori variance factor σ 2 0 : covariance matrix σ0 2 = v Pv (15) n obs n unk remember our residuals v are: v = Ax l Variance-Covariance Matrix Q xx R n unk n unk σ 2 x 1 σ x1,x 2... σ x1,x n Q xx = σ0n 2 1 σ x1,x = 2 σ 2 x 2... σ x2,x n (16) σ x1,x n σ x2,x n... σ 2 x n σx 2 i σ xi is variance of unknown parameter x i is standard deviation AVN school 32/ 19 (45)
34 Back to example Now we know how accurate our coordinates are P5(ˆx 5, ŷ 5 ) P4(ˆx 4, ŷ 4 ) AVN school 33/ 19 (45)
35 Concepts What are piecewise linear offsets? time dependent parameters Example: Troposphere (mf = mapping function) x 2 x 1 x L w (t) = mf w (t)x 1 + mf w (t) t t 1 (x 2 x 1 ) t 2 t 1 L w = mf w (t) mf w (t) t t 1 x 1 t 2 t 1 L w = mf w (t) t t 1 x 2 t 2 t 1 t 1 t t AVN school 34/ 19 (45)
36 Concepts What are constraints? pseudo observations help us with singularities useful concept A simple example: (think of troposphere between two time epochs) A c R nconst n unk l c R nconst A c = x j x i = 0 ± 5 x j 1 x i 1 = 0 ± 3 (± means standard deviation) ( ) ( ) 0 l c = 0 P c R nconst nconst ( 1 ) P c = AVN school / 19 (45)
37 Concepts What are constraints? N c R n unk n unk b c R nconst Combination with the real observations: N c = A c P c A c (17) b c = A c P c l c (18) N total = A PA + A c P c A c = N + N c (19) b total = A Pl + A c P c l c = b + b c (20) The only thing that changes is: σ 2 0,const = v Pv + v c P c v c n obs + n const n unk (21) AVN school 36/ 19 (45)
38 Concepts What is a geodetic datum? - Example 2 Think of example 2 This time no station has fixed coordinates unknowns: all coordinates observations: s ij we have a priori coordinates P2 P3 P5 P4 P AVN school 37/ 19 (45)
39 Concepts What is a geodetic datum? we know inner geometry very well we can not estimate absolute coordinates N matrix is singular is this correct? P2 P1 P4 P5 P AVN school 38/ 19 (45)
40 Concepts What is a geodetic datum? - Example 2 we know inner geometry very well we can not estimate absolute coordinates N matrix is singular or this? P2 P1 P4 P5 P AVN school 39/ 19 (45)
41 Concepts What is a geodetic datum? - Formulas We need datum stations with additional conditions for example in this case we could say: translation in x translation in y rotation around z scale dxi = 0 dyi = 0 (yi dx i x i dy i ) = 0 (xi dx i y i dy i ) = 0 This would lead to the following matrix: G R n datum n unk G = y 1 x 1... y n x n x 1 y 1... x n y n (22) AVN school 40/ 19 (45)
42 Concepts What is a geodetic datum? - Solution The solution is now: ( ) ( ) ( ) A PA G x A G = Pl 0 k 0 (23) What else is different: σ 2 0 = v Pv n obs + n datum n unk ; Geodetic datum in VLBI ( ) Qxx Q xk = σ0 2 Q xk Q kk ( A PA ) 1 G G 0 (24) usually we use NNT/NNR - No Net Rotation and No Net Translation conditions (3 NNT and 3 NNR) scale is fixed by observations AVN school 41/ 19 (45)
43 Concepts What is a global solution? combination of large number of sessions to estimate accurate common parameters combination is done at the normal equation level you need to make sure N matrix has same structure you want to get rid of some parameters ( ) ( ) ( ) N11 N Nx = b 12 x1 b1 = N 21 N 22 x 2 b 2 ( ) N 11 N 12 N22 1 N 21 }{{} N reduc 1 x 1 = b 1 N 12N22 b 2 (25) }{{} b reduc N REDUC = N reduc1 + N reduc N reducn (26) b REDUC = b reduc1 + b reduc b reducn (27) AVN school 42/ 19 (45)
44 Conclusion What is to remeber? we calculate our solutions using Least Squares Adjustment we minimize v Pv we need to model our observations! theoretical delay we need partial derivatives for our Jacobian matrix A we need a priori values X 0 (usually not a big deal) we estimate additions x to a priori values X AVN school 43/ 19 (45)
45 Link to VLBI AVN school 44/ 19 (45)
46 Matthias Schartner a, David Mayer a, a TU Wien, Department of Geodesy and Geoinformation
Introduction to geodetic VLBI
Introduction to geodetic VLBI David Mayer Andreas Hellerschmied Johannes Böhm Harald Schuh and Johannes Böhm, Very Long Baseline Interferometry for Geodesy and Astrometry, in Guochang Xu (editor): Sciences
More informationA priori gradients in the analysis of GPS and VLBI observations
A priori gradients in the analysis of GPS and VLBI observations Peter Steigenberger GeoForschungsZentrum Potsdam Volker Tesmer Deutsche Geodätisches Forschungsinstitut, München Johannes Boehm Institute
More informationVery Long Baseline Interferometry for Geodesy and Astrometry
VieVS User-Workshop 2017 Very Long Baseline Interferometry for Geodesy and Astrometry Johannes Böhm VLBI How does it work? VLBI a flowchart SINEX skd vex2 drudge drudge snp/prc snp/prc NGS Mark4 Mark4
More informationVery Long Baseline Interferometry for Geodesy and Astrometry
Very Long Baseline Interferometry for Geodesy and Astrometry Johannes Böhm Harald Schuh and Johannes Böhm, Very Long Baseline Interferometry for Geodesy and Astrometry, in Guochang Xu (editor): Sciences
More informationFormulation of L 1 Norm Minimization in Gauss-Markov Models
Formulation of L 1 Norm Minimization in Gauss-Markov Models AliReza Amiri-Simkooei 1 Abstract: L 1 norm minimization adjustment is a technique to detect outlier observations in geodetic networks. The usual
More informationMath 096--Quadratic Formula page 1
Math 096--Quadratic Formula page 1 A Quadratic Formula. Use the quadratic formula to solve quadratic equations ax + bx + c = 0 when the equations can t be factored. To use the quadratic formula, the equation
More informationE. Calais Purdue University - EAS Department Civil 3273
E. Calais urdue University - EAS Department Civil 3273 ecalais@purdue.edu Need for a Reference Frame 1. ositions and velocities from geodetic measurements: Are not direct observations, but estimated quantities
More informationGPS Geodesy - LAB 7. Neglecting the propagation, multipath, and receiver errors, eq.(1) becomes:
GPS Geodesy - LAB 7 GPS pseudorange position solution The pseudorange measurements j R i can be modeled as: j R i = j ρ i + c( j δ δ i + ΔI + ΔT + MP + ε (1 t = time of epoch j R i = pseudorange measurement
More informationBasic Linear Algebra in MATLAB
Basic Linear Algebra in MATLAB 9.29 Optional Lecture 2 In the last optional lecture we learned the the basic type in MATLAB is a matrix of double precision floating point numbers. You learned a number
More informationLeast squares: introduction to the network adjustment
Least squares: introduction to the network adjustment Experimental evidence and consequences Observations of the same quantity that have been performed at the highest possible accuracy provide different
More informationCONT08 - First Results and High-Frequency Earth Rotation
. EVGA-Meeting 29 (Bordeaux) CONT8 - First Results and High-Frequency Earth Rotation 1 Thomas Artz Sarah Böckmann Axel Nothnagel Institut of Geodesy and Geoinformation University of Bonn March 24, 29 CONT8
More informationPhysics 403. Segev BenZvi. Parameter Estimation, Correlations, and Error Bars. Department of Physics and Astronomy University of Rochester
Physics 403 Parameter Estimation, Correlations, and Error Bars Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Best Estimates and Reliability
More informationEuropean Sessions and Datum Definitions
, IVS 2012 General Meeting Proceedings, p.271 275 http://ivscc.gsfc.nasa.gov/publications/gm2012/tornatore.pdf Vincenza Tornatore 1, Hana Krásná (neé Spicakova) 2, Damiano Triglione 1 1) Politecnico di
More informationLeast squares problems
Least squares problems How to state and solve them, then evaluate their solutions Stéphane Mottelet Université de Technologie de Compiègne 30 septembre 2016 Stéphane Mottelet (UTC) Least squares 1 / 55
More informationOptimization and Calculus
Optimization and Calculus To begin, there is a close relationship between finding the roots to a function and optimizing a function. In the former case, we solve for x. In the latter, we solve: g(x) =
More informationInformation Dimension
Information Dimension Mina Karzand Massachusetts Institute of Technology November 16, 2011 1 / 26 2 / 26 Let X would be a real-valued random variable. For m N, the m point uniform quantized version of
More informationConvergence in Distribution
Convergence in Distribution Undergraduate version of central limit theorem: if X 1,..., X n are iid from a population with mean µ and standard deviation σ then n 1/2 ( X µ)/σ has approximately a normal
More informationSub-Daily Parameter Estimation in VLBI Data Analysis
Dissertation Sub-Daily Parameter Estimation in VLBI Data Analysis Supervisors o.univ.-prof. Dipl.-Ing. Dr.techn. Harald Schuh Priv.-Doz. Dr.-Ing. Axel Nothnagel Examiners Prof. Dr. Hansjörg Kutterer Dr.
More informationPrinciples of the Global Positioning System Lecture 11
12.540 Principles of the Global Positioning System Lecture 11 Prof. Thomas Herring http://geoweb.mit.edu/~tah/12.540 Statistical approach to estimation Summary Look at estimation from statistical point
More informationIntegration - Past Edexcel Exam Questions
Integration - Past Edexcel Exam Questions 1. (a) Given that y = 5x 2 + 7x + 3, find i. - ii. - (b) ( 1 + 3 ) x 1 x dx. [4] 2. Question 2b - January 2005 2. The gradient of the curve C is given by The point
More informationSimple Linear Regression for the MPG Data
Simple Linear Regression for the MPG Data 2000 2500 3000 3500 15 20 25 30 35 40 45 Wgt MPG What do we do with the data? y i = MPG of i th car x i = Weight of i th car i =1,...,n n = Sample Size Exploratory
More informationCh 4 Differentiation
Ch 1 Partial fractions Ch 6 Integration Ch 2 Coordinate geometry C4 Ch 5 Vectors Ch 3 The binomial expansion Ch 4 Differentiation Chapter 1 Partial fractions We can add (or take away) two fractions only
More informationINCREASING THE RELIABILITY OF GPS GEODETIC NETWORKS
Proceedings of the first international conference on satellite positioning systems, Alexandria, December 12-13, 1995. INCREASING THE RELIABILITY OF GPS GEODETIC NETWORKS ABSTRACT By Dr. Dalal S. Alnaggar
More informationMATH 115, SUMMER 2012 LECTURE 12
MATH 115, SUMMER 2012 LECTURE 12 JAMES MCIVOR - last time - we used hensel s lemma to go from roots of polynomial equations mod p to roots mod p 2, mod p 3, etc. - from there we can use CRT to construct
More informationCS4495/6495 Introduction to Computer Vision. 8B-L2 Principle Component Analysis (and its use in Computer Vision)
CS4495/6495 Introduction to Computer Vision 8B-L2 Principle Component Analysis (and its use in Computer Vision) Wavelength 2 Wavelength 2 Principal Components Principal components are all about the directions
More informationSimulations of VLBI observations with the Onsala Twin Telescope
Simulations of VLBI observations with the Onsala Twin Telescope Master of Science Thesis CAROLINE SCHÖNBERGER Department of Earth and Space Sciences Division of Space Geodesy and Geodynamics CHALMERS UNIVERSITY
More informationComputational Physics
Interpolation, Extrapolation & Polynomial Approximation Lectures based on course notes by Pablo Laguna and Kostas Kokkotas revamped by Deirdre Shoemaker Spring 2014 Introduction In many cases, a function
More informationNumerical Methods for Inverse Kinematics
Numerical Methods for Inverse Kinematics Niels Joubert, UC Berkeley, CS184 2008-11-25 Inverse Kinematics is used to pose models by specifying endpoints of segments rather than individual joint angles.
More informationGaussian and Linear Discriminant Analysis; Multiclass Classification
Gaussian and Linear Discriminant Analysis; Multiclass Classification Professor Ameet Talwalkar Slide Credit: Professor Fei Sha Professor Ameet Talwalkar CS260 Machine Learning Algorithms October 13, 2015
More informationJoint Probability Distributions and Random Samples (Devore Chapter Five)
Joint Probability Distributions and Random Samples (Devore Chapter Five) 1016-345-01: Probability and Statistics for Engineers Spring 2013 Contents 1 Joint Probability Distributions 2 1.1 Two Discrete
More informationIntelligent Systems Discriminative Learning, Neural Networks
Intelligent Systems Discriminative Learning, Neural Networks Carsten Rother, Dmitrij Schlesinger WS2014/2015, Outline 1. Discriminative learning 2. Neurons and linear classifiers: 1) Perceptron-Algorithm
More informationChallenges and perspectives for CRF and TRF determination
Challenges and perspectives for CRF and TRF determination J. Böhm, Z. Malkin, S. Lambert, C. Ma with contributions by H. Spicakova, L. Plank, and H. Schuh Consistency TRF EOP CRF ITRF2008 from VLBI/GNSS/SLR/DORIS
More informationCommon Realization of Terrestrial and Celestial Reference Frame
Common Realization of Terrestrial and Celestial Reference Frame M. Seitz, R. Heinkelmann, P. Steigenberger, T. Artz Abstract The realization of the International Celestial Reference System (ICRS) and the
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More informationImpact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames
Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames J. Böhm, H. Spicakova, L. Urquhart, P. Steigenberger, H. Schuh Abstract We compare the influence of two different a priori gradient
More informationSection Least Squares Regression
Section 2.3 - Least Squares Regression Statistics 104 Autumn 2004 Copyright c 2004 by Mark E. Irwin Regression Correlation gives us a strength of a linear relationship is, but it doesn t tell us what it
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationQ 2.0.2: If it s 5:30pm now, what time will it be in 4753 hours? Q 2.0.3: Today is Wednesday. What day of the week will it be in one year from today?
2 Mod math Modular arithmetic is the math you do when you talk about time on a clock. For example, if it s 9 o clock right now, then it ll be 1 o clock in 4 hours. Clearly, 9 + 4 1 in general. But on a
More informationLinear Regression. In this problem sheet, we consider the problem of linear regression with p predictors and one intercept,
Linear Regression In this problem sheet, we consider the problem of linear regression with p predictors and one intercept, y = Xβ + ɛ, where y t = (y 1,..., y n ) is the column vector of target values,
More informationIntroduction to Statistics and Data Analysis
Introduction to Statistics and Data Analysis RSI 2005 Staff July 15, 2005 Variation and Statistics Good experimental technique often requires repeated measurements of the same quantity These repeatedly
More informationIntroduction: the Abruzzo earthquake The network and the processing strategies. displacements estimation at earthquake epoch. horizontal displacements
The Abruzzo earthquake: temporal and spatial analysis of the first geodetic results L. Biagi, S. Caldera, D. Dominici, F. Sansò Politecnico di Milano Università degli studi de L Aquila Outline Introduction:
More informationLinear Regression Spring 2014
Linear Regression 18.05 Spring 2014 Agenda Fitting curves to bivariate data Measuring the goodness of fit The fit vs. complexity tradeoff Regression to the mean Multiple linear regression January 1, 2017
More informationGeodetic Network Adjustment Examples
Geodetic Network Adjustment Examples Friedhelm Krumm Geodätisches Institut Universität Stuttgart http//www.gis.uni-stuttgart.de krumm@gis.uni-stuttgart.de Rev. 3.1 June 13, 2018 Contents 1 Introduction
More informationAT760 Global Carbon Cycle. Assignment #3 Due Friday, May 4, 2007 Atmospheric Transport and Inverse Modeling of CO 2
AT760 Global Carbon Cycle Assignment 3 Due Friday, May 4, 2007 Atmospheric Transport and Inverse Modeling of CO 2 In this exercise you will develop a very simplified model of the mixing of the global atmosphere.
More informationA. Barbu, J. Laurent-Varin, F. Perosanz, F. Mercier and J. Marty. AVENUE project. June, 20
Efficient QR Sequential Least Square algorithm for high frequency GNSS Precise Point Positioning A. Barbu, J. Laurent-Varin, F. Perosanz, F. Mercier and J. Marty AVENUE project June, 20 A. Barbu, J. Laurent-Varin,
More informationVectors. Vector Practice Problems: Odd-numbered problems from
Vectors Vector Practice Problems: Odd-numbered problems from 3.1-3.21 After today, you should be able to: Understand vector notation Use basic trigonometry in order to find the x and y components of a
More informationTwin Telescopes at Onsala and Wettzell and their contribution to the VGOS System
Twin Telescopes at Onsala and Wettzell and their contribution to the VGOS System C. Schönberger, P. Gnilsen, J. Böhm, R. Haas Abstract During the last years the International VLBI Service for Geodesy and
More informationChapter 3 ALGEBRA. Overview. Algebra. 3.1 Linear Equations and Applications 3.2 More Linear Equations 3.3 Equations with Exponents. Section 3.
4 Chapter 3 ALGEBRA Overview Algebra 3.1 Linear Equations and Applications 3.2 More Linear Equations 3.3 Equations with Exponents 5 LinearEquations 3+ what = 7? If you have come through arithmetic, the
More informationBusiness Statistics. Lecture 9: Simple Regression
Business Statistics Lecture 9: Simple Regression 1 On to Model Building! Up to now, class was about descriptive and inferential statistics Numerical and graphical summaries of data Confidence intervals
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationMultiple Linear Regression for the Supervisor Data
for the Supervisor Data Rating 40 50 60 70 80 90 40 50 60 70 50 60 70 80 90 40 60 80 40 60 80 Complaints Privileges 30 50 70 40 60 Learn Raises 50 70 50 70 90 Critical 40 50 60 70 80 30 40 50 60 70 80
More informationLet X and Y denote two random variables. The joint distribution of these random
EE385 Class Notes 9/7/0 John Stensby Chapter 3: Multiple Random Variables Let X and Y denote two random variables. The joint distribution of these random variables is defined as F XY(x,y) = [X x,y y] P.
More informationGPS point positioning (using psuedorange from code). and earthquake location (using P or S travel times).
GPS point positioning (using psuedorange from code). and earthquake location (using P or S travel times). 15 17 From J. HOW, MIT, GPS-RJH5.pdf 18 P 1 = P 2 = P 3 = P 4 = ( X X 1 ) 2 + (Y Y 1 ) 2 + (Z
More informationImpact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames
Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames J. Böhm, H. Spicakova, L. Urquhart, P. Steigenberger, H. Schuh Abstract We compare the influence of two different a priori gradient
More informationDesign and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras
Design and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras Lecture No. # 10 Convergence Characteristics of Newton-Raphson Method
More informationSometimes the domains X and Z will be the same, so this might be written:
II. MULTIVARIATE CALCULUS The first lecture covered functions where a single input goes in, and a single output comes out. Most economic applications aren t so simple. In most cases, a number of variables
More informationLONG-TERM TRENDS IN THE AMOUNT OF ATMOSPHERIC WATER VAPOUR DERIVED FROM SPACE GEODETIC AND REMOTE SENSING TECHNIQUES
LONG-TERM TRENDS IN THE AMOUNT OF ATMOSPHERIC WATER VAPOUR DERIVED FROM SPACE GEODETIC AND REMOTE SENSING TECHNIQUES Rüdiger Haas, Tong Ning, and Gunnar Elgered Chalmers University of Technology, Onsala
More informationApproximation of ambiguity covariance matrix for integer de-correlation procedure in single-epoch GNSS positioning
he 9 th International Conference ENVIRONMENAL ENGINEERING 22 23 May 24, Vilnius, Lithuania SELECED PAPERS eissn 229-792 / eisbn 978-69-457-64-9 Available online at http://enviro.vgtu.lt Section: echnologies
More informationSolving Quadratic & Higher Degree Equations
Chapter 9 Solving Quadratic & Higher Degree Equations Sec 1. Zero Product Property Back in the third grade students were taught when they multiplied a number by zero, the product would be zero. In algebra,
More informationLinear Regression (continued)
Linear Regression (continued) Professor Ameet Talwalkar Professor Ameet Talwalkar CS260 Machine Learning Algorithms February 6, 2017 1 / 39 Outline 1 Administration 2 Review of last lecture 3 Linear regression
More informationThe Multivariate Normal Distribution 1
The Multivariate Normal Distribution 1 STA 302 Fall 2014 1 See last slide for copyright information. 1 / 37 Overview 1 Moment-generating Functions 2 Definition 3 Properties 4 χ 2 and t distributions 2
More informationPre-Algebra Notes Unit Three: Multi-Step Equations and Inequalities
Pre-Algebra Notes Unit Three: Multi-Step Equations and Inequalities A note to substitute teachers: pre-algebra teachers agree that all units of study are important, but understanding this unit seems to
More informationVelocity (Kalman Filter) Velocity (knots) Time (sec) Kalman Filter Accelerations 2.5. Acceleration (m/s2)
KALMAN FILERING 1 50 Velocity (Kalman Filter) 45 40 35 Velocity (nots) 30 25 20 15 10 5 0 0 10 20 30 40 50 60 70 80 90 100 ime (sec) 3 Kalman Filter Accelerations 2.5 Acceleration (m/s2) 2 1.5 1 0.5 0
More informationBindel, Spring 2016 Numerical Analysis (CS 4220) Notes for
Life beyond Newton Notes for 2016-04-08 Newton s method has many attractive properties, particularly when we combine it with a globalization strategy. Unfortunately, Newton steps are not cheap. At each
More information2010 GCE A Level H2 Maths Solution Paper 2 Section A: Pure Mathematics. 1i) x 2 6x + 34 = 0 6 ± x = 2
00 GCE A Level H Maths Solution Paper Section A: Pure Mathematics i) x 6x + 34 0 6 ± 36 36 x 6 ± 0i 3 ± 5i (ii) Since the coefficients are all real, another root of the equation is x i. [ x ( + i) ] [
More informationSolving Systems of Equations
Solving Systems of Equations Solving Systems of Equations What are systems of equations? Two or more equations that have the same variable(s) Solving Systems of Equations There are three ways to solve
More informationSYDE 112, LECTURE 7: Integration by Parts
SYDE 112, LECTURE 7: Integration by Parts 1 Integration By Parts Consider trying to take the integral of xe x dx. We could try to find a substitution but would quickly grow frustrated there is no substitution
More informationStat 206: Linear algebra
Stat 206: Linear algebra James Johndrow (adapted from Iain Johnstone s notes) 2016-11-02 Vectors We have already been working with vectors, but let s review a few more concepts. The inner product of two
More informationGaussian vectors and central limit theorem
Gaussian vectors and central limit theorem Samy Tindel Purdue University Probability Theory 2 - MA 539 Samy T. Gaussian vectors & CLT Probability Theory 1 / 86 Outline 1 Real Gaussian random variables
More informationA quadratic expression is a mathematical expression that can be written in the form 2
118 CHAPTER Algebra.6 FACTORING AND THE QUADRATIC EQUATION Textbook Reference Section 5. CLAST OBJECTIVES Factor a quadratic expression Find the roots of a quadratic equation A quadratic expression is
More informationNotes on Random Vectors and Multivariate Normal
MATH 590 Spring 06 Notes on Random Vectors and Multivariate Normal Properties of Random Vectors If X,, X n are random variables, then X = X,, X n ) is a random vector, with the cumulative distribution
More informationPhysics 403. Segev BenZvi. Propagation of Uncertainties. Department of Physics and Astronomy University of Rochester
Physics 403 Propagation of Uncertainties Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Maximum Likelihood and Minimum Least Squares Uncertainty Intervals
More informationEESC Geodesy with the Global Positioning System. Class 6: Point Positioning using Pseuduorange
EESC 9945 Geodesy with the Global Positioning System Class 6: Point Positioning using Pseuduorange GPS Positioning Solutions Point Positioning: Determination of the coordinates of a point with respect
More informationMath221: HW# 7 solutions
Math22: HW# 7 solutions Andy Royston November 7, 25.3.3 let x = e u. Then ln x = u, x2 = e 2u, and dx = e 2u du. Furthermore, when x =, u, and when x =, u =. Hence x 2 ln x) 3 dx = e 2u u 3 e u du) = e
More informationAtmospheric delay. X, Y, Z : satellite cartesian coordinates. Z : receiver cartesian coordinates. In the vacuum the signal speed c is constant
Atmospheric delay In the vacuum the signal speed c is constant c τ = ρ = ( X X ) + ( Y Y ) + ( Z Z ) S S S 2 S 2 S 2 X, Y, Z : receiver cartesian coordinates S S S X, Y, Z : satellite cartesian coordinates
More informationOriginally published as:
Originally published as: Seitz M., Steigenberger P., Artz T. (2xx) Consistent adjustment of combined terrestrial and celestial reference frames, Earth on the Edge: Science for a Sustainable Planet, IAG
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Positioning Sytems: Trilateration and Correlation
EECS 6A Designing Information Devices and Systems I Fall 08 Lecture Notes Note. Positioning Sytems: Trilateration and Correlation In this note, we ll introduce two concepts that are critical in our positioning
More informationMultiple Random Variables
Multiple Random Variables This Version: July 30, 2015 Multiple Random Variables 2 Now we consider models with more than one r.v. These are called multivariate models For instance: height and weight An
More informationLecture 13: Simple Linear Regression in Matrix Format. 1 Expectations and Variances with Vectors and Matrices
Lecture 3: Simple Linear Regression in Matrix Format To move beyond simple regression we need to use matrix algebra We ll start by re-expressing simple linear regression in matrix form Linear algebra is
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationLecture 22. r i+1 = b Ax i+1 = b A(x i + α i r i ) =(b Ax i ) α i Ar i = r i α i Ar i
8.409 An Algorithmist s oolkit December, 009 Lecturer: Jonathan Kelner Lecture Last time Last time, we reduced solving sparse systems of linear equations Ax = b where A is symmetric and positive definite
More informationNutation determination by means of GNSS
Nutation determination by means of GNSS - Comparison with VLBI Nicole Capitaine, Kunliang Yao SYRTE - Observatoire de Paris, CNRS/UPMC, France Introduction Space geodetic techniques cannot be used for
More informationReview for Final Exam, MATH , Fall 2010
Review for Final Exam, MATH 170-002, Fall 2010 The test will be on Wednesday December 15 in ILC 404 (usual class room), 8:00 a.m - 10:00 a.m. Please bring a non-graphing calculator for the test. No other
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationLinear Models Review
Linear Models Review Vectors in IR n will be written as ordered n-tuples which are understood to be column vectors, or n 1 matrices. A vector variable will be indicted with bold face, and the prime sign
More informationSolving Quadratic & Higher Degree Equations
Chapter 9 Solving Quadratic & Higher Degree Equations Sec 1. Zero Product Property Back in the third grade students were taught when they multiplied a number by zero, the product would be zero. In algebra,
More informationInteger-Valued Polynomials
Integer-Valued Polynomials LA Math Circle High School II Dillon Zhi October 11, 2015 1 Introduction Some polynomials take integer values p(x) for all integers x. The obvious examples are the ones where
More informationSection 20: Arrow Diagrams on the Integers
Section 0: Arrow Diagrams on the Integers Most of the material we have discussed so far concerns the idea and representations of functions. A function is a relationship between a set of inputs (the leave
More informationLecture 8 Analyzing the diffusion weighted signal. Room CSB 272 this week! Please install AFNI
Lecture 8 Analyzing the diffusion weighted signal Room CSB 272 this week! Please install AFNI http://afni.nimh.nih.gov/afni/ Next lecture, DTI For this lecture, think in terms of a single voxel We re still
More informationPhysics 411 Lecture 7. Tensors. Lecture 7. Physics 411 Classical Mechanics II
Physics 411 Lecture 7 Tensors Lecture 7 Physics 411 Classical Mechanics II September 12th 2007 In Electrodynamics, the implicit law governing the motion of particles is F α = m ẍ α. This is also true,
More informationLecture 5: September 12
10-725/36-725: Convex Optimization Fall 2015 Lecture 5: September 12 Lecturer: Lecturer: Ryan Tibshirani Scribes: Scribes: Barun Patra and Tyler Vuong Note: LaTeX template courtesy of UC Berkeley EECS
More informationImage Alignment and Mosaicing Feature Tracking and the Kalman Filter
Image Alignment and Mosaicing Feature Tracking and the Kalman Filter Image Alignment Applications Local alignment: Tracking Stereo Global alignment: Camera jitter elimination Image enhancement Panoramic
More informationOne sided tests. An example of a two sided alternative is what we ve been using for our two sample tests:
One sided tests So far all of our tests have been two sided. While this may be a bit easier to understand, this is often not the best way to do a hypothesis test. One simple thing that we can do to get
More informationA Probability Review
A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in
More informationECE 650 Lecture 4. Intro to Estimation Theory Random Vectors. ECE 650 D. Van Alphen 1
EE 650 Lecture 4 Intro to Estimation Theory Random Vectors EE 650 D. Van Alphen 1 Lecture Overview: Random Variables & Estimation Theory Functions of RV s (5.9) Introduction to Estimation Theory MMSE Estimation
More informationIntroduction to Gaussian Processes
Introduction to Gaussian Processes 1 Objectives to express prior knowledge/beliefs about model outputs using Gaussian process (GP) to sample functions from the probability measure defined by GP to build
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationPubl. Astron. Obs. Belgrade No. 91 (2012), REALIZATION OF ETRF2000 AS A NEW TERRESTRIAL REFERENCE FRAME IN REPUBLIC OF SERBIA
Publ. Astron. Obs. Belgrade No. 91 (2012), 199-206 Contributed paper REALIZATION OF ETRF2000 AS A NEW TERRESTRIAL REFERENCE FRAME IN REPUBLIC OF SERBIA D. BLAGOJEVIĆ and V. VASILIĆ Faculty of Civil Engineering,
More informationMachine Learning - MT & 5. Basis Expansion, Regularization, Validation
Machine Learning - MT 2016 4 & 5. Basis Expansion, Regularization, Validation Varun Kanade University of Oxford October 19 & 24, 2016 Outline Basis function expansion to capture non-linear relationships
More information1 Problem 1 Solution. 1.1 Common Mistakes. In order to show that the matrix L k is the inverse of the matrix M k, we need to show that
1 Problem 1 Solution In order to show that the matrix L k is the inverse of the matrix M k, we need to show that Since we need to show that Since L k M k = I (or M k L k = I). L k = I + m k e T k, M k
More information