优·matlab第四次课.pptVIP

* * * * 常微分方程的求解 统计作图 多项式拟合 随机试验 设微分方程初值问题： 常微分方程的求解 分为：符号解法和数值解法 符号解法： 命令形式1： dsolve(’eqution’,’var’) 命令形式2： dsolve(’eqution’ , ’cond1,cond2,…’ , ’var’) P158 初始条件 符号解法 dsolve(Dy=y^2,x) 自变量 大写 ?dsolve(Dy=y^2,x) ans = -1/(x-C1) 求解 ?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是初始状态列向量。 常 微 分 方 程 的 数 值 解 法 高阶微分方程需转化一阶微分方程组 注：只适用于一阶微分方程 例7-45：用数值积分的方法求解微分方程： 设初始时间 ；终止时间 初始条件 分析：求解 令： （化为一阶微分方程）即原微分方程化为： 写成矩阵形式为