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
|
pbig(1) Scilab Function pbig(1)
NAME
pbig - eigen-projection
CALLING SEQUENCE
[Q,M]=pbig(A,thres,flag)
PARAMETERS
A : real square matrix
thres : real number
flag : character string ('c' or 'd')
Q,M : real matrices
DESCRIPTION
Projection on eigen-subspace associated with eigenvalues with real part >=
thres (flag='c') or with magnitude >= thres (flag='d').
The projection is defined by Q*M, Q is full column rank, M is full row rank
and M*Q=eye.
If flag='c', the eigenvalues of M*A*Q = eigenvalues of A with real part >=
thres.
If flag='d', the eigenvalues of M*A*Q = eigenvalues of A with magnitude >=
thres.
If flag='c' and if [Q1,M1] = full rank factorization (fullrf) of eye-Q*M
then eigenvalues of M1*A*Q1 = eigenvalues of A with real part < thres.
If flag='d' and if [Q1,M1] = full rank factorization (fullrf) of eye-Q*M
then eigenvalues of M1*A*Q1 = eigenvalues of A with magnitude < thres.
EXAMPLE
A=diag([1,2,3]);X=rand(A);A=inv(X)*A*X;
[Q,M]=pbig(A,1.5,'d');
spec(M*A*Q)
[Q1,M1]=fullrf(eye-Q*M);
spec(M1*A*Q1)
SEE ALSO
psmall, projspec, fullrf
AUTHOR
F. D. (1988)
|
