系统非线性环节的仿真

本文格式为Word版，下载可任意编辑——系统非线性环节的仿真2.5系统非线性环节的仿真

在实际系统中，往往存在各种非线性特性，可将此当作非线性环节处理，这种环节的输入和输出之间关系是一种非线性函数关系，因此非线性环节的仿真就是用仿真语言来描述这些关系。本节介绍几种典型的非线性环节的仿真算法。2.5.1饱和环节

饱和环节在控制系统中较普遍，例如饱和放大器、限幅装置、伺服阀饱和特性等。饱和环节特性如下图。

图2.5-1饱和特性

该特性对应的数学表达式为：

?uu?c?y??cu?c??cu??c?(2.5-1)

式中，c为饱和环节特征参数，斜率为1，该环节特性可用MATLAB编程仿真，利用上面算法的编写的MATLAB函数SATURATION，调用格式为：

y?saturation(u,c)

其中，u为输入；c为饱和环节特征参数，y为饱和环节输出。Saturation.m;amp209.m

2.5.2死区环节

在控制装置中，放大器的不灵敏区，伺服阀和比例阀阀芯正遮羞特性，传动元件静摩擦等造成的死区特性。典型死区非线性环节特性如图2.5-2所示。

图2.5-2死区特性

可用下面数学关系来描述：

?0u?c?y??u?cu?c?u?cu??c?(2.5-2)

式中，c为死区特征参数，斜率为1。

该环节可根据上述算法编写MATLAB函数deadzone供调用，格式如下：

y?deadzone(u,c)

其中，u为环节输入；c为死区环节特征参数，y为死区环节输出。

Deadzone.m;amp210.m。

2.5.3齿隙非线性环节

齿轮传动副和丝杆螺母传动副中存在传动间隙都属这一类非线性因素，它对系统精度带来影响。齿隙非线性环节特性如图2.5-3所示。

图2.5-3齿隙特性

当输入u增加时，输出沿a?b?d线段变化；当输入u减小时，输出沿d?e?a线段变化。在线段bd上，输入增加时，当前输出值y(k)总是大于前一时刻的输出值y(k-1)。而在ea上，输入减小时，当前输出y(k)总是小于前一时刻的输出值y(k-1)。在ab段和de段，y(k)=y(k-1)。以上特性的数学描述如下：

?u(k)?c当u(k)?u(k?1)?0且y(k?1)?u(k)?c?y(k)??u(k)?c当u(k)?u(k?1)?0且y(k?1)?u(k)?c

?y(k?1)其余?（2.5-3）

式中，c为齿隙环节特征参数，斜率为1。

根据（2.5-3）算法编写的MATLAB函数backlash，调用格式如下：

y1?backlash(u1,u0,y0,c)

其中，u0,u1分别为前一时刻和当前时刻输入值；y0,y1分别为前一时刻和当前时刻输出值；c为齿隙特征参数。Backlash.m,amp211.m.

2.5.4继电非线性环节

继电非线性环节特性如图2.5-4所示，(b)为具有死区继电环节，(a)为(b)特例，即c=0。对于图(b)所示特性，可用下面数学关系描述：

图2.5-4继电环节特性

?0u(k)?c?y??1u(k)?c??1u(k)??c?(2.5-4)

根据（2.5-4）算法编写的MATLAB函数relaydead，调用格式如下：

y?relaydead(u,c)

其中，u为环节输入，c为环节输出，y为死区特征参数。Relaydead.m,amp212.m.

2.5.5具有滞环的继电环节

具有滞环的继电非线性环节特性如下图，可用下面数学关系描述：

图2.5-5带滞环继电特性

u(k)?u(k?1)且u(k)?c?h?y(k)???hu(k)?u(k?1)且u(k)??c

?y(k?1)其他?（2.5-5）

根据（2.5-5）算法编写的MATLAB函数relaydelay，调用格式如下：

y?relaydelay(u1,u0,y0,c,h)

其中，u0,u1分别为前一时刻和当前时刻环节输入值；y0,y1分别为前一时刻和当前时刻输出值；c,h为环节特征参数。Relaydelay.m,amp213.m.

2.5.6库仑—粘性摩擦力环节

图2.5-6库仑-粘性摩擦特性

机械摩擦具有非线性如下图，它由库仑摩擦和粘性摩擦两部分组成。当物体运动速度不为零时，出现的库仑摩擦力为和相对运动速度相反的、与速度大小无关的恒定力，而粘性摩擦力则速度大小成正比。上述特性可用数学关系来描述：

?cv(k)?0且v(k?1)?0?f(k)???cv(k)?0)且v(k?1)?0

?sign(u)?(G?u?c)其他?(2.5-6)

根据式（2.5-6）的算法，编写的MATLAB函数friction，调用格式为：

y?friction(u,u1,c,G)

其中，u,u1分别为当前时刻和下一时刻的输入值（速度）；y为输出值（摩擦力）；c为库仑摩擦力值；G为粘性摩擦系数。

2.6采样控制系统仿真

2.6.1采样控制系统的基本组成

采样控制系统是指系统一处或几处信号是以采样形式来传递的，而被控对象是连续的。典型的采样控制，是一种连续—离散

混合系统，系统如图2.6-1所示。

图2.6-1采样控制系统

随着计算机技术的发展，计算机控制系统成为采样控制系统的主要形式。系统（a）和系统（b）的主要区别在于系统闭环的输入信号类型，前者是模拟信号，后者是数字信号。系统中的A/D转换器是采样开关，将连续模拟量转变为离散的数字量并根据计算量的位数进行进行量化处理；D/A转换器是将离散的数字量转变为模拟量，同时它也是一个保持器，且一般具有零阶保持器特征。A/D和D/A转换器在模拟量和数字量之间存在确定的比例关系，相当于系统中的比例环节。2.6.2采样控制系统仿真特点

采样控制系统包含连续部分和离散部分。对于连续部分仿真方法和连续系统仿真一样，可采用数值积分法或离散相像法。若采用数值积分法则需要确定积分步长，若采用离散相像法则需要确定虚拟的采样周期。对于离散部分，A/D转换器和D/A转换

器是实际存在的，采样周期和保持器类型也均是实际存在的。因此采样控制系统仿真中，仿真步距或虚拟采样周期和系统实际采样周期之间存在同步问题。2.6.3仿真步长和采样周期

对采样控制系统进行仿真时，连续部分仿真步长的选择必需根据被控对象的动态特性、系统采样周期大小、仿真精度的要求等综合考虑。一般按下面两种状况处理：

（1）仿真步长h等于采样周期T；（2）仿真步长h小于采样周期T。

第一种方法适用于系统连续部分参数变化较缓慢或系统幅值穿越频率?c较小的系统。对于大多娄机电类采样控制系统，系统连续部分参数变化较快，常采用其次种方法，以保证系统中连续部分足够的仿真精度。

若仿真步长h小于采样周期T，为了全球仿真程序的实现，寻常取采样周期T恰好是仿真步长h的整数倍速，即h?TN，其

中N为正整数。采样系统仿真采用定步距，对于连续部分在每个计算点均作仿真运算，而对于离散部分（数字控制器）只有在采样时刻才执行仿真运算，在其他计算点不执行仿真运算。在仿真程序中，采样时刻点可借助仿真时间/采样周期为整数的关系来确定。采样控制系统仿真程序流程如图2.6-2所示。

2.6.4采样控制系统仿真方法

如前所述，采样控制系统分为连续部分和离散部分。对于连续部分仿真有不同的处理方法，采样控制系统仿真方法分为两种。

一、基于数值积分法

图2.6-2采样控制系统仿真程序流程图

对系统连续部分仿真采用数值积分法，图2.6-2是基于数值积

分的采样控制系统仿真程序框图。这种方法需选择连续部分仿真步长、仿真数值积分方法等。一般采用定步距，且仿真步长一般小于离散部分采样周期。离散部分仿真是基于递推法，十分简单。

二、基于离散相像法

系统的连续部分先进行z变换，若连续部分模型G(s)已知，则可借助MATLAB函数C2D将连续模型转换为离散模型G(z)，将G(z)和原系统离散部分模型D(z)合并后可求得采样控制系统的离散模型W(z)，这样就可以进行仿真运算。在连续部分开散化时，可优先选择虚拟的采样周期和系统实际采样周期一致。在仿真运算时可根据仿真速度和精度要求，采用