[工学]MATLAB数值计算.docVIP

第5章 MATLAB数值计算 MATLAB的数学计算＝数值计算＋符号计算 其中符号计算是指使用未定义的符号变量进行运算，而数值计算不允许使用未定义的变量。.1 矩阵和数组 MATLAB最基本也是最重要的功能就是进行实数或复数矩阵的运算。矩阵和数组的概念 在MATLAB的运算中，经常要使用标量、向量、矩阵和数组，这几个名称的定义如下： 标量：是指1×1的矩阵，即为只含一个数的矩阵。 向量：是指1×n或n×1的矩阵，即只有一行或者一列的矩阵。 矩阵：是一个矩形的数组，即二维数组，其中向量和标量都是矩阵的特例，0×0矩阵为空矩阵([])。 数组：是指n维的数组，为矩阵的延伸，其中矩阵和向量都是数组的特例。.1.1矩阵输入(1) 矩阵元素应用方括号([])括住； (2) 每行内的元素间用逗号或空格隔开； (3) 行与行之间用分号或回车键隔开； (4) 元素可以是数值或表达式。 1. 通过显式元素列表输入矩阵 c=[1 2;3 4;5 3*2] % [ ]表示构成矩阵,分号分隔行,空格分隔元素 c = 1 2 3 4 5 6 用回车键代替分号分隔行: c=[1 2 3 4 5 6] 1 2 3 4 5 6 2. 通过语句生成矩阵 (1) 使用from:step:to方式生成向量(2) 使用linspace和logspace函数生成向量 3. 由矩阵生成函数产生特殊矩阵 MATLAB提供了很多能够产生特殊矩阵的函数，各函数的功能如表5.2所示。 表5.2 矩阵生成函数 函数名 功能 例子 输入 结果 zeros(m,n) 产生m×n的全0矩阵 zeros(2,3) ans = 0 0 0 0 0 0 ones(m,n) 产生m×n的全1矩阵 ones(2,3) ans = 1 1 1 1 1 1 rand(m,n) 产生均匀分布的随机矩阵，元素取值范围0.0～1.0。 rand(2,3) ans = 0.9501 0.6068 0.8913 0.2311 0.4860 0.7621 randn(m,n) 产生正态分布的随机矩阵 randn(2,3) ans = -0.4326 0.1253 -1.1465 -1.6656 0.2877 1.1909 magic(N) 产生N阶魔方矩阵(矩阵的行、列和对角线上元素的和相等) magic(3) ans = 8 1 6 3 5 7 4 9 2 eye(m,n) 产生m×n的单位矩阵 eye(3) ans = 1 0 0 0 1 0 0 0 1 注意： zeros、ones、rand、randn和eye函数当只有一个参数n时，则为n×n的方阵; 当eye(m,n)函数的m和n参数不相等时则单位矩阵会出现全0行或列。 【例.3】查看eye函数的功能。 X1=eye(2,3) X1 = 1 0 0 0 1 0 X2=eye(3,2) X2 = 1 0 0 1 0 0 4. 通过MAT数据文件加载矩阵 通过“load”命令或选择菜单“File”→“Import Data”命令加载MAT数据文件来创建矩阵。 5. 在M文件中创建矩阵 M文件实际上是一种包含MATLAB代码的文本文件 通过在MATLAB命令窗口中运行M文件创建矩阵。 .1.2矩阵元素和操作 矩阵和多维数组都是由多个元素组成的，每个元素通过下标来标识。 1. 矩阵的下标 (1) 全下标方式 矩阵中的元素可以用全下标方式标识，即由行下标和列下标表示，一个m×n的矩阵的第i行第j列的元素表示为 (i,j)。 注意： 如果在提取矩阵元素值时，矩阵元素的下标行或列(i,j)大于矩阵的大小(m,n)，则MATLAB会提示出错； 而在给矩阵元素赋值时，如果行或列(i,j)超出矩阵的大小(m,n)，则MATLAB自动扩充矩阵，扩充部分以0填充。 a=[1 2;3 4;5 6] a = 1 2 3 4 5 6 a(3,