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="time_id">
15
<refname>time_id</refname>
16
<refpurpose> SISO 最小二乗同定</refpurpose>
20
<synopsis>[H [,err]]=time_id(n,u,y)</synopsis>
39
<para>システムへの入力ベクトル</para>
43
<term>"impuls" </term>
45
<para>y はインパルス応答</para>
51
<para>yはステップ応答.</para>
67
y(1)==0 の場合は分母がn次,分子がn-1次の有理関数,
68
y(1)<>0 の場合, 分母と分子が n次の有理関数.
76
<literal>||y - impuls(H,npt)||^2</literal>, ただし,
77
<literal>impuls(H,npt)</literal> は
78
<literal>H</literal>のインパルス応答の最初の
79
<literal>npt</literal>個の係数です.
89
<literal>y</literal>が厳密にプロパーの場合(<literal>y(1)=0</literal>) ,
90
<literal>time_id</literal>は線形方程式に関する以下の最小二乗解
91
を計算します: <literal>coeff(Den,n):=1</literal> を拘束条件とする
92
<literal>Den*y-Num*u=0</literal>.
93
<literal>y(1)~=0</literal>の場合,
94
アルゴリズムはまず解のプロパーな部分を計算した後,y(1)を
100
<programlisting role="example"><![CDATA[
15
<refname>time_id</refname>
16
<refpurpose> SISO 最小二乗同定</refpurpose>
20
<synopsis>[H [,err]]=time_id(n,u,y)</synopsis>
39
<para>システムへの入力ベクトル</para>
43
<term>"impuls" </term>
45
<para>y はインパルス応答</para>
51
<para>yはステップ応答.</para>
67
y(1)==0 の場合は分母がn次,分子がn-1次の有理関数,
68
y(1)<>0 の場合, 分母と分子が n次の有理関数.
76
<literal>||y - impuls(H,npt)||^2</literal>, ただし,
77
<literal>impuls(H,npt)</literal> は
78
<literal>H</literal>のインパルス応答の最初の
79
<literal>npt</literal>個の係数です.
89
<literal>y</literal>が厳密にプロパーの場合(<literal>y(1)=0</literal>) ,
90
<literal>time_id</literal>は線形方程式に関する以下の最小二乗解
91
を計算します: <literal>coeff(Den,n):=1</literal> を拘束条件とする
92
<literal>Den*y-Num*u=0</literal>.
93
<literal>y(1)~=0</literal>の場合,
94
アルゴリズムはまず解のプロパーな部分を計算した後,y(1)を
100
<programlisting role="example"><![CDATA[
102
102
h=(1-2*z)/(z^2-0.5*z+5)
103
103
rep=[0;ldiv(h('num'),h('den'),20)]; //インパルス応答
112
112
H=time_id(2,'step',rep)
113
113
H=time_id(3,u,rep) // uを入力として指定, 高次過ぎる系を指定
114
114
]]></programlisting>
116
<refsection role="see also">
118
<simplelist type="inline">
120
<link linkend="imrep2ss">imrep2ss</link>
123
<link linkend="arl2">arl2</link>
126
<link linkend="armax">armax</link>
129
<link linkend="frep2tf">frep2tf</link>
116
<refsection role="see also">
118
<simplelist type="inline">
120
<link linkend="imrep2ss">imrep2ss</link>
123
<link linkend="arl2">arl2</link>
126
<link linkend="armax">armax</link>
129
<link linkend="frep2tf">frep2tf</link>