用主元连乘法定义行列式.doc

用主元连乘法定义行列式工科线性代数现代化和大众化的思路之二西安电子科技大学 陈怀琛 高淑萍email： 摘要：提出了用主元连乘积法来定义行列式，可以把高斯消元法与行列式有机地联系起来，大大简化行列式的理论讲授难度，可以避开的许多困惑的概念，并且使行列式的计算和编程十分方便快捷。这 从1995到2004的十年间，为了把科学计算软件用于机械、电子、控制等类工科课程，我写了五本书，涉及十多门课程，解了超百道的线性代数问题，真正体验了后续课程和工程问题对线性代数的需要。发现工科线性代数现代化和大众化的首要手段是引进机算。其次，许多老师根本不知道工程上是如何应用线性代数的，凭想象选材。造成在内容上，“有用的不教，教了的没用”。本文只就行列式的用法和算法问题做些探讨。一、 工科后续课究竟如何用行列式的？工科学生遇到的主要线性代数问题是解方程组。其阶数通常为五阶以上，直到成百上千。求解的基本原理是高斯消元法，采用的工具是计算机，用手工笔算是不可能的。不过现在的大学数学就是偏偏只教笔算，只管解三阶以下的题目。“不教机算”的问题我已谈得很多，2009年高教司已立项“用MATLAB和建模实践改造工科线性代数”来解决这个问题，有19所大学，45000名学生参与了试点，深受师生欢迎，现在正在继续推广，本文就不多说了。在这里我要谈的是线性代数中的一个困难的理论问题行列式，应该如何教。大家知道，“行列式不为零”是判断线性方程组解存在的根据，从数学上看，必不可少。因此现在的线性代数老师教给学生的是这样一个概念：拿来一个线性方程组，先算它的系数行列式，如果它不等于零，再去求解；否则就不算了。实际情况远非如此，我解的上百道矩阵应用题目，都是直接求解方程，没有先算行列式“判解”这个步骤。其原因何在呢？1. 先判解，后求解，这样的操作顺序要有一个前提，那就是，判解应该比求解容易。如果求解快，判解慢，那有何必多此一举去“判解”呢？因此，计算复杂度是关键，当前线性代数中不考虑计算复杂度是一个大的缺陷。2. 方程组求解现在都用高斯消元法，数学上早有证明，那是最快的方法。行列式的计算法则随定义而定，现有三种定义方法1：显式法、代数余子式法和主元连乘积法。我国的现有教材从来都不讲第三种，而前两种定义下的计算复杂度极高。三阶及以上的系统，判解远慢于求解。阶数愈高，两者的差距愈大。因此，判解实际上是不进行的。 3. 消元法求解的过程中，已经经历了解的存在性的判断。只要消元所得的行阶梯矩阵的主对角线上的元素不为零，方程组就有解；主元出现零，方程就无解。其实，求出了行阶梯矩阵，用主元连乘积的定义，已经可以很方便地通过N-1次连乘，得到行列式的值，所以求解和求行列式几乎是同时完成的，并不需要多付出双倍以至千百倍的无用功来先判解的。尤其是使用计算机解题时，如果出现了某个零主元，计算机会发出出错告警，指出系数行列式接近或等于于零。二、 用主元连乘积法定义行列式1设二阶方阵开始，设二阶方阵的系数矩阵为，只用第三类初等变换来消元，所得的行阶梯矩阵为：，对角项（主元）的连乘积D不为零，得到 D=ad-bc0这个D=ad-bc就称为矩阵A的行列式。可见，二阶线性方程组主元的连乘积就等于行列式。 如果用矩阵的下标来标注元素，其行列式为 。2.三阶方程组的行列式设三阶方程组的系数矩阵为则只用第三类初等变换的高斯消元法求得其上三角矩阵如下：要求三个对角元素的连乘积不为零。结果为： 可见用主元连乘积也同样可定义三阶系数矩阵的行列式。为了使行列式的值具有唯一性，必须限定消元变换中不使用第一类（会改变正负号）和第二类（会改变乘数）初等变换。3. 向高阶行列式的演绎由此可以推想，将N阶系数方阵用高斯消元法变换为上三角方阵，其对角线上所有主元的连乘积就是该方阵的行列式，即。我们已经知道，消元法解方程时，若主元不等于零，解就存在并可用除法求得。所以，在这个定义下，判断行列式是否为零与判断诸主元是否为零是等价的。反过来说，用消元法如能求出方程组的解，则此方程组系数矩阵的主元，因而行列式必不等于零，再用行列式去判解是多此一举。许多数学书上对此定义方法早有定论。比如在1中明确地指出主元连乘积法是现有的行列式的三种定义方法之一。不知什么原因，得不到重视。中国直到2012年游宏教授的教材2中才首次见到其推导证明。三、 高阶行列式的三种定义方法可以认为三种高阶行列式定义方法都是从二、三阶行列式从形式上向上演绎而得出的。显式法是从各元素下标的排列组合规律向上演绎，代数余子式法是从矩阵按行展开的规则进行演绎，而主元连乘积法则是从消元法所得上三角矩阵向上演绎的。1显式法：按照这个定义，nn矩阵的行列式中每一项将由不同行不同列的n个矩阵元素乘积组成(即要做n-1次乘法)，这些项应该能覆盖所有可能的排列方式，根据排列理论，行列式将有N=n!项相加。即使n=10，N将达，而需要的乘法次数为n(n-1)!之多，而每项的正负号将由这n项下标排列的逆序数决定，还需要更多的计算量。所以显式法也称为大公式法。显式法中的逆序和排列计数对非数学系的大学新生往往是拦路虎。而它的运算量不仅超越了人们笔算可能性，也超越了计算机的能力。一个25阶的行列式若按这个大公式来算，用每秒1万亿次的超级计算机，也要算1200万年才能得出结果。这种现象在计算数学上称为“维数灾难”。所以它的主要用途是数学推理，搞数学的当然不可缺少，但在应用上并没有多少价值。2代数余子式法，其思路是将nn矩阵的行列式化为n个(n-1)(n-1)较小的行列式（考虑正负号后称为代数余子式）的线性组合。逐级分解，可以减少高阶行列式的计算量；逐级综合，就可由n-1阶行列式向上定义n阶行列式。因为二阶行列式要两次乘法，按照这个方法，三阶行列式要三个二阶行列式的线性组合，即要3+3*2=9次乘法，四阶行列式要4+4*9=40次乘法，依此类推，当n很大时，要算的是n!个二阶行列式的线性组合，近似为2n!次乘法。因此，这种方法的计算量与显式法相差不大，其主要好处是可以避开逆序定义和排列组合理论，但不可能成为有实用计算价值的方法。3主元连乘法，它的核心是高斯消元法，通过等价变换消元，将系数矩阵化为上三角阵，然后把主对角线上n个主元连乘；得到行列式。这种定义的运算量已经在消元法中讨论过，实现上三角矩阵所需的计算量约为Nn3/3，n=10时，N=333次，n=25时，N5200，用现有的微机可以在微秒级的时间内完成。求行列式只要做一个n元的连乘，其运算量可忽略不计。它的另一个好处是把方程组求解和求系数行列式在同一个运算过程用同一个程序来完成，也就是把判解（的存在唯一性）和求解统一起来，实际上所有数学软件计算行列式时都采用这种方法。用MATLAB为例，只用两条语句：L,U=lu(A); % 对方阵A做LU分解D=prod(diag(U)% U的主对角线元素连乘积即为A的行列式行列式的三种定义方法所需乘法次数列表如下阶数N23451025高斯消元法求主元N3/331121413335208消元法求行列式N3/3+ N-141324453425233代数余子式法求行列式2N!29402053.1022e+025显式法求行列式 (N-1)N!212724803.7227e+026从此表可以看出，只有N=2时，用显式法判解才比消元法方便。N=3时，两者的计算量基本相同。N3时，N每加一，用显式法定义的计算量成十倍地增长。代数余子式法计算量与显式法基本相同，只有消元法的计算量最小，而且不引进其他新概念，理应作为高阶行列式计算的首选。对于非数学类的学生，应该教他们走一条比较平坦好走的路走到目标，没必要选一条难走的悬崖峭壁让大家去攀爬，因而又得去学习各种攀登的技巧和工具，人为地给课程增加了难度。四、 行列式性质教法的改变采用主元连乘法讲行列式后，逆序数、排列组合、代数余子式、矩阵按行展开、伴随矩阵等许多概念都可以不讲，那样是不是会影响学生理解行列式的性质呢？初步的探索证明，比较容易用主元连乘法定义证明的性质有下面一些：*任意三角方阵、对角方阵的行列式等于其对角元素的连乘积；*第一、二、三类初等矩阵的行列式分别为-1,k和1；*从A的某行中加上另一行的倍数，其行列式不变；*方阵中任意i,j两行交换，行列式反号；*方阵中若有一个全零行，其行列式为零。*如果A中两行的元素相同或差同样倍数，行列式等于零。*如果A奇异，则det A=0，如果A可逆，则det A0，这已在行列式定义中证明。*方阵A与它的转置AT的行列式相等。即。也有一些性质则不易用主元连乘法证明，比如行列式按行展开等。哪些性质对学生重要，关键是要研究这些性质对学生未来工作有什么应用价值。行列式按行展开的作用是便于推理，这对数学系必不可少，另一个作用是为高阶行列式的手工计算服务，其实根据前面的分析，它并不能减少多少复杂度。对于工科学生，在采用主元连乘法和计算机软件来算行列式后，这类用处不大的性质,是可以不讲的。五、 结论高等教育的现代化和大众化是现代社会发展的需要，大学数学向真实世界数学靠拢是它改革的基本目标，工科线性代数的大众化改革是很有潜力的3。把行列式与消元法无缝连接可以大大减轻课程的难度并提高它的实践性，也提高了课程的内在逻辑性。老师可以把讲课的重点放在行列式的几何及物理意义上，无需花很多时间去讲它的性质与计算，其优点是无可比拟的。至于采取了这种讲法，原来为显式法或代数余子式法准备的那些概念应该怎样处理？哪些保留？哪些扬弃？肯定是一个极有争论性的问题。因为各个专业、各类学生的要求都会不同，传统的考研命题也会继续产生影响。各个学校、专业和各类学生也都