← Back to branch summary

~ubuntu-branches/ubuntu/hoary/scilab/hoary

« back to all changes in this revision

Viewing changes to man/control/projsl.cat

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2005-01-09 22:58:21 UTC
  • mfrom: (1.1.1 upstream)
  • Revision ID: james.westby@ubuntu.com-20050109225821-473xr8vhgugxxx5j
Tags: 3.0-12
http://bugs.debian.org/255869
changed configure.in to build scilab's own malloc.o, closes: #255869

Show diffs side-by-side

added added

removed removed

Lines of Context:
man/control/projsl.cat
1
 
projsl            Scilab Group            Scilab Function            projsl
2
 
NAME
3
 
   projsl - linear system projection
4
 
  
5
 
CALLING SEQUENCE
6
 
 [slp]=projsl(sl,Q,M)
7
 
PARAMETERS
8
 
 sl,slp          : syslin lists
9
 
                 
10
 
 Q,M             : matrices (projection factorization)
11
 
                 
12
 
DESCRIPTION
13
 
   slp= projected model of sl where Q*M is the full rank  factorization of
14
 
  the projection.
15
 
  
16
 
   If (A,B,C,D) is the representation of sl, the projected model is given by
17
 
  (M*A*Q,M*B,C*Q,D).
18
 
  
19
 
   Usually, the projection Q*M is obtained as the spectral projection of an
20
 
  appropriate auxiliary matrix W e.g. W = product of (weighted) gramians or
21
 
  product of Riccati equations.
22
 
  
23
 
EXAMPLE
24
 
 rand('seed',0);sl=ssrand(2,2,5);[A,B,C,D]=abcd(sl);poles=spec(A)
25
 
 [Q,M]=pbig(A,0,'c');  //keeping unstable poles
26
 
 slred=projsl(sl,Q,M);spec(slred('A'))
27
 
 sl('D')=rand(2,2);  //making proper system
28
 
 trzeros(sl)  //zeros of sl
29
 
 wi=inv(sl);  //wi=inverse in state-space
30
 
 [q,m]=psmall(wi('A'),2,'d');  //keeping small zeros (poles of wi) i.e. abs(z)<2
31
 
 slred2=projsl(sl,q,m);
32
 
 trzeros(slred2)  //zeros of slred2 = small zeros of sl
33
 
 //  Example keeping second order modes
34
 
 A=diag([-1,-2,-3]);
35
 
 sl=syslin('c',A,rand(3,2),rand(2,3));[nk2,W]=hankelsv(sl)
36
 
 [Q,M]=pbig(W,nk2(2)-%eps,'c');    //keeping 2 eigenvalues of W
37
 
 slr=projsl(sl,Q,M);  //reduced model
38
 
 hankelsv(slr)
39
 
SEE ALSO
40
 
   pbig
41
 
  
42
 
AUTHOR
43
 
   F. D. 
44