连续函数离散化——零阶保持
公式:
1.
h
t
=
A
‾
h
t
−
1
+
B
‾
x
t
h_t = \overline{A}h_{t-1} + \overline{B}x_t
ht=Aht−1+Bxt
2.
h
0
=
B
‾
x
0
h_0 = \overline{B}x_0
h0=Bx0
3.
y
t
=
C
‾
h
t
y_t = \overline{C}h_t
yt=Cht
详细解释:
这些公式描述了状态空间模型(SSM)中的离散化过程,特别是使用零阶保持方法对连续时间系统进行离散化。以下是对每个公式的详细解释:
1. 离散化的状态更新方程
h t = A ‾ h t − 1 + B ‾ x t h_t = \overline{A}h_{t-1} + \overline{B}x_t ht=Aht−1+Bxt
解释:
-
h
t
h_t
ht:表示在离散时间
t
t
t时刻的隐状态向量。
-
h
t
−
1
h_{t-1}
ht−1:表示在上一离散时间
t
−
1
t-1
t−1时刻的隐状态向量。
-
A
‾
\overline{A}
A:离散化后的状态转移矩阵。
-
x
t
x_t
xt:表示在离散时间
t
t
t时刻的输入向量。
-
B
‾
\overline{B}
B:离散化后的输入矩阵。
零阶保持:
- 零阶保持假设在每个离散时间间隔内输入保持不变。通过这种方法,连续时间的状态空间方程可以被离散化,以便应用在数字计算和离散时间系统中。
这个公式表示当前时刻的状态 h t h_t ht是由上一时刻的状态 h t − 1 h_{t-1} ht−1和当前输入 x t x_t xt通过各自的离散化矩阵 A ‾ \overline{A} A和 B ‾ \overline{B} B共同作用得到的。
2. 初始条件
h 0 = B ‾ x 0 h_0 = \overline{B}x_0 h0=Bx0
解释:
-
h
0
h_0
h0:初始状态。
-
x
0
x_0
x0:初始输入。
-
B
‾
\overline{B}
B:离散化后的输入矩阵。
这个公式表示初始状态 h 0 h_0 h0是由初始输入 x 0 x_0 x0通过离散化后的输入矩阵 B ‾ \overline{B} B直接得到的。这提供了系统在开始时的状态条件。
3. 离散化的输出方程
y t = C ‾ h t y_t = \overline{C}h_t yt=Cht
解释:
-
y
t
y_t
yt:表示在离散时间
t
t
t时刻的输出向量。
-
C
‾
\overline{C}
C:离散化后的输出矩阵。
-
h
t
h_t
ht:当前时刻的隐状态向量。
这个公式表示当前时刻的输出 y t y_t yt是由当前时刻的状态 h t h_t ht通过离散化后的输出矩阵 C ‾ \overline{C} C得到的。即,输出是状态的线性变换。
离散化过程中的注意事项:
-
状态转移矩阵的离散化: A ‾ \overline{A} A通常是通过连续时间状态转移矩阵 A A A的指数矩阵计算得到:
A ‾ = e A Δ t \overline{A} = e^{A\Delta t} A=eAΔt
其中, Δ t \Delta t Δt是离散化的时间步长。 -
输入矩阵的离散化: B ‾ \overline{B} B可以通过如下公式得到:
B ‾ = ( ∫ 0 Δ t e A τ d τ ) B \overline{B} = \left( \int_0^{\Delta t} e^{A\tau} d\tau \right) B B=(∫0ΔteAτdτ)B
这里使用了连续时间输入矩阵 B B B和时间步长 Δ t \Delta t Δt。 -
输出矩阵: C ‾ \overline{C} C通常与连续时间系统的输出矩阵 C C C相同,因为输出是状态的直接观测,离散化不改变这个关系。
通过这些离散化方法,可以将连续时间的状态空间模型转化为离散时间模型,使其适用于数字控制和离散时间信号处理。这种转换允许我们在计算机上模拟和控制连续时间系统。