Matlab计算与仿真38524.ppt

第五讲提纲 上機實習提示 交互式计算-III 上机实习提示 上机实习提示 - 时间/地点安排 单周二：12-13节 实验中心B座401 - 基本要求 准时出席/遵守实验室规定 完成实习指南要求的任务 认真做好实习报告或总结 交互式计算 (续) 矩阵计算进阶 - 算术运算 加减运算: 相同阶数或一个为标量 提示: A+b 等价于 [aij+b] 乘法运算: 相容阶数或一个为标量 提示: A*b 等价于 [aij*b] 除法运算: 相容阶数且除矩阵非奇异或非零标量 提示: A/b 等价于 [aij/b] 注意: 左右除的分别, 即 A/B 不等价于 B/A 交互式计算 (续) - 非算术运算 关系运算: = = 结果: 返回布尔值矩阵 条件: 相同阶数或一个为标量 提示: A==b 等价于 [aij==b] 逻辑运算: /|/~/xor 以逻辑和为例, 相同维数或至少一个为标量 C=A B; % cij=1 iff both aij and bij nonzeros cij=0 otherwise 交互式计算 (续) - 矩阵函数 数学函数: 函数名通常以m为结尾标志 提示: 矩阵函数 和 函数矩阵的差异 例如, expm(X)≠exp(X) logm(X)≠log(X) 逻辑函数 相同维数或至少一个为标量 AND/OR/XOR 提示: 与逻辑运算一致 交互式计算 (续) 交互式计算 (续) 交互式计算 (续) 交互式计算 (续) - 矩阵整形操作 重排操作 reshape 基本格式: reshape(A,m,n); 功能: 按列重排 举例: A=[1:4; 5:8; 9:12]’; B=reshape(A,2,6) B = 1 3 5 7 9 11 2 4 6 8 10 12 提示: 对二维信号处理特别有用 交互式计算 (续) - 矩阵整形操作 翻转操作 fliplr(水平翻转)/flipud(上下翻转)/rot90(绕转) 基本格式: fliplr(A); flipud(A); rot90(A); 功能: 按水平/上下翻转重排, 或逆时针旋转90度 举例: rot90(A) A=[1:3; 4:6]; ans = fliplr(A) flipud(A) 3 6 ans = ans = 2 5 3 2 1 4 5 6 1 4 6 4 5 1 2 3 交互式计算 (续) 交互式计算 (续) - 矩阵整形操作 块提取操作: diag(对角元素)/tril(下三角)/triu(上三角) 基本格式: diag(A); % 返回对角向量 tril(A); % 上三角元素清零 triu(A); % 下三角元素清零 功能:提取对角块/下三角/上三角块 举例: triu(A) A=[1:3; 4:6]; ans =