[發(fā)明專(zhuān)利]縱向飛行模型簇復(fù)合頻率控制器設(shè)計(jì)方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410070335.0
• 申請(qǐng)日: 2014-02-28
• 公開(kāi)（公告）號(hào): CN103809457A
• 公開(kāi)（公告）日: 2014-05-21
• 發(fā)明（設(shè)計(jì)）人: 史忠科
• 申請(qǐng)（專(zhuān)利權(quán)）人: 西安費(fèi)斯達(dá)自動(dòng)化工程有限公司
• 主分類(lèi)號(hào): G05B13/04
• 分類(lèi)號(hào): G05B13/04
• 代理公司: 暫無(wú)信息
• 代理人: 暫無(wú)信息

1.一種縱向飛行模型簇復(fù)合頻率控制器設(shè)計(jì)方法，其特點(diǎn)是包括以下步驟：

1)給定不同高度、馬赫數(shù)下通過(guò)掃頻飛行試驗(yàn)直接由允許飛行的全包線(xiàn)內(nèi)的幅頻和相頻特性構(gòu)成飛行器全包線(xiàn)內(nèi)的升降舵與飛行高度之模型簇，對(duì)應(yīng)的飛行器升降舵與飛行高度之間開(kāi)環(huán)傳遞函數(shù)簇描述為：

G(s)=e-σ(h,M)sK(h,M)A(h,M,s)B(h,M,s)]]>

其中

A(h,M,s)=s^m+a_m-1(h,M)s^m-1+a_m-2(h,M)s^m-2+…+a₁(h,M)s+a₀(h,M)、

B(h,M,s)=sⁿ+b_n-1(h,M)s^n-1+b_n-2(h,M)s^n-2+…+b₁(h,M)s+b₀(h,M)為多項(xiàng)式，s為傳遞函數(shù)中常用的拉普拉斯變化后的變量，h,M分別為飛行高度和馬赫數(shù)，σ(h,M)是俯仰回路的延遲時(shí)間，K(h,M)為隨h,M變化的增益，a_l(h,M),l=0,1,2,…,m-1為多項(xiàng)式A(h,M,s)中隨h,M變化的系數(shù)簇，b_i(h,M),i=0,1,2,…,n-1為多項(xiàng)式B(h,M,s)中隨h,M變化的系數(shù)簇；

2)根據(jù)飛行包線(xiàn)內(nèi)的幅頻特性直接確定開(kāi)環(huán)截止頻率區(qū)間確定方法為：

從|G(jω)|=1即中，得到開(kāi)環(huán)截止頻率ω_c解的最大值ω_cmax和最小值ω_cmin，開(kāi)環(huán)截止頻率ω_c區(qū)間為ω_cmin≤ω_c≤ω_cmax;

式中，jω為頻率特性中的變量，j為虛部表示，ω為角頻率；

3)根據(jù)飛行包線(xiàn)內(nèi)的相頻特性，計(jì)算包線(xiàn)內(nèi)最大相位裕度γ_max(ω_c)=max{180°-σ(h,M)ω_c+arg[A(h,M,jω_c)]-arg[B(h,M,ω_c)]},ω_cmin≤ω_c≤ω_cmax和包線(xiàn)內(nèi)最小相位裕度

γ_min(ω_c)=min{180°-σ(h,M)ω_c+arg[A(h,M,jω_c)]-arg[B(h,M,ω_c)]},ω_cmin≤ω_c≤ω_cmax直接確定與截止頻率區(qū)間所對(duì)應(yīng)的相位裕度區(qū)間為：γ_min(ω_c)≤γ(ω_c)≤γ_max(ω_c),ω_cmin≤ω_c≤ω_cmax；

4)候選多級(jí)串聯(lián)滯后—超前補(bǔ)償環(huán)節(jié)的傳遞函數(shù)為：

Gc(s)=kcΠi=1NTD1(i)s+1a(i)TL1(i)s+1·TD2(i)s/a(i)+1TL2(i)s+1]]>

式中，k_c為待確定的常數(shù)增益，N為整數(shù)，表示待確定的滯后—超前補(bǔ)償環(huán)節(jié)的級(jí)數(shù)，T_D1(i)、T_L1(i)、T_D2(i)、T_L2(i),i=1,2,…,N為待確定的時(shí)間常數(shù)，a(i)>1,i=1,2,…,N為待確定的參數(shù)；

加入多級(jí)串聯(lián)滯后—超前補(bǔ)償環(huán)節(jié)后，

從|G(jω)G_c(jω)|=1即|[kcΣi=1NTD1(i)s+1a(i)TL1(i)s+1·TD2(i)s/a(i)+1TL2(i)s+1e-σ(h,M)SK(h,M)A(h,M,s)B(h,M,s)]s=jω|=1]]>中，得到開(kāi)環(huán)截止頻率ω_fc解的最大值ω_fcmax和最小值ω_fcmin，開(kāi)環(huán)截止頻率ω_fc區(qū)間為ω_fcmin≤ω_fc≤ω_fcmax，

在飛行器全包線(xiàn)內(nèi)的相位裕度指標(biāo)γ^*給定情況下，加入多級(jí)串聯(lián)滯后—超前補(bǔ)償環(huán)節(jié)后系統(tǒng)的相位裕度γ_f(ω_fc)應(yīng)該滿(mǎn)足：

γ_f(ω_fc)=180°+arg[G(jω_fc)]+arg[G_c(jω_fc)]>γ^*,ω_fcmin≤ω_fc≤ω_fcmax

即滿(mǎn)足：180°-σ(h,M)ω_fc+arg[A(h,M,jω_fc)]-arg[B(h,M,ω_fc)]+arg[G_c(jω_fc)]>γ^*

ω_fcmin≤ω_fc≤ω_fcmax

在上述指標(biāo)和極大似然準(zhǔn)則或其它準(zhǔn)則共同約束下，可以根據(jù)系統(tǒng)模型結(jié)構(gòu)辨識(shí)中的極大似然方法或辨識(shí)方法確定滯后—超前補(bǔ)償環(huán)節(jié)的級(jí)數(shù)N、常數(shù)增益k_c、時(shí)間常數(shù)T_D1(i)、T_L1(i)、T_D2(i)、T_L2(i),i=1,2,…,N，待確定的參數(shù)a(i)>1,i=1,2,…,N；

5)在飛行器全包線(xiàn)內(nèi)的幅值裕度指標(biāo)L*分貝數(shù)給定情況下，

從20log₁₀G(jω)G_c(jω)=-L^*即

20log10|[kcΣi=1NTD1(i)s+1a(i)TL1(i)s+1·TD2(i)s/a(i)+1TL2(i)s+1e-σ(h,M)SK(h,M)A(h,M,s)B(h,M,s)]s=jω|=-L*]]>中，得到頻率ω_Lc解的最大值ω_Lcmax和最小值ω_Lcmin，ω_Lc區(qū)間為ω_Lcmin≤ω_Lc≤ω_Lcmax，判斷：

γ_L(ω_Lc)=180°+arg[G(jω_Lc)]+arg[G_c(jω_Lc)]>0,ω_Lcmin≤ω_Lc≤ω_Lcmax

即滿(mǎn)足：180°-σ(h,M)ω_Lc+arg[A(h,M,jω_Lc)]-arg[B(h,M,ω_Lc)]+arg[G_c(jω_Lc)]>0

ω_Lcmin≤ω_Lc≤ω_Lcmax

若滿(mǎn)足，則飛行控制器設(shè)計(jì)完成，若不滿(mǎn)足，再增加補(bǔ)償環(huán)節(jié)級(jí)數(shù)或減小常數(shù)增益k_c。