[发明专利]静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法

权利要求书

1.一种静不稳定飞行器等价气动伺服弹性鲁棒稳定性的飞行试验确定方法，其特征在于包括下述步骤：

(a)通过扫频飞行试验，获得静不稳定飞行器多回路系统含有不确定量的闭环传递函数频率特性矩阵Φ(jω)＝G(jω)[I+G(jω)]^-1；为了确定系统的相位裕度和幅值裕度，定义：Y_p(jω)＝K_pe^-τjω，并串联接至等价开环系统；

式中，j为虚数符号，ω表示频率，Φ(jω)为闭环传递函数频率特性矩阵，G(jω)为开环传递函数频率特性矩阵，I为单位矩阵，K_p为每条回路的附加的增益，τ为每条回路的附加时间滞后；

(b)当系统临界稳定时，复变量s的实部为零，令s＝jω为纯虚数，满足以下行列式关系|(K_pe^-τjω-1)Φ(jω)+I|＝0，

对Φ(jω)进行特征分解，得：Φ(jω)＝T(jω)D(jω)T^-1(jω)

式中，T(jω)为线性变换矩阵，D(jω)为约当阵；

(c)取τ＝0，选定颤振频率段[ω₀，ω₁]，求解满足以下等式：

Πi=1n|{1+(Kp-1)Re[di(jω)]}2+{(Kp-1)Im[di(jω)]}2|≤δ2]]>

时得到的最小ω值和回路幅值裕度K_p值，最小的K_p为整个系统在选定颤振频率段[ω₀，ω₁]的幅值裕度，δ＞0为不确定性影响的估计值，d_i(jω)为矩阵D(jω)的第i行第i列元素；

(d)相位裕度转化为K_p＝1时在选定颤振频率段[ω₀，ω₁]上计算

Πi=1n|{1+Re[(e-τjω-1)di(jω)]}2+Im{[(e-τjω-1)di(jω)]}2|≤λ2]]>

时得到的最小ω值和最小的τ，此时的ωτ为整体系统的相位裕度，λ＞0为不确定性影响的估计值。