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="glever"> inverse of matrix pencil</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>[Bfs,Bis,chis]=glever(E,A [,s]) </CALLING_SEQUENCE_ITEM>
15
<PARAM_NAME>E, A</PARAM_NAME>
17
<SP>: two real square matrices of same dimensions</SP>
21
<PARAM_NAME>s</PARAM_NAME>
23
<SP>: character string (default value '<VERB>s</VERB>')</SP>
27
<PARAM_NAME>Bfs,Bis</PARAM_NAME>
29
<SP>: two polynomial matrices</SP>
33
<PARAM_NAME>chis</PARAM_NAME>
46
by generalized Leverrier's algorithm for a matrix pencil.</P>
49
(s*E-A)^-1 = (Bfs/chis) - Bis.
52
<P><VERB>chis</VERB> = characteristic polynomial (up to a multiplicative constant).</P>
53
<P><VERB>Bfs</VERB> = numerator polynomial matrix.</P>
55
= polynomial matrix ( - expansion of <VERB>(s*E-A)^-1</VERB> at infinity).</P>
57
Note the - sign before <VERB>Bis</VERB>.</P>
59
<SECTION label="Caution">
61
This function uses <VERB>cleanp</VERB> to simplify <VERB>Bfs,Bis</VERB> and <VERB>chis</VERB>.</P>
65
s=%s;F=[-1,s,0,0;0,-1,0,0;0,0,s-2,0;0,0,0,s-1];
66
[Bfs,Bis,chis]=glever(F)
67
inv(F)-((Bfs/chis) - Bis)
90
<AUTHOR>F. D. (1988) </AUTHOR>