第6章常微分方程祥解.ppt

第六章 习题 ex1(1)(2)(3)(勘误: (2) -x’=2x+3y ) ex2 ex7 ex8 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * * * MATLAB数学实验 第六章 常微分方程 第六章 常微分方程 6.1 预备知识：常微分方程 6.2 解常微分方程的MATLAB指令 6.3 计算实验：Euler法和刚性方程组 6.4 建模实验: 导弹系统的改进 1.微分方程的概念 常微分方程: f (t,y,y’,y’’,…,y(n))=0 微分方程组: 联系一些未知函数的一组微 分方程 线性常微分方程: y(n) + a1 (t) y(n-1) + … + an-1 (t) y’ + an (t) y = b(t) 若ai (t) (i =1, …,n) 与t无关, 称为常系数的 若b(t)=0,称为齐次的 6.1 预备知识：常微分方程 2.初等积分法 分离变量法等 3.常系数线性微分方程 线性常微分方程的解为一个特解和相应 的齐次微分方程通解的叠加。 齐次微分方程的解可用特征根法求得 例1 求x’’+ 0.2 x’+3.92x = 0的通解 解 特征方程为?2 + 0.2? +3.92=0 ? roots([1 0.2 3.92] 求得共轭复根 ?+?i=-0.1?1.9774i, 通解为 x(t) = Ae?t cos(?t) +Be?t sin(?t) 4.初值问题数值解 数值解法:寻求解y(t)在一系列离散节点 t0 t1 … tn tf 上的近似值yk (k=0,1,…n)。 hk = tk+1 －tk 为步长，通常取为常量h 。 其中 高阶常微分方程初值问题可以化为一阶常微分方程组，已给一个n阶方程 y(n) = f (t, y, y’,…,y(n-1)) 设y1 = y, y2 = y’, …, y n= y(n-1), 化为一阶方程组 y0： 表示初值向量y0； t： 表示节点列向量 (t0 , t1 , … , tn)T； y： 数值解矩阵,每一列对应y的一个分量 若无输出参数，则作出图形 初值问题求解? 常用格式[t,y] = ode45(odefun, tspan, y0) odefun：表示f (t, y)的函数句柄或Inline函数 t是标量，y是标量或向量； tspan： 若为[t0, tf], 表示自变量初值t0和终值tf 若为[t0,t1,?, tn], 表示输出节点列向量 6.2解常微分方程MATLAB指令 完整格式 [t,y] = ode45(odefun, tspan, y0,options, p1,p2,?) options——为计算参数(如精度要求)设 置，默认可用空矩阵[]表示； p1,p2,? ——为附加传递参数，这时 odefun的表示为f (t, y, p1,p2, ?) odefun= inline(y-2*t/y,t,y); [t,y]=ode45(odefun,[0,4],1); [t,y] plot(t,y,o-) [t,y]=ode45(odefun,0:1:4,1);[t,y] 例2 解微分方程 y’ = y?2t/y, y(0)=1, 0t4 例3 解微分方程组 0t30 %M函数eg6_3fun.m function f=eg6_3fun(t,x) f(1)=-x(1)^3-x(2); f(2)=x(1)-x(2)^3; f=f(:); clear;[t,x]=ode45(@eg6_3fun,[0 30],[1;0.5]); subplot(1,2,1); plot(t,x(:,1),t,x(:,2),:); subplot(1,2,2); plo