现代化工计算技术求解线性方程组15求解线性方程组.pptVIP

单元1 数据处理 单元1 数据处理 * 单元1 数据处理 * 任务5 求解线性方程组 求解线性方程组 线性方程组的求解方法 手工求解线性方程组 采用自编软件求解线性方程组 采用Excel求解线性方程组 * 单元1 数据处理 * 任务5 求解线性方程组 一、线性方程组的求解方法 (一)消元法(高斯Gauss消去法 ) （1）消元过程 * 单元1 数据处理 * 任务5 求解线性方程组 一、线性方程组的求解方法 (一)消元法(高斯Gauss消去法 ) （1）消元过程 （2）回代过程 最后一个方程式直接解出: x4=b4(3)/a44(3) 将x4代入上一式解出x3 : x3=(b3(2)－a34(2) x4)/a33(2) 逐次往前计算,便可求出全部xi * 单元1 数据处理 * 任务5 求解线性方程组 一、线性方程组的求解方法 (二)高斯主元消去法 ☆消元过程需用系数矩阵的对角线元素ak, k(称为主元素)做除数。若主元素为零,通过行交换找到非零ak, k,但若主元素很小,由于舍入误差及有效数字损失,常有较大误差,以致使最终解极不准确。 ☆高斯列主元消去法,即在消元之前应对方程组的首行或首列元素进行检查,并将其中绝对值最大者调整到首行或首列。 ①遍查方程组第一式（行）中的所有元素(不包括常数项),找出其中绝对值最大者作为主元素,然后将主元素及其所在列的其他元素与第一列各对应元素互换位置。这种选取主元素的方法,称为行主元法。 ②遍查方程组第一列中的所有元素,找出其中绝对值最大者作为主元素,然后将主元素及其所在行的其他元素与第一式（行）的各对应元素互换位置。这种选取主元素的方法，称为列主元法。 ③同时采用行主元法与列主元法选取主元素,称为全主元法。 返回 * 单元1 数据处理 * 任务5 求解线性方程组 二、手工求解线性方程组 例1-9 采用高斯消元法求解三阶方程组 回代:将上三角形方程组自下而上逐步进行求解,从而得出 x3=－6、x2=－1、x1=9 返回 * 单元1 数据处理 * 任务5 求解线性方程组 三、采用Excel求解线性方程组 (一) 方法一——Gauss消去法与Excel单元格驱动相结合求解 返回 回代: x3=D16/C16 =－6 * 单元1 数据处理 * 任务5 求解线性方程组 三、 采用Excel求解线性方程组 (二) 方法二——采用Excel自带函数求解线性方程组 返回 第二步:求矩阵的逆 第三步:求方程的解 第一步:求系数行列式的值,判断方程组有无解 函数:①Mdeterm②Minverse③Mmult * 单元1 数据处理 * 任务5 求解线性方程组 三、 采用Excel求解线性方程组 (三)方法三——采用Excel规划求解法求解线性方程组 第二步:选中B33,点击“工具”,选择“规划求解”,调出“规划求解参数对话框”,如图所示。选择值为（V）,输入第一个方程的常数项值“1”。可变单元格输入框中；输入：$A$33:$A$35。 第一步:依次在单元格B33,B34,B35中输入3个方程的表达式,即: =2*A33-A34+3*A35;=4*A33+2*A34+5*A35;=A33+2*A34+0*A35 其中A33、A34、A35分别代表x1、x2、x3。 * 单元1 数据处理 * 任务5 求解线性方程组 三、 采用Excel求解线性方程组 (三)方法三——采用Excel规划求解法求解线性方程组 第二步:选中B33,点击“工具”,选择“规划求解”,调出“规划求解参数对话框”,如图所示。选择值为（V）,输入第一个方程的常数项值“1”。可变单元格输入框中；输入：$A$33:$A$35。 * 单元1 数据处理 * 任务5 求解线性方程组 三、 采用Excel求解线性方程组 (三)方法三——采用Excel规划求解法求解线性方程组 第三步:设置约束条件,在“约束”窗口中点击“添加”按钮,出现“添加约束”对话框,如图1所示。在单元格引用位置输入$B$34,约束值前选择“=”,其值输入4。以同样的方法添加另一约束条件。 * 单元1 数据处理 * 任务5 求解线性方程组 三、采用Excel求解线性方程组 (三)方法三——采用Excel规划求解法求解线性方程组 返回 第四步:求解计算,添加完约束条件后,在约束窗口就多了两个表达式,即$B$34=4、$B$35=7,如图所示。按“求解”按钮,得结果对话框，此时在A33∶A35位置出现方程组的解。 * 单元1 数据处理 * 任务5 求解线性方程组 四、 采用Delphi自编软件求解线性方程组 输入方程维数“3”?输入增广矩阵?“确定”