常微分方程课程设计课件（完整版）.ppt

Euler 折线法源程序 clear f=sym(y+2*x/y^2); a=0; b=2; h=0.4; n=(b-a)/h+1; % n=(b-a)/h; x=0; y=1; szj=[x,y]; for i=1:n-1 % i=1:n y=y+h*subs(f,{‘x’,‘y’},{x,y});%第一次代初值 x=x+h; szj=[szj;x,y]; end szj plot(szj(:,1),szj(:,2),or-) 最新.课件 * Euler折线法举例（续） 解析解： 解析解 近似解 最新.课件 * Runge-Kutta 方法 为了减小误差，可采用以下方法： 让步长 h 取得更小一些； 改用具有较高精度的数值方法： 龙格-库塔方法 Runge-Kutta (龙格-库塔) 方法 是一类求解常微分方程的数值方法 有多种不同的迭代格式 最新.课件 * Runge-Kutta 方法 用得较多的是 四阶R-K方法（教材第 79 页） 其中 最新.课件 * 四阶 R-K 方法源程序 clear; f=sym(y+2*x/y^2); a=0; b=2; h=0.4; n=(b-a)/h+1; % n=(b-a)/h; x=0; y=1; szj=[x,y]; for i=1:n-1 % i=1:n l1=subs(f,{x,y},{x,y}); l2=subs(f,{x,y},{x+h/2,y+l1*h/2}); l3=subs(f,{x,y},{x+h/2,y+l2*h/2}); l4=subs(f,{x,y},{x+h,y+l3*h}); y=y+h*(l1+2*l2+2*l3+l4)/6; x=x+h; szj=[szj;x,y]; end plot(szj(:,1),szj(:,2), dg-) 最新.课件 * Runge-Kutta 方法 最新.课件 * Euler 法与 R-K法误差比较 最新.课件 * Matlab 解初值问题 用 Maltab自带函数 解初值问题 求解析解：dsolve 求数值解：  ode45、ode23、 ode113、ode23t、ode15s、 ode23s、ode23tb 最新.课件 * dsolve 求解析解 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) y = 1/2*exp(-x^2)*x^2+exp(-x^2)*C1 最新.课件 * dsolve 的使用 几点说明 如果省略初值条件，则表示求通解； 如果省略自变量，则默认自变量为 t dsolve(Dy=2*x,x); ％ dy/dx = 2x dsolve(Dy=2*x); ％ dy/dt = 2x 若找不到解析解，则返回其积分形式。 微分方程中用 D 表示对 自变量 的导数，如： Dy y； D2y y； D3y y 最新.课件 * dsolve 举例 例 2：求微分方程 在初值条件 下的特解，并画出解函数的图形。 y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) ezplot(y); 最新.课件 * dsolve 举例 例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]); 最新.课件 * dsolve 举例 例： [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 是一个 结构类型