|
599 | 599 | @test mhe8.nint_u == [1, 1] |
600 | 600 | @test mhe8.nint_ym == [0, 0] |
601 | 601 |
|
602 | | - mhe9 = MovingHorizonEstimator(nonlinmodel, 5, 1:2, 0, [1, 1], I(6), I(6), I(2), Model(Ipopt.Optimizer)) |
| 602 | + I_6 = Matrix{Float64}(I, 6, 6) |
| 603 | + I_2 = Matrix{Float64}(I, 2, 2) |
| 604 | + mhe9 = MovingHorizonEstimator(nonlinmodel, 5, 1:2, 0, [1, 1], I_6, I_6, I_2, Model(Ipopt.Optimizer)) |
603 | 605 | @test mhe9.P̂0 ≈ I(6) |
604 | 606 | @test mhe9.Q̂ ≈ I(6) |
605 | 607 | @test mhe9.R̂ ≈ I(2) |
|
619 | 621 | linmodel1 = LinModel(sys,Ts,i_u=[1,2]) |
620 | 622 | f(x,u,_) = linmodel1.A*x + linmodel1.Bu*u |
621 | 623 | h(x,_) = linmodel1.C*x |
622 | | - nonlinmodel = setop!(NonLinModel(f, h, Ts, 2, 2, 2), uop=[10,50], yop=[50,30]) |
| 624 | + nonlinmodel = setop!(NonLinModel(f, h, 1000*Ts, 2, 2, 2), uop=[10,50], yop=[50,30]) |
623 | 625 | mhe1 = MovingHorizonEstimator(nonlinmodel, He=2) |
624 | 626 | @test updatestate!(mhe1, [10, 50], [50, 30]) ≈ zeros(4) atol=1e-9 |
625 | 627 | @test updatestate!(mhe1, [10, 50], [50, 30], Float64[]) ≈ zeros(4) atol=1e-9 |
|
0 commit comments