闭环系统可辨识性实验研究..docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
研1206 刘新菊 模式识别与智能系统 2012020176 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 1 实验三 闭环系统可辨识性实验研究 实验目的:通过实验掌握闭环系统可辨识性条件及其应用。 实验内容: 模拟如下闭环系统 前向通道方程 其中v(k)是服从正态分布N(0,1)的高斯白噪声。 反馈通道方程 二阶反馈且存在反馈通道噪声 其中w(k)是服从正态分布N(0,1)的高斯白噪声。 二阶反馈但不加反馈噪声 一阶反馈并加反馈噪声 一阶反馈不加反馈噪声 比例反馈加反馈噪声 比例反馈不加反馈噪声 在以上各种不同的闭环状态下,利用最小二乘法直接辨识方法估计方法估计前向通道的模型参数。 根据实验结果分析闭环系统的可辨识性条件。 实验步骤: 熟悉闭环系统的可辨识性概念及条件。 设计实验方案,编制试验程序,上机调剂,记录结果。 分析实验结果,完成书面实验报告,要求同时提供程序框图及其说明。 实验三 闭环系统可辨识性实验报告 一、熟悉闭环系统的可辨识性概念及条件 系统辨识中的闭环问题作为一类特殊的辨识问题,近年来越来越受到人们的关注。Soderstrom?指出,由于输出信号的干扰噪声通过反馈环节与输入信号相关,直接采用频谱分析法,辨识结果将是对象传递函数与反馈传递函数倒数间的一个加权平均值;Ljung?的理论分析表明,如果采用预报误差法进行开环辨识时,只要对象模型集包含真实对象的动态特性,即使噪声模型不足以描述噪声的真实动态特性,仍可获得对象的一致无偏估计,而将此方法直接用于闭环辨识时,只要噪声模型不能精确描述真实噪声,即使对象模型集包含了真实对象动态特性,得到的将是对象参数的有偏估计;Gustavsson?则举出了反馈环节是比例调节器时,直接采用预报误差法将导致对象不可辨识性的实例。? ?对上述问题的解决,存在三种经典的解决方法:直接法、间接法、联合输入输出法。? (1)直接辨识法。直接法是利用闭环条件下得到的数据,直接把对象当作在开环环境中进行辨识。但要保证系统的可辨识性及对象参数估计的一致无偏性,采用直接法时必须满足条件:①存在足够的外部激励信号;②控制器的阶数足够高;③控制器在不同的模式间切换;④控制器时变、非线性。? (2)间接辨识法。间接法的思想是先根据闭环系统的输入输出信号用开环辨识方法得到闭环系统的传递函数,再由已知的控制器传递函数推导出对象的辨识模型。? (3)联合输入输出法。将闭环系统的输入输出统一看作为白噪声驱动系统产生的输出信号。先辨识出整个系统的传递函数矩阵,然后利用传递函数矩阵中的分块矩阵与系统各环节传递函数的关系,求出各环节的传递函数。? 经典闭环辨识法在许多方面并不尽如人意:直接法虽然简洁,但必须满足较为苛刻的条件;间接法与联合输入输出法实验前要知道反馈控制器的传递函数,并且需要至少三个可检测的信号,此外,闭环系统与开环对象模型间的参数转换导致了繁重的计算工作量。 1、二阶反馈且存在反馈通道噪声 %******************************* %闭环系统 %最小二乘辨识的递推算法 %Z(k)=-1.4*Z(k-1)-0.45*Z(k-2)+u(k-1)+0.7*u(k-2)+v(k) %u(k)=0.33*z(k)+0.033*z(k-1)-0.4*z(k)+w(k) %******************************* clear clc %**************产生均值为0,方差为1 的正态分布噪声************** v=random(Normal,0,1,1,402); w=random(Normal,0,1,1,402); %************产生观测序列z******************* z=zeros(402,1); u=zeros(402,1); z(1)=-1; z(2)=0; u(1)=-1; u(2)=0; for i=3:402 z(i)=-1.4*z(i-1)-0.45*z(i-2)+u(i-1)+0.7*u(i-2)+v(i); u(i)=0.33*z(i)+0.033*z(i-1)-0.4*z(i-2)+w(i); end %************递推求解************ P=100*eye(4); %估计方差 Pstore=zeros(4,401); Pstore(:,1)=[P(1,1),P(2,2),P(3,3),P(4,4)]; Theta=zeros(4,401); %参数的估计值,存放中间过程估值 Theta(:,1)=[3;3;3;3]; % K=zeros(4,400); %增益矩阵 K=[10;10;10;10

您可能关注的文档

文档评论(0)

jiayou10 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

版权声明书
用户编号:8133070117000003

1亿VIP精品文档

相关文档