第2章 数值计算.ppt

2.3.5 矩阵元素及重排 1．矩阵元素 在MATLAB，也可以采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依此类推 。 例如： A =[ 1 2 3 4 5 6 7 8 9]; A(3) ans = 7 A(7) ans = 3 这与我们常用的取法A(3,1)，A(1,3)的结果相同。 2．矩阵重排 当向量的元素个数能表示成m×n形式，这时我们可以将其排为矩阵形式。 格式 A=reshape(x,m,n) %将向量x重新排成m×n的二 维矩阵A 例如： x=1:12 x = 1 2 3 4 5 6 7 8 9 10 11 12 A=reshape(x,3,4) A = 1 4 7 10 2 5 8 11 3 6 9 12 注 A(:)将矩阵A每一列元素堆叠起来，成为一个列向量，产生一个12×1的矩阵，等价于reshape(A,12,1)，其转置A(:)就是行向量x。 2.3.6 矩阵复制 当矩阵的阶次在4阶以下，矩阵元素的输入可以逐个写入。但是当矩阵阶次较大且结构相同时，用逐个输入的方法太费时间，这时可以用矩阵复制的办法来加快输入的速度，常用函数repmat来进行矩阵复制。 【例2-10】 已知 A=[1 3 6 9 12]，写出矩阵B，要求复制6行向量A。 A=[1 3 6 9 12] B=repmat(A,6,1) B = 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 C=repmat(A,6,2) C = 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 1 3 6 9 12 2.3.7 稀疏矩阵 当矩阵的大部分元素是零，只有少数元素为非零元素时，这种矩阵称为稀疏矩阵。 格式 S=sparse（A） %A为原全矩阵，返回稀疏矩阵S full(S) %将稀疏矩阵转化为全矩阵函数 nnz(S) %非零个数的计数函数 spy(S) %对应的稀释矩阵图 [i,j,v]=find(A) %查找矩阵非零数所在的行数i、列数j，及其对应的非零数值v 【例2-11】 已知，试求：（1）找出矩阵A的非零元素所在的行、列及其对应的值；（2）求出其稀疏矩阵B，再转换成全矩阵C，并画出稀疏矩阵B的标示图 A=[10 0 0 0; 20 30 0 0; 0 0 0 40] [i,j,v]=find(A) i‘ = %程序显示i，是一列,转置i表示 1 2 2 3 j = %程序显示j，是一列，转置j 1 1 2 4 v = %程序显示v，是一列，转置v 10 20 30 40 B=sparse(i,j,v) %稀疏矩阵的产生