矩阵及基本运算

矩阵及基本运算第一页，共三十七页，2022年，8月28日矩阵的建立：2直接输入以左方括号开始，右方括号结尾同一行中的元素用逗号或空格隔开行与行之间用分号或回车隔开>>a=[1,2,3;4,5,6;7,8,9]a=123456789>>a=[123;456;789]>>a=[123456789]这种方法适用于矩阵元素个数较少的情况第二页，共三十七页，2022年，8月28日3用矩阵编辑器输入在命令窗口创建变量双击它打开矩阵编辑器在编辑器中输入或修改，完毕后关闭即可第三页，共三十七页，2022年，8月28日4用矩阵函数创建矩阵zeros(m,n):m×n阶零矩阵eye(m,n):n阶单位矩阵ones(m,n):m×n阶元素全为1的矩阵rand(m,n):m×n阶均匀随机矩阵，取值0到1randn(m,n):m×n阶正态随机矩阵，取值0到1magic(n):n阶魔方矩阵(矩阵每行、每列及两对角线上的元素和相等)hilb(n):n阶Hilbert矩阵，即矩阵元素为1/(i+j-1)>>A=ones(2,3)A=111111>>eye(3)ans=100010001第四页，共三十七页，2022年，8月28日5A=diag(v,k)生成主对角线方向上的第k(整数)层元素为向量v的矩阵.规定：k=0时，表示矩阵的主对角线k>0时，表示主对角线的平行位置上方的第k层K<0时，表示主对角线的平行位置下方的第-k层v=diag(A,k)提取矩阵A中主对角线方向上第k(整数)层元素第五页，共三十七页，2022年，8月28日例利用函数生成矩阵6B=0100002000030000A=100020003>>v=[123]>>A=diag(v,0)>>B=diag(v,1)第六页，共三十七页，2022年，8月28日例随机生成数值在10到30之间的含5个元素的行向量7>>10+(30-10)*rand(1,5)ans=25.844129.189823.114810.714226.9826第七页，共三十七页，2022年，8月28日例生成三对角矩阵8A=120000112000021200003120000112000021>>a1=ones(1,6)>>a2=2*ones(1,5)>>a3=[12312]>>A=diag(a1,0)+diag(a2,1)+diag(a3,-1)第八页，共三十七页，2022年，8月28日9通过文件生成有时我们需要处理一些没有规律的数据，或数据量很大时，如在命令窗口输入，清除后再次使用需要重新输入，这就增加工作量.解决方案：直接把数据作为矩阵输入到M文件（Matlab的程序代码文件）中作为变量保存到MAT文件中SaveFilenamevar1var2var3LoadFilename第九页，共三十七页，2022年，8月28日数组的生成：10使用from:step:to生成数组

当step省略时，表示步长step=1当step为负数时，可以创建降序的数组除前面的矩阵输入外，还有>>a=1:0.5:2a=1.00001.50002.0000>>b=1:3b=123>>c=3:-1:1c=321第十页，共三十七页，2022年，8月28日11linspace(a,b,n)生成从a到b之间线性分布的n个元素的数组.如果n省略，则默认为100.logspace(a,b,n)生成从10^a到10^b之间按对数等分的n个元素的数组.如果n省略，则默认为50.>>a=linspace(1,9,5)a=13579>>b=logspace(1,5,3)b=101000100000logspace(a,b,n)=10.^linspace(a,b,n)第十一页，共三十七页，2022年，8月28日复数矩阵的输入12>>A=[1+2i3+4i;5+6i7+8i]>>A=[13;57]+i*[24;68]设A=1+2i3+4i5+6i7+8i第十二页，共三十七页，2022年，8月28日矩阵元素操作MATLAB利用下标来访问矩阵中的元素，下标可以是向量.13设A=123456789>>b=A(2,3)b=6>>c=A([12],[23])c=2356第十三页，共三十七页，2022年，8月28日如果在某个下标位置只有一个：号，则表示取对应的所有列或行.用end表示某一维数中的最大值.14设A=123456789>>b=A(2,:)b=456>>c=A(:,2:3)c=235689>>d=A(1:2,end)d=36第十四页，共三十七页，2022年，8月28日用单下标来表示

即将矩阵的所有列按先左后右的次序接成“一维长列”，然后再对元素位置进行编号15设A=123456789>>A(4)ans=2>>A(2:5)ans=4725第十五页，共三十七页，2022年，8月28日矩阵的合并16A=123456B=78910C=111213>>[A,B]ans=12378456910>>[A;C]ans=123456111213第十六页，共三十七页，2022年，8月28日通过空矩阵“[]”对矩阵元素进行删除17设A=123456789>>A(2,:)=[]A=123789>>A=[]A=[]>>isempty(A)ans=1第十七页，共三十七页，2022年，8月28日矩阵的运算基本运算加（+）：A+B,A+k(每个元素都加k)减（-）：A-B,A-k乘（*）：A*B左除（\）：A\B,即A-1B,A必须为方阵右除（/）：A/B,即AB-1,B必须为方阵乘幂（^）:A^n,A必须为方阵转置(‘):transpose(A)或A’点运算(.):其加、减、乘、除和乘方都是对两个相同维数的矩阵进行对应元素的运算.18第十八页，共三十七页，2022年，8月28日19A=123321125B=-29-15615-21-611-2例已知求出AX=B和XA=B的解>>A\Bans=2.00004.00003.50001.00001.0000-19.0000-2.00001.00006.5000>>inv(A)*Bans=2.00004.00003.50001.00001.0000-19.0000-2.00001.00006.5000>>B/Aans=25.2500-3.2500-17.500030.7500-0.7500-22.500026.2500-5.7500-15.0000>>B*inv(A)ans=25.2500-3.2500-17.500030.7500-0.7500-22.500026.2500-5.7500-15.0000第十九页，共三十七页，2022年，8月28日20A=123456B=78-9-1-3-2例已知>>A.*Bans=716-27-4-15-12>>A./Bans=0.14290.2500-0.3333-4.0000-1.6667-3.0000>>A.\Bans=7.00004.0000-3.0000-0.2500-0.6000-0.3333>>A.^3ans=182764125216第二十页，共三十七页，2022年，8月28日关系运算关系运算符：<(小于)、<=(小于等于)、>(大于)、>=(大于等于)、==(等于)、~=(不等于)关于运算在对应元素之间进行，结果是一个二值矩阵，其中0表示“假”，1表示“真”.21相关函数I=find(A)返回矩阵A的非零元素的线性指标[IJ]=find(A)返回矩阵A的非零元素的行、列指标，分别存于I和JA=1004>>find(A)ans=14>>[IJ]=find(A)I=12J=12第二十一页，共三十七页，2022年，8月28日22>>A=[1432];B=[5413];>>A>B%比较A与B矩阵的元素大小，输出逻辑数值ans=0010>>x=find(A==B)%找出A与B中对应元素相等的下标x=2>>A(find(A==B))%找出A和B中对应元素相等的元素ans=4例第二十二页，共三十七页，2022年，8月28日23例找出以下矩阵中绝对值大于3的所有元素>>A=zeros(2,5);>>A(:)=-4:5A=-4-2024-3-1135>>L=abs(A)>3L=1000100001A=-4-2024-3-1135>>XL=A(find(L))XL=-445第二十三页，共三十七页，2022年，8月28日逻辑运算逻辑运算符：&(与)、|(或)、~(非)、xor(异或)、&&(先决与)、||(先决或)24xyx&yxIy~xxor(x,y)000010010111100101111100&&:当该运算符的左边为1时，才执行右边的运算||:当该运算符的左边为0时，才执行右边的运算第二十四页，共三十七页，2022年，8月28日25A=1234B=4321例设>>(A>B)&(A==B)ans=0000>>(A>B)&(A>1)ans=0011>>(A>B)|(A>1)ans=0111>>xor(A>B,A>1)ans=0100第二十五页，共三十七页，2022年，8月28日其它一些常用矩阵函数size(A):返回矩阵的行数和列数size(A,1):返回行数size(A,2):返回列数length(a):返回向量a的长度reshape(A,m,n):重排矩阵A为m行n列矩阵，A必须为m*n个元素26A=123456a=12345>>size(A)ans=23>>size(A,2)ans=3>>reshape(A,3,2)ans=154326>>length(a)ans=5第二十六页，共三十七页，2022年，8月28日max(a):返回向量a的最大元素max(A,[],dim):返回矩阵沿着第dim维方向的最大元素向量max(A,B):返回矩阵A,B中对应位置上的较大元素min:返回最小元素，用法类似maxmean(a):返回向量a的平均值mean(A,dim):返回矩阵A第dim维方向上平均值median:返回中位数，用法类似meansum：求元素和，用法类似meanprod：求元素积，用法类似meancumsum:求元素累积和，用法类似meansort(a)：对向量a进行升序排序sort(A,dim,mode):对矩阵A沿着第dim维方向排序，mode:’ascend’(升序，默认)、’descend’(降序)27第二十七页，共三十七页，2022年，8月28日28a=12565931A=123456>>max(A,[],1)ans=456>>max(A,[],2)ans=36>>max(A)ans=456>>max(a)ans=65>>min(a)ans=1第二十八页，共三十七页，2022年，8月28日29a=12565931A=123456>>sum(a)ans=95>>cumsum(a)ans=121782919495>>sum(A,1)ans=579>>sum(A,2)ans=615>>sum(A)ans=579第二十九页，共三十七页，2022年，8月28日30a=12520931A=123456>>prod(a)ans=32400>>prod(A,1)ans=41018>>prod(A,2)ans=6120>>prod(A)ans=41018第三十页，共三十七页，2022年，8月28日31a=12520931A=123456>>sort(a)ans=13591220>>sort(a,'ascend')ans=13591220>>sort(a,'descend')ans=20129531>>sort(A,1)ans=123456>>sort(A,2)ans=123456>>sort(A)ans=123456第三十一页，共三十七页，2022年，8月28日例计算以下例子321+2+…+1001*2*…*1006+9+12+…+297+300>>sum(1:100)ans=5050>>prod(1:100)ans=9.3326e+157>>sum(6:3:300)ans=15147阶乘：factorial(N)=prod(1:N)>>factorial(100)ans=9.3326e+157第三十二页，共三十七页，2022年，8月28日33设a=12510931A=123456用多种方法求数组a的算术平均值，矩阵A每行的算术平均值>>mean(a)ans=6.6667>>sum(a)/length(a)ans=6.6667>>sum(a)/size(a,2)ans=6.6667>>mean(A,2)ans=25>>sum(A,2)/size(A,2)ans=25第三十三页，共三十七页，2022年，8月28日34例设a=[91,98,5