数学实验第三讲.pptVIP

3、Rossler微分方程组： 当固定参数b=2,c=4时，试讨论随参数a由小到大变化（如a∈(0,0.65])而方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状？ 数学实验之三 微分方程解的形式 ①　解析解 y = f(x) ②　数值解 (xi, yi) ③　图形解 x y o ①简单的微分方程。 ②复杂、大型的微分方程。 一阶微分方程: 获取解析解的方法归类: ① 分离变量法；如 dy/dx = x*y； ② 齐次方程的变换法； 如 dy/dx = f (y/x) ③ 线性方程的常数变易法或公式法. …… 解析解 MATLAB软件实现 解析解 dsolve(eqn1,eqn2, …, c1 ,…, var1 ,…) 微分方程组 初值条件 变量组 注意：① y Dy，y D2y ② 自变量名可以省略，默认变量名‘t’。 例① 输入：y=dsolve (Dy=1+y^2) y1=dsolve(Dy=1+y^2,y(0)=1,x) 输出：y= tan(t-C1) （通解，一簇曲线） y1= tan(x+1/4*pi)（特解，一条曲线） 例② 常系数的二阶微分方程 y=dsolve(D2y-2*Dy-3*y=0,x) y=dsolve(D2y-2*Dy-3*y=0,y(0)=1,Dy(0)=0,x) 输入 ： x=dsolve(D2x-(1-x^2)*Dx+x=0, x(0)=3,Dx(0)=0) 上述两例的计算结果怎样？由此得出什么结论？ 例③ 非常系数的二阶微分方程 例③无解析表达式！ x=dsolve((Dx)^2+x^2=1,x(0)=0) 例④ 非线性微分方程 x = [ sin(t)] [ -sin(t)] 若欲求解的某个数值解，如何求解？ t=pi/2; eval(x) 输入： [x,y]=dsolve(Dx=3*x+4*y,Dy=-4*x+3*y) [x,y]=dsolve(Dx=3*x+4*y,Dy=-4*x+3*y, x(0)=0,y(0)=1) 例④ 输出： （ li3.m） 数值解 1、欧拉法 2、龙格—库塔法 数值求解思想：(变量离散化) 引入自变量点列{xn} → {yn}， 在x0?x1?x2?…?xn?…上求y(xn)的近似值yn.通常取等步长 h，即xn = x0+ n×h，或 xn = xn-1+ h，（n=1,2,…)。 研究常微分方程的数值解法是十分必要的。 1) 向前欧拉公式: （y’= f (x, y) ） y (xn+1) ? y(xn) + h f(xn, y(xn)) (迭代式) yn+1 ? yn + h f(xn, yn) (近似式) 特点：f(x,y)取值于区间[xn, xn+1]的左端点. 1、欧拉方法 在小区间[xn, xn+1]上用差商代替微商(近似), yn+1 ? yn + h f(xn +1, yn +1) 特点：① f(x,y)取值于区间[xn, xn+1]的右端点. ② 非线性方程, 称‘隐式公式’。 yn+1 = yn + h f(xn, yn) 2) 向后欧拉公式 方法：迭代（ y’= f (x, y) ） x=[];y=[]; x(1)=x0; y(1)=y0; for n=1:k x(n+1)=x(n)+n*h; y(n+1) = y(n) + h *f(x(n), y(n)); (向前） end 例 1 观察向前欧拉、向后欧拉算法计算情况。与精确解进行比较。误差有多大？ 解：1） 解析解： y = x + e-x y=dsolve(Dy=-y+x+1,y(0)=1,x) 2） 向前欧拉法： yn+1 = yn + h(-yn + xn+ 1) = (1-h) yn + h xn+ h 3）向后欧拉法： yn+1 = yn + h(- yn +1+xn +1+1) 转化 yn+1 = (yn + h xn+1+ h )/(1+h) y’ = f(x,y) = -y + x +1; x1(1)=0;y1(1)=1;y2(1)=1;h=0.1;(died.m) for k=1:10 x1(k+1)=x1(k)+h; y1(k+1)=(1-h)*