数据结构 的案例5.docVIP

向量与矩阵 在许多科学和工程技术应用中，数组用于表示数学对象向量和矩阵。本节我们将学习如何在C中表示向量和矩阵，以及如何编写C函数实现对这些对象的操作。 向量的表示 向量是由一个数的序列（向量的分量）构成的数学对象。若向量由n个分量构成，则称向量是n维的。遗憾的是，“维”在C中的含义与向量的标准数学术语不一致。图8.8中的数组C是一维（在C的术语中）的，而它所表示的向量却是三维（在数学术语中）的。一个n一维向量在C中表示为一个长度为n的一维数组。也许你会注意到向量分量的符号标识与C中也不尽相同。向量X的第三个分量（4）在数学中记为x3，但它在C中被称为x[2],因为C数组的下标总是从0开始的。 C数组 x[0] x[1] x[2] 1 2 3 X=1,2,4 向量的 x1 x2 x3 分量 图8.8 向量在C中的表示：C数组x；向量X 标量积 分量个数相同的两个向量可以相乘，将对应分量的乘积求和，所得结果为标量（点或内）积。因此，若X=1,2,4且W=2,3,1,它们的标量积为 X·W=＝1*2+2*3+4*1=12 通常，若X=x1,x2,…,xn与W=w1,w2,…,,wn是两个n维的向量，则它们的标量积为 　　　　　　　　　　　　　 式中，xk与wk分别表示向量X和W的第k个分量。 编写计算两个n维向量X和W的程序很简单。需要一个循环计数器逐一遍历各个分量，变量sum_prod用于累加求和。计算标量积的C程序如下： sum_prod=0; for(k=0;kn;++k) sum_prod+=x[k]*w[k]; 下面我们将逐步编写的许多函数中都需要类似标量积这样的求和计算。表示求和的标记符号可以通过一个相当机械化的过程翻译为C的语句。图8.9中的求和符号如果按字面直译应当是 for(k=1;k=n;++k) 但考虑到向量的分量1至n是用数组的元素0到n-1表示的，我们必须对for循环的起点稍做改动，新的循环边界为 for(k=0;kn;++k) 　　X·W= 变量　　 for语句　　　和变量 sum_prod 　　及有关　　sum_prod 　保存和　　　变量　　　中的累加值 　　　　 　　　　　　　图8.9　求和符号与for循环参数 矩阵的表示 矩阵是由矩形数字阵列构成的一种数学对象，矩形阵列中的数称为矩阵的元素。一个m×n矩阵由m行n列构成，每行为一个n维向量，每列是一个m维向量。第i行第j列元素记为aij，在C语言中（设C中i和j的初值为0）表示为a[i][j]。故一个m×n矩阵在C语言中用一个m行n列的二维数组表示。 图8.10（a）是一个m×n（此处为4×3）矩阵，图8.10（b）是该矩阵的二维数组表示形式。 1 1 1 2 3 1 1 -1 -1 0 1 2 A= (a) (b) 图8.10 (a)　4×3矩阵 (b) 矩阵的二维数组表示 矩阵与向量的乘法 如果A是一个m×n的矩阵且X是一个n维向量，那么我们把A与X的乘积记为A*X，得到的是一个m维的向量V。一个由m行n列的矩阵只能在右边与n维向量相乘。如果A是一个m×n的矩阵且W是一个m维的向量，那么我们可以计算W和A的乘积，记为W*A，结果为一个n维的向量Z。一个m行n列的矩阵只能与一个m维的向量左乘。图8.11(a)是一个4×3矩阵A的右乘，图8.11(b) 是矩阵的左乘。下面我们对图8.11(a)进行详细讨论。 A * X = V W * A = Z *= *= 图8.11（a）中结果向量中的分量究竟是如何计算的呢？V中的第i个分量v[i] 是矩阵A的第i行与向量X的标量积。记住，在C中总是从零开始计数的。如v（l）(向量V的第二个分量)是A的行1（第二行）（可以视为一个向量）与向量X的标量积。这些向量在图8.11中以彩色标出。具体计算如下： v[1]=〈2，3，1〉*〈1，2，2=2*1+3*2+1*2=10 m*n矩阵A和n维向量X相乘计算vi 的数学通项公式是 vi=a i k* x k 将求和公式转化为C语句，并根据组下标从零开始的规