用C语言实现矩阵的运算-.docVIP

目 录 第一章 绪 论 1 1.1 矩阵的产生 1 1.2 几种特殊的矩阵及其性质 1 1.3 矩阵的应用 3 第二章 矩阵的几种运算 5 2.1 矩阵的加法与数乘 5 2.2 矩阵的乘法 6 2.3 矩阵的转置 8 2.4 矩阵的逆 8 2.5 矩阵的特征值 10 第三章 用C语言实现矩阵运算 13 3.1 算法设计分析 13 3.1.1 矩阵乘法 13 3.1.2 矩阵的逆 13 3.1.3 矩阵的三角分解 14 3.1.4 矩阵的特征值 15 3.1.5 稀疏矩阵迭代法 16 3.2 本章小结 17 结束语 19 致 谢 21 参考文献 23 附 录 25 第一章 绪 论 矩阵理论既是学习经典数学的基础，又是一门最有实用价值的数学理论。它不仅是数学的一个重要的分支，而且也已成为现代各科技领域处理大量有限维空间形式与数量关系的强有力的工具。特别是计算机的广泛应用，为矩阵理论的应用开辟了广阔的前景。 1.1 矩阵的产生 在解决众多理论研究和工程应用问题时，将其转化为线性代数的矩阵运算问题，通常是简洁高效的。有许多实际问题和数学研究对象常常可以用一张数表表示。因此，我们建立一个数学模型来统一深入的研究这种表格。 由m×n个数(i=1,2,…,m;j=1,2,…,n)排成一个m行n列的矩阵数表： 称为m×n矩阵，通常用大写字母如A，B，…或A，B，…表示，有时也记作A=或，其中(i=1,2,…,m;j=1,2,…,n)称为矩阵的第i行第j列元或(i,j)元。元都是实数的矩阵称为实矩阵：元都是复数的矩阵称为复矩阵。当m=n时，称它为n阶方阵或n阶矩阵。 1.2 几种特殊的矩阵及其性质 前面我们知道，由m×n个数(i=1,2,…,m;j=1,2,…,n)排成一个m行n列的矩阵数表： 就是m×n矩阵，元全为0的矩阵称为零矩阵，记为0；只有一行的矩阵： A= 称为行矩阵，也称为行向量；为了避免元间的混淆，行向量也记作： A= 只有一列的矩阵： B= 称为列矩阵，也称为列向量。 如果两个矩阵的行数和列数相等，则称它们为同型矩阵。 设A=与B=是同型矩阵且对应元相等，则称A与B相等，记作A=B。即 A=B, i=1,2,…,m;j=1,2,…,n 如果A是n阶方阵，从左上角到右下角的对角线，称为A的主对角线；从右上角到左下角的对角线，称为A的次对角线。 称主对角线以下的元全为零的方阵： A= 为上三角形矩阵，称主对角线以上的元全为零的方阵： B= 为下三角形矩阵。 如果方阵的主对角线以外的元全为零，即 A= 则称它为对角矩阵，记作或。在对角矩阵中，未写出的元表示零元。 主对角线上全为1的n阶对角矩阵，即 称为n阶单位阵，记作或。 1.3 矩阵的应用 矩阵的应用非常广泛，可以说我们日常生产生活中都会应用到矩阵。 在白酒工业，成品酒的勾兑这一工序尤为重要，勾兑的目的不仅要使成品酒达到规定的酒精度，更重要的是要使成品酒中影响酒体风味的几十种主要微量成分达到预先设计好的平衡比例和具体含量。因此，使用计算机配合相关勾兑软件，就显得非常重要。如基酒及调配液的配比计算模块中，按目标含量、实际含量、调配液密度、成品密度、调配比数N类参数，依照质量守恒原理，建立起质量平衡线性方程组；且编制有行列式计算模块，可对任意阶行列式进行计算，从而改变了传统的逐次迭代算法，可直接对线性方程组进行解的存在性的判定和求解计算；小到针对单一的酒精度勾兑的配比计算，大到同时针对30种目标成分、31种基酒及调配液的调配比进行计算。 我国高速公路网络中环形结构日益增多，发展形成了复杂的网状结构，如何实现准确的收费清分就成为不得不解决的问题。由于高速公路投资费用极高 ,在实际路网中相邻节点间有两条或以上的路段直接相连的情况一般不会出现，故高速路网的结构图一般都是简单图。因此，可以定义若干矩阵，实现对高速公路网(主线、 匝道、 交叉口等) 的参数化描述。如可通过连接矩阵、属性矩阵和标识站矩阵描述高速公路的物理结构，通过规则矩阵描述收费清分规则。 第二章 矩阵的几种运算 2.1 矩阵的加法与数乘 设有两个同型的m×n矩阵A=，B=，矩阵A与B的和记作A+B，规定为： A+B= 数k与矩阵A的乘积，简称数乘，记作kA或Ak，规定为： kA=Ak= 矩阵的加法与数乘统称为矩阵的线性运算。 矩阵A=的负矩阵，记为-A，定义为 -A=(-1)A=() 矩阵A与B的减法，记作A-B。定义为 A-B=A+(-B) 例 2.1 设 A=，B= 求2A-3B。 解 2A-3B=2-3 =+ = 2.2 矩阵的乘法 矩阵乘法是出于研究线性方程组以及线性