Matlab计算与仿真交互式计算(续).ppt

第五讲提纲 第四讲回顾 交互式计算(续) Matlab编程(I) 第四讲回顾 矩阵计算进阶 - 计算: 算术/非算术运算/矩阵函数 - 整形: 查询/重排/翻转/拼接 关于复数 - 内置操作函数 关于字符串 - 转换/查询/拼接 交互式计算 (续) 综合练习三: 基于函数eye创建矩阵 - 提示1: 借助矩阵翻转操作 回忆: lookfor flip??? - 提示2: 矩阵的算术运算 回忆: A*c; %=[aij*c] 交互式计算 (续) - 提示3: 矩阵左右翻转函数 回忆: fliplr(A); % 列指标逆序矩阵 - 求解 Tem=fliplr(eye(3)); % 单位阵的左右翻转 (-2)*Tem % ans = 0 0 -2 0 -2 0 -2 0 0 交互式计算 (续) 综合练习四: 求解下述线性方程组 - 求解: A=[2 1 -5 1; 1 -3 0 -6;… 0 2 -1 2; 1 4 -7 6]; b=[8; 9; -5; 0] x=b’/A % x=A\b 交互式计算 (续) 关于多项式 - 基本概念 数学表达式 求值计算 c=[c1,c2,c3,…,cn,cn+1]; polyval(c,x) 交互式计算 (续) 卷积计算 w=conv(u,v) % 返回多项式系数向量 重构操作 c=poly(v) % 返回多项式系数向量 微分操作 cdef=polyder(c) 拟合操作 c=polyfit(x,y,n) 交互式计算 (续) 使用实例 c = [3 2 1]; polyval(c,[5 7 9]) ans = 85 162 262 x = (0: 0.1: 2.5); y = erf(x); c = polyfit(x,y,6) c = 0.0084 -0.0983 0.4217 -0.7435 0.1471 1.1064 0.0004 交互式计算 (续) d=conv(c,[5 7 9]) ans = 15 31 46 25 9 dd=polyder(d) dd = 60 93 92 25 poly(roots(dd)) ans = 1.0000 1.5500 1.5333 0.4167 Matlab编程 程序设计概述 - 基本概念 程序: 数据结构+算法 注解: 数据结构定义操作对象 算法定义求解过程 - 必要性 问题求解需要复杂算法 交互计算方式难以满足 Matlab编程 (续) 程序结构 - 顺序结构 逐条顺序执行程序指令或模块 - 选择结构 依据条件选择执行程序指令或模块 if 语句 % 条件选择 switch 语句 % 多值切换 Matlab编程 (续) Matlab编程 (续) M-文件概述 m-编辑器: 创建以m为结尾标志的asccii文件 分类：脚本文件/函数文件 - 脚本文件(script) 执行批处理行命令 文件名：默认的宏命令 全局性：与工作空间共享变量名 适用性：实现小规模简单的问题求解 Matlab编程 (续) Matlab编程 (续) Matlab编程 (续) - 程序运行 程序定位查询 用法举例: which drawPlot D:\mat