1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
linfn Scilab Group Scilab Function linfn
NAME
linfn - infinity norm
CALLING SEQUENCE
[x,freq]=linfn(G,PREC,RELTOL,options);
PARAMETERS
G : is a syslin list
PREC : desired relative accuracy on the norm
RELTOL : relative threshold to decide when an eigenvalue can be
considered on the imaginary axis.
options : available options are 'trace' or 'cond'
x is the computed norm.
freq : vector
DESCRIPTION
Computes the Linf (or Hinf) norm of G This norm is well-defined as soon
as the realization G=(A,B,C,D) has no imaginary eigenvalue which is both
controllable and observable.
freq is a list of the frequencies for which ||G|| is attained,i.e., such
that ||G (j om)|| = ||G||.
If -1 is in the list, the norm is attained at infinity.
If -2 is in the list, G is all-pass in some direction so that ||G (j
omega)|| = ||G|| for all frequencies omega.
The algorithm follows the paper by G. Robel (AC-34 pp. 882-884, 1989).
The case D=0 is not treated separately due to superior accuracy of the
general method when (A,B,C) is nearly non minimal.
The 'trace' option traces each bisection step, i.e., displays the lower
and upper bounds and the current test point.
The 'cond' option estimates a confidence index on the computed value and
issues a warning if computations are ill-conditioned
In the general case (A neither stable nor anti-stable), no upper bound
is prespecified.
If by contrast A is stable or anti stable, lower and upper bounds are
computed using the associated Lyapunov solutions.
SEE ALSO
h_norm
AUTHOR
P. Gahinet
|
