@@ -102,10 +102,10 @@ def ensemble_kalman_filter_predict(
102
102
Y[t] ~ observation_fn(X[t])
103
103
```
104
104
105
- Given the ensemble `state.particles` sampled from `P(X[t-1] | Y[t-1])`, this
106
- function produces the predicted (a.k.a. forecast or background) ensemble
107
- sampled from `P(X[t] | Y[t-1])`. This is the predicted next state *before*
108
- assimilating the observation `Y[t]`.
105
+ Given the ensemble `state.particles` sampled from `P(X[t-1] | Y[t-1], ... )`,
106
+ this function produces the predicted (a.k.a. forecast or background) ensemble
107
+ sampled from `P(X[t] | Y[t-1], ... )`. This is the predicted next state
108
+ *before* assimilating the observation `Y[t]`.
109
109
110
110
Typically, with `F` some deterministic mapping, `transition_fn(X)` returns a
111
111
normal distribution centered at `F(X)`.
@@ -177,9 +177,9 @@ def ensemble_kalman_filter_update(
177
177
Y[t] ~ observation_fn(X[t])
178
178
```
179
179
180
- Given the ensemble `state.particles` sampled from `P(X[t] | Y[t-1])`, this
181
- function assimilates obervation `Y[t]` to produce the updated ensemble sampled
182
- from `P(X[t] | Y[t])`.
180
+ Given the ensemble `state.particles` sampled from `P(X[t] | Y[t-1], ... )`,
181
+ this function assimilates obervation `Y[t]` to produce the updated ensemble
182
+ sampled from `P(X[t] | Y[t], ... )`.
183
183
184
184
Typically, with `G` some deterministic observation mapping,
185
185
`observation_fn(X)` returns a normal distribution centered at `G(X)`.
0 commit comments