<< slicot_tb01id | SLICOT |
[A_OUT, E_OUT, B_OUT, C_OUT, LSCALE, RSCALE, INFO] = slicot_tg01ad(JOB, THRESH, A_IN, E_IN, B_IN, C_IN) |
= 'A': All matrices are involved in balancing; = 'B': B, A and E matrices are involved in balancing; = 'C': C, A and E matrices are involved in balancing; = 'N': B and C matrices are not involved in balancing.
Threshold value for magnitude of elements: elements with magnitude less than or equal to THRESH are ignored for balancing.
The leading L-by-N part of this array must contain the state dynamics matrix A.
The leading L-by-N part of this array must contain the descriptor matrix E.
The leading L-by-M part of this array must contain the input/state matrix B.
The leading P-by-N part of this array must contain the state/output matrix C.
The leading L-by-N part of this array contains the balanced matrix Dl*A*Dr.
The leading L-by-N part of this array contains the balanced matrix Dl*E*Dr.
The leading L-by-M part of this array contains the balanced matrix Dl*B.
The leading P-by-N part of this array contains the balanced matrix C*Dr.
The scaling factors applied to S from left.
The scaling factors applied to S from right.
= 0: successful exit.
To balance the matrices of the system pencil corresponding to the descriptor triple (A-lambda E,B,C), by balancing.
L = 4;
N = 4;
M = 2;
P = 2;
JOB = 'A';
THRESH = 0;
A_IN = [ -1 0 0 0.003;
0 0 0.1000 0.02;
100 10 0 0.4;
0 0 0 0.0];
E_IN = [1 0.2 0 0.0;
0 1 0 0.01;
300 90 6 0.3;
0 0 20 0.0];
B_IN = [10 0;
0 0;
0 1000;
10000 10000];
C_IN = [-0.1 0.0 0.001 0.0;
0.0 0.01 -0.001 0.0001];
[A_OUT, E_OUT, B_OUT, C_OUT, LSCALE, RSCALE, INFO] = slicot_tg01ad(JOB, THRESH, A_IN, E_IN, B_IN, C_IN)
Version | Description |
---|---|
1.0.0 | initial version |
SLICOT Documentation