第7章 微分方程问题的计算机求解.ppt

星期二, 2013-6- 18, 13:45:26 7.8 微分方程的框图求解 Simulink简介 Simulink相关模块 微分方程的Simulink建模与求解 星期二, 2013-6- 18, 13:45:26 7.8.1 Simulink简介 1990 年前后出现最早的 Simulink，当时名为SimuLAB，1992 年改为 Simulink Simulink 的名字有两重含义 仿真 (simu) 与模型/模块 连接 (link) 使用odegroup命令可以打开自定义模块集 或单击 图标启动 Simulink 星期二, 2013-6- 18, 13:45:26 7.8.2 Simulink相关模块 常用的模块： 输入输出端口(In1, out1) 时钟模块：产生时间变量t 常用输入模块： Sine模块产生正弦信号 Step模块产生阶跃信号 Constant模块产生恒值信号 积分器模块(Int)：描述一阶导数，搭建积分的动态系统 延迟模块(Transport Delay) 星期二, 2013-6- 18, 13:45:26 增益模块： Gain：主要用于信号的放大 Sliding Gain：滚动杆，实时地改变模块的增益 Matrix Gain：矩阵增益模块 数学运算模块 +, -, *, ÷ 数学函数模块 Sin, exp 信号向量化模块 Mux：将若干路信号混成向量型的信号 DeMux：将向量型信号解出单路的信号 星期二, 2013-6- 18, 13:45:26 7.8.3 微分方程的 Simulink建模与求解 建立起相应的微分方程的 Simulink 模型 可以用sim()函数对其模型直接求解，该函数的调用格式 得出微分方程的数值解 星期二, 2013-6- 18, 13:45:26 例 7.43 求解 Lorenz 方程的问题，其中，b = 8/3, r = 10, s = 28，各个状态变量的初值为 x1(0)= x2(0)= 0, x3(0)=10-10，试用Simulink搭建该模型并得出仿真结果 定义关键信号，引入积分器模块 星期二, 2013-6- 18, 13:45:26 搭建Simulink框图 星期二, 2013-6- 18, 13:45:26 创建文件c7mlor1b.mdl 求解系统 与函数ode45()的结果进行比较 仿真参数设置 状态变量初值的处理 仿真结果检验 星期二, 2013-6- 18, 13:45:26 例 7.44 延迟微分方程式 使用延迟模块定义关键信号 用Simulink搭建该微分方程模型 求解微分方程 星期二, 2013-6- 18, 13:45:26 搭建Simulink框图， 以文件c7mdde2.mdl存在MATLAB的工作空间下 星期二, 2013-6- 18, 13:45:26 设定新的状态变量： 得出一阶状态方程模型： 星期二, 2013-6- 18, 13:45:26 例 7.45 其中，t1=0.15, t2=0.5 用Simulink搭建系统模型，并进行仿真 星期二, 2013-6- 18, 13:45:26 搭建Simulink框图，以文件c7mdde3.mdl存在MATLAB的工作空间下 星期二, 2013-6- 18, 13:45:26 另一个模型 星期二, 2013-6- 18, 13:45:26 输入已知的矩阵： MATLAB求解语句： 星期二, 2013-6- 18, 13:45:26 例7-47 非线性随机输入系统仿真 建立仿真模型 星期二, 2013-6- 18, 13:45:26 仿真模型 星期二, 2013-6- 18, 13:45:26 本章内容简介 本章介绍了基于 MATLAB 符号运算工具箱 dsolve()函数的线性微分方程的解析解方法，并介绍基于该函数的特殊非线性微分方程的解析解。 对一般非线性微分方程来说，解析解是不存在的，只能依赖数值解的方法对其进行研究。 星期二, 2013-6- 18, 13:45:26 7.5.2 非线性方程边值问题的打靶算法 假定原始问题可以转换成下面的初值问题 则问题转换成求解 Newton 迭代法求取 m 星期二, 2013-6- 18, 13:45:26 式中 且可以由下面的微分方程初值问题来求解 星期二, 2013-6- 18, 13:45:26 其中，要求能显式地求出下式： 构造MATLAB函数实现上述算法： 星期二, 2013-6- 18, 13