[發明專利]一種基于高精度多項式運算的復雜系統傳遞函數計算方法

有效

• 申請號: 201310349264.3
• 申請日: 2013-08-12
• 公開（公告）號: CN103400043A
• 公開（公告）日: 2013-11-20
• 發明（設計）人: 陳宇;楊云飛;李君;程興;鄧舞燕;陳彬;賈大玲;李聃;何巍;牟宇
• 申請（專利權）人: 北京宇航系統工程研究所;中國運載火箭技術研究院
• 主分類號: G06F19/00
• 分類號: G06F19/00;F42B15/01
• 代理公司: 核工業專利中心 11007
• 代理人: 包海燕

1.一種基于高精度多項式運算的復雜系統傳遞函數計算方法，其特征在于：包括以下步驟：?

步驟1輸入數據前處理；?

步驟2狀態空間矩陣計算；?

步驟3采用FADEEVA法計算箭體的傳遞函數；?

步驟4擴充浮點數據信息，進行高精度多項式運算；將慣導平臺通道和速率陀螺通道的兩個傳遞函數相加，得到箭體加平臺加陀螺的傳遞函數；將箭體加平臺加陀螺的傳遞函數與舵機的傳遞函數相乘，得到箭體加平臺加陀螺加舵機的傳遞函數；?

步驟5將箭體加平臺加陀螺加舵機的傳遞函數與較正網絡的傳遞函數相乘，得到姿控系統開環傳遞函數。?

2.根據權利要求1所述的一種基于高精度多項式運算的復雜系統傳遞函數計算方法，其特征在于：步驟2包括以下步驟：?

火箭姿態控制系統狀態方程表述為如下狀態方程形式：?

式中：?

x=[x₁,x₂,...,x_m]^T為m維狀態變量；?

y=[y₁,y₂,...,y_m]^T為m維輸出變量；?

u=[u₁,u₂,...,u_r]^T為r維輸入變量；?

A、B、C、D為狀態矩陣。?

3.根據權利要求2所述的一種基于高精度多項式運算的復雜系統傳遞函?數計算方法，其特征在于：步驟3包括以下步驟：?

純箭體傳遞函數輸入包括：合成擺角芯級發動機擺角和助推器發動機擺角輸出包括：剛體姿態角信號和加速度信號敏感元件測量的姿態角信號和測量得到的加速度信號

對步驟2所述火箭姿態控制系統狀態方程做拉氏變換，得：?

y(s)=C(sI-A)^-1Bu(s)+Du(s)；?

輸出變量到輸入變量的傳遞函數矩陣為：?

式中，矩陣G(s)的第(i,j)個元素為第i個輸出變量到第j個輸入變量的傳遞函數：?

FADEEVA法將預解矩陣(sI-A)^-1表示為：?

則其分子的系數矩陣和分母的系數按如下FADEEVA遞推公式計算：?

其中tr(A)表示A的跡。?