ARTIFICIAL SATELLITES, Vol. 41, No DOI: /v
|
|
- Caren Jefferson
- 5 years ago
- Views:
Transcription
1 ARTIFICIAL SATELLITES, Vol. 41, No DOI: /v On InSAR AMBIGUITY RESOLUTION FOR DEFORMATION MONITORING P.J.G. Teunissen Department of Earth Observation and Space Sstems (DEOS) Delft Universit of Technolog Kluverweg 1, Delft, The Netherlands Fax: ; ABSTRACT Integer carrier phase ambiguit resolution is the ke to fast and highprecision satellite positioning and navigation. It applies to a great variet of current and future models of GPS, modernized GPS and Galileo. It also applies to stacked radar interferometr for deformation monitoring, see e.g. [Hanssen, et al, 001]. In this contribution we appl the integer least-squares principle to the rank defect model of stacked InSAR carrier phase data. We discuss two was of dealing with the rank defect for ambiguit resolution. One is based on the use of a priori data, the other is based on the use of an interval constraint on the deformation rate. 1 INTEGER LEAST-SQUARES (ILS) Consider the sstem of observation equations = Ax + Bz + e (1) where R m is the vector of observations, x R p and z Z n are the vectors of unknown parameters and e R m is the noise vector. Matrix (A, B) is given and assumed to be of full column rank. The structure of the sstem (1) is tpical for carrier phase GPS applications. In that case consists of the carrier phase and pseudo range data, z consists of the integer double difference ambiguities and x consists of the baseline components and possibl other real-valued parameters (e.g. atmospheric dela parameters). The model is however also applicable to radar interferometric phase ambiguit resolution problems, see e.g. [Bamler and Hartl, 1998], [Hanssen, 001], [Hanssen et al, 001], [Kampes and Hanssen, 004], [Kampes, 005]. In order to solve (1) in a least-squares sense, we follow the approach of [Teunissen, 1993]. Let Q be the variance matrix of and let. into a sum of three quadratic terms: e =(.) T (.). We first decompose e = Ax Bz = ê + ẑ z ẑ + ˆx(z) x ˆx z ()
2 18 in which ê = Aˆx Bẑ, ˆx and ẑ are the least-squares solutions for x and z, respectivel, assuming that z is real-valued instead of integer-valued, and ˆx(z) =(A T A) 1 A T ( Bz) (3) is the least-squares solution of x based on the assumption that z is known. In GPS terminolog, the solutions ˆx and ẑ are also referred to as the float solutions. Solution ˆx(z) is a conditional least-squares solution. Matrices Qẑ and Qˆx z are the variance matrices of ẑ and ˆx(z), respectivel. Note that the float solutions ˆx and ẑ would not be unique in case matrix (A, B) is rank defect. From () follows that min Ax + x R p,z Z Bz = ê n +min ẑ z +min ˆx(z) z Z n ẑ x R x p ˆx z This shows that the least-squares solutions for x R p and z Z n are given as (4) ž =argmin ẑ z Z z n ẑ and ˇx =ˆx(ž) (5) In [Teunissen, 1995] it is shown how ž can be computed in an efficient manner b means of the LAMBDA method. In [Teunissen, 1999] it is shown that the integer least-squares estimator ž is optimal in the sense that it maximizes the probabilit of correct integer estimation. For a brief review of statistical carrier phase ambiguit resolution, see e.g. [Teunissen, 001]. The probabilit of correct integer estimation (the so-called success-rate) is driven b the variance matrix of ẑ. The ADOP (Ambiguit Dilution of Precision), defined as ADOP = Qẑ 1 n (ccle) (6) is a simple scalar measure for the precision of ẑ. It is the geometric mean of the conditional standard deviations of ẑ. It has the advantage of being invariant for the class of admissible integer transformations. The ADOP can be used to get a first quick impression on the model s abilit to achieve successful ambiguit resolution. Fo this to be the case, the ADOP needs to be small enough (e.g. well below the ccle level). InSAR ILS WITH DEFORMATION DATA One of the simplest models of phase ambiguit resolution for stacked radar interferometr takes the form (see e.g. [Hanssen et al., 001]), i = a i x + z i + e i,i=1,...,n (7) with i the observed phase difference between two permanent scatterers (expressed in ccles), x the unknown deformation rate, z i the unknown integer ambiguit, e i the noise term, and a i = λ Δt i,whereλ denotes the wavelength (e.g. 5.8cm) andδt i is a time interval (expressed in ears) between the current time and a reference time. The model can be written in vector form as = ax + z + e (8)
3 This model appears to be a special case of (1), since A a, B I, p =1andm = n. Note however that the condition of having a coefficient matrix of full column rank is not met anmore, since matrix (a, I n ) has a rank defect of one. Hence, the sstem contains insufficient information to compute a unique float solution for x and z. This situation can be remedied b including data x o on the deformation rate x. This could be actual data (from earlier studies) or pseudo-data. As a result we obtain the full rank model [ ] = x o [ a 1 ] x + [ I 0 ] [ ] e z + e x with variance matrix [ Q 0 0 σ x in which σx is used to express the a priori uncertaint of the deformation rate. Now that the model is of full rank, the method of the previous section can be applied again. This is the approach which has been followed in [Hanssen et al, 001]. To determine the ADOP, we need the variance matrix of ẑ. It is given as Qẑ = Q + σx aat, from which the ADOP follows as ADOP = Qẑ 1 n = Q 1 n (1 + σ x a T a) 1 n (10) With Q = σφi n, a i = Δt λ i and small σ φ,weget ADOP σ φ ( λ σ x σ φ Δt i ) 1 n ] (9) (ccle) (11) This simple result can be used to get a quick impression of whether the model has enough strength for successful ambiguit resolution. 3 InSAR ILS WITH DEFORMATION CONSTRAINT Instead of working with a priori data on the deformation rate x, one ma also consider using a constraint in the form of an interval α x β, withα and β given. In practice one usuall has a fair idea over which range one can expect the deformation to occur. In that case it is not difficult to choose α and β. Hence, the ILS-problem now becomes min ax x R,z Z z subject to α x β (1) n Note that we could have added the constraint also to the model based on the a priori deformation data. In the following, however, we will work with the model (1). The procedure for solving it can also be used for the model based on the a priori deformation data. Unconstrained minimization Let us first consider minimizing the least-squares objective function without the constraint. With the use of the orthogonal projectors P a = a(a T a) 1 a T and Pa = I n P a,wehave min x R,z Z n ax z = min x R,z Z n Pa ( z) + P a ( z) ax = min z Z n Pa ( z) +min x R P a ( z) ax = min z Z n Pa ( z) +min x R ˆx(z) x a T a (13) 19
4 0 with ˆx(z) =(a T a) 1 a T ( z). This shows that the solution is given as ž =argmin P z Z n a ( z) Q 1 and ˇx =ˆx(ž) (14) This result seems to impl that in order to find ž, we need to search over the complete space of integers Z n. Clearl, this is not practicall feasible. Fortunatel, we can confine the search to a subset of Z n. The relevant integer vectors Let S : R n Z n be the integer map of componentwise integer rounding. Then all R n which are mapped b S to z Z n form the set S z = R n S() =z (15) This set is the n-dimensional unit-cube centred at z. Note,sinceQ is assumed to be a diagonal matrix, that S has the propert min z Z n z = S().Wema therefore decompose the least-squares objective function also as min ax x R,z Z z n =min min ax x R z Z z =min ax S( n x R ax) (16) This shows that the relevant integer vectors are those which are the centres of the unitcubes through which the line ax passes. The form the set Thus instead of (14) we ma now write Ω z = z Z n S( ax) =z, x R (17) ž =argmin Pa ( z) z Ω z and ˇx =ˆx(ž) (18) However, since the range of x has not et been bounded, the set Ω z still contains an infinite number of integer vectors. The finite set Ω z We will now make use of the constraint α x β. This allows us to reduce the infinite set Ω z to the finite set Ω z = z Zn S( ax) =z, α x β (19) Hence, we are now onl considering a line segment of length a T a(β α) of the line ax. Note that the length of the line segment is an approximate upper bound for the number of integer vectors contained in Ω z. Finding the integer vectors z Ω z To find all integer vectors that satisf S( ax) =z for α x β, we first make the problem a bit more tractable. Let x = x α, z = S( aα), z = z z,and = aα z.thenz = S( ax) =S( ax )+z and thus z = S( ax ) with 0 x γ =(β α) and S 0 (0) Hence, x is nonnegative and all components of lie within the interval [ 1, + 1 ]. Note that z = S( ax )isequivalenttothen intervals z i 1 i a i x z i + 1, i =1,...,n.
5 To find the elements of Ω z, one can now proceed as follows. Consider for each i the graph of the line ɛ i (x) = i a ix as function of x. The line starts at the point (0, i ), where for each i, i [ 1, + 1 ]. The line ma be located in either the first or the fourth quadrant, depending on whether a i 0ora i 0. Starting with x = 0, the first candidate integer vector is z = 0. The next candidate integer vector follows from which of the n lines ɛ i (x) first crosses one of the levels ± 1 modulus 1. For instance, if it is the line ɛ j (x) which first crosses the level ± 1, then the second candidate integer vector is given as z =(...,0, 1, 0,...) T,withthe1inslotj. In this wa one can proceed collecting all candidate integer vectors z for the range 0 x γ. From them the elements of Ω z are obtained as z = z + S( aα). Computing the final solution Now that the finite set Ω z is given, we compute 1 ž =argmin P z Ω a ( z) z and ˇx =ˆx(ž) (1) On this result one last check has to be performed, namel whether the constraint α ˇx =ˆx(ž) β is satisfied or not. If it is, then (1) is the solution sought. If it is not, then the solution is given b one of the two endpoints of the line segment, x = α or x = β. The solution is x = α if aα S( aα) < aβ S( aβ). Note that in practice one ma question the choice of the constraint if one of the two endpoints turns out to be the solution. The purpose of the constraint in our present application is namel to be able to construct a finite integer search space Ω z, but it is not chosen so as to influence the final solution. Thus if one of the two endpoints turns out to be the solution, it is likel that the constraint has been chosen too strict. 4 A POSTERIORI PROBABILITY OF THE INTEGER SOLUTION Once the integer vector ž has been computed, one would like to know whether one can have an trust in the solution. Intuitivel, this boils down to a comparison of the values taken b Pa ( z) for z Ω z. One would then have more trust in the solution if Pa ( ž) is significantl the smallest value, than when P a ( ž) is comparable to one or more of the other values. As shown below this is what is measured b the conditional distribution of z. The following assumptions are made: x, z N(ax + z, Q ), with conditional probabilit densit function (pdf) f( x, z). Thus f( x, z) exp 1 ax z (where stands for proportional to). x is uniforml distributed over the interval [α, β]. Its pdf g(x) is therefore proportional to the indicator function δ α,β (x) of the interval, g(x) δ α,β (x). For the distribution of z, we assume a noninformative prior. x and z are assumed to be independent and the distribution of z is denoted as h(z). We therefore have for the joint distribution f(x,, z) = f( x, z)g(x)h(z), from which the a posteriori distribution of z follows as h(z ) = f( x, z)g(x)dx z Z n f( x, z)g(x)dx = exp 1 P a ( z) k(z) z Z n exp 1 P a ( z) k(z) ()
6 where k(z) = 1 β σ exp 1 ( ) x ˆx(z) dx with σ 1 = π α σ a T a (3) For sufficientl small σ, k(z) will be close to one for ˆx(z) [α, β] and close to zero outside this range. Under this assumption () ma be approximated as h(z ) exp 1 P a ( z) z Ω z exp 1 P a ( z) (4) Note that ž is the maximizer of the right hand side. The a posteriori probabilit of the integer solution (i.e. the probabilit conditioned on the data) is given as h(ž ). The closer this value is to one, the more trust one will have in the computed solution. 5 REFERENCES Bamler, R. and Hartl, P. (1998): Snthetic aperture radar interferometr. Problems, 14, R1-R54. Inverse Hanssen, R.F. (001): Radar Interferometr: Data Interpretation and Error Analsis. Kluwer Academic Publishers, Dordrecht. Hanssen, R.F., P.J.G. Teunissen and P. Joosten (001): Phase ambiguit resolution for stacked radar interferometric data. In: Proc. KIS001, International Smposium on Kinematic Sstems in Geodes, Geomatics and Navigation, Banff, Canada, pp Kampes, B.M. and R.F. Hanssen (004): Ambiguit resolution for permanent scatterer interferometr. IEEE Transactions on Geoscience and Remote Sensing, 4(11), Kampes, B.M. (005): Displacement Parameter Estimation using Permanent Scatterer Interferometr. Netherlands Geodetic Commission. Publications on Geodes, 168 p. Teunissen, P.J.G. (1993): Least-squares estimation of the integer GPS ambiguities. Invited Lecture, Section IV Theor and Methodolog, IAG General Meeting, Beijing, China, August Also in: LGR Series, No. 6, Delft Geodetic Computing Centre. Teunissen, P.J.G. (1995): The least-squares ambiguit decorrelation adjustment: a method for fast GPS integer ambiguit estimation. Journal of Geodes, 70: Teunissen, P.J.G. (1999): An optimalit propert of the integer least- squares estimator. Journal of Geodes, 73: Teunissen, P.J.G. (001): Statistical GNSS carrier phase ambiguit resolution: a review. IEEE, /01, pp Received: , Reviewed: , b W. Pachelski, Accepted:
New Global Navigation Satellite System Ambiguity Resolution Method Compared to Existing Approaches
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 29, No. 4, Jul August 2006 New Global Navigation Satellite Sstem Ambiguit Resolution Method Compared to Existing Approaches Sandra Verhagen and Peter J.
More informationIt is the purpose of this contribution to present the probability distribution of the 'xed' baseline. This will be done for the case the ambiguities a
THE DISTRIBUTION OF THE GPS BASELINE IN CASE OF INTEGER LEAST-SQUARES AMBIGUITY ESTIMATION P.J.G. Teunissen Delft Geodetic Computing Centre (LGR) Faculty of Geodesy Delft University of Technology Thijsseweg
More informationNoise Characteristics in High Precision GPS Positioning
Noise Characteristics in High Precision GPS Positioning A.R. Amiri-Simkooei, C.C.J.M. Tiberius, P.J.G. Teunissen, Delft Institute of Earth Observation and Space systems (DEOS), Delft University of Technology,
More informationIonosphere influence on success rate of GPS ambiguity resolution in a satellite formation flying
Journal of Physics: Conference Series PAPER OPEN ACCESS Ionosphere influence on success rate of GPS ambiguity resolution in a satellite formation flying To cite this article: Leandro Baroni 2015 J. Phys.:
More informationInfluence of ambiguity precision on the success rate of GNSS integer ambiguity bootstrapping
J Geod (007) 8:35 358 DOI 0007/s0090-006-0-3 ORIGINAL ARTICLE Influence of ambiguity precision on the success rate of GNSS integer ambiguity bootstrapping P J G Teunissen Received: March 006 / Accepted:
More informationLeast-Squares Estimation of. 1 Introduction
Abstract Least-Squares Estimation of the Integer GPS Ambiguities PJG Teunissen Delft Geodetic Computing Centre (LGR) Department of the Geodetic Engineering Delft University of Technology Thijsseweg, 69
More informationOn the Spectrum of the GPS DD-Ambiguities
On the Spectrum of the GPS DD-Ambiguities Dr. Peter Teunissen, Paul de Jonge, and Christiaan Tiberius, Delft University of Technology (in: Proceedings ION GPS-94, Salt Lake City, Utah, USA, September 20-23,
More informationARTIFICIAL SATELLITES, Vol. 41, No DOI: /v THE LAMBDA METHOD FOR THE GNSS COMPASS
ARTIFICIAL SATELLITES, Vol. 41, No. 3-2006 DOI: 10.2478/v10018-007-0009-1 THE LAMBDA METHOD FOR THE GNSS COMPASS P.J.G. Teunissen Delft Institute for Earth Observation and Space Systems (DEOS Delft University
More informationDeformation measurement using SAR interferometry: quantitative aspects
Deformation measurement using SAR interferometry: quantitative aspects Michele Crosetto (1), Erlinda Biescas (1), Ismael Fernández (1), Ivan Torrobella (1), Bruno Crippa (2) (1) (2) Institute of Geomatics,
More informationInteger least-squares theory for the GNSS compass
J Geod 2010 84:433 447 DOI 10.1007/s00190-010-0380-8 ORIGINAL ARTICLE Integer least-squares theory for the GNSS compass P. J. G. Teunissen Received: 31 December 2009 / Accepted: 11 March 2010 / Published
More informationCarrier-phase Ambiguity Success Rates for Integrated GPS-Galileo Satellite Navigation
Proceedings Space, Aeronautical and Navigational Electronics Symposium SANE2, The Institute of Electronics, Information and Communication Engineers (IEICE), Japan, Vol., No. 2, pp. 3- Carrier-phase Ambiguity
More informationCarrier Phase Integer Ambiguity Resolution Recent results and open issues
Carrier Phase Integer Ambiguity Resolution Recent results and open issues Sandra Verhagen DEOS, Delft University of Technology, The Netherlands Kolloquium Satellitennavigation, TU München, 9 February 2010
More informationAutocorrelation Functions in GPS Data Processing: Modeling Aspects
Autocorrelation Functions in GPS Data Processing: Modeling Aspects Kai Borre, Aalborg University Gilbert Strang, Massachusetts Institute of Technology Consider a process that is actually random walk but
More informationInteger carrier-phase ambiguity resolution is the key to fast. Integer Aperture Estimation. A Framework for GNSS Ambiguity Acceptance Testing
working papers Integer Aperture Estimation A Framework for GNSS Ambiguity Acceptance Testing Tracking the carrier phase of GNSS signals has evolved into a widespread practice for achieving rapid and very
More informationLeast-Squares Variance Component Estimation:
Least-Squares Variance Component Estimation: Theory and GPS Applications A.R. Amiri-Simkooei Least-Squares Variance Component Estimation: Theory and GPS Applications PROEFSCHRIFT ter verkrijging van de
More informationResearch Article Chaotic Attractor Generation via a Simple Linear Time-Varying System
Discrete Dnamics in Nature and Societ Volume, Article ID 836, 8 pages doi:.//836 Research Article Chaotic Attractor Generation via a Simple Linear Time-Varing Sstem Baiu Ou and Desheng Liu Department of
More informationOn the realistic stochastic model of GPS observables: Implementation and Performance
he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov
More informationA. Derivation of regularized ERM duality
A. Derivation of regularized ERM dualit For completeness, in this section we derive the dual 5 to the problem of computing proximal operator for the ERM objective 3. We can rewrite the primal problem as
More informationSubband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014
Subband Coding and Wavelets National Chiao Tung Universit Chun-Jen Tsai /4/4 Concept of Subband Coding In transform coding, we use N (or N N) samples as the data transform unit Transform coefficients are
More informationMODELING INTERFEROGRAM STACKS FOR SENTINEL - 1
MODELING INTERFEROGRAM STACKS FOR SENTINEL - 1 Fabio Rocca (1) (1) Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italy, Email: rocca@elet.polimi.it ABSTRACT The dispersion of the optimal estimate
More informationLecture 4 Propagation of errors
Introduction Lecture 4 Propagation of errors Example: we measure the current (I and resistance (R of a resistor. Ohm's law: V = IR If we know the uncertainties (e.g. standard deviations in I and R, what
More informationNoise covariance model for time-series InSAR analysis. Piyush Agram, Mark Simons
Noise covariance model for time-series InSAR analysis Piyush Agram, Mark Simons Motivation Noise covariance in individual interferograms has been well studied. (Hanssen, 2001) Extend the noise model to
More informationUSING THE INTEGER DECORRELATION PROCEDURE TO INCREASE OF THE EFFICIENCY OF THE MAFA METHOD
ARIFICIAL SAELLIES, Vol. 46, No. 3 2 DOI:.2478/v8-2-2- USING HE INEGER DECORRELAION PROCEDURE O INCREASE OF HE EFFICIENCY OF HE MAFA MEHOD S. Cellmer Institute of Geodesy University of Warmia and Mazury
More informationExperimental Uncertainty Review. Abstract. References. Measurement Uncertainties and Uncertainty Propagation
Experimental Uncertaint Review Abstract This is intended as a brief summar of the basic elements of uncertaint analsis, and a hand reference for laborator use. It provides some elementar "rules-of-thumb"
More informationEntropy, Relative Entropy and Mutual Information Exercises
Entrop, Relative Entrop and Mutual Information Exercises Exercise 2.: Coin Flips. A fair coin is flipped until the first head occurs. Let X denote the number of flips required. (a) Find the entrop H(X)
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 information3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES
3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES 3.1 Introduction In this chapter we will review the concepts of probabilit, rom variables rom processes. We begin b reviewing some of the definitions
More informationODIF FOR MORPHOLOGICAL FILTERS. Sari Peltonen and Pauli Kuosmanen
ODIF FOR MORPHOLOGICAL FILTERS Sari Peltonen Pauli Kuosmanen Signal Processing Laborator, Digital Media Institute, Tampere Universit of Technolog P.O. Box 553, FIN-331 TAMPERE, FINLAND e-mail: sari@cs.tut.fi
More informationUsing Modern Satellite Data to Constrain Climate Models, and Bayesian Climate Projection
Using Modern Satellite Data to Constrain Climate Models, and Baesian Climate Projection Stephen Lero, Yi Huang, Richard Good, James Anderson (Harvard Universit) International Detection and Attribution
More informationEvaluation of Two Integer Ambiguity Resolution Methods for Real Time GPS Positioning
Evaluation of Two Integer Ambiguity Resolution Methods for Real Time GPS Positioning LEANDRO BARONI, HÉLIO KOITI KUGA Space Mechanics and Control Division National Institute for Space Research 12227-010,
More informationCovariance and Correlation Class 7, Jeremy Orloff and Jonathan Bloom
1 Learning Goals Covariance and Correlation Class 7, 18.05 Jerem Orloff and Jonathan Bloom 1. Understand the meaning of covariance and correlation. 2. Be able to compute the covariance and correlation
More informationy = f(x + 4) a) Example: A repeating X by using two linear equations y = ±x. b) Example: y = f(x - 3). The translation is
Answers Chapter Function Transformations. Horizontal and Vertical Translations, pages to. a h, k h, k - c h -, k d h 7, k - e h -, k. a A (-,, B (-,, C (-,, D (,, E (, A (-, -, B (-,, C (,, D (, -, E (,
More informationCarrier Phase Techniques
Carrier Phase Techniques Michael P Vitus August 30, 2005 Introduction This paper provides an explanation of the methods as well as the wor performed on the GPS project under Professor Ng ding the summer
More informationERAD Water vapor observations with SAR, microwave radiometer and GPS: comparison of scaling characteristics
Proceedings of ERAD (2002): 190 194 c Copernicus GmbH 2002 ERAD 2002 Water vapor observations with SAR, microwave radiometer and GPS: comparison of scaling characteristics D. N. Moisseev 1, R. F. Hanssen
More informationMAT 1275: Introduction to Mathematical Analysis. Graphs and Simplest Equations for Basic Trigonometric Functions. y=sin( x) Function
MAT 275: Introduction to Mathematical Analsis Dr. A. Rozenblum Graphs and Simplest Equations for Basic Trigonometric Functions We consider here three basic functions: sine, cosine and tangent. For them,
More informationUniversity of California, Los Angeles Department of Statistics. Joint probability distributions
Universit of California, Los Angeles Department of Statistics Statistics 100A Instructor: Nicolas Christou Joint probabilit distributions So far we have considered onl distributions with one random variable.
More informationELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications
ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications Professor M. Chiang Electrical Engineering Department, Princeton University March
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 information1 Cricket chirps: an example
Notes for 2016-09-26 1 Cricket chirps: an example Did you know that you can estimate the temperature by listening to the rate of chirps? The data set in Table 1 1. represents measurements of the number
More informationBACKWARD FOKKER-PLANCK EQUATION FOR DETERMINATION OF MODEL PREDICTABILITY WITH UNCERTAIN INITIAL ERRORS
BACKWARD FOKKER-PLANCK EQUATION FOR DETERMINATION OF MODEL PREDICTABILITY WITH UNCERTAIN INITIAL ERRORS. INTRODUCTION It is widel recognized that uncertaint in atmospheric and oceanic models can be traced
More informationSolutions Week 2. D-ARCH Mathematics Fall Linear dependence & matrix multiplication. 1. A bit of vector algebra.
D-ARCH Mathematics Fall 4 Solutions Week Linear dependence & matrix multiplication A bit of vector algebra (a) Find all vectors perpendicular to Drawing a sketch might help! Solution : All vectors in the
More informationCHAPTER 5. Jointly Probability Mass Function for Two Discrete Distributed Random Variables:
CHAPTER 5 Jointl Distributed Random Variable There are some situations that experiment contains more than one variable and researcher interested in to stud joint behavior of several variables at the same
More informationStability Analysis for Linear Systems under State Constraints
Stabilit Analsis for Linear Sstems under State Constraints Haijun Fang Abstract This paper revisits the problem of stabilit analsis for linear sstems under state constraints New and less conservative sufficient
More informationSubfield value sets over finite fields
Subfield value sets over finite fields Wun-Seng Chou Institute of Mathematics Academia Sinica, and/or Department of Mathematical Sciences National Chengchi Universit Taipei, Taiwan, ROC Email: macws@mathsinicaedutw
More informationBest prediction in linear models with mixed integer/real unknowns: theory and application
J Geod (2007) 81:759 780 DOI 10.1007/s00190-007-0140-6 ORIGINAL ARTICLE Best prediction in linear models with mixed integer/real unknowns: theory and application P. J. G. Teunissen Received: 19 September
More information1 GSW Gaussian Elimination
Gaussian elimination is probabl the simplest technique for solving a set of simultaneous linear equations, such as: = A x + A x + A x +... + A x,,,, n n = A x + A x + A x +... + A x,,,, n n... m = Am,x
More informationMultiscale Principal Components Analysis for Image Local Orientation Estimation
Multiscale Principal Components Analsis for Image Local Orientation Estimation XiaoGuang Feng Peman Milanfar Computer Engineering Electrical Engineering Universit of California, Santa Cruz Universit of
More information2 Random Variable Generation
2 Random Variable Generation Most Monte Carlo computations require, as a starting point, a sequence of i.i.d. random variables with given marginal distribution. We describe here some of the basic methods
More informationSupplementary Material: A Robust Approach to Sequential Information Theoretic Planning
Supplementar aterial: A Robust Approach to Sequential Information Theoretic Planning Sue Zheng Jason Pacheco John W. Fisher, III. Proofs of Estimator Properties.. Proof of Prop. Here we show how the bias
More informationInversion of Earthquake Rupture Process:Theory and Applications
Inversion of Earthquake Rupture Process:Theory and Applications Yun-tai CHEN 12 * Yong ZHANG 12 Li-sheng XU 2 1School of the Earth and Space Sciences, Peking University, Beijing 100871 2Institute of Geophysics,
More information0.24 adults 2. (c) Prove that, regardless of the possible values of and, the covariance between X and Y is equal to zero. Show all work.
1 A socioeconomic stud analzes two discrete random variables in a certain population of households = number of adult residents and = number of child residents It is found that their joint probabilit mass
More informationCOMPRESSIVE (CS) [1] is an emerging framework,
1 An Arithmetic Coding Scheme for Blocked-based Compressive Sensing of Images Min Gao arxiv:1604.06983v1 [cs.it] Apr 2016 Abstract Differential pulse-code modulation (DPCM) is recentl coupled with uniform
More informationLecture 7. Econ August 18
Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar
More informationPICONE S IDENTITY FOR A SYSTEM OF FIRST-ORDER NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS
Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 143, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu PICONE S IDENTITY
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More informationConstrained Optimization and Lagrangian Duality
CIS 520: Machine Learning Oct 02, 2017 Constrained Optimization and Lagrangian Duality Lecturer: Shivani Agarwal Disclaimer: These notes are designed to be a supplement to the lecture. They may or may
More informationOn the Exploitation of Target Statistics for SAR Interferometry Applications
On the Exploitation of Target Statistics for SAR Interferometry Applications A. Monti Guarnieri, S. Tebaldini Politecnico di Milano Abstract This paper focuses on multi-image Synthetic Aperture Radar Interferometry
More information2938 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 46, NO. 11, NOVEMBER Integer Parameter Estimation in Linear Models with Applications to GPS
2938 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 46, NO. 11, NOVEMBER 1998 Integer Parameter Estimation in Linear Models with Applications to GPS Arash Hassibi Stephen Boyd Abstract We consider parameter
More informationFrom Bayes to Extended Kalman Filter
From Bayes to Extended Kalman Filter Michal Reinštein Czech Technical University in Prague Faculty of Electrical Engineering, Department of Cybernetics Center for Machine Perception http://cmp.felk.cvut.cz/
More informationModeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation
Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung
More informationQuantum Computer Algorithms: basics
Quantum Computer Algorithms: basics Michele Mosca Canada Research Chair in Quantum Computation SQUINT Summer School on Quantum Information Processing June 3 Overview A quantum computing model Basis changes
More informationLecture: Quadratic optimization
Lecture: Quadratic optimization 1. Positive definite och semidefinite matrices 2. LDL T factorization 3. Quadratic optimization without constraints 4. Quadratic optimization with constraints 5. Least-squares
More informationLet X denote a random variable, and z = h(x) a function of x. Consider the
EE385 Class Notes 11/13/01 John Stensb Chapter 5 Moments and Conditional Statistics Let denote a random variable, and z = h(x) a function of x. Consider the transformation Z = h(). We saw that we could
More informationMechanical Spring Reliability Assessments Based on FEA Generated Fatigue Stresses and Monte Carlo Simulated Stress/Strength Distributions
Purdue Universit Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 4 Mechanical Spring Reliabilit Assessments Based on FEA Generated Fatigue Stresses and Monte
More information2: Distributions of Several Variables, Error Propagation
: Distributions of Several Variables, Error Propagation Distribution of several variables. variables The joint probabilit distribution function of two variables and can be genericall written f(, with the
More informationMetrol. Meas. Syst., Vol. XVIII (2011), No. 2, pp METROLOGY AND MEASUREMENT SYSTEMS. Index , ISSN
Metrol. Meas. Sst., Vol. XVIII (11), No., pp. 335 34 METROLOGY AND MEASUREMENT SYSTEMS Index 3393, ISSN 86-89 www.metrolog.pg.gda.pl INVESTIGATION OF THE STATISTICAL METHOD OF TIME DELAY ESTIMATION BASED
More informationChapter 7. Extremal Problems. 7.1 Extrema and Local Extrema
Chapter 7 Extremal Problems No matter in theoretical context or in applications many problems can be formulated as problems of finding the maximum or minimum of a function. Whenever this is the case, advanced
More informationTHE UNIVERSITY OF HONG KONG DEPARTMENT OF MATHEMATICS
THE UNIVERSITY OF HONG KONG DEPARTMENT OF MATHEMATICS MATH853: Linear Algebra, Probability and Statistics May 5, 05 9:30a.m. :30p.m. Only approved calculators as announced by the Examinations Secretary
More informationStability Analysis of a Geometrically Imperfect Structure using a Random Field Model
Stabilit Analsis of a Geometricall Imperfect Structure using a Random Field Model JAN VALEŠ, ZDENĚK KALA Department of Structural Mechanics Brno Universit of Technolog, Facult of Civil Engineering Veveří
More informationBayesian decision theory Introduction to Pattern Recognition. Lectures 4 and 5: Bayesian decision theory
Bayesian decision theory 8001652 Introduction to Pattern Recognition. Lectures 4 and 5: Bayesian decision theory Jussi Tohka jussi.tohka@tut.fi Institute of Signal Processing Tampere University of Technology
More informationVISION TRACKING PREDICTION
VISION RACKING PREDICION Eri Cuevas,2, Daniel Zaldivar,2, and Raul Rojas Institut für Informati, Freie Universität Berlin, ausstr 9, D-495 Berlin, German el 0049-30-83852485 2 División de Electrónica Computación,
More informationInterpreting Regression Results
Interpreting Regression Results Carlo Favero Favero () Interpreting Regression Results 1 / 42 Interpreting Regression Results Interpreting regression results is not a simple exercise. We propose to split
More informationMATH 2070 Test 3 (Sections , , & )
Multiple Choice: Use a # pencil and completel fill in each bubble on our scantron to indicate the answer to each question. Each question has one correct answer. If ou indicate more than one answer, or
More informationx. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).
.8.6 µ =, σ = 1 µ = 1, σ = 1 / µ =, σ =.. 3 1 1 3 x Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ ). The Gaussian distribution Probably the most-important distribution in all of statistics
More informationConstructing high-resolution, absolute maps of atmospheric water vapor by combining InSAR and GNSS observations
Constructing high-resolution, absolute maps of atmospheric water vapor by combining InSAR and GNSS observations Fadwa Alshawaf, Stefan Hinz, Michael Mayer, Franz J. Meyer fadwa.alshawaf@kit.edu INSTITUTE
More informationMILES: MATLAB package for solving Mixed Integer LEast Squares problems Theory and Algorithms
MILES: MATLAB package for solving Mixed Integer LEast Squares problems Theory and Algorithms Xiao-Wen Chang and Tianyang Zhou Scientific Computing Laboratory School of Computer Science McGill University
More informationLecture 7. Floating point arithmetic and stability
Lecture 7 Floating point arithmetic and stability 2.5 Machine representation of numbers Scientific notation: 23 }{{} }{{} } 3.14159265 {{} }{{} 10 sign mantissa base exponent (significand) s m β e A floating
More informationLecture 8: Channel Capacity, Continuous Random Variables
EE376A/STATS376A Information Theory Lecture 8-02/0/208 Lecture 8: Channel Capacity, Continuous Random Variables Lecturer: Tsachy Weissman Scribe: Augustine Chemparathy, Adithya Ganesh, Philip Hwang Channel
More information15. Conic optimization
L. Vandenberghe EE236C (Spring 216) 15. Conic optimization conic linear program examples modeling duality 15-1 Generalized (conic) inequalities Conic inequality: a constraint x K where K is a convex cone
More informationELE539A: Optimization of Communication Systems Lecture 6: Quadratic Programming, Geometric Programming, and Applications
ELE539A: Optimization of Communication Systems Lecture 6: Quadratic Programming, Geometric Programming, and Applications Professor M. Chiang Electrical Engineering Department, Princeton University February
More informationSupplementary information. Analytical Techniques and Measurement Uncertainty
Supplementary information Analytical Techniques and Measurement Uncertainty Interferomeric SAR (InSAR) allows measurement of millimetre-level surface displacements, including land subsidence, in the radar
More informationOrdinal One-Switch Utility Functions
Ordinal One-Switch Utilit Functions Ali E. Abbas Universit of Southern California, Los Angeles, California 90089, aliabbas@usc.edu David E. Bell Harvard Business School, Boston, Massachusetts 0163, dbell@hbs.edu
More informationAtmospheric Phase Screen (APS) estimation and modeling for radar interferometry
Atmospheric Phase Screen (APS) estimation and modeling for radar interferometry Ramon Hanssen¹, Alessandro Ferretti², Marco Bianchi², Rossen Grebenitcharsky¹, Frank Kleijer¹, Ayman Elawar¹ ¹ Delft University
More informationLinear programming: Theory
Division of the Humanities and Social Sciences Ec 181 KC Border Convex Analsis and Economic Theor Winter 2018 Topic 28: Linear programming: Theor 28.1 The saddlepoint theorem for linear programming The
More informationEEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as
L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationDynamic data processing recursive least-squares
Dynamic data processing recursive least-squares Dynamic data processing recursive least-squares Delft University of Technology Department of Mathematical Geodesy and Positioning VSSD Series on Mathematical
More informationestimator is a random vector. Conceptually we may write a = F (y), where a denotes the integer ambiguity vector, y the vector of observables, and F (:
ON THE INTEGER NORMAL DISTRIBUTION OF THE GPS AMBIGUITIES P.J.G. Teunissen Delft Geodetic Computing Centre (LGR) Faculty of Geodesy Delft University of Technology Thijsseweg 2629 JA Delft, The Netherlands
More information6.451 Principles of Digital Communication II Wednesday, May 4, 2005 MIT, Spring 2005 Handout #22. Problem Set 9 Solutions
6.45 Principles of Digital Communication II Wednesda, Ma 4, 25 MIT, Spring 25 Hand #22 Problem Set 9 Solutions Problem 8.3 (revised) (BCJR (sum-product) decoding of SPC codes) As shown in Problem 6.4 or
More information7-6. nth Roots. Vocabulary. Geometric Sequences in Music. Lesson. Mental Math
Lesson 7-6 nth Roots Vocabular cube root n th root BIG IDEA If is the nth power of, then is an nth root of. Real numbers ma have 0, 1, or 2 real nth roots. Geometric Sequences in Music A piano tuner adjusts
More informationIntegration. Tuesday, December 3, 13
4 Integration 4.3 Riemann Sums and Definite Integrals Objectives n Understand the definition of a Riemann sum. n Evaluate a definite integral using properties of definite integrals. 3 Riemann Sums 4 Riemann
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace
More informationreview math0410 (1-174) and math 0320 ( ) aafinm mg
Eam Name review math04 (1-174) and math 0320 (17-243) 03201700aafinm0424300 mg MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Simplif. 1) 7 2-3 A)
More informationSCORE BOOSTER JAMB PREPARATION SERIES II
BOOST YOUR JAMB SCORE WITH PAST Polynomials QUESTIONS Part II ALGEBRA by H. O. Aliu J. K. Adewole, PhD (Editor) 1) If 9x 2 + 6xy + 4y 2 is a factor of 27x 3 8y 3, find the other factor. (UTME 2014) 3x
More informationPrimer on Mathematical Geodesy CTB3310 / CTB3425
Delft Universit of Technolog Primer on Mathematical Geodes CTB3310 / CTB3425 Christian Tiberius ii The front cover shows a two dimensional normal Probabilit Densit Function. C.C.J.M. Tiberius Lecture notes
More informationAnn-Brith Strömberg. Lecture 4 Linear and Integer Optimization with Applications 1/10
MVE165/MMG631 Linear and Integer Optimization with Applications Lecture 4 Linear programming: degeneracy; unbounded solution; infeasibility; starting solutions Ann-Brith Strömberg 2017 03 28 Lecture 4
More informationComputation of Csiszár s Mutual Information of Order α
Computation of Csiszár s Mutual Information of Order Damianos Karakos, Sanjeev Khudanpur and Care E. Priebe Department of Electrical and Computer Engineering and Center for Language and Speech Processing
More informationThe Multivariate Gaussian Distribution
The Multivariate Gaussian Distribution Chuong B. Do October, 8 A vector-valued random variable X = T X X n is said to have a multivariate normal or Gaussian) distribution with mean µ R n and covariance
More informationConvex envelopes, cardinality constrained optimization and LASSO. An application in supervised learning: support vector machines (SVMs)
ORF 523 Lecture 8 Princeton University Instructor: A.A. Ahmadi Scribe: G. Hall Any typos should be emailed to a a a@princeton.edu. 1 Outline Convexity-preserving operations Convex envelopes, cardinality
More informationGrothendieck s Inequality
Grothendieck s Inequality Leqi Zhu 1 Introduction Let A = (A ij ) R m n be an m n matrix. Then A defines a linear operator between normed spaces (R m, p ) and (R n, q ), for 1 p, q. The (p q)-norm of A
More information