计算机对化工的影响PPT,化工计算机数据与图形处理PPT-3.ppt

化工计算机数据与图形处理PPT-3.ppt

(55页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

19.90 积分

第三章 EXCEL解方程组 本章主要内容3.1 Excel矩阵计算3.2 解线性联立方程3.2.1 行列式法3.2.2 矩阵法3.3 Newton-Raphson迭代法 解非线性方程组3.4 方程组的规划求解 3.1 Excel 矩阵运算1. Excel的矩阵格式矩阵是数据(数值和变量)按行和列排成矩形的阵列： Excel工作表中用数组表示矩阵，单元格的位置对应矩阵元素的位置，但标记顺序相反: 矩阵元素aij是先行(i)后列(j)，而单元格A3是先列(A)后行(3)。3.1 Excel 矩阵运算1. Excel的矩阵格式 列矩阵或行矩阵在Excel工作表中放在同一列或同一行，称为一维数组。 m×n矩阵元素依次放在工作表的一个区域中，称为二维数组。3.1 Excel 矩阵运算2. 数组加(减)法 进行数组加(减)法运算时，一般要求数组中单元格的个数相同。示例：建立二维数组A1:D3，一维数组F1:I1。一维数组加(减)法： 1)A1:D1一行数据与F1:I1另一行数组相加，得到的仍是含四个单元格的一行数组，因此必须选定加和结果数组的位置，本例选F3:I3。3.1 Excel 矩阵运算2. 数组加(减)法2)键入“＝”，用鼠标选定A1:D1区域，被选区域出现闪动的虚线框，同时编辑框的等号后出现A1:D1。3)键入“＋”，选中F1:I1区域。同样被选区域出现闪动虚线框，加号后出现F1:I1。4)按Ctrl＋Shift＋Enter键。F3:I3区域出现加和结果，编辑框的公式被大括号括住:{=A1:D1+F1:I1}。3.1 Excel 矩阵运算注意：1)若开始未选定区域F3:I3，而只选中了F3单元格，则完成上述各步计算后在F3出现数值2。2)数组公式是一个整体，无论是编辑、移动、清除都必须选取公式包括对所有单元格。3)大括号不用输入，按Ctrl＋Shift＋Enter后自动产生，表明这是数组公式。如果输入公式时自己加上，Excel将认为输入的是文字。4)若相加的两个数组单元格数目不等，则出现错误信号#N/A。3.1 Excel 矩阵运算注意：5)若加一常数到一数组，Excel会自动将该常数加到数组的每个单元格中。 二维数组加法类似于一维数组，此外还有如下特点：二维数组可以与行单元格或列单元格数目相同的一维数组相加。3.1 Excel 矩阵运算示例结果：3.1 Excel 矩阵运算3. 数组乘以一常数 数组乘以一常数，得到与原数组大小相同的新数组，新数组的每个单元格的值等于原数组单元格的值乘以此常数。3.1 Excel 矩阵运算4. 矩阵乘法 矩阵A(m×n)与B(n×p)相乘，得一矩阵C(m×p)。仅当第一个矩阵的列数等于第二个矩阵的行数时，两矩阵才能相乘。 矩阵相乘时，各矩阵元素按一定规则相乘再相加。Excel有矩阵相乘的工作表函数MMULT(array1，array2)。3.1 Excel 矩阵运算4. 矩阵乘法计算规则：其中i为行数，j为列数。例1. MMULT({1,3;7,2},{2,0;0,2})例2. MMULT({3,0;2,0}, {2,0;0,2}) 例3. MMULT({1,3,0;7,2,0;1,0,0}, {2,0;0,2}) 3.1 Excel 矩阵运算注意：1)输入MMULT函数前要根据乘积矩阵大小(m×p)选定它的单元格区域。2)函数名前有等号，输入完毕同时按Ctrl＋Shift＋Enter。3)数组乘法不等于矩阵乘法，它是两个数组对应单元格各自相乘，数组相乘用*，公式写作＝数组1*数组2。3.1 Excel 矩阵运算5. 逆矩阵 行数等于列数的矩阵称为方阵。若一方阵的主对角元素aii均为1，其他位置上的元素均为0，则该矩阵称为单位矩阵I 。 设有一n阶方阵A和另一n阶方阵B，它们的乘积为：AB=BA=I 则B为A的逆阵，可记为A-1。3.1 Excel 矩阵运算5. 逆矩阵 Excel求逆矩阵的工作表函数为 MINVERSE(array)例1. MINVERSE({4,-1;2,0})例2. MINVERSE({1,2,1;3,4,-1;0,2,0})例3. MINVERSE()为一矩阵，可参与其他矩阵运算。3.1 Excel 矩阵运算6. 行列式行列式是在研究线性方程组的解法中产生的概念。n阶行列式可表示为： 3.1 Excel 矩阵运算6. 行列式 行列式通常用D表示，其元素用竖直线而非圆括号括住。 方阵与行列式是两个不同的概念：n阶方阵是n2个数按一定方式排列成的数据表；n阶行列式是按一定的运算规则得到的一个确定值。3.1 Excel 矩阵运算6. 行列式Excel计算行列式值的函数： MDETERM(array)例1 MDETERM({1,3,8,5;1,3,6,1;1,1,1,0;7,3,10,2})3.1 Excel 矩阵运算7. 转置矩阵 矩阵A的行和列的元素交换后得到它的转置矩阵A’：3.1 Excel 矩阵运算7. 转置矩阵1)Excel转置矩阵函数： TRANSPOSE(数组)注：要先根据原矩阵选定转置矩阵区域，且数组函数输入同时按Ctrl＋Shift＋Enter键。2)选择性粘贴－转置注：粘贴时只需选定区域的左上单元格即可。3.1 Excel 矩阵运算7. 转置矩阵3.2 解线性联立方程线性联立方程： 设一组含n个未知数的n个方程： 其中xi(x1, x2, x3,…) 是实验未知数，ci是实验测量值，aij是系数。这些方程必须是线性独立的，上式的矩阵表示为： AX=C3.2 解线性联立方程A是系数矩阵：X和C分别为解矩阵和常数矩阵：3.2.1 行列式法根据Cramer法则，线性联立方程有唯一解的条件是其系数行列式为非零值：此时方程的解为：3.2.1 行列式法 是方程组的系数行列式，但其中的第i列元素被常数列阵C所取代。例如 是用常数C列阵取代系数矩阵的第一列元素所得的行列式：3.2.1 行列式法示例：Excel步骤：1)输入注释文字和系数数值。2)在E2输入: =MDETERM(A2:C4)，得到系数行列式的值－9，说明联立方程有唯一解。3.2.1 行列式法3)复制联立方程的系数矩阵A2:C4，粘贴到下边A6:C8、A10:C12和A14:C16。再复制常数列数组D2:D4，将其分别粘贴到A6:A8、B10:B12、C14:C16。得到分子行列式。4)分别解出三个分子行列式的值，再根据公式解出x，y，z的值。3.2.1 行列式法3.2.2 矩阵法线性方程的矩阵表示经数学变换得： X＝A-1C即解矩阵(X)等于系数矩阵的逆矩阵(A-1)乘于常数矩阵(。省略部分。on-Raphson迭代法原理： 对方程组的每一个方程在其近似解 处用Teller级数展开，只取线性项得：3.3 Newton-Raphson迭代法原理： 其中 是方程 在 处的一阶偏导数。上式为△xi的线性方程组：3.3 Newton-Raphson迭代法原理： 将初始值 代入函数fi和偏导数 ，然后用解线性方程组的方法计算得到 ，从而得到自变量xi的第一次迭代值：继续迭代，可得：3.3 Newton-Raphson迭代法原理： 直至： 为一小正数。 即为满足指定精度的原方程的解。3.3 Newton-Raphson迭代法示例： 方程组： 每个方程对每一自变量的偏导数为：3.3 Newton-Raphson迭代法示例： 于是有线性方程组(1)：3.3 Newton-Raphson迭代法Excel解方程组步骤： 1)A列标记迭代次数。在A2、A3、A4输入x1:0、 x2:0、 x3:0后，选定A2:A4区域，拖拽填充柄到适当位置。冒号后的数字表示迭代次数。2)B2、B3、B4分别输入x1、x2、x3的初始值0.5。3)C2:E4单元格区域输入线性方程组(1)的系数矩阵。Δxi的系数是含xi的数学表达式。3.3 Newton-Raphson迭代法Excel解方程组步骤： 4)F2:F4区域输入常数列矩阵。根据式(1)，F2输入=-(B2^2+B3^2+B4^2-1)/2，F3单元格输入=-(2*B2^2+B3^2-4*B4)/2，F4输入 =-(3*B2^2-4*B3+B4^2)/2。5)用矩阵法解线性方程式(1)：选定单元格G2:G4，输入公式：=MMULT(MINVERSE(C2:E4),F2:F4)，同时按Ctrl+Shift+Enter，得到第一次迭代改变量Δxi1 。3.3 Newton-Raphson迭代法Excel解方程组步骤： 6)用矩阵加法在B5:B7填入第一次迭代后的自变量xi1 ：选定B5:B7，输入公式=B2:B4+G2:G4，按Ctrl＋Shift＋Enter。7)重复解线性方程组过程：选取C2:G4，将填充柄向下拖3行。在G4:G7得到第二次迭代改变量Δxi1 。8)重复迭代：选B5:G7，向下拖填充柄直至两次迭代改变量小于10－15。3.3 Newton-Raphson迭代法3.4 规划求解规划求解解方程组是单变量求解解一元方程的延伸。目标单元格是方程组中任一方程所在，目标值是该方程的常数。可变单元格是方程组的自变量所在。约束为方程组中其余方程。“规划求解”改变可变单元格的值，使得目标值与约束均得到满足。3.4 规划求解3.4.1 规划求解解线性方程组 以前面的线性方程组为例： x，y，z是可变单元格变量，将F1的常数作为目标单元格的数值，F2和F3的常数作为约束条件。任意设置可变单元格的值，目标是寻找满足约束条件(6,3)，并使目标值为0的x，y，z。3.4 规划求解Excel步骤：1)B2:D2输入未知数初始值(任意数)1,1,1。2)根据F1式，在B4单元格输入公式：＝2*B2+C2-D2；根据F2，在C4输入：＝B2-C2+D2；根据F3，在D4输入：=B2+2*C2+D2。3)打开“工具”菜单，单击“规划求解”，出现“规划求解参数”对话框。 3.4 规划求解Excel步骤：4)“设置目标单元格”输入方程F1所在位置：$B$4，“等于”栏设置目标值为0。“可变单元格”编辑框输入未知数所在单元格：$B$2:$D$2。3.4 规划求解Excel步骤：5)单击“约束”选项区的“添加”按钮，出现添加约束对话框，在“单元格引用位置”处输入方程F2的位置：$C$4，用下拉箭头选取等号“＝”，在“约束值”处输入F2的常数6。单击“添加”，给出新的添加约束对话框，继续增“单元格引用位置”处输入F3的位置：$D$4，选“＝”，约束值为3。单击“确定”，回到“规划求解参数”对话框。 3.4 规划求解3.4 规划求解6)约束条件将出现在“约束”列表区：$C$4=6和$D$4=3。若约束条件设置有误，可通过“更改”或“删除”对约束条件进行改变。7)所有参数设置好后，按“求解”，得到相应结果：3.4 规划求解3.4.2 解非线性方程组 多元一次线性方程比较容易，对初始值设定无严格要求。若方程为非线性，则必须设置不同初始值，才能保证解的可信性。例：在Cd(ClO4)2和NaCN的水溶液中有如下平衡：3.4 规划求解 根据上述三个解离平衡，由物料平衡和电荷平衡可得到含三个未知数[CN]、[Cd2+]和[H+]的三个方程： 3.4 规划求解 CL是加入的NaCN总浓度，CM是加入的Cd(ClO4)2总浓度，KW是水的离子积(＝10-14)。 三个方程均为未知浓度的高次函数，用“规划求解”步骤如下：1)在B1:B3单元格输入Cd2＋总浓度CM ，CN－的总浓度CL，HCN的酸解离常数Ka。2)D2:F2是要求解的未知浓度[CN]、[Cd2+]和[H+]所在位置。3.4 规划求解3)三个方程均有求和项，先将其分别算出。B6:B9是各级积累稳定常数β i。在C6、D6、E6分别输入公式：=B6*$D$2^A6、=A6*C6、=(2-A6)*C6。用自动填充得各项的值，然后用自动求和，在C10、D10、E10分别得到F1、F2、F3方程中的求和项。4)在D4:F4单元格，根据三个方程用D2:F2单元格的值计算，分别输入公式： ＝E2+E2*C10-B1 =D2+D2*F2/B3+E2*D10-B2 =F2+B2-1.0E-14/F2-2*B1-D2+2*E2+E2*E103.4 规划求解5)解非线性方程组，初始值的设置很重要。由于Cd2＋与CN－形成的配合物很稳定，加上HCN的解离常数很小，可以设[CN-]的初始值：0.01-3*0.001＝0.007。 Cd2＋几乎全部与CN－配位，游离的很少，可以设初始[Cd2＋]＝10－12。HCN酸性很弱，CN-在水溶液中显碱性，初始[H＋]＝10－106)以D4为目标，E4和F4为约束，用“规划求解”改变D2:F2单元格的值，使目标值和约束值均为0。按“求解”后，显示找不到有用的解。3.4 规划求解7)改变设置，初始[Cd2＋]＝10－13，仍然找不到有用的解。8)再改变设置，初始[H＋]＝10－11，“规划求解结果”显示可找到一解，满足所有的约束及最优状态。9)结果很好，F1、F2、F3接近于0。3.4 规划求解 关 键 词： 数据 化工 ppt 处理 图形 计算机

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。