第四章连续系统的离散化方法.pptVIP

第四章 连续系统的离散化方法 4.1 常微分方程的数值解法 其一般公式为 解：因欧拉法的递推公式为 4.1.3 龙格－库塔法 将 2、 其意义如下：用欧拉法以斜率先求取一点， 再由此点求得另一斜率 四阶龙格－库塔法的计算公式为： 其中， 四阶/五阶龙格－库塔法系数表 4.1.4 微分方程数值解的MATLAB实现 1、 微分方程组中的方程个数必须等于初始条件数，这是求微分方程特解所必须的条件。 6、输入参数options表示选项参数（包括tol，trace），可缺省，即取默认值，tol是控制结果精度的选项对ode23( )函数取 ，对ode45( )函数取 。trace为输出形式控制变量，如果trace不为0，则 例 4－1 求解常微分方程 首先编制M文件，并且函数名和M文件名相同。 function xdot=wffc_1(t,x) %定义t,x,xdot和文件名 xdot=[0 1;-1 0]*x+[0;1]*(2+t^2/pi); %状态方程的表示形式 编制M文件，并且函数名和M文件名相同。 function xdot=vdpl(t,x) xdot=zeros(2,1); ％赋初值，并规定向量的维数。 xdot(1)=x(2);