对控制系统随机性能分析进行性能评价怎么估计过程的滞后时间

拒绝访问 | www.bzfxw.com | 百度云加速
请打开cookies.
此网站 (www.bzfxw.com) 的管理员禁止了您的访问。原因是您的访问包含了非浏览器特征(3ef47d-ua98).
重新安装浏览器,或使用别的浏览器如何评价自动控制系统的性能?_百度知道
如何评价自动控制系统的性能?
我有更好的答案
稳定状态下和期望值的误差,越小越好。响应时间:系统对于输入信号的反应时间,响应时间越快越好。超调:在过渡过程中和期望值最大的差值。超调量越小越好。稳态误差通常来说评价自动控制系统的性能有下面三个指标
采纳率:37%
那就要看你试验的是什么
和你的各项指标
为您推荐:
其他类似问题
自动控制系统的相关知识
换一换
回答问题,赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。火电厂热控系统控制性能评价方法的研究_甜梦文库
火电厂热控系统控制性能评价方法的研究
华北电力大学(保定) 硕士学位论文 火电厂热控系统控制性能评价方法的研究 姓名:张彩 申请学位级别:硕士 专业:控制理论与控制工程 指导教师:王印松
华北电力大学硕士学位论文摘要摘要以Swanda提出的控制器跟踪性能评价方法为基础,研究了基于无量纲性能指 标的单回路PID反馈控制系统的设定值跟踪性能评价方法以及热工模型中纯迟延 时间的递推最t]、--乘辨识算法;在一阶惯性加纯迟延模型的基础上,以二阶惯性加 纯迟延模型为被控对象,以内模控制原理来整定PID控制器参数为基础,得到了基 于实际阶跃响应数据的无量纲性能指标基准,对控制器的性能划分了若干等级;采 用这种方法对蒸汽温度和压力控制系统的性能进行评价,仅利用闭环阶跃响应数据 计算动态性能指标并辨识得到纯迟延时间,利用纯迟延时间对性能指标无量纲化, 得到性能评价指标,对照本文制定的性能评价等级给出控制器性能的优劣。关键词:性能评价,无量纲性能指标,最小二乘辨识,内模控制ABSTRACTBasedassessORthe methodology Swanda proposedadevelopedmethodologyis studiedon atothe perfoI芏】啮nce tracking setpoint of PID controller,which is basedsetofdimensionless performance indicesand锄on-lineidentification basedto aOnstandard recursiveleast squares,for the estimation of pure time delay.With regardsctpointareramp signalchange from closed-loop response,the dimensionless performance indices model of firstorbasedOnasecond order rational transfer function plusatime delay.The satisfied to define performanceperformance criteria got from internal mode classescontroltuning Rreusedandgive the band of the indices for every class.Performance of steam temperatureand pressure calculatedsystem wasasses,溉],fromclosed-loop data performance indices Wereis used to make traditional indicescarlandpure time delay wasidentified,whichdimensionless.According to the performance band in the paper,performancebe assessed.ZhangCai(Control Theory and Control Engineering)by Prof.Wang YinsongDirectedKEYWORDS:performanceassessment,dimensionlessperformanceindices,recursive leastsquares identification,intemal mode control 华北电力大学硕士学位论文摘要摘要以Swanda提出的控制器跟踪性能评价方法为基础,研究了基于无量纲性能指 标的单回路PID反馈控制系统的设定值跟踪性能评价方法以及热工模型中纯迟延 时间的递推最小二乘辨识算法;在一阶惯性加纯迟延模型的基础上,以二阶惯性加 纯迟延模型为被控对象,以内模控制原理来整定PID控制器参数为基础,得到了基 于实际阶跃响应数据的无量纲性能指标基准,对控制器的性能划分了若于等级;采 用这种方法对蒸汽温度和压力控制系统的性能进行评价,仅利用闭环阶跃响应数据 计算动态性能指标并辨识得到纯迟延时间,利用纯迟延时间对性能指标无量纲化, 得到性能评价指标,对照本文制定的性能评价等级给出控制器性能的优劣。关键词:性能评价,无量纲性能指标,最小二乘辨识,内模控制ABSTRACTBasedassess011the methodology Swanda Woposedadeveloped methodology is studied011 atothe perfoanance tracking sctpoint of PID controller,which is basedsecofdimensionless performance indicesand锄onqincidentification basedaOilstandard rccursiveleast squares,for the estimation of pure time delay.With regard to change from closed-loop rcsponse.tlc dimensionlesssetpoint ramp signalafeperformanceaindicesbased011 amodel of first 01"second order rational transfer function plus performance criteria got矗Dm internal mode controlclasses and give the and pressuretime delay.The satisfied to defme performancetuning黜usedbandof the mdices for every class.Performance of steam temperature closed-loop data performancesystem WaS assessed,fromindiceswerccalculatexl and pure time delay was identified,which isusedtomake traditional indicesdimensionless.According to the performance band in the paper,performance Can be assesse正zllallgCai(Control Theory and Control Engineering)byDirectedProf.Wang YinsongKEYWORDS:performanceassessment,dimensionlessverformanceindices,recursive leastsquares identification,intemal mode control 声明本人郑重声明:此处所提交的硕士学位论文‘火电厂热控系统控制性能评价方法的 研究》,是本人在华北电力大学攻读硕士学位期间,在导师指导下进行的研究工作和取 得的研究成果。据本人所知,除了文中特别加以标注和致谢之处外,论文中不包含其他 人已经发表或撰写过的研究成果,也不包含为获得华北电力大学或其他教育机构的学位 或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均己在论文中作了明确的说明并表示了谢意。学位论文作者签名:.缢丝日期;a缒』:丝关于学位论文使用授权的说明本人完全了解华北电力大学有关保留、使用学位论文的规定,即:①学校有权保管、 并向有关部门送交学位论文的原件与复印件;②学校可以采用影印、缩印或其它复制手 段复制并保存学位论文;③学校可允许学位论文被查阅或借阅:④学校可以学术交流为 目的,复制赠送和交换学位论文;⑤同意学校可以用不同方式在不同媒体上发表、传播学位论文的全部或部分内容。(涉密的学位论文在解密后遵守此规定)作者签名:裂塞El导师签名:至翌丝日期:丝:£f£期:丑点丛 华北电力大学硕士学位论文第一章引言1.1选题背景和意义随着工业自动化水平的提高,工业控制系统的规模越来越大,一个大型的工业 控制系统有成百套控制回路,过程控制技术面临的重要任务之一就是对控制系统的 监视和维护。在控制系统的使用过程中,有很多原因会导致控制器性能下降,比如设备的老化、非线性的作用等,为了确保控制器能在设计的性能下工作,就需要对大量的过程数据进行分析来判断控制器在工作过程中性能如何。一个大的自动化系统控制回路众多,人工来分析数据并判断控制器性能的好坏工作量相当大,这种人工监视评价控制器的性能的方法难以适应当前的高度自动化水平。研究控制系统的 性能评价方法,形成与当前自动化水平相适应的控制系统性能评价体系,对实际的 控制系统进行监视,确保控制器在工作过程中性能良好有很重要的意义。 控制系统性能评价的主要目的就是对控制系统的性能作出定量的评定,评定出它的性能等级:优、良、中、可、差等,并对控制系统存在的问题提出早期的识别,它可以实时地指导过程控制工程师和工艺操作人员对控制性能的潜在问题采取各种应对措施,做到防患于未然,减少经济损失增强安全性。缺少有效的控制器性能评价与诊断工具,很可能要在问题暴露之后才能采取补救措施。因此虽然控制器性能评价这一领域的理论研究才刚刚起步,但却倍受工业界的关淀【l】。在这一理论的工业应用方面取得了一定的成果,并且针对一些特定的工业过程出现了一些性能评 价软件。 目前,已有许多成熟的控制理论和技术成功地应用于计算机闭环控制系统中, 如;模糊控制、PID控制、预测控制,鲁棒控制等等。对自动控制系统性能的主要 要求在时域中一般为:稳定。并有一定的裕量:符含要求的瞬态响应;符合要求的 稳态响应。当给定一个控制系统时,根据什么来判断是不是一个好的系统,与其它 控制系统怎样比较才好,以目前用得最多的PID控制为例,在PID控制中,PID参 数的整定是一个关键的问题,PID参数直接影响控制效果。应如何选择最合理的PID 控制参数,在选择合理的PID参数之前,必须判断现有的PID控制器的参数是否 选择合理,对控制系统的性能进行评价。确定PID参数传统的方法也是在获取对象数学模型的基础上,根据一定的原则来确定。因此在获得PID控制的参数,进行PID的整定和自整定之前,必须知道被控对象的模型,对控制系统的性能进行评价。 所以控制系统在线闭环辨识和性能评价方法的研究对于PID参数的整定和自整定 也是一个必不可少的环节. 目前在工业过程控制中普遍采用PID控制回路,在生产过程自动控制的发展历l 华北电力大学硕士学位论文程中,PID控制是历史最久,生命力最强,应用最广泛的基本控制方式。这是由于 PID控制器算法有其自身的优点:原理简单、使用方便、适应性强、鲁棒性强等, 然而据统计有60%左右存在着性能缺陷【2.3】,其中有些问题可通过参数的调整来解 决,而另外一些问题只能通过采用新的控制策略或改造硬件设备来改善控制器的性能。在火电厂的控制系统有其自身的特点,在自动化水平比较高的火电厂中,采用 DCS控制系统,运行数据很容易存储和获取,为我们利用历史数据进行分析和计算提供了丰富的资源;电厂热工系统中有上百套自动控制回路,控制品质的好坏影响着系统的安全性和经济性;控制系统基本上采用PI和PID控制器。针对PID控制 的热工系统提出量化的性能评价指标,给出系统的性能评价方法,利用丰富的数据 资源来得到实际控制系统的性能状况,这对电厂的安全经济运行,减少工程师和操 作人员的工作量以及迸一步提高电厂的自动化水平有很重要的意义。 本文主要研究了基于无量纲性能指标的控制系统性能评价方法。首先针对具有 纯迟延的一阶惯性环节的单回路反馈PI控制系统,采用内模控制原理整定控制器 参数,在此基础上利用纯迟延时间对阶跃响应的动态性能指标无量纲化,得到PI 控制器的设定值跟踪性能指标的基准值以及系统无量纲动态性能指标和频域指标 的关系,据此设置了量化的性能等级,给出各等级性能中指标的范围。并且在一阶 系统的基础上借助于MATLAB工具进一步得到二阶系统和输入实际阶跃信号(斜率 为有限值的阶跃信号)下的性能评价指标和性能等级划分。在实际应用中利用系统 的闭环响应数据辨识得到系统的纯迟延时间,同时计算得到所需要的性能指标,从 而得到无量纲化的性能评价指标,根据文中给出的PID控制系统性能等级划分原则得到相应的性能优劣级剐。1.2课题研究现状控制系统性能评价主要有两个分支,一个是基于最小方差控制的随机性能评价 基准的研究,另一个是基于设定值响应的传统性能指标为基准的评价方法的研究, 先前的研究主要集中在最小方差基准的研究上,然而两种评价方法各有自身的优缺点。1.2.1最小方差性能评价基准最小方差控制律就是使目标函数(输出时间序列的方差)最小。基于最小方差 控制律的性能评价方法的原理图1.1和实现方法如下: (1)需要估计得到系统的纯滞后时间d;(2)对于闭环控制回路输出y,迸行时间序列分析,建立从扰动6/,到输出只的时间序列模型乃=如+地。(3)估计得2 华北电力大学硕士学位论文到过程输出的最小方差仃知(4)利用系统输出数据计算过程输出的实际方差盯;(5)计算得到母q)=仃知,叮2?图1-1最小方差性能评价方法原理图尽管at是未知的,这里可以对输出”进行时间序列分析得到白噪声的估计值磊,。仃:是控制系统的实际输出方差,玖d)在【o'1】之间,且露@)的值越大意味着更好,更强的控制,但是,其它性能指标如:鲁棒性等可能不能满足要求。因此,需要迸 一步考察控制系统的其它性能指标如:满足约束条件、闭环带宽、调节时间和上升时间等。Harris[3】首先提出了用最小方差控制来评价单回路控制器性能的思想,并用时间 序列分析的方法找到反馈控制器的非时变形式,然后用这个非时变形式作为评价控制系统性能的一个准则。随后,Desborough和Harris[41,Kozub和Gareia[s]等利用统计方法提出了无量纲的性能指标,即系统的最小方差与系统的实际方差之比。1yler 和Momti[61把研究扩充到单输入单输出的非最小相位和有不稳定极点的情况。 Huangt7】推广到多输入多输出控制器的性能评价,包括非最小相位的情况。Ko和 Edgarfsl把最小方差准则应用到串级控制回路中。ThornhiU等f9J将这一思想应用到对 单输入单输出控制系统设定点跟踪的性能评价中。Kozub[1。】从实践的角度论述了这 种性能评价方法的工程应用经验。Qinlll】在总结前人研究成果的基础上探讨了随机 最优性能(最小方差性能),确定性性能和系统鲁棒性性能相结合的可能性。1.2.2基于传统指标的性能评价基准最小方差准则主要是评价控制系统的随机性性能,,而对于控制系统的确定性 性能,主要指如衰减比(及衰减率),最大动态偏差(及超调量),残余偏差,调节 时间(及振荡频率)等动态响应指标。Shinskey提出用最小化系统的绝对值误差积 分值,来作为一阶线性加纯迟延系统的性能评价基准;Astrom提出了一种对PID 控制器性能进行评价的基准,利用系统的阶跃响应的性能指标与Z-N整定法得到的 性能指标进行对比;Huang和Shah[121借助于用户定义的特定闭环动态特性,如调节 时间,超调等提出了一种相对于基准的实际控制回路性能评价技术。Ju和Chiut”1针 对SISO和MIMO系统提出了一种基于最大闭环对数指标厶一的监控方法;Wang3 华北电力大学硕士学位论文和Cldu!141在此基础上进一步提出用设定点变化代替延时反馈的测试方法,并针对 误差信号利用快速傅立叶变换为基础的频率响应辨识技术从频率信息中计算L。? Kammer等[15】针对需完成闭环测试的线性二次型(LQ)性能评价提出了无模型方法。 Swanda和Sebor9116l提出用设定点响应数据来估计标量的性能指标。标量选用 的是调节时间和误差的绝对积分时间。首先选用一阶惯性加纯迟延的模型来等效实 际的被控制对象,他提出对于调整好的控制系统各性能指标与纯迟延时间的关系很 相近而与实际的被控对象关系不大,因此仅利用有限的过程知识来对PID控制器的 性能进行评价。用纯迟延时间对系统的动态性能指标无量纲化,并划分了性能等级, 来对实际的系统进行性能评价,这种性能评价方法主要针对在实际中应用最广泛的 PID控制器。这种用设定值响应数据来评价系统的性能有以下优点:阶跃信号在实 际的系统操作中容易获取并且便于分析和计算,并且操作人员对于阶跃响应的各项 性能指标的意义理解比较深刻,便于快速对问题做出处理。1.3本文研究内容首先,本文研究了Swanda提出的基于闭环理想阶跃响应数据的PID控制器的 跟踪性能评价方法。利用一阶惯性加纯迟延模型在内模控制下,得到在理想阶跃响 应下的无量纲性能指标的基准值,并给出了一阶惯性加纯迟延模型在PID控制下性 能等级的划分方法。 接着,在一阶惯性加纯迟延模型的基础上,研究给出了二阶惯性加纯迟延模型 在内模控制下的性能评价指标的基准值.这里借助于MATLAB工具,做了大量的 实验,分析了性能指标之间的函数关系以及性能指标与控制器和被控对象参数之间 的函数关系,得到性能评价指标的基准值,并划分了性能等级;进一步研究了,当 输入信号为实际阶跃信号(斜率为有限值)时的PID控制系统的性能评价,通过实 验给出了,经斜率修正的性能评价指标,并给出了相应的性能等级。 性能评价指标是通过模型的纯迟延时间来无量纲化的,本文研究了一阶和二阶 惯性加纯迟延时间模型的递推最小二乘辨识方法,通过辨识可以得到连续模型的具 有小数部分的纯迟延时间。 利用所研究的控制系统性能评价方法和纯迟延时间的辨识算法,对火电厂喷水 减温控制的过热汽温系统以及燃料量控制的主蒸汽压力系统的主调节器跟踪性能 进行评价,利用闭环响应数据计算性能指标并辨识出纯迟延时间,进一步得到用于 性能评价的无量纲性能指标,得到控制系统的性能优劣。4 华北电力大学硕士学位论文第二章基于无量纲指标的控制性能评价基准对控制系统性能评价的研究主要集中在基于最小方差控制的性能评价基准,而 对传统的阶跃响应性能指标关注较少,本章主要介绍了基于IMC.PI控制(根据内 模控制原理整定PID控制器)的无量纲性能指标的PI控制器性能评价基准。被控 对象采用广泛的一阶惯性加纯迟延模型,在内模控制的基础上,整定出PID控制器 参数,将控制系统的性能指标用纯迟延时间无量纲化,利用无量纲的调整时间和绝 对值误差积分以及超调量作为性能评价的基准,为了评价控制器性能的好坏,利用时域指标和频域指标之间的制约关系划分了性能等级。通过这种方法,仅利用闭环阶跃响应数据以及对象的纯迟延时间这些有限知识,可以对控制器的性能进行评 价。2.1无量纲性能指标,传统的性能指标,如调整时间、绝对值误差积分等这些时域指标反映了整个 系统的性能好坏,并不能很好的反映出控制器的性能。系统的性能还受到对象本身 的参数限制,如果对象的特性不同,即使在最优的控制器参数下,系统所能达到的 传统的性能指标也有很大的差别,这样就很难用统一的标准来对不同系统的控制性 能进行评价。在内模型控制中被控对象的纯迟延时间是一个重要的影响因素【"】。本文采用无量纲的调整时间、绝对值误差积分以及超调量来评价控制器的性能。调整时间是超调量和衰减率的综合反映。另外,对于不稳定系统,调整时间是 无穷大的,这就表示控制器的性能不好。调整时间用纯迟延时间来无量纲化 乙=‘Iz,Z比较小就表示控制器的性能好,比较大就证明控制器的性能差。调整 时间描述了响应曲线上一个点,绝对值误差积分则反映了整个时间响应过程的特性丝彩=髓E/眈l『)。lrol为阶跃信号的大小.调整时间主要反映了控制的速度,而绝对值误差积分和超调量主要反映了控制的质量。上面利用纯迟延时间定义了无量纲的时域动态性能指标,在无量纲化的过程中 已经消除了被控对象影响因素,反映了控制器本身的性能优劣,可以作为性能评价 的指标。本文利用反映动态性能的无量纲化的调整时间Z、绝对值误差积分尉E和 超调量0S以及反映系统稳定裕度的幅值裕量以、相位裕量丸来对控制系统性能进 行评价。对于PID控制控制系统,性能指标基准值的获取可以采用多种控制器参数 整定方法,Z―N法、c―C法、IMC―PI等等,这里以内模控制为基准,得到各无量纲 性能指标值,作为性能评价的基准值。 华北电力大学硕士学位论文2.2被控对象描述用传递函数来表示实际的-r业对象是一个近似的过程,大多数有自平衡能力的 对象一般可以用一阶或二阶惯性加纯迟延的形式来表示。如下形式的六个模型,代 表了一大部分稳定过程对象。模型11Go)=焉模型21G(J)=磊丽ke-tpo_25is+万1)厂e-一e'41模型2 模型4G。)=丽e-J G鼢丽1模型3 G(砖=―(I模型5∞)=等挚模型2到模型5这四个模型,在形式上与模型11和模型21有很大差别,模型 4和模型5没有纯迟延时闻,模型3和模型5具有零点,其中模型3有超调而模型 5有反调,但是都可以等效为一阶或二阶惯性加纯迟延的形式。用辨识的方法可以 得到这后四个模型的近似形式,用形如模型11和模型2l的一阶和二阶惯性加纯迟 延的模型表示。参数辨识可以从简单阶跃响应曲线利用最d,--乘法得到(豫擤20l。 利用辨识程序计算得到模型2到模型5的等效模型传递函数如下: 模型12 模型13模型14 模型15G(力=五蒜G(J)2至云e丽--G◇2i蟊eG(D5蒜i-鬲-模型22模型23G(s)=五夏历季干丽 模型24G(力2i西石―e万---j丽 模型25G(s)2百话嚣i了丽G(D。石磊磊≯e了-西丽下面给出了模型2到模型5以及其各自的一阶、二阶惯性加纯迟延模型的动态 特性和频率特性的对比,如图:6 华北电力大学硕士学位论文图2-1模型2到5(一)、模型12到15(一~)以及模型22到25(一)的单位阶跃响应曲线对比Frequency(rad/蛳)modeI 5图2-2模型2到5(一)、模型12到15(一一)、模型22到25(一)的频率响应曲线对比 华北电力大学硕士学位论文图2.1给出了模型2到模型5以及它们各自的等效模型1l和模型21的单位阶 跃响应曲线,从图中得出模型2到模型5的动态特性可以由其等效模型11和模型 2l近似表示,并且模型2l能更准确地描述原模型的动态特性。图2.2给出了模型 2到模型5和各自的等效模型的频域特性对比,在穿越频率附近每个模型与其近似 形式有相似的频率特性。在小于穿越频率时,模型1l和模型21都可以很准确地来 估计出真实模型的频率特性;在高频阶段模型2、3、4、5与各自等效模型ll的特 性有很大的差别,但是其二阶模型2l可以比较准确地来估计真实模型的高频特性, 另外二阶模型也可以表示对象的振荡特性,所以二阶惯性加纯迟延模型具有更广泛的代表性,本文就采用这一模型来研究控制系统的性能评价方法。2.3内模控制原理PID控制器参数整定方法有很多种,本文针对一阶惯性加纯迟延和二阶惯性加 纯迟延对象,采用内模控制原理来整定控制器,在这种原则下得到,性能指标,如调整时间t、绝对值误差积分伍以、超调量OS、幅值裕量彳。和相位裕量≯。与控制器参数之间的函数关系,以及性能指标之间的相互关系,获取了性能指标的基准值和性能等级的划分。用内模控制原理来整定PID控制器参数(文中简称IMC-PID)的方法很流行,并且在实际应用中能够很好地调整控制器来满足工业性能。IMC-PID整定方法是基于一阶惯性加纯迟延模型的,并且闭环传递函数也可以近似为一阶惯性加纯迟延的 形式,所以闭环响应速度可以由IMC中的滤波常数来决定,这个滤波常数也是闭环 传递函数的惯性时间。性能评价中选用IMC―PID整定方法有以下原因,首先,整定 方法类似与Lambda整定方法,而lambda整定在实际中有很广泛的应用121l;其次, 控制系统的跟踪性能和鲁棒性是控制性能很重要的两方面,内模整定方法能综合考 虑跟踪性能和鲁棒性能,内模控制中的可调参数滤波时间常数直接决定了系统的响 应速度和反馈控制回路的鲁棒性;最后,选择合适的f,可以使系统阶跃响应的绝对值误差积分最小122]。内模控制能克服被控对象增益和时间常数变化对控制性能的影响[23,36】,介绍了 内模控制的原理,以及根据内模控制原理来整定PlD控制器参数的方法,给出了 PID控制器参数与内模控制器参数的关系,并且分析了无量纲性能指标和内模控制 器参数之间的关系。下面给出了内模型控制的原理框图,其中,G(s)为被控对象,0(s)为对象的模型,Gc(s)为控制器。将否(s)分为两部分,G(D=盈0短(D,其中G+(s)=e一,瑾O)是不包含纯迟延B1部分的传递函数。q(s)=Q(DFO),其中F(D是滤波器,通常采用一阶或二阶滤波器:以s)=÷或F∞=-―二―万。1 1 S 华北电力大学硕士学位论文图2-3内模控制系统图£图2-4内模控制系统等效结构图 当广义被控对象内部模型确定后,低通滤波器时间常数f。成为控制器唯一的可 调整参数,f。增加,系统的响应减慢,t减小,系统的响应加快。这样参数整定简 单,便于工程应用【24J。纯迟延采用的近似一般为:泰勒近似:F1*l一嚣;Pade近似口1≈―1-m―/2。 l+嚣/22.3.1一阶纯滞后对象的IMC.PID控制当被控对象为一阶惯性加纯迟延形式,纯而后耶节用一阶黍勒展采j丘似,耀凝 器采用一阶,从而单位反馈控制器的传递函数可以用一个理想的PI控制器来等效,滤波器的PI控制器来实现,c(砖=i石G丽As)数为GOLc④=面‘T而s+l=巧(1+毒);如果采用二阶滤波器,单位负反馈控制器可以用一个带 CG=石j面1廿-口,则闭环传递函数为羔r=否南,分母中譬迟r=≤≥‘竺≠?系统的开环传递函延采用泰勒一阶近似P,:1一岱,得到羔:一e--".,则J,(f):r(1一口―-(t_-t))。由此计算tS+1系统的动态性能指标。控制器的典型形式∞)=%器,开环传递函数为GD工=瓦vs而/2+li,,闭环传递函数为羔:,如果纯迟延采用一阶pade近似,采用一阶滤波器,单位反馈控制器具有PID(rs/2+1)e一。分母中纯迟延采用一阶pade近9 华北电力大学硕士学位论文似得到闭环传递函数詈=等1一则等-1+秀咄呻町),~?,f』+,Zf.2.3.2二阶纯滞后对象的IMC.PID控制如果被控对象为二阶惯性加纯迟延的形式,同理可得Gc(D=堡鼍譬铲,相应的单位负反馈控制器c(蛳传递函数为c鲈南?盟学业,可以用一个带滤波器的PID控制器来实现,系统的开环传递函数为G&0)=―s(r―d二+―2r。一+O系统的闭环传递函数为等=南则y妒?一e掣c-+等。2.4本壹小结在上述推导的基础上,我们得到一阶滞后,二阶滞后对象在内模整定PI、PID 控制下的开环传递函数和闭环传递函数,由此并根据性能指标的定义得出最优无量 纲性能指标,作为控制系统性能评价的基准值。进一步利用这些指标和控制器以及 被控对象参数的函数关系,来标定控制器性能评价指标并划分性能等级。本章主要介绍了基于无量纲性能指标的控制系统性能评价方法所用到的基本 理论知识,包括控制系统传统的动态性能指标和频域指标的定义,内模控制原理和 基于内模控制原理PID控制器的参数整定方法,以及被控对象的具有纯迟延的一 阶、二阶模型。给出了IMC.PID控制器中PID控制器参数和内模控制器滤波常数 之问的关系,以及阶跃响应时间序列与滤波常数之间的关系,为后续的性能指标计 算提供了依据,本文所研究的控制系统性能评价基准以及性能等级的划分都是在本章的理论基础上进行公式推导或者编程计算得到的。lO 华北电力大学硕士学位论文第三章基于理想阶跃响应数据的PID跟踪性能评价根据第二章中的基础理论知识,针对PI控制的一阶惯性加纯迟延模型的单回 路反馈控制系统,利用阶跃响应的无量纲动态性能指标,给出单回路PI控制系统 的设定值跟踪性能的评价方法;并推广到基于阶跃响应的二阶惯性加纯迟延模型的PID反馈控制系统跟踪性能评价。在利用内模控制原理整定PID控制器参数的基础上,给出控制器设定值跟踪性能的无量纲指标基准值,并设置了五个量化的性能等 级,给出各等级性能中指标的范围。对参数设置不同的PI控制性能进行评价,褥出正确的评价结果。不仅得出量化的性能等级,而且能够给出各个等级中的PI控 制器参数范围。3.1基于一阶模型的PI跟踪性毹评价3.1.1无量纲时域性能指标考=兰三,可以计算出系统的动态性能指标。调整时间‘和误差绝对值积分厨E:t t5’1‘=鬣::篡三翟m=』Iy(f)一,(t)ldt=,(‘+f)@-,(3.2)名:≥:{3≥+1(△=±5徇 14L-"+1(△:垃%)覃f。。.,, ”’7觑易2带2争+?协4,同样由单位反馈开环传递函数G。(s)推导出频率特性,模型11在内模控制下魄=∞2赢矿幅值和相位为:艘=IG&u国)I2石了i而1,o=甜g[60L(jco)]=一三一锄 华北电力大学硕士学位论文由arg【G&(’,眈)】=万,得到相位穿越频率为:纹2云,从而幅值裕量为:4:上:要盘+1) 4….4R(coD i哼¨)门。’ (3.5)由l魄媳)I-1,得到幅值穿越频率为:%2夏i万1面,从而相位裕量为:叱=万+。(%)=争刁1鬲幅值裕量和相位裕量的关系可以表示为:西m 3.1.3时域/翱域指标之闻的笑藜(3.6)2三(1一去?可以看出基于闭环阶跃响应的无量纲化性能指标仅是f。,f的函数,o越小系 统的动态性能越好,但是在推导过程中采取了一定的近似,计算一般要求其最小值 为被控对象的纯迟延时间f,由于在内模控制中,滤波常数f,的选择受被控对象的 影响,所以有一定的范围。Thomasson[25】把f,2T+f作为内模控制器的调整规则, 然而在被控对象的惯性时间远大于纯迟延时间时,这个条件对于控制器的调整过于 保守,使系统的响应速度比较慢,Fruehauf等1261指出L应该大于或者等于过程的纯 迟延时问f,Rivera等【22】认为滤波常数最小值为f,此时的等效PI控制系统的绝对 值误差积分值最小并且系统的超调量很小或者等于零,在本文中取f,≥f。 L=3三£+1(△=±5%)ff。≥f剧目=至+lffc≥f厶=考哮钉)fc≥f九=争赤姆=≥+o.667L≥f同时给出返儿个尢量绷即住li邑于|『栋相詈襻住宿林隅僵裕置、干甚位诠重天糸? z≥4(3-7)如=三睁+o.667).z≥4 九:互一――L 九一2―+0―.667协8,(3_9) (3.9)】J 24 £≥4由上述推导可以看出。参数不同的一阶惯性加纯迟延模型在IMC,PI控制器作 用下.所得到的最优菁量纲件能指标都是相同的。这磐指标仅反映了粹制器的性能 华北电力大学硕士学位论文好坏,在内模整定的最优PI控制器参数下,对于参数不同的被控对象都可以得到最小的无量纲调整时闻£=4。3.1.4性能等级划分 以无量纲的性能指标作为判断PI控制器性能的基准,为了评定PI控制器距离 所能达到的最好性能有多远,将控制器的性能量化为五个等级,根据[27]中提出的 等级划分原则,这里利用I以丸以及超调量将控制器的跟踪性能划分为优、良、 中、可、差五个等级。根据C来划分优良级和良中级,然后利用厶,办来划分良和 中,在利用超调量划分中和差。控制中鲁捧性和动态性能是相互制约的,过高的鲁棒性能会导致动态性能变差,实际中A。>10和九>80。也是没有必要的,当A,=10,丸=80时,根据公式可以得到相应的t=17.1,M易=6.4,这两个值用来划分中和 可两个级剐,用超调量10%来划分可和差级别。表3-1 PI控制器性能分类表笠堡优 良 中 可 差互£<=4.04.0<瓦<=8.4 8.4<Z<=17.1竺幸 宰 聿丝旦髓易<=2 2<三4幺<=3.47 3.47(髓易<26.37 厨易>6.37 觑日>6.37£>17.1 Z>17.1∞<=lO% ∞>10%表3-i中给出了各等级对应的性能指标的双值范围,对于某一个等级各个性能 指标必须满足此等级的范围。在实际的控制系统中,利用无量纲的性能指标对控制 器性能进行评价时,仅需要知道系统的t,05以及已把d这三个指标的值,当这三 个指标都满足等级划分中的范围时,就可以得到相应的性能等级。3.1.5仿真研究利用表3―1中的性能等级划分标准,针 对两个被控对象,在单回路PI控制下, 对控制系统进行性能评价。㈣=黑 ,G2∽=筹3.1.5.1Z-N和IMC-PI整定性能对比分别用Z-N法和IMC-PI整定方法,得到PI控制器的参数,在这两种最优的PI 控制下,利用单位阶跃信号的闭环响应数据计算得到无量纲的时域指标来评价控制 器的性能,利用开环传递函数计算出系统的相位裕量和幅值裕量作为验证指标,来 判断设定值跟踪性能较优的情况下,鲁棒性能是否满足工业要求。两个具有差异的 被控对象在同一种控制器参数整定方法下所得到的无量纲的性能指标很接近,这就 华北电力大学硕士学位论文说明了用对象的纯迟延时间来无量纲化的性能指标可以作为控制器的性能评价基准。表3-2 IMC―PI和Z-N参数整定的性能对比两个被控对象在同一种控制器参数整定方法下所得到的无量纲的性能指标很接近,这就说明了用对象的纯迟延时间来无量纲化的性能指标可以作为控制器的性 能评价基准。另外可以看出由IMc―PI整定的控制器具有较好的动态特性:超调量 小、调整速度快、绝对值误差积分小并且鲁棒性能好,有合适的相位裕量和幅值裕 量,这就充分体现了用最优IMc―PI控制得到的无量纲性能指标作为控制器性能评 价的基准值是准确、合理的。从表中无量纲的性能指标可以判断Z-N法整定的PI 控制器的设定值跟踪性能较IMC-PI整定法差。 3.1.5.2PI控制器性能评价 固定结构的控制器参数设置不合理,或者被控对象的参数发生了变化都会导致 控制性能下降。上述两个被控制对象在控制器参数发生变化时,利用MATLAB仿真 得到其闭环阶跃响应数据,根据数据计算得到无量纲的时域指标,来判断不同参数 下的控制系统跟踪性能优劣,参照性能等级划分原则得到相应的性能等级。下面给 出了两个对象在PI控制下的性能评价。表3-3对象I PI控制性能评价lI 华北电力大学硕士学位论文表3-4对象2 PI控制性能评价利用闭环阶跃响应数据计算得到无量纲的时域指标作为控制性能指标来进行 性能评价和等级判定,并对照系统的阶跃响应曲线进行对比,分析这种评价是否符 合实际。图3-1和图3-2是对象1和对象2对应的阶跃响应曲线,图中的曲线编号 与表中的编号相对应。t/s图3-1对象1在PI控制下的阶跃响应在对系统的动态性能进行分析时,经常采用阶跃响应曲线得到其动态性能指标 来判断这一系统的性能好坏。由阶跃响应曲线可以直观地看出曲线1超调量小响应 迅速性能较好;曲线2、4、7超调量递增,震荡较剧烈,响应速度减慢;曲线3、5、t,l图3―2对象2在PI控制下的阶跃响应6为慢爬行,响应速度递减。由无量纲性能评价指标得到的性能等级与由阶跃响应曲线分析得到的性能判断结果相符合。 华北电力大学硕士学位论文3.1.5.3性能指标与控制器参数 对象l和对象2以X。和Z为坐标轴。运用MATLAB工具画出各个性能指标的等 商线,置。和Z的取值范围以Z-N和IMC整定的参数值为边界。对象1取:0.1<K。<0.25,300<霉<600;对象2取:o.1<恐<1.6,5<写<7t5。这样不仅能够对参数发生变化的控制器的性能进行评价,而且可以给出当某个或某 几个性能指标达到最优时的PI参数。比如由图3-3可以得到Z<5,5<£<10等 区间内的置。和正的取值范围。.如果综合各个性能指标就可以求出不同性能等级中 的K。和正的取值。岛 模型I f的等高线Contour鄞 模型1lEA。的等高线1模型1 oS的等高线lines。of£for modelContour lines of皿■。for model 1Contourlinc¥of DSformodel lI, 模型2 r的等高线 Contour抽∞甜L formodel2I'Contour模型2lEAd的等高线 lines of上EHJ for model 2模型2 05的等高线 Contour lines of os for model 2图3-3 PI控制器参数和性能指标关系等高线对象l和对象2相对应无量纲性能指标的等高线的形状、变化趋势和取值范围 都很接近,说明这种用纯迟延时间来无量纲化的指标以及等级的划分可以作为统一 的基准来评价PI控制器的性能,不仅可以评价出控制系统性能的好坏,而且可以给出不同性能等级的控制器参数范围。 华北电力大学硕士学位论文3.2基于一阶模型的PID控制器性能评价3.2.1无量纲时域性能指标 由第二章可知,一阶惯性加纯迟延模型,在内模控制下,如果采用一阶滤波器, 纯迟延采用一阶pade近似,内模控制器的等效单位反馈控制器可以用一个理想的 PID控制器来表示,并且系统闭环传递函数为:r手=等e一,则孚小秀咄却-f)饥fJ+lr Zf.i g- g- ZZc如果f。≥f,2,则阶跃响应没有超调,Rivera指出在IMC.PID参数整定中一般取‘≥f。则无量纲的调整时间很容易得到£=蔓=3.o生+丘ln(1一≠-)+1,得到调整时间的近似值:Z=3.0÷+o.5朋El≤L,彳≤4(3-10)2呱舢呱l(1-寺e-O-r)l*‘卜m…5¨弘三一=导+o.5f0,f21(3.11)3.2.2频域性能指标 与IMC.PI类似,IMC.PID控制开环的幅值裕量和相位裕量也可以表示为f,/,-的函数。模型11在IMC.PID控制下的开环传递函数为: 二s+1GoL=CG=。L_:一F“(f。+≥s幅值裕量和相位裕量分别为:A。=1.56二生+O.776 (3-12)(3‘ (3.13)妒。~u―1.21―46 丸~2―2r。―+l 可以得到幅值裕量和相位裕量的关系近似为九=三(1一号等与?3.2.3时域/频域指标之间的关系在上述给出的各个性能指标和L,f的函数关系,进一步推导得到其他性能指标 和无量纲调整时间的关系,由于用到f。,f21,所以这里L≥3.5。 华北电力大学硕士学位论文M易=子+o.333厶=O.52I",+O.516辅.5瓦23.5(3.14) (3.15) (3.16)九:i7/"一萼l萼L-rv3.2.4性能等级划分咖丁霸c玎5 t孔5与PI控制器性能等级的划分原则相同,我们得到PID控制器性能等级,采用 互=3.5,C--8.4,£=18.2作为优、良、中、可等级的分界点,并利用公式计算出分界 点目岛=1.5,尉易--3.13,乜易=6.4,可与差等级仍采用超调量来划分。表3-5 PID控制器性能分类表3.3基于二阶模型的PID控制性能评价在传统的工业过程中,二阶惯性加纯迟延模型能更准确地反映热工系统的特 性,并且用二阶惯性加纯迟延来近似高阶对象精度很高131】,并且这一点在第二章的 仿真研究中进行了证明。对于二阶惯性加纯迟延的模型,在PID控制下的设定值跟 踪性能评价,仍采用IMC整定PID控制器参数,来得到无量纲性能指标和控制器 参数之间的函数关系。只是由于二阶系统在公式推导上比较复杂,无量纲性能指标 和控制器参数之间的函数关系借助于MATLAB仿真实验来获取。3.3.1曲线的最小二乘拟合12列在科学实验的统计方法研究中,往往要从一组实验数据《z。,y;)中寻找出自变量工和因变量Y之间的函数关系Y=.厂(x)。由于观测数据往往不准确,因此并不要求Y=I厂(∞经过所有的点(r。,Y1),而只要求在给定点五上误差4=,(一)一弘按照某种标准达到最小,通常采用欧氏范数例12作为误差量度的标准。这就是所谓的最d,--乘法。在MATLAB中实现最小二乘法拟和通常可以采用如下两种途径:利用polyfit 函数进行多项式拟合;利用矩阵除法解决复杂型函数的拟合。 3.3.2一阶系统性能评价方法的程序实现 利用公式推导的方法,得到一阶系统时域和频域性能指标与控制器和被控对象 华北电力大学硕士学位论文参数的函数关系,现在利用仿真来获取闭环阶跃响应数据,通过最小二乘法拟合得 到f。,f作为自变量和性能指标的值作为因变量的函数关系。 下面是一阶PI、一阶PID和二阶PID系统的单回路闭环控制系统,其中的单 位反馈控制器是内模控制器等效的PI和PID控制器。图3-4一阶IMC-PI控制等效框图图3-5一阶IMC.PID控制等效框图图3-6二阶IMC.PID控制等效框图.性能指标£、尉易、以与%,f的函数关系,以及上4目、以与L的函数关系可以表示为y=甜+6的形式,丸一o,f以及九一L的关系可以表示为y=i旨+b的形式,其中L/f是自变量,性能指标是因变量,具体函数关系如下所示:L=q等+b-(A=±5哟;妣=口:詈+6:;以=码睾+毛;九=z等并+b?以L为自变量,其他性能指标为因变量的函数关系如下:己“乞=a5正+以;以=a6Ta+b6;丸=;;}了+b7 』,下1下面是一阶系统IMC?PI控制下的各个性能评价指标C,/A目,4,o。与‘,f以及£的函数关系曲线。19 华北电力大学硕士学位论文乙表3-6参数拟合结果4。6. 1.∞毛0.32a,屯1.00口,6|1.58 玛?108.0口.3.001.OO1.586.89.996,0.874. 0.536.1.37岛86.26图3?7 IMC―PI无量纲性能指真值(-?)和近似拟合值对比(一)其中(…)是通过程序计算得到的函数关系,相当于上一章中运用公式推导得到 的对应关系,(.)是利用最/b--乘法拟和得到的函数关系。表3-6给出了上述公式 中的各个参数,与第三章中推导的参数进行对比可以发现拟合结果很精确。 华北电力大学硕士学位论文同样,得到IMC.PID控制下的一阶惯性加纯迟延模型的无量纲性能指标: Z,M岛,4,①。与fc,f以及t的函数关系,曲线对比关系与上图类似,下面给出参数的拟合结果。 表3.7参数拟合结果 竺!3.00垒0.45垒1.00塾0.55竺1.58垒0.79竺1-52.9塾’墨0.32塾0.87堡0.53堡垒苎91.2594.530.57-101.33.3.3二阶IMC.PID控制系统性能评价 当被控对象为二阶惯性加纯迟延模型,在IMC.PID控制下无量纲性能指标的计 算相当复杂,这里就需要采用上述的性能指标计算方法和最小二乘拟合方法,得到性能指标与参数r。/r的函数关系,以及各个性能指标之间的函数关系? 华北电力大学硕士学位论文表3.8参数拟和结果口I4.75岛0.54 吼 0.41以1.94岛1.51 口‘ 0.42呜 1.98趣1.84 岛 ?70.96q -35.8屯80.686,1.29屯1.63岛76.18图3-8IMC.PID无量纲性能指标真值(一)和近似拟合值对比卜)一阶系统仿真曲线和推导的函数关系曲线能很好的吻合,结果表明,利用仿真程序得到的性能指标和f。,f的函数关系精确度比较高,可以用来求取二阶系统在不同的L,f下的性能指标,并进~步利用最小二乘法拟合得到函数关系。由图3-8可以看出对于二阶系统这种性能指标与参数之间函数关系的求取方法仍然适用,二阶 系统PID控制器性能指标函数关系如下所示:£=4.75等+o.54 舻1.98等+1.84/AEd:O.41r,+1.29九=80.68一羔彳。:O.42T,+1.63埘岛=1.94等+1.51丸:76.18一i7百0.96iJ-r1与一阶系统PI控制器性能等级的划分原则相同,我们得到二阶系统的PID控 制器性能等级,采用Z=5.3,C=lO.2,Z=20.1作为优、良、中、可等级的分界点, 并利用公式计算出分界点髓易=2.68,Ⅸ易=4.69,姐易--8.75,可与差等级仍采用超 调量来划分,具体的等级划分方法如表3-9所示。表3-9 PID控制器性能分类表3.3.4仿真研究针对二阶惯性加纯迟延模型G。(s)=K,P“其中%=l,乙l=5,‘:=10,乙=10。采用几种不同参数的PID控制器,仅利用闭环的单位阶跃响应数‘据得到控制器性能评价结果。22 华北电力大学硕士学位论文 表3-10 PID控制器性能评估l/see图3-9 PI控制回路性能评估阶跃响应利用MATLAB仿真得到相应的单位阶跃响应曲线,可以直观判断控制性能的 优劣依次为曲线2,l,6,5,7,3,4所表示的控制系统。与上述采用IMC?PID控制为基准的无量纲性能指标评价方法很好地吻合。3.4本章小结本章利用PI反馈控制系统的理想阶跃响应数据,来对PI控制器的设定值跟 踪性能进行评价,这种性能评价方法主要基于无量纲的调整时间和绝对值误差积分 这两个性能指标,对控制器的性能等级进行划分,在实际被控对象未知的情况下, 这种评价方法仅需要系统的闭环阶跃响应数据,辨识得到被控对象的纯迟延时间, 并计算出无量纲动态性能指标,与等级划分中的性能指标范围进行对比,判断出这 些指标所属的等级,从而得出控制器量化的性能等级,用来判断控制器距离调整好的控制器性能相差多远。 华北电力大学硕士学位论文第四章基于实际阶跃信号的PID控制系统性能评价上述几章在输入理想阶跃信号的基础上,研究了PID控制系统的设定值跟踪性 能评价的方法,给出了基于一阶、二阶惯性加纯迟延模型的性能评价指标的基准值, 给出量化的性能等级。然而在实际控制系统的性能测试或工况改变时。输入信号更 多地采用具有一定斜率的阶跃信号,本文称为实际阶跃信号,本章在输入理想阶跃 信号的基础上,研究了基于实际阶跃响应数据的PID反馈控制系统的跟踪性能评价 方法。4.1实际阶跃信号下的一阶系统性能指标利用对象GO)=壬,输入实际阶跃信号的斜率为k的实际阶跃信号,Ls十1k+=k/r,其中k’为标准化斜率取0到l之间的实数,根据系统实际阶跃响应输出数 据计算得到过渡过程时间&和误差绝对值积分HE。, t和烈易是当输入为标准 单位阶跃信号时系统的过渡过程时间和误差绝对值积分,经过大量的实验得到,以七。为自变量,乙和髓%为因变量的函数关系可以近似地表示为:Y=口+G一1)+b+ye属1(4-1)其中a着llb为待估计参数,Y表示实际阶跃信号下的无量纲性能指标,儿为理想阶跃信号下的无量纲性能指标,则乙、三t‰与k‘的函数关系为:珏=at‘(-专一1)+A+Z缸.1IA互☆=巴+(古一1)+62+髓易?J‘根据不同七‘对应的性能指标,采用最/b--乘法拟合得到待估计参数值,下面给出了瓦与k‘,髓‰与k’的实际(…)函数关系曲线以及由最小二乘法拟合(一)得到的近似函数关系曲线。k’ 图4_lL与七.的函数关系圈睨坂与矿的函数关系k‘ 华北电力大学硕士学位论文其中al=0.81,bl=0.32;a2=0.47,屯=0.48。函数关系可以近似的表示为:疋=z+o.75+嘻1一1)(4?2)“如=尉乜+o,45‘嘻一1)4.2实际阶跃信号下的二阶系统性能指标(4?3)对于二阶系统选择对象GQ)=j―圭一,采用式(4一1)所示的函数关系,利用 J’+S+I同样的参数求取方法得到aI=0.16,bl=--0.24;口2=O.11,b2=O.06;tl,3=0.18,63=_o.01。函数关系可以近似简化表示为:砭=乃4-0?75+咕一1),Ⅸ%=殂蜀+0?45‘嘻一1),OS,20S+O?17(k‘一妒丘 厍(4-4)k’||}’ /AE女与k‘的函数关系k’气与k‘的函数关系OSI与矿韵函数关系图4.3性能指标与k。的函数关系曲线(~实际,一拟合) 实际的性能指标与标准化斜率的函数关系曲线与由上述三个式子所表示的性 能指标与斜率的近似的函数关系曲线对比,如上图所示。由图可以看出当标准化斜 率由0到l时间,拟合得到的的函数关系对实际的曲线近似程度很好,不管是一阶 系统还是二阶系统输入信号的斜率与性能指标的函数关系都可以用相同的式子表示。4.3实际阶跃信号下一阶IMC.PI控制器性能评价在IMC.PI控制的一阶惯性加纯迟延的系统,如果设定值为斜率为k的实际阶 跃信号,则无量纲调整时间和误差绝对值积分与f√f的函数关系,根据上述公式进 行修正,将这种输入信号斜率和指标之间的函数关系,推广应用到无量纲性能指标 的相互关系中,下面给出具有k。值修正的性能指标与f,,f的函数关系计算公式。图4-4对比了在k’;0.1,0.2,0.5,0.8。l,1.5,2时,是性能指标乙一f。,f,以及埘艮一%/f的实际曲线(一)和利用修正公式拟合曲线(~)对比. 华北电力大学硕士学位论文F=o.I时瓦对比r哪.1时屿对比F=o.2lied对比F=o.2时L对比矿Io.8IAga对l:gFlo.s对瓦对比F=l对t对比矿tl/AEa对比F=I。5时L对比Ftlj时z4繇对比F=2对瓦对比F=2时五tEt对比图4.4不同斜率的实际阶跃信号下的性能指标与控制器参数关系曲线可以看出七’值修正公式对实际函数曲线的拟合关系很好,这种修正关系同样可 以应用到一阶IMC―PID控制系统和二阶IMC-PID控制系统性能评价中去。当 矿=l,I.5,2时,取实际阶跃响应性能指标的计算公式与理想阶跃响应性能指标的计算公式相同,其中幅值裕量和相位裕量保持不变与第三章相同。乙=3二‘+1+o.75(÷一1) 上t巨&=二£-+1+o.45(÷一1)由上述四式得到性能指标间相互关系:(4―5) (4-6)zd点I:―T,k―+―0.6=0―/k一"-I)+O.667(4.7)4=引型警业+o.舐,I九~瓦丽獗而_沥丽丽 九=三一酝面面万12面丽丽(4-s)‘4。9’ (4-9) 华北电力大学硕士学位论文4.4性能等级划分利用实际阶跃响应的输出数据计算系统的无量纲性能评价指标疋、/'肛at,根据信号的斜率矿值,将性能指标标准化为理想阶跃响应的无量纲指标,对照第三章中一阶IMC.PI、IMC.PID控制系统以及二阶IMC.PID控制系统性能等级的划分,就可以得到基于实际阶跃信号的IMC.PI、IMC.PID控制系统和二阶IMC?PID控制系统性能等级的划分原则,如表4一l表4-2表4.3所示.与基于理想阶跃输入信号的控制系统性能评价对比。基于实际阶跃信号的性能评价中,除了需要根据输出数 据计算性能指标和辨识出纯迟延时间外,还需要利用输入数据近似计算出输入信号的斜率。 表4-1一阶IMC.PI实际阶跃输入性能评价等级划分竺丝盘二旦:!!l!!生:二坠丝血二壁:塑!!墨:二里旦墨二Q:!!l墨:二12=优良五<=6.0 6.0<矗<=12 12<毛<=11.7五>11.7DIEd<=2.6 2.6<‘t日(=4.7 4.7</AEd<=6.4 /ffJ2d>6.4 IAEd>6.4车 + +中可 差08<_lo%08>1096矗>11.7表4-2一阶IMC.PID实际阶跃输入性能评价等级分类表笠堡互t二Q:!i丛丝:=坠丝墨.止=Q:堡!i!!墨:=堡垒旦二垒:12f!!墨:=壁优中 可 差五<=4.35 /AEd<=1.78 2.6<三t%<=3.23 3.23</AEd<=6.4 Lt妇)6.4 /AEd>6.4¥ 4 木良4.35<瓦<=8.78.7<£(=18.2 毛>18.2 毛>18.2OS<=10% 03>1096表4-3二阶IMC.PID实际阶跃输入性能评价等级分类表箜堡L二!:!i(!丝:=!)丝生’正=Q:垒壑!丛:二!)鱼}互二Q:!!(!』墨:=!):优良 中 可 差弓(=7.95 7.95<马<-lO.6 10.6<毛<=19.9 矗>19.9 /AEd(=5.25 5.25<‘t幻<=7.0 7.O</AEd<=8.2牛 + ?地>8.205<=1096 OS>10%矗>19.9‘锄>8.24.5仿真研究利用本章所褥到的基于实际阶跃信号的性能评价方法,我们针对典型的二阶惯 性加纯迟延对象。对控制器性能进行评价,这里借助于MATLAB仿真工具进行研 究。选取不同的PID控制器参数对其进行控制,输入分别采用单位理想阶跃信号和 实际阶跃信号,利用MATLAB工具得到响应曲线,并根据输出数据和性能指标的 华北电力大学硕士学位论文定义计算得到两种不同信号下的无量纲的性能评价指标,对输入实际阶跃信号的系 统进行性能评价时,将实际阶跃信号下的性能指标转换成性能评价等级中所需要的 指标值。G(s)=(10-―s+坚1)(15一s+1)表4.4基于理想阶跃响应PID控制器性能评估图4-5阶跃响应曲线 下面给出直接利用实际阶跃响应数据计算得到的无量纲性能指标,其中无量纲 性能评价指标是根据实际阶跃响响应指以及标准化斜率之间的关系转化得到的阶 跃信号输入下的无量纲性能指标。 袭4-5基于斜率F=O.5实际阶跃响应PID控制器性能评估 华北电力大学硕士学位论文表4-6基于斜率k‘=O.4实际阶跃响应PID控制器性能评估…印嚣。。.£堡叁型墨叁墼玉量塑丝鳇堂握L4F_dkOSk烈Tsk玉墨塑丝丝竖堕堂堡 傩 Ts执瞄性能 等级霉毫萼T型i翌Td―F』sk型IAE燮dk堡-OSk―jTa塑燮IAEd堂L05等级号砀1”表4-9斜率k’=l性能评价等级划分由以上数据可以得到,可以得出根据实际阶跃响应对系统的性能评价与根据理想阶 跃响应进行的性能评价相一致,这进一步证明拟合得到的基于阶跃响应的性能指标 和基于斜坡响应的性能指标之问的函数关系精度很高。 通过仿真得到实际阶跃信号下的动态响应曲线,将上述评价结果与曲线对比, 可以得出基于无量纲性能指标评价结果与直观的响应曲线相符合。 华北电力大学硕士学位论文k‘=D.5时的实际阶跃响应矿=O.4时的实际阶跃响应旷=O.2时的实际阶跃响应F=O.1时的实际阶跃响应圈4?6不同实际阶跃信号的动态响应曲线4.6本章小结在实际的控制系统中。为了检验控制系统性能好坏,很多时候输入信号采用斜 率为有限值的阶跃信号,或者在工况发生改变时,一些扰动信号也是类似于这样的 输入的信号,本章在前几章研究的基础上,对实际阶跃输入时的控制器性能评价方 法进行了研究,用仿真实验的方法给出了实际阶跃信号的斜率对性能指标以及性能 评价的影响,对性能评价公式进行了修正,对基于无量纲往能指标的控制器性能评 价基准做了进一步完善。在应用中只需要根据输出数据计算出性能指标值,再利用 纯迟延时间和实际阶跃信号的斜率转化为相当于理想阶跃响应的无量纲性能评价 指标,对应性能等级的划分得到控制器的优劣程度。 华北电力大学硕士学位论文第五章纯迟延时间的估计算法在利用无量纲的性能指标进行控制系统性能评价时,一方面需要利用闭环响应 数据计算出传统的性能指标,另一方面需要辨识出模型的纯迟延时间,利用纯迟延 时间对传统的性能指标无量纲化.所以纯迟延时间辨识算法的准确性直接影响着控制系统性能评价的结果。本章针对典型的热工模型――二阶惯性加纯迟延模型的参数进行在线辨识,介绍并推导了基于递推最小二乘在线辨识算法,对带有小数部分的纯迟延时间进行参数估计。5.1模型的变回归估计算法本章介绍了将二阶惯性加纯迟延的连续时间模型离散化,通过对离散模型进行 辨识来得到连续模型的参数。对于给定的采样周期,连续系统的小数部分迟延在离 散化时将会产生一个负实零点,因此可以离散零点引起的相位增益来估计小数部分 迟延,这就要求离散的零点只与小数部分迟延有关,但是连续系统的有限的零点在离散化过程中也将产生离散零点。这里我们需要辨识的典型的热工对象的离散零点只与小数部分迟延有关。对于离散化后的模型采用固定模型的变回归估计(fixedmodel variable regressor estimation)。FMVRE方法是Elnaggar【29】提出来的一种适合估计工业工程的一阶惯性加延迟 模型的纯迟延时间的估计算法。FMVRE方法比起其它的延迟的估计方法有许多优 点:应用起来非常简单和收敛非常快。FMVRE法是从估计其它的过程延迟参数里 推导出来的这使得它收敛起来非常快也很容易实现。具体的辨识算法如下:(I)直 接应用递归最/b-乘法和最小化损失函数辨识离散模型参数;(2)由辨识出的离散 模型得到连续系统模型,其中连续模型的小数迟延由离散零点计算得到;(3)此算 法可以将连续系统的有理传递函数以及整数和小数迟延一起辨识得到。5.2延迟的小数部分估计【30】对于一个采样数字系统,一阶惯性环节为如下形式:G(s);:鲁P1,T>015十I(5-1)其中,f=砌+巩h是采样周期,d是非负的整数,对于0<占<l,我们可以很容易的得到u和Y之间的脉冲传递函数掣:z一∥生堡!生坚,其中,p:exp(-h/D,utz) z―Pza-=p(p一一1)/(i―p),这样小数部分数延迟崩在z平面里可以产生一个在原点的极点和一个在负实轴的零点q=吨,(1一口),(0<口<1).同样,对于一l<占<0,我们3l 华北电力大学硕士学位论文有墨婴=z-dpl2【(1一∥弦+用其中∥=(p-l―p)/(1一p)和卢哼(1+占)当|Il。o,这样材IzJ z―P小数部分延迟动在z平面只产生一个在负实轴的零点叠;―∥馋一历’(o<夕<黟,现假设 延迟时问f是不确定的:f=砌±动,一1<F.<l,那么整个采样控制系统对于任何的 占∈卜l,1】可以表示成如下的脉冲传递函数:H(z1。z-(d“)垒±121±刍兰:。(5.2)z―p对于小的采样间隔,式(5-2)的零点(z‘,z‘)的关系如图5.1所示。群 鬻辜:叶.斗-J:叶.Jma m6 m4m2图5-1系统的零点z‘和2。的值在实际情况下,=’和:。总是在负实轴上,有如下的性质:仃净西+万。【1一功仃净了_三。+了j【_。;(1一D其中仃是频率响应◇扣一:。)0扣一2。)的幅角在p=0处的级数展开的第一项。 由上面的可以得到估计延迟的小数部分的递推辨识算法如下:‘5。3’ (5.3)a)用流程图所实现的RLS算法估计参数63,62,巩,p以及d的估计值寞,占2,蠢,声,0:”利用估计值63,62,丘通过式(5?2)计算得到z‘,z。,由(5.3)式可以计算出善=l一盯;c)如果一l<F<1回到a),否则停止算法。这时步骤c)可以被以下步骤d)取代;d)如果一l<善<1回到a),否则增加(减少)延迟的整数部分d的估计值和相应得减少(增加)善几个单位值以使I爿<l。这样修正估计参数蠡,龟,丘完成关系式(5.2)维持被估计模型的增益;最后回到a)。.系统的零点的求取利用方程毛z2+b2z+玩=O,所以在计算P时,我们不必求 出离散系统的零点,而利用参数bt,b2,以的值以及一元二次方程根与系数之间的关 系,直接计算得到F=(毛-bt)/b,+ba+b.),在参数估计中,小数部分迟延离散化 后可以统一设为6lz2+b2z+b3的形式,这样就可以在模型阶次已知的情况下进行离散模型的参数辨识,再根据F=慨一61),如+6:+61)得到连续传递函数的小数部分迟延。 华北电力大学硕士学位论文5.3整数部分纯迟延估计算法具有纯迟延时间系统的参数辨识,这里采用基于离散系统的最小二乘辨识算法。带有纯迟延时间的连续系统传递函数:G(D=G0(曲P1f=dhd=o,1,2,…其中G0(s)是关于s的分式多项式,f为纯迟延时间,h为采样周期,d为非负整数.对上述系统进行采样离散化,得到其脉冲传递函数G∽=舞玎,=1blz-(d矿+Db蕊2Z-(d+≮2)+?.鬲oq-bnzr-{d+n)A(z、7叫)l+珥z叫+口,z叫+…+口.。孑一”最/b--乘形式为:y(t)=妒(r)口+P(r),各个参数的估计值:乡(f)=≯(f归。其中,驴(f)=卜y(t-1),…一(f一甩),“(,一d-1),…uq―d―m)lr,0=∞t,口2,…a。,bl,b2,…6,r矿(f)=【-y(t-1),…,―】,(r一片),“0一d一1)’…uq-d一删7,0=[鼠,a2,…a。,岛,62,…kr损失函数采用以下误差的平方和的累积形式,各个参数的估计值就是使损失函数达到最小。e(t)=y(f)一多(,),‘,<口,D=∑p(f)】2,‘,矽,0)=nfm(g(e,d))面’当己知纯迟延时问在某一范围内【d。,d。】时,可以采用最小化损失函数的方法将延迟时间0。和其他参数向量百一起进行辨识。可以通过一个给定的迟延,反 复进行最小二乘估计最好的迟延就是能使损失函数最小的值。具体实现方法可以分 为两步来进行。 1)假设纯迟延时间已知,利用最小二乘法对其饱参数进行估计 t,(口,d,)=nfm(J(e,d))4∈【d。m,d。。】 2)对关于纯迟延时间的损失函数进行优化,得到纯迟延时间的估计值 ‘,(口,d)=min(2(O,4)) 这种两步算法在具体的递推程序实现过程中,只需在标准递推最小二乘算法其 中加入最小化损失函数的求取【3孤。首先假设上一步迟延的估计值是正确的,计算出输入输出的观测向量乒(七),利用标准最小二乘递推公式(5.4,5-5,5-6)得出其他参数估计值蚕(.|});然后假设其他参数的估计值毋(七)是正确的,让d从d。变化到d。,从而由式(5-7)得到d。一缸+1个观测向量,进而由式(5―8)得到d。一‰+1个损失函数值,第k步d的估计值就是使损失函数值最小的那一个,这个值又作为 定值用于下一步的参数估计中。下面给出改进的具有最小化损失函数的递推最小二乘算法的递推公式【虬I:≯(f)=卜寸仲一1)'...,..贝七一,磅,u(k一0一1),…u(k一0一m)f百(_|})=否(意一1)+x(七)【),(|j})一≯(.|})7舀(七一1)】(5.4) (5.5) (5.6)K(D=以.i}一1)乒(.j})眵7(七)P(.|}一D≯(七)+明一1 华北电力大学硕士学位论文P(七)=【,一足(七){;7(.|})】P(_|}一D/2 J(k,矗)=M(k一1,d)+眇<露)一多(毛矗妒(素)】2S(k,0)=nlinJ(k,力(5.7) (5,8) (5.9)w∈[d血,d。】方程(5-4)是输入输出的观测向量,方程(5-5,5-6,5-7)是标准的递推最 小二乘公式,在一设定的纯迟延时间的情况下对其他参数进行估计,后两个方程 (5-8,5―9)当纯迟延时间在某一个范围内变化时用来最小化损失函数并且求取使 损失函数最小的d的值,进而得到最佳的迟延时间。5.4热工对象的辨识方法5.4.I辨识算法 连续传递函数在离散化后,小数部分迟延将会产生在负实轴的零点。在进行参数估计时需要先利用输入输出数据进行脉冲传递函数的参数估计,再利用脉冲传递函数的零点计算出连续传递函数的小数部分迟延.这种参数估计的方法要求原连续 系统不带迟延部分的传递函数在离散化后不会产生零点,所以原传递函数应该具有标准形式:石(s)=e-。/DR(D,对于二阶级惯性加纯迟延模型,t=dh+玩h是采样周期,d是非负整数这样离散化后的脉冲传递函数为:酢阳㈨=号瓮掣G(s)=石丽IS e~丽Ke-d:从而得到最小二乘形式为:yO)=一aly(t―1)一a2y(t一2)+blu(t―d―1)+b2u(t―d一2)+b3u(t―d一3) 根据上述推导,对文献【30】中提到的算法进行了修改,我们得到带有小数部分 纯迟延时间的传递函数标准形的参数佶计算法步骤如下;1)初始化被辨识参数d,五。,…,a:,a.,d∈【d曲,‰】并且使占的估计值苦等于零,也就是初始化毛=丘:0,龟≠0;2)’利用输入输出数据,通过递推最小二乘法以及最小化损失函数对未知参数 %,..?,a2,q,毛,如,6I以及d进行修正,采用上述的整数部分纯迟延估计算法;3)计算系统的零点z'着llz。,并且根据公式计算善=瞳一茜)/瞧+龟+毛);4)如果一1<害<l,则转到7),否则转到5);5)反复进行以下算法直到渊<1:如果善<一1,则令(0一‘孑一l且吾.--h善+1),如果善>l,则令(d--hd+坦荸---h吾一1); 华北电力大学硕士学位论文6)令,=gj+龟+蠢,害≥o时,设丘=o,意=(1一舍)扣毛=耖,占so时,设毛=o,b2=(1+锄6l=匆那么于=h(a4-句;7)将离散模型参数转化成连续模型参数,转2)。 5.4.2仿真实验分析 选取以下两个模型作为研究对象,利用上述的辨识算法,来估计模型中的参数。G№,=毒‰,晓伊羔采用MATLAB编程来实现以上算法,对模型的参数进行辨识,具体的辨识过 程可以简单的描述为;首先进行数据的采集,输入为伪随机信号,输出为上述模型 的输出,取采样周期为一秒;接下来利用参考模型输入数据和输出数据,采用上述 辨识算法1).7)步,得到离散模型的参数,将离散模型的参数转化为连续模型的 参数;再进行下一次数据的采集,如此一直到参数收敛,结束运算。通过辨识得到离散传递函数和连续传递函数中各个参数的收敛情况。t/st,s图5-2离散模型一参数收敛曲线图5-3连续模型一参数收敛曲线图5-4离散模型二参数收敛曲线圈5-5连续模型一参数收敛曲线 华北电力大学硕士学位论文表5-1模型参数的收敛终值图5.2,图5.3分别为模型一的参数的收敛曲线,图5-4,图5-5分别为模型二 的参数的收敛曲线图中的al,a2’b1,b2,b3,d是辨识得到的(5-15)式中离散模型 的参数,e是小数部分纯迟延时阅;a,b,de,K是计算得到的(5-16)式中连续模型的 参数。由图中参数收敛曲线,我们可以看出在递推算法中,纯迟延时间的估计值不 会随着其他参数估计的变化而变化,也就是说在其他参数逐渐得到修正的每一步递 推中,纯迟延没有一直变化。随着辨识数据的采集和递推算法对参数的修正。纯迟 延时间首先达到某一定值不再发生变化,随后参数估计值的变化越来越小并最终收 敛到真值。经过60次采样计算便可以得到各个参数的收敛值如表5-1。得到两个模型的传递函数为:饥,=意豢‰,会产生误差。觎加意‰模型参数收敛值和真值间的误差较小,其中小数部分纯迟延时间误差相对大一 些,主要是小数部分纯迟延的计算公式会带来一定的误差,另外模型转换过程中也5.5本章小结本章介绍并推导了基于递推最/bZ-乘在线辨识算法,对带有纯迟延时间系统进 行参数估计。连续系统小于采样周期的纯迟延时间,在相应的采样系统中就会产生 一个副实零点,利用这个零点引起的相位变化来对小于采样周期的迟延进行估计; 采用最小化损失函数的两步辨识法,对系统的其他参数和采样周期整数倍的迟延进 行估计。利用这种递推辨识算法对热工模型进行参数辨识,得到模型的离散传递函 数进而求得其对应的连续传递函数。结果表明这种方法对典型二阶惯性加纯迟延对 象进行参数估计时精度比较高,在控制系统性能评价中我们就采用这种方法辨识得 到模型的纯迟延时间。在参数辨识中,输入信号必须充分激励,本文通过观察的方 法来选取合适的数据;另外,通过估计得到的相邻若干参数的绝对值误差和来判断 参数的收敛情况,给出何时停止递推. 华北电力大学硕士学位论文第六章基于设定值响应的蒸汽温度/压力系统控制性能评价针对串级喷水减温控制的过热汽温系统和燃料量控制的主蒸汽压力系统,以内回路和惰性区组成的广义被控制对象来作为本章控制器性能评价的研究对象,随着 运行工况发生变化,过热汽温和主蒸汽压力的动态特性有很大的变化,对于整定好 的主调节器和副调节器的控制性能也会随着对象特性的变化而变差。本章分别仿真 研究了四种典型工况下的过热汽温和主蒸汽压力控制系统,利用上述几章的性能评 价方法对不同工况下主调PID控制器的性能进行了评价,给出了量化的性能等级。6.1串级过热汽温控制系统锅炉蒸汽温度控制直接影响到全厂的热效率和设备的安全运行,因此汽温控制 系统是锅炉的重要控系统之一。目前广泛采用喷水来控制蒸汽温度,对于这种控制 方式,喷水量扰动就是基本扰动。由于过热器是具有分布参数的对象,可以把管内 的蒸汽和金属管壁看作是无穷多个单容对象串联组成的多容对象。当喷水量发生变 化后,需要通过这些串联单容对象,最终引起出口蒸汽温度的变化,因此出口汽温 的响应具有很大的迟延。减温器离过热器出口温度越远,迟延越大。 串级汽温调节系统具有两个闭环控制回路,可以简单的描述为:由被调对象的 导前区和副调节器组成内回路(或称副回路);由被调对象的惰性区、主调节器和内回路组成外回路(或称主回路)。串级汽温调节系统基本框图,如图6.1所示。 其中%。(s)是导前区对象传递函数;Wo:(s)是惰性区对象传递函数;日(J)是副调传 递函数,D,(j)是主调传递函数。一麓逐;37图6.1串级控制系统框图导前汽温能够快速反映扰动,尤其是减温水的自发性扰动,副调节器根据导前 汽温的变化改变减温水量,对主汽温起粗调作用。副回路任务是快速消除内扰,要 求调节过程的持续时间较短,但不要求无差;主调节器则通过对副回路的校正,对 主汽温起细调作用,主调的任务是维持主汽温恒定,当主汽温偏离给定值时,主调 输出校正信号,使副调不断调节减温水量,直到主汽温恢复到给定值为止。稳定时, 导前汽温可能稳定在与原来不同的数值上,但主汽温则一定等于给定值。弼一 华北电力大学硕士学位论文广义被控制对象的单回路PID反馈控制系统G(D等价于虚线框内内回路与惰性区对象的串联,G(D可以利用运行数据“(r)和y(f)通过最小二乘方法辨识得到。串级汽温调节系统框图可以进一步等效简化为如下形式,串级单位反馈控制回路。图6.2等效的单回路控制框图表6.1为某超临界600MW直流锅炉高温过热器在减温水流量扰动下二级减温 控制对象动态特性的数学模型,其动态特性如表所示,通过4种典型负荷工况点的 试验,分别给出了高温过热器导前区和惰性区的传递函数‘341。 表6-1某超临界600MW直流锅炉高温过热器的动态特性利用上述对象作为实际的过热汽温控制对象来仿真产生性能评价所需要的闭 环响应数据。这里副回路pID控制器参数:DI(s)=k。+1/Tjs,其中k。=8’互=25。这时候在四中工况下副回路都有比较好的设定值响应特性。将内回路和惰性区对象 串连作为广义的被控对象,我们针对广义的过热汽温对象,对主调节器的跟踪性能进行评价。当设定值发生阶跃变化时,利用单位阶跃响应的输出数据计算得到性能 指标,过渡过程时间,绝对值误差积分,超调量;同时利用最小二乘法辨识得到一 阶惯性加纯迟延形式的广义被控对象模型或者二阶形式{利用纯迟延时间对性能指 标进行无量纲化,得到无量纲的性能评价指标。6.2广义被控对象的辨识主调节器采用el控制器,形式为c(s)=k,+11r,s 调节器的参数。 37%时,kpj=0.36,瓦=588;50%el-,k,2=0.35,正2=400; 在不同的工况下整定得到主75%时,‘3=o.38,死=286;100%时,k,.;o.33,死=182?38 华北电力大学硕士学位论文利用第五章中的热工模型纯迟延时间的递推辨识算法,在输入阶跃信号的基础 上,我们得到四种不同工况下广义被控对象的二阶惯性加纯迟延模型。.(啬裴荷k#s)吒(D=赢矿1.0丽42试验工况表6-2过熟汽温广义被控对象模型 广义被控对象的二阶模型eⅢ知.(d=5。2%42负荷.2kg/s)吃(力=而丽1.1而18丽8跏.(d写臻毒。)吃(曲=丽瓦1.面2丽e棚“.(0燃.Sk荷g/。)丸(s)=丽万1.2而73鬲e~,对比广义被控对象的二阶惯性加纯迟延模型的动态特性,可以看出模型能很好 地复现原被控对象的动态特性,一方面说明实际的热工对象能很好地由惯性加纯迟 延模型来近似,另一方面表明第五章所研究的热工模型的辨识算法精确度比较高, 能够用来对实际的热工对象进行闭环辨识。图6-3二阶模型(~)和实际对象(.一)的阶跃响应曲线对比 由动态曲线可以看出,二阶惯性加纯迟延形式模型能很好地描述实际对象的动 态特性。在辨识对象模型的同时,利用输出数据计算得到这种最优PID参数控制器 参数下的四种不同工况的性能指标,利用纯迟延时间对性能指标无量纲化得到无量 纲的性能评价指标,得到控制器性能的优劣。 华北电力大学硕士学位论文表6.3不同工况过热汽温系统最优性能指标6.3主调节器跟踪性能评价如果运行工况发生变化,对象的特性有很大改变,在控制器参数不变的情况下, 控制系统的控制性能有变差的趋势,这就需要采集阶跃响应数据,计算性能评价指 标对控制器的性能进行评价。下面分别采用四种工况下整定好的控制器(37%(a), 50%(b),75(e),t00(d))在运行工况发生变化时,对主调节器的控制性能进行评价。,表6-4过热汽温系统主调节器控制性能评价乃尉Ed2.7919 4,0285 5.6864 8.6354OS(a) 50% 75%评价(b) 37% 50% 75% 100% (d) 37% S0% 75%乃11.6289 6.5867 10.0926 18.5458朋E玎0s3.5727 3.0297 3.93 14 5.9540 0.2807 0.0501 0 00s评价甲 昆 艮37%4.4337 10.0875 17.60740.0463 00优良 中 司loo%29.1235 (c) 37% 50%00s甲 评价j莹 司 民乃14.1398 11.8066埘Ed5.4810 3.6368 2.98024.3305评价中缸∞ 23.4852 10.0405朋Edoo0.5196 0.2730 0.0266 0一5.0048 0.6465 0.3336中7.4400 3.8049 2.979875%4.8589 100% 10.1530优艮L00%4.92030.0217优从表中可以看出,随负荷(或主蒸汽流量)的变化,模型参数变化很大,特别是惰性区的时间常数T,以及导前区的静态增益K。用一套固定的控制器参数不能保证各种工况下控制器的品质,因此,在控制系统设计时,应考虑变参数控制方案。 下面仿真得到四种PI控制器参数下,不同工况时的单位阶跃响应,将上述评 价结果与响应曲线进行对比,可以直观地验证评价方法的合理性和结果的正确性。 表中(a,b,c,d)与下图(a,b,c,d)相对应,表示四种工况。S(a)kp--0.36,TI=588四种工况的阶跃响应应40 华北电力大学硕士学位论文(c)七p=O.38,正=286四种工况的阶跃响应(d)七p=0.33,瓦=182四种工况的阶跃响应图中曲线(.一)表示37%负荷下阶跃响应,r)表示50%负荷下阶跃响应,(…-)表示75%负荷下阶跃响应,(一.一.一.)表示100%负荷下阶跃响应。图6-4不同工况下过热汽温控制系统阶跃响应曲线6.4主汽压力控制系统跟踪性能评价6.4.1主汽压力控制系统 火电厂的主汽压力一直通过改变燃料量和送风量来进行调节,主汽压力相对于 燃料量变化是一个大惯性、大滞后的对象。燃料在炉内的燃烧过程是非常迅速的, 一般燃料在炉内停留约为几秒钟,相对于大滞后的主汽压力来讲是一个快速的过 程,因此由火焰图象处理技术得到的炉膛辐射能信号能够快速反映进入炉膛的燃料 量变化。文献[351采用炉膛辐射能信号作为中间被调量,构成了串级PID控制系统。同样,针对串级PID主汽压控制系统,如图6.1所示,图中D’(曲是主调节器;D1(s)是副调节器;%(砖是燃料量与妒膛辐射能信号之间的传递函数;Wo。④是炉膛辐射能与主汽压力之间的传递函数;F(D为燃料量扰动信号;艺(s)为主汽压力 给定值;P(J)为主汽压力测量值;“s)为虚线框内的广义被控对象;识s)为以s)的内部模型。 袭6-5某超临界600MW直流锅炉主汽压力的动态特性41 华北电力大学硕士学位论文6.4.2主汽压力模型辨识 副回路调节器参数取.j}。=8,霉=20副回路具有比较好的给定值响应特性,我们 将广义的被控对象模型作为研究的对象,对主调节器进行性能评价。副调对象和PI控制器构成单回路系统,按随动系统整定调节器参数,将整定好的单回路系统等效为~个环节,与主对象共同形成广义的被控对象以J),在此基础上建立广义被控对 象坝s)的内部模型,然后设计外部控制器并整定其参数。表6-6主汽压力广义被控对象模型这笙三况44%负荷广义被控对象的三丽硬亘哦。)=虿丽语≮1.i66夏i;石e“jJ3||、,Stl厶己.6s+I62%负荷88%负荷哌∽=历矛蕊1.83瓦/鬲。“‘二)‘JtoQ.‘5十l碗(s)=面丽72i.0面88丽8Ⅲj,1竺竺荷吭o)=历瓦百了2,4五56了三鬲em“6.4.3主调节器跟踪性毹评价对于不同工况整定得到控制器参数,仿真得到四种PID控制器参数下,不同工况时的单位阶跃响应。44%K,=0.2501,K』=O.0051,Kd=1.9702 62%Kp=0.2922,Kj=0.0047,Ka=2.9700外回路PID控制器参数按照IMc理论设计,采用Do)=Kp+等+匕J DI(s):0.2501+旦塑+1.9702s D2(s):0,2922+旦:业+2.9700sS J88%Kp=0.3330,Ks=0.0044,Kd=4.1461 100%K,=0.3677,Kl=0.0036,Ka=6.0939D3(s):0.3330+.0.0044+4.146JJD4(s):0。3677+竺堂+6.0939sS计算得到在每种工况下,控制系统的动态性指标,并利用MATLAB仿真得到 阶跃响应曲线。 华北电力大学硕士学位论文图6-5最优PID下的阶跃响应曲线运行工况发生变化时,则被控对象的特性将发生很大的变化,如果控制器参数 不变,控制器的性能有变差的趋势。下面分别采用44%(a)和100%(d)工况下整定好 的控制器,在运行工况发生变化时,对主调节器的控制性能进行评价。表6.3主蒸汽压力系统主调节器跟踪性能评价在工况发生改变时,被控对象的特性发生变化,由上述性能评价可以看出,?在 某一工况下整定的PID控制器参数,基本上能够适应工况的变化,在不同工况下仍 有比较好的设定值跟踪性能。表中(a,d)与下图中(a,d)相对应,表示两种工况,一表示37%负荷下阶跃响应,――表示62%负荷下阶跃响应,….表示88%负荷下阶跃响应,....表示100%负荷下阶跃响应。(a)Dl(s)控制下四种工况的阶跃响应(d)D4(s)控制下四种工况的阶跃响应圈6-6不同工况下主汽压力控制系统阶跃响应益线 华北电力大学硕士学位论文6.5本章小结以串级过热汽温系统和串级主蒸汽压力控制系统的广义被控对象作为性能评 价的研究对象,利用两个系统的典型模型来仿真产生性能评价所需要的阶跃响应数 据,利用仿真得到的数据一方面辨识得N--阶惯性加纯迟延模型,另一方面计算出 稳定时间,绝对值误差积分和超调量,进而得到无量纲的性能评价指标,从评价的 结果以及对应仿真产生的阶跃响应曲线都可以看出本文所采用的性能评价方法与 直观的曲线特性相符合。对于实际的控制系统,当有阶跃信号产生时,获取闭环响 应数据,进行计算得到性能指标,并辨识得到被控对象模型的纯迟延时间,利用无量纲性能指标对应跟踪性能评价等级,得

我要回帖

更多关于 自动控制系统性能指标 的文章

 

随机推荐