微分方程的定.docVIP

微分方程的定义：凡含有未知函数的导数（偏导数）或微分的方程叫微分方程。（记号：在表达微分方程时，用字母D表示求微分，D2，D3等表示求高阶微分，任何D后所跟的字母为因变量，自变量可以指定或有系统规则选定，例如用表示微分方程。 ：）。微分方程中出现的未知函数最高阶导数的阶数，称为微分方程的阶。定义式如下： F(x, y, y￠, ...., y(n)) = 0 常微分方程：在很多学科领域内有着重要的应用，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。这些问题都可以化为求常微分方程的解，或者化为研究解的性质的问题。 常微分方程实例 下下列方程都是微分方程 (其中 y, v, q 均为未知函数). (1) y￠= kx, k 为常数； (2) ( y - 2xy) dx + x2 dy = 0； (3) mv￠(t) = mg - kv(t)； 一阶微分方程的形式及解法 （1）一阶微分方程的普遍形式 一般形式：F（x,y,y）=0 标准形式：y=f(x,y) 微分方程的解： 任何代入微分方程后使其成为恒等式的函数，都叫做该方程的解.若微分方程的解中含有任意常数的个数与方程的阶数相同，且任意常数之间不能合并，则称此解为该方程的通解(或一般解).当通解中的各任意常数都取特定值时所得到的解，称为方程的特解。 一般地说，n 阶微分方程的解含有 n个任意常数。也就是说，微分方程的解中含有任意常数的个数和方程的阶数相同，这种解叫做微分方程的通解。通解构成一个函数族。 如果根据实际问题要求出其中满足某种指定条件的解来，那么求这种解的问题叫做定解问题，对于一个常微分方程的满足定解条件的解叫做特解。对于高阶微分方程可以引入新的未知函数，把它化为多个一阶微分方程组。 微分方程求解： 1.可分离变量：形如y’=f(x)*g(y)解法 两边同时积分为 2.齐次方程：形如y’=令u=y/x有u+u’x=y’代入原方程有u+u’x=f(u)[可分离变量] 3.全微分方程：形如其中 4.线性微分方程 微分方程求解方法：解析解、数值解法、定性理论 解析解： Matalb中求微分方程解析解的函数是dsolve，其调用格式为： dsolve(‘eq1,eq2,…’, ‘cond1,cond2,…’, ‘v’)该函数求解常微分方程组eq1,…,eqn在初值条件cond1,…,condn下的特解，若不给出初值条件，则求方程组的通解，v给出指定的自变量，如果不给出，默认的自变量为。 例如，求常微分方程的通解。 输入： dsolve(Dy=2) 其默认的独立变量为。 例如，求常微分方程的通解。 输入：dsolve(Dy=1/(x+y),x) 输出结果为： ans = -lambertw(-C1*exp(-x-1))-x-1 这里Y=lambertw(X)表示：Y*exp(Y)=X。 例如，求常微分方程的通解。 输入： dsolve(D2y-2*Dy+5*y=exp(x)*sin(2*x),x) ans = exp(x)*sin(2*x)*C2+exp(x)*cos(2*x)*C1-1/4*exp(x)*cos(2*x)*x 例如，求常微分方程满足初始条件的特解。 输入： dsolve(D2y-y=4*x*exp(x),y(0)=0,Dy(0)=1,x) ans = exp(x)-exp(-x)+(-1+x)*x*exp(x) 例如，求常微分方程组的通解。 输入：[x,y]=dsolve(Dx=y+1,Dy=x+1,t) x =exp(t)*C2-exp(-t)*C1-1 y =exp(t)*C2+exp(-t)*C1-1 数值解： 微分方程数值解法的命令格式：[t,x]= solver (xfun,[t0 tf],y0,tol) (1) 在解n个未知函数的方程组时，x0和x均为n维向量，m-文件中的待解方程组应以x的分量形式写成； (2) 使用Matlab软件求数值解时，高阶微分方程必须等价地变换成一阶微分方程组 例如，求微分方程在区间上的数值解。 解： 建立m文件example0707.m如下： function dy=example0707(t,y) dy=t; 输入命令： [t, y]=ode45(example0707,[-1,1],1); plot(t,y,*) 一阶微分方程： y’=f(x,y) (1) 对称形式：p(x,y)dx+q(x,