[发明专利]静不稳定飞行器等价稳定裕度的飞行试验鲁棒确定方法

权利要求书

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

(c)取τ＝0，有：

|(K_p-1)Φ(jω)+I|＝0

对Φ(jω)进行特征分解，得：

Φ(jω)＝T(jω)D(jω)T^-1(jω)

则幅值裕度可以近似为计算行列式det[(K_p-1)D(jω)+I]的模值：

|det[(K_p-1)D(jω)+I]|≤δ，

或者Π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列元素；

式中，det为行列式符号，T(jω)为线性变换矩阵，D(jω)为约当阵；

(d)相位裕度转化为计算行列式det[I+(e^-τjω-1)D(jω)]的模值：

|det[I+(e^-τjω-1)D(jω)]|≤λ

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

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