第9章 常微分方程初值问题的数值解.pptVIP

内蒙古大学自动化系 第9章 常微分方程初值问题的数值解 《数值计算与MATLAB》 内蒙古大学自动化系 * 内蒙古大学自动化系 * 数值计算与MATLAB * 第9章 常微分方程初值问题的数值解 G(x,y,y,…,y (m) ) = 0 或 y (m) = f(x,y,y…, y (m-1) ) 通解：含有待定常数 特解：不含待定常数 （给出初始条件） 解析解：解析表达式 数值解：表格法或图示法 * 1 求解常微分方程的MATLAB符号法 1.1 常微分方程的MATLAB符号表示法 （1）用“Dmy”表示函数y=f(x)的m阶导数y(m)=f(m)(x)。 （2）初始条件的写法和上述规定完全一致。 （3）不特别界定时，通常默认小写字母“t”为函数的自变量。 例如： Dy+2*x*y=1，y(0)=0 * 1.2 求解常微分方程的符号法指令dsolve [y1, y2，…, y12]=dsolve(a1, a2, …, a12) （1）每个输入参数a1，a2，…, a12； （2）“初始条件” ； （3）“界定的自变量” ； （4）参量ai可以多达12个； （5）输出参量。 举例。 * 2 常微分方程数值解的基本原理 2.1 求常微分方程数值解的基本原理 1. 高阶微分方程与一阶微分方程组 y1(x0)=y(x0)=y0，y2