北航数学实验matlab2 幻灯片.pptVIP

矩阵的数组运算 数组运算：对应元素进行运算 点与算术运算符之间不能有空格！ 数组运算包括：点乘、点除、点幂 相应的数组运算符为： “.* ” ， “./ ” ， “.\ ” 和“ .^ ” 参与运算的对象必须具有相同的形状！ 例： A=[1 2 3; 4 5 6]; B=[3 2 1; 6 5 4]; C=A.*B; D=A./B; E=A.\B; F=A.^B; 函数取值 设 x 是变量， f 是一个函数 当 x = a 是标量时，f(x) = f(a)也是一个标量 当 x = [a, b, … , c] 是向量时，f(x)= [f(a), f(b), … , f(c)] 函数作用在矩阵上的取值 若 A 是矩阵，则 f(A) 是一个与 A 同形状的矩阵 f 作用在 x 的每个分量上 函数取值 怎样计算 eA ? 例： x=[0:pi/4:pi]; A=[1 2 3; 4 5 6]; y1=sin(x); y2=exp(A); y3=sqrt(A); 例： 矩阵的超越函数 Matlab 提供了三种矩阵函数：expm、sqrtm、logm 详情参见联机帮助（help expm / sqrtm / logm ） 更一般的矩阵函数： funm funm(A,@fun) 参数 fun 的可以是 exp,，log，cos，sin，cosh，sinh 矩阵与数的运算 加减：矩阵的每个元素都与数作加减运算 数乘：矩阵的每个元素都与数作乘法运算 矩阵除以一个数：每个元素都除以这个数 点幂： 底为矩阵，指数为标量 底为标量，指数为矩阵 数与数组的点幂 x.^y =[1^4,2^5,3^6]=[1,32,729] x.^2 =[1^2,2^2,3^2]=[1,4,9] 2 .^x = ? .^ 前面留个空格 例：x=[1 2 3]; y=[4 5 6]; 2 .^[x;y]= ? Matlab中的所有 标点符号必须在 英文状态下输入 Matlab中常见数学函数 sin、cos、tan、cot、sec、csc、… asin、acos、atan、acot、asec、acsc、… exp、log、log2、log10、sqrt abs、conj、real、imag、sign fix、floor、ceil、round、mod、rem max、min、sum、mean、sort、… det、inv、eig、rank、… …… log 是自然对数，即以 e 为底数 mod(x,y) 结果与 y 同号，rem(x,y) 则与 x 同号 max 等函数的参数是矩阵时，是作用在矩阵各列上 练习题 试分别生成 5 阶的单位阵、8 阶均匀分布的随机矩阵及其下三角矩阵 生产列向量 x=[1, 3, 5, 7, 9, … , 29] 生成以 x 的前 5 个元素为对角线的矩阵 A 生成一个与 A 同阶的正态分布的随机矩阵 B 输出 A 与 B 的 kronecker 乘积矩阵 C 生成由 B 与其同阶魔方矩阵点乘得到的矩阵 D 生成由 D 的第 2、4、5 行和第4、1 列组成的子矩阵 E 求出矩阵 E 的最大元素 数学软件 向量与矩阵运算 向量与矩阵的生成 向量与矩阵运算 向量的生成 直接输入: a=[1,2,3,4] 冒号运算符 a=[1:4] == a=[1, 2, 3, 4] b=[0:pi/3:pi] == b=[0, 1.0472, 2.0944, 3.1416] c=[6:-2:0] == c = [6, 4, 2, 0] 例： 从矩阵中抽取行或列 向量与矩阵的生成（续） 向量与矩阵运算 矩阵的生成 直接输入: A=[1, 2, 3; 4, 5, 6; 7, 8, 9] 由向量生成 由函数生成 通过编写m文件生成 例： x=[1,2,3];y=[2,3,4]; A=[x,y], B=[x;y] 例： C=magic(3) 常见矩阵生成函数 zeros(m,n) 生成一个 m 行 n 列的零矩阵，m=n 时可简写为 zeros(n) ones(m,n) 生成一个 m 行 n 列的元素全为 1 的矩阵,  m=n 时可写为 ones(n) eye(m,n) 生成一个主对角线全为 1 的 m 行 n 列矩阵, m=n 时可简写为 eye(n)，即为 n 维单位矩阵 diag(X) 若 X 是矩阵，则 diag(X) 为 X 的主对角线向量 若 X 是向量，diag(X) 产生以 X 为主对角线的对角矩阵 tril(A) 提取一个矩阵的下三角部分 triu(A) 提取一个矩阵的上三角部分 rand(m,n) 产生 0～1 间均匀分布的随机矩阵 m=n