系统辨识习题解答
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
。
由于 , ,
所以
由Markov估计, ,其中矩阵C为
,
取加权阵 。
P532/4:解:(1)由参数估计值偏差的估计式:
我们有:
(A)
由于 为独立同分布,均值为零的不相关随机变量,因此有:
对(A)式两边求期望值,我们有:
由此递归式子,可得: (B)
(2)由指标的估计式:
将初值 和(B)式代入,两边取对数,有:
证毕。
则噪声模型也化成最小二乘格式:
数据向量he(k)包含着不可测的噪声量,这可用相应的估计值代替:
其中,
则可写出利用增广最小二乘法得到的递推算法:
θ可表示成:
4-19、考虑如下模型
其中,u(k)和z(k)分别为模型的输入和输出变量,它们是可测的;v(k)是零均值白噪声,它是不可测的。试从Markov估计概念出发,证明该模型的参数向量 的估计值 可以写成如下加权最小二乘算法的形式 ,
解:根据表示定理,在一定条件下,有色噪声e(k)可以看成是由白噪声v(k)驱动的线性环节的输出,该线性环节称为成形滤波器,其脉冲传递函数可写成
即
其中
根据其结构,噪声模型可区分为以下三Байду номын сангаас:
自回归模型(AR模型):
平均滑动模型(MA模型):
自回归平均滑去模型(ARMA模型):
3-4、根据离散Wiener-Hopf方程,证明
(1)
(2) ,
(3) ,
(4) ,
解:(1)由于
,
所以
(2)由于
,
及
(3)由于
,
所以
(4)由于
,
所以
4-18、考虑如下模型
其中,u(k)和z(k)是模型的输入输出变量,v(k)是零均值白噪声。定义参数向量
请利用增广最小二乘思想,写出模型参数 的递推辨识算法。
解:令
及
则模型化成最小二乘格式:
令 ,及
式中, 为数据矩阵, 为输出向量,加权矩阵取 ,其中矩阵C为
解:令
及
则模型化成最小二乘格式:
准则函数取 ,其中 为加权因子,对所有的k, 都必须大于零。
对于 (L为数据长度),可以构成线性方程组
式中
则 ,式中 为加权矩阵,它是正定的对角阵,由加权因子 构成
,
设 使得J(θ)最小,则有:
从而:
,
当 是正则矩阵时,模型的加权最小二乘解为
系统辨识习题解答
1-14、若一个过程的输入、输出关系可以用MA模型描述,请将该过程的输入输出模型写成最小二乘格式。
提示:①MA模型
② 定义
解:因为MA模型 ,其中
,从而
所以当定义 ,则有最小二乘格式:
,
其中e(k)是误差项。
2-3、设 是一个平稳的有色噪声序列,为了考虑这种噪声对辨识的影响,需要用一种模型来描述它。请解释如何用白噪声和表示定理把 表示成AR模型、MA模型和ARMA模型。
解:由于M序列是循环周期为 , , 为 序列移位脉冲周期,自相关函数近似于 函数, 为 序列的幅度。设数据的采样时间等于 ,则离散Wiener-Hopf方程为:
当 序列的循环周期 大于过程的过渡过程时间时,即 充分大时,离散Wiener-Hopf方程可写成:
由于M序列的自相关函数为
,
代入上式得
4- 证明:
由于 , ,
所以
由Markov估计, ,其中矩阵C为
,
取加权阵 。
P532/4:解:(1)由参数估计值偏差的估计式:
我们有:
(A)
由于 为独立同分布,均值为零的不相关随机变量,因此有:
对(A)式两边求期望值,我们有:
由此递归式子,可得: (B)
(2)由指标的估计式:
将初值 和(B)式代入,两边取对数,有:
证毕。
则噪声模型也化成最小二乘格式:
数据向量he(k)包含着不可测的噪声量,这可用相应的估计值代替:
其中,
则可写出利用增广最小二乘法得到的递推算法:
θ可表示成:
4-19、考虑如下模型
其中,u(k)和z(k)分别为模型的输入和输出变量,它们是可测的;v(k)是零均值白噪声,它是不可测的。试从Markov估计概念出发,证明该模型的参数向量 的估计值 可以写成如下加权最小二乘算法的形式 ,
解:根据表示定理,在一定条件下,有色噪声e(k)可以看成是由白噪声v(k)驱动的线性环节的输出,该线性环节称为成形滤波器,其脉冲传递函数可写成
即
其中
根据其结构,噪声模型可区分为以下三Байду номын сангаас:
自回归模型(AR模型):
平均滑动模型(MA模型):
自回归平均滑去模型(ARMA模型):
3-4、根据离散Wiener-Hopf方程,证明
(1)
(2) ,
(3) ,
(4) ,
解:(1)由于
,
所以
(2)由于
,
及
(3)由于
,
所以
(4)由于
,
所以
4-18、考虑如下模型
其中,u(k)和z(k)是模型的输入输出变量,v(k)是零均值白噪声。定义参数向量
请利用增广最小二乘思想,写出模型参数 的递推辨识算法。
解:令
及
则模型化成最小二乘格式:
令 ,及
式中, 为数据矩阵, 为输出向量,加权矩阵取 ,其中矩阵C为
解:令
及
则模型化成最小二乘格式:
准则函数取 ,其中 为加权因子,对所有的k, 都必须大于零。
对于 (L为数据长度),可以构成线性方程组
式中
则 ,式中 为加权矩阵,它是正定的对角阵,由加权因子 构成
,
设 使得J(θ)最小,则有:
从而:
,
当 是正则矩阵时,模型的加权最小二乘解为
系统辨识习题解答
1-14、若一个过程的输入、输出关系可以用MA模型描述,请将该过程的输入输出模型写成最小二乘格式。
提示:①MA模型
② 定义
解:因为MA模型 ,其中
,从而
所以当定义 ,则有最小二乘格式:
,
其中e(k)是误差项。
2-3、设 是一个平稳的有色噪声序列,为了考虑这种噪声对辨识的影响,需要用一种模型来描述它。请解释如何用白噪声和表示定理把 表示成AR模型、MA模型和ARMA模型。
解:由于M序列是循环周期为 , , 为 序列移位脉冲周期,自相关函数近似于 函数, 为 序列的幅度。设数据的采样时间等于 ,则离散Wiener-Hopf方程为:
当 序列的循环周期 大于过程的过渡过程时间时,即 充分大时,离散Wiener-Hopf方程可写成:
由于M序列的自相关函数为
,
代入上式得
4- 证明: