8---微分方程(组)求解.pptVIP

数学实验 微分方程(组)求解 微分方程的解 代入微分方程后能使方程两端恒等的函数称为微分方程的解。 通解：若微分方程的解中所含任意常数的个数恰好等于微分方程的阶数，则称此解为微分方程的通解。 特解：通解中的任意常数得到确定后，称其为微分方程的特解。 数值解：求得的特解函数是用“列表法”表示的函数（非解析式）。 符号解(解析解)：解函数以解析式表示。 定解条件 定解条件：确定通解中的常数，使其成为特解所需的额外附加条件称为定解条件。 初始条件：定解条件由系统在某一瞬间所处的状态给出时（t0时刻的函数值及各阶导数值），称其为初始条件。 定解问题：带定解条件的微分方程问题。 初值问题：带初始条件的微分方程问题。 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Matlab函数数值求解 Matlab提供的ODE求解器 参数说明 数值求解举例 数值求解举例 再说定义函数文件 再说定义函数文件 举例说明 举例说明（单个方程） 举例说明（方程组） 思考 说明 例 高阶常微分方程 思考 对于高阶微分方程组的数值求解。应用如何处理？ 微分代数方程(DAE) DAE 举例 下面以实例说明，求解方程组的数值解 (2)编写Matlab代码 x0=[0.8;0.1;0.1];%初值 [t,x]=ode15s(odefun,[0 20],x0,options);%这里好像不能使用ode45 plot(t,x);legend(x1(t),x2(t),x3(t)) * * 主要内容 解析解 数值解 微分代数方程求解 含有未知函数的导数或微分的方程称微分方程。 微分方程中出现的各阶导数（微分）的最高阶数称为微分方程的阶。 未知函数为一元函数的微分方程称常微分方程。 未知函数为多元函数的微分方程称偏微分方程。 微分方程相关概念 常微分方程的一般形式为: dsolve 的使用 y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) 其中 y 为输出， eq1、eq2、...为微分方程，cond1、cond2、...为初值条件，v 为自变量。 例 1：求微分方程 的通解，并验证。 y=dsolve(Dy+2*x*y=x*exp(-x^2),x) syms x; diff(y)+2*x*y - x*exp(-x^2) 几点说明 如果省略初值条件，则表示求通解； 如果省略自变量，则默认自变量为 t dsolve(Dy=2*x,x); % dy/dx = 2x dsolve(Dy=2*x); % dy/dt = 2x 若找不到解析解，则返回其积分形式。 微分方程中用 D 表示对 自变量 的导数，如： Dy y； D2y y； D3y y 例 2：求微分方程 在初值条件 下的特解，并画出解函数的图形。 y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) ezplot(y); 例3：求微分方程组 在初值条件 下的特解，并画出解函数的图形。 [x,y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=0, ... x(0)=1, y(0)=0, t) ezplot(x,y,[0,1.3]); 注：解微分方程组时，如果所给的输出个数与方程个数相同，则方程组的解按词典顺序输出；如果只给一个输出，则输出的是一个包含解的结构（structure）类型的数据。 例： [x,y]=dsolve(Dx+5*x=0,Dy-3*y=0, ... x(0)=1, y(0)=1,t) r = dsolve(Dx+5*x=0,Dy-3*y=0, ... x(0)=1, y(0)=1,t) 这里返回的 r 是一个 结构类型 的数据 r.x %查看解函数 x(t) r.y %查看解函数 y(t) 只有很少一部分微分方程（组）能求出解析解。大部分微分方程（组）只能利用数值方法求数值解。 dsolve的输出个数只能为一个 或 与方程个数相等。 [T,Y] = solver(odefun,tspan,y0) 其中 y0 为初值条件，tspan为求解区间；Matlab在数值求解时自动对求解区间进行分割，T (向量) 中返回的是分割点的值(自变量)，Y (向量) 中返回的是解函数在这些分割点上的