第四讲 线性方程组与最小二乘法.pptVIP

《测量程序设计》 第四讲 线性方程组与最小二乘法 成都理工大学地球科学学院测绘工程系测绘工程专业 主讲教师： 杨容浩 联系电话： 电子邮箱： yangronghao@cdut.cn * 一、线性方程组及其求解方法概述 二、一般线性方程组求解 第四讲 线性方程组与最小二乘法 * 4.1 线性方程组及其求解方法概述 线性方程组的一般形式： （1）当m=n时，若A非奇异，则线性方程组有唯 一解， 否则，有无穷多组解； （2）当mn时，线性方程组有无穷多组解； （3）当mn时，若A的秩tr(A)=n，则线性方程组有唯 一 解；若tr(A)n，有无穷多组解；若tr(A)n，方程 不存在精确解，可以求得近似解。 * 4.1 线性方程组及其求解方法概述 （1）当tr(A)=n时，A非奇异 ① 直接法：将结构上比较复杂的原始方程组，通过等价变换化成结构简单的方程组，使之变成易于求解的形式，然后再通过求解结构简单的方程组来得到原始方程组的解。通过有限次的算术运算，若计算过程中没有舍入误差，可以求出精确解的方法。 （G通常是对角矩阵、三角矩阵或者一些简单结构的矩阵。） ② 间接法：把原方程组改写成一个等价方程组，对这个等价方程组建立迭代格式。是用某种极限过程去逐次逼近方程组解的方法。 * 4.1 线性方程组及其求解方法概述 （2）当tr(A)n时 无法求解精确解，只能求解近似解。 最小二乘解：指当方程组 无法求得精确解， 而以欧式距离下最接近的解作为近似，所求得的解。 最小二乘的正规方程： 为最小二乘解。 * 4.2 一般线性方程组求解 （1）高斯消元法 对于线性方程组 ，写成分量形式为 若记 则有 * 4.2 一般线性方程组求解 （1）高斯消元法 第一次消元：设 由第2,3，…，n个方程减去 第一个方程乘以 则有 其中， 消除第一个未知数x1后，与原方程组同解的方程组 * 4.2 一般线性方程组求解 （1）高斯消元法 第二次消元：设 由第3,4，…，n个方程减去 第二个方程乘以 则有 消除第二个未知数x2后，与原方程组同解的方程组 其中， * 4.2 一般线性方程组求解 （1）高斯消元法 继续前述过程，经过n-1次消元后，得到与原方程组同解的方程组 其中， * 4.2 一般线性方程组求解 （1）高斯消元法 通过回代过程，可以求得每一个未知数的值： 高斯消元法总的计算量约为 高斯消元法始终假设主元 若 则消元 过程无法进行；实际上，即使 但 很小 时，用它作除数对计算结果也是很不利的。 * 4.2 一般线性方程组求解 （2）全选主元高斯消元法 对于如下所示的方程组，有性质： 任意交换线性方程组两行系数和常量，线性方程 组解不变； 任意交换线性方程组两列，并同时交换与列号对 应的未知数的位置，线性方程组解不变； 线性方程组一行系数加上另一行系数的m倍（m 为非零实数），线性方程组解不变。 * 4.2 一般线性方程组求解 （2）全选主元高斯消元法 输 入：方程组系数矩阵A（1:n,1:n），常量数组b[1:n]。 输 出：方程组系数上三角矩阵A（1:n,1:n）， 列交换记录 数组T（1:n）。 过程名：CPGElimination(A,b) begin T（1:n）← 1:n; for k=1:n-1 begin 确定满足k≤μ≤n和k≤λ≤n，使得 |A(μ,λ)|=max{A(i,j)|i=k:n, j=k:n} A(k,1:n) ? A(μ,1:n) ; A(1:n,k) ? A(1:n,λ) ; T(k) ? T(λ); if A(k,k)≠0 begin for i=k+1:n begin m=A(i,k)/A(k,k); for j=k+1:n A(i,j)=A(i,j)-m*A(k,j);