一种动中通天线系统辨识方法技术方案

本发明专利技术公开了一种动中通天线系统辨识方法，属于动中通天线伺服控制系统辨识技术领域。其首先生成伪随机序列信号作为天线辨识输入，控制天线转动，采集辨识过程中给定速度、天线轴角；然后对采集数据进行处理，得到辨识的输入输出矩阵，根据离散状态方程递推关系，由输入输出矩阵得到Hankel矩阵，提取子矩阵，进行奇异值分解，确定辨识模型的阶次；最后，根据模型阶次构建新的奇异值矩阵，求解待辨识天线的离散状态方程。该方法实时性高，辨识简便，不需要反复试凑确定模型阶次，辨识结果为离散模型，辨识结果受噪声影响小，模型准确度高。

An identification method for moving through antenna system

The invention discloses a method for identifying a moving through antenna system, which belongs to the technical field of identifying a servo control system of the moving through antenna. Firstly, pseudo-random sequence signals are generated as input of antenna identification, which controls the rotation of antenna and collects the given speed and antenna axis angle in the identification process. Then, the input-output matrix of identification is obtained by processing the collected data. According to the recursive relation of discrete state equation, the Hankel matrix is obtained from the input-output matrix, and the extractor is obtained from the input-output matrix. The singular value decomposition is used to determine the order of the identification model. Finally, a new singular value matrix is constructed according to the order of the model to solve the discrete state equation of the antenna to be identified. The method has the advantages of high real-time performance, simple identification, no need to repeatedly try to determine the order of the model, and the identification result is discrete model. The identification result is less affected by noise and the accuracy of the model is high.

【技术实现步骤摘要】
一种动中通天线系统辨识方法
本专利技术涉及动中通天线伺服控制系统辨识
，特别是指一种动中通天线系统辨识方法。
技术介绍
动中通天线伺服控制系统主要包括电机、驱动器、减速机、齿轮传动、天线面及馈源网络，由于负载转动惯量、传动链阻尼、齿隙等无法准确获取，根据伺服系统工作原理直接建模得到准确的数学模型非常困难，因此现有技术中通常采用系统辨识方法进行建模。目前，动中通天线系统辨识过程分为生成辨识信号、编写程序、采集数据、选择模型形式、参数拟合、模型误差对比、模型修正等过程，反复迭代，直到模型误差满足要求。采用上述方式的辨识过程存在如下不足：1)辨识过程不连续，信号生成、数据采集、数据处理在多个软件平台之间切换，工作量大。2)模型形式、模型阶次需要预设，系统辨识主要辨识的是预设传递函数的参数，无法体现模型其他信息。3)根据辨识误差，反复修正模型阶次，直到误差满足要求，工作效率不高。4)仅适用于较低阶次模型的系统辨识，模型阶次较高时辨识误差较大。
技术实现思路
有鉴于此，本专利技术的目的在于避免上述
技术介绍
中的不足之处而提供的一种动中通天线系统辨识方法，其具有使用操作方便，辨识模型阶次准确、辨识精度高的特点。为了实现上述目的，本专利技术所要解决的技术问题是由以下技术方案实现的：一种动中通天线系统辨识方法，该方法用于获得由状态空间描述的动中通天线的数学模型，包括以下步骤：(1)给定辨识激励信号；具体方式为，生成伪随机序列，将伪随机序列作为半实物仿真平台的辨识输入；根据电机减速比和齿轮传动比，计算电机驱动器速度，将电机驱动器速度作为给定速度通过CAN总线发送到天线的电机驱动器，控制天线电机按照该给定速度转动；(2)辨识数据采集；具体方式为，采集并记录天线的给定速度、方位角度和俯仰角度；(3)数据预处理；具体方式为，从辨识采集的数据中剔除野值，根据辨识输入输出矩阵定义，将预处理后的输入输出数据代入到矩阵中，得到辨识输入矩阵U和输出矩阵Y；(4)确定模型阶次；具体方式为，根据状态空间离散状态方程的伪随机序列响应关系Y＝HU，通过矩阵运算求解Hankel矩阵H；选择期望的模型最大阶次，构建子矩阵H1，对子矩阵H1进行奇异值分解：H1＝VΓ2UT，Γ为对角矩阵，其对角元素是H1的奇异值，V、U为正交矩阵；根据奇异值的大小确定辨识模型阶次k；(5)确定状态方程矩阵；具体方式为，选择矩阵V的前k行、矩阵U的前k列，矩阵Γ的k行k列，构构建新的矩阵V、U、Γ，得到能观矩阵P、能控矩阵Q：P＝VΓ,Q＝ΓUT，选择矩阵H的2～k+1行、2～k+1列，构建子矩阵H2，得到状态方程矩阵A、B、C，其中，A＝P-1H2Q-1，B为矩阵Q的第一列，C为矩阵P的第一行；(6)对动中通天线进行建模，天线的状态空间表示如下：xi+1＝Axi+BUiyi＝Cxi其中，xi为状态变量，yi为输出变量；完成动中通天线系统辨识建模。可选的，所述伪随机序列由Matlab软件生成，所述伪随机序列的参数设置为：时钟周期0.01秒，信号长度2047，重复周期4，信号幅度2°/s。可选的，所述步骤(4)和步骤(5)中奇异值分解的方式为采用Matlab软件中的SVD函数进行分解。本专利技术与
技术介绍
相比具有如下优点：1)本专利技术利用半实物仿真平台进行天线系统辨识，可以结合Matlab和半实物仿真平台的接口，便于利用Matlab强大的数学功能，具有编程简单、操作使用方便的优点。2)本专利技术通过输入输出数据构建Hankel矩阵，通过矩阵分解得到的奇异值从而判定模型阶次，不需要反复试凑确定模型阶次。3)本专利技术采用半实物仿真平台构成高速CAN总线数据采集系统，实时性高，保证了数据的准确性。附图说明图1是本专利技术实施例中动中通天线系统辨识平台的结构框图；图2是本专利技术实施例中动中通天线系统辨识过程的流程图。具体实施方式一种动中通天线系统辨识方法，该方法用于获得由状态空间表示的天线的数学模型。状态空间描述包括状态方程和输出方程，状态方程是描述系统状态变量和输入变量之间关系的一阶向量微分方程或差分方程，输出方程是描述系统输出变量与系统状态变量和输入变量之间函数关系的代数方程。状态空间表示如下：xi+1＝Axi+BUiyi＝Cxi，其中，矩阵A为系统矩阵或状态矩阵，B为控制矩阵或输入矩阵，矩阵C为输出矩阵或者观测矩阵。该辨识方法的主要任务就是完成对矩阵A、B、C的求解。该方法包括以下步骤：①给定辨识激励信号；具体方式为，选取伪随机信号作为辨识输入，伪随机信号由Matlab生成；根据减速机减速比和齿轮传动比，计算电机转速，通过CAN总线发送到电机驱动器，控制天线按照该期望速度转动；②辨识数据采集；具体方式为，采集并记录天线的给定速度、方位角度和俯仰角度；③数据预处理；具体方式为，从辨识采集的数据中剔除野值，根据辨识输入输出矩阵定义，计算输入矩阵、输出矩阵内部元素数值，构建辨识的输入矩阵U和输出矩阵Y；④确定模型阶次；具体方式为，离散状态方程零初始状态下的脉冲响应可以表示成Y＝HU，利用步骤③中的输入矩阵U和输出矩阵Y，通过矩阵运算求解Hankel矩阵H；接着，选择期望的模型最大阶次，可选的期望模型阶次不超过40，构建子矩阵，并对子矩阵进行奇异值分解，根据奇异值大小确定辨识模型阶次；⑤确定状态方程模型；具体方式为，状态方程为平衡实现形式，根据模型阶次构建新的Hankel矩阵，进行奇异值分解，得到能观矩阵P、能控矩阵Q，根据平衡实现和Hankel矩阵关系，求解矩阵A、B、C，完成动中通天线系统辨识建模。具体来说，如图1和2所示，一种动中通天线系统辨识方法，其针对的动中通天线可以为AE座架、AEC座架或者ACE座架天线，各轴系相互正交，耦合度低，辨识过程可以只考虑单个轴系，将天线系统作为单输入单输出(SISO)系统。该方法采用半实物仿真机AND5436、Matlab软件及Simulink工具搭建系统辨识平台，其中，由Matlab软件编制伪随机序列信号生成、数据处理、辨识算法等，Simulink工具搭建天线运动控制及数据采集模块，全部软件在AND5436仿真机内运行。AND5436仿真机运行程序，完成辨识数据采集后，由Matlab进行数据处理，构造输入、输出矩阵。下一步需要确定模型阶次，具体方式为：根据离散状态方程伪随机序列信号响应，由输入矩阵和输出矩阵解算Hankel矩阵H，进一步构造新的Hankel矩阵H1，对矩阵H1进行奇异值分解，确定重要奇异值数量，从而确定模型阶次。最后，根据模型阶次，构建新的奇异值分解子矩阵，根据平衡实现和Hankel矩阵关系，解算状态方程中的待定矩阵A、B、C，从而最终确认待辨识天线的状态方程。具体步骤如下：①给定激励信号系统辨识时，首先给定天线激励信号，天线按照给定速度转动。本例选取伪随机序列信号作为辨识输入，伪随机序列具有良好的随机性且其自相关性和功率谱是对白噪声的近似，能够很好地削减测量噪声对辨识的影响，是用于辨识的理想激励。伪随机信号参数有时钟周期，信号长度，重复周期数以及信号幅度等。动中通天线口径一般在0.3米～1.2米之间，谐振频率在10Hz～45Hz之间，根据动中通天线的结构和伺服系统特性，伪随机序列信号的参数选择为时钟周期0.01s，信号长度2047，重本文档来自技高网...

【技术保护点】
1.一种动中通天线系统辨识方法，其特征在于，用于获得由状态空间描述的动中通天线的数学模型，包括以下步骤：(1)给定辨识激励信号；具体方式为，生成伪随机序列，将伪随机序列作为半实物仿真平台的辨识输入；根据电机减速比和齿轮传动比，计算电机驱动器速度，将电机驱动器速度作为给定速度通过CAN总线发送到天线的电机驱动器，控制天线电机按照该给定速度转动；(2)辨识数据采集；具体方式为，采集并记录天线的给定速度、方位角度和俯仰角度；(3)数据预处理；具体方式为，从辨识采集的数据中剔除野值，根据辨识输入输出矩阵定义，将预处理后的输入输出数据代入到矩阵中，得到辨识输入矩阵U和输出矩阵Y；(4)确定模型阶次；具体方式为，根据状态空间离散状态方程的伪随机序列响应关系Y＝HU，通过矩阵运算求解Hankel矩阵H；选择期望的模型最大阶次，构建子矩阵H1，对子矩阵H1进行奇异值分解：H1＝VΓ2UT，Γ为对角矩阵，其对角元素是H1的奇异值，V、U为正交矩阵；根据奇异值的大小确定辨识模型阶次k；(5)确定状态方程矩阵；具体方式为，选择矩阵V的前k行、矩阵U的前k列，矩阵Γ的k行k列，构构建新的矩阵V、U、Γ，得到能观矩阵P、能控矩阵Q：P＝VΓ,Q＝ΓUT，选择矩阵H的2～k+1行、2～k+1列，构建子矩阵H2，得到状态方程矩阵A、B、C，其中，A＝P‑1H2Q‑1，B为矩阵Q的第一列，C为矩阵P的第一行；(6)对动中通天线进行建模，天线的状态空间表示如下：xi+1＝Axi+BUiyi＝Cxi其中，xi为状态变量，yi为输出变量；完成动中通天线系统辨识建模。...

【技术特征摘要】
1.一种动中通天线系统辨识方法，其特征在于，用于获得由状态空间描述的动中通天线的数学模型，包括以下步骤：(1)给定辨识激励信号；具体方式为，生成伪随机序列，将伪随机序列作为半实物仿真平台的辨识输入；根据电机减速比和齿轮传动比，计算电机驱动器速度，将电机驱动器速度作为给定速度通过CAN总线发送到天线的电机驱动器，控制天线电机按照该给定速度转动；(2)辨识数据采集；具体方式为，采集并记录天线的给定速度、方位角度和俯仰角度；(3)数据预处理；具体方式为，从辨识采集的数据中剔除野值，根据辨识输入输出矩阵定义，将预处理后的输入输出数据代入到矩阵中，得到辨识输入矩阵U和输出矩阵Y；(4)确定模型阶次；具体方式为，根据状态空间离散状态方程的伪随机序列响应关系Y＝HU，通过矩阵运算求解Hankel矩阵H；选择期望的模型最大阶次，构建子矩阵H1，对子矩阵H1进行奇异值分解：H1＝VΓ2UT，Γ为对角矩阵，其对角元素是H1的奇异值，V、U为正交矩阵；根据...

【专利技术属性】
技术研发人员：闫少雄，杜要锋，秦超，孙孟林，
申请(专利权)人：中国电子科技集团公司第五十四研究所，
类型：发明
国别省市：河北,13