个人上机仿真辨识作业.docVIP

个人上机仿真辨识作业 问题 CARMA模型的参数辨识。 其中： del 为延时，未知； na 为A多项式的阶次，未知； nb 为B多项式的阶次，未知； nc 为C多项式的阶次，未知； u 501*1的输入数据向量，已知； y 501*1的输出数据向量，已知。 要求画出辨识程序框图，编写系统结构辨识和参数辨识的程序，给出辨识结果图和数据，并比较分析，完成报告 辨识程序框图 由于题目中给定的是延时未知、阶次未知的模型，所以我编写的程序是先使用输入和输出数据的前301组数据进行模型阶和时滞的辨识。辨识程序的框图如图1所示。在辨识过程中，本题我采用的RELS算法和残差平方总和确定模型阶和时滞的方法。 图1 模型阶和时滞辨识程序框图 用后200组数据对所辨识出来的模型进行验证。验证的程序框图如图2所示。 图2 模型验证程序框图 辨识结果和数据 结果 模型阶和时滞的确定 本次作业题中，我使用的是残差平方总和的方法来进行模型阶和时滞的辨识的。其结果如图3所示。从图中可以看出，当d=3时，J(d)最小，可确定出模型的时滞为d=3；由于模型阶n≥d，又由J(n)曲线可以看出，n=3时，其曲线开始平滑，故取n=4。因此，通过辨识程序，以及对J(d)、J(n)曲线的分析，最后确定本题模型的阶为n=4，时滞为d=3。