第10章 数学问题的非传统解法解析.ppt

星期二, 2013-5- 21, 05:50:51 例10-51 框图求解零初值线性分数阶微分方程 引入 变换后 星期二, 2013-5- 21, 05:50:51 星期二, 2013-5- 21, 05:50:51 例10-52 非线性分数阶微分方程（零初值） 变换形式 星期二, 2013-5- 21, 05:50:51 星期二, 2013-5- 21, 05:50:51 非零初值的 Caputo 微分方程数值解法 需要对初值补偿 构造封装模块 星期二, 2013-5- 21, 05:50:51 取 定义辅助信号 由poly2caputo() 填写Fcn模块的补偿函数 可以封装 Caputo 微分算子模块 仿照前面例子建模 星期二, 2013-5- 21, 05:50:51 例10.53 非线性分数阶微分方程 已知初值 解析解 辅助信号 等效于 星期二, 2013-5- 21, 05:50:51 微分方程变换 直接得出 化简 新微分方程 星期二, 2013-5- 21, 05:50:51 星期二, 2013-5- 21, 05:50:51 10.6.6 对象编程实例——分数阶传递函数模型 分数阶传递函数（FOTF）对象 FOTF对象的创建 FOTF短信的重载函数编写 面向对象的分析与设计方法 星期二, 2013-5- 21, 05:50:51 FOTF对象定义 FOTF类用于描述分数阶传递函数，其作用类似于MATLAB控制系统工具箱中的“tf” 建立类的最基本步骤是建立一个@FOTF目录 在目录中编写两个函数fotf()和display()，前者描述该类的内容，后者描述该类的显示方法 还可以选择性地建立其他函数用于函数重载 plus(),minus(),mtimes(),feedback() unique(), bode() 等 星期二, 2013-5- 21, 05:50:51 建立起fotf.m文件 星期二, 2013-5- 21, 05:50:51 显示函数 星期二, 2013-5- 21, 05:50:51 例 10.54 给定 试将分数阶传递函数输入到MATLAB的工作空间 MATLAB求解语句： 星期二, 2013-5- 21, 05:50:51 用Grünwald-Letnikov定义求解分数阶微分 数学描述 其中， 星期二, 2013-5- 21, 05:50:51 构造一个MATLAB函数 函数调用格式 分数阶微分数值计算 星期二, 2013-5- 21, 05:50:51 例10-42 阶跃函数的导数和积分是什么？ 整数阶的积分与微分 对不同的阶次 星期二, 2013-5- 21, 05:50:51 例 10.43 考虑一个初值非零的函数 试求出其分数阶导数。 MATLAB求解语句： 星期二, 2013-5- 21, 05:50:51 对不同的 g 选值，可以调用下面的语句绘制出分数阶导函数的三维图 星期二, 2013-5- 21, 05:50:51 例 10.44 试用不同定义求取函数为 f (t) = sin(3t + 1) 的0.75阶微分，并比较得出的结果。 Cauchy公式结果 MATLAB求解语句： 星期二, 2013-5- 21, 05:50:51 Caputo 微积分定义的数值计算 Caputo 分数阶积分与 Grünwald-Letnikov 定义完全一致 Caputo / Grünwald-Letnikov 微分关系 MATLAB 函数实现 星期二, 2013-5- 21, 05:50:51 MATLAB 代码 星期二, 2013-5- 21, 05:50:51 例10-45 函数为 f (t) = sin(3t + 1) ，求 0.3、1.3、2.3 阶的Caputo导数 0.3阶两个定义之差为 绘制 0.3 阶导数 星期二, 2013-5- 21, 05:50:51 其余两个导数曲线绘制 星期二, 2013-5- 21, 05:50:51 Oustaloup滤波算法 如果函数 f (t) 的表达式不是事先已知 连续滤波器传递函数模型：Oustaloup算法 感兴趣区间 星期二, 2013-5- 21, 05:50:51 构造MATLAB函数，设计连续滤波器 函数调用格式 Oustaloup 滤波器求解 星期二, 2013-5- 21, 05:50