matlab教学第二章微分方程讲义.pptVIP

1. 微分方程 例2.6 求解一阶微分方程 ??? (1) 求通解 输入： dsolve(Dy=1+y^2) 输出： ans = tan(t+C1) （2）求特解 输入： dsolve(Dy=1+y^2,y(0)=1,x) 指定初值为1，自变量为x 输出： ans = tan(x+1/4*pi) * * 例2.7 求解二阶微分方程 输入: dsolve(D2y+(1/x)*Dy+(1-(1/2)^2/x^2)*y=0,y(pi/2)=2,Dy(pi/2)=-2/pi,x) ans =2^(1/2)*pi^(1/2)/x^(1/2)*sin(x) 化简输出结果，输入： pretty(ans) 1/2 1/2 2 pi sin(x) ----------------- 1/2 x 即： * * 2．微分方程组 例2.8 求解???? df/dx=3f+4g;? dg/dx=-4f+3g。 （1）通解:?? [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g) f = exp(3*t)*(C1*sin(4*t)+C2*cos(4*t)) g = exp(3*t)*(C1*cos(4*t)-C2*sin(4*t)) 特解： [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g(0)=1) f = exp(3*t)*sin(4*t) g = exp(3*t)*cos(4*t) * * 2.4.2 数值解 在微分方程(组)难以获得解析解的情况下，可以用Matlab方便地求出数值解。格式为: [t,y] = ode23(‘F’,ts,y0)， [t,y] = ode45(‘F’,ts,y0) 注意： 微分方程的形式：y = F(t, y)，t为自变量，y为因变量（可以是多个，如微分方程组）； [t, y]为输出矩阵，分别表示自变量和因变量的取值； F代表微分方程组的函数名（m文件，必须返回一个列向量）； ts的取法有几种，（1）ts=[t0, tf] 表示自变量的取值范围，（2）ts=[t0,t1,t2,…,tf],则输出在指定时刻t0,t1,t2,…,tf处给出，（3）ts=t0:k:tf,则输出在区间[t0,tf]的等分点给出； y0为初值条件。 * * 第二章 微分方程模型 * * 本章学习目的： 学习微分方程模型的建立、求解方法、分析结果及解决实际问题的全过程。 熟练掌握使用MATLAB软件的函数求微分方程的解析解、数值解和图形解 。 * * 2.1 引例 对于圆柱形状容器壁上的容积刻度，可以利用圆柱体体积公式： ，其中容器的直径D为常数，体积V与相对于容器底部的任意高度H成正比，因此在容器壁上可以方便地标出容积刻度。 而对于几何形状不规则的容器，比如“倒葫芦形状”的容器壁上如何标出容积刻度呢？ x * * 建立坐标系，由微元法分析可知： 其中x表示高度，直径是高度的函数，记为D(x) 可得微分方程： 如果该方程中的函数D(x)无解析表达式，只给出D(x)的部分测试数据，如何求解此微分方程呢？ * * h=0.2; d=[0.04,0.11,0.26,0.56,1.04,1.17]; x(1)=0;v(1)=0; for k=1:5 x(k+1)=x(k)+h; v(k+1)=v(k)+(h/2)*(pi/4)*(d(k)^2+d(k+1)^2); end x=x(1:6),v=v(1:6), plot(x,v) x = Columns 1 through 5 0 0.2000 0.4000 0.6000 0.8000 Column 6 1.0000 v = Columns 1 through 5 0 0.0011 0.0073 0.0373 0.1469 Column 6 0.3393 * * 2.2 微分方程模型的建立 在工程实际问题中，“改变”、“变化”、“增加”、“减少”等关键词提示我们注意什么量在变化，关键词“速率”、“增长”、“衰变”、“边际的”等常涉及到导数。 我们熟悉的速度公式： 就是一个简单的一阶微分方程