@@ -344,21 +344,13 @@ each control period ``k``, see [`initpred!`](@ref) and [`linconstraint!`](@ref).
344344 ```math
345345 \b egin{aligned}
346346 \m athbf{e_x̂} &= \b egin{bmatrix}
347- \m athbf{W}(H_p-1) \m athbf{B̂_u} &
348- \m athbf{W}(H_p-n_1-1) \m athbf{B̂_u} &
349- \m athbf{W}(H_p-n_2-1) \m athbf{B̂_u} &
350- \c dots &
351- \m athbf{W}(H_p-n_{H_c-1}-1) \m athbf{B̂_u} \e nd{bmatrix} \\
347+ \m athbf{W}(H_p-1)\m athbf{B̂_u} & \m athbf{W}(H_p-n_1-1)\m athbf{B̂_u} & \c dots & \m athbf{W}(H_p-n_{H_c-1}-1)\m athbf{B̂_u} \e nd{bmatrix} \\
352348 \m athbf{g_x̂} &= \m athbf{Â}^{H_p-1} \m athbf{B̂_d} \\
353349 \m athbf{j_x̂} &= \b egin{bmatrix}
354- \m athbf{Â}^{H_p-2} \m athbf{B̂_d} &
355- \m athbf{Â}^{H_p-3} \m athbf{B̂_d} &
356- \c dots &
357- \m athbf{0}
358- \e nd{bmatrix} \\
350+ \m athbf{Â}^{H_p-2}\m athbf{B̂_d} & \m athbf{Â}^{H_p-3}\m athbf{B̂_d} & \c dots & \m athbf{0} \e nd{bmatrix} \\
359351 \m athbf{k_x̂} &= \m athbf{Â}^{H_p} \\
360352 \m athbf{v_x̂} &= \m athbf{W}(H_p-1)\m athbf{B̂_u} \\
361- \m athbf{b_x̂} &= \m athbf{W}(H_p-1) \m athbf{\b ig(f̂_{op} - x̂_{op}\b ig)}
353+ \m athbf{b_x̂} &= \m athbf{W}(H_p-1)\m athbf{\b ig(f̂_{op} - x̂_{op}\b ig)}
362354 \e nd{aligned}
363355 ```
364356"""
@@ -392,14 +384,11 @@ function init_predmat(
392384 for j= 1 : Hp
393385 iRow = (1 : ny) .+ ny* (j- 1 )
394386 K[iRow,:] = Ĉ* getpower (Âpow, j)
395- end
387+ end
396388 # --- previous manipulated inputs lastu0 ---
397389 vx̂ = W (Hp- 1 )* B̂u
398390 V = Matrix {NT} (undef, Hp* ny, nu)
399- for j= 1 : Hp
400- iRow = (1 : ny) .+ ny* (j- 1 )
401- V[iRow,:] = Ĉ* W (j- 1 )* B̂u
402- end
391+ Q! (V, 0 , Hp)
403392 # --- decision variables Z ---
404393 nZ = get_nZ (estim, transcription, Hp, Hc)
405394 ex̂ = Matrix {NT} (undef, nx̂, nZ)
0 commit comments