采样系统鲁棒稳定性分析的新方法
鲁棒性优化的原理、评估方法及应用-放射医学论文-基础医学论文-医学

鲁棒性优化的原理、评估方法及应用-放射医学论文-基础医学论文-医学鲁棒性优化的原理、评估方法及应用放射医学论文基础医学论文医学放射医学作为一门重要的医学分支,应用广泛且发展迅猛。
在放射医学的实践中,为了保证诊断结果的准确性和稳定性,提高影像质量和疾病诊断的可信度,鲁棒性优化成为一种重要的手段。
本论文将着重探讨鲁棒性优化的原理、评估方法以及其在放射医学中的应用。
一、鲁棒性优化原理鲁棒性优化是指在实际应用中,通过在系统中引入一定程度的冗余,使得系统对各种干扰因素和不确定性具有强健性。
在放射医学领域中,鲁棒性优化的原理主要包括以下几个方面。
1. 信号处理技术鲁棒性优化中的信号处理技术主要针对图像数据的处理。
比如在辐射剂量计算中,为了减小各种因素对剂量计算结果的影响,可以基于模型订正或者增加剂量分配的冗余,提高系统的鲁棒性。
2. 特征提取与选择特征提取与选择是鲁棒性优化的关键环节。
通过合理选择影像中的关键特征,可以减少噪声和其他干扰因素对诊断结果的影响。
比如在肿瘤检测中,可以通过计算形状特征、纹理特征等来提高肿瘤检测的准确性和鲁棒性。
3. 算法优化算法优化是鲁棒性优化的重要手段。
通过改进或设计新的算法,可以提高系统对各种噪声和变化的适应能力。
例如,对于放射源和探测器位置的微小变化,可以采用基于机器学习的方法来优化图像重建算法,从而提高图像质量和诊断准确性。
二、鲁棒性优化的评估方法为了评估鲁棒性优化的效果,我们需要选择合适的评估方法和指标。
以下是几种常用的评估方法。
1. 灵敏度分析灵敏度分析是评估系统对输入参数变化的鲁棒性的一种方法。
通过改变系统参数或输入数据的扰动幅度,观察输出结果的变化情况,可以评估系统在不同干扰因素下的鲁棒性。
2. 参数估计参数估计是通过对输入参数进行统计分析,估计系统对参数变化的鲁棒性。
通过观察参数估计结果的方差、置信区间等指标,可以评估系统在不同干扰条件下对参数的稳定性和可信度。
工业过程中系统鲁棒性研究及仿真分析

工业过程中系统鲁棒性研究及仿真分析随着工业化进程的推进,各种工业生产系统也得到了迅猛发展。
在这个过程中,除了注重生产效率、成本控制等指标之外,还需要考虑生产过程的稳定性和可靠性,以保证产品的质量和安全。
系统鲁棒性研究及仿真分析因此成为了工业领域的一个重要课题。
一、什么是系统鲁棒性所谓系统鲁棒性,通俗地说,就是指系统在面对外部环境变化、内部异常等因素时的稳定性和适应能力。
不同的系统在不同的环境下,面对各种不同的干扰或者异常,其表现和反应也会有所不同。
在工业生产过程中,系统鲁棒性的研究对于确保生产过程的稳定性和性能优化至关重要。
例如,在计算机制造过程中,如果生产系统的鲁棒性不够强,那么在组装、测试等过程中出现的异常情况可能会导致计算机产品存在缺陷或者无法正常工作。
因此,对于提高工业生产系统生产效率、可靠性和维修成本控制等方面具有至关重要的意义。
二、系统鲁棒性的研究方法如何研究系统的鲁棒性呢?从数学上讲,系统鲁棒性可以用一些参数来描述。
例如,对于一个控制系统或者人工智能算法,可以描述其对于外界参数、初始条件和噪声的容忍程度。
为了测量系统的鲁棒性,可以通过仿真分析等方式进行研究。
仿真分析可以模拟出工业生产过程中各种异常情况,从而测试系统的适应能力和弹性。
经过仿真分析研究后,可以得到一个系统鲁棒性指标。
这个指标通常由以下几个方面构成:1. 参数不确定性:即系统中参数变化或者干扰的容忍度。
2. 初始条件不确定性:即系统对于初始输入或者状态的容忍度。
3. 物理噪声:指物理环境中可能会对系统造成影响的因素。
4. 模型不确定性:指模型误差可能引起的不确定性。
通过对这些指标进行分析,可以评估系统的稳定性和可靠性,从而确定如何对系统进行改进。
三、系统鲁棒性的应用系统鲁棒性研究的应用领域非常广泛,除了工业生产系统之外,还包括金融、航空航天、医疗等领域。
具体应用如下:1. 工业生产领域在制造业、物流、供应链和质量管理等方面,系统鲁棒性的研究和应用都有很多的实例。
单级倒立摆LQR控制方法的鲁棒稳定性分析

单级倒立摆LQR控制方法的鲁棒稳定性分析刘微微;张静【摘要】鲁棒稳定性是控制系统的一个重要指标,针对单级倒立摆系统,采用LQR 最优控制方法,通过MATLAB实验环境下的仿真,得到了良好的控制效果.重点对该控制方法的鲁棒稳定性做了详细分析,通过增加系统自身的扰动及LQR控制器中加权阵R的改变考察该控制方法的鲁棒稳定性,对比仿真结果表明该控制方法鲁棒稳定性良好,进一步验证了LQR控制方法在对于倒立摆这样一个非线性极强的控制系统上所具有的控制优势.【期刊名称】《黑龙江大学工程学报》【年(卷),期】2010(037)002【总页数】4页(P105-108)【关键词】鲁棒稳定性;倒立摆;LQR;最优控制【作者】刘微微;张静【作者单位】哈尔滨理工大学,自动化学院,哈尔滨,150080;哈尔滨理工大学,自动化学院,哈尔滨,150080【正文语种】中文【中图分类】TP130 引言对于单级倒立摆系统,目前已有多种控制方法可对其实现稳摆控制。
典型的有线性PID控制[1-2]、常规PID控制[3]、LQR控制、智能控制[4]等。
其中,LQR控制方法的优势在于其控制方案简单,超调量小,且响应速度快,该方法不仅对单级倒立摆系统能够进行有效控制,且已经成功地应用于直线双倒立摆[5]和双足机器人[6]的控制。
鲁棒稳定性是衡量控制系统控制品质优劣的一个重要指标[7],鉴于LQR控制方法已有的控制优势,有必要对该控制方法的鲁棒稳定性做详细分析,以进一步验证LQR控制方法的控制优势。
本文针对单级倒立摆系统,完成了系统建模及LQR控制的MATLAB仿真,通过增加系统自身的扰动及LQR控制器中加权阵R的改变考察该控制方法的鲁棒稳定性,对比仿真结果表明该控制方法鲁棒稳定性良好。
在较大的系统参数摄动下系统达到稳态的时间变化较小。
1 单级倒立摆系统建模单级倒立摆的控制目标是:由倒摆和小车组成的倒立摆在适当的控制力作用下,在有限长度的导轨上,受到干扰后,倒立摆仍然能够竖直立稳[8],即θ≈0,且小车位移x≈0(图1)。
机床网络分布式伺服系统的稳定性与鲁棒性分析与优化策略

机床网络分布式伺服系统的稳定性与鲁棒性分析与优化策略摘要:机床是现代制造业中关键的工具,而网络分布式伺服系统的稳定性与鲁棒性对机床性能的影响至关重要。
本文旨在分析机床网络分布式伺服系统的稳定性与鲁棒性,并提出相应的优化策略。
首先介绍了机床网络分布式伺服系统的基本原理和结构,并分析了系统中可能存在的稳定性和鲁棒性问题。
接着,通过数学建模和仿真实验,对机床网络分布式伺服系统的稳定性和鲁棒性进行了详细分析。
最后,提出了一些优化策略,如控制算法优化、故障诊断与容错措施等,以提高机床网络分布式伺服系统的稳定性和鲁棒性。
1. 引言机床网络分布式伺服系统是由多个伺服节点组成,并通过网络进行通信和控制的一种先进控制系统。
该系统架构具有分布式、协作和实时性强的特点。
然而,由于网络通信延迟、信号干扰、节点故障等因素的存在,机床网络分布式伺服系统的稳定性和鲁棒性面临一些挑战。
因此,对该系统的稳定性和鲁棒性进行分析与优化具有重要意义。
2. 机床网络分布式伺服系统的稳定性分析机床网络分布式伺服系统稳定性是指系统在一定工作条件下是否能保持稳定的运行状态。
在稳定性分析中,需要考虑网络通信延迟对系统响应的影响、节点之间的通信协议、控制算法的选择等因素。
通过数学建模和仿真实验,可以得到系统的传输函数和稳定域,进而评估系统的稳定性。
另外,还需考虑系统中可能出现的饱和效应、奇异性问题等,并采取适当的措施进行补偿或优化。
3. 机床网络分布式伺服系统的鲁棒性分析机床网络分布式伺服系统的鲁棒性是指系统对于外部扰动和不确定性的抵抗能力。
在鲁棒性分析中,需要考虑节点之间通信的可靠性、节点故障的容错机制、外部干扰对系统性能的影响等因素。
通过故障诊断和容错措施,可以提高系统的鲁棒性。
此外,还需采用适当的控制算法和滤波器,以降低系统对于噪声和干扰的敏感性,提高系统的鲁棒性。
4. 机床网络分布式伺服系统的优化策略为了提高机床网络分布式伺服系统的稳定性和鲁棒性,可以采取一系列优化策略。
帆板控制系统的鲁棒性分析及改进方法研究

帆板控制系统的鲁棒性分析及改进方法研究简介:帆板控制系统是指用于调整帆板角度以准确捕捉风能的控制系统。
在帆板能源利用领域,鲁棒性是一个重要的指标,旨在保证系统在各种外部扰动下的稳定性和可靠性。
本文将对帆板控制系统的鲁棒性进行分析,并研究改进方法,以提升系统的稳定性和可靠性。
一、鲁棒性分析1. 外部扰动的分析:首先,对帆板系统中可能遇到的外部扰动进行详细分析。
这些扰动可能包括:风速变化、风向变化、船体运动等。
2. 系统响应的分析:通过数学模型建立系统的状态空间方程,并分析系统对于不同外部扰动的响应情况,考虑到系统的跟踪误差和稳定性。
二、鲁棒性改进方法研究1. 鲁棒控制设计:基于鲁棒控制理论,设计出一种对外部扰动具有强鲁棒性的控制器。
具体包括:a. H∞控制方法:利用H∞控制方法将系统的鲁棒性分析转化为一个优化问题,设计出具有强稳定性和鲁棒性能的控制器。
b. μ合成控制方法:利用μ合成控制方法对帆板系统进行频域分析,并设计出一个具有强鲁棒性的控制器。
2. 鲁棒估计器设计:针对帆板系统中存在的不确定性,设计出一种鲁棒估计器来对系统进行状态估计和鲁棒性优化。
具体包括:a. 鲁棒滤波器设计:采用鲁棒滤波器对传感器测量信号进行滤波和融合,以提高测量的准确性和可靠性。
b. 鲁棒辨识算法:利用鲁棒辨识算法对系统的参数进行估计和辨识,以提升系统的鲁棒性和准确性。
3. 鲁棒策略优化:通过优化策略,对帆板系统的鲁棒性进行进一步改进。
具体包括:a. 高鲁棒性控制策略:通过改进控制策略,增强系统对外部扰动的抵抗能力,提升鲁棒性和稳定性。
b. 多模型控制策略:利用多模型控制策略,将帆板系统分成不同的模型区域,并分别设计控制器,以提高系统的稳定性和鲁棒性。
总结:帆板控制系统的鲁棒性分析及改进方法的研究对于提升系统的稳定性和可靠性具有重要意义。
通过对外部扰动的分析,建立系统的数学模型,并设计合适的控制策略和估计器,可以提高系统对外部扰动的鲁棒性。
采样控制系统的稳定性分析

z 或 w 特征方程的系数,按照下述方法构造(2n-3)行、(n+1)列朱利阵列,见表8-2:
w 1 其中比较常用的代数判据就是劳斯判据。
式中z和w均为复数,分别把它们表示成实部和虚部相加的形式,即
其模等于1,与频率ω无关;其相角为ωT,随频率ω 而改变。
可见,S平面上的虚轴映射到Z平面上,为以原 点为圆心的单位圆。
当s位于S平面虚轴的左边时,σ为负数, z eT
小于1。反之,当s位于s平面虚轴的右半平面时,为 正数,z eT 大于1。s平面的左、右半平面在z平 面上的映像为单位圆的内、外部区域。
z 1 z 1
当动点z在Z平面的单位圆上和单位圆之内时,应满足:
式中z和w均为复数,分别把它们表示成实部和虚部 |a0|< an, |b0|>|bn-1|, |c0|>|cn-2|
可见,S平面上的虚轴映射到Z平面上,为以原点为圆心的单位圆。 082, 满足|b0|>|b3|
根据给定的D(z)知:
相加的形式,即 1(s),试求系统稳定时k的变化范围。
(1)当采样周期T分别为1(s),0.
z x jy
w u jv
因此,必须采用一种新的变换,使z平面上的单位圆,在新的坐标系中的映象为虚轴。
化简后,得W域特征方程
显然,闭环系统特征方程的根λ1、λ2、…λn即是闭环脉冲传递函数的极点。
因此,必须采用一种新的变换,使z平面上的单位圆,在新的坐标系中的映象为虚轴。
线性采样系统稳定的充要条件 图8-21:线性采样系统结构图
线性采样系统如图8-21所示。 其特征方程为
D(z) 1 GH(z) 0
机器学习中的模型稳定性与鲁棒性分析方法(八)

机器学习中的模型稳定性与鲁棒性分析方法在机器学习领域,模型的稳定性和鲁棒性分析是非常重要的,它们关乎模型的可靠性和泛化能力。
本文将探讨机器学习中模型稳定性和鲁棒性的分析方法,并探讨如何应对模型不稳定和不鲁棒的情况。
一、什么是模型稳定性?首先,我们来看看什么是模型稳定性。
在机器学习中,模型稳定性指的是当输入数据发生微小变化时,模型的输出结果是否会有显著变化。
如果模型对输入数据的微小变化非常敏感,那么这个模型就是不稳定的。
模型稳定性是评估模型可靠性的重要指标,稳定的模型更有利于在真实场景中得到准确的预测结果。
二、模型稳定性分析方法了解了模型稳定性的概念后,我们来看看如何进行模型稳定性分析。
一种常见的方法是使用bootstrap法。
Bootstrap法通过从原始数据集中随机抽取一定数量的样本,并使用这些样本来重新训练模型,然后观察模型输出结果的变化情况。
如果模型的输出结果在不同的bootstrap样本上变化较小,那么这个模型就可以认为是稳定的。
另一种常见的方法是使用交叉验证。
交叉验证通过将原始数据集划分为多个子集,然后使用不同的子集来训练模型,并使用剩余的子集来测试模型。
通过观察不同训练集上模型输出结果的一致性来评估模型的稳定性。
除了这些方法之外,还有一些统计学方法和图形方法可以用来分析模型的稳定性。
三、什么是模型鲁棒性?除了模型稳定性,模型鲁棒性也是非常重要的。
模型鲁棒性指的是当输入数据包含噪声或异常值时,模型是否能够产生合理的输出结果。
一个鲁棒的模型能够忽略输入数据中的噪声和异常值,从而得到更加稳健的预测结果。
四、模型鲁棒性分析方法针对模型鲁棒性,我们可以使用一些方法来进行分析。
一种常见的方法是使用敏感性分析。
敏感性分析通过在输入数据中引入噪声或异常值,来观察模型输出结果的变化情况。
如果模型对噪声或异常值的影响较小,那么可以认为这个模型是鲁棒的。
另一种方法是使用对抗性训练。
对抗性训练通过在训练过程中引入对抗样本,来提高模型对抗攻击的能力,从而提高模型的鲁棒性。
时滞不确定采样控制系统的鲁棒稳定性

时滞不确定采样控制系统的鲁棒稳定性刘彦文;王广雄;綦志刚;许保同【摘要】本文给出了一种可定量分析采样控制系统的时滞鲁棒稳定性的方法.因为采样系统的对象是连续时间的,所以对象中的时滞也应该是按连续时间来处理.文中指出,一个整数倍时滞是稳定的采样系统,可能会因为有并不很大的连续时间时滞而失稳.定义了一个新的变量w(t),用来描述这个不确定连续时间时滞带来的动特性.将w(t)的反馈回路分成与时滞无关和有关的两个部分,并提出了一种用频率响应来确定是否存在由不确定时滞引起的周期解的方法.用修正z-变换法和仿真验证了这个由图解解析所求得的解.本方法既可用于采样系统,也可用于一般的连续时间系统.%We propose a quantitative method for analyzing the robust stability of sampled-data systems with uncertain time-delays. Because the sampled-data systems are obtained from continuous-time systems by sampling, the time-delay in the sampled-data system must also be treated in the continuous-time system. It is pointed out that a stable sampled-data system with a time-delay equal to the integer-multiple of the sampling period may be destabilized by a small continuous time-delay. A new variable w(t) is defined to describe the dynamic response caused by the uncertain continuous time-delay. The feedback loop of w(t) is then divided into two parts. One depends on the uncertain time-delay, and the other is independent of the time-delay. A special frequency response method is proposed to determine the existence of the periodic solution of the system caused by the uncertain time-delay. The graphic-analytical solution is then verified by the modified z-transform method and by simulation. Theproposed method can also be used for robust stability analysis of continuous-time systems with time-delays.【期刊名称】《控制理论与应用》【年(卷),期】2013(030)002【总页数】5页(P238-242)【关键词】鲁棒稳定性;采样控制系统;连续时间时滞;频率响应;时滞不确定性【作者】刘彦文;王广雄;綦志刚;许保同【作者单位】哈尔滨工程大学自动化学院,黑龙江哈尔滨150001【正文语种】中文【中图分类】TP2731 引言(Introduction)采样控制系统是指用一离散时间的控制器,例如数字计算机,控制一个连续时间的对象,系统中既包含有连续的动特性,也包含有离散时间的动特性.这里讨论的时滞是指对象中存在的时间上的滞后现象,这是在过程控制领域中经常存在的一种信号或能量传递滞后的现象.目前大多数的理论工作都是将这时间上的滞后用离散系统的概念来处理,即将时滞τ看成采样周期h的倍数.但是过程控制中对象的时滞并不一定等于h的整数倍,且带有一定的不确定性.文献[1]指出,整数倍时滞稳定的离散系统,当实际的时滞τ与整数倍有差别时,实际的采样系统有可能是不稳定的.这里的问题是这个τ是一种模拟量之间的滞后关系,不是离散的信号之间的关系.文献[1]首次指出了这种采样系统中连续时间时滞不确定性的鲁棒性问题,不过文献[1]是用小增益定理来进行处理的.小增益定理因为没有包含相位信息,所以具有一定的保守性.文献[1]的主要贡献是提出了一种使权函数尽量贴近不确定性的界的做法以减少保守性. 从研究方法来说,离散系统的时滞分析一直是许多文献关注的热点,例见文献[2-3]及其后所附的文献.但是采样控制系统与离散系统不同,采样系统的对象中的信号都是连续时间的,因此采样系统的性能与纯离散(时间)系统是不一样的.对于含有连续信号的采样系统的性能分析,曾提出过提升技术[4-5].但是提升技术在控制系统设计中是有局限性的[6-7].为了研究采样系统中连续信号的性能,文献[8-9]提出了采样系统频率特性的概念.频率响应法在控制系统稳定性分析中占有重要地位.这是一种图解解析法,基本上不受阶次的限制.本文将根据文献[10]中采样系统的连续信号之间频率特性的概念,直接用连续信号之间的频率特性来分析这种时滞不确定系统的鲁棒性,方法简单实用.这里要说明的是,修正z变换也可用来分析采样系统中的时间滞后.采样系统中两个信号之间如果存在时间滞后δ(δ<h),那么滞后信号采样所对应的z变换称为修正z 变换[11].理论上修正z变换是可用来分析非整数倍时间滞后的,但计算复杂,一般只用于低阶系统,例如二阶系统的分析[1].本文中将采用修正z变换来进行配合性的说明和验证.2 时滞不确定系统的稳定性分析(Stability analysis of the system with uncertain timedelay)设所考虑的连续对象为式中:P0(s)为一稳定的有理传递函数,τ为滞后时间,图1所示就是所研究的时滞采样控制系统.这里将时滞τ按采样周期h的整数倍部分和余下部分分开,整数倍时滞与连续部分的P0合为一个对象特性P,而余下部分e-τu s代表了模拟量信号滞后时间的不确定性.图中:H为保持器,S为采样器,为离散的数字控制器,前后的开关是为了强调控制器前后的信号都是离散的,u 为连续对象的输入,z为对象时滞部分的输入,y为对象的连续输出,w信号表示时滞环节前后的信号差.图1 时滞采样控制系统Fig.1 Sampled-data control system with time-delay 现将图1的系统分拆成图2的形式.这是将时滞部分分成并行的两个通道:一个直通通道和一个包含时滞特性的通道(e-τu s-1).这样,系统就由上下两部分组成,下半部的输入信号为w,输出信号为z.注意到在采样控制系统中w和z都是连续信号,不过在w到z的连续信号的回路中还包含有离散信号.图2中w信号的含义可以从图1上来说明.图1中的虚线部分表示了这个w信号的组成,可见这个w信号就是时滞环节前后的信号差.图2中将系统在w和z处进行分开,所考察的就是这个时滞前后信号差(w)的动态特性.如果这个w(t)能收敛到稳态值0,就表明这个时滞系统是稳定的.由于时滞环节这个连续时间的动特性与系统其他动态部分的相互作用,时滞环节前后的信号差有可能是发散的,即w可能发散.这就是采样系统中连续时间时滞不确定性影响系统稳定性的原因.图2 时滞采样控制系统分析用的框图Fig.2 Block diagram used for analysis of sampled-data control systems with time-delay下面的推导中要用到一些标准的表达式:用*号表示采样信号表示采样信号的拉氏变换:其中:Y*(jω)表示其频谱,h为采样周期,ωs=2π/h.现在来计算图2中从w到z的频率响应特性[10-11].这一回路中各信号的拉氏变换式为式中H(s)为保持器,式(7)中括号部分(P H)*表示P(s)和H(s)相乘后再离散化.根据式(6)-(8)可得输出信号的拉氏变换式设输入信号是一正弦函数w(t)=exp(jω0 t)[9],并设ω0<π/h.这种函数也称为复数正弦(phasor),其频谱为注意到这里ω0<ωs/2,并不存在频率混叠现象,当只研究主频段(-ωs/2,ωs/2)上的特性时,W*(jω)与W(jω)是一样的.因此对正弦输入信号来说,可以将式(10)中的W*换成W,从而得输出对输入的频率响应特性为式(12)中的负号反映了控制器的负反馈作用(见图2).现将该负号单独提出,并用T zw(jω)来定义这部分的频率响应,即定义这里的分析中要求T zw(jω)是稳定的.而从图2可以看到从w→z的这一回路是整数倍时滞采样系统,其稳定性很容易用常规的离散化设计来保证.在正弦信号的假设下,图2可简化成图3的形式,图中D(jω)为图3系统的反馈连接Fig.3 The feedback connection of the system图3 是一种负反馈连接.根据频率法可知,如果即系统的Nyquist图线经过-1点,这时系统就是临界稳定的.将式(15)改写如下:式(16)的左侧与时滞τu无关,而右侧则只与时滞的参数τu有关,根据二者的相对关系就可以判断系统在此时滞下的稳定性.注意到-1/D(jω)的图形是非常简单的.根据式(14)可写得式(17)表明,D的负倒特性是一条在实轴0.5处平行于虚轴的直线,实轴以下的一段直线对应于ωτu从0→ π.由此可见,式(16)左右两项的交点在第2象限.-1/D的直线是由下往上,而频率特性T zw(jω)的走向(ω增加方向)是由右向左.设二者相交时的时滞为τuc,频率为ωc.当系统的时滞τu>τuc时,-1/D上的点将处于频率特性频率增加方向的右侧.这个-1/D相当于频率法中的-1点,当-1点在频率特性的右侧时,系统不稳定.也就是说,若系统的时滞大于τuc时,该采样系统是不稳定的.τuc就是鲁棒稳定的上限.如果系统的频率特性T zw(jω)在进入第2象限时其实数部分已小于0.5,那就不会与-1/D线相交,就不会因为有时滞τu而不稳定.这就是时滞无关稳定性.当然式(16)只是正弦周期解的条件,如果不稳定时的波形与正弦型出入较大,那么计算结果是会有误差的(见算例).3 算例(Examples)算例1 设一单位负反馈的采样控制系统,其连续对象为[1]式中:τ为时滞时间,见式(2)所示,ς为阻尼比.本例中设采样周期为如果τ为整数倍时滞,τ=vh,此时对象特性为根据常规的离散化方法可以知道,当采样周期h为式(19)时,式(20)的z传递函数P(z)=0[1].这相当于系统开路,但因为对象是稳定的,在单位负反馈控制下(图1),这个系统显然是稳定的.而且这个系统在任意整数v下都是稳定的.即在离散(时间)的概念下,这个单位负反馈的闭环系统是时滞无关稳定的.但是如果τ与整数倍时滞vh 有差别,这个采样系统就有可能是不稳定的.本例中式(19)的采样和P(z)=0,属于病态采样[1],但是因为式(20)比较简单,可以用解析的方法来进行分析,所以文献[1]用这个例子来说明采样系统中存在时滞的鲁棒性问题.但因为是病态采样,文献[1]的方法最后并没有用于这个例子.本文的方法则不受病态采样的限制,而且这个例子确有其特殊之处,通过这个例子还可进一步说明本文方法的适用条件.具体计算时,本例中设整数v=1,即对象的时滞为即图1中含有整数倍时滞的对象为根据式(19),设本例中的采样周期h=3.3 s,ς=0.3061.本例为单位负反馈,即=1.将式(22)代入式(13),并注意到本例中的P(z)=0,即(P H)*=0,得图4所示即为所得的频率响应T zw(jω).当ω=0.5 rad/s时,T zw(jω)与-1/D线相交.根据交点处的座标,从式(14)可得对应的τu=1.12 s,说明该采样系统的时滞当超出采样周期h的值达到1.12 s时就会失去稳定性.图4 时滞采样系统的稳定性判别Fig.4 Stability test of the sampled-data system with time-delay现在对所得结果进行验算.利用修正的z变换(modified z-transform)公式[11],根据式(1)-(2)可得图1系统中对象的z传递函数为本例中ς=0.3061,h=3.3 s,v=1.根据式(24)可写得本例中的闭环系统的特征方程式.当τu=1.28 s时,得该特征方程式为式(25)表明,z平面上的一对特征根正好超出单位圆,说明时滞大于1.28 s时系统就不稳定了.图5就是τu=1.28 s时,按图1的系统结构所得的混合仿真结果.仿真时的初始条件是z(0)=1,系统在这组参数下刚开始要发散,与式(25)的特征方程式的分析结果是一致的.上面图4用图解解析法求得的系统鲁棒稳定的时滞上限是τu=1.12 s,而实际的上限是1.28 s(式(25)).误差的原因是这个病态采样系统w(t)的波形与正弦型有一定差别(见图5).这说明,如果波形较差,上面的图解解析法可以提供一个定性分析的结果,如果波形接近正弦型,那么这个方法就可给出一个定量的结果.图5 时滞τ=3.3+1.28时的调节过程Fig.5 The transient response for time-delayτ=3.3+1.28算例2 本例是一个正常的采样系统.设图1中的连续对象为并设采样周期h=1 s.与算例1类似,根据式(26)可算得在单位反馈(=1)作用下的频率响应T zw(jω)(式(13)). 该T zw(jω)曲线与-1/D线相交处的参数为ω=1.77 rad/s,τu=0.56 s.式(26)是比较简单的,故可求得其修正的z变换式,并进而求得在这个摄动值τu=0.56 s时闭环系统的特征方程式为式(27)表明,该系统的一对特征根正好超出单位圆,与上面图解解析法所得的结果是一致的.图6所示就是该系统在这组参数下的仿真曲线.这里只给出w(t)在开始要发散的前40 s的图形.由于波形接近正弦,所以分析的结果比较正确.本例属于正常设计,表明本文的图解法可用于采样系统时滞鲁棒性的定量分析.图6 算例2的响应曲线w(t)Fig.6 The response of w(t)for example 24 结论(Conclusions)对于有时滞的采样控制系统,如果采用离散化设计,只能保证整数倍时滞时的稳定性.而连续对象中的时滞并不正好等于整数倍的采样周期,实际上的时滞值相对于整数倍的摄动会使系统失去稳定性.文中定义的信号w(t)反映了这个摄动影响的动态特性.利用采样系统的频率响应特性可以定量地确定鲁棒稳定性的时滞不确定性的上限.即使当周期解的波形偏离正弦型时,本方法仍可以对系统的鲁棒性做出相当有效的定性判断.本文提出的用w(t)信号来分析时滞不确定性的方法,对于有时滞的一般连续时间系统也是适用的.参考文献(References):【相关文献】[1]ALTERMAN I,MIRKIN L.On the robustness of sampled-data systems to uncertainty in continuous-time delays[J].IEEE Transactions on Automatic Control,2011,56(3):686-692. [2]LIX,GAO H.A new model transformation of discrete-time systems with time-varying delay and its application to stability analysis[J].IEEE Transactions on Automatic Control,2011,56(9):2172-2178.[3]LIU J,ZHANG J,HE M,et al.New results on robust H-infinity control for discrete-time systems with interval time-varying delays[J].Journal of Control Theory and Applications,2011,9(4):611-616.[4]BAMIEH B A,PEARSON J B.A general framework for linear periodic systemswith applications to H∞sampled-data control[J].IEEE Transactions on AutomaticControl,1992,37(4):418-435.[5]CHEN T,FRANCIS B A.H∞-optimal sampled-data control:computation anddesign[J].Automatica,1996,32(2):223-228.[6]WANGG X,LIUY W,HE Z,etal.The lifting technique for sampleddata systems:useful or useless?[J].Acta Automatica Sinica,2005,31(3):491-494.[7]WANG G X,LIU Y W,HE Z.H∞design for sam pled-data systems via lifting technique:conditions and limitation[J].Acta Automatica Sinica,2006,32(5):791-795. [8]刘彦文,王广雄,何朕.采样系统的频率响应和L 2诱导范数[J].控制与决策,2005,20(10):1133-1136.(LIU Yanwen,WANG Guangxiong,HE Zhen.Frequency response and the L 2-induced norm of sampled-data systems[J].Control and Descision,2005,20(10):1133-1136.)[9]YAMAMOTO Y.KHARGONEKAR P P.Frequency response of sampled-datasystems[J].IEEE Transactions on Automatic Control,1996,41(2):166-175.[10]王广雄,何朕.应用H∞控制[M].哈尔滨:哈尔滨工业大学出版社,2010.(WANG Guangxiong,HE Zhen.Applied H∞Control[M].Harbin:Harbin Institute of Technology Press,2010.)[11]FRANKLIN G F,POWELL J D,WORKMAN M.Digital Control of Dynamic Systems[M].3rd editon.Beijing:Tsinghua University Press,2001.。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
o. 4
WANG Guang-Xiong et al.: A New Approach to Robust Stability Analysis of · · ·
511
(3)
or, in operator notation[1,2] z ˆk (t) = Φc xk + Φ11 w ˆk where Φ11 is the convolution operator, Φ11 : L2 [0, τ ] → L2 [0, τ ], and Φc is the state transition operator, Φc : Rx → L2 [0, τ ]. (The notation Rx means x is the dimension of the signal x). Φb and Φ12 in (2) are also the related operators, and B2d maps u(k) to the discrete-time state xk , where u(k) is produced by the zero-order hold. It is the B2 matrix that is formed from the discretization with hold, and is written with a subscript d. In control problems, (1) represents the generalized plant, the second input of the generalized plant uk and the second (sampled) output yk are connected with the discrete controller K (Z ). Let ˆ and K is K (z ) = Ck (zI − Ak )−1 Bk ; then the corresponding closed-loop system formed from G ˆ K) = Fl (G,
∞
ˆ f k
k=0
2 L2 [0,τ ]
<∞
Consider a time-continuous system as follows. ˙ (t) = Ax(t) + B1 w(t) + B2 u(t), x z (t ) = C 1 x (t ), y (t ) = C 2 x (t ) (1)
The relationship between the discrete-time state of the system and the lifted signals should be described by operator equations. Realization of the lifted system in operator form is given as ˆ= G
Acl Bcl Ccl 0
=
Ad B2d Ck Φb Bk C2 Ak 0 Φc Φ12 Ck 0
(4)
where Acl is a matrix, Ad = eAτ is the discretized state matrix of the plant, and Bcl and Ccl are operators. Notice that in forming (4), the operator Φ11 is assumed to be zero. It is true for the robust stability problem (see later). (4) shows that in the lifted system, the operator Φb maps the lifted input signal {w ˆk } to the ˆ . This discrete-time state and the discrete-time output of the discrete-time state xg (kτ ) of the plant G controller uk are mapped to the lifted output {z ˆk } by operators Φc and Φ12 , respectively. 2 2 ˆ : lL Let the system operator be Σ : L2 [0, ∞) → L2 [0, ∞), and its lifted be Σ → lL . It is 2 [0,τ ] 2 [0,τ ] [1] ˆ . proved that the system norm is preserved after lifting, i.e., Σ = Σ ˆ into a matrix one The last step of the lifting design is to transform the operator realization G Gdd , which is also called the H∞ discretization (it means that the H∞ norm is preserved)[3] . Then the general method for discrete-time system can be used to design the system. Now consider the robust stability problem of Fig. 1, where P is the plant, K is the discrete controller, H is the zero-order hold, S is the sampler, F is the antialiasing filter, and W is the weighting function of the multiplicative uncertainty.
eAτ Φc C2 Φb Φ11 0 B2d Φ12 0
(2)
The following equation involving Φ11 is given here as an example of these operator equations. Consider only the input w. The relationship between the lifted output {z ˆk } of system (1) and the lifted input {w ˆk } can be given as z ˆk (t) = C1 eAt xk + C1
(Department of Control Science and Engineering, Harbin Institute of Technology, Harbin 150001) (E-mail: gxwang@) Abstract The lifting technique is now the most popular tool for dealing with sampled-data control systems. However, for the robust stability problem the system norm is not preserved by the lifting as expected. And the result is generally conservative under the small gain condition. The reason for the norm difference by the lifting is that the state transition operator in the lifted system is zero in this case. A new approach to the robust stability analysis is proposed. It is to use an equivalent discrete-time uncertainty to replace the continuous-time uncertainty. Then the general discretized method can be used for the robust stability problem, and it is not conservative. Examples are given in the paper. Key words Sampled-data system, lifting technique, robust stability, small gain theorem
ˆ The sequences {f k } are discrete-time signals which take values in the function space L2 [0, τ ]. Let 2 ˆ {fk } ∈ lL2 [0,τ ] , i.e., an L2 [0, τ ]-valued sequence whose norm sequence is square integrable[1] ,
0 t
eA(t−s) B1 w ˆk (s)ds
1) Supported by the Harbin Institute of Technology Fund for the Key Subjects (54100179) Received May 14, 2004; in revised form December 28, 2004
1
Introduction The controller in a sampled-data system is discrete, but the input and output signals of the plant are continuous-time signals. So the lifting technique has become the first choice for sampled-data system analysis and design in recent years. As for the robust stability problem, since the perturbation of the continuous-time plant is also continuous-time, the system norm between the corresponding continuoustime signals should be considered when using the small gain theorem. The robust stability problem is also one of the two reasons for introducing the lifting technique[1] . But it is just this robust stability problem where the system norm is not preserved by the lifting. In this paper the problems and the conservativeness of the lifting technique are discussed, and a new and not conservative method is proposed. 2 Problems of lifting design The lifting can be visualized as breaking up at each sampling time the continuous-time signal f (t) ˆ into an infinite number of consecutive pieces f k (t ) ˆ f k (t ) = f (τ k + t ), 0 t τ