arxiv: v2 [math.na] 27 Dec 2016
|
|
- Mervin Chambers
- 5 years ago
- Views:
Transcription
1 An algorithm for constructing Equiangular vectors Azim rivaz a,, Danial Sadeghi a a Department of Mathematics, Shahid Bahonar University of Kerman, Kerman , IRAN arxiv: v2 [math.na] 27 Dec 2016 Abstract In this paper a generalization of the Gram-Schmidt Algorithm is presented. Actually we provide an algorithm to construct a set of equiangular vectors with the prescribed angle θ (0, π/2) using a set of independent vectors in R n. Then a new type of matrix decomposition is derived. Keywords: Gram-Schmidt Algorithm, Equiangular matrix, Equiangular vectors, Law of cosines, Matrix Decomposition 2010 MSC: 15B99, 15A03, 15A23, 15A30 1. Introduction We know that the Gram-Schmidt (GS) Algorithm has numerous applications in applied mathematics such as the theory of least squares, Projection methods and Eigenvalue problems. This algorithm takes a set of linearly independent vectors {a 1, a 2,..., a k } of a n dimensional space (k n), and then generates a set of orthogonal vectors {q 1, q 2,..., q k } successively step-by-step so that they span the same subspaces [1, 2, 6, 7, 8]. If this algorithm is applied to the column vectors of a nonsingular n n matrix A = [a 1,..., a n ], then the obtained vectors q 1, q 2,..., q n are orthogonal and make the orthogonal matrix Q column-wise. The outcome of GS algorithm is the QR decomposition A = QR, where R is an upper triangular matrix. The purpose of this paper is to generalize the GS algorithm by providing an algorithm that generate linearly independent vectors s 1, s 2,..., s n using a 1, a 2,..., a n, provided that s T i s j is equal to α (0, 1) if i j and 1 if i = j. It means the s i s are equiangular with angle θ that θ = arccos α (0, π/2). In this paper we deal with subspaces of R n with the standard inner product and Euclidean norm. Throughout this paper, x denotes the Euclidean norm of vector x and α = cos θ. Also e = [1,..., 1] T R n is the vector of ones and {e 1, e 2,..., e n } is the standard orthogonal basis of R n. Moreover, All matrices in this paper are real. Proposition 1.1. Suppose that three vectors v 1, v 2, v 3 in R 3 are constructed by linear combination of the vector e and the basis {e 1, e 2, e 3 } as v 1 = [1, x, x] T, v 2 = [x, 1, x] T and v 3 = [x, x, 1] T. Then the vectors s i = v i / v i for i = 1, 2, 3 are equiangular with α = cos θ (0, 1), where x = ( (1 α)(1 + 2α) 1)/(1 2α). Proof. The proof is clear by obtaining s T i s j = (x 2 + 2x)/(2x 2 + 1) = α for i j. A well-ordering representation of an equiangular basis in R n analogous to the basis {e 1, e 2,..., e n } is introduced in the following lemma. Proposition 1.2. Let v i = xe + (1 x)e i for i = 1,..., n. Then the set of {s 1,..., s n } where s i = v i / v i (1 α)(1+(n 1)α) 1 are equiangular vectors with α (0, 1) and x = (1 α)(n 1) 1. Moreover if the vectors s 1,..., s n are 1 embedded into the positive coordinate axes so that n 1 entries of all s i are equal to t (0, n ) and the last one is s = 1 (n 1)t 2. Then s = (1 α)(n 1)(n 2)+n±2(n 1) (1 α)(1+(n 1)α) αn+2(1 α (1 α)(1+(n 1)α)) n, t = n, (1.1) so that the plus sign in the formula of s must be selected and minus sign in t. Corresponding author. arivaz@uk.ac.ir. address: dl.sadeghi@math.uk.ac.ir (Danial Sadeghi) Preprint submitted to Elsevier September 3, 2018
2 Proof. We can check that s T i s 2x + (n 2)x2 j = 1 + (n 1)x 2 = α for i j. For the next part we must have s2 +(n 1)t 2 = 1 and α = 2ts + (n 2)t 2. Then we can verify the Equation (1.1). On the other hand, the scalar t tends to zero as α tends to zero, which justifies the minus sign in Equation (1.1) and vise versa in s. Proposition 1.3. Any set of equiangular vectors {s 1,..., s n } in R n with α ( 1 n 1, 1) is linearly independent. Proof. let x 1 s x n s n = 0, for arbitrary scalars x i. Premultiplying by s T i s gives a system of linear equations Gx = 0, where x = [x 1, x 2,..., x n ] T and G has ones on the main diagonal and cos θ on all off-diagonals. G is a Gram matrix of matrix S = [s 1,..., s n ] [4, 5], and G has no zero eigenvalue, so it is nonsingular and x 0. In this case we have θ (0, arccos( 1 1 n 1 )) and α ( n 1, 1). If the upper bound of α holds, then G = eet which is rank-one and singular. If the lower bound of α holds, then G is singular with the eigenpair (0, e) with algebraic multiplicity one at 0. So the rank of G and S is n 1 and then s 1,..., s n are equiangular but linearly dependent. Actually they are equiangular lines [3]. This paper is organized as follows. In section 2, using the basis {e 1, e 2..., e n }, we construct the so-called standard equiangular basis {s 1, s 2..., s n } with the angle θ. Next we construct the equiangular vectors of n- dimensional case using a set of linearly independent vectors. Then we have A = S.R so that the column vectors of A and S are the input and output vectors of the algorithm and R is an upper triangular matrix. Lemma 1.4. Suppose that s 1, s 2,..., s k are equiangular normalized vectors in R n (k n) with α arccos( 1 n 1, 1). Then k s i = k (1 + (k 1)α). Proof. The proof is clear by obtaining k s i 2 = ( k s i) T ( k s i) = k +k(k 1)α = k(1+(k 1)α). The next problem is finding the relationship between the two angles ϕ and η, where ϕ is the angle of the vectors k s i and s k+1, and η is the angle of the vectors k s i and s j (1 j k < n). An interesting triangular relationship is obtained as follows ( k ) ( s T k+1 s k ) i cos ϕ = s k+1 k s i = kα s T j s i, cos η = k (1 + (k 1)α) s j k s i = 1 + (k 1)α. (1.2) k (1 + (k 1)α) Multiplying the above equalities together implies that cos ϕ cos η = cos θ = α. (1.3) Indeed ϕ, η and θ are the angles of a trihedral angle, in such a way that two faces corresponding to ϕ and η are perpendicular together. Actually (1.3) represents the first law of cosine in the spherical pyramid [9]: Lemma 1.5. Let a, b and c are three linearly independent vectors in R n (n 3) which make a trihedral angle with the angles ϕ, η and θ and also the dihedral angle between two faces opposite to θ is π/2 as illustrated in Figure 1. Then cos θ = cos ϕ cos η. Proof. cos θ = at c a c = (proj ba + (a proj b a)) T (proj b c + (c proj b c)) = (proj ba)(proj b c) = cos ϕ cos η. (1.4) a c a c 2. An algorithm for producing a set of equiangular vectors Suppose that {a 1, a 2,..., a n } is a set of linearly independent vectors in R n. In this section we drive a set of n equiangular vectors with angle θ using the given vectors. First, we normalize a 1 as s 1 = a 1 1 a 1. Next, suppose that vector a 2 makes an acute angle with a 1. Let u 2 = a 2 s 1 s T 1 a 2, b = cot θ u 2 s 1 and v 2 = a 2 (s 1 s T 1 a 2 b). Then v 2 = u 2 + cot θ u 2 s 1 (Fig. 2). Without loss of generality we can write v 2 = q 2 + cot θ s 1, (2.1) where q 2 = u 2 1 u 2 is orthogonal to s 1. So the next normalized equiangular vector spanned by {a 1, a 2 } is s 2 = v 2 / v 2. For the case of obtuse angles between a 1 and a 2 we can use the same method with some slight differences. We now ready to present the main theorem as follows 2
3 Figure 1: Three vectors a, b and c form a pyramid, so that two faces which corresponding to ϕ and η are perpendicular. Theorem 2.1. Suppose that the vectors a 1,..., a n in R n are linearly independent. For an angle θ the equiangular vectors s 1, s 2 which are generated from a 1, a 2 by (2.1) are given. If q k+1 span a 1,..., a k+1 is a normalized orthogonal vector to s i s (i = 1,..., k). Then for 2 k n 1 the (k + 1)th equiangular vector is obtained as s k+1 = v k+1 / v k+1 where v k+1 = q k+1 + ( ) k k s i (sec θ 1)(sec θ+k) k s i. (2.2) Proof. We define a subspace named P 2 span s 1, s 2, a 3 as P 2 = {p p span s 1, s 2, a 3, p T s 1 = p T s 2 }. So each vector p of P 2 has the same angle θ p with s 1 and s 2. Indeed P 2 and the subspace spanned by s 1 and s 2 are perpendicular and their intersection is a line spanned by s 1 + s 2. Since p T s 1 = p T s 2, then p T (s 1 s 2 ) = 0 and s 1 s 2 is the orthogonal complement of P 2 so dim(p 2 ) = 2. Let ā 3 be the projection of a 3 onto P 2, thus it makes the same angle with s 1 and s 2. We can write ā 3 = proj P2 a 3 = a 3 proj s1 s 2 a 3 (2.3) From (1.3) the vector named v 3 which makes angle θ with both s 1, s 2, is in subspace P 2 (v 3 P 2 ). Suppose that ϕ is the angle between v 3 and s 1 + s 2. We have ϕ < θ and from the definition of η in (1.3), η = θ/2. Moreover, cot θ is replaced by cot ϕ as follows cos ϕ = cos θ 2 cos θ cos η = = cos θ 1 + cos θ cos(θ/2), sin ϕ = cos2 (θ/2) cos 2 θ 2 cos θ, then cot ϕ = cos(θ/2) 1 + cos θ 2 cos2 θ. The purpose is to compute v 3 using ā 3. In comparison with the previous step (2.1), ā 3 and s 1 + s 2 can be considered as a 2 and s 1, respectively. Let u 3 = ā 3 (s 1 + s 2 )(s 1 + s 2 ) T ā 3 s 1 + s 2 2. Then from (2.3) u 3 = a 3 proj s1 s 2 a 3 (s 1 + s 2 )(s 1 + s 2 ) T (a 3 proj s1 s 2 a 3 ) s 1 + s 2 2 = a 3 proj s1 s 2 a 3 (s 1 + s 2 )(s 1 + s 2 ) T a 3 s 1 + s 2 2 = a 3 (proj s1 s 2 a 3 + proj s1+s 2 a 3 ). (2.4) Note that u 3 is orthogonal to s 1, s 2 because s 1 s 2 s 1 + s 2. Hence from (2.1) we have v 3 = u 3 + cot ϕ u 3 s 1 + s 2. Without loss of generality we can take s 1 + s 2 v 3 = q 3 + cot ϕ s 1 + s 2 s 1 + s 2, (2.5) where q 3 = u 3 1 u 3 which is orthogonal to s 1, s 2. So the next normalized equiangular vector to s 1, s 2 spanned by {a 1, a 2, a 3 } is s 3 = v 3 / v 3. After k steps the equiangular vectors s 1,..., s k will be produced (k 3). For obtaining s k+1 we begin with defining P k as P k = P ij = { p p span s1,..., s k, a k+1, p (s i s j ) ; 1 i < j k }. (2.6) 3
4 Figure 2: Vector v 2 makes the angle θ with s 1. In this case we compute the vector ā k+1 which is the projection of the vector v k+1 onto the subspace P k and then the unit vector s k+1 which makes the same angle θ with the previous s i s will be obtained using ā k+1. For simplicity, we project a k+1 onto the orthogonal subspace of P k, say W k = {s i s j } then using (2.3) in case of (k + 1)-dimensional, implies a k+1 = proj span a1,...,a k+1 a k+1 = proj Wk P k a k+1 = proj Wk a k+1 + proj Pk a k+1 = proj Wk a k+1 + ā k+1, (2.7) and then ā k+1 = a k+1 proj Wk a k+1. (2.8) It remains to obtain proj Wk a k+1. For this we will find an orthogonal basis B Wk = {B k,i } of W k. We claim that W k has a basis of equiangular vectors as S Wk = {S k,i } = {s 1 s k, s 2 s k,..., s s k }. To check the linearly independence take x i(s i s k ) = 0 so x is i ( x i)s k = 0. Since s 1,..., s k are linearly independent, then x 1 = = x = 0. To prove W k = span S Wk it is enough to show that W k span S Wk. Let w W k, then from the definition of W k w = k x i,j (s i s j ) = i<j k x i,j [(s i s k ) (s j s k )] = y i (s i s k ), (2.9) i<j so w span S Wk. This argument yields that dim(w k ) = k 1 so dim(p k ) = k + 1 (k 1) = 2 and P k is a plane. Now we present a method that takes a set of equiangular vectors such as S Wk and generates an orthogonal basis B Wk = {B k,1,..., B k, }. If k = 2 it is trivial: B W2 = S W2 = {s 1 s 2 }. If k = 3 then S W3 = {S 3,1, S 3,2 } = {s 1 s 3, s 2 s 3 }. Since S 3,1 and S 3,2 are equiangular with the same norm, then we choose B W3 as B 3,1 = S 3,1 S 3,2 = s 1 s 2 and B 3,2 = S 3,1 + S 3,2 = s 1 + s 2 2s 3. So B W3 = {s 1 s 2, s 1 + s 2 2s 3 }. If k = 4 then S W4 = {s 1 s 4, s 2 s 4, s 3 s 4 }. Let B 4,3 = S 4,1 + S 4,2 + S 4,3 = s 1 + s 2 + s 3 3s 4. We define C 4,1 = S 4,1 S 4,3 and C 4,2 = S 4,2 S 4,3. Both C 4,1 and C 4,2 are orthogonal to B 4,3, but C 4,1 C 4,2. Now we choose two linear combinations of them which are orthogonal, similar to the case k = 3. So B 4,1 = C 4,1 C 4,2 = S 4,1 S 4,2 = s 1 s 2 and B 4,2 = C 4,1 + C 4,2 = S 4,1 + S 4,2 2S 4,3 = s 1 + s 2 2s 3 and then B W4 = {s 1 s 2, s 1 + s 2 2s 3, s 1 + s 2 + s 3 3s 4 }. For k > 4, B Wk is obtained recurrently, e.g., if k = 5, then B 5,1 = s 1 s 2, B 5,2 = s 1 + s 2 2s 3, B 5,3 = s 1 + s 2 + s 3 3s 4 and B 5,4 = s 1 + s 2 + s 3 + s 4 4s 5. The general form of B Wk for k = 2,..., n 1 is illustrated as follows B Wk = { s 1 s 2, s 1 + s 2 2s 3,..., } s j (k 1)s k. (2.10) j=1 Since B 2,i = B 3,i =... = B k,i for i = 1,..., k 1, then for simplicity we let B i = B j,i, where B i = s s i is i+1. From (2.8) ā k+1 = a k+1 proj B i a k+1. Moreover using (1.3) we have cos ϕ = cos θ k cos θ 1 + (k 1) cos θ k cos η = cos2 θ and sin ϕ = and then 1 + (k 1)(cos θ) 1 + (k 1)(cos θ) k cos θ cot ϕ = 1 + (k 1) cos θ k cos2 θ. ( Let u k+1 = a k+1 proj k sia k+1 + ) proj B i a k+1. Note that {B 1,..., B, k s i} is a set of orthogonal vectors and span B 1,..., k s i =span s 1,..., s k. Then u k+1 is orthogonal to s 1,..., s k. Now we 4
5 substitute k s i and proj B i a k+1 instead of s 1 + s 2 and proj s1 s 2 a 3 in (2.5), respectively to get v k+1. k We have v k+1 = u k+1 + cot ϕ u k+1 s i k and without loss of generality s i k v k+1 = q k+1 + cot ϕ s i k, (2.11) where q k+1 = u k+1 1 u k+1 is orthogonal to s 1,..., s k. So the (k +1)th normalized equiangular vector spanned by {a 1,..., a k } is s k+1 = v k+1 / v k+1. We can proceed this method to obtain n equiangular vectors in R n. The formula of cot ϕ can be simplified as 1+() k cos θ = k cos θ k cos 2 θ sec 2 θ+() sec θ k = k (sec θ 1)(sec θ+k). s i Algorithm 1. (Equiangular Algorithm). This algorithm produces equiangular vectors s i R n with the angle θ via the linearly independent vectors a k, for k = 1 : n 1: s 1 = a 1 / a 1 ; s = zeros(n, 1) ; r = 0 ; 2: for k = 2 : n do 3: s = s + s ; s = s/ s ; 4: q = v k ( ss T a k + r ) 1 (a k (ss T a k + r) ) ( ; k ) 5: v k = q + 1/ ((sec θ 1)(sec θ + k 1)) s ; 6: s k = v k / v ( k ; r = 0 ; ) 7: B = j=1 s j (k 1)s k ; 8: for i = 1 : k 1 do 9: r = r + B i B T i a k / B i 2 ; 10: end for 11: end for We present an algorithm named Equiangular Algorithm (EA) for this method. We call the outcome matrix of this method S = [s 1,..., s m ] of size n m, Equiangular matrix. We denote the set of all full-rank equiangular matrices of size n m with cos θ = α by EMα n m and also the corresponding nonsingular square ones by EMα n. Orthogonal vectors. The formula (2.11) reduces to the orthogonalization process if α = 0. In this case at the step k+1, we have a k+1 = a k+1 (proj B1 a k proj B a k+1 + proj Bk a k+1 ), k = 1, 2,..., n 1, (2.12) where B k = k s i. It can be said that the normalized orthogonal vector s k+1 in case of α = 0, is constructed using the set of orthogonal vectors {B 1,, B, B k } instead of the previous set of orthonormal vectors {s 1,..., s k } which is unlike the Gram-Schmidt algorithm. Below is an example for outcome of this algorithm using MATLAB. Example 2.2. For the Vandermonde matrix A = gives S 1 = and S 2 =, applying EA with θ 1 = π/3 and θ 2 = π/ , respectively. We can check that the condition number of A is κ(a) = e + 03 but κ(s 1 ) = and κ(s 2 ) = It means that the transformation of the illconditioned matrix A into S 1 and S 2 reduces the condition number. The process of producing equiangular vectors provides a matrix factorization by generating an upper triangular matrix within EA, like QR decomposition. Theorem 2.3. Suppose A R n m has full-rank (m n). For given α (0, 1) there is an equiangular matrix S EMα n m R n m and an upper triangular R R m m so that A = SR. 5
6 Proof. In step k (k < m) of EA which is applied to the column vectors of A = [a 1,..., a m ], we derive a representation of v k followed by (2.11) as v k s k = a k + d is i. So a k = d i s i + v k s k = [s 1,..., s k ] d 1. d v k = S kr k. (2.13) By EA S k can be extended as S = [ S k, S ] m k, where Sm k is a n (m k) equiangular matrix and then a k = [ S k, S ] [ ] [ ] R k Rk m k = S. (2.14) 0 0 This is the kth column vector in both sides of A = SR Example 2.4. For the Minij matrix A = , Theorem 2.3 results A = S 1 R 1 by θ 1 = π/6 where S 1 = and R = and also results A = S 2 R 2 by θ 2 = 3π/8 where S 2 = and R 2 = Corollary 2.5. Any symmetric positive definite (spd) matrix A can be factorized as a Cholesky factorization plus a rank-one matrix. Proof. Suppose A R n n is a real spd matrix, then A can be written as B T B, where B is a nonsingular matrix. SR factorization by Theorem 2.3 implies 1 α A =B T B = R T S T SR = R T R = R T [ (1 α)i n + αee T ] R α 1 =(1 α)r T R + α(r T e)(r T e) T = Cholesky + rank-one. (2.15) It is notable that this factorization is not unique with respect to a fixed angle. Actually, the number of SR factorization of a matrix A R n m of full rank, is 2 n 1. Because by altering the sign of vector q 2 in (2.1) there are two choices for constructing s 2 in step 2. Also by altering the sign of vector q 3 in (2.5) there are two choices for constructing s 3 in step 3. Therefore we have two choices in each step until the nth vector s n is obtained. References References [1] A. Björck, Numerical Methods for Least Squares Problems, SIAM, Philadelphia, PA, [2] A. Björck, Solving linear least squares problems by Gram-Schmidt orthogonalization, BIT 7, 1-21, [3] P. Casazza and J. Kovacevic, Equal-norm tight frames with erasures, Adv. Comp. Math, 18, , [4] R. A. Horn and C. R. Johnson, Matrix Analysis, Cambridge University Press, Cambridge, [5] G. R. G. Lanckriet, N. Cristianini, P. Bartlett, L. E. Ghaoui, M. I. Jordan, Learning the kernel matrix with semidefinite programming, Journal of Machine Learning Research, 5, 2772, p. 29, [6] E. Schmidt, Über die Auflösung linearer Gleichungen mit unendlich vielen Unbekannten, Rend. Circ. Math. Palermo. Ser. 1, 25, 53-77,
7 [7] E. Solivrez, E. Gagliano, Orthonormalization on the plane: a geometric approach, Mex. J. Phys, 31, no. 4, , [8] L. N. Trefethen, and D. Bau, Numerical Linear Algebra III, Philadelphia, PA: SIAM, [9] I. Todhunter, Spherical Trigonometry for the use of college and schools, with numerous examples,, fifth edition, Macmillan and Co., London, 1886; an excellent source on spherical geometry may be found at Nov 12,
Linear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationMAT Linear Algebra Collection of sample exams
MAT 342 - Linear Algebra Collection of sample exams A-x. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationOrthonormal Bases; Gram-Schmidt Process; QR-Decomposition
Orthonormal Bases; Gram-Schmidt Process; QR-Decomposition MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 205 Motivation When working with an inner product space, the most
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More information2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian
FE661 - Statistical Methods for Financial Engineering 2. Linear algebra Jitkomut Songsiri matrices and vectors linear equations range and nullspace of matrices function of vectors, gradient and Hessian
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
More informationRounding error analysis of the classical Gram-Schmidt orthogonalization process
Cerfacs Technical report TR-PA-04-77 submitted to Numerische Mathematik manuscript No. 5271 Rounding error analysis of the classical Gram-Schmidt orthogonalization process Luc Giraud 1, Julien Langou 2,
More informationLecture 6, Sci. Comp. for DPhil Students
Lecture 6, Sci. Comp. for DPhil Students Nick Trefethen, Thursday 1.11.18 Today II.3 QR factorization II.4 Computation of the QR factorization II.5 Linear least-squares Handouts Quiz 4 Householder s 4-page
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationChapter 6: Orthogonality
Chapter 6: Orthogonality (Last Updated: November 7, 7) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). A few theorems have been moved around.. Inner products
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
More informationFrame Diagonalization of Matrices
Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)
More informationOn the loss of orthogonality in the Gram-Schmidt orthogonalization process
CERFACS Technical Report No. TR/PA/03/25 Luc Giraud Julien Langou Miroslav Rozložník On the loss of orthogonality in the Gram-Schmidt orthogonalization process Abstract. In this paper we study numerical
More information2. Every linear system with the same number of equations as unknowns has a unique solution.
1. For matrices A, B, C, A + B = A + C if and only if A = B. 2. Every linear system with the same number of equations as unknowns has a unique solution. 3. Every linear system with the same number of equations
More informationMath 407: Linear Optimization
Math 407: Linear Optimization Lecture 16: The Linear Least Squares Problem II Math Dept, University of Washington February 28, 2018 Lecture 16: The Linear Least Squares Problem II (Math Dept, University
More informationMATH 167: APPLIED LINEAR ALGEBRA Least-Squares
MATH 167: APPLIED LINEAR ALGEBRA Least-Squares October 30, 2014 Least Squares We do a series of experiments, collecting data. We wish to see patterns!! We expect the output b to be a linear function of
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More informationMATH Linear Algebra
MATH 304 - Linear Algebra In the previous note we learned an important algorithm to produce orthogonal sequences of vectors called the Gramm-Schmidt orthogonalization process. Gramm-Schmidt orthogonalization
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationAssignment 1 Math 5341 Linear Algebra Review. Give complete answers to each of the following questions. Show all of your work.
Assignment 1 Math 5341 Linear Algebra Review Give complete answers to each of the following questions Show all of your work Note: You might struggle with some of these questions, either because it has
More informationMath 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Exam 2 will be held on Tuesday, April 8, 7-8pm in 117 MacMillan What will be covered The exam will cover material from the lectures
More informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
More informationLecture 4 Orthonormal vectors and QR factorization
Orthonormal vectors and QR factorization 4 1 Lecture 4 Orthonormal vectors and QR factorization EE263 Autumn 2004 orthonormal vectors Gram-Schmidt procedure, QR factorization orthogonal decomposition induced
More information1 Last time: least-squares problems
MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that
More informationSolutions to Review Problems for Chapter 6 ( ), 7.1
Solutions to Review Problems for Chapter (-, 7 The Final Exam is on Thursday, June,, : AM : AM at NESBITT Final Exam Breakdown Sections % -,7-9,- - % -9,-,7,-,-7 - % -, 7 - % Let u u and v Let x x x x,
More informationMATH 22A: LINEAR ALGEBRA Chapter 4
MATH 22A: LINEAR ALGEBRA Chapter 4 Jesús De Loera, UC Davis November 30, 2012 Orthogonality and Least Squares Approximation QUESTION: Suppose Ax = b has no solution!! Then what to do? Can we find an Approximate
More informationDot Products. K. Behrend. April 3, Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem.
Dot Products K. Behrend April 3, 008 Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem. Contents The dot product 3. Length of a vector........................
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationMATH 240 Spring, Chapter 1: Linear Equations and Matrices
MATH 240 Spring, 2006 Chapter Summaries for Kolman / Hill, Elementary Linear Algebra, 8th Ed. Sections 1.1 1.6, 2.1 2.2, 3.2 3.8, 4.3 4.5, 5.1 5.3, 5.5, 6.1 6.5, 7.1 7.2, 7.4 DEFINITIONS Chapter 1: Linear
More informationMAT2342 : Introduction to Applied Linear Algebra Mike Newman, fall Projections. introduction
MAT4 : Introduction to Applied Linear Algebra Mike Newman fall 7 9. Projections introduction One reason to consider projections is to understand approximate solutions to linear systems. A common example
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationIntroduction to Linear Algebra, Second Edition, Serge Lange
Introduction to Linear Algebra, Second Edition, Serge Lange Chapter I: Vectors R n defined. Addition and scalar multiplication in R n. Two geometric interpretations for a vector: point and displacement.
More informationWorksheet for Lecture 25 Section 6.4 Gram-Schmidt Process
Worksheet for Lecture Name: Section.4 Gram-Schmidt Process Goal For a subspace W = Span{v,..., v n }, we want to find an orthonormal basis of W. Example Let W = Span{x, x } with x = and x =. Give an orthogonal
More information(v, w) = arccos( < v, w >
MA322 Sathaye Notes on Inner Products Notes on Chapter 6 Inner product. Given a real vector space V, an inner product is defined to be a bilinear map F : V V R such that the following holds: For all v
More informationPractice Exam. 2x 1 + 4x 2 + 2x 3 = 4 x 1 + 2x 2 + 3x 3 = 1 2x 1 + 3x 2 + 4x 3 = 5
Practice Exam. Solve the linear system using an augmented matrix. State whether the solution is unique, there are no solutions or whether there are infinitely many solutions. If the solution is unique,
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationSection 6.4. The Gram Schmidt Process
Section 6.4 The Gram Schmidt Process Motivation The procedures in 6 start with an orthogonal basis {u, u,..., u m}. Find the B-coordinates of a vector x using dot products: x = m i= x u i u i u i u i Find
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationMath 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook.
Math 443 Differential Geometry Spring 2013 Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook. Endomorphisms of a Vector Space This handout discusses
More informationRecall the convention that, for us, all vectors are column vectors.
Some linear algebra Recall the convention that, for us, all vectors are column vectors. 1. Symmetric matrices Let A be a real matrix. Recall that a complex number λ is an eigenvalue of A if there exists
More informationorthogonal relations between vectors and subspaces Then we study some applications in vector spaces and linear systems, including Orthonormal Basis,
5 Orthogonality Goals: We use scalar products to find the length of a vector, the angle between 2 vectors, projections, orthogonal relations between vectors and subspaces Then we study some applications
More informationMath 321: Linear Algebra
Math 32: Linear Algebra T. Kapitula Department of Mathematics and Statistics University of New Mexico September 8, 24 Textbook: Linear Algebra,by J. Hefferon E-mail: kapitula@math.unm.edu Prof. Kapitula,
More informationMath Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT
Math Camp II Basic Linear Algebra Yiqing Xu MIT Aug 26, 2014 1 Solving Systems of Linear Equations 2 Vectors and Vector Spaces 3 Matrices 4 Least Squares Systems of Linear Equations Definition A linear
More informationLINEAR ALGEBRA SUMMARY SHEET.
LINEAR ALGEBRA SUMMARY SHEET RADON ROSBOROUGH https://intuitiveexplanationscom/linear-algebra-summary-sheet/ This document is a concise collection of many of the important theorems of linear algebra, organized
More informationMATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.
MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationa s 1.3 Matrix Multiplication. Know how to multiply two matrices and be able to write down the formula
Syllabus for Math 308, Paul Smith Book: Kolman-Hill Chapter 1. Linear Equations and Matrices 1.1 Systems of Linear Equations Definition of a linear equation and a solution to a linear equations. Meaning
More informationLecture 3: QR-Factorization
Lecture 3: QR-Factorization This lecture introduces the Gram Schmidt orthonormalization process and the associated QR-factorization of matrices It also outlines some applications of this factorization
More informationLinear Models Review
Linear Models Review Vectors in IR n will be written as ordered n-tuples which are understood to be column vectors, or n 1 matrices. A vector variable will be indicted with bold face, and the prime sign
More informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
More informationLINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS
LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,
More informationTBP MATH33A Review Sheet. November 24, 2018
TBP MATH33A Review Sheet November 24, 2018 General Transformation Matrices: Function Scaling by k Orthogonal projection onto line L Implementation If we want to scale I 2 by k, we use the following: [
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationCheat Sheet for MATH461
Cheat Sheet for MATH46 Here is the stuff you really need to remember for the exams Linear systems Ax = b Problem: We consider a linear system of m equations for n unknowns x,,x n : For a given matrix A
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More information2. Review of Linear Algebra
2. Review of Linear Algebra ECE 83, Spring 217 In this course we will represent signals as vectors and operators (e.g., filters, transforms, etc) as matrices. This lecture reviews basic concepts from linear
More informationChapter 4 Euclid Space
Chapter 4 Euclid Space Inner Product Spaces Definition.. Let V be a real vector space over IR. A real inner product on V is a real valued function on V V, denoted by (, ), which satisfies () (x, y) = (y,
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 informationAM 205: lecture 8. Last time: Cholesky factorization, QR factorization Today: how to compute the QR factorization, the Singular Value Decomposition
AM 205: lecture 8 Last time: Cholesky factorization, QR factorization Today: how to compute the QR factorization, the Singular Value Decomposition QR Factorization A matrix A R m n, m n, can be factorized
More informationInner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold:
Inner products Definition: An inner product on a real vector space V is an operation (function) that assigns to each pair of vectors ( u, v) in V a scalar u, v satisfying the following axioms: 1. u, v
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More informationMathematical Foundations: Intro
Mathematical Foundations: Intro Graphics relies on 3 basic objects: 1. Scalars 2. Vectors 3. Points Mathematically defined in terms of spaces: 1. Vector space 2. Affine space 3. Euclidean space Math required:
More informationON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH
ON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH V. FABER, J. LIESEN, AND P. TICHÝ Abstract. Numerous algorithms in numerical linear algebra are based on the reduction of a given matrix
More informationMath 416, Spring 2010 Gram-Schmidt, the QR-factorization, Orthogonal Matrices March 4, 2010 GRAM-SCHMIDT, THE QR-FACTORIZATION, ORTHOGONAL MATRICES
Math 46, Spring 00 Gram-Schmidt, the QR-factorization, Orthogonal Matrices March 4, 00 GRAM-SCHMIDT, THE QR-FACTORIZATION, ORTHOGONAL MATRICES Recap Yesterday we talked about several new, important concepts
More informationLecture 1: Review of linear algebra
Lecture 1: Review of linear algebra Linear functions and linearization Inverse matrix, least-squares and least-norm solutions Subspaces, basis, and dimension Change of basis and similarity transformations
More informationMath 261 Lecture Notes: Sections 6.1, 6.2, 6.3 and 6.4 Orthogonal Sets and Projections
Math 6 Lecture Notes: Sections 6., 6., 6. and 6. Orthogonal Sets and Projections We will not cover general inner product spaces. We will, however, focus on a particular inner product space the inner product
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationApplied Mathematics 205. Unit II: Numerical Linear Algebra. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit II: Numerical Linear Algebra Lecturer: Dr. David Knezevic Unit II: Numerical Linear Algebra Chapter II.3: QR Factorization, SVD 2 / 66 QR Factorization 3 / 66 QR Factorization
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
More informationDot product and linear least squares problems
Dot product and linear least squares problems Dot Product For vectors u,v R n we define the dot product Note that we can also write this as u v = u,,u n u v = u v + + u n v n v v n = u v + + u n v n The
More informationReview of similarity transformation and Singular Value Decomposition
Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm
More informationMath Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p
Math Bootcamp 2012 1 Review of matrix algebra 1.1 Vectors and rules of operations An p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the
More information(v, w) = arccos( < v, w >
MA322 F all206 Notes on Inner Products Notes on Chapter 6 Inner product. Given a real vector space V, an inner product is defined to be a bilinear map F : V V R such that the following holds: Commutativity:
More information6. Orthogonality and Least-Squares
Linear Algebra 6. Orthogonality and Least-Squares CSIE NCU 1 6. Orthogonality and Least-Squares 6.1 Inner product, length, and orthogonality. 2 6.2 Orthogonal sets... 8 6.3 Orthogonal projections... 13
More informationLinear Algebra. Paul Yiu. 6D: 2-planes in R 4. Department of Mathematics Florida Atlantic University. Fall 2011
Linear Algebra Paul Yiu Department of Mathematics Florida Atlantic University Fall 2011 6D: 2-planes in R 4 The angle between a vector and a plane The angle between a vector v R n and a subspace V is the
More informationLinear Algebra Highlights
Linear Algebra Highlights Chapter 1 A linear equation in n variables is of the form a 1 x 1 + a 2 x 2 + + a n x n. We can have m equations in n variables, a system of linear equations, which we want to
More information1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det
What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix
More informationMatrix Theory, Math6304 Lecture Notes from September 27, 2012 taken by Tasadduk Chowdhury
Matrix Theory, Math634 Lecture Notes from September 27, 212 taken by Tasadduk Chowdhury Last Time (9/25/12): QR factorization: any matrix A M n has a QR factorization: A = QR, whereq is unitary and R is
More informationThe Singular Value Decomposition and Least Squares Problems
The Singular Value Decomposition and Least Squares Problems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo September 27, 2009 Applications of SVD solving
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationApplied Linear Algebra
Applied Linear Algebra Gábor P. Nagy and Viktor Vígh University of Szeged Bolyai Institute Winter 2014 1 / 262 Table of contents I 1 Introduction, review Complex numbers Vectors and matrices Determinants
More informationMATH 1120 (LINEAR ALGEBRA 1), FINAL EXAM FALL 2011 SOLUTIONS TO PRACTICE VERSION
MATH (LINEAR ALGEBRA ) FINAL EXAM FALL SOLUTIONS TO PRACTICE VERSION Problem (a) For each matrix below (i) find a basis for its column space (ii) find a basis for its row space (iii) determine whether
More informationMATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003
MATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003 1. True or False (28 points, 2 each) T or F If V is a vector space
More informationMTH 2032 SemesterII
MTH 202 SemesterII 2010-11 Linear Algebra Worked Examples Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education December 28, 2011 ii Contents Table of Contents
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationLinear Algebra in Actuarial Science: Slides to the lecture
Linear Algebra in Actuarial Science: Slides to the lecture Fall Semester 2010/2011 Linear Algebra is a Tool-Box Linear Equation Systems Discretization of differential equations: solving linear equations
More informationLinear Algebra (Review) Volker Tresp 2017
Linear Algebra (Review) Volker Tresp 2017 1 Vectors k is a scalar (a number) c is a column vector. Thus in two dimensions, c = ( c1 c 2 ) (Advanced: More precisely, a vector is defined in a vector space.
More informationENGG5781 Matrix Analysis and Computations Lecture 8: QR Decomposition
ENGG5781 Matrix Analysis and Computations Lecture 8: QR Decomposition Wing-Kin (Ken) Ma 2017 2018 Term 2 Department of Electronic Engineering The Chinese University of Hong Kong Lecture 8: QR Decomposition
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationMATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)
MATH 20F: LINEAR ALGEBRA LECTURE B00 (T KEMP) Definition 01 If T (x) = Ax is a linear transformation from R n to R m then Nul (T ) = {x R n : T (x) = 0} = Nul (A) Ran (T ) = {Ax R m : x R n } = {b R m
More informationLinear Algebra, part 3 QR and SVD
Linear Algebra, part 3 QR and SVD Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2012 Going back to least squares (Section 1.4 from Strang, now also see section 5.2). We
More informationMATH 304 Linear Algebra Lecture 34: Review for Test 2.
MATH 304 Linear Algebra Lecture 34: Review for Test 2. Topics for Test 2 Linear transformations (Leon 4.1 4.3) Matrix transformations Matrix of a linear mapping Similar matrices Orthogonality (Leon 5.1
More informationELE/MCE 503 Linear Algebra Facts Fall 2018
ELE/MCE 503 Linear Algebra Facts Fall 2018 Fact N.1 A set of vectors is linearly independent if and only if none of the vectors in the set can be written as a linear combination of the others. Fact N.2
More information(v, w) = arccos( < v, w >
MA322 F all203 Notes on Inner Products Notes on Chapter 6 Inner product. Given a real vector space V, an inner product is defined to be a bilinear map F : V V R such that the following holds: For all v,
More informationBackground Mathematics (2/2) 1. David Barber
Background Mathematics (2/2) 1 David Barber University College London Modified by Samson Cheung (sccheung@ieee.org) 1 These slides accompany the book Bayesian Reasoning and Machine Learning. The book and
More information