13
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="equil1">
15
<refname>equil1</refname>
16
<refpurpose>行列の平衡 (非負定) の組</refpurpose>
20
<synopsis>[T [,siz]]=equil1(P,Q [,tol])</synopsis>
40
<para>3要素の整数ベクトル</para>
54
<literal>equil1</literal> は次のような <literal>t</literal> を計算します:
57
<literal>P1=T*P*T'</literal> および <literal>Q1=inv(T)'*Q*inv(T)</literal>
61
<literal>P1 = diag(S1,S2,0,0)</literal> および <literal>Q1 = diag(S1,0,S3,0)</literal>
63
<literal>S1,S2,S3</literal> は各次元が<literal>siz=[n1,n2,n3]</literal>の正定対角行列です.
66
<literal>tol</literal> はSVDのランク定義の閾値です.
71
<programlisting role="example"><![CDATA[
15
<refname>equil1</refname>
16
<refpurpose>行列の平衡 (非負定) の組</refpurpose>
20
<synopsis>[T [,siz]]=equil1(P,Q [,tol])</synopsis>
40
<para>3要素の整数ベクトル</para>
54
<literal>equil1</literal> は次のような <literal>t</literal> を計算します:
57
<literal>P1=T*P*T'</literal> および <literal>Q1=inv(T)'*Q*inv(T)</literal>
61
<literal>P1 = diag(S1,S2,0,0)</literal> および <literal>Q1 = diag(S1,0,S3,0)</literal>
63
<literal>S1,S2,S3</literal> は各次元が<literal>siz=[n1,n2,n3]</literal>の正定対角行列です.
66
<literal>tol</literal> はSVDのランク定義の閾値です.
71
<programlisting role="example"><![CDATA[
72
72
S1=rand(2,2);S1=S1*S1';
73
73
S2=rand(2,2);S2=S2*S2';
74
74
S3=rand(2,2);S3=S3*S3';
81
81
Q1=clean(inv(T)'*Q*inv(T))
82
82
]]></programlisting>
84
<refsection role="see also">
86
<simplelist type="inline">
88
<link linkend="balreal">balreal</link>
91
<link linkend="minreal">minreal</link>
94
<link linkend="equil">equil</link>
97
<link linkend="hankelsv">hankelsv</link>
84
<refsection role="see also">
86
<simplelist type="inline">
88
<link linkend="balreal">balreal</link>
91
<link linkend="minreal">minreal</link>
94
<link linkend="equil">equil</link>
97
<link linkend="hankelsv">hankelsv</link>