将(1)中的第一个方程除以系数a1(它叫做“主导”元素)并令 (2) 则得一个新的方程 (3) 再由方程(3)及(1)中的后面三个方程消去x,这样便得到了一个辅助方程组,它包括具 有三个末知数的三个方程,此种消去法易于施行,只须顺次将方程(3)乘以an,11(也 就是乘以第二、第三和第四行的“主导”元素),再由(1)中的对应方程减去此式即可,消去 一个末知数以后所得的新方程组其系数用a代表 an1b1(i,j≥2 其次将新方程组中的第一式除以它的“主导”元素a,则得方程 (5) 其中 然后仿照前面的方法继续进行,我们便得到了一组具有两个未知数的两个方程,它们的系 数呈如次的形式 ain-an,byn(i,i≥3) 将这组方程的第一式除以主导元素a3,并令 4;2 b3i (j>3) 则得方程 x3 T burrs b3s, 2 最后再做一步,即可得出一个方程,它只含一个未知数而其系数为a4s将这个方程 除以a43则得 x4= b45 将具有系数b-1(G>i的一切方程合并,便得到一个三角形的方程组它与原有的方 程组等价;它的解就是原有方程组的解,我们要注意,上述方法只有当所有的“主导”元紊 都不等于零时才能使用 我们把求三角形方程组的系数的手续称为正面过程,而把求三角形方程组的解的手 续称为反面过程(参看附表) 我们还要讲一下验算的方法,用代换x一x+1,则我们得到一组以x为变数的方 程组,它的系数与原来的方程相同,而它的常数项等于原方程的系数与常数项之和,我们 可以同时计算这两个方程组,求出解彩并视其是否等于x1+1,这就是验算方法 现在简单地说明一下附表 正面过程是用如下方法来施行的,写出矩阵系数(包括常数项与核验和)将第一行除 以主导元素,并将结果写成矩阵最末一行,再求出第一个辅助系数ain(i≥2):从已知 矩阵任取一个元素,由它减去一个乘积—一就是上述元素所在的那一行的主导元素与上 述元素所在的那一列的最末元素的乘积,重复施行这种手续,当我们得出了仅含一行的矩 阵时,正面过程便完成了
附表1 a1,0010.4210.540.60.3 0421.00 32 0.44 0540.321 220.7 + 0660.44 3.22 0.54 0.66 0.3 82300.093200.16280.374001.45360 a31, t33,1 a34,,,tale, 093200.70840-0.136400.538001,20320 0.70201.29280 0.113160.197670.454101.76193 069785-0,154820.495681.03871 0221850.70301.4881 04977073804|1.23591 1.482402,48240 x3 1.039172.03916 1.25780-0.257 在反面过程中,我们利用包含1的各行而由最末一行开始,精确地说,在这些行的最 后一行里,我们从常数项的一列中得到了最后一个未知量的值,而在核验列中得到了核验 值然后可以逐次得出各个未知量的值只要由倒数第二列的元素减去对应系数b与前面 所得未知量值的乘积即可,在表格的末尾写出1字,可以帮助我们找出在所要各行中对 应于已知x的系数,例如 0.45410—0.I1316×1.03917-0.19767×1.48240=0.04348 最后,我们还要指出用这种方法解n个变数的线性方程组所需的乘法与除法的运算 次数为n(n2+3n-1) §3.消去法的几何解释 先看两个变数的情况 l by I: ax+by=c 在平面上各表示一条直线两个直线有一交点:消去y得出仅有x的方程,这是这个交 点在x轴上的投影 也可以这样看:第一、二方程各表示一条直线l与.由方程 a
,定义出一族直线,这些直线由x+n表示,这些直线有一个重要性质,就是通过l与 的交点,不难证明:反过来,凡是通过4与l的交点的直线也在这族之中,在这族直线 中有一条平行于y轴的,这条直线便是消去y后的方程 再看三个变数的情况 ax十by+cx=d, x+by+cs-d ax+by+cz=d 这表示三个平面,平面族 a ul 代表通过1与以的交线的所有的平面。由!与 中消去x而得出的方程可以看成为:它代表通过 x交线而平行于x轴的平面,也可以看成为:这条 交线在(y,z)平面上的投影,就是y,x平面上的 条直线,再从l,ψ中消去x,又得y,z平面上 图1 的一条直线 因此,消去r,可以看作把三维空间的三个平面求交点的问题变为在y,r平面上求 两条直线的交点的问题这两条直线,正是两条空间直线(平面的交线)的投影 般讲来:一个 十…+ 可以看成为n维空间的超乎面.消去法便是把n从空间m个超平面求交点的问题化为 n-1维空间n-1个超乎面求交点的问题 §4.消去法的力学解释 在一条两端固定的弦线上取n点P1,……,Pn,在这n点各加一重物,也就是在这些点 各有一向下的力F1,……,F我们来研究这些点的垂度y,…,y, 我们假定弦线上的力适合于“线性叠加原 则 1 1°,两组力叠加,其对应的垂度也相加 2°.所有力都乘以同一实数,则所有的垂 度也乘上这一个相同的数 以表示当在P点上作用一个单位力时 点P的垂度,这样,力F1,…,F。的联合作用 后的垂度yx……y等于 图2 2 解线性方程组的问题,也就是给了垂度y,…,y,要求出力F,…,Fn的问题了 在P1点加一个反作用力R,这样单位力作用于P时,P点的垂度等于 十R
考虑把弦线固定于P的情况,也就是 R 也就是在P点加一个单位力,如果要在P加一个力使P固定,这个力是-a/4au,这 时 从(1)式消去F1得 Ca: i-ayiai/an) F,- y; -an1/4m. 这就是加支点后的平衡方程在P加了支点,在P作用一个单位力,b就是P的垂度 逐步消去,就是逐步加支点的过程 55经济平衡 假定有n种生产品P1,…,P。生产一个单位P需要a单位P,如果各产品的数量 是x1,…,x,为了生产这些产品,P类产品的总消耗是 ∑a;,1≤i≤n 能够供给市场的数量是 a计x 因此知道了市场需要b1…,b反回来考虑给各工业的生产指标x1, 也是一个 解线性方程的问题 这类方程当然可以用消去法解,但更好是用叠代法解,关于叠代法将来再谈 §6.线性迴归分析 某一变量ξ决定于n个因素 我们已经做了N次实验得出的实验数据是 我们考虑线性关系 问题是怎样的线性关系,差方和最小,也就是,如果使 求怎样的a使
最小即求 1) 的极小值 ank, b k 我们现在证明 k 的解答x〓;使(1)取最小值 我们现在来证明这一点:如果a1;……,a并不适合于(2)。例如;有一个廴使 1;一b 我们考虑 F(4……·ak-1,a+s,a+1,…,an) ∑a,mP+εn 2 a2)的+8∑( am)+28(b 4:)+e2a 8"d 凑方得 F(4,…,ak-1:+,;…,an F(a,a1,…,an)+a数(6+ 如果a≠0,则F(a1…,a)不是最小值,因为在(3)式中取E akL aRKE FCa 的数值小于F(a1,…,a)的数值了 因此:求迴归平面的问题一变而为解线性方程的问题了 致于要证明,适合于(2)的解一定使F取最小值,这一点的证明不难,如果(2)仅有 个解,当然亳无问题,因为由(3)可知不适合(2)的都不可能使F极小(读者自证 (2)一定有解,并处理(2)有不止一个解的情况) 方程组(2)当然可以用消去法来解但是这是一个有对称系数的联立方程式即 关于这样的方程组我们另有较好的计算方法 以上的证明的优点之一,也许有人会指出,它避开了微积分,直接用初等的“凑方”法 来处理了,实际上,更好的优点在于这个方法介绍了计算数学上的一个重要方法—松弛 法