第二章matlab02数值运算功能1.ppt

1）四舍五入取整：round(A) %将矩阵A中的元素按最近的整数取整,即四舍五入取整 2）按离0近的方向取整： fix(A) %将矩阵A中元素按离0近的方向取整 3）floor(x)向负无穷方向取整,即不足取整 4）ceil(x)向正无穷方向取整,即过剩取整 (二)取整函数 v=[-2.6 -2.5 -2.4 4.2 4.5 4.8]; n=[round(v); fix(v); floor(v); ceil(v)] n = -3 -3 -2 4 5 5 -2 -2 -2 4 4 4 -3 -3 -3 4 4 4 -2 -2 -2 5 5 5 例：A=[1 0.2 -1 -0.3 0]； x=sign(A) x = 1 1 -1 -1 0 xnew=(sign(A)+1)/2 xnew = 1.0000 1.0000 0 0 0.5000 xnew=floor((sign(A)+1)/2) xnew = 1 1 0 0 0 (三)逻辑函数 1）all函数 所有向量元素（对于矩阵来说则指列向量）非零时，输出为1；否则输出为0 all(1:9); 输出为1 all([1:9 0]); 输出为0 all([1:3;1:3;2:4]);输出为[1 1 1 ] all([1:3;1:3;0:2]);输出为［0 1 1 ] 2）any函数: 向量中有非零元素（对于矩阵来说则指列向量）时，输出为1；否则，输出为0。 any([zeros(1,9) 1]);输出为1 any(zeros(1,9) );输出为0 x=find(a= =k) ％产生a中等于k的元素的索引值，没有则生成空矩阵 [i,j]=find(a) %产生a中非零元素的行列位置 (4)元素检索函数：find [i,j]=find(a= =k) ％产生a中等于k的元素的行列位置，没有则生成空矩阵 a=[2 0 3 1 0 0]; [i,j]=find(a= =1) i = 2 j = 1 a(find(a==1)) a=1 a=[1 2 0;3 0 5;7 8 9］ a =1 2 0 3 0 5 7 8 9 a(3,3)=0 a =1 2 0 3 0 5 7 8 0 可用find函数修改 可以利用下标对矩阵元素进行修改。 例:矩阵的修改 a(find(a==9))=0 4.矩阵的一些特殊操作 矩阵的变维 a=[1:12]; b=reshape(a,3,4); a(:)%列向量 c=zeros(3,4);c(:)=a(:) 矩阵的变向 rot90:旋转; fliplr:左右翻; flipud:上下翻 矩阵的抽取 diag:抽取矩阵主对角线;或以向量为主对角 tril: 抽取主下三角; triu:抽取主上三角 矩阵的扩展 【例1】 reshape与 diag的使用演示。 a=-4:4 % 产生一维数组 a =-4 -3 -2 -1 0 1 2 3 4A=reshape(a,3,3) ％注意元素个数不能变 % 把一维数组 a 重排成3×3矩阵 A = -4 -1 2 -3 0 3 -2 1 4 按索引顺序排 A=reshape(a,length(a),1) %行向量变为列向量 D =1 2 34 1 56 7 1 diag(diag(D)) % 内diag 取 D 的对角元素，外 diag 利用向量生成对角阵ans =1 0 00 1 00 0 1 diag(D)ans =111 X=diag(v) %以向量v为主对角线元素，其余元素为0构成矩阵X。 diag(D) % 取 D 矩阵的主对角元素 v =[1 1 1] a1=diag(A,1) % 取 A 阵“第一上对角线”的元素 a1 = -1 3 A1=diag(a1,-1) % 产生以 a1 数组元素为“第一下对角线”元素的二维数组A1 = 0 0