File: testnormopt

package info (click to toggle)
scilab 2.6-4
  • links: PTS
  • area: non-free
  • in suites: woody
  • size: 54,632 kB
  • ctags: 40,267
  • sloc: ansic: 267,851; fortran: 166,549; sh: 10,005; makefile: 4,119; tcl: 1,070; cpp: 233; csh: 143; asm: 135; perl: 130; java: 39
file content (22 lines) | stat: -rw-r--r-- 476 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
// Copyright INRIA
getf(SCI+'/demos/lmitool/normopt.sci','c')
n=5;nu=2;ny=2;
//rand('seed',0)
A=rand(n,n);
B=rand(n,nu);
C=rand(ny,n);
D=rand(nu,ny);
//find H_inf norm of linear system A,B,C,D
//   min gama 
//   such that
//   X-X'=0
//  -[A*X+X*A',B,X*C';B',-gama*Ib,D';C*X,D,-gama*Ic]>0
[X,gamaopt]=normopt(A,B,C,D)
h_norm(syslin('c',A,B,C,D))


getf(SCI+'/demos/lmitool/normoptd.sci','c')
[X,g2opt]=normoptd(A,B,C,D);
sqrt(g2opt)
dhnorm(syslin('d',A,B,C,D),0.00001,1000)