第四章空间行列式运算与空间权重矩阵设定.doc

空间计量经济学导论（詹姆斯.勒沙杰）课件 范 巧 fanqmn@ 重庆科技学院经济系 小范经济工作室 在经济学的边缘上 拟讲授的主要内容 ? 空间行列式的作用：多维空间的平面转换 ? 空间行列式的基本运算法则 ? 空间权重矩阵和空间最近邻矩阵的设定 I ? ?W ? 空间行列式 的设定与转换 n ( ) I ? ?W ? 1 ? 逆矩阵 的基本计算 n ( ) I ? ?W ?1 ? 空间模型参数效应计算中 的矩阵转化 n ln I ? ?W ? 对数行列式 的基本计算 n ln I ? ?W ? 对数行列式 的蒙特卡罗模拟 n ln I ? ?W ? 对数行列式 的复杂运算 n ln I ? ?W ? 对数行列式 的误差对参数估计的影响 n ? 空间计量经济模型参数估计的闭式解 1.1 多维空间的平面转换：转换矩阵的利用 i?i ? 转换矩阵：设主对角线元素为1，其余元素为C的 矩阵，为i阶转换 矩阵。 ? ? 1 C C ? ? ? C 1 C ? ? ? ? T ,其中，i 2, 3, ,n ? ? ? ? ? ? ? ? ? ? C C 1 ? ? ? i?i ? 转换矩阵的作用：将多维的空间体系转换成为平行的n面体。 ? 转换矩阵的运用案例：正方形的平面拉伸 ?0 1 0 1? 1 C ?0 1 0 1? S ? ? ??????????????T S ? ? ? ? ? ? ? ? ? 左乘2阶转换矩阵 [ ] 0 0 1 1 C 1 0 0 1 1 1.1 多维空间的平面转换：转换矩阵的利用(续) ? 转换矩阵的运用案例（续）； ? ? ? ? ? ? ? 0 1 0.9 1.9 0 1 0.5 1.5 T S S; T S ; T S ? ? ??? ? ? ? ? ??? ? ? ? ? 0 0.9 1 1.9 0 0.5 1 1.5 ? ? ? ? C?0 C?0.9 C?0.5 ? C的不同取值对平面拉伸影响的图示分析： (1.9,1.9) 行列式S代表的面积为1，则行列 式转换导致TS的面积取决于|T|。 (0.5,1) (0,1) (0.9,1) (1, 1) (1, 0.9) (1.5,1.5) T S ? (1?C2 )S ? S 转换矩阵的引入将导致四边形 面积缩小。 (1, 0.5) (0, 0) (1, 0) 1.2 转换矩阵的病态取值 ? 转换矩阵病态：转换矩阵的行列式取值为0时导致转换失败的空间维 度变换过程。 ? 转换矩阵病态的案例 ： Y ? X ? ?? 左乘 ? 的转换矩阵 例1. , T 0 ? ? 时也能完全拟合，模型失效 =0 ? ?T 例2. , , Y ? X ? ?? 其中，Y ? u v w X ?? 3 ?1 1 1? ? ? ? ?T TY ? u ?v ? w u ?v ? w u ?v ? w 令 ? ? ? ? ，则 T 1 1 1 0 ? ? 1 1 1 ? ? ? ? ? ? ? 时，能完全拟合，模型失效 u v w 1.3 平面转换的收缩效应与分布调整 ? 平面转换的收缩效应：如PPT1.1正方形的拉伸面积缩小，更一般如书 上P58. ? 平面转换收缩效应条件下保持单位体积的分布调整： 2 1 1 ?1 2 N u ? ? ?n ? u ? ? T ?? u ? ? 二元正态分布： ( ) (2 ) exp( ( ) ( )) 2 n n 1 2 (u ) (u ) 2 ? ? ? ? ? ? ?? ? ? n?? T 1 ? ? ? 1 2 ? ? ? ? ? ? T分布：T(u) ( ) ( ) (n ) (1 ) 2 2 n 2.1 矩阵运算的几种主要算法 ? 转换运算的主要算法： 高斯消元法：将某一行的倍数加到另一行得到三角矩阵； LU法：基于单位下三角矩阵L、上三角矩阵U、置换矩阵P，结合 A ? L U P PA=LU，得到 ； Cholesky分解： 2 A ? RT R ? A ? R RT ? R 对于对称正定矩阵， ； 对于非对称矩阵A，则： n n A abs(U ) A A abs A abs U ? ? ? ? ? ln T 2ln ( ) ? 2 ln ( ) ii ii i 1 1 ? i? 2.1 矩阵运算的几种主要算法(续) ? 转换运算的主要算法（续）： 频谱分解法： ?0 1? A I aB,其中，B ,求解 A ? ? ? ? ? ? 2 1 0 ? ? ? 构造矩阵B的频谱分解等式： B ?V? V B T ?? ? ?? ? ? 1 0 0.7071 0.7071 ? ? ?? ?? ?? ? ??? ?? ?? ? ??? ? ? B ; V ; VVT In; V VT 1 0 1 0.70