第三讲 典型函数的使用.ppt

第三讲 典型函数的使用 ode45函数（求解常微分方程组） Ordinary Differential Equation(常微分方程） Partial Differential Equations (偏微分方程） fsolve函数（求解非线性方程组） 问题的提出—如何求解下面振动问题？ 如下图所示弹簧—质量—阻尼器系统中，振体质量为m，弹簧的线性系数为k，非线性系数为α，阻尼器的阻尼系数为c。振体受一外力F(t)=X0kf(t)作用。 根据已有知识求解上述问题 列运动方程为 该方程为研究力学系统的一种典型模型，研究范围可以包括洗衣机到车辆的振动问题。 常微分方程 常微分方程的概念 含有自变量，未知函数及其导数（微分或偏导数）的方程称为微分方程，如果未知函数只含有一个变量，称为常微分方程，如果未知函数含有若干个变量，则称为偏微分方程，微分方程中未知函数的导数或偏导数的最高阶次称为微分方程的阶。 常微分方程数值解 实际中，很多问题的数学模型都是微分方程。我们可以研究它们的一些性质。但是，只有极少数特殊的方程有解析解。对于绝大部分的微分方程是没有解析解的。常微分方程作为微分方程的基本类型之一，在自然界与工程界有很广泛的应用。很多问题的数学表述都可以归结为常微分方程的定解问题。很多偏微分方程问题，也可以化为常微分方程问题来近似求解。 常微分方程数值解法 对于一个常微分方程： 通常会有无穷个解。如 因此，我们要加入一个限定条件。通常会在端点处给出，如下面的初值问题： 常微分方程数值解法 常微分方程的解是一个函数，但是，计算机没有办法对函数进行运算。因此，常微分方程的数值解并不是求函数的近似，而是求解函数在某些节点的近似值。 常微分方程数值解法（初值问题的Euler方法） 初值问题的Euler方法 对Euler方法的简单讨论 数值方法，主要研究如何建立差分方程，并研究差分方程的性质。 为了考察数值方法提供的数值解，是否有实用价值，需要知道如下几个结论： ① 步长充分小时，所得到的数值解能否逼近问题的真解；即收敛性问题 ② 误差估计 ③ 产生的舍入误差，在以后的各步计算中，是否会无限制扩大；稳定性问题 n个一阶常微分方程组的标准形式 ode45函数的调用方式 函数名为FunctionName的函数的创建 ode45函数应用举例1 某机械运动系统如下图所示，物体质量为M=1kg,弹簧的虎克系数为常数K=20N/m,摩擦系数为B=5N/m/s。t=0时，施加外力f(t)=30N,试给出系统的运动方程，画出10秒内机械系统位移、速度随时间变化的坐标曲线。 求解过程步骤1 首先根据牛顿运动定理，建立系统的运动方程为 求解过程步骤2 创建函数如下： function s=ForcingFunction(t,w) s=[w(2);-5*w(2)-20*w(1)+30]; 调用形式如下： [tt,yy]=ode45(@ForcingFunction,[0 10],[0 0],[ ]); plot(tt,yy(:,1))%位移曲线 hold on plot(tt,yy(:,2))%速度曲线 ode45函数应用举例2（有p1,p2等参数的传递） 求解单变量二阶常微分方程 求解步骤1—创建函数 function s=ForcingFunction(t,w,x,c) switch c case 1 s=[w(2);-2*x*w(2)-w(1)+1]; case 2 s=[w(2);-2*x*w(2)-w(1)]; case 3 h=sin(pi*t/5).*(t=5); s=[w(2);-2*x*w(2)-w(1)+h]; end 其中 求解步骤2—直接调用ode45函数获得问题的解 第1种情况， 程序如下 [tt,yy]=ode45(ForcingFunction,[0 35],[0 0],[],0.15,1); plot(tt,yy(:,1)) 结果 情况2和3的调用语句 [tt,yy]=ode45(ForcingFunction,[0 35],[1 0],[],0.15,2); plot(tt,yy(:,1)) [tt,yy]=ode45(ForcingFunction,linspace(0,35,200),[0 0],[],0.15,3); plot(tt,yy(:,1)) 非线性方程的数值求解 n元非线性方程组的标准形式 fsolve的调用格式 fsolve(FunctionNa