资源目录
压缩包内文档预览:
编号:10288720
类型:共享资源
大小:1.66MB
格式:RAR
上传时间:2018-07-14
上传人:hon****an
认证信息
个人认证
丁**(实名认证)
江苏
IP属地:江苏
2.4
积分
- 关 键 词:
-
系统
辨识
讲稿
- 资源描述:
-
系统辨识讲稿,系统,辨识,讲稿
- 内容简介:
-
1系统辨识第 12 讲要点第 6 章 模型阶次的辨识(结构辨识)6.1 引言各种模型参数辨识方法一般需要假定模型的结构已知,实际上在多数情况下这是不现实的。当没有模型结构的先验知识时,需要利用系统的输入输出数据来确定模型的结构。这就是所谓的模型结构辨识问题。对单输入单输出(SISO) 系统来说,模型结构辨识也就是模型的阶次辨识。下面给出各种模型结构辨识方法(线性系统)。模型结构辨识的过程: 、 、 WinerVolta3216.2 根据 Hankel 矩阵的秩来估计模型的阶次 脉冲响应值的确定: 利用经典辨识方法中的相关分析法确定; 利用最小二乘类方法确定;设线性系统的输入输出序列分别为: ,脉Lkzku,21),(,冲响应值序列为 ,它们之间由卷积表示为:Nig,210),( )()0wiigkzi 2其中: 是过程的输出白噪声。)(kw令: )(,)2,1(,LwwNggzzL )()1()0)NLuuLH 则有: LLwgz利用最小二乘法可以获得脉冲响应的估计值: SzH1)( 利用脉冲响应值确定模型的阶次: 定义 Henkel 矩阵: lXkglkglkgll )2()()1(2)1()()(),( H其中, 为系统的脉冲响应值, 决定 Henkel 矩阵的维数。)(g注意: 的取值范围, 。k)2(lL Henkel 矩阵的性质:(此为定阶的依据) knlkl,ran00H 定阶准则(弱噪声的情形): lLklkl lD21),(det,3 强噪声下的 Henkel 矩阵:(强噪声的情形) lXklklkll )2()()1(2)1()()(),( H其中, , 。)0(gRkLigig16.3 根据行列式比确定模型的阶次6.3.1 无噪声情形假设模型为:(阶次参数为 )nnii ikubkya11 )()()(令: ,其中,nnUYH )()2()1( )()()( LyLnyny )()2()1( 1LuLnunuU 如果输入是持续激励信号,则有 ,关于矩阵 ,有以下性质:rakUnH,mi0Hn其中: 为过程模型的真实阶次。0n记:4nHL1)(则有以下结论:(判别模型阶次的依据) 0,)(det为了提高判别阶次的精度,一般用以下的判别式子: )1(det)(nHDR若 较 有显著增加时,则取阶次估计值为 。)(nDR)1( n06.3.2 白(弱)噪声情形假设模型为:(阶次参数为 )n)()()(11 kviubikzaznii 令: ,其中,nnUZH )()2()1( )1()()( LzLnznz zzz )()2()1( 1LuLnunuU 定义判别式子(行列式比): )1(det)(HDR其中: 1)(,1)( nnLLH5若 较 有显著增加时,则取阶次估计值为 。)(nDR)1( n06.3.3 有色噪声情形假设模型为: )()()(11 kvidikeceubzazmimi nini其中有两个阶次参数 。n,记: ,)(,)2(1,2121dcbadcaaLnvvvzzzmnn 及: ,其中:, mnmnVEUZH )()2()1( )1()()( LzLnznz zzz )()2()1( 1LuLnunuU 6 )()2()1( 2)1()()( LmneLneneEm )()2()1( 21LnvLnvnvVm 由此可得: nmnHz,依据最小二乘法估计参数 ,得:nnz,1,)(利用以下的估计式代替 中得元素,迭代计算得 。mVE, mVE,dvceba)()()(kkvzeu其中: )(,)2(,1)( , mkvkvkeenuuzz ve因此得到: ,nmnVEUZH记:7 )()( 11),(211,mHn VEUVZEZLmmnnnmnnnmn 定义判别式子(行列式比): )1(det)(nHnDR和 )(et)(2m即可对参数 进行估计。mn,为了减少计算量,可以用以下的判别式子: )1(det)(nDRt2m其中: )()(1221Hn6.4 利用残差的方差估计模型的阶次6.4.1 残差方差分析考虑如下模型 )()()(11 kvuzBkzA式中 u(k)和 z(k) 分别为模型输入和输出变量; v(k) 是均值为零、方差8为 不相关随机噪声; 和 为迟延算子多项式,记作2v)(1zA)(1zBaazBbbn() 211其中 n 为假设模型的阶次参数。模型的最小二乘格式可写成 zkvkn()()hT式中,数据向量和参数向量定义为 hnnuznuab(),),1 TT运用最小二乘原理,可获得模型参数 的最小二乘估计为n()nnHzTT1式中,数据矩阵和输出向量定义为 hznnLzLTT T()(),(),()1212其中, L 为数据长度。模型阶次为 n 时,输出残差向量可写成 nnvHz式中 (),()nnnvvL012 T且残差的方差具有如下性质 PP+PTTTlim()lilimliLLnLnLbnVn1v它会呈现如下的变化特性96.4.2 白噪声过程的 F-Test 定阶法通过观察残差方差 的变化情况,可确定模型的阶次,具体步骤:)(nV 阶次 n 逐一增加,当 n 增加至 时,若 呈现显著下降,)(nV则确认模型阶次为 ; 引进统计量 tVnLnL()()(),)1221212Ft 服从 F 分布,自由度为 2(n2-n1) 和 L-2n2; 设零假设 ,取风险水平 ,对应的阀值00:H%5。),(221Lnt 按下式判断模型阶次若 ,拒绝 ,若 ,接受 ,tt, 0:ntt1,( 001:nH其中 为风险水平 下的阀值。这时模型的阶次估计值可取 。6.4.3 有色噪声过程的 F-Test 定阶法考虑如下模型V(n)nv2n0显著下降点10 )()()( 11 ikvdkieckeiubizazmimi nini式中 u(k)和 z(k) 分别为模型输入和输出变量; v(k) 是均值为零、方差为 不相关随机噪声; n 和 m 为模型阶次。v2 引进统计量 )2,(2)1,()1,( ,)(),(, nLFLnVmt nn 设零假设 ,取风险水平0)(0)0 1: mHHm,对应的阀值 。%5,2(Ft 按下式判断模
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。