3Matlab矩阵.ppt

3 Matlab 矩阵运算 3.1 MATLAB矩阵 3.1.1 矩阵的建立 1．直接输入法 直接输入需遵循以下基本规则： 整个矩阵应以“ [ ]”为首尾，即整个输入矩阵必须包含在方括号中； 矩阵中，行与行之间必须用分号“ ；”或 Enter 键（ 按 Enter 键）符分隔； 每行中的元素用逗号“ ，”或空格分隔； 矩阵中的元素可以是数字或表达式，但表达式中不可包含未知的变量，MATLAB用表达式的值为该位置的矩阵元素赋值。 例3-1： A=[1,2,3,4;5,6,7,8;9,10,11,12;13,14,15,16] 2．利用M文件建立矩阵 对于比较大且比较复杂的矩阵，可以为它专门建立一个M文件。 启动编辑器（输入edit命令，或者使用按钮） 输入待建立矩阵，输入格式要求和直接输入法类似，只不过不同行元素必须使用分号作为间隔 3．利用冒号表达式建立一个向量 冒号表达式可以产生一个行向量，一般格式是： e1:e2:e3 其中e1为初始值，e2为步长，e3为终止值。 在MATLAB中，还可以用linspace函数产生行向量。其调用格式为： linspace(a,b,n) 其中a和b是生成向量的第一个和最后一个元素，n是元素总数。 x=0:0.5:2; x = 0 0.5000 1.0000 1.5000 2.0000 y=linspace(0,2,7) y = ? 0 0.3333 0.6667 1.0000 1.3333 1.6667 2.0000 z=[-1 x 3] z = -1.0000 0 0.5000 1.0000 1.5000 2.0000 3.0000 u=[y;z] u = 0 0.3333 0.6667 1.0000 1.3333 1.6667 2.0000 -1.0000 0 0.5000 1.0000 1.5000 2.0000 3.0000 3.1.2 矩阵的拆分 1．矩阵元素 通过下标引用矩阵的元素 例如: A(3,2)=200 采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。 在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依次类推。例如 A=[1,2,3;4,5,6]; A(3) ans = 2 2．矩阵拆分 (1) 利用冒号表达式获得子矩阵 ① A(:,j)表示取A矩阵的第j列全部元素； A(i,:)表示A矩阵第i行的全部元素； A(i,j)表示取A矩阵第i行、第j列的元素。 ② A(i:i+m,:)表示取A矩阵第i～i+m行的全部元素； A(:,k:k+m)表示取A矩阵第k～k+m列的全部元素， A(i:i+m,k:k+m)表示取A矩阵第i～i+m行内，并在第k～k+m列中的所有元素。 此外，还可利用一般向量和end运算符来表示矩阵下标，从而获得子矩阵。end表示某一维的末尾元素下标。 (2) 利用空矩阵删除矩阵的元素 在MATLAB中，定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。 可以通过给指定矩阵的行列赋以空矩阵来删除指定的行列。 3.1.3 特殊矩阵1．通用的特殊矩阵常用的产生通用特殊矩阵的函数有：zeros：产生全0矩阵(零矩阵)。ones：产生全1矩阵(幺矩阵)。eye：产生单位矩阵。rand：产生0～1间均匀分布的随机矩阵。randn：产生均值为0，方差为1的标准正态分布随机矩阵。 例3-3 分别建立3×3、3×2和与矩阵A同样大小的零矩阵 (1) 建立一个3×3零矩阵。 zeros(3) (2) 建立一个3×2零矩阵。 zeros(3,2) (3)设A为2×3矩阵，则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。 A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵A zeros(size(A)) %产生一个与矩阵A同样大小的零矩阵 2．用于专门学科的特殊矩阵 (1) 魔方矩阵魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵，其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n)，其功能是生成一个n阶魔方阵。 例3-5 将101~125