|
30 | 30 |
|
31 | 31 | # visualize the results |
32 | 32 | plot_cmp([ri_without_bound, ri_with_bound], [0, 1], cs=["#FF5722", "#303F9F"]) |
33 | | - |
34 | | -# if __name__ == '__main__': |
35 | | -# # init dynamic system |
36 | | -# system = NonLinSys(Model(synchronous_machine, [2, 1])) |
37 | | -# |
38 | | -# # settings for the computation |
39 | | -# options = ASB2008CDC.Options() |
40 | | -# options.t_end = 3 |
41 | | -# options.step = 0.01 |
42 | | -# options.tensor_order = 3 |
43 | | -# options.taylor_terms = 4 |
44 | | -# |
45 | | -# options.u = Zonotope.zero(1, 1) |
46 | | -# options.u_trans = np.zeros(1) |
47 | | -# |
48 | | -# # settings for the using geometry |
49 | | -# Zonotope.REDUCE_METHOD = Zonotope.REDUCE_METHOD.GIRARD |
50 | | -# Zonotope.ORDER = 50 |
51 | | -# |
52 | | -# # z = Zonotope([0, 3], np.diag([0.2, 0.2])) |
53 | | -# z = Interval([-0.2, 2.8], [0.2, 3.2]) |
54 | | -# |
55 | | -# # reachable sets computation without boundary analysis |
56 | | -# options.r0 = [cvt2(z, Geometry.TYPE.ZONOTOPE)] |
57 | | -# ti_whole, tp_whole, _, _ = ASB2008CDC.reach(system, options) |
58 | | -# |
59 | | -# # reachable sets computation with boundary analysis |
60 | | -# options.r0 = boundary(z, 1, Geometry.TYPE.ZONOTOPE) |
61 | | -# ti_bound, tp_bound, _, _ = ASB2008CDC.reach(system, options) |
62 | | -# |
63 | | -# # visualize the results |
64 | | -# plot_cmp([tp_whole, tp_bound], [0, 1], cs=["#FF5722", "#303F9F"]) |
0 commit comments