matlab第四次课42433.pptVIP

* * 常微分方程的求解 数据的统计分析 曲线拟合 随机试验 设微分方程初值问题： 命令形式1： dsolve（’eqution’,’var’） 命令形式2： dsolve（’eqution’ , ’cond1,cond2,…’ , ’var’） 常微分方程的求解 初始条件 符号解法 dsolve(Dy=y/(x^2),x) 自变量 大写 ?dsolve(Dy=y/(x^2),x) ans = C1*exp(-1/x) 求解 ?dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) 大写 对应求导阶数 condition ans =-1/3*x^3+125/468+31/468*x^4 求解常微分方程 [t,y]=ode23(‘fun’,tspan,y0) 需要把微分方程化为一阶微分方程组 tspan=[t0 tf]表示积分的起始值和终止值； fun是定义函数的文件名； [t,y]=ode45(‘fun’,tspan,y0) y0是初始状态列向量。 常 微 分 方 程 的 数 值 解 法 例：用数值积分的方法求解微分方程： 设初始时间 ；终止时间 初始条件 分析：求解 令： （化为一阶微分方程）即原微分方程化为： 写成矩阵形式为 （化为一阶微分方程） xdot= 放入函数exf.m中 u function xdot=exf(t,x