第八章 微分方程1.pptVIP

第八章 微分方程 8.1 引例 8.2 微分方程的基本概念 微分方程的解 8.3 常微分方程数值解法介绍 8.4 微分方程求解的MATLAB实现 例8.1 求微分方程dy/dx=1+y2的通解 2. 求微分方程的特解 3. 求解微分方程组 高阶微分方程化为一阶微分方程组的方法 对任意高阶微分方程 x(n)=f(t,x,x, ···,x(n-1)) 一阶微分方程组的向量形式表示 一阶微分方程组的向量形式表示 一阶线性微分方程组的向量形式表示 求一阶微分方程组数值解的MATLAB命令ode 例8.5 求解线性微分方程组X=AX，其中 例8.6 求解Van der Pol方程x+(x2-1)x+x=0 积分曲线 由y1=x可知，原方程的解是输出矩阵Y的第一列 Van der Pol方程 x+(x2-1)x+ x = 0满足初始条件x(0)=0,x(0)=1的解的轨迹可由下面的命令给出 plot(T,Y(:,1)) 向量微分方程组dy/dt=f(t,y)的ode输出 [T,Y] = ode45 (fun, Ts, Y0, options ) x1 = x2 x2 = 4x1+3x2-4x3 x3 = x1+2x2+x3 解：建立函数M文件exam1： function dx=exam1(t,x) dx=zeros(3,1); dx(1)=x(2); dx(2)=4*x(1)+3*x(2)-4*x(3); dx(3)=x(1)+2*x(2)+x(3); 取初始条件x(0)=[1;2;3], Ts=[0,20],利用ODE45命令求解如下 [T,Y]=ode45(exam1, [0,20], [1;2;3]); 利用命令plot(T,Y(:,1))可得积分曲线. 解：令y1 = x, y2 = y1=x, 则方程化为方程组 y1 =y2 , y2 = (1- y12)y2 - y1 建立函数M文件vdp： function dy=vdp(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=(1-y(1)^2)*y(2)-y(1); 取初始条件y(0)=[0;1], Ts=[0,20],利用ODE45命令求解如下 Ts =[0,20]; y0=[0;1]; [T,Y]=ode45(vdp,Ts,y0); * * 厘 厘 解: 设所求曲线方程为y=y(x), 则 又y(x)应满足条件x=1时y=2，因此可解得C=1 两端对x积分，得 例1 求xOy平面上过点(1,2)的一条曲线方程，使其在任意点处的斜率为2x。 故所求曲线方程为y=x2+1。 例2 列车上以20米/秒的速度行驶，当制动时，列车获得加速度-0.4米/秒2，问开始制动后多长时间能停住，以及列车在这段时间行驶了多少路程？ 等式两边同时积分，得 解: 设列车开始制动后t秒内行驶了s米，则 再积分得 利用条件t=0,s=0,ds/dt=20，可以解得C1=20, C2=0,因此 利用条件列车停住时速度为零，由 解得:s=500(米) 代入 及C1=20可得: t=50(秒) 微分方程: 例如 含有未知函数的导数或微分的方程 常微分方程(ODE) (Ordinary Differential Equations) 仅含一个自变量的微分方程. 微分方程的阶: 微分方程中出现的未知函数的最高阶导数的阶数. 线性与非线性微分方程:方程中关于未知函数及其各阶导数均是一次的,则称为线性微分方程. 线性微分方程 非线性微分方程 微分方程的解:代入微分方程能使方程成为恒等式的函数称为微分方程的解。 微分方程解的分类： (1)通解: 微分方程的解中含有任意常数,且任意常数的个数与微分方程的阶数相同。 (2)特解: 确定了通解中任意常数以后的解。 解的图象(积分曲线): 微分方程的解曲线。 初始条件: 用来确定任意常数的条件。 考虑引例1，对于微分方程 函数y=x2+C是其通解，y=x2+1是其满足初始条件x0=1，y(x0)=2的特解。 微分方程的(部分)积分曲线(0≤C≤2)见右图，其中红色曲线为满足初始条件的特解曲线。 考虑引例2，对于微分方程 函数s=-0.2t2+C1t+C2是其通解, s=-0.2t2+20t是其满足条件t0=0，s(t0)=0, ds/dt (t0)=20的特解. 满足初始条件的积分(特解)曲线见下图。 路程曲线 速度曲线 ?微分方程是研究函数变化规律的有力工具，在科技、工程技术、经济管理以及生态、环境、人口、交通等各个领域中有着十分广泛的应用。 ?很多实际问题的数学模型都可以用微分方程来表示，但建立微分方程只是解决问题的第一步，通常在此基础上还需要求出微分方程的解并对实际问题进行