准滑动模态控制matlab仿真实例

准滑动模态控制准滑动模态控制 2 8 1 准滑动模态控制准滑动模态控制 在滑动模态控制系统中 如果控制结构的切换具有理想的开关特性 则能在 切换面上形成理想的滑动模态 这是一种光滑的运动 渐进趋近于原点 但在 实际工程中 由于存在时间上的延迟和空间上的滞后等原因 使得滑动模态呈 抖振形式 在光滑的滑动上叠加了抖振 理想的滑动模态是不存在的 现实中 的滑动模态控制均伴随有抖振 抖振问题是影响滑动模态看控制广泛应用的主 要障碍 所谓准滑动模态 是指系统的运动轨迹被限制在理想滑动模态的某一领 域内的模态 从相轨迹方面来说 具有理想滑动模态的控制是使一定范围内的 状态点均被吸引至切换面 而准滑动模态控制则是使一定范围内的状态点均被 吸引至切换面的某一领域内 通常称此领域为滑动模态切换面的边界层 在边界层内 准滑动模态不要求满足滑动模态的存在条件 因此准滑动模 态不要求在切换面上进行控制结构的切换 它可以在边界层上进行结构变换的 控制系统 也可以根本不进行结构变换的连续状态反馈控制系统 准滑动模态 控制在实现上的这种差别 使它从根本上避免或削弱了抖振 从而在实际中得 到了广泛的应用 在连续系统中 常用的准滑动模态控制有以下两种方法 1 用饱和函数代替理想滑动模态中的符号函数 sat ssgn s 2 46 1 1 s sat skss s 1 k 其中称为 边界层 饱和函数如图 2 26 所示 饱和函数的本质 sat s 为 在边界层外 采用切换控制 在边界之内 采用线性化反馈控制 2 将继电特性连续化 用连续函数取代 s sgn s 2 47 s s s 式中是很小的正常数 2 8 2 仿真实例仿真实例 对象为二阶传递函数 2 p b Gs sas 其中 25 133ab 将传递函数描述为位置状态方程的形式 xxu A 其中 01 0a A 0 b 设位置指令信号为 则r 1 2 1 2 sc rxdrx sc rxrxslaw 采用指数趋近律 sgn slawsks 将代入上式 得控制律为 2 2 xaxbu 1 2 2 uc rxrslawax b 系统的初始状态为 在仿真程序中 指令为正弦信号 0 50 为指数趋近律 为采用饱和函数的指数趋近律 为采用继电1M 2M 3M 特性的指数趋近律 取 进行如下几种仿真 51C 15 30k 1 采用指数趋近律 取 正弦跟踪结果如图所示 1M 00 511 522 533 544 55 0 8 0 6 0 4 0 2 0 0 2 0 4 0 6 time s r yout 正弦跟踪 M 1 00 511 522 533 544 55 0 1 0 0 1 0 2 0 3 0 4 0 5 0 6 time s error 误差曲线 M 1 0 100 10 20 30 40 50 6 3 2 1 0 1 2 3 4 e de 相轨迹 M 1 00 511 522 533 544 55 1 0 1 2 3 4 5 6 time s s 切换函数 s M 1 00 511 522 533 544 55 0 4 0 2 0 0 2 0 4 0 6 0 8 1 1 2 1 4 1 6 time s u 控制器输出 M 1 2 验证准滑动模态控制 为采用饱和函数的方法 取 仿真结2M 5 0 果如图所示 00 511 522 533 544 55 0 8 0 6 0 4 0 2 0 0 2 0 4 0 6 time s r yout 正弦信号跟踪 M 2 00 511 522 533 544 55 0 1 0 0 1 0 2 0 3 0 4 0 5 0 6 time s error 误差曲线 M 2 0 100 10 20 30 40 50 6 3 2 1 0 1 2 3 4 e de 相轨迹 M 2 00 511 522 533 544 55 1 0 1 2 3 4 5 6 time s s 切换函数 s M 2 00 511 522 533 544 55 1 0 5 0 0 5 1 1 5 2 time s u 控制器输出 M 2 3 为采用继电特性的方法 取 仿真结果如图所示3M 0 05 00 511 522 533 544 55 0 8 0 6 0 4 0 2 0 0 2 0 4 0 6 time s r yout 正弦信号跟踪 M 3 00 511 522 533 544 55 0 1 0 0 1 0 2 0 3 0 4 0 5 0 6 time s error 误差曲线 M 3 0 100 10 20 30 40 50 6 3 2 1 0 1 2 3 4 e de 相轨迹 M 3 00 511 522 533 544 55 1 0 1 2 3 4 5 6 time s s 切换函数 s M 3 00 511 522 533 544 55 1 0 5 0 0 5 1 1 5 2 time s u 控制器输出 M 3 仿真程序如下仿真程序如下 1 主程序 figure2 8 m clear all close all global a b c A F M ep k delta ts 0 001 T 5 TimeSet 0 ts T c 5 0 para t x ode45 figure2 8eq TimeSet 0 5 0 para x1 x 1 x2 x 2 r A sin 2 pi F t dr A 2 pi F cos 2 pi F t ddr A 2 pi F 2 sin 2 pi F t s c r x 1 dr x 2 if M 1 slaw ep sign s k s Exponential velocity trending law u 1 b c dr x 2 ddr slaw a x 2 elseif M 2 kk 1 delta for i 1 1 T ts 1 if s i delta sats i 1 elseif abs s i delta sats i kk s i elseif s i delta sats 1 elseif abs s delta sats kk s elseif s delta sats 1 end slaw ep sats k s else