matla第2篇 矩阵运算基础.pptVIP

第二章 矩阵运算基础 计算机及信息工程学院 MATLAB的运算事实上是以数组 (array) 及矩阵 (matrix) 方式在做运算，而这二者在MATLAB的基本运算性质不同，数组强调元素对元素的运算，而矩阵则采用线性代数的运算方式。 2.1 矩阵的创建 而宣告一变量为数组或是矩阵时，如果是要个别键入元素，须用中括号[ ] 将元素置于其中。数组为一维元素 所构成，而矩阵为多维元素所组成，例如 x = [1 2 3] % 一维 1x3 阵列 矩阵的同行元素之间用空格或逗号分隔； x = [1 2 3; 4 5 6] % 二维 2x3 矩阵，以;区隔各列的元素 x = [1 2 3 % 二维 2x3 矩阵，各列的元素分二行键入 4 5 6] ％矩阵的尺寸不必预先定义 假设要计算 y = sin (x), 0=x=π而x = 0, 0.2π, 0.4π,...,π，即可用数组方式运算，例如 x = [0 0.2*pi 0.4*pi 0.6*pi 0.8*pi pi] % 注意数组内也可作运算 x = 0 0.6283 1.2566 1.8850 2.5133 3.1416 y=sin(x) y = 0 0.5878 0.9511 0.9511 0.5878 0.0000 要找出数组的某个元素或数个元素，可参考以下的例子 x(3) % 第三个x的元素 ans = 1.2566 y(5) % 第五个y的元素 ans = 0.5878 x(1:5) % 列出第一到第五个x的元素 ans = 0 0.6283 1.2566 1.8850 2.5133 y(3:-1:1) % 列出第三到第一个y的元素，3为起始值，1为终止值，-1为增量 ans = 0.9511 0.5878 0 x(2:2:6) % 列出第二到第六个x的元素，2为起始值，6为终止值，2为增量 ans = 0.6283 1.8850 3.1416 y([4 2 5 1]) % 列出y元素，排列元素依序为原来y数组的4,2,5,1个 ans = 0.9511 0.5878 0.5878 0 (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=[]。 注意，X=[]与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间，只是维数为0。 将某些元素从矩阵中删除，采用将其置为空矩阵的方法就是一种有效的方法。 3. 建立大矩阵 大矩阵可由方括号中的小矩阵建立起来。例如 A=[1 2 3;4 5 6;7 8 9]; C=[A,eye(size(A)); ones(size(A)),A] C = 1 2 3 1 0 0 4 5 6 0 1 0 7 8 9 0 0 1 1 1 1 1 2 3 1 1 1 4 5 6 1 1 1 7 8 9 2.1.1 命令窗口直接输入 【例2-1】在命令窗口创建简单的数值矩阵 在命令窗口输入： x= [1 2 3 ； 4 5 6 ； 7 8 9] x= 1 2 3 4 5 6 7 8 9 2.1.2 通过M文件创建矩阵 当矩阵的尺寸较大时，直接在命令窗口输入矩阵元素，容易出错误且不便修改。为了解决此问题，可以先将矩阵按创建原则写入一个M文件中，在MATLAB的命令窗口或程序中直接执行该M文件，即将矩阵调入工作空间。 2.1.3 通过函数创建矩阵 【例2-3】 创建0到2*π间的正弦函数矩阵。 x=0:pi/4:2*pi; y=sin(x) y = Columns 1 through 5 0 0.7071