1
<?xml version="1.0" encoding="UTF-8"?>
3
* Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
4
* Copyright (C) 2008 - INRIA
6
* This file must be used under the terms of the CeCILL.
7
* This source file is licensed as described in the file COPYING, which
8
* you should have received as part of this distribution. The terms
9
* are also available at
10
* http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
13
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:lang="ja" xml:id="kroneck">
15
<refname>kroneck</refname>
16
<refpurpose>行列ペンシルのクロネッカー形式</refpurpose>
20
<synopsis>[Q,Z,Qd,Zd,numbeps,numbeta]=kroneck(F)
21
[Q,Z,Qd,Zd,numbeps,numbeta]=kroneck(E,A)
31
実数行列ペンシル <literal>F=s*E-A</literal>
38
<para>同じ次元の実数行列</para>
54
<term>numbeps,numeta</term>
64
行列ペンシルのクロネッカー形式: <literal>kroneck</literal> は,
65
ペンシル<literal>F=s*E -A</literal>を以下のような上三角形式に変換する
66
2つの直交行列<literal>Q, Z</literal>を計算します:
68
<programlisting role=""><![CDATA[
69
| sE(eps)-A(eps) | X | X | X |
70
|----------------|----------------|------------|---------------|
71
| O | sE(inf)-A(inf) | X | X |
72
Q(sE-A)Z = |---------------------------------|----------------------------|
74
| 0 | 0 | sE(f)-A(f) | X |
75
|--------------------------------------------------------------|
77
| 0 | 0 | 0 | sE(eta)-A(eta)|
80
4個のブロックの次元は以下のように指定されます:
83
<literal>eps=Qd(1) x Zd(1)</literal>, <literal>inf=Qd(2) x Zd(2)</literal>,
84
<literal>f = Qd(3) x Zd(3)</literal>, <literal>eta=Qd(4)xZd(4)</literal>
87
<literal>inf</literal>ブロックにはペンシルの無限大モードが含まれます.
90
<literal>f</literal> ブロックにはペンシルの有限モードが含まれます.
93
イプシロンとetaブロックの構造は以下のように指定されます:
96
<literal>numbeps(1)</literal> = 大きさ 0 x 1のepsブロックの番号
99
<literal>numbeps(2)</literal> = 大きさ 1 x 2のepsブロックの番号
102
<literal>numbeps(3)</literal> = 大きさ 2 x 3のepsブロックの番号 etc...
105
<literal>numbeta(1)</literal> = 大きさ 1 x 0のetaブロックの番号
108
<literal>numbeta(2)</literal> = 大きさ 2 x 1のetaブロックの番号
111
<literal>numbeta(3)</literal> = 大きさ 3 x 2のetaブロックの番号 etc...
114
このコードはT. Beelen (Slicot-WGS group)によるものです.
119
<programlisting role="example"><![CDATA[
120
F=randpencil([1,1,2],[2,3],[-1,3,1],[0,3]);
121
Q=rand(17,17);Z=rand(18,18);F=Q*F*Z;
122
//random pencil with eps1=1,eps2=1,eps3=1; 2 J-blocks @ infty
123
//with dimensions 2 and 3
124
//3 finite eigenvalues at -1,3,1 and eta1=0,eta2=3
125
[Q,Z,Qd,Zd,numbeps,numbeta]=kroneck(F);
126
[Qd(1),Zd(1)] //eps. part is sum(epsi) x (sum(epsi) + number of epsi)
127
[Qd(2),Zd(2)] //infinity part
128
[Qd(3),Zd(3)] //finite part
129
[Qd(4),Zd(4)] //eta part is (sum(etai) + number(eta1)) x sum(etai)
134
<refsection role="see also">
136
<simplelist type="inline">
138
<link linkend="gschur">gschur</link>
141
<link linkend="gspec">gspec</link>
144
<link linkend="systmat">systmat</link>
147
<link linkend="pencan">pencan</link>
150
<link linkend="randpencil">randpencil</link>
153
<link linkend="trzeros">trzeros</link>