8-常微分方程数值解法.pptVIP

常微分方程的数值解法 引入 微分方程的数值解法是动态系统仿真的基础。 思考：数值分析课程-计算机求解数学问题？ 仿真软件的实现（具体执行步骤） ？ 常微分方程的数值解法 Euler法 Runge-Kutta法 Adams算法 Gear算法 Matlab下的常微分方程求解函数 二阶、三阶龙格库塔法ode23() 四阶、五阶龙格库塔法ode45() 自适应变步长求解法 Matlab下的常微分方程求解函数 问题描述： 调用格式： [t,x]=ode23(方程函数名,tspan,x0,选项,附加参数) [t,x]=ode45(方程函数名,tspan,x0,选项,附加参数) 选项可以通过odeget(),odeset()函数来设置，通常采用默认值即可。 Matlab下的常微分方程求解函数 调用格式： [t,x]=ode23(方程函数名,tspan,x0,选项,附加参数) [t,x]=ode45(方程函数名,tspan,x0,选项,附加参数) 说明： 方程函数名:为描述系统状态的M函数的名称,用单引号括住 tspan:[t0,tf]分别为起始、终止时间； x0: 状态变量的初始值。 t:求解的时间变量 x:求得的状态变量 方程函数名的编写 编写格式固定 function xdot = 方程函数名(t,x,flag,附加参数) t:时间变量; x:方程的状态变量; xdot:状态变量的导数; flag:运行标志位,系统变量; 附加参数:用逗号分隔; An example 设著名的Lorenz模型的状态方程表示为 若令其初值为 试用Matlab求解该方程的数值解。 参见：lorenzeq，main_lorenzeq ex2: 已知著名的Van der Pol方程 选择状态变量 ,则原方程可以变换成 这里 为一可变参数，使用附加参数的方法对其进行传递。 参见函数：vdp main_vdp Ex3 考虑著名的Rossler微分方程组， 选定a=b=0.2,c=5.7,且x1(0)=x2(0)=x3(0)=0,求解该微分方程。 求职应注意的礼仪 求职时最礼貌的修饰是淡妆 面试时最关键的神情是郑重 无论站还是坐，不能摇动和抖动 对话时目光不能游弋不定 要控制小动作 不要为掩饰紧张情绪而散淡 最优雅的礼仪修养是体现自然 以一种修养面对两种结果 必须首先学会面对的一种结果----被拒绝 仍然感谢这次机会，因为被拒绝是面试后的两种结果之一。 被拒绝是招聘单位对我们综合考虑的结果，因为我们最关心的是自己什么地方与用人要求不一致，而不仅仅是面试中的表现。 不要欺骗自己，说“我本来就不想去”等等。 认真考虑是否有必要再做努力。 必须学会欣然面对的一种结果----被接纳 以具体的形式感谢招聘单位的接纳，如邮件、短信 考虑怎样使自己的知识能力更适应工作需要 把走进工作岗位当作职业生涯的重要的第一步，认真思考如何为以后的发展开好头。 Thank you * * * *