[工学]MATLAB B第九讲.pdf

《MATLAB数据处理与应用》2010-2011 学年选修课 第九讲 第九讲 MATLAB 的微分与微分方程数值分 MATLAB 的微分与微分方程数值分 析 析 王文健 王文健 wwj527@163.com wwj527@163.com Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY 主要内容 数值积分 数值微分 一阶常微分方程(ODE) 的数值解法 高阶ODE 和OED 组的的数值解法 傅里叶变换与泰勒级数展开 2 西南交通大学摩擦学研究所 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY MATLAB的微积分与微分方程 简介 微分与积分是数学分析中的重要概念，虽然很多微积 分问题可以用解析的方式得到精确答案 但有更多的微积分问题没有解析解存在，必须借助数 值方法求解 MATLAB在数值积分与微分数值解法存在很大的优点 3 西南交通大学摩擦学研究所 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY MATLAB的微积分与微分方程 数值积分 对于积分运算有梯形法trapz(x,y)和二次函数法quad() 两种运算函数 梯形法trapz(x,y) 是MATLAB提供的最简单的积分函 数，式中的x 、y 分别代表数目相同的矩阵或数组，y 是 x 的函数 注意：trapz(x,y) 只针对离散型的数据做积分 4 西南交通大学摩擦学研究所 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY MATLAB的微积分与微分方程 数值积分— 梯形法 举例：求函数 K ∫π sin(x)dx 在区间[0,pi]上的积分值 0 对于此积分很容易得到精确值为2 x=0:pi/100:pi;