数组和向量学习教案.pptx

《数组和向量学习教案.pptx》由会员分享，可在线阅读，更多相关《数组和向量学习教案.pptx（24页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1数组和向量数组和向量(xingling)第一页，共24页。2023年2月7日第 2 页教学(jio xu)重点n n数组的创建(chungjin)方法n n数组的操作n n向量的创建(chungjin)方法n n向量的操作第2页/共24页第二页，共24页。2023年2月7日第 3 页教学内容n nMATLAB数组n n创建数组n n数组操作(cozu)n nMATLAB向量第3页/共24页第三页，共24页。2023年2月7日第 4 页3.1 MATLAB数组数组n nMATLAB的一个重要功能是进行向量(xingling)和矩阵的运算，而向量(xingling)和矩阵主要由数组来表示。

2、n n数组是MATLAB中的基本构件。n n数组的概念：（略）n n一个矢量可以只有一列，或者也可以只有一行。因此，矢量是数组的一个特殊情况。通常，数组可以有多列和多行。第4页/共24页第四页，共24页。2023年2月7日第 5 页数组的性质数组的性质(xngzh)(1/4)n n数组中的单个数据项称为元素。这些元素具有取值和位置组合而成的唯一属性。在一个二维数组中，位置是指元素所在的行号与列号(按顺序(shnx)。一般而言，n维数组中元素的位置是一个包含n个索引值的向量。第5页/共24页第五页，共24页。2023年2月7日第 6 页数组的性质数组的性质(xngzh)(2/4)n n当应用(y

3、ngyng)于n维数组A时，函数size()按下面两种方式中的一种返回信息：n n如果按照只有一个返回值的形式调用，如sz=size(A)，将会返回一个向量（n维），其中包含了数组A每一维的大小。n n如果按照具有多个返回值的形式调用，如rows,cols=size(A)，将返回具有所要求数量的数组A的各维大小。n n为了避免出错，给定返回变量的个数应当与数组维数相同。第6页/共24页第六页，共24页。2023年2月7日第 7 页数组的性质数组的性质(xngzh)(3/4)n nlength()函数返回数组的最大维值，因此如果我们创建一个(y)具有283大小的3维数组，size(A)将会返回2

4、 8 3，而length(A)将返回8。n n对于mn数组(实数)的转置（该操作由置于数组标识符(变量)后面的撇号字符()来表明），将返回一个(y)nm的数组，其行与列发生了交换。第7页/共24页第七页，共24页。2023年2月7日第 8 页数组的性质数组的性质(xngzh)(4/4)一些值得注意的特殊情况一些值得注意的特殊情况(qngkung)(qngkung)如下：如下：当一个二维矩阵具有相同的行数和列数时，称其为方当一个二维矩阵具有相同的行数和列数时，称其为方阵阵(square)(square)。当数组中的非零值仅出现在行号和列号相同的位置时，当数组中的非零值仅出现在行号和列号相同的位置

5、时，称该数组为对角的称该数组为对角的(diagonal)(diagonal)。如果仅有一行，则称其为行向量如果仅有一行，则称其为行向量(row vector)(row vector)，或者如，或者如前一样称为向量。如果仅有一列，则称其为列向量前一样称为向量。如果仅有一列，则称其为列向量(column vector)(column vector)，它是行向量的转置。，它是行向量的转置。如果一个数组如果一个数组(矩阵矩阵)没有任何元素，则称为空数组没有任何元素，则称为空数组(empty array)(empty array)(空矩阵空矩阵)。第8页/共24页第八页，共24页。2023年2月7日第

6、9 页3.2 创建创建(chungjin)数数组（组（1/2）可以通过直接输入数值来创建数组，也可以通过MATLAB内置函 数 来 创 建 具 有 某 一 特 点(tdin)的数组。直接在数组中输入数值创建数组时，使用“；”或者另起一行表明一行的结束。第9页/共24页第九页，共24页。2023年2月7日第 10 页创建创建(chungjin)数组数组（2/2）n n通过通过(tnggu)(tnggu)函数创建数组函数创建数组n n函数函数zeros(m,n)zeros(m,n)与与ones(m,n)ones(m,n)分别产生填充值全部为分别产生填充值全部为0 0或或1 1的的mm行行n n列的

7、数组。列的数组。n n函数函数rand(m,n)rand(m,n)产生取值为产生取值为0101之间的均匀分布随机之间的均匀分布随机数的数组，数的数组，randn(m,n)randn(m,n)产生均值为产生均值为0 0的正态分布随机的正态分布随机数的数组。数的数组。n n函数函数diag()diag()具有多种形式。具有多种形式。diag(A)diag(A)，其中，其中A A为一数组，为一数组，它返回一个向量，将它返回一个向量，将A A的对角线上的元素作为向量的对角线上的元素作为向量元素；元素；diag(V)diag(V)，其中，其中V V为一向量，它返回一个对角为一向量，它返回一个对角方阵，对

8、角线元素为方阵，对角线元素为V V的元素。的元素。n nMATLABMATLAB还提供了还提供了magic(m)magic(m)函数，它产生一个填充函数，它产生一个填充1 1到到m2m2之间数字的数组，其组织方式使得每一行、每之间数字的数组，其组织方式使得每一行、每一列以及对角线上的元素分别加起来等于相同的数。一列以及对角线上的元素分别加起来等于相同的数。第10页/共24页第十页，共24页。2023年2月7日第 11 页3.3 数组操作数组操作(cozu)（1/2）n n获取数组中的元素 n n可以通过将需要获取元素的索引值用括号括起来的方式对数组中的元素进行寻址，其中第一个值为行索引值，第二

9、个值为列索引值。n n也可以在数组中写入(存储)数值，如A(2,3)=0，若位置超出范围，则自动(zdng)扩展数值。n n从数组中移除元素n n可以使用空数组 的一种特殊用法来完成此项任务。但数组元素必须整行或整列移除。如A(3,:)=或A(:,2)=等。第11页/共24页第十一页，共24页。2023年2月7日第 12 页数组操作数组操作(cozu)（2/2）n n数组其它数组其它(qt)(qt)操作操作n n数组算术操作数组算术操作(P/42-43)(P/42-43)n n数组逻辑运算数组逻辑运算(P/43-44)(P/43-44)n n使用库函数使用库函数(P/44)(P/44)n n数

10、组连接数组连接(P/45)(P/45)n n数组切片数组切片(P/45)(P/45)n n数组重排数组重排(P/45-46)(P/45-46)n n线性化数组线性化数组(P/46-49)(P/46-49)第12页/共24页第十二页，共24页。2023年2月7日第 13 页3.4 应用应用(yngyng)实例实例计算土壤体积计算土壤体积(1/3)n n估计为一座建筑准备地基所需要移走泥土的数量。估计为一座建筑准备地基所需要移走泥土的数量。n n第一步是测量待修建筑物的地面，将得到一个矩第一步是测量待修建筑物的地面，将得到一个矩形的方格点，每一个方格点上定义了该点的高度形的方格点，每一个方格点上定

11、义了该点的高度 。n n下一步是考虑下一步是考虑(k(k ol ol)建筑地基的建筑图建筑地基的建筑图 。其中。其中阴影区域是为了建造建筑地基而必须将泥土移走阴影区域是为了建造建筑地基而必须将泥土移走的地方。从建筑图中可以估计出每一个测量的方的地方。从建筑图中可以估计出每一个测量的方块中必须被移走泥土的比例块中必须被移走泥土的比例(对于这个实例来说，对于这个实例来说，它是一个介于它是一个介于0101之间的数之间的数)。n n需要被移走的泥土的总量就等于每一小方块的深需要被移走的泥土的总量就等于每一小方块的深度与该方块中需要移除面积的乘积之和。度与该方块中需要移除面积的乘积之和。第13页/共24

12、页第十三页，共24页。2023年2月7日第 14 页应用应用(yngyng)实例实例 计算土壤体计算土壤体积（积（2/3）nnclearclearn nclcclcn n%soil depth data for each square produced by the survey%soil depth data for each square produced by the surveyn ndpth=dpth=n n8 8 9 8 8 8 8 8 7 8 7 7 7 7 8 8 8 78 8 9 8 8 8 8 8 7 8 7 7 7 7 8 8 8 7n n8 8 8 8 8 8 8 7

13、7 7 7 7 8 7 8 8 8 78 8 8 8 8 8 8 7 7 7 7 7 8 7 8 8 8 7n n8 8 8 8 7 7 8 7 8 8 8 8 8 7 8 8 8 88 8 8 8 7 7 8 7 8 8 8 8 8 7 8 8 8 8n n7 7 7 8 7 8 8 8 8 8 8 8 7 6 7 7 7 77 7 7 8 7 8 8 8 8 8 8 8 7 6 7 7 7 7nn8 8 8 8 8 8 8 8 7 7 7 7 7 6 6 7 7 88 8 8 8 8 8 8 8 7 7 7 7 7 6 6 7 7 8n n8 7 7 8 7 7 8 7 7 7 7 7 7

14、 7 7 7 7 88 7 7 8 7 7 8 7 7 7 7 7 7 7 7 7 7 8nn9 8 8 9 8 7 8 7 7 7 7 7 6 7 6 7 7 89 8 8 9 8 7 8 7 7 7 7 7 6 7 6 7 7 8n n8 8 8 9 9 8 8 8 7 6 6 6 6 7 7 8 7 88 8 8 9 9 8 8 8 7 6 6 6 6 7 7 8 7 8n n9 8 8 7 7 7 7 7 7 6 6 7 7 7 8 8 7 89 8 8 7 7 7 7 7 7 6 6 7 7 7 8 8 7 8n n9 8 8 7 7 7 6 7 7 6 6 8 8 8 9 9 7

15、89 8 8 7 7 7 6 7 7 6 6 8 8 8 9 9 7 8n n9 9 8 8 8 8 7 7 7 7 7 8 8 9 9 9 8 89 9 8 8 8 8 7 7 7 7 7 8 8 9 9 9 8 8nn9 8 8 7 7 8 7 7 7 7 8 8 9 9 9 8 7 8;9 8 8 7 7 8 7 7 7 7 8 8 9 9 9 8 7 8;l%estimated proportion of each square that should be excavatedlarea=l1 1 1 1 1 1 1 1 1 1 .3 0 0 0 0 0 0 0l1 1 1 1 1 1

16、 1 1 1 1 .7 0 0 0 0 0 0 0l1 1 1 1 1 1 1 1 1 1 1 .8 .4 0 0 0 0 0l1 1 1 1 1 1 1 1 1 1 1 1 1 .8 .3 0 0 0l1 1 1 1 1 1 1 1 1 1 1 1 1 1 .7 .2 0 0l1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .6 0 0l0 0 0 .7 1 1 1 1 1 1 1 1 1 1 1 .8 0 0l0 0 0 .7 1 1 1 1 1 1 1 1 1 1 1 .7 0 0l0 0 0 .4 1 1 1 1 1 1 1 1 1 1 1 .6 0 0l0 0 0 .1

17、.8 1 1 1 1 1 1 1 1 1 1 1.4 0l0 0 0 0 .2.7 1 1 1 1 1 1 1 1 1 1.9 .1l0 0 0 0 0 0 .4 .8 .9 1 1 1 1 1 1 1 1.6;lsquare_volume=dpth.*area;ltotal_soil=sum(sum(square_volume);代码代码(di m)清单清单 计算泥土总量的脚本计算泥土总量的脚本第14页/共24页第十四页，共24页。2023年2月7日第 15 页应用应用(yngyng)实例实例 计算土壤计算土壤体积（体积（3/3）地形地形(dxng)测量测量 计算泥土计算泥土(nt)体积体积

18、 运行该程序，得到结果为1117.5立方单位 第15页/共24页第十五页，共24页。2023年2月7日第 16 页3.5 MATLAB向量向量(xingling)n n向量是对相似数据项的集合进行分组的最简单方式，向量是数据的一维分组。n n向量中单个的数据项通常称为元素(element).向量元素的两个独立且截然不同的属性(数值及在向量中的位置(wi zhi)决定了其在某一特定向量中的唯一性。第16页/共24页第十六页，共24页。2023年2月7日第 17 页创建创建(chungjin)向量向量n n创建向量有两种方法：n n使用一系列常值来创建向量。n n直接输入。n n使用冒号输入某个范

19、围内的数值。n n使用linspace()函数创建；如linspace(0,10,7)，第一个参数为下限，第二个参数为上限，第三个参数为向量中的数值个数。n n使用zeros(1,n)、ones(1,n)、rand(1,n)、randn(1,n)等函数来创建。n n通过对已有的向量进行操作(cozu)产生新的向量。第17页/共24页第十七页，共24页。2023年2月7日第 18 页向量向量(xingling)尺寸尺寸n n向量有一个专门的属性：向量的长度。在多数情况(qngkung)下，该长度在向量创建时就已固定。n nMATLAB能够通过插入或删除元素来增加或减小向量的长度。n n函数len

20、gth(V)返回数组大小中的最大值n n 对于向量，该数值即为其长度。第18页/共24页第十八页，共24页。2023年2月7日第 19 页向量向量(xingling)索引索引n n数值(shz)索引 n nn n逻辑索引 n n如由Mask=true false false true得n n mask=1 0 0 1n n（注意：命令窗口中显示的逻辑变量值是用1代表真，用0代表假）n n.第19页/共24页第十九页，共24页。2023年2月7日第 20 页缩短缩短(sudun)向量向量n n从向量中移除元素可以缩短(sudun)向量。可以使用空向量 的一种特殊用法来完成此项任务。n n 注意：

21、实际中很少将缩短(sudun)向量作为对问题的正确求解方法，而且这样做可能会导致逻辑混乱。有可能的话，要使用索引来对需要保留的元素进行复制。第20页/共24页第二十页，共24页。2023年2月7日第 21 页向量向量(xingling)操作操作n nMATLABMATLAB语言的本质核心在于具有丰富的数组和向语言的本质核心在于具有丰富的数组和向量操作工具。量操作工具。n n有有3 3种方法直接从对标量的操作推广而来：种方法直接从对标量的操作推广而来：n n（1 1）算术操作。）算术操作。n n（2 2）逻辑操作。）逻辑操作。n n（3 3）应用库函数。）应用库函数。n n有两种方法是针对一般有

22、两种方法是针对一般(ybn)(ybn)情况下的数组和特殊情况下的数组和特殊情况下的向量的：情况下的向量的：n n（1 1）连接。）连接。n n（2 2）切片）切片(广义索引广义索引)。第21页/共24页第二十一页，共24页。2023年2月7日第 22 页3.6 应用应用(yngyng)实例实例力与力矩（力与力矩（1/2）n n考虑作用于物体上考虑作用于物体上P P点的两个力，点的两个力，计算计算P P点的合力、该合力方向上点的合力、该合力方向上的单位向量以及该力关于的单位向量以及该力关于MM点的点的力矩。力矩。n n可以使用可以使用MATLABMATLAB中的向量来表中的向量来表示该问题中的各

23、个向量，它们示该问题中的各个向量，它们(t men)(t men)均具有均具有3 3个元素：向量个元素：向量的的x x、y y、z z方向上的值。方向上的值。第22页/共24页第二十二页，共24页。2023年2月7日第 23 页应用实例应用实例(shl)力力与力矩（与力矩（2/2）n n1.clear1.clearnn2.clc2.clcnn3.PA=0 1 13.PA=0 1 1nn4.PB=1 1 04.PB=1 1 0n n5.P=2 1 15.P=2 1 1nn6.M=4 0 16.M=4 0 1nn7.%find the resultant of PA and PB7.%find t

24、he resultant of PA and PBn n8.PC=PA+PB8.PC=PA+PBn n9.%find the unit vector in the direction of PC9.%find the unit vector in the direction of PCn n10.mag=sqrt(sum(PC.2)10.mag=sqrt(sum(PC.2)nn11.unit_vector=PC/mag11.unit_vector=PC/magn n12.%find the moment of the force PC about M12.%find the moment of

25、 the force PC about Mn n13.%this is the cross product of MP and PC13.%this is the cross product of MP and PCn n14.MP=P-M14.MP=P-Mnn15.moment=cross(MP,PC)15.moment=cross(MP,PC)当运行该代码后，返回的结果如下：当运行该代码后，返回的结果如下：PC=1 2 2unit_vector=0.4082 0.8165 0.4082moment=1 2-5 为了检查力矩计算结果，想象为了检查力矩计算结果，想象PC绕绕M点旋转，应用右手准

26、则点旋转，应用右手准则(zhnz)得到该力矩的旋转轴。大致说来，右手准则得到该力矩的旋转轴。大致说来，右手准则(zhnz)规定力矩的方向与在规定力矩的方向与在M点的一个标准右旋螺钉在该力的作用下旋转的方向相同。如果不用太精确的话，我们可以推断该力矩旋转轴大致沿点的一个标准右旋螺钉在该力的作用下旋转的方向相同。如果不用太精确的话，我们可以推断该力矩旋转轴大致沿z轴的负方向。显示的结果证实了这一估计是正确的。轴的负方向。显示的结果证实了这一估计是正确的。代码清单代码清单 求解向量问题求解向量问题(wnt)的脚本的脚本 第23页/共24页第二十三页，共24页。2023年2月7日SCHOOL OF PHYSICS AND TECHNOLOGY N.N.U.第 24 页感谢您的观看感谢您的观看(gunkn)！第24页/共24页第二十四页，共24页。