1
<?xml version="1.0" encoding="ISO-8859-1" standalone="no"?>
2
<!DOCTYPE MAN SYSTEM "../../manrev.dtd">
4
<LANGUAGE>eng</LANGUAGE>
6
<TYPE>Scilab Function</TYPE>
7
<DATE>April 1993</DATE>
8
<SHORT_DESCRIPTION name="randpencil"> random pencil</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>F=randpencil(eps,infi,fin,eta) </CALLING_SEQUENCE_ITEM>
15
<PARAM_NAME>eps</PARAM_NAME>
17
<SP>: vector of integers</SP>
21
<PARAM_NAME>infi</PARAM_NAME>
23
<SP>: vector of integers</SP>
27
<PARAM_NAME>fin</PARAM_NAME>
29
<SP>: real vector, or monic polynomial, or vector of monic polynomial</SP>
33
<PARAM_NAME>eta</PARAM_NAME>
35
<SP>: vector of integers</SP>
39
<PARAM_NAME>F</PARAM_NAME>
41
<SP>: real matrix pencil <VERB>F=s*E-A</VERB> (<VERB>s=poly(0,'s')</VERB>)</SP>
49
<VERB>F=randpencil(eps,infi,fin,eta)</VERB> returns a random pencil <VERB>F</VERB>
50
with given Kronecker structure. The structure is given by:
51
<VERB>eps=[eps1,...,epsk]</VERB>: structure of epsilon blocks (size eps1x(eps1+1),....)
52
<VERB>fin=[l1,...,ln]</VERB> set of finite eigenvalues (assumed real) (possibly [])
53
<VERB>infi=[k1,...,kp]</VERB> size of J-blocks at infinity
54
<VERB>ki>=1</VERB> (infi=[] if no J blocks).
55
<VERB>eta=[eta1,...,etap]</VERB>: structure ofeta blocks (size eta1+1)xeta1,...)</P>
56
<P><VERB>epsi</VERB>'s should be >=0, <VERB>etai</VERB>'s should be >=0, <VERB>infi</VERB>'s should
59
If <VERB>fin</VERB> is a (monic) polynomial, the finite block admits the roots of
60
<VERB>fin</VERB> as eigenvalues.</P>
62
If <VERB>fin</VERB> is a vector of polynomial, they are the finite elementary
63
divisors of <VERB>F</VERB> i.e. the roots of <VERB>p(i)</VERB> are finite
64
eigenvalues of <VERB>F</VERB>.</P>
68
F=randpencil([0,1],[2],[-1,0,1],[3]);
69
[Q,Z,Qd,Zd,numbeps,numbeta]=kroneck(F);
72
F=randpencil([],[1,2],s^3-2,[]); //regular pencil