Matlab数组和向量的介绍.ppt

2019年4月29日,第1页,第3章 数组和向量,MATLAB R2010a 基础教程 清华大学出版社,教学目标 教学重点 教学内容,2019年4月29日,第2页,教学目标,掌握数组和向量的概念与性质 掌握数组与向量的操作和运算方法 了解数组与向量的实际应用,2019年4月29日,第3页,教学重点,数组的创建方法 数组的操作 向量的创建方法 向量的操作,2019年4月29日,第4页,教学内容,MATLAB数组 创建数组 数组操作 数组应用实例计算土壤体积 MATLAB向量 向量应用实例力与力矩,2019年4月29日,第5页,MATLAB数组,数组是MATLAB中的基本构件。 数组的概念 一个矢量可以只有一列，或者也可以只有一行。因此，矢量是数组的一个特殊情况。通常，数组可以有多列，并且有多行。,2019年4月29日,第6页,数组的性质(1/4),数组中的单个数据项称为元素。这些元素也同样具有取值和位置组合而成的唯一属性。在一个二维数组中，位置是指元素所在的行号与列号(按顺序)。一般而言，n维数组中元素的位置是一个包含n个索引值的向量。,2019年4月29日,第7页,数组的性质(2/4),当应用于n维数组A时，函数size()按下面两种方式中的一种返回信息： 如果按照只有一个返回值的形式调用，如sz=size(A)，将会返回一个n维向量，其中包含了数组每一维的大小。 如果按照具有多个返回值的形式调用，如rows,cols=size(A)，将返回具有所要求数量的数组A的各维大小。为了避免出错，应当给定与数组维数同的变量个数。,2019年4月29日,第8页,数组的性质(3/4),length()函数返回数组的最大维值，因此如果我们创建一个具有283大小的3维数组，size(A)将会返回2 8 3，而length(A)将返回8。 对于mn数组的转置(该操作由置于数组标识符后面的撇号字符()来表明)，将返回一nm的数组，其行与列的数值发生了交换。,2019年4月29日,第9页,数组的性质(4/4),一些值得注意的特殊情况如下： 当一个二维矩阵具有相同的行数和列数时，称其为方阵(square)。 当数组中的非零值仅出现在行号和列号相同的位置时，称该数组为对角的(diagonal)。 如果仅有一行，则称其为行向量(row vector)，或者如前一样称为向量。 如果仅有一列，则称其为列向量(column vector)，它是行向量的转置。,2019年4月29日,第10页,创建数组（1/2）,可以通过直接输入数值来创建数组，也可以通过MATLAB内置函数来创建具有某一特点的数组。 可以直接在数组中输入数值，使用“；”或者另起一行表明一行的结束。,2019年4月29日,第11页,创建数组（2/2）,函数zeros(m,n)与ones(m,n)分别产生填充值为0或1的m行n列的数组。 函数rand(m,n)与randn(m,n)分别产生取值为01之间的随机数的数组。 函数diag()具有多种形式，最常用的是diag(A)，其中A为一数组，它将A的对角线上的元素作为向量返回；diag(V)，其中V为一向量，它返回一个对角线为V的方阵。 MATLAB还提供了magic(m)函数，它产生一个填充1到m2之间数字的数组，其组织方式使得每一行、每一列以及对角线上的元素分别加起来等于相同的数。,2019年4月29日,第12页,数组操作（1/2）,获取数组中的元素 可以通过将需要获取元素的索引值用括号括起来的方式对数组中的元素进行寻址，其中第一个值为行索引值，第二个值为列索引值。 从数组中移除元素 可以使用空数组 的一种特殊用法来完成此项任务。,2019年4月29日,第13页,数组操作（2/2）,数组算术操作 数组逻辑运算 使用库函数 数组连接 数组切片 数组重排 线性化数组,2019年4月29日,第14页,应用实例 计算土壤体积（1/3）,估计为一座建筑准备地基所需要移走泥土的数量。 第一步是测量待修建筑物的地面，将得到一个矩形的方格点，每一个方格点上定义了该点的高度 。 下一步是考虑建筑地基的建筑图 。其中阴影区域是为了建造建筑地基而必须将泥土移走的地方。从建筑图中可以估计出每一个测量的方块中必须被移走泥土的比例(对于这个实例来说，它是一个介于01之间的数)。 需要被移走的泥土的总量就等于每一小方块的深度与该方块中需要移除面积的乘积之和。,2019年4月29日,第15页,应用实例 计算土壤体积（2/3）,clear clc % soil depth data for each square produced by the survey dpth = 8 8 9 8 8 8 8 8 7 8 7 7 7 7 8 8 8 7 8 8 8 8 8 8 8 7 7 7 7 7 8 7 8 8 8 7 8 8 8 8 7 7 8 7 8 8 8 8 8 7 8 8 8 8 7 7 7 8 7 8 8 8 8 8 8 8 7 6 7 7 7 7 8 8 8 8 8 8 8 8 7 7 7 7 7 6 6 7 7 8 8 7 7 8 7 7 8 7 7 7 7 7 7 7 7 7 7 8 9 8 8 9 8 7 8 7 7 7 7 7 6 7 6 7 7 8 8 8 8 9 9 8 8 8 7 6 6 6 6 7 7 8 7 8 9 8 8 7 7 7 7 7 7 6 6 7 7 7 8 8 7 8 9 8 8 7 7 7 6 7 7 6 6 8 8 8 9 9 7 8 9 9 8 8 8 8 7 7 7 7 7 8 8 9 9 9 8 8 9 8 8 7 7 8 7 7 7 7 8 8 9 9 9 8 7 8;,% estimated proportion of each square that should be excavated area = 1 1 1 1 1 1 1 1 1 1 .3 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 .7 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 .8 .4 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 .8 .3 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .7 .2 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .6 0 0 0 0 0 .7 1 1 1 1 1 1 1 1 1 1 1 .8 0 0 0 0 0 .7 1 1 1 1 1 1 1 1 1 1 1 .7 0 0 0 0 0 .4 1 1 1 1 1 1 1 1 1 1 1 .6 0 0 0 0 0 .1 .8 1 1 1 1 1 1 1 1 1 1 1 .4 0 0 0 0 0 .2 .7 1 1 1 1 1 1 1 1 1 1 .9 .1 0 0 0 0 0 0 .4 .8 .9 1 1 1 1 1 1 1 1 .6; square_volume = dpth .* area; total_soil = sum(sum(square_volume);,代码清单 计算泥土总量的脚本,2019年4月29日,第16页,应用实例 计算土壤体积（3/3）,地形测量,计算泥土体积,运行该程序，得到结果为1117.5立方单位。,2019年4月29日,第17页,MATLAB向量,向量是对相似数据项的集合进行分组的最简单方式，向量是数据的一维分组。 向量中单个的数据项通常称为元素(element)。向量元素的两个独立且截然不同的属性(数值及在向量中的位置)决定了其在某一特定向量中的唯一性。,2019年4月29日,第18页,创建向量,创建向量有两种方法： 使用一系列常值来创建向量。 通过对已有的向量进行操作产生新的向量。,2019年4月29日,第19页,向量尺寸,向量有一个专门的属性：向量的长度。在多数情况下，该长度在向量创建时就已固定。 MATLAB能够通过插入或删除元素来增加或减小向量的长度。 函数length(V)返回数组大小中的最大值 对于向量，该数值即为其长度。,2019年4月29日,第20页,向量索引,数值索引 逻辑索引 注意：命令窗口中显示的逻辑变量值是用1代表真，用0代表假。,2019年4月29日,第21页,缩短向量,注意：实际中很少将缩短向量作为对问题的正确求解方法，而且这样做可能会导致逻辑混乱。有可能的话，要使用索引来对需要保留的元素进行复制。,2019年4月29日,第22页,向量操作,MATLAB语言的本质核心在于具有丰富的数组和向量操作工具。 有3种方法直接从对标量的操作推广而来： （1）算术操作。 （2）逻辑操作。 （3）使用库函数。 有两种方法是针对一般情况下的数组和特殊情况下的向量的： （1）连接。 （2）切片(广义索引)。,2019年4月29日,第23页,应用实例力与力矩（1/2）,考虑作用于物体上P点的两个力，计算P点的合力、该合力方向上的单位向量以及该力关于M点的力矩。 可以使用MATLAB中的向量来表示该问题中的各个向量，它们均具有3个元素：向量的x、y、z方向上的值。,2019年4月29日,第24页,应用实例力与力矩（2/2）,1. clear 2. clc 3. PA = 0 1 1 4. PB = 1 1 0 5. P = 2 1 1 6. M = 4 0 1 7. % find the resultant of PA and PB 8. PC = PA + PB 9. % find the unit vector in the direction of PC 10. mag = sqrt(sum(PC.2) 11. unit_vector = PC/mag 12. % find the moment of the force PC about M 13. % this is the cross product of MP and PC 14. MP = P - M 15. moment = cross( MP, PC ),当运行该