系统LH1.docVIP

课程实验报告 实 验 名 称 线性系统的经典辨识方法 所 在 学 院 控制工程学院 专 业 指 导 教 师 年月一、 实验名称：线性系统的经典辨识方法 实验目的：1、掌握系统辨识的基本步骤和主要内容。 2、熟悉Matlab下系统辨识的编程（M文件）。 3、实验一种经典辨识方法。 三、 实验设备：计算机、软件MATLAB 7.1 四、 实验原理：本实验采用脉冲响应法对连续系统的传递函数进行系统辨识，根据系统n阶差分方程当输入为脉冲响应后可得到输出脉冲响应的值而系统的传递函数 差分方程： 传递函数： 经拉氏逆变换得 ①式 同时得： ②式 ③式 把①②③式带入差分方程中得： 欲使等式等于零，不等于零只有当时才能成立 令则有成立 1、加入脉冲响应后得到各个采样输出值g(0)…g(t) 2、根据系统n阶差分方程，联络方程 由此可推出a1...an 3、将a1...an带入方程中可解出的值 4、将带入中可导出的值 5.根据已得的的值联络方程 其中可解出的值。 将的值和的值带入 即得出辨识的方程。 五、 实验结果：（含程序、解释、图形，结果和分析） %传递函数G(s)=/(s+3)*(s+7),展开为Gs=c1/(s-s1)-c2/(s-s2) t=[0:0.1:]; %设置采样间隔和采样范围 Gs=tf(12,[1 10 21]); %给出打算辨识的传递函数 M=impulse(Gs,t); %得到采样时间点的数值；将脉冲响应的值存入数组y中 A=[M(2) M(3);M(3) M(4)]; %设置矩阵A，其值为采样值的g2，g2，g3 B=[-M(1);-M(2)]; %设置矩阵B，其值为采样值的-g1,-g2 a=A\B; %得到a1，a2的值 N=[a(2) a(1) 1]; %设置矩阵N，其值为a2，a1，1 x=roots(N); %求解方程a2*x^2 + a1*x + 1 = 0，并将辨识的x1，x2放于x中 s1=10*log(x(1)); %求解s1的值 s2=10*log(x(2)); %求解s2的值 C=[1 1 ;x(1) x(2)]; %设置矩阵C，其值为g(0) = c1 + c2 中的系数，与g(1)=c1*x1+c2*x2中的系数 D=[-M(1);-M(2)]; %设置矩阵H，其值为采样值的-g1,-g2 c=C\D; %设置矩阵c，并将辨识的c1，c2的值放入其中 impulse(Gs,t); %绘制辨识的传递函数图形 grid; %在绘制的传递函数的图形中加入网格 图形： 数据： a1=-3.3636, a2=2.7183; s1=-3, s2=-7; c1=3, c2=3 2.实验结果： Transfer function: 12 G(s)= ------------- s^2 + 10 s + 21 辨识后的传递函数：G(s)=3/(s+3) - 3/(s+7) 3.结果分析：定义的待辨识的传递函与辨识出的传递函数相同故而实验方法正确，但是这是在未考虑信噪比没有误差的情况下进行的。所以实验结果无误差，需要改进的地方则为应加入误差，即信噪比的干扰。