缉私艇追走私船模型实验.pptxVIP

会计学;用MATLAB软件求微分方程解析解、数值解理论. 编程求解微分方程数值解. 微分方程模型实验:缉私艇追赶走私船. ;1. 微分方程的解; 用matlab软件求微分方程的解析解;例1 (1);微分方程数值解的计算 欧拉方法、欧拉两步法、改进的欧拉方法;所以，过初始点A0(x0,y0)，做切线;这样依次类推得到上述微分方程的解在x=xn+1处 函数值的近似计算公式，;特别地，如果节点之间是等分的，则;将上两式算术平均即得梯形公式：;2）欧拉两步法;3）改进的欧拉公式;步骤2 对预测值再用梯形公式校正一次得yn+1，称???校正值;再解释：;编程计算微分方程数值解;解：（1）直接在命令窗执行命令 dsolve(Dy=3/x*y+x^3*(exp(x)+cos(x))-2*x, y(pi)=(exp(pi)+2/pi)*pi^3,x) ans= x^3*exp(x)+x^3*sin(x)+2*x^2 ;clc;clf; szy_eu=[]; y=(exp(pi)+2/pi)*pi^3; for x=pi:0.1:2*pi y=y+0.1*(3/x*y+x^3*(exp(x)+cos(x))-2*x); szy_eu=[szy_eu,y]; end szy_eu t=pi:0.1:2*pi; f=t.^3.*exp(t)+t.^3.*sin(t)+2*t.^2; plot