Créer une présentation
Télécharger la présentation

Download

Download Presentation

Earth Science Applications of Space Based Geodesy DES-7355 Tu-Th 9:40-11:05 Seminar Room in 3892 Central

110 Vues
Download Presentation

Télécharger la présentation
## Earth Science Applications of Space Based Geodesy DES-7355 Tu-Th 9:40-11:05 Seminar Room in 3892 Central

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -

**Earth Science Applications of Space Based Geodesy**DES-7355 Tu-Th 9:40-11:05 Seminar Room in 3892 Central Ave. (Long building) Bob Smalley Office: 3892 Central Ave, Room 103 678-4929 Office Hours – Wed 14:00-16:00 or if I’m in my office. http://www.ceri.memphis.edu/people/smalley/ESCI7355/ESCI_7355_Applications_of_Space_Based_Geodesy.html Class 7**More inversion pitfalls**Bill and Ted's misadventure. Bill and Ted are geo-chemists who wish to measure the number of grams of each of three different minerals A,B,C held in a single rock sample. Let a be the number of grams of A, b be the number of grams of B, c be the number of grams of C d be the number of grams in the sample. From Todd Will**By performing complicated experiments Bill and Ted are able**to measure four relationships between a,b,c,d which they record in the matrix below: Now we have more equations than we need What to do? From Todd Will**One thing to do is throw out one of the equations**(in reality only a Mathematician is naïve enough to think that three equations is sufficient to solve for three unknowns – but lets try it anyway). So throw out one - leaving (different A and b from before) From Todd Will**Remembering some of their linear algebra they know that the**matrix is not invertible if the determinant is zero, so they check that OK so far (or “fat, dumb and happy”) From Todd Will**So now we can compute**So now we’re done. From Todd Will**Or are we?**From Todd Will**Next they realize that the measurements are really only good**to 0.1 So they round to 0.1 and do it again From Todd Will**Now they notice a small problem –**They get a very different answer (and they don’t notice they have a bigger physical problem in that they have negative weights/amounts!) From Todd Will**So what’s the problem?**First find the SVD of A. Since there are three non-zero values on the diagonal A is invertible From Todd Will**BUT, one of the singular values is much, much less than the**others So the matrix is “almost” rank 2 (which would be non-invertible) From Todd Will**We can also calculate the SVD of A-1**From Todd Will**So now we can see what happened**(why the two answers were so different) Let y be the first version of b Let y’ be the second version of b (to 0.1) So A-1 stretches vectors parallel to h3 and a3 by a factor of 5000. From Todd Will**We have 4 unknowns (xR,yR,zR and tR)**And 4 (nonlinear) equations (later we will allow more satellites) So we can solve for the unknowns**Again, we cannot solve this directly**Will solve interatively by 1) Assuming a location 2) Linearizing the range equations 3) Use least squares to compute new (better) location 4) Go back to 1 using location from 3 We do this till some convergence criteria is met (if we’re lucky) Blewitt, Basics of GPS in “Geodetic Applications of GPS”**linearize**So - for one satellite we have Blewitt, Basics of GPS in “Geodetic Applications of GPS”**Linearize (first two terms of Taylor Series)**Blewitt, Basics of GPS in “Geodetic Applications of GPS”**Residual**Difference between observed and calculated (linearized) Blewitt, Basics of GPS in “Geodetic Applications of GPS”**So we have the following for one satellite**Which we can recast in matrix form Blewitt, Basics of GPS in “Geodetic Applications of GPS”**For m satellites (where m≥4)**Which is usually written as Blewitt, Basics of GPS in “Geodetic Applications of GPS”**Calculate the derivatives**Blewitt, Basics of GPS in “Geodetic Applications of GPS”**So we get**Is function of direction to satellite Note last column is a constant Blewitt, Basics of GPS in “Geodetic Applications of GPS”**Consider some candidate solution x’**Then we can write b are the observations are the residuals We would like to find the x’ that minimizes Blewitt, Basics of GPS in “Geodetic Applications of GPS”**So the question now is how to find this x’**One way, and the way we will do it, Least Squares Blewitt, Basics of GPS in “Geodetic Applications of GPS”**Since we have already done this – we’ll go fast**Use solution to linearized form of observation equations to write estimated residuals Vary value of x to minimize Blewitt, Basics of GPS in “Geodetic Applications of GPS”**Normal equations**Solution to normal equations**Assumes**Inverse exists (m greater than or equal to 4, necessary but not sufficient condition) Can have problems similar to earthquake locating (two satellites in “same” direction for example – has effect of reducing rank by one)**GPS tutorial Signals and Data**http://www.unav-micro.com/about_gps.htm**GPS tutorial Signals and Data**http://www.unav-micro.com/about_gps.htm**Elementary Concepts**Variables: things that we measure, control, or manipulate in research. They differ in many respects, most notably in the role they are given in our research and in the type of measures that can be applied to them. From G. Mattioli**Observational vs. experimental research.**Most empirical research belongs clearly to one of those two general categories. In observational research we do not (or at least try not to) influence any variables but only measure them and look for relations (correlations) between some set of variables. In experimental research, we manipulate some variables and then measure the effects of this manipulation on other variables. From G. Mattioli**Observational vs. experimental research.**Dependent vs. independent variables. Independent variables are those that are manipulated whereas dependent variables are only measured or registered. From G. Mattioli**Variable Types and Information Content**Measurement scales. Variables differ in "how well" they can be measured. Measurement error involved in every measurement, which determines the "amount of information” obtained. Another factor is the variable’s "type of measurement scale." From G. Mattioli**Variable Types and Information Content**Nominal variables allow for only qualitative classification. That is, they can be measured only in terms of whether the individual items belong to some distinctively different categories, but we cannot quantify or even rank order those categories. Typical examples of nominal variables are gender, race, color, city, etc. From G. Mattioli**Variable Types and Information Content**Ordinal variables allow us to rank order the items we measure in terms of which has less and which has more of the quality represented by the variable, but still they do not allow us to say "how much more.” A typical example of an ordinal variable is the socioeconomic status of families. From G. Mattioli**Variable Types and Information Content**Interval variables allow us not only to rank order the items that are measured, but also to quantify and compare the sizes of differences between them. For example, temperature, as measured in degrees Fahrenheit or Celsius, constitutes an interval scale. From G. Mattioli**Variable Types and Information Content**Ratio variables are very similar to interval variables; in addition to all the properties of interval variables, they feature an identifiable absolute zero point, thus they allow for statements such as x is two times more than y. Typical examples of ratio scales are measures of time or space. From G. Mattioli**Systematic and Random Errors**Error: Defined as the difference between a calculated or observed value and the “true” value From G. Mattioli**Systematic and Random Errors**Blunders: Usually apparent either as obviously incorrect data points or results that are not reasonably close to the expected value. Easy to detect (usually). Easy to fix (throw out data). From G. Mattioli**Systematic and Random Errors**Systematic Errors: Errors that occur reproducibly from faulty calibration of equipment or observer bias. Statistical analysis in generally not useful, but rather corrections must be made based on experimental conditions. From G. Mattioli**Systematic and Random Errors**Random Errors: Errors that result from the fluctuations in observations. Requires that experiments be repeated a sufficient number of time to establish the precision of measurement. (statistics useful here) From G. Mattioli**Accuracy vs. Precision**From G. Mattioli**Accuracy vs. Precision**Accuracy: A measure of how close an experimental result is to the true value. Precision: A measure of how exactly the result is determined. It is also a measure of how reproducible the result is. From G. Mattioli**Accuracy vs. Precision**Absolute precision: indicates the uncertainty in the same units as the observation Relative precision: indicates the uncertainty in terms of a fraction of the value of the result From G. Mattioli**Uncertainties**In most cases, cannot know what the “true” value is unless there is an independent determination (i.e. different measurement technique). From G. Mattioli**Uncertainties**Only can consider estimates of the error. Discrepancy is the difference between two or more observations. This gives rise to uncertainty. Probable Error: Indicates the magnitude of the error we estimate to have made in the measurements. Means that if we make a measurement that we will be wrong by that amount “on average”. From G. Mattioli**Parent vs. Sample Populations**Parent population: Hypothetical probability distribution if we were to make an infinite number of measurements of some variable or set of variables. From G. Mattioli**Parent vs. Sample Populations**• Sample population: • Actual set of experimental observations or measurements of some variable or set of variables. • In General: • (Parent Parameter) = limit (Sample Parameter) • When the number of observations, N, goes to infinity. From G. Mattioli**some univariate statistical terms:**mode: value that occurs most frequently in a distribution (usually the highest point of curve) may have more than one mode (eg. Bimodal – example later) in a dataset From G. Mattioli