一种风洞多模态实时系统辨识方法
【专利摘要】一种风洞多模态实时系统辨识方法，对当前风洞工况进行辨识，得到当前工况下的辨识模型，检测当前工况下的辨识模型能否对当前风洞工况进行预测和逼近；若不能，则重新采集数据，若能则实现对当前风洞工况的预测和逼近；当风洞工况变化后，实时检测当前风洞工况下的辨识模型能否对新风洞工况实现预测和逼近，若当前风洞工况下的辨识模型不能对新风洞工况实现预测和逼近，则重复步骤1)；随风洞工况不断变化，重复步骤2)实现整个风洞系统工况的预测和逼近。本发明在不需要断开运行的风洞系统的反馈环就可以进行系统输入输出数据采集，并对当前整个风洞工况进行系统辨识，所辨识系统能够反映风洞系统的整体特性，更有利于对当前风洞工况的分析。
【专利说明】
-种风洞多模态实时系统辨识方法
技术领域
[0001] 本发明属于系统辨识领域，尤其设及一种风桐多模态实时系统辨识方法。
【背景技术】
[0002] 常溫连续式跨声速风桐是一种W干燥空气作为介质的连续式跨声速空气动力实 验设备。该风桐主要用于解决我国新型战斗机、大型运输机、轰炸机、机动弹头、中远程导弹 和大型运载火箭W及可重复使用的天地往返运输系统的跨声速气动问题。其中参与流场控 制的系统包括:压缩机系统，进/排气系统，热交换器系统，二喉道中屯、调解片和再倒入段调 节片等。流场参数具有典型的多变量、非线性、强禪合等特点，与现有的低速风桐和亚音速 风桐W及生产型风桐相比其系统更为复杂。因此，得到此类风桐的有效且实用的控制模型 成为一个亟待解决的问题。
[0003] 近年来，余崇稽提出了运用有限元理论对跨声速风桐测力实验模型建立高维有限 元模型，并在此基础上建立了状态空间描述模型；冯凯提出运用机理建模结合参数辨识对 高超声速燃烧风桐加热器供气子系统建立了传递函数形式的模型，上述两种有限元模型和 传递函数形式的模型对风桐建模有一定的借鉴作用，但是考虑到常溫跨声速风桐系统的复 杂性，采用机理建模方法需要大量且合理的近似和简化才能完成。由于近似和简化存在，通 过运种方法所建模型与真实风桐系统误差可能较大，无法进行实际应用，因此运种方法不 适用于复杂系统建模。王锋，贺伟和毛鹏飞等对风桐系统的测力模型运用BP神经网络进行 建模。他们提出的风桐系统的辨识方法存在W下缺点：
[0004] (1)仅对风桐某局部系统或部分参数进行辨识，运种方法不足W反映风桐系统的 整体特性，如风桐系统是否具有稳定性W及是否具有能控能观测特性；
[0005] (2)实际常溫连续式跨声速风桐系统可能需要工作在不同工况，目前提出的方法 无法实时检测当前风桐工况辨识模型否能对新风桐工况进行预测和逼近，并在不能进行预 测和逼近时进行辨识模型切换。
【发明内容】
[0006] 针对【背景技术】所指出的风桐系统辨识存在的问题，本发明的目的在于提供一种风 桐多模态实时系统辨识方法，该方法具体采用DSR闭环系统辨识方法能够辨识出风桐不同 工况下的系统模型，此模型能够全面反映风桐系统的整体特性，并且能够保证系统辨识时 风桐实际运行的安全性和经济性；同时，实时检测当前风桐工况辨识模型否能对新风桐工 况进行预测和逼近，并在不能进行预测和逼近时进行辨识模型切换。
[0007] 为达到上述目的，本发明采用如下技术方案予W实现：
[000引一种风桐多模态实时系统辨识方法，包括W下步骤：
[0009] 1)对当前风桐工况采用DSR闭环子空间辨识方法进行辨识，得到当前工况下的辨 识模型，检测当前工况下的辨识模型能否对当前风桐工况进行预测和逼近;若不能，则重新 采集数据，若能则实现对当前风桐工况的预测和逼近；
[0010] 2)当风桐工况变化后，实时检测当前风桐工况下的辨识模型能否对新风桐工况实 现预测和逼近，若当前风桐工况下的辨识模型不能对新风桐工况实现预测和逼近，则重复 步骤1);
[0011] 随风桐工况不断变化，重复步骤2)实现整个风桐系统工况的预测和逼近。
[0012] 对当前风桐工况采用DSR闭环子空间辨识方法进行辨识，得到当前工况下的辨识 模型的具体过程如下：
[0013] (1)将白噪声信号输入风桐闭环系统，用信号采集仪采集白噪声输入信号和相应 的输出信号，在信号采样频率为10化，采样数据长度为2000的条件下，得到输入数据和输出 数据；
[0014] (2)对输入数据和输出数据进行奇异值剔除，然后数据归一化处理，得到处理后的 数据，根据处理后的数据得到输出信号的自相关曲线和功率谱密度曲线；
[0015] (3)对输出信号的自相关曲线和功率谱密度曲线进行检测，具体过程为:用输入信 号的自相关曲线和功率谱密度曲线与输出信号的自相关曲线和功率谱密度曲线做对比，检 测输出信号是否具有白噪声特性，如果具有白噪声的特性，则输出信号为白噪声信号；
[0016] (4)运用系统辨识方法对风桐当前工况进行辨识，得到当前工况下的辨识模型。
[0017] 检测当前工况下的辨识模型能否对当前风桐工况进行预测和逼近的具体过程如 下:运用残差序列白化检验方法对已辨识的系统进行预测和逼近能力测试，得到测试结果， 若辨识模型的残差序列不具有白噪声特性，则此辨识算法的辨识模型不能对当前风桐工况 进行预测和逼真;若模型的残差序列具有白噪声特性，则对应的辨识模型能够对当前风桐 工况进行预测和逼真。
[0018] 所述残差序列白化检验方法包括自相关曲线检验方法W及拟合残差序列的概率 密度方法。
[0019] 奇异值剔除具体采用拉依达奇异值剔除法。
[0020] 归一化处理具体采用Z-Score标准化数据归一化方法。
[0021] 实时检测当前风桐工况下的辨识模型能否对新风桐工况实现预测和逼近的具体 过程为：
[0022] (1)检测辨识系统输出序列值与真实系统输出序列值之间的匹配度是否大于 70%;
[0023] (2)检测辨识系统输出序列值与真实系统输出序列值的之间的残差序列是否白噪 声；
[0024] 若满足匹配度大于70%且残差序列是白噪声，则能够对对新风桐工况实现预测和 逼近，若不满足，则不能够对新风桐工况实现预测和逼近。
[0025] 与现有技术相比较，本发明的有益效果在于：
[0026] 1.本发明运用风桐输入输出数据直接对整个风桐工况进行系统辨识，系统辨识模 型具体为状态空间描述。由于常溫连续式跨声速风桐是个闭环系统如果采用开环辨识方 法，需要断开运行的风桐系统的反馈环才可进行系统输入输出数据采集，存在较大的不安 全因素及带来经济性问题。本发明运用DSR闭环子空间辨识方法在不需要断开运行的风桐 系统的反馈环就可W进行系统输入输出数据采集，并对当前整个风桐工况进行系统辨识， 由于是对整个风桐工况进行辨识，所W所辨识系统能够反映风桐系统的整体特性，更有利 于对当前风桐工况的分析。
[0027] 2.检测当前风桐工况辨识模型能否对当前风桐工况进行预测和逼近时，如果不能 当前辨识模型不能对当前风桐工况进行预测和逼近，可能是由于采样数据受到高频干扰， 需对采样数据进行低通滤波或进行数据重采样，并再次运用所述系统辨识方法进行风桐当 前工况辨识。
[00%] 3.实时检测当前风桐工况辨识模型能否对新风桐工况进行预测和逼近。如果当前 风桐工况辨识模型不能对新风桐工况进行预测和逼近，则需要进行辨识模型切换，重新利 用所述辨识方法进行新风桐工况下的系统辨识。根据不断变化的风桐工况运用所述系统辨 识方法和辨识系统检测方法对复杂风桐系统实现多模态辨识，更加符合实际应用。
[0029] 4.运用DSR闭环子空间辨识算法得到的辨识模型精度高于其他子空间辨识算法辨 识的模型。
【附图说明】
[0030] 图1是本发明单模态系统辨识流程图；
[0031 ]图2是本发明多模态系统辨识切换流程图；
[0032] 图3是DSR辨识系统与真实系统输出匹配度比较图；其中，（a)为Zl输出序列，（b)为 Z2输出序列，（C)为Z3输出序列。
[0033] 图4是MSID辨识系统与真实系统输出匹配度比较图；其中，（a)为Zl输出序列，（b) 为Z2输出序列，（C)为Z3输出序列。
[0034] 图5是PEM辨识系统与真实系统输出匹配度比较图；其中，（a)为Zl输出序列，（b)为 Z2输出序列，（C)为Z3输出序列。
[00对图6是DSR;输出残差序列自相关性检测曲线;其中，（a)为Zl，（b)为Z2,（C)为Z3。
[0036] 图7是MSIDS输出残差序列自相关性检测曲线；其中，（a)为Zl，（b)为Z2,（C)为 Z3。
[0037] 图8是阳输出残差序列自相关性检测曲线;其中，（a)为Zl，（b)为Z2,（C)为Z3。
[003引图9是DSRS输出残差序列概率密度估计曲线;其中，（a)为Zl，（b)为Z2,（C)为Z3。
[0039] 图10是MSIDS输出残差序列概率密度估计曲线；其中，（a)为Zl，（b)为Z2,（C)为 Z3。
[0040] 图11是阳MS输出残差序列概率密度估计曲线。其中，（a)为Zl，（b)为Z2,（C)为Z3。
【具体实施方式】
[0041 ]下面结合附图及【具体实施方式】对本发明作进一步详细说明。
[0042] 参见图1和图2,本发明一种风桐多模态实时系统辨识方法:风桐工况具体为风桐 模态，不同工况即为不同的模态，一种工况对应于一种模态，工况变化即模态发生变化;一 种风桐模态对应一种风桐小扰动模型。本发明中辨识系统用辨识模型进行表示。
[0043] 1)对当前风桐工况采用DSR闭环子空间辨识方法进行辨识，得到当前工况下的辨 识模型，检测当前工况下的辨识模型能否对当前风桐工况进行预测和逼近;若不能，则重新 采集数据，若能则实现对当前风桐工况的预测和逼近；
[0044] 2)当风桐工况变化后，实时检测当前风桐工况下的辨识模型能否对新风桐工况实 现预测和逼近，若当前风桐工况下的辨识模型不能对新风桐工况实现预测和逼近，则重复 步骤1);
[0045] 随风桐工况不断变化，重复步骤2)实现整个风桐系统工况的预测和逼近。
[0046] 其具体实施过程包括下述步骤：
[0047] 1.运用DSR闭环子空间辨识方法对整个风桐工况进行辨识，得到测试结果，具体包 括：
[0048] (1)将白噪声信号输入风桐闭环系统，用信号采集仪采集白噪声输入信号和相应 的输出信号，在信号采样频率为10化，采样数据长度为2000的条件下，得到输入数据和输出 数据；
[0049] (2)首先对输入数据和输出数据进行奇异值剔除，具体采用拉依达奇异值剔除法， 然后数据归一化处理，具体采用Z-score标准化数据归一化方法，得到处理后的数据，根据 处理后的数据得到输出信号的自相关曲线和功率谱密度曲线。
[0050] (3)对输出信号的自相关曲线和功率谱密度曲线进行检测。具体过程为:用输入信 号的自相关曲线和功率谱密度曲线与输出信号的自相关曲线和功率谱密度曲线做对比，如 果输出信号具有白噪声特性，则输出信号为白噪声信号；
[0051] (4)系统辨识。具体为运用S种系统辨识方法DSR，N4SID和TOM分别对风桐当前工 况进行辨识，得到当前工况下的辨识模型。
[0052] (5)对风桐当前工况下的辨识模型测试，具体过程为：
[0053] 对步骤(4)中=种辨识模型输出值和真实系统输出值进行匹配度（吻合度）比较。 为了能够清楚观察辨识系统与真实系统输出值之间的接近程度，现截取原2000组输出数据 中的50组数据(50~100)进行绘图，分别如图3~图5W及具体匹配度数值如表1所示。
[0054] 由图3~5可知，系统辨识算法MSID所辨识的模型性能远低于其他两种辨识算法 辨识模型，尤其体现在对输出量Z3的预测上；同时可知，DSR和阳M两种辨识算法所辨识的模 型能够对真实系统的输出进行很好的预测。
[0055] 表1 =种辨识系统输出值与真实系统输出值匹配度表
[0化6]
[0化7]
[0058] 由表1中S种算法辨识模型匹配度具体数值可知，可W得到与图3~图5相同结论， 同时，可知DSR算法的S输出匹配度比PEM算法的相应值要大，即DSR算法对所述风桐系统的 预测比阳M算法的精度要高。
[0059] 检测辨识系统输出序列值与真实系统输出序列值的之间的残差序列是否为白噪 声;若满足残差序列是白噪声，则当前风桐工况辨识模型能够对当前风桐工况进行预测和 逼近，若当前风桐工况辨识模型不能对当前风桐工况进行预测和逼近，可能是由于采样数 据受到高频干扰，需对采样数据进行低通滤波或进行数据重采样，并再次运用所述系统辨 识方法进行风桐当前工况辨识并检测。
[0060] 运用残差序列白化检验方法，具体为自相关曲线检验方法，W及拟合残差序列的 概率密度方法对已辨识的系统进行预测和逼近能力测试，得到测试结果，具体如图6~图 11，其中图6~图8为S输出残差序列自相关曲线，图9~图11为S输出残差序列概率密度估 计曲线。
[0061] 由图6~图8可得；种辨识模型输出残差序列的自相关曲线，其中由图7的自相关 曲线可知在横轴约[-10,10]时，对应的幅值都非零且非常大，由此可知MSID辨识算法辨识 模型的残差序列不具有白噪声特性，即此算法辨识的模型不能对原系统进行预测；同时，由 图6和图8可得，曲线仅在0时刻具有自相关特性，在其他时刻均接近零，即运两种辨识算法 辨识的模型的残差序列具有白噪声特性，也即对应的辨识模型能够对原始系统进行充分预 测。
[0062] 由图9~图11可得S种辨识模型输出残差序列的概率密度估计曲线，其中由图10 可知，N4SID辨识算法对应的残差序列概率密度曲线为非正态分布概率密度曲线，而由图9 和图11可知，其余两种辨识算法的相应残差序列的概率密度曲线均为正态分布概率密度曲 线，即运两种辨识算法辨识的模型的残差序列具有白噪声特性，也即对应的辨识模型能够 对原始系统进行充分预测。
[0063] 2.当风桐工况变化后，实时检测当前风桐工况下的辨识模型能否对新风桐工况实 现预测和逼近，具体过程如下：
[0064] (1)检测辨识系统输出序列值与真实系统输出序列值之间的匹配度是否大于 70%;
[0(?日]匹配度公式为
[0066]
(%) (1)
[0067] 其中Y为真实系统输出序列，f为辨识系统输出序列，norm( g)为求范数，mean (g) 为求均值。
[0068] 根据大量工程实践所得，一般情况下当匹配度大于70%时认为辨识模型合格。
[0069] (2)检测辨识模型输出序列值与真实系统输出序列值的之间的残差序列是否为白 噪声。
[0070] 运用残差序列白化检验方法对残差序列是否是白噪声进行检测，残差序列白化检 验方法包括自相关曲线检验方法W及拟合残差序列的概率密度方法。
[0071] 若满足匹配度大于70%且残差序列是白噪声，则当前风桐工况辨识模型能够对新 风桐工况进行预测和逼近;若不满足，则需要对新风桐工况进行系统辨识。
[0072] 随风桐工况不断变化，重复辨识和检测过程，实现整个风桐系统工况的预测和逼 近。
【主权项】
1. 一种风洞多模态实时系统辨识方法，其特征在于，包括以下步骤： 1) 对当前风洞工况采用DSR闭环子空间辨识方法进行辨识，得到当前工况下的辨识模 型，检测当前工况下的辨识模型能否对当前风洞工况进行预测和逼近;若不能，则重新采集 数据，若能则实现对当前风洞工况的预测和逼近； 2) 当风洞工况变化后，实时检测当前风洞工况下的辨识模型能否对新风洞工况实现预 测和逼近，若当前风洞工况下的辨识模型不能对新风洞工况实现预测和逼近，则重复步骤 1)； 随风洞工况不断变化，重复步骤2)实现整个风洞系统工况的预测和逼近。2. 根据权利要求1所述的一种风洞多模态实时系统辨识方法，其特征在于，对当前风洞 工况采用DSR闭环子空间辨识方法进行辨识，得到当前工况下的辨识模型的具体过程如下： (1) 将白噪声信号输入风洞闭环系统，用信号采集仪采集白噪声输入信号和相应的输 出信号，在信号采样频率为10Hz，采样数据长度为2000的条件下，得到输入数据和输出数 据； (2) 对输入数据和输出数据进行奇异值剔除，然后数据归一化处理，得到处理后的数 据，根据处理后的数据得到输出信号的自相关曲线和功率谱密度曲线； (3) 对输出信号的自相关曲线和功率谱密度曲线进行检测，具体过程为：用输入信号的 自相关曲线和功率谱密度曲线与输出信号的自相关曲线和功率谱密度曲线做对比，检测输 出信号是否具有白噪声特性，如果具有白噪声的特性，则输出信号为白噪声信号； (4) 运用系统辨识方法对风洞当前工况进行辨识，得到当前工况下的辨识模型。3. 根据权利要求1所述的一种风洞多模态实时系统辨识方法，其特征在于，检测当前工 况下的辨识模型能否对当前风洞工况进行预测和逼近的具体过程如下:运用残差序列白化 检验方法对已辨识的系统进行预测和逼近能力测试，得到测试结果，若辨识模型的残差序 列不具有白噪声特性，则此辨识算法的辨识模型不能对当前风洞工况进行预测和逼真;若 模型的残差序列具有白噪声特性，则对应的辨识模型能够对当前风洞工况进行预测和逼 真。4. 根据权利要求3所述的一种风洞多模态实时系统辨识方法，其特征在于，所述残差序 列白化检验方法包括自相关曲线检验方法以及拟合残差序列的概率密度方法。5. 根据权利要求3所述的一种风洞多模态实时系统辨识方法，其特征在于，奇异值剔除 具体采用拉依达奇异值剔除法。6. 根据权利要求3所述的一种风洞多模态实时系统辨识方法，其特征在于，归一化处理 具体采用Z-score标准化数据归一化方法。7. 根据权利要求1所述的一种风洞多模态实时系统辨识方法，其特征在于，实时检测当 前风洞工况下的辨识模型能否对新风洞工况实现预测和逼近的具体过程为： (1) 检测辨识系统输出序列值与真实系统输出序列值之间的匹配度是否大于70%; (2) 检测辨识系统输出序列值与真实系统输出序列值的之间的残差序列是否白噪声； 若满足匹配度大于70%且残差序列是白噪声，则能够对对新风洞工况实现预测和逼 近，若不满足，则不能够对新风洞工况实现预测和逼近。
【文档编号】G06F17/50GK106021693SQ201610323773
【公开日】2016年10月12日
【申请日】2016年5月16日
【发明人】蔡远利, 李红霞
【申请人】西安交通大学