5
5
#include < vector>
6
6
7
7
int main ()
8
- {
9
-
10
- double Td = 0.11 ; // time delay in seconds.
11
- size_t order = 3 ; // order of the Pade approximation.
12
-
13
- auto tf_delay = ns_control_toolbox::pade (Td, order);
14
- tf_delay.print ();
15
-
16
- // Define Ts
17
- double Ts{ 0.1 };
18
- auto ss_sys = ns_control_toolbox::tf2ss (tf_delay);
19
- ss_sys.print ();
20
-
21
- // Test discretization and compare with Matlab.
22
-
23
- ss_sys.discretisize (Ts);
24
- ss_sys.print_discrete_system ();
8
+ {
9
+
10
+ double Td = 0.11 ; // time delay in seconds.
11
+ size_t order = 3 ; // order of the Pade approximation.
12
+
13
+ auto tf_delay = ns_control_toolbox::pade (Td, order);
14
+ tf_delay.print ();
15
+
16
+ // Define Ts
17
+ double Ts{ 0.1 };
18
+ auto ss_sys = ns_control_toolbox::tf2ss (tf_delay);
19
+ ss_sys.print ();
20
+
21
+ // Test discretization and compare with Matlab.
22
+
23
+ ss_sys.discretisize (Ts);
24
+ ss_sys.print_discrete_system ();
25
25
26
26
// ns_control_toolbox::balance(ss_sys.A_);
27
27
// ns_utils::print("Balanced matrix A_");
28
28
// ns_eigen_utils::printEigenMat(ss_sys.A_);
29
-
30
- // Concat
31
- auto AB = ns_eigen_utils::hstack<double >(ss_sys.A_ , ss_sys.B_ );
32
- auto CD = ns_eigen_utils::hstack<double >(ss_sys.C_ , ss_sys.D_ );
33
- auto ss_system = ns_eigen_utils::vstack<double >(AB, CD);
34
-
35
- ns_utils::print (" Matrix to be balanced " );
36
- ns_eigen_utils::printEigenMat (ss_system);
37
-
38
- ns_utils::print (" Balanced matrix for the system " );
39
- ns_control_toolbox::balance (ss_system);
40
- ns_eigen_utils::printEigenMat (ss_system);
41
-
42
- // PERMUTATION TEST
43
-
44
-
45
-
46
- return 0 ;
47
- }
29
+
30
+ // Concat
31
+ auto AB = ns_eigen_utils::hstack<double >(ss_sys.A () , ss_sys.B () );
32
+ auto CD = ns_eigen_utils::hstack<double >(ss_sys.C () , ss_sys.D () );
33
+ auto ss_system = ns_eigen_utils::vstack<double >(AB, CD);
34
+
35
+ ns_utils::print (" Matrix to be balanced " );
36
+ ns_eigen_utils::printEigenMat (ss_system);
37
+
38
+ ns_utils::print (" Balanced matrix for the system " );
39
+ ns_control_toolbox::balance (ss_system);
40
+ ns_eigen_utils::printEigenMat (ss_system);
41
+
42
+ // PERMUTATION TEST
43
+
44
+
45
+
46
+ return 0 ;
47
+ }
0 commit comments