CN108170637B - 一种具有微分特性过程的传递函数模型辨识方法 - Google Patents
一种具有微分特性过程的传递函数模型辨识方法 Download PDFInfo
- Publication number
- CN108170637B CN108170637B CN201711274324.4A CN201711274324A CN108170637B CN 108170637 B CN108170637 B CN 108170637B CN 201711274324 A CN201711274324 A CN 201711274324A CN 108170637 B CN108170637 B CN 108170637B
- Authority
- CN
- China
- Prior art keywords
- time
- transfer function
- data
- step response
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Feedback Control In General (AREA)
- Complex Calculations (AREA)
Abstract
Description
技术领域
本发明属于自动控制技术领域,具体涉及一种具有微分特性过程的传递函数模型辨识方法。
背景技术
发明内容
发明目的:针对上述具有微分特性过程,提出一种通过计算过程阶跃响应曲线上的几个特征参数,获得过程相应传递函数模型的方法。
技术方案:一种具有微分特性过程的传递函数模型辨识方法,所述的模型辨识方法包括如下步骤:
(1)进行现场动态特性试验,获取过程阶跃响应曲线数据;
(2)根据过程阶跃响应曲线数据,求取阶跃响应曲线的特征参数;
(3)根据过程阶跃响应曲线的特征参数,求取过程的传递函数。
进一步的,步骤(1)具体包括在过程输入u和过程输出y保持稳定不变的条件下,记录y的稳态值y0,并将u阶跃变化Δu,并以T秒为采样周期,采集过程输出y在各采样时刻的数据y'(k),用y'(k)减去y0,得到过程的阶跃响应数据序列y(k),其中k为采样时刻,k=1,2,…,N,N为采样数据个数,满足N*T大于过程输出从变化到稳定所需的时间,T取1-5秒;
步骤(2)中过程阶跃响应曲线的特征参数为曲线极值点对应的时间T0和相应的极值Km、T0时刻之前曲线的拐点对应的时间Tq,特征参数的具体计算步骤如下:
(2.1)以y的数据序列{y(i)}及y(i)对应的采样时刻时间t数据序列{i·T}为样本,其中i=k0-M,k0-M+1,...,k0,k0+1,...,k0+M,k0为y(k)中绝对值最大的数据所对应的采样时刻,M<min(k0,N-k0),min为取小运算,M取2-5,通过最小二乘法拟合,得到二次多项式y=at2+bt+c;
(2.3)以y的数据序列{y(j)}及y(j)对应的采样时刻时间t数据序列{j·T}为样本,其中j=N1,N1+1,...,N2,N1=int(0.1k0),N2=int(0.9k0),int为取整运算,通过最小二乘法拟合,得到三次多项式y=At3+Bt2+Ct+D;
步骤(3)求取过程的传递函数方法为:
将步骤2中计算得到的过程阶跃响应曲线特征参数T0、Km和Tq代入到如下计算式:
计算得到n、Tc和K,则具有微分特性过程的传递函数模型为其中n为模型阶次,Tc为时间常数,K为模型传递系数;如果计算得到的n是小数,n=n0+α,其中n0为n的整数部分,α为n的小数部分,则过程的传递为
更进一步的,所述的模型辨识方法适用于模型阶次n>1的过程。
本发明根据传递函数通过反拉氏变换,得到其阶跃响应曲线y(t),通过求该曲线的极值和拐点,推导出传递函数三个参数n、Tc和K与其阶跃曲线的三个特征参数Tq、T0和Km的关系式,其中Km是曲线的极值,T0为极值点对应的时间,Tq为T0时刻之前曲线的拐点对应的时间。这样,通过现场试验获取过程阶跃响应曲线,计算阶跃响应曲线的三个特征参数Tq、T0和Km,就可以得到过程的传递函数模型。
有益效果:本发明相比现有技术的显著效果在于采用该方法求解具有微分特性过程的传递函数,可直接根据过程的阶跃响应曲线,获得过程的传递函数模型,其计算过程简单,便于工程应用。
附图说明
图1为过程阶跃响应曲线及其特征参数。
具体实施方式
为了详细的说明本发明公开的技术方案,下面结合说明书附图和具体实施例做进一步的阐述。
步骤1:进行现场动态特性试验,在稳态工况下,将过程输入u阶跃变化Δu,并以T秒为采样周期,即过程输入u和过程输出y保持稳定不变的条件下,记录y的稳态值y0,并将u阶跃变化Δu,采集过程输出y在各采样时刻的数据y'(k),用y'(k)减去过程输入u阶跃变化前y的稳态值,得到过程的阶跃响应数据序列y(k),其中k为采样时刻,k=1,2,…,N,N为采样数据个数,满足N*T大于过程输出从变化到稳定所需的时间;
在此处过程输入阶跃幅度取Δu=1,采样数据个数取N=60,阶跃响应数据y(k)如下:
步骤2:找出y(k)中绝对值最大的数据所对应的采样时刻k0,以y的数据序列{y(i)}及y(i)对应的采样时刻时间t数据序列{i·T}为样本,其中:i=k0-M,k0-M+1,...,k0,k0+1,...,k0+M,M<min(k0,N-k0),一般M取2-5,通过最小二乘法拟合,得到二次多项式y=at2+bt+c;
通过上表得到k0为11,取M为3,可得到拟合的二次多项式为y=-0.00002161x2+0.0022x-0.0023;
步骤4:以y的数据序列{y(j)}及y(j)对应的采样时刻时间t数据序列{j·T}为样本,其中j=N1,N1+1,...,N2,N1=int(0.1k0),N2=int(0.9k0),通过最小二乘法拟合,得到三次多项式y=At3+Bt2+Ct+D,易查得N1=3,N2=8,拟合得到多项式为y=-0.0000006953x3+0.00002549x2+0.0015-0.0064;
步骤6:将步骤3和步骤5中得到的T0、Km和Tq代入到如下计算式:
Claims (2)
1.一种具有微分特性过程的传递函数模型辨识方法,其特征在于,包括如下步骤:
(1)对以汽机调门开度为输入,发电功率为输出的过程进行现场动态特性试验,获取过程阶跃响应曲线数据,具体包括:在过程输入汽机调门开度u和过程输出发电功率y保持稳定不变的条件下,记录y的稳态值y0,并将u阶跃变化Δu,并以T秒为采样周期,采集过程输出y在各采样时刻的数据y'(k),用y'(k)减去y0,得到过程的阶跃响应数据序列y(k),其中k为采样时刻,k=1,2,…,N,N为采样数据个数,满足N*T大于过程输出从变化到稳定所需的时间,T取1-5秒;
(2)根据过程阶跃响应曲线数据,求取阶跃响应曲线的特征参数,过程阶跃响应曲线的特征参数为曲线极值点对应的时间T0和相应的极值Km、T0时刻之前曲线的拐点对应的时间Tq,特征参数的具体计算步骤如下:
(2.1)以y的数据序列{y(i)}及y(i)对应的采样时刻时间t数据序列{i·T}为样本,其中i=k0-M,k0-M+1,...,k0,k0+1,...,k0+M,k0为y(k)中绝对值最大的数据所对应的采样时刻,M<min(k0,N-k0),min为取小运算,M取2-5,通过最小二乘法拟合,得到二次多项式y=at2+bt+c,a、b和c为二次多项式的系数;
(2.3)以y的数据序列{y(j)}及y(j)对应的采样时刻时间t数据序列{j·T}为样本,其中j=N1,N1+1,...,N2,N1=int(0.1k0),N2=int(0.9k0),int为取整运算,通过最小二乘法拟合,得到三次多项式y=At3+Bt2+Ct+D,A、B、C和D为三次多项式的系数;
(3)根据过程阶跃响应曲线的特征参数,求取过程的传递函数,求取方法为:
将步骤(2)中计算得到的过程阶跃响应曲线特征参数T0、Km和Tq代入到如下计算式:
2.根据权利要求1所述的一种具有微分特性过程的传递函数模型辨识方法,其特征在于:所述的模型辨识方法适用于模型阶次n>1的过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711274324.4A CN108170637B (zh) | 2017-12-06 | 2017-12-06 | 一种具有微分特性过程的传递函数模型辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711274324.4A CN108170637B (zh) | 2017-12-06 | 2017-12-06 | 一种具有微分特性过程的传递函数模型辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108170637A CN108170637A (zh) | 2018-06-15 |
CN108170637B true CN108170637B (zh) | 2021-11-12 |
Family
ID=62525225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711274324.4A Active CN108170637B (zh) | 2017-12-06 | 2017-12-06 | 一种具有微分特性过程的传递函数模型辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108170637B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110887943B (zh) * | 2019-12-27 | 2022-07-08 | 中国科学院合肥物质科学研究院 | 一种氧气浓度传感器响应速度提升方法 |
CN112180801B (zh) * | 2020-10-13 | 2021-07-09 | 郑州大学 | 一种过程控制系统的控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103744286A (zh) * | 2013-12-31 | 2014-04-23 | 广东电网公司电力科学研究院 | 一种火力发电系统的控制器的设计方法和装置 |
CN105204336A (zh) * | 2015-11-02 | 2015-12-30 | 江西洪都航空工业集团有限责任公司 | 一种飞机运动模态辨识方法 |
CN105929683A (zh) * | 2016-06-23 | 2016-09-07 | 东南大学 | 一种微分可调pid控制器参数工程整定方法 |
CN106054596A (zh) * | 2016-06-23 | 2016-10-26 | 东南大学 | 一种基于改进性能指标的pid控制器参数优化整定方法 |
CN106054610A (zh) * | 2016-06-23 | 2016-10-26 | 东南大学 | 一种优化的pi控制器参数工程整定方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6697767B2 (en) * | 2000-10-18 | 2004-02-24 | The National University Of Singapore | Robust process identification and auto-tuning control |
JP6017232B2 (ja) * | 2012-09-07 | 2016-10-26 | ルネサスエレクトロニクス株式会社 | 送信装置及び通信システム |
-
2017
- 2017-12-06 CN CN201711274324.4A patent/CN108170637B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103744286A (zh) * | 2013-12-31 | 2014-04-23 | 广东电网公司电力科学研究院 | 一种火力发电系统的控制器的设计方法和装置 |
CN105204336A (zh) * | 2015-11-02 | 2015-12-30 | 江西洪都航空工业集团有限责任公司 | 一种飞机运动模态辨识方法 |
CN105929683A (zh) * | 2016-06-23 | 2016-09-07 | 东南大学 | 一种微分可调pid控制器参数工程整定方法 |
CN106054596A (zh) * | 2016-06-23 | 2016-10-26 | 东南大学 | 一种基于改进性能指标的pid控制器参数优化整定方法 |
CN106054610A (zh) * | 2016-06-23 | 2016-10-26 | 东南大学 | 一种优化的pi控制器参数工程整定方法 |
Non-Patent Citations (1)
Title |
---|
基于阶跃辨识的PID自整定算法研究及其应用;邱亮;《中国优秀硕士学位论文全文数据库 信息科技辑》;20130715(第2013年7期);I140-399 * |
Also Published As
Publication number | Publication date |
---|---|
CN108170637A (zh) | 2018-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110609973B (zh) | 一种用于流量测量的卡尔曼滤波方法 | |
CN108170637B (zh) | 一种具有微分特性过程的传递函数模型辨识方法 | |
JP2015076024A (ja) | プラント制御装置の制御ゲイン最適化システム | |
CN104199993B (zh) | 一种针对间歇过程的二维动态卡尔曼滤波器设计方法 | |
CN117236084B (zh) | 一种木工机械加工生产智能管理方法及系统 | |
CN113960923B (zh) | 基于离散扩展状态观测器的无模型自适应滑模控制方法 | |
CN101997471B (zh) | 基于pid预测函数的励磁控制方法 | |
CN110705042A (zh) | 基于t-s模糊模型的汽轮机阀门流量特性辨识方法、系统及介质 | |
CN103684351A (zh) | 具有pid指数因子的自适应惯性滤波方法 | |
Cvejn | The design of PID controller for non-oscillating time-delayed plants with guaranteed stability margin based on the modulus optimum criterion | |
CN111854109A (zh) | 一种室内温度控制方法、装置、电子设备及存储介质 | |
CN107102553B (zh) | 基于时滞鲁棒模型算法的并网逆变器控制方法 | |
CN108549211A (zh) | 一种时滞系统的分数阶pid控制器设计方法 | |
CN116881665A (zh) | 一种基于CMOA优化的TimesNet-BiLSTM光伏出力预测方法 | |
CN105865621A (zh) | 激光功率的校准方法和装置 | |
CN106569406B (zh) | 基于尺度变换的有刷直流电机模型参数快速辨识方法 | |
Du et al. | Predicting the solar maximum with the rising rate | |
JP2020064439A (ja) | 制御モデル同定方法、制御モデル同定装置及びプログラム | |
CN103701738B (zh) | 一种预失真参数确定方法、设备及数字预失真处理系统 | |
CN110262221B (zh) | 一种热工过程中对象的pid控制器参数控制方法 | |
Alkhafaji et al. | A novel PID robotic for speed controller using optimization based tune technique | |
Sariyildiz et al. | Robust stability and performance analysis of the control systems with higher order disturbance observer: Frequency approach | |
CN108932483B (zh) | 一种基于机器嗅觉的气味识别方法及装置 | |
CN112256067B (zh) | 温度控制方法、装置、电子设备及存储介质 | |
CN114711638A (zh) | 烘焙设备的控制方法、装置、烘焙设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |