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 * Copyright (C) 2012 - 2016 - Scilab Enterprises
8 * This file is hereby licensed under the terms of the GNU GPL v2.0,
9 * pursuant to article 5.3.4 of the CeCILL v.2.1.
10 * This file was originally licensed under the terms of the CeCILL v2.1,
11 * and continues to be available under such terms.
12 * For more information, see the COPYING file which you should have received
13 * along with this program.
16 <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" xmlns:scilab="http://www.scilab.org" xml:lang="ja" xml:id="quaskro">
18 <refname>quaskro</refname>
19 <refpurpose>準クロネッカー形式</refpurpose>
24 [Q, Z, Qd, Zd, numbeps, numbeta] = quaskro(F)
25 [Q, Z, Qd, Zd, numbeps, numbeta] = quaskro(E,A)
26 [Q, Z, Qd, Zd, numbeps, numbeta] = quaskro(F,tol)
27 [Q, Z, Qd, Zd, numbeps, numbeta] = quaskro(E,A,tol)
30 <refsection role="parameters">
37 実数行列ペンシル <literal>F=s*E-A</literal> (<literal>s=poly(0,'s')</literal>)
44 <para>同じ次元の実数行列</para>
50 <para>実数 (許容誤差,デフォルト値=1.d-10)</para>
73 <refsection role="description">
77 <literal>quaskro</literal>は,ペンシル<literal>F=s*E -A</literal>を上三角行列形式に変換する
78 直交行列 <literal>Q, Z</literal>を計算します:
81 | sE(eps)-A(eps) | X | X |
82 |----------------|----------------|------------|
83 | O | sE(inf)-A(inf) | X |
84 Q(sE-A)Z = |=================================|============|
92 <literal>eps=Qd(1) x Zd(1)</literal>, <literal>inf=Qd(2) x Zd(2)</literal>,
93 <literal>r = Qd(3) x Zd(3)</literal>
96 <literal>inf</literal> ブロックには, ペンシルの無限大モードが含まれます.
99 <literal>f</literal> ブロックには, ペンシルの有限モードが含まれます.
102 epsilonブロックの構造は次のように指定されます:
105 <literal>numbeps(1)</literal> = 大きさ 0 x 1のepsブロックの数
108 <literal>numbeps(2)</literal> = 大きさ 1 x 2のepsブロックの数
111 <literal>numbeps(3)</literal> = 大きさ 2 x 3のepsブロックの数 etc...
114 完全な(4ブロックの)クロネッカー形式は,
115 (pertransposed)ペンシル<literal>sE(r)-A(r)</literal>を指定して
116 <literal>quaskro</literal>をコールする
117 関数<literal>kroneck</literal>により指定されます.
120 このコード T. Beelenによるものです.
123 <refsection role="see also">
125 <simplelist type="inline">
127 <link linkend="kroneck">kroneck</link>
130 <link linkend="schur">schur</link>
133 <link linkend="spec">spec</link>