CN111832115B - 对称迟后-超前校正的一体化频率设计方法 - Google Patents
对称迟后-超前校正的一体化频率设计方法 Download PDFInfo
- Publication number
- CN111832115B CN111832115B CN202010489269.6A CN202010489269A CN111832115B CN 111832115 B CN111832115 B CN 111832115B CN 202010489269 A CN202010489269 A CN 202010489269A CN 111832115 B CN111832115 B CN 111832115B
- Authority
- CN
- China
- Prior art keywords
- correction
- phase angle
- determining
- frequency
- link
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 77
- 238000013461 design Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000010008 shearing Methods 0.000 claims abstract description 12
- 238000012546 transfer Methods 0.000 claims description 13
- 230000005251 gamma ray Effects 0.000 claims description 2
- 230000010363 phase shift Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了对称迟后‑超前校正的一体化频率设计方法,该方法包括(S1)根据所要求的稳态性能指标确定系统的开环增益;(S2)利用已知的开环增益,绘制未校正系统G0(s)的Bode图,并计算未校正系统的剪切频率ωc0,相角裕度γ0和增益裕度Lg0;(S3)根据相角裕度的要求确定超前‑迟后校正环节的α,为使相角裕度r0达到要求值,计算超前环节所需提供的超前相角(S4)确定校正后系统的剪切频率ωc;(S5)确定超前校正环节;(S6)确定迟后校正环节;(S7)检验是否满足系统的性能指标,若不满足要求,可增大附加相角Δ的值,从步骤(S3)重新计算。通过上述方案,本发明达到了一次性处理迟后与超前的目的,具有很高的实用价值和推广价值。
Description
技术领域
本发明属于飞行器的液压伺服控制系统频域技术领域,具体地讲,是涉及 对称迟后-超前校正的一体化频率设计方法。
背景技术
本发明考虑无线电控制的飞行器上伺服液压伺服系统的设计。液压伺服系 统在无线电控制的飞机上被用来控制尾舵和船舵。由液压伺服系统构成的舵机 是保证飞行器性能的关键执行器件,用于测试舵机的可靠性,对于提高飞机的 安全性具有非常重要的意义。本发明提供了一种对称迟后-超前校正的一体化频 率设计方法对该液压伺服系统性能进行改进,以增强系统的稳定性,提高系统 的响应速度。
频率校正法是根据频域指标确定校正装置的参数。迟后-超前校正环节的传 递函数为:
其中,和/>分别超前校正环节和迟后校正环节的传递函数。迟 后-超前校正的传统方法是将迟后环节和超前环节分开进行设计,没有充分利用 迟后-超前校正整体环节的特性。采用迟后-超前分离校正方法往往需要反复校正 多次才能达到期望的性能。从而传统的迟后-超前校正方法存在耗时长、效率低 等问题,因此如何改变传统的迟后-超前校正方法的缺陷是本领域技术人员亟需 解决的问题。
本发明针对具有如下形式的迟后-超前校正环节:
给出一体化的设计方法,其中,α,τ,T是需要设计的参数。本发明利用该迟 后-超前校正环节的整体频率特性进行设计,该方法的有益效果是:该方法充分 利用在校正后系统剪切频率附近,校正环节具有正的超前相角和负的幅频特性 这一特点,一体化地对校正参数进行设计。具体说来,就是利用在校正后剪切 频率附件,通过校正环节提供超前相角;通过该迟后-超前校正环节在剪切频率 处具有的负频率特性10lgα将未校正系统的正幅值拉回0分贝线。
发明内容
本发明的目的在于提供对称迟后-超前校正的一体化频率设计方法,主要解 决现有技术中存在的迟后与超前分开设计存在耗时长、效率低等问题。
为了实现上述目的,本发明采用的技术方案如下:
对称迟后-超前校正的一体化频率设计方法,基于如下形式的对称迟后-超前 校正环节进行校正:
其中,α,τ,T是待设计参数;Gc(s)在ω=ωm附近,且时,其相角近 似为其幅值为-10lgα,其中/>Δ1∈(6°~12°);
具体校正步骤如下:
(S1)根据所要求的稳态性能指标确定系统的开环增益;
(S2)利用已知的开环增益,绘制未校正系统G0(s)的Bode图,并计算未校 正系统的剪切频率ωc0,相角裕度γ0和增益裕度Lg0;
(S3)根据相角裕度的要求确定超前-迟后校正环节的参数α,为使相角裕 度r0达到要求值,计算超前环节所需提供的超前相角即:
其中,γ为要求的相角裕度;γ0c为期望未校正系统提供的相位储备,满足 γ0c<180°+∠G0(jωcL),式中ωcL为所要求频率的下界;Δ为附加相角;
其中,迟后-超前校正环节的参数α由下式确定:
(S4)确定校正后系统的剪切频率ωc;
(S5)确定超前校正环节;
(S6)确定迟后校正环节;
(S7)检验是否满足系统的性能指标,若不满足要求,可增大附加相角Δ的 值,从步骤(S3)重新计算。
进一步地,所述步骤(S4)中在校正后的剪切频率处ωc,未校正系统的正 幅值由迟后-超前校正环节拉回到0dB线,因此校正后系统的剪切频率ωc由下式 确定:
20lg|G0(ωc)|=10lgα (3)
如果求得的剪切频率ωc满足要求,则进行下一步;否则返回步骤(S3)中 调整Δ和γ0c。
进一步地,所述步骤(S5)中确定超前校正环节是让超前校正环节取得最 大值的频率对准校正后系统的剪切频率,即ωm=ωc,根据确定校正装 置的转时间常数,进而确定超前环节的传递函数。
具体地,所述步骤(S6)中确定迟后校正环节是为了减小串联迟后校正对 系统相位裕度的影响,要求校正环节在剪切频率ωc处的迟后相移在6°~12°以下, 应选择:
进而确定迟后-超前校正装置的传递函数。
其中,在上面的设计方法中,附加相角Δ是考虑到迟后环节的相位延迟Δd和校正装置的近似最优频率与剪切频率ωc不一致而产生的附件相角Δ0两方面的因素。 如果迟后环节的转折频率在/>之间,相位延迟Δd在6°~12°。另外,由于 校正后的剪切频率ωc比未校正系统的剪切频率ωc0小,一般来说,未校正系统在 剪切频率ωc附近不会产生大的相位衰减,附加相角Δ0不用太大,甚至是0都可以, 综合两方面的因素,附加相角Δ一般取Δ=10°~15°即可。如果选择/>可选 取Δ=10°。
与现有技术相比,本发明具有以下有益效果:
本发明充分利用了迟后-超前校正环节作为一个整体的频率特性,结合迟后 环节与超前环节的频率特性进行一体化设计,有效解决了传统技术中存在的迟 后与超前分开设计存在耗时长、效率低等问题。本发明由于将迟后环节与超前 环节进行了一体化设计,因此在校正的过程中通常一次就能达到性能指标,具 有耗时短、处理效率高等优点。
附图说明
图1为本发明的系统流程图。
图2为本发明对称迟后-超前校正一体化设计方法示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明,本发明的实施方式包括但 不限于下列实施例。
实施例
如图1与图2所示,对称迟后-超前校正的一体化频率设计方法,
对具有如下形式的对称迟后-超前校正环节进行校正:
其中,α,τ,T是待设计参数,Gc(s)在ω=ωm附近,且时,其相角近 似为其幅值为-10lgα,其中/>Δ1∈(6°~12°)。
通过以下例子来说明本发明方法。
例:对于一个飞行器上的电动液压伺服系统,要设计一个校正补偿器使该 伺服系统满足下述控制指标,其中该系统的开环传递函数为:
要求满足下列性能指标:
(1)单位斜坡输入r(t)=t时的稳态误差(2)相角裕度γ≥35°;(3) 剪切频率ωc≥20rad/s。
试用频率法(相位裕度法)法设计对称迟后-超前校正装置的传递函数Gc(s)。
(S1)根据所要求的稳态误差有取K=126。未校正的开环传 递函数为:
(S2)在K=126的情形下,可以确定未校正系统的剪切频率为
相角裕度为
γ0=-14.9° (8)
未校正系统是最小相位的,而且相角裕度小于零,因此是不稳定的。
(S3)根据设计要求,校正后系统的剪切频率的下界是ωcL=20rad/s,在此 频率处相位储备为180°+∠G0(jωcL)=8.2°。取期望未校正系统提供的相位储备 γ0c=7°。另外,取裕量Δ=10°。
为使校正后系统的相角裕度达到要求,超前环节应提供的超前相角为:
因此迟后-超前校正环节的参数α为:
(S4)按下式确定校正后系统的剪切频率ωc
20lg|G0(ωc)|=10lgα (3)
采用渐近幅频特性简单求解得ωc=24.80。满足设计要求ωc>20°。
(S5)取ωm=ωc进行设计,超前校正环节的转折频率为:
则超前校正环节为:
(S6)确定迟后校正环节。取则τ=0.40。此时,迟后校正环节为:
由前两步可得所设计的迟后-超前校正装置的传递函数为:
(S7)校正后系统的开环传函为:
按照设计过程,估算剪切频率在至/>之间,按下式求取剪切频率ωc:
即:
简单求解该方程可得ωc=24.90rad/s。
校正后系统的相位裕度为:
γ=32.68° (18)。
由于γ=32.68°小于初始设定的性能指标γ≥35°,因此不满足设计要求,需要 返回步骤(S3),重新设计。
返回步骤(S3),取裕量Δ=12°。
为使校正后系统的相角裕度达到要求,超前环节应提供的超前相角为:
因此迟后-超前校正环节的参数α为:
(S4)按下式确定校正后系统的剪切频率ωc:
20lg|G0(ωc)|=10lgα (3)
采用渐近幅频特性简单求解得ωc=23.4。满足设计要求ωc>20°。
(S5)取ωm=ωc进行设计,超前校正环节的转折频率为:
则超前校正环节为:
(S6)确定迟后校正环节。取则τ=0.43。此时,迟后校正环节为:
由前两步可得所设计的迟后-超前校正装置的传递函数为:
(S7)校正后系统的开环传函为:
按照设计过程,估算剪切频率在至/>之间,按下式求取剪切频率ωc:
即:
简单求解该方程可得ωc=23.5rad/s。
校正后系统的相位裕度为:
γ=39.6° (27)。
由于γ=39.6°>35°,因此满足设计要求。
由于满足设计要求,因此本发明实现了充分利用了对称迟后-超前环节作为 一个整体的频率特性,实现了一体化迟后-超前环节的一体化设计。
上述实施例仅为本发明的优选实施例,并非对本发明保护范围的限制,但 凡采用本发明的设计原理,以及在此基础上进行非创造性劳动而做出的变化, 均应属于本发明的保护范围之内。
Claims (1)
1.对称迟后-超前校正的一体化频率设计方法,其特征在于,基于如下形式的对称迟后-超前校正环节进行校正:
其中,α,τ,T是待设计参数;Gc(s)在ω=ωm附近,且时,其相角近似为其幅值为-10lgα,其中/>Δ1∈(6°~12°);
具体校正步骤如下:
(S1)根据所要求的稳态性能指标确定系统的开环增益;
(S2)利用已知的开环增益,绘制未校正系统G0(s)的Bode图,并计算未校正系统的剪切频率ωc0,相角裕度γ0和增益裕度Lg0;
(S3)根据相角裕度的要求确定超前-迟后校正环节的参数α,为使相角裕度r0达到要求值,计算超前环节所需提供的超前相角即:
其中,γ为要求的相角裕度;γ0c为期望未校正系统提供的相位储备,满足γ0c<180°+∠G0(jωcL),式中ωcL为所要求频率的下界;Δ为附加相角;
其中,迟后-超前校正环节的参数α由下式确定:
(S4)确定校正后系统的剪切频率ωc;在校正后的剪切频率处ωc,未校正系统的正幅值由迟后-超前校正环节拉回到0dB线,因此校正后系统的剪切频率ωc由下式确定:
20lg|G0(jωc)|=10lgα (3)
如果求得的剪切频率ωc满足要求,则进行下一步;否则返回步骤(S3)中调整Δ和γ0c;
(S5)确定超前校正环节;确定超前校正环节是让超前校正环节取得最大值的频率对准校正后系统的剪切频率,即ωm=ωc,根据确定校正装置的转时间常数,进而确定超前环节的传递函数;
(S6)确定迟后校正环节;确定迟后校正环节是为了减小串联迟后校正对系统相位裕度的影响,要求校正环节在剪切频率ωc处的迟后相移在6°~12°以下,应选择:
进而确定迟后-超前校正装置的传递函数;
(S7)检验是否满足系统的性能指标,若不满足要求,可增大附加相角Δ的值,从步骤(S3)重新计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010489269.6A CN111832115B (zh) | 2020-06-02 | 2020-06-02 | 对称迟后-超前校正的一体化频率设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010489269.6A CN111832115B (zh) | 2020-06-02 | 2020-06-02 | 对称迟后-超前校正的一体化频率设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111832115A CN111832115A (zh) | 2020-10-27 |
CN111832115B true CN111832115B (zh) | 2023-12-19 |
Family
ID=72899069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010489269.6A Active CN111832115B (zh) | 2020-06-02 | 2020-06-02 | 对称迟后-超前校正的一体化频率设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111832115B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101867096A (zh) * | 2009-04-10 | 2010-10-20 | 霍尼韦尔国际公司 | 产生用于系统相位校正的参考信号的系统和方法 |
CN103076806A (zh) * | 2011-10-26 | 2013-05-01 | 北京航天长征飞行器研究所 | 三回路自动驾驶仪控制参数一体化解析整定方法 |
CN106909761A (zh) * | 2017-04-12 | 2017-06-30 | 北华航天工业学院 | 设计相位超前校正环节的频域性能曲线图法 |
CN106919051A (zh) * | 2017-04-12 | 2017-07-04 | 北华航天工业学院 | 设计相位滞后‑超前校正环节的性能曲线图法 |
-
2020
- 2020-06-02 CN CN202010489269.6A patent/CN111832115B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101867096A (zh) * | 2009-04-10 | 2010-10-20 | 霍尼韦尔国际公司 | 产生用于系统相位校正的参考信号的系统和方法 |
CN103076806A (zh) * | 2011-10-26 | 2013-05-01 | 北京航天长征飞行器研究所 | 三回路自动驾驶仪控制参数一体化解析整定方法 |
CN106909761A (zh) * | 2017-04-12 | 2017-06-30 | 北华航天工业学院 | 设计相位超前校正环节的频域性能曲线图法 |
CN106919051A (zh) * | 2017-04-12 | 2017-07-04 | 北华航天工业学院 | 设计相位滞后‑超前校正环节的性能曲线图法 |
Non-Patent Citations (1)
Title |
---|
对称迟后-超前校正的一体化设计方法;张颖等;《电气电子教学学报》;第43卷(第5期);第127-133、161页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111832115A (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103076806B (zh) | 三回路自动驾驶仪控制参数一体化解析整定方法 | |
US9821837B2 (en) | Electric power steering apparatus | |
CN104252133B (zh) | 一种无人机纵向控制律平滑切换方法 | |
EP0037159B1 (en) | Control system for aircraft vertical path guidance | |
CN109508027B (zh) | 基于鲁棒控制理论的火箭助推发射的控制方法 | |
CN100404372C (zh) | 一种旋翼飞行器偏航控制的优化方法和实施该方法的系统 | |
US8275496B2 (en) | Longitudinal and vertical gust feed forward compensation using lateral control surfaces | |
US9671788B2 (en) | Vertical path control augmentation using lateral control surfaces | |
GB1374103A (en) | Computers | |
CN103587681A (zh) | 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法 | |
CN107390708A (zh) | 一种基于鲁棒伺服控制无人机起飞拉起的方法 | |
CN112000127B (zh) | 一种基于反步法的飞行器横侧向联合控制方法 | |
CN111832115B (zh) | 对称迟后-超前校正的一体化频率设计方法 | |
JPH0738128B2 (ja) | 制御装置 | |
CN114942649B (zh) | 一种基于反步法的飞机俯仰姿态与航迹角解耦控制方法 | |
CN110597281A (zh) | 一种自动着舰纵向飞行控制系统参数获取方法 | |
CN113428219A (zh) | 一种基于传递函数快速响应控制汽车安全的系统和方法 | |
CN112325710B (zh) | 一种运载火箭大推力直接入轨高精度姿态控制方法和系统 | |
CN102935862A (zh) | 电动式助力转向系统的转向控制装置 | |
US2567922A (en) | Automatic pilot for aircraft | |
CN112372639B (zh) | 一种机械臂系统的二重相位超前补偿迭代学习控制方法 | |
CN113485394B (zh) | 一种高精度固定时间收敛的相对姿态容错跟踪控制方法 | |
JP3028888B2 (ja) | オートパイロット装置 | |
CN109653885B (zh) | 一种发动机实际转速均值配平方法及其控制系统 | |
CN113848960B (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 |