《控制系统计算机仿真课件》ch2.1 2.2 连续系统仿真.pptVIP

图2.18 (exam2_7.mdl) * 然后编制M函数内容如下： % saturation_ zone function function [uo] = saturation_ zone (ui) if ui=1 uo=1; elseif ui= -1 uo= -1; else uo=ui; end 启动仿真后，得到如右图所示仿真结果。 * 从图中可以看出，解析解单调下降并迅速收敛到0。 当 时，欧拉法和RK4法的解曲线均发散，数值积分算法的解是错误的。 当 时，欧拉法的解曲线仍然发散，对应的解是错误的；RK4法的解曲线单调下降并收敛到0，对应的解是正确的。 当 时，欧拉法和RK4法的解均收敛到0（虽然欧拉法的解曲线是振荡收敛的）。如果只要求得到 处的 的解，则两种数值积分算法的解都可以认为是正确的。 * 欧拉法 RK4法 解析解 为什么会出现上述情况呢？这是因为数值积分算法只是一种近似积分方法，在反复的递推计算中会引进误差。如果误差的累积越来越大，将使计算出现不稳定，从而得到错误的结果。 系统的稳定性与计算稳定性是两个不同的概念。前者用原系统的微分方程、传递函数来讨论，后者用逼