北航自控原理课件7(英文版)

非线性系统分析AnalysisofNonlinearControlSystems北京航空航天大学自动控制原理第七章Contents目录非线性系统分析的核心理论框架与方法体系01非线性系统概述与典型特性02相平面法基本概念与绘制03相平面法在系统分析中的应用04描述函数法原理与应用05两种方法的综合对比与总结CHAPTER01非线性系统概述OverviewofNonlinearSystems—理解非线性现象的本质与工程意义CHAPTER07·NONLINEARSYSTEMS非线性系统的基本定义非线性系统是指包含一个或多个非线性元件或环节的控制系统，其核心特征是不满足叠加原理。由于几乎所有实际物理系统都含有某种程度的非线性，因此掌握非线性分析方法对于工程实践具有不可替代的重要意义。工程中的非线性控制系统实例DEFINITION系统中至少包含一个非线性环节，输入输出关系不能用线性微分方程完整描述，数学模型通常含有变量的乘积、幂次或分段函数。KEYDIFFERENCE不满足叠加原理：两个输入同时作用的总响应不等于各自单独作用响应之和，拉普拉斯变换和传递函数方法不再普遍适用。WHYITMATTERS几乎所有实际物理系统都存在固有非线性——摩擦、间隙、饱和等，忽略这些因素可能导致分析结果严重偏离实际行为。NonlinearCharacteristics四种典型非线性特性工程中绝大多数非线性问题可以归纳为饱和、死区、间隙和继电器四种典型特性——理解它们是进行非线性系统分析的基石。饱和特性Saturation输入超过阈值后输出保持恒定，等效增益随输入增大而下降，大信号时响应变慢运放输出限幅、电机最大转矩、液压阀最大流量限制等效增益↓死区特性DeadZone输入在零附近某一范围内输出为零，超过死区后才产生线性响应静摩擦力、液压阀重叠量、传感器最小分辨率阈值，造成稳态误差稳态误差↑间隙特性Backlash主动件改变运动方向时从件存在延迟响应，输入输出呈滞环形状齿轮传动间隙、丝杠螺母副间隙，导致精度下降和自激振荡滞环延迟继电器特性Relay输出仅有有限个离散状态（±1），输入超阈值时切换，典型强非线性温控开关、电磁继电器、比较器电路，用于bang-bang最优控制Bang-BangCHAPTER07·NONLINEARDYNAMICS非线性系统的独特现象非线性系统展现出多种线性系统中不可能出现的独特行为，包括自激振荡（极限环）、跳跃谐振和谐波响应等。这些现象源于非线性方程的复杂动力学特性，是驱动我们发展专门分析方法的根本原因。01极限环系统在无外部周期激励下可自发产生稳定的等幅振荡，振荡频率和幅值由系统本身参数决定，这是非线性系统最典型的特征之一。02滞回激励频率缓慢连续变化时，稳态振幅和相位可能发生突变跳跃，正向扫频和反向扫频的跳跃点不同，形成滞回现象。03谐波畸变对正弦输入的稳态响应不再是同频率正弦，而是包含基波和多个高次谐波分量，使得传统频率响应的概念需要推广和修正。04分岔非线性系统可能存在多个平衡状态，系统行为随参数变化发生质的突变，这是线性系统完全不具备的复杂动力学特征。自激振荡的典型工程实例——飞机机翼颤振现象NONLINEARSTABILITY非线性系统的稳定性特点非线性系统的稳定性分析远比线性系统复杂：系统可能存在多个平衡点，每个平衡点具有不同的稳定性特征，且稳定性通常是局部的而非全局的。多平衡点特征非线性系统可能拥有多个平衡状态，如单摆有下垂稳定和倒立不稳定两个平衡点，每个平衡点需要独立进行稳定性分析。这与线性系统仅有一个平衡点的特性形成鲜明对比。局部稳定与全局稳定非线性系统的渐近稳定性通常只在平衡点附近的某一区域内成立，超出该区域系统可能趋于其他平衡点或发散。只有在特定条件下才能实现全局渐近稳定。李雅普诺夫稳定性定义平衡点x*称为渐近稳定的，若从其附近出发的轨迹随时间趋于无穷时收敛于x*。这与BIBO稳定是从不同角度定义的概念，二者在非线性系统中并不等价。稳定性的参数依赖性稳定性不仅取决于系统结构，还受初始条件和参数值的影响。参数变化可能导致稳定性质的突变，即分岔现象，这使得非线性系统的预测和控制更具挑战性。CHAPTER02相平面法基本概念与绘制PhasePlaneMethod—FundamentalsandConstructionTechniquesNONLINEARSYSTEMS·基本概念相平面法的基本概念相平面法是一种适用于二阶系统的几何分析方法，通过在状态变量构成的平面上绘制运动轨迹来直观展示系统的动态行为。它将时间响应信息压缩到二维图形中，无需解析求解微分方程即可判断系统的稳定性与运动趋势。01相平面PhasePlane以二阶系统的两个状态变量x₁和x₂为坐标轴构成的二维平面，系统中每个状态对应平面上的一个点。状态空间02相轨迹PhaseTrajectory系统状态随时间演化在相平面上描绘出的曲线，消除时间变量t，直接展示x₂与x₁之间的函数关系。消去时间t03相图PhasePortrait从不同初始条件出发的多条相轨迹的集合，完整描绘系统在所有可能初始条件下的全局动态行为。全局动态04方法适用范围严格适用于二阶自治系统（不显含时间t），也可通过分段线性化近似处理某些特定的非自治系统。二阶自治Chapter7·非线性系统分析奇点分类与特征值判据奇点（平衡点）是非线性系统相平面分析的核心，通过在奇点处线性化系统并分析雅可比矩阵的特征值，可将奇点分为六种典型类型。FOCUS&NODE焦点与节点特征值为复数或同号实数相轨迹螺旋收敛/发散或直接趋向/偏离SADDLE&CENTER鞍点与中心特征值异号实数或纯虚数流形分离或封闭周期轨道01·稳定焦点特征值为一对负实部复数，相轨迹呈螺旋线收敛于奇点，系统表现为衰减振荡02·不稳定焦点特征值为一对正实部复数，相轨迹呈螺旋线远离奇点，系统表现为增幅振荡03·稳定节点特征值为两个负实数，相轨迹直接收敛于奇点，系统无振荡地趋于平衡04·不稳定节点特征值为两个正实数，相轨迹直接远离奇点，系统无振荡地偏离平衡05·鞍点特征值一正一负两个实数，存在稳定流形与不稳定流形，整体不稳定06·中心点特征值为一对纯虚数，相轨迹为封闭椭圆曲线，系统做等幅周期运动，临界稳定非线性系统分析·相平面方法相轨迹的绘制方法相轨迹的绘制主要有解析法和等倾线法两种途径。解析法通过积分获得精确轨迹方程但适用范围有限；等倾线法通过构建等斜率曲线族并连接短线段来近似绘制轨迹，是处理非线性系统最实用的定性分析工具。解析法AnalyticalMethod从系统状态方程出发，消去时间变量t，得到dx₂/dx₁=f(x₁,x₂)/g(x₁,x₂)的微分方程对该方程积分求解，得到相轨迹的解析表达式x₂=F(x₁,C)，其中C为由初始条件确定的积分常数适用范围有限，仅当微分方程可积时使用，大多数非线性系统无法获得解析解等倾线法IsoclineMethod令相轨迹斜率dx₂/dx₁=α（常数），解出α对应等倾线方程，在相平面上绘制一族等倾线在每条等倾线上标出斜率为α的短线段（方向场），从给定初始条件出发依次连接这些线段得到近似轨迹该方法虽为近似方法但适用范围广，可处理任意形式的非线性系统，是工程分析中最常用的绘图手段PhaseTrajectory·TimeDomain从相轨迹获取时间信息相轨迹本身不包含显式时间信息，但可以通过数值积分等方法从相轨迹中恢复出时间响应。利用状态变量之间的微分关系，将轨迹分段积分可获得各点对应的时刻，进而计算振荡周期、调节时间等关键时域指标。01数值积分法利用dx₁/dt=x₂的关系，沿相轨迹进行数值积分：t=∫(dx₁/x₂)，将轨迹分成若干小段逐段累加得到对应时刻。该方法适用于任意形状的相轨迹，精度取决于积分步长的选取。02环路积分法对于接近极限环的封闭轨迹，可通过计算一个完整周期的环路积分T=∮(dx₁/x₂)来获得自激振荡的周期。此方法避免了逐点计算的繁琐，直接得到系统的固有振荡特性。03增量法近似在轨迹上取相邻两点，用Δt≈Δx₁/ẋ₁或Δt≈Δx₂/ẋ₂估算时间间隔，取平均值可提高精度。这是一种简便的近似方法，适用于快速估算或实时分析场景。04时域指标提取时间信息的恢复使得我们可以从相图中读取调节时间、峰值时间、振荡周期等工程所需的时域性能指标。这为相平面分析法与经典控制理论的衔接提供了桥梁。PhaseTrajectoryAnalysis线性系统的相轨迹特征二阶线性系统ẍ+2ζωₙẋ+ωₙ²x=0的相轨迹形状完全由阻尼比ζ决定，呈现稳定节点、稳定焦点、中心点和不稳定焦点四种典型形态。理解这些基本模式是分析复杂非线性系统相平面行为的必要基础。过阻尼特征值为两个负实数，相轨迹从任意初始点单调收敛于原点，对应稳定节点，无振荡过程ζ>1欠阻尼特征值为一对负实部共轭复数，相轨迹呈螺旋线收敛于原点，对应稳定焦点，振荡逐渐衰减0<ζ<1无阻尼特征值为一对纯虚数，相轨迹为围绕原点的封闭椭圆族，对应中心点，系统等幅振荡永不停止ζ=0负阻尼特征值实部为正，相轨迹从原点附近螺旋向外发散，对应不稳定焦点，振幅不断增大直至系统失稳ζ<0Chapter03相平面法的系统分析应用PhasePlaneAnalysisApplications—FromTheorytoEngineeringPracticePhasePlane·PiecewiseLinear分段线性系统的相平面分析分段线性化是处理实际非线性系统最常用的相平面分析策略。通过将非线性特性划分为多个线性区域，在每个区域内独立绘制相轨迹并在边界处进行拼接，可以将复杂的非线性问题转化为多个简单线性问题的组合求解。01分区原则根据非线性特性的分段点（如饱和阈值、死区边界、继电器切换面），将相平面划分为若干个线性工作区域。每个区域对应系统的一种线性化工作状态，便于独立分析各区域内的动力学行为。02区域独立分析在每个区域内，系统退化为线性系统，可求出该区域的奇点位置和类型，绘制对应的相轨迹族。利用线性系统的成熟理论，快速确定局部稳定性和动态响应特性。03边界拼接规则相轨迹到达区域分界线时，该点即为相邻区域轨迹的起始点，通过逐段拼接得到完整的全局相轨迹。保证轨迹在边界处的连续性，形成跨区域的系统整体运动图景。04切换面上的特殊行为轨迹可能在切换面上发生滑动运动（SlidingMode），沿切换面滑向某一方向，这是分段线性系统特有的动力学现象。滑动模态的分析需要考察切换面两侧向量场的相对关系。PhasePlaneAnalysis饱和非线性系统的相平面分析带饱和特性的二阶系统是分段线性相平面分析的经典案例。饱和使大信号区域等效增益下降、系统响应变慢，但通常不会改变系统的整体稳定性。线性工作区LinearRegion当系统误差e的绝对值小于饱和阈值时，系统按线性模型运行，奇点位于原点若线性区设计为欠阻尼（0<ζ<1），相轨迹呈螺旋线收敛于原点，系统表现为衰减振荡线性区的动态性能指标（超调量、调节时间等）可用标准二阶系统公式直接计算饱和工作区SaturationRegion当|e|超过饱和阈值，控制量被限幅，系统等效增益下降，表现为阻尼增大、响应变慢饱和区相轨迹弯曲程度减小，收敛速度降低，从饱和区进入线性区的切换点决定了后续振荡幅度饱和特性的主要影响是增大了大偏差情况下的调节时间，但一般不会导致系统失稳PhasePlaneAnalysis继电器系统的极限环分析继电器控制系统因其输出的离散切换特性，在相平面分析中经常产生极限环（自激振荡）。通过在切换面两侧分别绘制相轨迹并分析拼接结果，可以判定极限环的存在性、振荡幅值和频率，以及极限环的稳定性。01相平面分区：继电器系统将相平面分为e>0和e<0两个区域，每个区域内系统按不同线性方程运动，切换面为e=0直线。02封闭环路形成：当两侧区域的相轨迹在切换面上首尾相接形成封闭环路时，系统产生极限环，对应持续的等幅自激振荡。03稳定极限环：附近轨迹趋向于该封闭环，表明系统在扰动后能恢复到等幅振荡状态，是工程中最常遇到的情形。04不稳定极限环：附近轨迹远离封闭环，表明等幅振荡状态不能维持，系统要么收敛到平衡点要么发散到更大振幅。05振荡参数计算：极限环的振荡周期T可通过沿封闭轨迹积分计算：T=∮dt，振荡幅值由封闭环的最大坐标值确定。PhasePlaneAnalysis死区非线性系统的相平面分析死区特性使控制器在小误差范围内失去作用，导致系统进入死区后无法继续收敛到原点，从而产生不可消除的稳态误差。相平面分析清晰地揭示了死区宽度与稳态误差之间的定量关系，为工程设计中权衡死区补偿策略提供了理论依据。01三区划分结构死区将相平面分为三个区域：|e|<Δ的死区（无控制作用）和|e|>Δ的两个线性控制区，死区内系统仅靠惯性自由运动。|e|<Δ02死区内轨迹特征在死区区域内，由于控制力为零，相轨迹为水平直线（ẋ恒定）或抛物线，系统状态不再向原点收敛。惯性运动03稳态误差定量关系轨迹进入死区后停止运动的位置即为系统最终稳态值，该稳态误差的大小由死区宽度Δ和进入死区时的状态决定。Δ→稳态值04工程补偿策略减小死区对稳态精度影响的工程策略包括：提高执行机构灵敏度、引入积分补偿或在死区边界处增加高增益区间。积分补偿NONLINEARSYSTEMANALYSIS间隙（回差）非线性系统分析间隙非线性因传动机构中主动件与从动件之间的空隙而产生，输入输出关系呈现滞环特征。它不仅降低系统定位精度，更可能在闭环系统中引入额外相位延迟并导致自激振荡。相平面分析可揭示间隙宽度与振荡幅值的定量关系。精密齿轮传动与消隙机构工程场景01滞环特性：主动件反向运动时，从件在间隙宽度b内无响应，导致信号传递存在延迟和失真02相平面表现：系统轨迹在切换面附近反复穿越，收敛速度明显减慢，严重时收敛为封闭极限环03自激振荡：振荡频率和幅值与间隙宽度b、系统增益和固有频率密切相关，b越大振荡越剧烈04工程措施：采用高精度消隙齿轮、双片齿轮预紧、软件补偿算法以及提高传动链刚度Chapter04描述函数法原理与应用DescribingFunctionMethod—PrinciplesandEngineeringApplicationsNonlinearSystemAnalysis描述函数的定义与基本思想描述函数法的核心思想是用非线性环节输出基波分量与正弦输入的复数比来定义等效增益N(A)，从而将非线性系统近似为幅值依赖的线性系统进行分析。该方法成立的前提是线性部分具有良好的低通滤波特性，能有效衰减高次谐波分量。01描述函数定义当非线性环节输入正弦信号x(t)=Asin(ωt)时，N(A)=Y₁/A·ejφ₁，其中Y₁和φ₁分别为输出基波的幅值和相位。02物理意义N(A)相当于非线性环节的"等效频率响应增益"，但它仅依赖于输入幅值A而非频率ω，这是与线性传递函数的关键区别。03基本假设条件系统的线性部分G(jω)具有低通滤波特性，能有效衰减非线性环节输出的高次谐波分量，使闭环信号近似为正弦。04计算步骤将输出y(t)展开为傅里叶级数，提取基波分量的系数A₁和B₁，则N(A)=(B₁+jA₁)/A。NONLINEARANALYSIS·DESCRIBINGFUNCTIONS典型非线性的描述函数计算不同非线性特性的描述函数N(A)具有不同的数学表达和变化规律。掌握这些典型结果是运用描述函数法分析系统的基础。理想继电器01输出为±M的方波信号，傅里叶展开后基波幅值B₁=4M/π，描述函数N(A)=4M/(πA)为纯实数02N(A)与输入幅值A成反比，A越大等效增益越小；无相位偏移，因为继电器是单值奇函数非线性03常用于开关控制与bang-bang调节，是研究非线性系统极限环特性的典型模型IDEALRELAY饱和特性01当A≤a（阈值）时系统工作在线性区，N(A)=k；当A>a时N(A)<k且随A增大单调递减02描述函数表达式含arcsin项，极限情况A→∞时N(A)→0，反映饱和对大信号的强衰减作用03广泛存在于放大器、执行机构中，是限制系统输出幅值、防止过载的核心非线性环节SATURATION死区特性01当A≤Δ（死区宽度）时无输出，N(A)=0；当A>Δ时N(A)随A增大单调递增并趋向于线性增益k02死区与饱和的描述函数互补：二者之和恒等于线性增益k，可用于简化复合非线性分析03常见于阀门、电机启动与传感器阈值，导致系统对小信号不响应、存在调节静差DEADZONESTABILITYANALYSIS描述函数法的稳定性分析描述函数法将非线性系统的稳定性分析转化为奈奎斯特判据的推广：通过判断G(jω)曲线与-1/N(A)曲线在复平面上是否存在交点来预测自激振荡。交点处的A值和ω值分别对应振荡的幅值和频率，而交点处两条曲线的相对走向决定了极限环的稳定性。01自激振荡条件G(jω)·N(A)=-1，等价于G(jω)曲线与-1/N(A)曲线在复平面上相交，交点即为可能的振荡状态。G·N=-102交点参数读取交点处对应的频率ω即为自激振荡频率，对应的幅值A即为自激振荡幅值，两者可直接从图中获得。ω与A03稳定性判定规则若-1/N(A)曲线随A增大的方向从G(jω)曲线的不稳定侧穿向稳定侧，则该交点对应稳定极限环。极限环稳定04无交点的情况若G(jω)曲线不包围-1/N(A)曲线的任何部分，则系统在任何初始条件下均不产生自激振荡，系统全局稳定。全局稳定STABILITYANALYSIS·四步法描述函数法分析实例通过完整的四步分析流程——计算描述函数、绘制-1/N(A)曲线、画出G(jω)奈奎斯特图、判断交点与稳定性——可以系统地预测非线性系统是否产生自激振荡及其参数。01计算N(A)根据非线性特性类型求出描述函数表达式。例如理想继电器N(A)=4M/(πA)，则−1/N(A)=−πA/(4M)，为负实轴上的射线。N(A)02绘制曲线在复平面上同时画出G(jω)的奈奎斯特图（ω:0→∞）和−1/N(A)曲线（A:0→∞），标注方向箭头以明确参数变化趋势。G(jω)03寻找交点观察两条曲线是否存在交点。若存在则记录交点处的ω值和A值，分别对应自激振荡的频率和幅值。ω,A04判定稳定性在交点处，若−1/N(A)沿A增大方向从G(jω)包围区域内部穿出到外部，则该极限环稳定；反之则不稳定。稳定判据Limitations描述函数法的局限性与适用条件描述函数法是一种近似分析工具，其精度依赖于线性部分的低通滤波特性和系统结构的规范性。理解其适用条件和局限性对于避免误用、正确评估分析结果至关重要。在实际工程中应与仿真验证相结合以确保结论的可靠性。01近似性局限描述函数法忽略高次谐波的影响，当线性部分G(jω)的高频衰减不够充分时，分析结果与实际系统行为可能产生显著偏差。02结构限制标准描述函数法要求系统可化为一个非线性环节与一个线性环节串联的典型结构，多个非线性环节需先等效合并。03信息有限该方法仅能预测自激振荡的存在性、幅值和频率，无法提供系统的完整瞬态响应和收敛过程的详细信息。04工程建议分析结果应作为初步判断，最终结论需结合数值仿真或实验验证，特别是对于高阶和复杂非线性系统。Chapter05方法综合对比与课程总结ComparativeSummary—PhasePlanevs.DescribingFunctionMethodNONLINEARANALYSIS相平面法与描述函数法对比相平面法和描述函数法是非线性系统分析的两大核心工具，各有所长：相平面法精确但限于二阶系统，描述函数法可处理高阶系统但为近似方法。对比维度相平面法描述函数法适用阶次仅适用于二阶系统可处理任意阶次的系统分析精度精确方法（无近似假设）近似方法（忽略高次谐波）分析内容给出完整的动态行为全景图主要预测自激振荡的存在与参数数学基础基于微分方程的几何理论基于谐波平衡的频域近似图形工具相平面上的轨迹与奇点分析G(jω)与−1/N(A)的奈奎斯特图工程适用性适合低阶系统的定性分析适合高阶系统的稳定性预测信息丰富度提供全局动态行为信息仅提供稳态振荡信息两种方法各有优劣，工程中应根据系统阶次和分析目标选择合适的方法Chapter07·Summary本章核心知识点总结本章系统学习了非线性系统分析的理论基础与两大核心方法。从非线性特性的识别到相平面的精确分析，再到描述函数的频域近似，形成了完整的非线性系统分析知识体系，为后续现代控制理论的学习奠定了坚实基础。01非线性基本概念掌握本质非线性与人为非线性的分类，理解饱和、死区、间隙、继电器四种典型特性的输入输出关系及其对系统性能的影响。能够识别实际系统中的非线性环节，分析其对稳定性、响应速度和稳态精度的作用机理。02相平面法核心技能能判断奇点类型（焦点、节点、鞍点、中心），掌握解析法和等倾线法绘制相轨迹，会对分段线性系统进行分区分析和边界拼接。能够根据相图判断系统的稳定性、瞬态响应特性及是否存在极限环振荡。03描述函数法核心技能能计算典型非线性的描述函数N(A)，会在复平面上绘制G(jω)和−1/N(A)曲线，能判断交点处极限环的稳定性和振荡参数。掌握描述函数法的频域分析思路，理解谐