24-25-2-系统辨识与建模-期末

说明:所有计算题需写出详细公式和步骤,结果保留三位小数。

一、简答题(共5小问,每题7分)

1.

说明什么是系统辨识并写出系统辨识的三要素。

2.

什么是白噪声和有色噪声,并说明有色噪声如何白化。

3.

解释“持续激励信号”的定义,并说明其在系统辨识中的作用。

4.

数据病态会导致最小二乘估计失效,分析其产生原因及解决方法。

5.

卡尔曼滤波中过程噪声 wkw_{k} 和观测噪声 vkv_{k} 的协方差矩阵 QkQ_{k}RkR_{k} 如何影响滤波性能?

二、(10 分)

实验测得某物体温度 yy(℃)随时间 tt(小时)的变化数据如下:

tt1234
yy57911

假设模型为 y=at+by = at + b

  1. 使用标准最小二乘法估计参数 aabb。(8 分)
  2. 计算残差平方和。(2 分)

三、(15 分)

三个传感器测量同一物理量 xx,建立观测模型 y=Hx+vy = Hx + v,测量数据如下:

  • 传感器1:y1=8.2y_1=8.2,测量噪声方差 σ12=0.25\sigma_1^2=0.25
  • 传感器2:y2=8.5y_2=8.5,测量噪声方差 σ22=0.16\sigma_2^2=0.16
  • 传感器3:y3=7.9y_3=7.9,测量噪声方差 σ32=0.64\sigma_3^2=0.64
  1. 给定权重矩阵 Σ=diag(1,2,3)\Sigma=\operatorname{diag}(1,2,3),求 xx 的加权最小二乘估计值。(7 分)
  2. 利用测量噪声协方差矩阵 R=diag(0.25,0.16,0.64)R = \operatorname{diag}(0.25, 0.16, 0.64),求 xx 的马尔科夫估计值。(8 分)

四、(25 分)

考虑连续时间系统:

G(s)=Ks+a(K=2,a=1)G ( s ) = \frac { K } { s + a } \quad ( K = 2 , \, a = 1 )

使用双线性变换法(s=2Tz1z+1s=\frac{2}{T}\cdot\frac{z-1}{z+1}),采样周期 T=1T=1 秒。

1. 系统离散化(10分)

G(s)G(s) 转换为离散传递函数 G(z)G(z),推导其表达式并表示为:

G(z)=b0+b1z11+a1z1G ( z ) = \frac { b _ { 0 } + b _ { 1 } z ^ { - 1 } } { 1 + a _ { 1 } z ^ { - 1 } }

a1a_{1}b0b_{0}b1b_{1} 的数值表达式。

2. 递推最小二乘参数估计(15 分)

实验测得输入输出数据如下:

kk0123
u(k)u(k)1.20.60.30.3
y(k)y(k)0.01.20.90.7

假设模型结构为:

y(k)=a1y(k1)+b0u(k)+b1u(k1)+e(k)y ( k ) = - a _ { 1 } y ( k - 1 ) + b _ { 0 } u ( k ) + b _ { 1 } u ( k - 1 ) + e ( k )

初始估计为 θ^0=[0,0,0]\hat { \theta } _ { 0 } = [ 0 , 0 , 0 ] ^ { \top },协方差矩阵 P0=1000IP _ { 0 } = 1 0 0 0 I

要求:

  1. 写出递推最小二乘的 KkK_{k}θ^k\hat{\theta}_{k}PkP_{k} 更新公式。(5 分)
  2. 递推计算 k=1,2,3k=1,2,3 时的参数估计 θ^k\hat{\theta}_{k}。(7 分)
  3. 验证最终估计值是否与(1)中理论离散化参数是否一致,若不一致,可能原因有哪些。(3分)

五、(15 分)

离散系统动态方程为:

xk=[1101]xk1+wk1,wkN(0,Q),Q=[0.1000.1]x _ { k } = \left[ \begin{array} { l l } { 1 } & { 1 } \\ { 0 } & { 1 } \end{array} \right] x _ { k - 1 } + w _ { k - 1 } , \quad w _ { k } \sim \mathcal { N } ( 0 , Q ) , \quad Q = \left[ \begin{array} { l l } { 0 . 1 } & { 0 } \\ { 0 } & { 0 . 1 } \end{array} \right]

观测方程为:

yk=[10]xk+vk,vkN(0,R),R=0.5y _ { k } = [ 1 \quad 0 ] x _ { k } + v _ { k } , \quad v _ { k } \sim \mathcal { N } ( 0 , R ) , \quad R = 0 . 5

已知初始状态 x^00=[00],P00=Ic\hat{x}_{0|0}=\left[\begin{array}{l}0\\0\end{array}\right], P_{0|0}=I_{c}

1.

计算一步预测 x^10\hat{x}_{1|0} 和预测协方差 P10P_{1|0};(5 分)

2.

计算卡尔曼增益 K1K_{1};(5 分)

3.

y1=2.0y_{1}=2.0,更新状态估计 x^11\hat{x}_{1|1} 和协方差 P11P_{1|1}。(5 分)