第2章--矩阵与MATLAB.docVIP

第2章 矩阵与MATLAB 2.1 概述 MATLAB是一门语言，它的运算指令和语法基于一系列基本的矩阵运算以及它们的扩展运算。 2.2 矩阵和向量 矩阵a的大小可用如下命令获得: size(a)----size可显示出两个值:第一个值为行数(m),第二个值为列数(n)。 [m,n]=size(a)----[m,n]表示矩阵的行赋给m，列数赋给n。 2.2.1 方阵 当m=n时,a称为方阵。 2.2.2 对角阵 当aij=0,i≠j,且m=n时,得对角阵: ……(nn) 如果a=[1 3 8],那么, a=[1 3 8]； z=diag(a) 产生如下对角阵：z= 1 0 0 0 3 0 0 0 8 如果b是一个方阵，那么， diag(b)-----返回b的对角阵元素。 单位矩阵I的MATLAB表达式为： eye(n)-----其中n是方阵的大小。 c=eye(3) 产生：c= 1 0 0 0 1 0 0 0 1 2.2.3 列矩阵和行矩阵(向量) 当aij=ai1(即只有1列时),a称为列矩阵(列向量),记做: 当aij=a1j时,即矩阵只有1行,则矩阵称为行矩阵或行向量,记做: a=[a11 a12 … a1n]=[a1 a2 … a 在MATLAB中,这是向量的默认定义。 2.2.4 矩阵和向量的转置 矩阵和向量的转置用(＇)表示。 向量a的长度，即元素的个数：L=length(a) 或 L=size(a) 矩阵的大小只能用命令size获得。使用size的好处是不需要事先知道a是向量还是矩阵。但对于向量来说，length命令更方便。 2.3 向量的创建 在MATLAB中给出了几种方法用来对向量和矩阵中的元素进行赋值。第一种方法是使用冒号来指定数值范围和相邻值的步长；第二种方法是指定数值的范围和期望的数值个数。在第一种方法中，步长最为重要。在第二种方法中，数值的个数最为重要。 第一种方法： x=s:d:f 其中，s=初值，d=增量或减量值，f=结束值或终值。可产生如下的行向量x: x =[s+d,s+2d,…,s+nd] 其中，s+nd =f 。当d省略时，MATLAB默认d=1。即： x=s:f 产生向量： x =[s+1,s+2,…,s+n] 第二种方法：（指定n为从s到f的等间隔值） s=linspace(s,f,n) 其中，增量或减量值d由MATLAB通过下式计算得出： 如果没有指定n的值，MATLAB则默认n=100。 如果需要用对数等间隔，则： x=logspace(s:f:n) 产生向量： x =[10s,10s+d,10s+2d,…10f 如果没有指定n的值，MATLAB则默认n=50。 例1：若想创建[-2,1,3,5,7,9,10]的向量x。 x=[-2,1:2:9,10] 或 x=[-2,1,3,5,7,9,10] 若 z=x-1 → z=[-3,0,2,4,6,8,9] z=x；z(2)=z(2)/2；（括号里2代表某列的值） → z=[-2,0.5,3,5,7,9,10] z=x；z(2)=z(2)/2；z(3:4)=z(3:4)*3-1；（第三个和第四个数乘三减一）→ z=[-2,0.5,8,14,7,9,10] 例2：如果， y=[-1,6,15,-7,31,2,-4,-5] x=y(3:5) → x=[15,-7,31] x=[y(1),y(2),y(7),y(8)] 或 x=y([1,2,7,8]) 或 index=[1 ,2,7,8]；x=y(index) → x=[-1,6,-4,-5] 例3：如果y=[-1,6,15,-7,31,2,-4,-5]；z=[10,20,30,40,50,60,70,80] 用sort函数把向量y按升序排列，然后根据向量y中的新元素顺序重新排列向量z。 函数sort的调用格式为： [ynew,index]=sort(y) 其中，ynew是重新排列后的向量；index是一个由元素位置所构成的向量，它表示排列后向量y中元素的原始位置号。 y=[-1,6,15,-7,31,2,-4,-5]；z=[10,20,30,40,50,60,70,80]； [ynew,index]=sort(y),znew=z(index) 执行结果为： ynew → [-7 -5 -4 -1 2 6 15 31] index → [4 8 7 1 6 2 3 5] znew