File: testnormopt

package info (click to toggle)
scilab 2.4-1
  • links: PTS
  • area: non-free
  • in suites: potato, slink
  • size: 55,196 kB
  • ctags: 38,019
  • sloc: ansic: 231,970; fortran: 148,976; tcl: 7,099; makefile: 4,585; sh: 2,978; csh: 154; cpp: 101; asm: 39; sed: 5
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)