1.一種歸一化的人機(jī)系統(tǒng)飛行品質(zhì)預(yù)測方法,其特征在于:具體包括如下步驟,
第一步,構(gòu)建增廣被控對象;
其中,xs是帶有延遲的增廣狀態(tài)向量,為xs的一階導(dǎo)數(shù),As、Bs、Cs、Ds、Es是增廣系數(shù)矩陣,up是飛行員操縱量,y是輸出狀態(tài)量,w是外界擾動向量。Cs=[C DCd],Cs=[C DCd],Ds=D,x為飛機(jī)小擾動方程狀態(tài)向量,xd為帶有時間延遲的增廣狀態(tài)向量,Ad、Bd、Cd分別為帶有時間延遲的增廣系數(shù)矩陣,A、B、C、D、E分別為飛機(jī)小擾動狀態(tài)方程特性系數(shù)矩陣;
飛行員實際感知向量yobs為
yobs=Csxs+Dsup+vy (2)
其中,vy是觀測噪聲,觀測噪聲強(qiáng)度Vy為:
其中,ρy是感知噪聲信噪比,fy是注意力分配系數(shù),為觀測噪聲方差;
第二步,根據(jù)飛行任務(wù),構(gòu)建飛行員最優(yōu)控制模型指標(biāo)函數(shù),給定初始加權(quán)系數(shù)值Qy和ru;
第三步,計算飛行員最優(yōu)控制增益;
第四步,求解觀測噪聲方差和操縱噪聲方差,計算Kalman濾波增益;
第五步,計算指標(biāo)函數(shù)J;
第六步,以指標(biāo)函數(shù)J為目標(biāo)函數(shù),以加權(quán)系數(shù)為自變量,進(jìn)行尋優(yōu)計算,直至找到令J取得最小值的加權(quán)系數(shù);
第七步,根據(jù)最優(yōu)注意力分配加權(quán)系數(shù),計算指標(biāo)函數(shù),進(jìn)行人機(jī)系統(tǒng)飛行品質(zhì)預(yù)測模型擬合;
根據(jù)飛行試驗結(jié)果對形如下式的飛行品質(zhì)預(yù)測模型進(jìn)行擬合:
其中的參數(shù)A和B為模型系數(shù),采用現(xiàn)有飛行品質(zhì)試驗數(shù)據(jù)擬合得到,
第八步,采用上述模型進(jìn)行飛行品質(zhì)預(yù)測。
2.根據(jù)權(quán)利要求1所述的一種歸一化的人機(jī)系統(tǒng)飛行品質(zhì)預(yù)測方法,其特征在于:第三步中所述的計算飛行員最優(yōu)控制增益,具體為:
通過最優(yōu)控制理論得到控制關(guān)系為:
其中,是帶有操縱量的狀態(tài)向量X的估計值,是飛行員最優(yōu)操縱量,Gp是調(diào)節(jié)器增益向量,K是由下列Riccati方程確定的唯一解
0=(Ao)TK+KAo+Qo-KBof-1(Bo)TK (6)
其中
將X=[xs up]T=[x xd up]T代入(5)式,則
其中Gn是狀態(tài)向量估計值增益向量,Gn1是的增益向量;令
則
Ip即為飛行員最優(yōu)控制增益,因此,式(8)寫為
令引入操縱噪聲vu,則,
其中vu是強(qiáng)度為Vu的零均值高斯白噪聲,其中ρu是操縱噪聲信噪比系數(shù),是操縱噪聲方差。
3.根據(jù)權(quán)利要求1或2所述的一種歸一化的人機(jī)系統(tǒng)飛行品質(zhì)預(yù)測方法,其特征在于:第四步具體為:
聯(lián)立(1)和(12)式得到:
其中,為狀態(tài)向量X的一階導(dǎo)數(shù),C1=[C DCd D],w為外界擾動向量,vu為操縱噪聲,vy為觀測噪聲;
上式中狀態(tài)向量X的估計值由Kalman濾波得到,其中濾波增益:
F=Σ1(C1)T(Vy)-1 (14)
其中估計誤差矩陣Σ1是由下列Riccati方程確定的唯一解:
0=A1Σ1+Σ1(A1)T+W1-Σ1(C1)T(Vy)-1(C1)Σ1 (15)
其中W1=diag(W,Vu),W為外界擾動強(qiáng)度,Vu為操縱噪聲強(qiáng)度。帶有狀態(tài)估計的人機(jī)閉環(huán)狀態(tài)方程為:
其中I1=[Ip,0],Cδ=[0 Cd 1],F(xiàn)為Kalman濾波增益矩陣,δ為飛機(jī)舵面偏轉(zhuǎn)量;
協(xié)方差矩陣Xcov是下列Lyapunov方程的解,
其中Qlyp=diag(Vw,Vy,Vu),Vy是觀測噪聲強(qiáng)度;
則,輸出協(xié)方差矩陣,
由此得到觀測噪聲方差和操縱噪聲方差分別為:
4.根據(jù)權(quán)利要求1所述的一種歸一化的人機(jī)系統(tǒng)飛行品質(zhì)預(yù)測方法,其特征在于:第五步中所述的指標(biāo)函數(shù)J:
其中,Je=Y(jié)cov(1,1),Ju=Y(jié)cov(2,2),row_u=row_X-row_x-row_xd,row_X是向量X的行數(shù),row_x是向量x的行數(shù),row_xd是向量xd的行數(shù)。