系统辨识大作业.pdfVIP

《系统辨识》大作业

学号：

班级：自动化1班

姓名：张青

信息与控制工程学院自动化系

2015-07-11

第一题

模仿indcx2,建对象，由相关分析法，获得脉冲响应序列反伏），由后供），参照讲

义,获得系统的脉冲传递函数Gz（）和传递函数G）；应用最小二乘辨识,获得脉冲响应序

列以幻；同图显示两种方法的辨识效果图：应用相关最小二乘法,拟合对象的差分方程模型;

构建时变对象，用最小二乘法和带遗忘因子的最小二乘法，可（以用辨识工具箱）辨识模型

的参数，比较两种方法的辨识效果差异；

答：根据index2建结构框图：

相关分析法：利用结构框图得到UY和tout

VariableEditor-UY田口,XWorkspace-*,n?X

屯蒲啕出■Base▼Selectdata...▼口一Xffl]r：画㈤Sel...▼»

◎小

田UY20x2doubleName▲Value

23456用UY20

ffitout00

15-60■

25-34.2558

35-4.494

45-2.347

5-54.7736

6-53.3802

7-5-2.4902

-

其中的U就是题目中要求得出的M序列，根据结构框图可知序列的周期是

N.=2-=24-=5。

在commandwindow中输入下列指令，既可以得到脉冲相应序列gk（）

aa=5;NNPP=15;ts=2;

RR=ones(15)+eye(15);

fori=15:-l:l

UU(16-i,:)=UY(16+i:30+i,l);

en

YY=[UY(31:45,2)];

GG=RR*UU*YY/[aa*aa*(NNPP+l)*ts];

plot(0:2:29,GG)

holon

stem(0:2:29,GG,fille)

最小二乘法建模的响应序列

由于是二阶水箱系统，可以假设系统的传递函数为G(s)=b/屏,