说明:所有计算题需写出详细公式和步骤,结果保留三位小数。
一、简答题(共5小问,每题7分)
1.
说明什么是系统辨识并写出系统辨识的三要素。
2.
什么是白噪声和有色噪声,并说明有色噪声如何白化。
3.
解释“持续激励信号”的定义,并说明其在系统辨识中的作用。
4.
数据病态会导致最小二乘估计失效,分析其产生原因及解决方法。
5.
卡尔曼滤波中过程噪声 wk 和观测噪声 vk 的协方差矩阵 Qk 和 Rk 如何影响滤波性能?
二、(10 分)
实验测得某物体温度 y(℃)随时间 t(小时)的变化数据如下:
假设模型为 y=at+b。
- 使用标准最小二乘法估计参数 a 和 b。(8 分)
- 计算残差平方和。(2 分)
三、(15 分)
三个传感器测量同一物理量 x,建立观测模型 y=Hx+v,测量数据如下:
- 传感器1:y1=8.2,测量噪声方差 σ12=0.25
- 传感器2:y2=8.5,测量噪声方差 σ22=0.16
- 传感器3:y3=7.9,测量噪声方差 σ32=0.64
- 给定权重矩阵 Σ=diag(1,2,3),求 x 的加权最小二乘估计值。(7 分)
- 利用测量噪声协方差矩阵 R=diag(0.25,0.16,0.64),求 x 的马尔科夫估计值。(8 分)
四、(25 分)
考虑连续时间系统:
G(s)=s+aK(K=2,a=1)
使用双线性变换法(s=T2⋅z+1z−1),采样周期 T=1 秒。
1. 系统离散化(10分)
将 G(s) 转换为离散传递函数 G(z),推导其表达式并表示为:
G(z)=1+a1z−1b0+b1z−1
求 a1、b0、b1 的数值表达式。
2. 递推最小二乘参数估计(15 分)
实验测得输入输出数据如下:
| k | 0 | 1 | 2 | 3 |
|---|
| u(k) | 1.2 | 0.6 | 0.3 | 0.3 |
| y(k) | 0.0 | 1.2 | 0.9 | 0.7 |
假设模型结构为:
y(k)=−a1y(k−1)+b0u(k)+b1u(k−1)+e(k)
初始估计为 θ^0=[0,0,0]⊤,协方差矩阵 P0=1000I。
要求:
- 写出递推最小二乘的 Kk 、 θ^k 和 Pk 更新公式。(5 分)
- 递推计算 k=1,2,3 时的参数估计 θ^k。(7 分)
- 验证最终估计值是否与(1)中理论离散化参数是否一致,若不一致,可能原因有哪些。(3分)
五、(15 分)
离散系统动态方程为:
xk=[1011]xk−1+wk−1,wk∼N(0,Q),Q=[0.1000.1]
观测方程为:
yk=[10]xk+vk,vk∼N(0,R),R=0.5
已知初始状态 x^0∣0=[00],P0∣0=Ic
1.
计算一步预测 x^1∣0 和预测协方差 P1∣0;(5 分)
2.
计算卡尔曼增益 K1;(5 分)
3.
若 y1=2.0,更新状态估计 x^1∣1 和协方差 P1∣1。(5 分)