Ch线性代数应用.ppt

⑷相应的MatLab函数 三角分解 三角分解 正定阵的 分解 矩阵或向量的 范数 矩阵或向量的 范数 矩阵或向量的 范数 矩阵或向量的 条件数 矩阵或向量的 条件数 矩阵或向量的 条件数 对前面的例子运行下面程序, 结果为 和前面计算结果一致. 例 对 阶 矩阵及 矩阵, 计算相应的条 件数. 程序为: 结果为: 进一步地有: 结果为 误差分析 设方程 给 以微小扰动 由此引起的 的 扰动记为 即有 从而有 再由范数性质得 又 即有 两式做商得: 由此可见: 对 的敏感程度依赖于 回到前例, 方程为 此时 2.迭代法 对大型的稀疏矩阵 求解矩阵方程 通常使 用迭代的方法加以求解. 以下面的例子来说明迭代过程. 例 求解方程 解 将方程改写成 上式表明, 对每一组 代入上 式的右端, 即可得到另一组 若这两组 相等, 即得到 方程之解. 程序如下: 结果为: 该方程的精确解为 计算结果表明, 经 过 次迭代, 所得到的向量与精确解相当接近. ⑴雅克比迭代 设 其中 若对角阵 非奇异, 则方程 转化为 记: 则方程的迭代形式 可写成 由此可见, 若序列 收敛于 则 必是方程之解. 该迭代即称为雅克比迭代. ⑵高斯—赛德尔迭代 在迭代法的引例中, 可以发现迭代可以进一步修改为: 即 若对角阵 非奇异, 则迭代式转化为 记 则 该迭代即称为高斯—赛德尔迭代. 程序为: 下三角阵 下三角阵 两种迭代法结果比较: 雅克比迭代 赛德尔迭代 ⑶收敛性 两种迭代公式的一般形式为 设 是方程的解, 即 两式相减后有 将后式代入前式则有 因此 设 的特征值为 记 称其为矩阵 的谱半径, 由此得到: 结论1 若矩阵 严格对角占优, 即 则雅克比迭代与高斯—赛德尔迭代均收敛. 结论2 若矩阵 为对称正定阵, 则高斯—赛德尔迭代 收敛. 在上例中, 矩阵 此时 为严格对角占优的. 因而雅克比迭代收敛. 六、应 用 1.投入和产出模型 国民经济各个部门之间存在着相互依存的关系, 每个 部门在运转中将其它部门的产品或半成品（称为投入） 经过加工变为自己的产品（称为产出）, 如何根据各部 门间的投入产出关系, 确定各部门的产出水平, 以满足 社会需要, 是投入产出分析中的课题. 投入产出表描述了国民经济各部门之间的生产和消耗, 投入和产出的数量关系. 下表是一张简化的中国2002年投入产出表. 表中国民经济由农业、工业、建筑业、运输邮电业、 批零餐饮业和其他服务业6个中间部门构成, 数字表示它 们之间的投入产出关系、外部需求、初始投入等. 5414 2341 1570 2875 16814 2918 总投 3093 1218 703 659 4815 1663 初始 5414 2725 542 219 225 272 1285 146 服务 2341 927 273 130 43 140 749 79 批零 1570 477 146 67 163 128 527 62 运输 2875 2691 124 23 20 3 9 5 建筑 16814 4083 1223 557 403 1444 8605 499 工业 2918 1284 13 127 13 229 788 464 农业 产出 外部 服务 批零 运输 建筑 工业 农业 投入 产出 投入产出表 数据表的说明: 行表示该部门的产品分配给各部门作生产使用的数量. 列表示该部门生产中消耗的各部门产品的数量. 消耗系数 定义 表示第 个部门1个单位的产出对第 个部门的 直接消耗. 消耗系数由投入产出表可直接计算. 例如对 此时总产出 而工业部门 对其投入为 因而消耗系数为 同理可以计算其它系数. 由此得消耗系数表. 0.100 0.094 0.143 0.095 0.076 0.050 服务 0.050 0.056 0.027 0.049 0.045 0.027 批零 0.027 0.029 0.104 0.045 0.031 0.021 运输 0.