(*^ ::[ Information = "This is a Mathematica Notebook file. It contains ASCII text, and can be transferred by email, ftp, or other text-file transfer utility. It should be read or edited using a copy of Mathematica or MathReader. If you received this as email, use your mail application or copy/paste to save everything from the line containing (*^ down to the line containing ^*) into a plain text file. On some systems you may have to give the file a name ending with ".ma" to allow Mathematica to recognize it as a Notebook. The line below identifies what version of Mathematica created this file, but it can be opened using any other version as well."; FrontEndVersion = "NeXT Mathematica Notebook Front End Version 2.2"; NeXTStandardFontEncoding; fontset = title, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e8, 24, "Times"; ; fontset = subtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e6, 18, "Times"; ; fontset = subsubtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, italic, L1, e6, 14, "Times"; ; fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, bold, L1, a20, 18, "Times"; ; fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, bold, L1, a15, 14, "Times"; ; fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, L1, a12, 12, "Times"; ; fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 10, "Times"; ; fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, L1, 12, "Courier"; ; fontset = output, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier"; ; fontset = message, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, L1, 12, "Courier"; ; fontset = name, inactive, noPageBreakInGroup, nohscroll, preserveAspect, M7, italic, B65535, L1, 10, "Times"; ; fontset = header, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, L1, 12, "Times"; ; fontset = leftheader, 12; fontset = footer, inactive, nohscroll, noKeepOnOnePage, preserveAspect, center, M7, italic, L1, 12, "Times"; ; fontset = leftfooter, 12; fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12, "Courier"; ; fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; paletteColors = 128; currentKernel; ] :[font = title; inactive; Cclosed; preserveAspect; startGroup] Subspaces :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Review :[font = text; inactive; preserveAspect] Let us recall the definition of an (abstract) vector space: Definition: A VECTOR SPACE (v.s.) over a field of scalars is a set V of elements called vectors, together with two binary operations, scalar multiplication * (usually suppressed) and vector addition + , subject to the following vector space laws: For all u, v, w Î V and scalars a and b, (1) (Closure of vector additon) u + v Î V. (2) (Commutativity of addition) u + v = v + u. (3) (Associativity of addition) u + (v + w) = (u + v) + w. (4) (Additivive Identity) There exists an element 0 Î V such that u + 0 = u = 0 + u. (5) (Additive Inverse) There exists, for each u Î V, an element -u Î in V such that u + (-u) = 0 = (-u) + u. (6) (Closure of scalar multiplication) a*v Î V. (7) (Distributive Law) a*(u+v) = a*u + a*v. (8) (Distributive Law) (a + b)*u = a*u + b*u. (9) (Associative Law) (a b) * u = a*(b*u). (10) (Monoidal Law) 1*u = u. NOTATION: u + (-v) = u - v. ALSO: we often refer to the ªvector space Vº, without identifying explicitly the binary operations. They are taken for granted. SIMPLE FACTS: (deducible from the vector space properties alone): For all vectors u, 0*u = 0, where the latter 0 means the zero vector. For all vectors u, (-1)*u = -u. ;[s] 13:0,0;331,1;332,2;403,3;404,4;575,5;577,6;659,7;660,8;678,9;679,10;771,11;772,12;1299,-1; 13:1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = text; inactive; preserveAspect; endGroup] EXAMPLES OF VECTOR SPACES: (1) The spaces Rn and Cn, n >=1, over the fields R and C, respectively. Also, C is a vector space over R. (2) The set M(2,2) of all 2x2 real matrices, together with the usual matrix addition and scalar multiplication. (3) The set V = C[0,1] = {f(x) | f(x) is a continuous function on the interval [0,1] }, together with the usual function addition and scalar multiplication: (f+g)(x) = f(x) + g(x) (a*f)(x) = a f(x). This is a very non-standard vector space. It seems kind of odd, but we are really thinking of functions f(x), 0 <= x <= 1, as ªvectorsº. Of course, one has to check that all the laws are satisfied. This is routine: they are. What plays the part of the zero vector is the zero constant function. Now, what is a subspace of a given vector space V? DEFINITION: A subspace of the vector space V is a subset W of V such that W, together with the binary operations it inherits from V , forms a vector space (over the same field of scalars) in its own right. EXAMPLES OF SUBSPACES: (1) Let V = R3 and W = {[x,y,x] | x, y Î R}. (2) Let V = M(2,2) and W = {{{a,b},{0,a}} | a , b Î R } (3) Let V = C[0,1] and W = P_2 = {f(x) = a + b x + c x2 | a, b, c in R} In every case we have to verify that laws (1) - (10) are valid for each subset W of the v.s. V. This is unpleasant!! Do we have to do all this work?? No, and here is why: SUBSPACE TEST: Let W be a nonempty subset of the v.s. V. Then W is a subspace of V if and only if (1) For all u, v Î W, u + v Î W (this is really closure of addition on elements of W). (2) For all u Î W and scalars a, a*u Î W (this is really closure of scalar multiplication on elements of W). It isn't too hard to see what the Subspace Test works: certainly, any subspace W of V must satisfy these conditions since it is a vector space in its own right with the operations inherited from V. Conversely, if the nonempty subset W of V satisfies these conditions, then vector space properties (1) and (6) follow. Properties (2)-(3) and (7)-(10) follow simply because elements of W are elements of V and they hold for elements of V. The only properties to be checked are (4) and (5). Since W is nonempty, there is an element v Î W. Now one of the Simple Facts above is that (-1)v = -v, the additive inverse of v as an element of V. By condition (2) of the Subspace Test, (-1)v Î W. Therefore -v Î W. By condition (1) of the Subspace Test, v + (-v) = 0 (the zero vector of the vector space V) also belongs to W. It follows that W contains an additive identity, namely, the very same one as V has, and that every element of W has an additive inverse in W. This disposes of vector space laws (4) and (5) and shows that W is a vector space in its own right, hence a subspace of V. APPLICATION: A very common occurrence is a linear combination (l.c.) of vectors v1, v2, ... , vn, by which we mean an expression of the form, for suitable scalars a1, a2, ...an, a1 v1 + a2 v2 + ... + an vn. The set of all possible linear combinations of the vectors v1, v2, ... ,vn has a name: the linear span of v1, v2, ..., vn. DEFINITION: The linear span of vectors v1, v2, ...., vn is the set of all vectors span{v1, v2, ..., vn} = { a1 v1 + a2 v2 + ... + an vn | a1, a2, ..., an are scalars}. Note that the span depends on the field of scalars in use. Here is the application we referred to: FACT: If v1, v2, ..., vn are vector in the v.s. V then lin{v1, v2, ..., vn} is a subspace of V. This is a fact that you should be able to prove. It is also in your text. ANOTHER IDEA: if V is a vector space such that V = lin{v1, v2, ..., vn}, then the vectors v1, v2, ..., vn are said to SPAN the vector space V. EXAMPLE: {1,0,0}, {0,1,0}, and {0,0,1} span the vector space R3. As a matter of fact, so does the set of vectors {1,0,0}, {0, 1,0},{0,0,1}, {1,1,2}. Why do you think that one would prefer the first spanning set? DEFINITION: A set of vectors v1, v2, ... , vn is a BASIS of the vector space V if this set spans V and no subset of v1, v2, ... vn spans V. :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Subspaces and Matrices :[font = text; inactive; preserveAspect] Suppose we are given a matrix A. For the sake of definiteness, we will work with the following matrix A: :[font = input; preserveAspect] A={{1,-2,0,0,3},{2,-5,-3,-2,6},{0,5,15,10,0},{2,6,18,8,6}} :[font = text; inactive; preserveAspect] Here are three very important vector spaces that are associated with the matrix A. Let us say that A = [C1, C2, ..., Cn] where the Ai's are the columns of A, and that A = [R1, R2, ..., Rm]T where the Ri's are the rows of A. Then: The column space of A: C(A) = lin{C1, C2, ..., Cn} The row space of A: R(A) = lin{R1, R2, ..., Rn} The null space of A: N(A) = { v Î Rn ½ A.v = 0} (Note: replace R by C in the case of complex scalars in all that follows. We have to use one symbol, so we default to R.) R(A) and C(A) are definitely vector spaces. As a matter of fact, R(A) is a subspace of Rm and C(A) is a subspace of Rn. (How come?) EXAMPLE: Describe the row space, column space and null space of the simple matrix of the following cell in simple geometrical terms. Where do these subspaces live? For the null space, by the way, it might help to click on the cells after the definition of ;[s] 41:0,0;111,1;112,2;115,3;116,4;124,5;125,6;139,7;140,8;183,9;184,10;187,11;188,12;196,13;197,14;198,15;199,16;212,17;213,18;284,19;285,20;288,21;289,22;297,23;298,24;336,25;337,26;340,27;341,28;349,29;350,30;388,31;389,32;391,33;392,34;394,35;395,36;622,37;623,38;653,39;654,40;928,-1; 41:1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,32,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,64,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,32,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,32,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,32,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = input; preserveAspect] B = {1, 3, 1} :[font = input; preserveAspect] v = {x,y,z} :[font = input; preserveAspect] B.v :[font = text; inactive; preserveAspect] Now let's take a look at a harder example, namely, the matrix A we defined above. How would you define the row space and column space of A? Write this out on a separate piece of paper. It might help to use the following comman :[font = input; preserveAspect] MatrixForm[A] :[font = input; preserveAspect] NullSpace[A] :[font = text; inactive; preserveAspect] Now for another idea about spanning sets. Remember the examples of the last section (look them up). There was something undesirable about one spanning set: it had unnecessarily many vectors in the list. A spanning set which does not have too many vectors is called a basis of the vector space it spans. Thus: DEFINITION: A basis for the vector space V is a minimal spanning set, that is, a spanning set which has the property that no smaller set will span the space. Give an example of a basis for R3. Now for a key question: how can we spot redundancies in the list of rows or columns of a given matrix like our A? I claim the answer lies in the following cell: ;[s] 5:0,0;333,1;338,2;514,3;515,4;682,-1; 5:1,11,8,Times,0,12,0,0,0;1,10,8,Times,2,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,11,8,Times,32,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = input; preserveAspect; endGroup] RowReduce[A] :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Exercises :[font = text; inactive; preserveAspect] Find the row space, the column space, and the null space of the following matrices. Are there any relationships between the dimensions of (the number of independent vectors in) these spaces ? :[font = input; preserveAspect] æ1 0 2 7 5ö 1) A = ÷ 3 3 -1 -3 0÷ è2 1 -2 0 5ø ;[s] 12:0,0;15,1;16,2;36,3;37,4;49,5;51,6;70,7;71,8;87,9;88,10;108,11;109,-1; 12:1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0; :[font = input; preserveAspect] æ3 2 10 5ö 2) B = ç1 1 0 2÷ ç4 -2 8 0÷ è-3 5 -2 5ø ;[s] 16:0,0;15,1;16,2;31,3;32,4;44,5;45,6;61,7;62,8;78,9;79,10;94,11;95,12;111,13;112,14;127,15;129,-1; 16:1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0; :[font = input; preserveAspect; endGroup; endGroup] æ 2 2 1 0ö ç 1 0 2 1÷ 3) C = ç 3 -1 -1 0÷ ç -2 0 0 2÷ è 4 -4 2 1ø ;[s] 20:0,0;14,1;15,2;33,3;34,4;48,5;49,6;67,7;68,8;78,9;79,10;97,11;98,12;112,13;113,14;131,15;132,16;146,17;147,18;165,19;166,-1; 20:1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,9,Helvetica,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0; ^*)