1
<?xml version="1.0" encoding="UTF-8"?>
2
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns4="http://www.w3.org/1999/xhtml" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:id="cspect" xml:lang="ja">
4
<refname>cspect</refname>
13
[sm [,cwp]]=cspect(nlags,npoints,wtype,x [,y] [,wpar])
14
[sm [,cwp]]=cspect(nlags,npoints,wtype,nx [,ny] [,wpar])
23
<para>ベクトル, 最初の信号のデータ.</para>
30
ベクトル, 二番目の信号のデータ. <literal>y</literal>を
31
省略した場合,<literal>x</literal>に等しいと仮定されます(自己相関).
32
指定された場合,<literal>x</literal>と同じ数の要素を有している必要が
41
スカラー : <literal>x</literal>信号の点の数. この場合,
42
信号xのセグメントは<literal>getx</literal>という名前のユーザ定義の関数
51
<literal>y</literal>信号の点の数. この場合,
52
<literal>y</literal>のセグメントが <literal>gety</literal>
53
という名前のユーザ定義関数により読み込まれます (下記参照).
54
指定された場合,<literal>ny</literal>は<literal>nx</literal>に
62
<para>相関遅延の数s (正の整数)</para>
68
<para>変換する点の数 (正の整数)</para>
78
<literal>'re'</literal>: 矩形
83
<literal>'tr'</literal>: 三角形
88
<literal>'hm'</literal>: ハミング
93
<literal>'hn'</literal> : ハニング
98
<literal>'kr'</literal>: カイザー,この場合,wpar
104
<literal>'ch'</literal>: チェビシェフ, この場合, wpar
115
<literal>Kaiser および Chebyshev
122
'kr': <literal>wpar は厳密に正の数である必要があります</literal>
127
'ch': <literal>wpar</literal> は要素数2のベクトル
128
<literal> [main_lobe_width,side_lobe_height]である必要があります.
131
<literal>0<main_lobe_width<.5</literal>, および
132
<literal>side_lobe_height>0</literal>
142
正規化された周波数の間隔<literal>[0,1]</literal>において
144
大きさ<literal>npoints</literal>の行配列です.
154
チェビシェフウインドウの場合は,チェビシェフウインドウのパラメータ
156
指定されていないチェビシェフウインドウのパラメータ
165
2つの信号<literal>x</literal> および <literal>y</literal>の
166
信号の相互スペクトル,そうでない場合は <literal>x</literal>を計算します.
179
<mml:mi mml:fontstyle="italic">xy</mml:mi>
183
<mml:mo mml:stretchy="false">(</mml:mo>
184
<mml:mo mml:stretchy="false">ω</mml:mo>
185
<mml:mo mml:stretchy="false">)</mml:mo>
187
<mml:mo mml:stretchy="false">=</mml:mo>
194
<mml:mo mml:stretchy="false">(</mml:mo>
198
<mml:mo mml:stretchy="false">∑</mml:mo>
201
<mml:mo mml:stretchy="false">=</mml:mo>
206
<mml:mo mml:stretchy="false">−</mml:mo>
213
<mml:mo mml:stretchy="false">(</mml:mo>
215
<mml:mo mml:stretchy="false">)</mml:mo>
221
<mml:mo mml:stretchy="false">−</mml:mo>
224
<mml:mo mml:stretchy="false">ω</mml:mo>
229
<mml:mo mml:stretchy="false">)</mml:mo>
232
<mml:mo mml:stretchy="false">(</mml:mo>
236
<mml:mo mml:stretchy="false">∑</mml:mo>
239
<mml:mo mml:stretchy="false">=</mml:mo>
244
<mml:mo mml:stretchy="false">−</mml:mo>
248
<mml:mover mml:accent="true">
250
<mml:mo mml:stretchy="false">ˉ</mml:mo>
254
<mml:mo mml:stretchy="false">(</mml:mo>
256
<mml:mo mml:stretchy="false">)</mml:mo>
262
<mml:mo mml:stretchy="false">ω</mml:mo>
267
<mml:mo mml:stretchy="false">)</mml:mo>
270
<mml:annotation mml:encoding="StarMath 5.0"> S_{xy}(%omega)={1}
271
over {N} (sum from{n=0} to{N-1} x(n)e^{-i %omega n}) ( sum
272
from{n=0} to{N-1} bar y(n)e^{i %omega n})
280
<para>The correlation method calculates the spectral estimate as the
281
Fourier transform of a modified estimate of the auto/cross correlation
282
function. This auto/cross correlation modified estimate consist of
283
repeatedly calculating estimates of the autocorrelation function from
284
overlapping sub-segments if the data, and then averaging these estimates
285
to obtain the result.
287
<para>The number of points of the window is
288
<literal>2*nlags-1.</literal>
291
For batch processing, the<literal> x</literal> and
292
<literal>y</literal> data may be read segment by segment using the
293
<literal>getx </literal>and <literal>gety</literal> user defined
294
functions. These functions have the following calling sequence:
297
<literal>xk=getx(ns,offset)</literal> and
298
<literal>yk=gety(ns,offset)</literal> where <literal>ns</literal> is the
299
segment size and <literal>offset</literal> is the index of the first
300
element of the segment in the full signal.
305
<para>5.0.2までのバージョンのScilabでは,返り値は
311
<para>Oppenheim, A.V., and R.W. Schafer. Discrete-Time Signal Processing,
312
Upper Saddle River, NJ: Prentice-Hall, 1999
317
<programlisting role="example"><![CDATA[
318
rand('normal');rand('seed',0);
321
nf=33;bedge=[0 .1;.125 .5];des=[1 0];wate=[1 1];
322
h=eqfir(nf,bedge,des,wate);
323
//有色のデータを得るために白色のデータをフィルタ処理
324
h1=[h 0*ones(1:max(size(x))-1)];
325
x1=[x 0*ones(1:max(size(h))-1)];
326
hf=fft(h1,-1); xf=fft(x1,-1);yf=hf.*xf;y=real(fft(yf,1));
327
sm=cspect(100,200,'tr',y);
328
smsize=max(size(sm));fr=(1:smsize)/smsize;
332
<refsection role="see also">
334
<simplelist type="inline">
336
<link linkend="pspect">pspect</link>
339
<link linkend="mese">mese</link>
342
<link linkend="corr">corr</link>