上课材料之二:第二章 数学基础 (Mathematics)第一节 矩阵(Matrix)及其二次型(Quadratic Forms)第二节 分布函数(Distribution Function),数学期望(Expectation)及方差(Variance) 第三节 数理统计(Mathematical Statistics ) 第一节 矩阵及其二次型(Matrix and its Quadratic Forms)2.1 矩阵的基本概念与运算 一个m ×n 矩阵可表示为:v a a a a a aa a a a A mn m m n n ij ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡== 212222111211][矩阵的加法较为简单,若C=A +B ,c ij =a ij +b ij但矩阵的乘法的定义比较特殊,若A 是一个m ×n 1的矩阵,B 是一个n 1×n 的矩阵,则C =AB 是一个m ×n 的矩阵,而且∑==nk kj ikij b ac 1,一般来讲,AB ≠BA ,但如下运算是成立的:● 结合律(Associative Law ) (AB )C =A (BC ) ● 分配律(Distributive Law ) A (B +C )=AB +AC 问题:(A+B)2=A 2+2AB+B 2是否成立?向量(Vector )是一个有序的数组,既可以按行,也可以按列排列。
行向量(row ve ctor)是只有一行的向量,列向量(column vector)只有一列的向量。
如果α是一个标量,则αA =[αa ij ]。
矩阵A 的转置矩阵(transpose matrix)记为A ',是通过把A 的行向量变成相应的列向量而得到。
显然(A ')′=A ,而且(A +B )′=A '+B ',● 乘积的转置(Transpose of a production ) A B AB ''=')(,A B C ABC '''=')(。
● 可逆矩阵(inverse matrix ),如果n 级方阵(square matrix)A 和B ,满足AB=BA=I 。
则称A 、B 是可逆矩阵,显然1-=B A ,1-=A B 。
如下结果是成立的:1111111)()()()(-------='='=A B AB A A AA 。
2.2 特殊矩阵1)恒等矩阵(identity matrix)对角线上元素全为1,其余全为0,可记为I ; 2)标量矩阵(scalar matrix) 即形如αI 的矩阵,其中α是标量; 3)幂等矩阵(idempotent matrix)如果矩阵A 具有性质A A A A ==⋅2,这样的矩阵称为幂等矩阵。
定理:幂等矩阵的特征根要么是1,要么是零。
4)正定矩阵(positive definite )和负定矩阵(negative definite ),非负定矩阵(nonnegative ) 或 半正定矩阵(positive semi-definite ),非正定矩阵(nonpositive definite) 或 半负定矩阵(negative semi-definite );对于任意的非零向量x ,如有x A x '>0(<0),则称A 是正(负)定矩阵;如有xA x'≥0(≤0),非负(非正)定矩阵。
如果A 是非负定的,则记为A ≥0;如果是正定的,则记为A >0。
协方差矩阵∑是半正定矩阵,几个结论:a )恒等矩阵或单位矩阵是正定的;b )如果A 是正定的,则1-A 也是正定的;c )如果A 是正定的,B 是可逆矩阵,则AB B '是正定的;d )如果A 是一个n ×m 矩阵,且n >m ,m A r =)(,则A A '是正定的,A A '是非负定矩阵。
5)对称矩阵(symmetric matrix ); 如果A =A ′,则A 称为对称矩阵。
2.3 矩阵的迹(trace )一个n ×n 矩阵的迹被定义为它的对角线上的元素之和,记为)(A tr ,则∑==ni iiaA tr 1)(,如下结论是显然的。
1))()(A tr A tr αα= (α是标量) 特例n I tr =)(2))()(A tr A tr ='3))()()(B tr A tr B A tr +=+ 4))()(BA tr AB tr =,特例211)(ij nj n i a A A tr ∑∑==='5)循环排列原则 tr(ABCD)=tr(BCDA)=tr(CDAB)=tr(DABC) 定理:实对称矩阵A 的迹等于它的特征根之和。
因为A 是实对称矩阵,故有在矩阵C ,使得⎪⎪⎪⎭⎫⎝⎛=Λ='n AC C λλ 1,其中I C C =',所以,∑==='='=Λ=ni iA tr AI tr C C A tr AC C tr tr 1)()()()()(λ。
2.4 矩阵的秩(rank)一个矩阵A 的行秩和列秩一定相等,一个矩阵的秩就可以定义为它的行秩或列秩,记为r(A),不加证明,我们给出如下结果:1))()(A r A r '=≤min (行数、列数)2)1)()(n B r A r -+≤)(AB r ≤min ))(),((B r A r ,其中A 、B 分别为m ×n 1、n 1×n 矩阵,特例:如果A 、B 为n ×n 矩阵,而且AB=0,则)()(B r A r +≤n3))()()(A A r A A r A r '='=,其中A 是n ×n 的方阵 4))(B A r +≤)()(B r A r +5)设A 是n ×n 矩阵,且I A =2,则n I A r I A r =-++)()( 6)设A 是n ×n 矩阵,且A A =2,则n I A r A r =-+)()( 2.5 统计量的矩阵表示向量可理解为特殊的矩阵。
i 是一个其元素都为1的n 维列向量,即i'=(1,1,…,1),如果我们再假定),,,(21n x x x x=',计量经济模型中的许多统计量就可以用矩阵的形式表示出来,很方便进行数学推导。
显而易见,∑=⋅'=n i i x i x 1,∑=⋅'=n i i x x x 12,样本的均值与方差的矩阵表示如下:1)样本均值矩阵表示;事实上n i i =' 即11='i i n,而⎪⎪⎪⎪⎭⎫ ⎝⎛='111111111 i i ,x i n x n x n i i ⋅'==∑=111;2)样本方差矩阵表示易知:x i i n x i n i x i x x '=⋅'⋅⋅==⎪⎪⎪⎭⎫⎝⎛11。
其中矩阵i i n '1是一个每个元素都为n 1的n 阶方阵,从而x M x i i n I x i i n x x i x x x x x x x n 021)1()1()(∆'⋅-='-=-=⎪⎪⎪⎪⎭⎫ ⎝⎛---。
定理:矩阵0M 是幂等矩阵。
矩阵0M 的对角线上的元素为)11(n -,非对角线的元素为n1-,是一个对称矩阵。
故样本方差:)()(1)(1122x x x x nx x n S n i i -'-=-=∑=x M x nx M x n x M M x n02000111'=='⋅=。
2.6 矩阵的二次型与多元正态分布1)矩阵的二次型(Quadratic Forms )和线性变换(linear transferring ) 设P 是一数域,一个系数在数域P 中的n x x x ,,,21 的二次齐次多项式n n n x x a x x a x a x x x f 11211221112122),,,(+++= n n x x a x a 2222222+++……………………………2n nn x a + (1)称为数域P 上的一个n 元二次型,或者,在不致引起混淆时简称二次型。
例如2332223*********x x x x x x x x x +++++就是有理数域上的一个三元二次型,为了以后讨论上的方便,在(1)中,i x x j i (<)j 的系数写在ij a 2。
而不简单地写成ij a 。
和在几何中一样,在处理许多其它问题时也常常希望通过变量的线性替换简化有关的二次型,为此,我们引入定义1 设.n x x ,,1 ;n y y ,,1 是两组文字,系数在数域...........P .中的一级关系式.......⎪⎪⎩⎪⎪⎨⎧+++=+++=+++=nnn n n n nn n n y c y c y c x y c y c y c x y c y c y c x 22112222121212121111 (2) 称为由...n x x ,,1 ,n x 到.n y y ,,1 的一个线性替换,或简称线性替换,如果系数行列式.......................0≠ij c那么线性替换......(2)就称为非退化的.......。
在讨论二次型时,矩阵是一个有力的工具,因此我们先把二次型与线性替换用矩阵来表示。
令ij ji a a =, i <j由于i j j i x x x x =所以二次型(1)可以写成n n n x x a x x a x a x x x f 112112211121),,,(+++=n n x x a x a x x a 2222221221++++……………………………………22211n nn n n n n x a x x a x x a ++++∑∑===n i nj j i ijx x a11(3)把(3)的系数排成一个n ×n 矩阵⎪⎪⎪⎪⎭⎫ ⎝⎛=nn n n n n a a a a a a a a a A 212222111211 (4) 它就称为二次型(3)的矩阵,因为ji ij a a =,i ,,,,1n j =所以A A '=我们把这样的矩阵称为对称矩阵,因此,二次型的矩阵都是对称的...........。
令⎪⎪⎪⎪⎭⎫⎝⎛=n x x x X 21 于是,二次型可以用矩阵的乘积表示出来,AX X '⎪⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎪⎭⎫ ⎝⎛=n nn n n n n n x x x a a a a a a a a a x x x 2121222211121121),,,( ⎪⎪⎪⎪⎭⎫ ⎝⎛+++++++++=n nn n n n n n n n x a x a x a x x x a x a x a x a x a x x x 22112222121121211121),,,( ∑∑===ni nj j i ij x x a 11故 AX X x x x f n '=),,,(21 应该看到,二次型(1)的矩阵A 的元素ji ij a a =正是它的j i x x 项的系数的一半,因此二次型和它的矩阵是相互唯一决定的,由此还能得到,若二次型BX X AX X x x x f n '='=),,,(21且A A =',B B =',则B A =。