16 Conformal Mapping. Overview
|
|
- Evangeline Carpenter
- 5 years ago
- Views:
Transcription
1 Algebraic Geodesy and Geoinformatics PART II APPLICATIONS 16 Conformal Mapping Overview First, the 3- point problem is discussed. A preliminary elimination of the translation vector reduces the size of system to 4 polynomial equations. This system can be solved in symbolic way, either with Dixon resultant or with reduced Groebner basis, both methods result in the same quartic univariate polynomial. Concerning numerical solution Extended NewtonRaphson method can be employed to solve all of the 9 equations in least square sense, utilizing the symbolic solution of any 7 equations as initial guess. For N- point problem Gauss- Jacobi solution improved by Extended Newton- Raphson method is very attractive. However, the General Procrustes algorithm is the fastest and precise as the direct global minimization technique requiring about 5 times longer computation time, than the General Procrustes. All numerical data of the examples are from the text book Problem definition Let us consider coordinates given in two systems A and B. The coordinates of the same physical point, Pi in system A are (Xi, Yi, Zi ), while its corresponding coordinates in system B are (xi, yi, zi ). We suppose that the relation between the two systems can be described by conformal mapping, namely Xi X0 xi yi = s R Yi + Y0 zi Zi Z0 This formula represents 3 elementary transformations, - scaling, with positive, real s, X0 - translation, with vector Y0, Z0 - rotation, with matrix R. The rotation matrix, R can be expressed by the skew matrix, S R = HI3 - SL-1 HI3 + SL where Clear@"Global *"D
2 2 ConformalMapping_16.nb S 0 c b c 0 a b a 0 ; and I 3 IdentityMatrix 3 ; Then the rotation matrix, R Inverse I 3 S. I 3 S Simplify; MatrixForm R 1 a 2 b 2 c 2 1 a 2 b 2 c 2 2 a b 2 c 1 a 2 b 2 c 2 2 a b c 1 a 2 b 2 c 2 1 a 2 b 2 c 2 1 a 2 b 2 c 2 2 b a c 1 a 2 b 2 c 2 2 b a c 1 a 2 b 2 c 2 2 a b c 1 a 2 b 2 c 2 2 a b c 1 a 2 b 2 c 2 1 a 2 b 2 c 2 1 a 2 b 2 c 2 which satisfies the following relation, I 3 R.Transpose R Simplify True The transformation has 7 parameters, (a, b, c, X 0, Y 0, Z 0, s) and to determine them, we need the coordinates of minimum 3 points in both systems (3 -Point Problem). The prototype equation for a point, P i, f 3 i 2 f 3 i 1 I 3 S. f 3 i x i y i z i s I 3 S. X i Y i Z i I 3 S. X 0 Y 0 Z 0 Expand; MatrixForm f 3 i 2 f 3 i 1 f 3 i Then for i = 1 x i X 0 s X i c y i c Y 0 c s Y i b z i b Z 0 b s Z i c x i c X 0 c s X i y i Y 0 s Y i a z i a Z 0 a s Z i b x i b X 0 b s X i a y i a Y 0 a s Y i z i Z 0 s Z i f 1 f 2 f 3 f 3 i 2 f 3 i 1 f 3 i. i 1 Expand; MatrixForm f 1 f 2 f 3 for i = 2 x 1 X 0 s X 1 c y 1 c Y 0 c s Y 1 b z 1 b Z 0 b s Z 1 c x 1 c X 0 c s X 1 y 1 Y 0 s Y 1 a z 1 a Z 0 a s Z 1 b x 1 b X 0 b s X 1 a y 1 a Y 0 a s Y 1 z 1 Z 0 s Z 1 f 4 f 5 f 6 f 3 i 2 f 3 i 1 f 3 i. i 2 Expand; MatrixForm f 4 f 5 f 6 for i = 3 x 2 X 0 s X 2 c y 2 c Y 0 c s Y 2 b z 2 b Z 0 b s Z 2 c x 2 c X 0 c s X 2 y 2 Y 0 s Y 2 a z 2 a Z 0 a s Z 2 b x 2 b X 0 b s X 2 a y 2 a Y 0 a s Y 2 z 2 Z 0 s Z 2
3 ConformalMapping_16.nb f 7 f 8 f 9 f 3 i 2 f 3 i 1 f 3 i. i 3 Expand; MatrixForm f 7 f 8 f 9 x 3 X 0 s X 3 c y 3 c Y 0 c s Y 3 b z 3 b Z 0 b s Z 3 c x 3 c X 0 c s X 3 y 3 Y 0 s Y 3 a z 3 a Z 0 a s Z 3 b x 3 b X 0 b s X 3 a y 3 a Y 0 a s Y 3 z 3 Z 0 s Z 3 Translation parameters, (X 0, Y 0, Z 0 ) can be eliminated by differencing, f 14 f 1 f 4 Simplify x 1 x 2 s X 1 s X 2 c y 1 c y 2 c s Y 1 c s Y 2 b z 1 b z 2 b s Z 1 b s Z 2 f 25 f 2 f 5 Simplify c x 1 c x 2 c s X 1 c s X 2 y 1 y 2 s Y 1 s Y 2 a z 1 a z 2 a s Z 1 a s Z 2 f 39 f 3 f 9 Simplify b x 1 b x 3 b s X 1 b s X 3 a y 1 a y 3 a s Y 1 a s Y 3 z 1 z 3 s Z 1 s Z 3 f 69 f 6 f 9 Simplify b x 2 b x 3 b s X 2 b s X 3 a y 2 a y 3 a s Y 2 a s Y 3 z 2 z 3 s Z 2 s Z 3 Now, we have four equations and four unknown parameters (a, b, c, s). The nonlinearity is represented by the variable s only. Let us introduce new variables, newvars x12 x 1 x 2, x13 x 1 x 3, x23 x 2 x 3, y12 y 1 y 2, y13 y 1 y 3, y23 y 2 y 3, z12 z 1 z 2, z13 z 1 z 3, z23 z 2 z 3, X12 X 1 X 2, X13 X 1 X 3, X23 X 2 X 3, Y12 Y 1 Y 2, Y13 Y 1 Y 3, Y23 Y 2 Y 3, Z12 Z 1 Z 2, Z13 Z 1 Z 3, Z23 Z 2 Z 3 ; Then our system becomes, sys x12 s X12 c y12 c s Y12 b z12 b s Z12, c x12 c s X12 y12 s Y12 a z12 a s Z12, b x13 b s X13 a y13 a s Y13 z13 s Z13, b x23 b s X23 a y23 a s Y23 z23 s Z23 ; Let us check it, f 14, f 25, f 39, f 69 sys. newvars Simplify 0, 0, 0, Symbolic solution Dixon Resultant Resultant Dixon We eliminate the linear parameters, a, b, and c, in order to get an univariate polynomial of s,
4 4 ConformalMapping_16.nb AbsoluteTiming solsdx DixonResultant sys, a, b, c, U, V, W Simplify; , Null solsdx x12 2 x23 y13 s Y13 x13 y23 s X23 y13 X13 y23 x13 Y23 s 2 X23 Y13 X13 Y23 y12 x23 y12 y13 z12 z13 x13 y12 y23 z12 z23 s X13 y12 2 y23 x13 y12 2 Y23 X12 y23 z12 z13 X23 y12 y12 y13 z12 z13 x23 y12 2 Y13 Y12 z12 z13 y12 Z12 z13 y12 z12 Z13 X13 y12 z12 z23 x13 Y12 z12 z23 X12 y13 z12 z23 x13 y12 Z12 z23 x13 y12 z12 Z23 x12 z12 s Z12 y23 z13 s Z13 y13 z23 s Y23 z13 Y13 z23 y13 Z23 s 2 Y23 Z13 Y13 Z23 s 2 X23 y12 2 Y13 x13 Y12 2 y23 X12 2 x23 y13 x13 y23 X13 y12 2 Y23 X23 Y12 z12 z13 X23 y12 Z12 z13 X23 y12 z12 Z13 x23 Y12 2 y13 Y12 Z12 z13 Y12 z12 Z13 y12 Z12 Z13 X13 Y12 z12 z23 X13 y12 Z12 z23 x13 Y12 Z12 z23 X13 y12 z12 Z23 x13 Y12 z12 Z23 x13 y12 Z12 Z23 X12 Y23 z12 z13 y23 Z12 z13 y23 z12 Z13 Y13 z12 z23 y13 Z12 z23 y13 z12 Z23 s 3 x23 Y12 2 Y13 X13 Y12 2 y23 x13 Y12 2 Y23 X12 2 X23 y13 x23 Y13 X13 y23 x13 Y23 x23 Y12 Z12 Z13 X23 Y12 2 y13 Y12 Z12 z13 Y12 z12 Z13 y12 Z12 Z13 X13 Y12 Z12 z23 X13 Y12 z12 Z23 X13 y12 Z12 Z23 x13 Y12 Z12 Z23 X12 Y23 Z12 z13 Y23 z12 Z13 y23 Z12 Z13 Y13 Z12 z23 Y13 z12 Z23 y13 Z12 Z23 s 4 X12 2 X23 Y13 X13 Y23 X12 Y23 Z12 Z13 Y13 Z12 Z23 Y12 X23 Y12 Y13 Z12 Z13 X13 Y12 Y23 Z12 Z23 The result is a quartic polynomial of s, Exponent[solsdx, {s,a,b,c}] 4, 0, 0, 0 The coeffients, q0s Coefficient solsdx, s, 0 Simplify x12 2 x23 y13 x13 y23 x12 z12 y23 z13 y13 z23 y12 x23 y12 y13 z12 z13 x13 y12 y23 z12 z23 q1s Coefficient solsdx, s Simplify x23 y12 2 Y13 X13 y12 2 y23 x13 y12 2 Y23 x12 2 X23 y13 x23 Y13 X13 y23 x13 Y23 x23 Y12 z12 z13 X12 y23 z12 z13 x23 y12 Z12 z13 X23 y12 y12 y13 z12 z13 x23 y12 z12 Z13 X13 y12 z12 z23 x13 Y12 z12 z23 X12 y13 z12 z23 x13 y12 Z12 z23 x13 y12 z12 Z23 x12 Y23 z12 z13 y23 Z12 z13 y23 z12 Z13 Y13 z12 z23 y13 Z12 z23 y13 z12 Z23 q2s Coefficient solsdx, s 2 Simplify x12 2 X23 Y13 X23 y12 2 Y13 x13 Y12 2 y23 X12 2 x23 y13 x13 y23 x12 2 X13 Y23 X13 y12 2 Y23 X23 Y12 z12 z13 X23 y12 Z12 z13 x12 Y23 Z12 z13 X23 y12 z12 Z13 x12 Y23 z12 Z13 x12 y23 Z12 Z13 x23 Y12 2 y13 Y12 Z12 z13 Y12 z12 Z13 y12 Z12 Z13 X13 Y12 z12 z23 X13 y12 Z12 z23 x13 Y12 Z12 z23 x12 Y13 Z12 z23 X13 y12 z12 Z23 x13 Y12 z12 Z23 x12 Y13 z12 Z23 x13 y12 Z12 Z23 x12 y13 Z12 Z23 X12 Y23 z12 z13 y23 Z12 z13 y23 z12 Z13 Y13 z12 z23 y13 Z12 z23 y13 z12 Z23 q3s Coefficient solsdx, s 3 Simplify x23 Y12 2 Y13 X13 Y12 2 y23 x13 Y12 2 Y23 X12 2 X23 y13 x23 Y13 X13 y23 x13 Y23 x23 Y12 Z12 Z13 x12 Y23 Z12 Z13 X23 Y12 2 y13 Y12 Z12 z13 Y12 z12 Z13 y12 Z12 Z13 X13 Y12 Z12 z23 X13 Y12 z12 Z23 X13 y12 Z12 Z23 x13 Y12 Z12 Z23 x12 Y13 Z12 Z23 X12 Y23 Z12 z13 Y23 z12 Z13 y23 Z12 Z13 Y13 Z12 z23 Y13 z12 Z23 y13 Z12 Z23
5 ConformalMapping_16.nb q4s Coefficient solsdx, s 4 Simplify X12 2 X23 Y13 X13 Y23 X12 Y23 Z12 Z13 Y13 Z12 Z23 Y12 X23 Y12 Y13 Z12 Z13 X13 Y12 Y23 Z12 Z23 Let us check the result, solsdx q4s s 4 q3s s 3 q2s s 2 q1s s q0s Simplify True The other parameters can be computed as function of the nonlinear parameter s. The first equation does not contain the parameter a, Map Exponent, a, b, c &, sys 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0 Therefore, we do not take it into consideration when computing a, soladx DixonResultant Drop sys, 1, 1, b, c, V, W Simplify x12 s X12 a s X23 y13 s 2 X23 Y13 x23 y13 s Y13 x13 y23 s X13 y23 s x13 Y23 s 2 X13 Y23 x23 z13 s X23 z13 s x23 Z13 s 2 X23 Z13 x13 z23 s X13 z23 s x13 Z23 s 2 X13 Z23 Similarly, we leave out equations, which does not contain the corresponding parameter in case of b, and c, solbdx DixonResultant Drop sys, 2, 2, a, c, U, W Simplify y12 s Y12 b s X23 y13 s 2 X23 Y13 x23 y13 s Y13 x13 y23 s X13 y23 s x13 Y23 s 2 X13 Y23 y23 z13 s Y23 z13 s y23 Z13 s 2 Y23 Z13 y13 z23 s Y13 z23 s y13 Z23 s 2 Y13 Z23 solcdx DixonResultant Drop sys, 3, 3, a, b, U, V Simplify z12 s Z12 c x23 y12 y12 y23 x12 x23 s X23 c y23 c s Y23 z12 z23 s c X23 y12 c x23 Y12 Y12 y23 X12 x23 c y23 y12 Y23 Z12 z23 z12 Z23 s 2 c X23 Y12 Y12 Y23 X12 X23 c Y23 Z12 Z Reduced Groebner Basis The similar result can be achieved by reduced Groebner basis. AbsoluteTiming solsgb GroebnerBasis sys, s, a, b, c, a, b, c Simplify; , Null
6 6 ConformalMapping_16.nb solsgb x12 2 s X23 y13 s 2 X23 Y13 x23 y13 s Y13 x13 y23 s X13 y23 s x13 Y23 s 2 X13 Y23 y12 x23 y12 y13 z12 z13 x13 y12 y23 z12 z23 s X13 y12 2 y23 x13 y12 2 Y23 X12 y23 z12 z13 X23 y12 y12 y13 z12 z13 x23 y12 2 Y13 Y12 z12 z13 y12 Z12 z13 y12 z12 Z13 X13 y12 z12 z23 x13 Y12 z12 z23 X12 y13 z12 z23 x13 y12 Z12 z23 x13 y12 z12 Z23 x12 z12 s Z12 y23 z13 s Z13 y13 z23 s Y23 z13 Y13 z23 y13 Z23 s 2 Y23 Z13 Y13 Z23 s 2 X23 y12 2 Y13 x13 Y12 2 y23 X12 2 x23 y13 x13 y23 X13 y12 2 Y23 X23 Y12 z12 z13 X23 y12 Z12 z13 X23 y12 z12 Z13 x23 Y12 2 y13 Y12 Z12 z13 Y12 z12 Z13 y12 Z12 Z13 X13 Y12 z12 z23 X13 y12 Z12 z23 x13 Y12 Z12 z23 X13 y12 z12 Z23 x13 Y12 z12 Z23 x13 y12 Z12 Z23 X12 Y23 z12 z13 y23 Z12 z13 y23 z12 Z13 Y13 z12 z23 y13 Z12 z23 y13 z12 Z23 s 3 x23 Y12 2 Y13 X13 Y12 2 y23 x13 Y12 2 Y23 X12 2 X23 y13 x23 Y13 X13 y23 x13 Y23 x23 Y12 Z12 Z13 X23 Y12 2 y13 Y12 Z12 z13 Y12 z12 Z13 y12 Z12 Z13 X13 Y12 Z12 z23 X13 Y12 z12 Z23 X13 y12 Z12 Z23 x13 Y12 Z12 Z23 X12 Y23 Z12 z13 Y23 z12 Z13 y23 Z12 Z13 Y13 Z12 z23 Y13 z12 Z23 y13 Z12 Z23 s 4 X12 2 X23 Y13 X13 Y23 X12 Z12 Y23 Z13 Y13 Z23 Y12 X23 Y12 Y13 Z12 Z13 X13 Y12 Y23 Z12 Z23 This is the same result as the result of the Dixon resultant, but the computation time is considerably longer. solsgb solsdx Simplify True The determination of the other parameters are similar. Again, we consider the parameter s as a constant parameter, solagb GroebnerBasis Drop sys, 1, 1, a, b, c, b, c a x23 y13 a s X23 y13 a s x23 Y13 a s 2 X23 Y13 a x13 y23 a s X13 y23 a s x13 Y23 a s 2 X13 Y23 x23 z13 s X23 z13 s x23 Z13 s 2 X23 Z13 x13 z23 s X13 z23 s x13 Z23 s 2 X13 Z23 The Dixon solution has two factors. The reduced Groebner basis gives the second one, solagb soladx 2 Simplify True Similarly solbgb GroebnerBasis Drop sys, 2, 2, a, b, c, a, c b x23 y13 b s X23 y13 b s x23 Y13 b s 2 X23 Y13 b x13 y23 b s X13 y23 b s x13 Y23 b s 2 X13 Y23 y23 z13 s Y23 z13 s y23 Z13 s 2 Y23 Z13 y13 z23 s Y13 z23 s y13 Z23 s 2 Y13 Z23 solbgb solbdx 2 Simplify True solcgb GroebnerBasis Drop sys, 3, 3, a, b, c, a, b x12 x23 s X12 x23 s x12 X23 s 2 X12 X23 c x23 y12 c s X23 y12 c s x23 Y12 c s 2 X23 Y12 c x12 y23 c s X12 y23 y12 y23 s Y12 y23 c s x12 Y23 c s 2 X12 Y23 s y12 Y23 s 2 Y12 Y23 z12 z23 s Z12 z23 s z12 Z23 s 2 Z12 Z23
7 ConformalMapping_16.nb solcgb solcdx 2 Simplify True Computation of the translation vector The translation vector can be computed from the original system consisting of 9 equations, syso Table f i, i, 1, 9 x 1 X 0 s X 1 c y 1 c Y 0 c s Y 1 b z 1 b Z 0 b s Z 1, c x 1 c X 0 c s X 1 y 1 Y 0 s Y 1 a z 1 a Z 0 a s Z 1, b x 1 b X 0 b s X 1 a y 1 a Y 0 a s Y 1 z 1 Z 0 s Z 1, x 2 X 0 s X 2 c y 2 c Y 0 c s Y 2 b z 2 b Z 0 b s Z 2, c x 2 c X 0 c s X 2 y 2 Y 0 s Y 2 a z 2 a Z 0 a s Z 2, b x 2 b X 0 b s X 2 a y 2 a Y 0 a s Y 2 z 2 Z 0 s Z 2, x 3 X 0 s X 3 c y 3 c Y 0 c s Y 3 b z 3 b Z 0 b s Z 3, c x 3 c X 0 c s X 3 y 3 Y 0 s Y 3 a z 3 a Z 0 a s Z 3, b x 3 b X 0 b s X 3 a y 3 a Y 0 a s Y 3 z 3 Z 0 s Z 3 Concerning the translation parameters, we have only 3 unknown parameters, but 9 equations. The translation vector, for i = 1, 2, 3. Therefore the coefficient matrix has special structure, namely Flatten Table IdentityMatrix 3, 3, 1 ; MatrixForm the pseudoinverze of,
8 8 ConformalMapping_16.nb pi PseudoInverse ; pi MatrixForm Therefore, the least square solution is a simple averaging, see Gauss-Jacobi combinatorial solution, see Chapter 7. The solution, pi. Α 1 Β 1 Γ 1 Α 2 Β 2 Γ 2 Α 3 Β 3 Γ 3 Α 1 3 Α 2 3 Α 3 3, Β 1 3 Β 2 3 Β 3 3, Γ 1 3 Γ 2 3 Γ 3 3 or in detailed form, solxyz0 1 3 x 1 y 1 z 1 s R. X 1 Y 1 Z 1 x 2 y 2 z 2 s R. X 2 Y 2 Z 2 x 3 y 3 z 3 s R. X 3 Y 3 Z 3 Simplify 1 1 a 2 b 2 c 2 x 1 1 a 2 b 2 c 2 x 2 x a 2 b 2 c 2 a 2 x 3 b 2 x 3 c 2 x 3 s X 1 a 2 s X 1 b 2 s X 1 c 2 s X 1 s X 2 a 2 s X 2 b 2 s X 2 c 2 s X 2 s X 3 a 2 s X 3 b 2 s X 3 c 2 s X 3 2 a b s Y 1 2 c s Y 1 2 a b s Y 2 2 c s Y 2 2 a b s Y 3 2 c s Y 3 2 b s Z 1 2 a c s Z 1 2 b s Z 2 2 a c s Z 2 2 b s Z 3 2 a c s Z 3, 1 2 a b c s X 1 2 a b c s X 2 2 a b s X 3 2 c s X 3 y a 2 b 2 c 2 a 2 y 1 b 2 y 1 c 2 y 1 y 2 a 2 y 2 b 2 y 2 c 2 y 2 y 3 a 2 y 3 b 2 y 3 c 2 y 3 s Y 1 a 2 s Y 1 b 2 s Y 1 c 2 s Y 1 s Y 2 a 2 s Y 2 b 2 s Y 2 c 2 s Y 2 s Y 3 a 2 s Y 3 b 2 s Y 3 c 2 s Y 3 2 a s Z 1 2 b c s Z 1 2 a s Z 2 2 b c s Z 2 2 a s Z 3 2 b c s Z 3, 1 2 b a c s X 1 2 b a c s X 2 2 b s X 3 2 a c s X 3 2 a s Y a 2 b 2 c 2 2 b c s Y 1 2 a s Y 2 2 b c s Y 2 2 a s Y 3 2 b c s Y 3 z 1 a 2 z 1 b 2 z 1 c 2 z 1 z 2 a 2 z 2 b 2 z 2 c 2 z 2 z 3 a 2 z 3 b 2 z 3 c 2 z 3 s Z 1 a 2 s Z 1 b 2 s Z 1 c 2 s Z 1 s Z 2 a 2 s Z 2 b 2 s Z 2 c 2 s Z 2 s Z 3 a 2 s Z 3 b 2 s Z 3 c 2 s Z Numerical Example for the 3-Point Problem Let us consider the following data of 3 physical points, datac3 X , Y , Z , X , Y , Z , X , Y , Z , x , y , z , x , y , z , x , y , z ;
9 ConformalMapping_16.nb Let us compute the parameter s. The quartic polynomial for s, eqs s 4 q4s s 3 q3s s 2 q2s s q1s q0s. newvars. datac s s s s 4 Let us normalize it, eqs eqs Coefficient eqs, s, 4 Expand s s s 3 1. s 4 The roots, sols NSolve eqs, s Flatten s , s , s , s The admissible solution should be positive, real, s0 Select sols, Im & s The elements of the skew matrix can be directly computed from the symbolic result, for example let us take the result of the Dixon resultant, then a0 Solve soladx. newvars 0, a Simplify a s X 2 z 1 s X 3 z 1 x 1 z 2 s X 1 z 2 s X 3 z 2 x 1 z 3 s X 1 z 3 s X 2 z 3 s 2 X 2 Z 1 s 2 X 3 Z 1 s x 1 Z 2 s 2 X 1 Z 2 s 2 X 3 Z 2 x 3 z 1 z 2 s Z 1 s Z 2 s x 1 Z 3 s 2 X 1 Z 3 s 2 X 2 Z 3 x 2 z 1 z 3 s Z 1 s Z 3 s X 2 y 1 s X 3 y 1 x 1 y 2 s X 1 y 2 s X 3 y 2 x 1 y 3 s X 1 y 3 s X 2 y 3 s 2 X 2 Y 1 s 2 X 3 Y 1 s x 1 Y 2 s 2 X 1 Y 2 s 2 X 3 Y 2 x 3 y 1 y 2 s Y 1 s Y 2 s x 1 Y 3 s 2 X 1 Y 3 s 2 X 2 Y 3 x 2 y 1 y 3 s Y 1 s Y 3 a0 a0. s0. datac3 Flatten a similarly b0 Solve solbdx. newvars 0, b Simplify b s Y 2 z 1 s Y 3 z 1 y 1 z 2 s Y 1 z 2 s Y 3 z 2 y 1 z 3 s Y 1 z 3 s Y 2 z 3 s 2 Y 2 Z 1 s 2 Y 3 Z 1 s y 1 Z 2 s 2 Y 1 Z 2 s 2 Y 3 Z 2 y 3 z 1 z 2 s Z 1 s Z 2 s y 1 Z 3 s 2 Y 1 Z 3 s 2 Y 2 Z 3 y 2 z 1 z 3 s Z 1 s Z 3 s X 2 y 1 s X 3 y 1 x 1 y 2 s X 1 y 2 s X 3 y 2 x 1 y 3 s X 1 y 3 s X 2 y 3 s 2 X 2 Y 1 s 2 X 3 Y 1 s x 1 Y 2 s 2 X 1 Y 2 s 2 X 3 Y 2 x 3 y 1 y 2 s Y 1 s Y 2 s x 1 Y 3 s 2 X 1 Y 3 s 2 X 2 Y 3 x 2 y 1 y 3 s Y 1 s Y 3 b0 b0. s0. datac3 Flatten b
10 10 ConformalMapping_16.nb c0 Solve solcdx. newvars 0, c Simplify c x 2 2 s x 3 X 1 s x 3 X 2 s 2 X 1 X 2 s 2 X 2 2 s 2 X 1 X 3 s 2 X 2 X 3 x 2 x 3 s X 1 s X 3 x 1 x 2 x 3 s X 2 s X 3 y 1 y 2 y 2 2 y 1 y 3 y 2 y 3 s y 2 Y 1 s y 3 Y 1 s y 1 Y 2 s y 3 Y 2 s 2 Y 1 Y 2 s 2 Y 2 2 s y 1 Y 3 s y 2 Y 3 s 2 Y 1 Y 3 s 2 Y 2 Y 3 z 1 z 2 z 2 2 z 1 z 3 z 2 z 3 s z 2 Z 1 s z 3 Z 1 s z 1 Z 2 s z 3 Z 2 s 2 Z 1 Z 2 s 2 Z 2 2 s z 1 Z 3 s z 2 Z 3 s 2 Z 1 Z 3 s 2 Z 2 Z 3 s X 2 y 1 s X 3 y 1 x 1 y 2 s X 1 y 2 s X 3 y 2 x 1 y 3 s X 1 y 3 s X 2 y 3 s 2 X 2 Y 1 s 2 X 3 Y 1 s x 1 Y 2 s 2 X 1 Y 2 s 2 X 3 Y 2 x 3 y 1 y 2 s Y 1 s Y 2 s x 1 Y 3 s 2 X 1 Y 3 s 2 X 2 Y 3 x 2 y 1 y 3 s Y 1 s Y 3 c0 c0. s0. datac3 Flatten c The translation vector, XYZ0 MapThread 1 2 &, X 0, Y 0, Z 0, solxyz0. datac3. s0. a0. b0. c0 Flatten X , Y , Z The residium of the equations of the original system, rs syso. datac3. s0. a0. b0. c0. XYZ , , , , , , , , The residium of the solution, Norm This method is implemented as a Mathematica function, Conform3DV7, in the GeoAlgebra package, GeoAlgebra Conform3DV7? Conform3DV7 Solves the 7 Parameter Datum Transformation Problem computing the best fitting parameters of the linear transform, between systems X,Y,Z x,y,z in form x,y,z s R a, b, c X,Y,Z X0,Y0,Z0. The inputs: xyz x1,y1,z1, x2,y2,z2, x3,y3,z3 XYZ X1,Y1,Z1, X1,Y1,Z1, X1,Y1,Z1. The result: s, a, b, c, X0, Y0, Z0 scale,s, elements of the skew matrix a, b, c, translation vector, X0, Y0, Z0 Let us employ the function,
11 ConformalMapping_16.nb AbsoluteTiming sol Conform3DV , , , , , , , , , , , , , , , , , NumberForm, 12 &; 0., Null sol , , , , , , The function is very fast, thanks for the symbolic solution! 16-4 Numerical Solutions First, we start with the Global Numerical Solver. Using it as a global method, any 7 equations can be solved from the 9 ones. The system of the 9 equations, syso x 1 X 0 s X 1 c y 1 c Y 0 c s Y 1 b z 1 b Z 0 b s Z 1, c x 1 c X 0 c s X 1 y 1 Y 0 s Y 1 a z 1 a Z 0 a s Z 1, b x 1 b X 0 b s X 1 a y 1 a Y 0 a s Y 1 z 1 Z 0 s Z 1, x 2 X 0 s X 2 c y 2 c Y 0 c s Y 2 b z 2 b Z 0 b s Z 2, c x 2 c X 0 c s X 2 y 2 Y 0 s Y 2 a z 2 a Z 0 a s Z 2, b x 2 b X 0 b s X 2 a y 2 a Y 0 a s Y 2 z 2 Z 0 s Z 2, x 3 X 0 s X 3 c y 3 c Y 0 c s Y 3 b z 3 b Z 0 b s Z 3, c x 3 c X 0 c s X 3 y 3 Y 0 s Y 3 a z 3 a Z 0 a s Z 3, b x 3 b X 0 b s X 3 a y 3 a Y 0 a s Y 3 z 3 Z 0 s Z 3 Take the first 7 equations, sysor Take syso, 1, 7 x 1 X 0 s X 1 c y 1 c Y 0 c s Y 1 b z 1 b Z 0 b s Z 1, c x 1 c X 0 c s X 1 y 1 Y 0 s Y 1 a z 1 a Z 0 a s Z 1, b x 1 b X 0 b s X 1 a y 1 a Y 0 a s Y 1 z 1 Z 0 s Z 1, x 2 X 0 s X 2 c y 2 c Y 0 c s Y 2 b z 2 b Z 0 b s Z 2, c x 2 c X 0 c s X 2 y 2 Y 0 s Y 2 a z 2 a Z 0 a s Z 2, b x 2 b X 0 b s X 2 a y 2 a Y 0 a s Y 2 z 2 Z 0 s Z 2, x 3 X 0 s X 3 c y 3 c Y 0 c s Y 3 b z 3 b Z 0 b s Z 3 The variable list a, b, c, s, X 0, Y 0, Z 0 ; We have two solutions and the first one has considerable "error", solabcs NSolve sysor. newvars. datac3, a , b , c , s , X , Y , Z , a , b , c , s , X , Y , Z
12 12 ConformalMapping_16.nb Remark : Do not mix the error of the method and the error of the technique, which means how to use it! The reason of the "error" here, is the fact that the system of the 7 equations mathematically has two solutions! Would be the model perfect and the data without error, then the additional two equations were redundant! However, it is generally not true, therefore these solutions do not represent the least square solution of the overdetermined system of 9 equations. In addition again, we can consider the first solution as a parasitic solution, but not as an error of the computation or the method. This phenomenon can be considered as a side effect of the algebraic solution! Now, let us use these results as initial values for the Extended Newton- Raphson method employed for solving the overdetermined system, GeoAlgebra NewtonExtended? NewtonExtended Computes the solution of an overdetermined nonlinear system. Input parameters: f list of functions of the system, x list of variables, x0 list of the initial values, eps error limit for the iteration, default value: 10^ 12 n maximum number of the iterations, default value: 100. Output: list of the iterative solutions We select the second solution as initial values, 0 Map 2 &, solabcs , , , , , , The result is somewhat different from the symbolic solution, AbsoluteTiming soln MapThread 1 2 &,, NewtonExtended syso. datac3,, 0 Last ; , Null soln a , b , c , s , X , Y , Z However, the residium, rn syso. datac3. soln , , , , , , , , is considerably smaller, see Fig.16.1, Norm and it is distributed nearly uniformly,
13 ConformalMapping_16.nb ListPlot rs, rn, Joined True, PlotRange All, Frame True Fig.16.1 Distribution of the residiums in case of symbolic (blue) and numeric solution of the 9 equations in least square sense (maroon) Remark: Probably, the best strategy is to employ symbolic solution first then to improve it with Extended Newton- Raphson method applied to the overdetermined system, to 9 equations. Let us do it! Employing the solution of the symbolic method as initial values, 0 Map 2 &, Join s0, a0, b0, c0, XYZ , , , , , , then applying Extended Newton- Raphson method, soln MapThread 1 2 &,, NewtonExtended syso. datac3,, 0 Last a , b , c , s , X , Y , Z Remark: Generally the best strategy is to combine symbolic and robust local numeric methods to get unique, precise solution without quessing initial values and with just a few iterations! 16-5 N-Point Problem The system of equations and data structures In this case, there are data for more than 3 points at our disposal. The prototype equation, e I 3 S. x i y i z i s I 3 S. X i Y i Z i I 3 S. X 0 Y 0 Z 0 x i X 0 s X i c y i c Y 0 c s Y i b z i b Z 0 b s Z i, c x i c X 0 c s X i y i Y 0 s Y i a z i a Z 0 a s Z i, b x i b X 0 b s X i a y i a Y 0 a s Y i z i Z 0 s Z i Expand Flatten In our illlustrative example there are seven points. The system in case of these 7 points,
14 14 ConformalMapping_16.nb sys Table e, i, 1, 7 Flatten x 1 X 0 s X 1 c y 1 c Y 0 c s Y 1 b z 1 b Z 0 b s Z 1, c x 1 c X 0 c s X 1 y 1 Y 0 s Y 1 a z 1 a Z 0 a s Z 1, b x 1 b X 0 b s X 1 a y 1 a Y 0 a s Y 1 z 1 Z 0 s Z 1, x 2 X 0 s X 2 c y 2 c Y 0 c s Y 2 b z 2 b Z 0 b s Z 2, c x 2 c X 0 c s X 2 y 2 Y 0 s Y 2 a z 2 a Z 0 a s Z 2, b x 2 b X 0 b s X 2 a y 2 a Y 0 a s Y 2 z 2 Z 0 s Z 2, x 3 X 0 s X 3 c y 3 c Y 0 c s Y 3 b z 3 b Z 0 b s Z 3, c x 3 c X 0 c s X 3 y 3 Y 0 s Y 3 a z 3 a Z 0 a s Z 3, b x 3 b X 0 b s X 3 a y 3 a Y 0 a s Y 3 z 3 Z 0 s Z 3, x 4 X 0 s X 4 c y 4 c Y 0 c s Y 4 b z 4 b Z 0 b s Z 4, c x 4 c X 0 c s X 4 y 4 Y 0 s Y 4 a z 4 a Z 0 a s Z 4, b x 4 b X 0 b s X 4 a y 4 a Y 0 a s Y 4 z 4 Z 0 s Z 4, x 5 X 0 s X 5 c y 5 c Y 0 c s Y 5 b z 5 b Z 0 b s Z 5, c x 5 c X 0 c s X 5 y 5 Y 0 s Y 5 a z 5 a Z 0 a s Z 5, b x 5 b X 0 b s X 5 a y 5 a Y 0 a s Y 5 z 5 Z 0 s Z 5, x 6 X 0 s X 6 c y 6 c Y 0 c s Y 6 b z 6 b Z 0 b s Z 6, c x 6 c X 0 c s X 6 y 6 Y 0 s Y 6 a z 6 a Z 0 a s Z 6, b x 6 b X 0 b s X 6 a y 6 a Y 0 a s Y 6 z 6 Z 0 s Z 6, x 7 X 0 s X 7 c y 7 c Y 0 c s Y 7 b z 7 b Z 0 b s Z 7, c x 7 c X 0 c s X 7 y 7 Y 0 s Y 7 a z 7 a Z 0 a s Z 7, b x 7 b X 0 b s X 7 a y 7 a Y 0 a s Y 7 z 7 Z 0 s Z 7 The numerical data are, xyz ; and XYZ ; where xyz are the coordinates for the WGS-84 system, while XYZ are the coordinates for the local system. In rule form, xyzr MapThread x 1 2 1, y 1 2 2, z &, Range 7, xyz Flatten x , y , z , x , y , z , x , y , z , x , y , z , x , y , z , x , y , z , x , y , z XYZR MapThread X 1 2 1, Y 1 2 2, Z &, Range 7, XYZ Flatten X , Y , Z , X , Y , Z , X , Y , Z , X , Y , Z , X , Y , Z , X , Y , Z , X , Y , Z
15 ConformalMapping_16.nb Global Minimization The global minimization is the most simple and robust method, but quite time consuming. The system equations in numerical form, sysn sys. xyzr. XYZR; Short sysn, b c s b s c s X 0 c Y 0 b Z 0, a c s a s c s c X 0 Y 0 a Z 0, a b s a s b s b X 0 a Y 0 Z 0, 16, a c s a s c s c X 0 Y 0 a Z 0, a b s a s b s b X 0 a Y 0 Z 0 The objective function, The result, obj Apply Plus, Map ^2 &, sysn ; Short obj, a b s a s b s b X 0 a Y 0 Z a b s a s b s b X 0 a Y 0 Z b c s b s c s X 0 c Y 0 b Z b c s b s c s X 0 c Y 0 b Z 0 2 AbsoluteTiming solgm NMinimize obj, ; , Null solgm , a , b , c , s , X , Y , Z The rotation matrix, Rn R. solgm 2 ; MatrixForm Rn Gauss-Jacobi solution Now, we have 7 points and any 3 of them form a subset, n 3; m 7; The number of the subsets
16 16 ConformalMapping_16.nb mn Binomial m, n 35 This is not a big number, therefore it is reasonable to use combinatorial solution, especially because the solution of the 3- Point Problem is very fast. These subsets are, qs Partition Map &, Flatten Subsets Range m, n, n 1, 2, 3, 1, 2, 4, 1, 2, 5, 1, 2, 6, 1, 2, 7, 1, 3, 4, 1, 3, 5, 1, 3, 6, 1, 3, 7, 1, 4, 5, 1, 4, 6, 1, 4, 7, 1, 5, 6, 1, 5, 7, 1, 6, 7, 2, 3, 4, 2, 3, 5, 2, 3, 6, 2, 3, 7, 2, 4, 5, 2, 4, 6, 2, 4, 7, 2, 5, 6, 2, 5, 7, 2, 6, 7, 3, 4, 5, 3, 4, 6, 3, 4, 7, 3, 5, 6, 3, 5, 7, 3, 6, 7, 4, 5, 6, 4, 5, 7, 4, 6, 7, 5, 6, 7 The data values for the subsets can be generated as it follows, datagj Map Transpose &, Table xyz qs i, j, XYZ qs i, j, i, 1, mn, j, 1, n ; Short datagj, , , , , , , , , , , , , , , , , , , 33, , , , , , , , , , , , , , , , , , Loading the symbolic solution of the 3 - point problem, GeoAlgebra Conform3DV7 ; The solution of the 35 subsets, AbsoluteTiming solgj Map Conform3DV7 1, 2 &, datagj ; , Null
17 ConformalMapping_16.nb The average NumberForm solgj TableForm, solgjavg Map Mean &, Transpose solgj NumberForm, 12 & , , , , , , which is a quite bad result. Now, instead of using weighting technique, we use a different method to improve this result. Let us compute the value of objective function for every subset solution, objgj Map obj. s 1, a 2, b 3, c 4, X 0 5, Y 0 6, Z 0 7 &, solgj , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,
18 18 ConformalMapping_16.nb Select the solution, which has the minimal residual, solgjs solgj First Flatten Position objgj, Min objgj ; solgjs NumberForm, 12 & , , , , , , and do improve it with the Extended Newton- Raphson method, using it as initial guess, Extended Newton- Raphson method GeoAlgebra NewtonExtended AbsoluteTiming solne NewtonExtended sysn,, solgjs Last; , Null solne , , , , , , In order to demonstrate the robustness of the method, let us select the worst result, solgjs solgj First Flatten Position objgj, Max objgj ; solgjs NumberForm, 12 & , , , , , , Again, the solution fast and precise, AbsoluteTiming solne NewtonExtended sysn,, solgjs Last; , Null solne , , , , , , General Procrustes method Another numerical solution technique, the General Procrustes method is also a good candidate for solving the problem, see Chapter 9. We also implemented it as a function in the GeoAlgebra package, GeoAlgebra GeneralProcrustes? GeneralProcrustes Solves the 7 Parameter Datum Transformation Problem computing the scale parameter, s, the translation vector, X0, Y0, Z0 and the rotation matrix, R, as well as the norm of the error matrix, nel. The input data are: Y1 matrix n 3, the coordinates of the image points, xi, yi, zi, Y2 matrix n 3, the coordinates of the object points, Xi, Yi, Zi, W weight matrix n n. Here n is the number of the pairs of points. The output is a list, R, s, X0, Y0, Z0, nel.
19 ConformalMapping_16.nb Now we employ identity matrix for the weigth matrix in order to compare the result with that of the other methods. W IdentityMatrix m 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1 AbsoluteTiming solgp GeneralProcrustes xyz, XYZ, W ; , Null The rotation matrix, solgp 1 MatrixForm The scaling parameter, solgp 2 NumberForm, 12 & The translation vector, The residual, solgp 3 MatrixForm solgp This is the same error norm as was computed with direct global minimization, (Section ), but the computation time is considerably less in that case. Conclusions In case of 3- point problem, Dixon resultant and reduced Groebner basis give the same symbolic result, although the Dixon method is faster than the Groebner method. The computation of the translation vector in least square sense - 3 unknowns and 9 equations - can be achieved by simple averaging because of the linearity of the problem. Employing symbolic solution, a very fast Mathematica function, Conform3DV7 was implemented in the GeoAlgebra package to solve conformal mapping in case of 3-points. Direct numerical solution with Extended Newton- Raphson method is also fast and precise, but it needs initial guess, which can be provided by the Global Numerical Solver, solving any 7 equations from the 9 ones. In all, the best strategy to solve 3-point problem, is the symbolic solution, Conform3DV7 and considering its result as a guess value for the function NewtonExtended. In case of N-point problem, direct global minimization in least square sense is alwasy a good choice, but it can be time consuming. Gauss-Jacobi combinatorial solution is reasonable, if the number of the points - therefore the number of the triplets- is not too high and the technique to solve a triplet is fast. However, the results of the different triplets can differ from each other very considerably. It seems to be the a good strategy to solve a triplet selected randomly, then the result can provide an initial guess for the Extended Newton- Raphson method. Employing General Procrustes algorithm is probably the best choice. It is precise and at least faster than any other methods.
17 Affine Mapping. Overview
Algebraic Geodesy and Geoinformatics - 2009 - PART II APPLICATIONS 17 Affine Mapping Overview The number of equations of the 3 - point problem can be reduced to 6 equations by eliminating the 3 translation
More information13 Positioning by Photogrammetric Resection
Algebraic Geodesy and Geoinformatics - 2009 - PART II APPLICATIONS 13 Positioning by Photogrammetric Resection Overview First the 3- point photogrammetric resection model is solved following the Grafarend-
More informationNow let us consider the solutions of the combinatorial pairs. the weights are the square of the corresponding determinants,
Algebraic Geodesy and Geoinformatics - 2009 - PART I METHODS 7 Gauss- Jacobi Combinatorial Algorithm 7-1 Linear model Another technique to solve overdetermined system is proposed by Gauss and Jacobi. The
More information11 Ranging by Global Navigation Satellite Systems (GNSS)
Algebraic Geodesy and Geoinformatics - 2009 - PART II APPLICATIONS 11 Ranging by Global Navigation Satellite Systems (GNSS) Overview First the observation equations are developed for implicit and explicit
More informationThe original expression can be written as, 4-2 Greatest common divisor of univariate polynomials. Let us consider the following two polynomials.
Algebraic Geodesy and Geoinformatics 2009 PART I - METHODS 4 Groebner Basis 4- Greatest common divisor of integers Let us consider the following integer numbers 2, 20 and 8. Factorize them, Clear@"Global
More information5 Linear Homotopy. 5-1 Introduction
Algebraic Geodesy and Geoinformatics - 2009 PART I METHODS 5 Linear Homotopy 5-1 Introduction Most often, there exist a fundamental task of solving systems of equations in geodesy. In such cases, many
More informationComputational Study of 3D Affine Coordinate Transformation
Computational Study of D Affine Coordinate Transformation Part I. -point Problem Bela Palancz 1, Robert H. Lewis 2, Piroska Zaletnyik and Joseph Awange 4 1 Department of Photogrammetryand Geoinformatics
More informationAlgorithms. Shanks square forms algorithm Williams p+1 Quadratic Sieve Dixon s Random Squares Algorithm
Alex Sundling Algorithms Shanks square forms algorithm Williams p+1 Quadratic Sieve Dixon s Random Squares Algorithm Shanks Square Forms Created by Daniel Shanks as an improvement on Fermat s factorization
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationSolving Algebraic Systems of Equations
Solving Algebraic Systems of Equations Daniel Lichtblau Wolfram Research, Inc. July 2000 Presented at: SCI2000, Orlando Overview Systems of polynomial equations with finitely many solutions arise in many
More informationLINEAR ALGEBRA: NUMERICAL METHODS. Version: August 12,
LINEAR ALGEBRA: NUMERICAL METHODS. Version: August 12, 2000 74 6 Summary Here we summarize the most important information about theoretical and numerical linear algebra. MORALS OF THE STORY: I. Theoretically
More information2 Basics of Polynomials
Algebraic Geodesy and Geoinformatics - 2009 PART I - METHODS 2 Basics of Polynomials 2-1 Representations of Polynomials 2-1- 1 List of monomials In Geodesy and Geoinformatics, most observations are related
More informationUNIT - 2 Unit-02/Lecture-01
UNIT - 2 Unit-02/Lecture-01 Solution of algebraic & transcendental equations by regula falsi method Unit-02/Lecture-01 [RGPV DEC(2013)] [7] Unit-02/Lecture-01 [RGPV JUNE(2014)] [7] Unit-02/Lecture-01 S.NO
More informationAssignment 2.1. Exponent Properties: The Product Rule
Assignment.1 NAME: Exponent Properties: The Product Rule What is the difference between x and x? Explain in complete sentences and with examples. Product Repeated Multiplication Power of the form a b 5
More informationI. Numerical Computing
I. Numerical Computing A. Lectures 1-3: Foundations of Numerical Computing Lecture 1 Intro to numerical computing Understand difference and pros/cons of analytical versus numerical solutions Lecture 2
More informationSomething that can have different values at different times. A variable is usually represented by a letter in algebraic expressions.
Lesson Objectives: Students will be able to define, recognize and use the following terms in the context of polynomials: o Constant o Variable o Monomial o Binomial o Trinomial o Polynomial o Numerical
More informationp 1 p 0 (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- cant method.
80 CHAP. 2 SOLUTION OF NONLINEAR EQUATIONS f (x) = 0 y y = f(x) (p, 0) p 2 p 1 p 0 x (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- Figure 2.16 cant method. Secant Method The
More information1 Number Systems and Errors 1
Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........
More informationData Fitting and Uncertainty
TiloStrutz Data Fitting and Uncertainty A practical introduction to weighted least squares and beyond With 124 figures, 23 tables and 71 test questions and examples VIEWEG+ TEUBNER IX Contents I Framework
More informationEngineering. Mathematics. GATE 2019 and ESE 2019 Prelims. For. Comprehensive Theory with Solved Examples
Thoroughly Revised and Updated Engineering Mathematics For GATE 2019 and ESE 2019 Prelims Comprehensive Theory with Solved Examples Including Previous Solved Questions of GATE (2003-2018) and ESE-Prelims
More informationEE 581 Power Systems. Admittance Matrix: Development, Direct and Iterative solutions
EE 581 Power Systems Admittance Matrix: Development, Direct and Iterative solutions Overview and HW # 8 Chapter 2.4 Chapter 6.4 Chapter 6.1-6.3 Homework: Special Problem 1 and 2 (see handout) Overview
More informationLinear Algebraic Equations
Linear Algebraic Equations Linear Equations: a + a + a + a +... + a = c 11 1 12 2 13 3 14 4 1n n 1 a + a + a + a +... + a = c 21 2 2 23 3 24 4 2n n 2 a + a + a + a +... + a = c 31 1 32 2 33 3 34 4 3n n
More informationQR Decomposition. When solving an overdetermined system by projection (or a least squares solution) often the following method is used:
(In practice not Gram-Schmidt, but another process Householder Transformations are used.) QR Decomposition When solving an overdetermined system by projection (or a least squares solution) often the following
More informationA Classical Introduction to Modern Number Theory
Kenneth Ireland Michael Rosen A Classical Introduction to Modern Number Theory Second Edition Springer Contents Preface to the Second Edition Preface v vii CHAPTER 1 Unique Factorization 1 1 Unique Factorization
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationChapter 9 Factorisation and Discrete Logarithms Using a Factor Base
Chapter 9 Factorisation and Discrete Logarithms Using a Factor Base February 15, 2010 9 The two intractable problems which are at the heart of public key cryptosystems, are the infeasibility of factorising
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationMath 314 Lecture Notes Section 006 Fall 2006
Math 314 Lecture Notes Section 006 Fall 2006 CHAPTER 1 Linear Systems of Equations First Day: (1) Welcome (2) Pass out information sheets (3) Take roll (4) Open up home page and have students do same
More informationGrade AM108R 7 + Mastering the Standards ALGEBRA. By Murney R. Bell
Hayes AM108R Mastering the Standards ALGEBRA By Murney R. Bell Grade 7 + Mastering the Standards Algebra By Murney R. Bell Illustrated by Reneé Yates Copyright 2008, Hayes School Publishing Co., Inc.,
More informationFitting. PHY 688: Numerical Methods for (Astro)Physics
Fitting Fitting Data We get experimental/observational data as a sequence of times (or positions) and associate values N points: (x i, y i ) Often we have errors in our measurements at each of these values:
More informationMathematica Project, Math21b Spring 2008
Mathematica Project, Math21b Spring 2008 Oliver Knill, Harvard University, Spring 2008 Support Welcome to the Mathematica computer project! In case problems with this assignment, please email Oliver at
More informationIntroduction to Applied Linear Algebra with MATLAB
Sigam Series in Applied Mathematics Volume 7 Rizwan Butt Introduction to Applied Linear Algebra with MATLAB Heldermann Verlag Contents Number Systems and Errors 1 1.1 Introduction 1 1.2 Number Representation
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationGATE Engineering Mathematics SAMPLE STUDY MATERIAL. Postal Correspondence Course GATE. Engineering. Mathematics GATE ENGINEERING MATHEMATICS
SAMPLE STUDY MATERIAL Postal Correspondence Course GATE Engineering Mathematics GATE ENGINEERING MATHEMATICS ENGINEERING MATHEMATICS GATE Syllabus CIVIL ENGINEERING CE CHEMICAL ENGINEERING CH MECHANICAL
More informationAlgebra 1 Seamless Curriculum Guide
QUALITY STANDARD #1: REAL NUMBERS AND THEIR PROPERTIES 1.1 The student will understand the properties of real numbers. o Identify the subsets of real numbers o Addition- commutative, associative, identity,
More informationApplied Mathematics 205. Unit I: Data Fitting. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit I: Data Fitting Lecturer: Dr. David Knezevic Unit I: Data Fitting Chapter I.4: Nonlinear Least Squares 2 / 25 Nonlinear Least Squares So far we have looked at finding a best
More informationMATHEMATICS FOR COMPUTER VISION WEEK 2 LINEAR SYSTEMS. Dr Fabio Cuzzolin MSc in Computer Vision Oxford Brookes University Year
1 MATHEMATICS FOR COMPUTER VISION WEEK 2 LINEAR SYSTEMS Dr Fabio Cuzzolin MSc in Computer Vision Oxford Brookes University Year 2013-14 OUTLINE OF WEEK 2 Linear Systems and solutions Systems of linear
More informationComputational Methods. Least Squares Approximation/Optimization
Computational Methods Least Squares Approximation/Optimization Manfred Huber 2011 1 Least Squares Least squares methods are aimed at finding approximate solutions when no precise solution exists Find the
More informationNon-polynomial Least-squares fitting
Applied Math 205 Last time: piecewise polynomial interpolation, least-squares fitting Today: underdetermined least squares, nonlinear least squares Homework 1 (and subsequent homeworks) have several parts
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation
ECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation Prof. om Overbye Dept. of Electrical and Computer Engineering exas A&M University overbye@tamu.edu Announcements
More informationLeast Squares Fitting of Data by Linear or Quadratic Structures
Least Squares Fitting of Data by Linear or Quadratic Structures David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution
More informationPose estimation from point and line correspondences
Pose estimation from point and line correspondences Giorgio Panin October 17, 008 1 Problem formulation Estimate (in a LSE sense) the pose of an object from N correspondences between known object points
More informationSpanning and Independence Properties of Finite Frames
Chapter 1 Spanning and Independence Properties of Finite Frames Peter G. Casazza and Darrin Speegle Abstract The fundamental notion of frame theory is redundancy. It is this property which makes frames
More informationMATH 213 Linear Algebra and ODEs Spring 2015 Study Sheet for Midterm Exam. Topics
MATH 213 Linear Algebra and ODEs Spring 2015 Study Sheet for Midterm Exam This study sheet will not be allowed during the test Books and notes will not be allowed during the test Calculators and cell phones
More informationThe Perrin Conjugate and the Laguerre Orthogonal Polynomial
The Perrin Conjugate and the Laguerre Orthogonal Polynomial In a previous chapter I defined the conjugate of a cubic polynomial G(x) = x 3 - Bx Cx - D as G(x)c = x 3 + Bx Cx + D. By multiplying the polynomial
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationCAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS
CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS Preliminaries Round-off errors and computer arithmetic, algorithms and convergence Solutions of Equations in One Variable Bisection method, fixed-point
More informationCS 221 Lecture 9. Tuesday, 1 November 2011
CS 221 Lecture 9 Tuesday, 1 November 2011 Some slides in this lecture are from the publisher s slides for Engineering Computation: An Introduction Using MATLAB and Excel 2009 McGraw-Hill Today s Agenda
More informationIntroduction to Mobile Robotics Compact Course on Linear Algebra. Wolfram Burgard, Bastian Steder
Introduction to Mobile Robotics Compact Course on Linear Algebra Wolfram Burgard, Bastian Steder Reference Book Thrun, Burgard, and Fox: Probabilistic Robotics Vectors Arrays of numbers Vectors represent
More informationDevelopment of an algorithm for the problem of the least-squares method: Preliminary Numerical Experience
Development of an algorithm for the problem of the least-squares method: Preliminary Numerical Experience Sergey Yu. Kamensky 1, Vladimir F. Boykov 2, Zakhary N. Khutorovsky 3, Terry K. Alfriend 4 Abstract
More informationProfessor Terje Haukaas University of British Columbia, Vancouver Notation
Notation This document establishes the notation that is employed throughout these notes. It is intended as a look-up source during the study of other documents and software on this website. As a general
More informationParameter-Dependent Eigencomputations and MEMS Applications
Parameter-Dependent Eigencomputations and MEMS Applications David Bindel UC Berkeley, CS Division Parameter-Dependent Eigencomputations and MEMS Applications p.1/36 Outline Some applications Mathematical
More informationProcess Model Formulation and Solution, 3E4
Process Model Formulation and Solution, 3E4 Section B: Linear Algebraic Equations Instructor: Kevin Dunn dunnkg@mcmasterca Department of Chemical Engineering Course notes: Dr Benoît Chachuat 06 October
More informationEngineering Mathematics
Thoroughly Revised and Updated Engineering Mathematics For GATE 2017 and ESE 2017 Prelims Note: ESE Mains Electrical Engineering also covered Publications Publications MADE EASY Publications Corporate
More informationThe residual again. The residual is our method of judging how good a potential solution x! of a system A x = b actually is. We compute. r = b - A x!
The residual again The residual is our method of judging how good a potential solution x! of a system A x = b actually is. We compute r = b - A x! which gives us a measure of how good or bad x! is as a
More informationJACOBI S ITERATION METHOD
ITERATION METHODS These are methods which compute a sequence of progressively accurate iterates to approximate the solution of Ax = b. We need such methods for solving many large linear systems. Sometimes
More informationNewton's method for complex polynomials
Newton's method for complex polynomials A preprint version of a Mathematical graphics column from Mathematica in Education and Research. Mark McClure mcmcclure@unca.edu Department of Mathematics University
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
More information3.2 Iterative Solution Methods for Solving Linear
22 CHAPTER 3. NUMERICAL LINEAR ALGEBRA 3.2 Iterative Solution Methods for Solving Linear Systems 3.2.1 Introduction We continue looking how to solve linear systems of the form Ax = b where A = (a ij is
More informationThe Normal Equations. For A R m n with m > n, A T A is singular if and only if A is rank-deficient. 1 Proof:
Applied Math 205 Homework 1 now posted. Due 5 PM on September 26. Last time: piecewise polynomial interpolation, least-squares fitting Today: least-squares, nonlinear least-squares The Normal Equations
More informationThe Jacobian. Jesse van den Kieboom
The Jacobian Jesse van den Kieboom jesse.vandenkieboom@epfl.ch 1 Introduction 1 1 Introduction The Jacobian is an important concept in robotics. Although the general concept of the Jacobian in robotics
More informationNumerical solution of Least Squares Problems 1/32
Numerical solution of Least Squares Problems 1/32 Linear Least Squares Problems Suppose that we have a matrix A of the size m n and the vector b of the size m 1. The linear least square problem is to find
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationWeighted Activity Selection
Weighted Activity Selection Problem This problem is a generalization of the activity selection problem that we solvd with a greedy algorithm. Given a set of activities A = {[l, r ], [l, r ],..., [l n,
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationLecture 44. Better and successive approximations x2, x3,, xn to the root are obtained from
Lecture 44 Solution of Non-Linear Equations Regula-Falsi Method Method of iteration Newton - Raphson Method Muller s Method Graeffe s Root Squaring Method Newton -Raphson Method An approximation to the
More informationChapter 2. Solving Systems of Equations. 2.1 Gaussian elimination
Chapter 2 Solving Systems of Equations A large number of real life applications which are resolved through mathematical modeling will end up taking the form of the following very simple looking matrix
More information1. Let r, s, t, v be the homogeneous relations defined on the set M = {2, 3, 4, 5, 6} by
Seminar 1 1. Which ones of the usual symbols of addition, subtraction, multiplication and division define an operation (composition law) on the numerical sets N, Z, Q, R, C? 2. Let A = {a 1, a 2, a 3 }.
More informationIterative Solvers. Lab 6. Iterative Methods
Lab 6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationNumerical methods part 2
Numerical methods part 2 Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE6103: Week 6 Numerical methods part 2 1/33 Content (week 6) 1 Solution of an
More informationChapter Six. Polynomials. Properties of Exponents Algebraic Expressions Addition, Subtraction, and Multiplication Factoring Solving by Factoring
Chapter Six Polynomials Properties of Exponents Algebraic Expressions Addition, Subtraction, and Multiplication Factoring Solving by Factoring Properties of Exponents The properties below form the basis
More informationDecoding linear codes via systems solving: complexity issues
Decoding linear codes via systems solving: complexity issues Stanislav Bulygin (joint work with Ruud Pellikaan) University of Kaiserslautern June 19, 2008 Outline Outline of the talk Introduction: codes
More informationAn Accelerated Block-Parallel Newton Method via Overlapped Partitioning
An Accelerated Block-Parallel Newton Method via Overlapped Partitioning Yurong Chen Lab. of Parallel Computing, Institute of Software, CAS (http://www.rdcps.ac.cn/~ychen/english.htm) Summary. This paper
More informationBootstrapping, Randomization, 2B-PLS
Bootstrapping, Randomization, 2B-PLS Statistics, Tests, and Bootstrapping Statistic a measure that summarizes some feature of a set of data (e.g., mean, standard deviation, skew, coefficient of variation,
More informationarxiv: v1 [hep-lat] 13 Oct 2017
arxiv:1710.04828v1 [hep-lat] 13 Oct 2017 Satisfying positivity requirement in the Beyond Complex Langevin approach Adam Wyrzykowski 1,, and Błażej Ruba 1, 1 M. Smoluchowski Institute of Physics, Jagiellonian
More informationOnline Courses for High School Students
Online Courses for High School Students 1-888-972-6237 Algebra I Course Description: Students explore the tools of algebra and learn to identify the structure and properties of the real number system;
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationPreliminaries and Complexity Theory
Preliminaries and Complexity Theory Oleksandr Romanko CAS 746 - Advanced Topics in Combinatorial Optimization McMaster University, January 16, 2006 Introduction Book structure: 2 Part I Linear Algebra
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 3 Chapter 10 LU Factorization PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
More informationQuantum Mechanics using Matrix Methods
Quantum Mechanics using Matrix Methods Introduction and the simple harmonic oscillator In this notebook we study some problems in quantum mechanics using matrix methods. We know that we can solve quantum
More informationThe Newton-Raphson method accelerated by using a line search - comparison between energy functional and residual minimization
Physics Electricity & Magnetism fields Okayama University Year 2004 The Newton-Raphson method accelerated by using a line search - comparison between energy functional and residual minimization Koji Fujiwara
More informationHot-Starting NLP Solvers
Hot-Starting NLP Solvers Andreas Wächter Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu 204 Mixed Integer Programming Workshop Ohio
More informationSolving of logic functions systems using genetic algorithm
Solving of logic functions systems using genetic algorithm V G Kurbanov,2 and M V Burakov Chair of control system of Saint-Petersburg State University of Aerospace Instrumentation, Bolshaya Morskaya, 67,
More informationcha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO
cha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO ROOTS OF EQUATIONS PT2.1 MOTIVATION Years ago, you learned to use the quadratic formula x = b ± b 2 4ac 2a to solve f(x) = ax 2 + bx + c = 0 (PT2.1) (PT2.2)
More informationLinear Systems of Differential Equations
Chapter 5 Linear Systems of Differential Equations Project 5. Automatic Solution of Linear Systems Calculations with numerical matrices of order greater than 3 are most frequently carried out with the
More informationKim Bowman and Zach Cochran Clemson University and University of Georgia
Clemson University, Conference October 2004 1/31 Linear Dependency and the Quadratic Sieve Kim Bowman and Zach Cochran Clemson University and University of Georgia and Neil Calkin, Tim Flowers, Kevin James
More informationEPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7
Introduction to Generalized Univariate Models: Models for Binary Outcomes EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7 EPSY 905: Intro to Generalized In This Lecture A short review
More informationMathematica for Calculus II (Version 9.0)
Mathematica for Calculus II (Version 9.0) C. G. Melles Mathematics Department United States Naval Academy December 31, 013 Contents 1. Introduction. Volumes of revolution 3. Solving systems of equations
More informationMATHEMATICS. Course Syllabus. Section A: Linear Algebra. Subject Code: MA. Course Structure. Ordinary Differential Equations
MATHEMATICS Subject Code: MA Course Structure Sections/Units Section A Section B Section C Linear Algebra Complex Analysis Real Analysis Topics Section D Section E Section F Section G Section H Section
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationnonrobust estimation The n measurement vectors taken together give the vector X R N. The unknown parameter vector is P R M.
Introduction to nonlinear LS estimation R. I. Hartley and A. Zisserman: Multiple View Geometry in Computer Vision. Cambridge University Press, 2ed., 2004. After Chapter 5 and Appendix 6. We will use x
More informationSymbolic solution to expansion of natural frequencies method for determining e At
Home PDF version of this note Mathematica notebook Symbolic solution to expansion of natural frequencies method for determining e At Written by Nasser M. Abbasi, based on lecture given by Professor Barmish,
More informationMCS 563 Spring 2014 Analytic Symbolic Computation Monday 14 April. Binomial Ideals
Binomial Ideals Binomial ideals offer an interesting class of examples. Because they occur so frequently in various applications, the development methods for binomial ideals is justified. 1 Binomial Ideals
More informationTopics. The CG Algorithm Algorithmic Options CG s Two Main Convergence Theorems
Topics The CG Algorithm Algorithmic Options CG s Two Main Convergence Theorems What about non-spd systems? Methods requiring small history Methods requiring large history Summary of solvers 1 / 52 Conjugate
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationA Review of Matrix Analysis
Matrix Notation Part Matrix Operations Matrices are simply rectangular arrays of quantities Each quantity in the array is called an element of the matrix and an element can be either a numerical value
More information