1常微分方程常见解法.pptVIP

一、 向量场;曲线（称为积分曲线），且 ;它所确定的向量场中的一条曲线，该曲线所经过的;例1.3.1 在区域 ;Maple指令：;回车后Maple就在 ;所谓图解法就是不用微分方程解的具体表达式，直接根据右端函数的结构和向量场作出积分曲线的大致图形。 图解法只是定性的，只反映积分曲线的一部分主要特征。 该方法的思想却十分重要。因为能够用初等方法求解的方程极少，用图解法来分析积分曲线的性态对了解该方程所反映的实际现象的变化规律就有很重要的指导意义。;三、一阶常微分方程的解法;初值问题 ;Bernoulli方程;例 湖泊的污染;因此有;积分得;当 , ;齐次函数: 函数;可化为齐次方程的方程;此时二元方程组;(2) 若;4. 对特殊方程;例 求方程;还原后得原方程通解为;例:雪球融化问题;引入新常数 ;中连续且有连续的一阶偏导数，则 ;例：验证方程;由公式(2.3.4)得: ;(2)偏积分法;而 ;例：验证方程;方程的通解为： ;四、微分方程的近似解法;能得到解析解的方程：; 对初始值问题 构造迭代序列 该序列一致收敛到解，故迭代一定次数后就可以作为一个近似;……;迭代的误差 (|x|n);例：求初始值问题解的迭代序列的前三项; 利用Maple软件可以求出更多的项 y[0]:=1; for j from 1 to 4 do y[j]:=1+int(x^2+y[j-1]^2,x=0..x); end do;;近似程度的显示;Taylor级数法;例： 用Tailor级数法求初始值 问题的近似解。 解：计算解函数在x=0点的函数值和各阶导数值得 所以，该初始值问题的近似解为; 用Maple处理 restart: ode1:=diff(y(x),x)-x^2-y(x)^2=0; for j from 1 to 7 do Order:=j*2: dsolve({ode1,y(0)=1},y(x),type=series); sol[j]:=rhs(%): end do: for j from 1 to 7 do y[j]:=sol[j]; end do;;运行后的结果;用Maple 处理并用图形显示 restart: ode1:=diff(y(x),x)-x^2-y(x)^2=0; for j from 1 to 7 do Order:=j*2; convert(dsolve({ode1,y(0)=1},y(x), type=series), polynom ): sol[j]:=rhs(%): end do: plot({sol[1],sol[2],sol[3],sol[4],sol[5],sol[6],sol[7]},x=-1..1);;在区间[-1,1]的近似情况;在区间[-0.5,0.5]内的近似情况;设初始值问题的解在 可以展开为幂级数 代入初始条件 方程后得 展开后比较两端同次幂的系数确定;例：用待定系数法求 ; 计算出解函数 在一系列节点 处的近似值, 节点间距;Euler折线法;;利用Taylor公式;改进的Euler折线法; 利用计算机编程 给出步长和初始值 循环计算各点上函数的近似值 显示结果 例 求初始值问题的数值解 ;printlev1:=0: h:=0.1: x[0]:=0: y[0]:=0.5: z[0]:=0.5: f1:=(x,y)-1+(y-x)^2; f2:=(x,y)-2*(x-y)+2*(y-x)*(1+(y-x)^2); for n from 0 to 9 do x[n+1]:=h*(n+1); y[n+1]:=y[n]+h*f1(x[n],y[n]); z[n+1]:=z[n]+h*f1(x[n],z[n])+h^2*f2(x[n],z[n])/2; u[n+1]:=x[n+1]+1/(2-x[n+1]); print (x[n+1],y[n+1],z[n+1],u[n+1]); od: 可以改变步长和增加分点来观察计算精度的变化情况;;对于常微分方程的边值问题;----------(2);低阶Runge-Kutta方法;(由(4)式);-----------(5);高阶Runge-Kutta方法;令;取;还可构造四阶(经典)Runge-Kutta方法;例 求初始值问题的数值解 ;printlev1:=0: h:=0.1: x[0]:=0: y[0]:=0.5: f:=(x,y)-1+(y-x)^