~ubuntu-branches/ubuntu/raring/scilab/raring-proposed

« back to all changes in this revision

Viewing changes to modules/signal_processing/help/ja_JP/pspect.xml

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2012-08-30 14:42:38 UTC
  • mfrom: (1.4.7)
  • Revision ID: package-import@ubuntu.com-20120830144238-c1y2og7dbm7m9nig
Tags: 5.4.0-beta-3-1~exp1
* New upstream release
* Update the scirenderer dep
* Get ride of libjhdf5-java dependency

Show diffs side-by-side

added added

removed removed

Lines of Context:
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="pspect" xml:lang="ja">
3
 
  <refnamediv>
4
 
    <refname>pspect</refname>
5
 
    <refpurpose>
6
 
      Welchの平均ペリオドグラム法による
7
 
      2つの離散時間信号の間の両面相互スペクトル推定.
8
 
    </refpurpose>
9
 
  </refnamediv>
10
 
  <refsynopsisdiv>
11
 
    <title>呼び出し手順</title>
12
 
    <synopsis>
13
 
      [sm [,cwp]]=pspect(sec_step,sec_leng,wtype,x [,y] [,wpar])
14
 
      [sm [,cwp]]=pspect(sec_step,sec_leng,wtype,nx [,ny] [,wpar])
15
 
    </synopsis>
16
 
  </refsynopsisdiv>
17
 
  <refsection>
18
 
    <title>パラメータ</title>
19
 
    <variablelist>
20
 
      <varlistentry>
21
 
        <term>x</term>
22
 
        <listitem>
23
 
          <para>ベクトル, 最初の信号の時間領域標本.</para>
24
 
        </listitem>
25
 
      </varlistentry>
26
 
      <varlistentry>
27
 
        <term>y</term>
28
 
        <listitem>
29
 
          <para>
30
 
            ベクトル, 2番目の信号の時間領域標本. <literal>y</literal>が
31
 
            省略された場合,<literal>x</literal>に等しいと仮定されます(自己相関).
32
 
            指定された場合, <literal>x</literal>と同じ数の要素を有している必要があります.
33
 
          </para>
34
 
        </listitem>
35
 
      </varlistentry>
36
 
      <varlistentry>
37
 
        <term>nx</term>
38
 
        <listitem>
39
 
          <para>
40
 
            スカラー : <literal>x</literal>信号の標本数.
41
 
            この場合,<literal>x</literal>信号のセグメントは<literal>getx</literal>という名前の
42
 
            ユーザ定義の関数によりロードされます (下記参照).
43
 
          </para>
44
 
        </listitem>
45
 
      </varlistentry>
46
 
      <varlistentry>
47
 
        <term>ny</term>
48
 
        <listitem>
49
 
          <para>
50
 
            スカラー : <literal>y</literal>信号の標本数.
51
 
            この場合,<literal>x</literal>信号のセグメントは<literal>gety</literal>という名前の
52
 
            ユーザ定義の関数によりロードされます (下記参照).
53
 
            指定された場合, <literal>ny</literal>は<literal>nx</literal>と等しい必要があります.
54
 
          </para>
55
 
        </listitem>
56
 
      </varlistentry>
57
 
      <varlistentry>
58
 
        <term>sec_step</term>
59
 
        <listitem>
60
 
          <para>
61
 
            各データウインドウのオフセット.
62
 
            オーバーラップ <literal>D</literal> はi <literal>sec_leng -sec_step</literal>により指定されます.
63
 
            <literal>sec_step==sec_leng/2</literal>の場合, 50% オーバーラップとなります. 
64
 
          </para>
65
 
        </listitem>
66
 
      </varlistentry>
67
 
      <varlistentry>
68
 
        <term>sec_leng</term>
69
 
        <listitem>
70
 
          <para>ウインドウの点の数.</para>
71
 
        </listitem>
72
 
      </varlistentry>
73
 
      <varlistentry>
74
 
        <term>wtype</term>
75
 
        <listitem>
76
 
          <para>ウインドウの種類</para>
77
 
          <itemizedlist>
78
 
            <listitem>
79
 
              <para>
80
 
                <literal>'re'</literal>: 矩形
81
 
              </para>
82
 
            </listitem>
83
 
            <listitem>
84
 
              <para>
85
 
                <literal>'tr'</literal>: 三角形
86
 
              </para>
87
 
            </listitem>
88
 
            <listitem>
89
 
              <para>
90
 
                <literal>'hm'</literal>: ハミング
91
 
              </para>
92
 
            </listitem>
93
 
            <listitem>
94
 
              <para>
95
 
                <literal>'hn'</literal> : ハニング
96
 
              </para>
97
 
            </listitem>
98
 
            <listitem>
99
 
              <para>
100
 
                <literal>'kr'</literal>: カイザー,この場合,wpar引数を指定する必要があります
101
 
              </para>
102
 
            </listitem>
103
 
            <listitem>
104
 
              <para>
105
 
                <literal>'ch'</literal>: チェビシェフ, この場合,wpar引数を指定する必要があります
106
 
              </para>
107
 
            </listitem>
108
 
          </itemizedlist>
109
 
        </listitem>
110
 
      </varlistentry>
111
 
      <varlistentry>
112
 
        <term>wpar</term>
113
 
        <listitem>
114
 
          <para>
115
 
            <literal>カイザーおよびチェビシェフウインドウ</literal>に関する
116
 
            オプションのパラメータ:
117
 
          </para>
118
 
          <itemizedlist>
119
 
            <listitem>
120
 
              <para>
121
 
                'kr': <literal>wpar</literal>は厳密に正の数とする必要があります
122
 
              </para>
123
 
            </listitem>
124
 
            <listitem>
125
 
              <para>
126
 
                'ch': <literal>wpar</literal> は要素数2のベクトル
127
 
                <literal> [main_lobe_width,side_lobe_height]</literal>
128
 
                とする必要があります.
129
 
                ただし,<literal>0&lt;main_lobe_width&lt;.5</literal>, および
130
 
                <literal>side_lobe_height&gt;0</literal>です.
131
 
              </para>
132
 
            </listitem>
133
 
          </itemizedlist>
134
 
        </listitem>
135
 
      </varlistentry>
136
 
      <varlistentry>
137
 
        <term>sm</term>
138
 
        <listitem>
139
 
          <para>
140
 
            正規化された周波数の範囲 <literal>[0,1]</literal>における両面スペクトル推定.
141
 
            これは,<literal>sec_len</literal>個の要素を有する行配列です.
142
 
            この配列は自己相関の場合には実数,相互相関の場合には複素数となります.
143
 
          </para>
144
 
          <para>
145
 
            関係する正規化された周波数配列は
146
 
            <literal>linspace(0,1,sec_len)</literal>です.
147
 
          </para>
148
 
        </listitem>
149
 
      </varlistentry>
150
 
      <varlistentry>
151
 
        <term>cwp</term>
152
 
        <listitem>
153
 
          <para>
154
 
            チェビシェフウインドの場合は
155
 
            チェビシェフウインドウの未指定のパラメータ,
156
 
            または空の行列.
157
 
          </para>
158
 
        </listitem>
159
 
      </varlistentry>
160
 
    </variablelist>
161
 
  </refsection>
162
 
  <refsection>
163
 
    <title>説明</title>
164
 
    <para>
165
 
      2つの信号<literal>x</literal> および <literal>y</literal>が両方共指定
166
 
      された場合,相互スペクトル推定値,
167
 
      そうでない場合に自己相関スペクトルを計算します.
168
 
      修正ペリオドグラム法により得られたスペクトル推定値.
169
 
    </para>
170
 
    <para>
171
 
      2つの信号<literal>x</literal> および  <literal>y</literal>の
172
 
      相互スペクトルは以下のように定義されます
173
 
    </para>
174
 
    <para>
175
 
      <inlinemediaobject>
176
 
        <imageobject>
177
 
          <imagedata>
178
 
            <mml:math>
179
 
              <mml:semantics>
180
 
                <mml:mrow>
181
 
                  <mml:msub>
182
 
                    <mml:mi>S</mml:mi>
183
 
                    <mml:mi mml:fontstyle="italic">xy</mml:mi>
184
 
                  </mml:msub>
185
 
                  <mml:mrow>
186
 
                    <mml:mrow>
187
 
                      <mml:mo mml:stretchy="false">(</mml:mo>
188
 
                      <mml:mo mml:stretchy="false">ω</mml:mo>
189
 
                      <mml:mo mml:stretchy="false">)</mml:mo>
190
 
                    </mml:mrow>
191
 
                    <mml:mo mml:stretchy="false">=</mml:mo>
192
 
                    <mml:mfrac>
193
 
                      <mml:mn>1</mml:mn>
194
 
                      <mml:mi>N</mml:mi>
195
 
                    </mml:mfrac>
196
 
                  </mml:mrow>
197
 
                  <mml:mrow>
198
 
                    <mml:mo mml:stretchy="false">(</mml:mo>
199
 
                    <mml:mrow>
200
 
                      <mml:mrow>
201
 
                        <mml:munderover>
202
 
                          <mml:mo mml:stretchy="false">∑</mml:mo>
203
 
                          <mml:mrow>
204
 
                            <mml:mi>n</mml:mi>
205
 
                            <mml:mo mml:stretchy="false">=</mml:mo>
206
 
                            <mml:mn>0</mml:mn>
207
 
                          </mml:mrow>
208
 
                          <mml:mrow>
209
 
                            <mml:mi>N</mml:mi>
210
 
                            <mml:mo mml:stretchy="false">−</mml:mo>
211
 
                            <mml:mn>1</mml:mn>
212
 
                          </mml:mrow>
213
 
                        </mml:munderover>
214
 
                        <mml:mi>x</mml:mi>
215
 
                      </mml:mrow>
216
 
                      <mml:mrow>
217
 
                        <mml:mo mml:stretchy="false">(</mml:mo>
218
 
                        <mml:mi>n</mml:mi>
219
 
                        <mml:mo mml:stretchy="false">)</mml:mo>
220
 
                      </mml:mrow>
221
 
                      <mml:msup>
222
 
                        <mml:mi>e</mml:mi>
223
 
                        <mml:mrow>
224
 
                          <mml:mrow>
225
 
                            <mml:mo mml:stretchy="false">−</mml:mo>
226
 
                            <mml:mi>i</mml:mi>
227
 
                          </mml:mrow>
228
 
                          <mml:mo mml:stretchy="false">ω</mml:mo>
229
 
                          <mml:mi>n</mml:mi>
230
 
                        </mml:mrow>
231
 
                      </mml:msup>
232
 
                    </mml:mrow>
233
 
                    <mml:mo mml:stretchy="false">)</mml:mo>
234
 
                  </mml:mrow>
235
 
                  <mml:mrow>
236
 
                    <mml:mo mml:stretchy="false">(</mml:mo>
237
 
                    <mml:mrow>
238
 
                      <mml:mrow>
239
 
                        <mml:munderover>
240
 
                          <mml:mo mml:stretchy="false">∑</mml:mo>
241
 
                          <mml:mrow>
242
 
                            <mml:mi>n</mml:mi>
243
 
                            <mml:mo mml:stretchy="false">=</mml:mo>
244
 
                            <mml:mn>0</mml:mn>
245
 
                          </mml:mrow>
246
 
                          <mml:mrow>
247
 
                            <mml:mi>N</mml:mi>
248
 
                            <mml:mo mml:stretchy="false">−</mml:mo>
249
 
                            <mml:mn>1</mml:mn>
250
 
                          </mml:mrow>
251
 
                        </mml:munderover>
252
 
                        <mml:mover mml:accent="true">
253
 
                          <mml:mi>y</mml:mi>
254
 
                          <mml:mo mml:stretchy="false">ˉ</mml:mo>
255
 
                        </mml:mover>
256
 
                      </mml:mrow>
257
 
                      <mml:mrow>
258
 
                        <mml:mo mml:stretchy="false">(</mml:mo>
259
 
                        <mml:mi>n</mml:mi>
260
 
                        <mml:mo mml:stretchy="false">)</mml:mo>
261
 
                      </mml:mrow>
262
 
                      <mml:msup>
263
 
                        <mml:mi>e</mml:mi>
264
 
                        <mml:mrow>
265
 
                          <mml:mi>i</mml:mi>
266
 
                          <mml:mo mml:stretchy="false">ω</mml:mo>
267
 
                          <mml:mi>n</mml:mi>
268
 
                        </mml:mrow>
269
 
                      </mml:msup>
270
 
                    </mml:mrow>
271
 
                    <mml:mo mml:stretchy="false">)</mml:mo>
272
 
                  </mml:mrow>
273
 
                </mml:mrow>
274
 
                <mml:annotation mml:encoding="StarMath 5.0"> S_{xy}(%omega)={1}
275
 
                  over {N} (sum from{n=0} to{N-1} x(n)e^{-i %omega n}) ( sum
276
 
                  from{n=0} to{N-1} bar y(n)e^{i %omega n})
277
 
                </mml:annotation>
278
 
              </mml:semantics>
279
 
            </mml:math>
280
 
          </imagedata>
281
 
        </imageobject>
282
 
      </inlinemediaobject>
283
 
    </para>
284
 
    <para>
285
 
      スペクトル推定の修正ペリオドグラム法は,
286
 
      <literal>x</literal> および <literal>y</literal>に含まれる
287
 
      データのウインドウ処理されたサブセクションのペリオドグラムを
288
 
      繰り返し計算します.
289
 
      これらのペリオドグラムは,最終的なスペクトル推定値を得るために平均化され,
290
 
      適当な定数により正規化されます.
291
 
      これは推定値の分散を減らすための平均化処理です.
292
 
    </para>
293
 
    <para>
294
 
      バッチ処理の場合, <literal>x</literal> および
295
 
      <literal>y</literal> データはユーザ定義関数
296
 
      <literal>getx</literal> および <literal>gety</literal>により
297
 
      セグメント毎に読まれます.
298
 
      これらの関数の呼び出し手順は以下のようになります:
299
 
    </para>
300
 
    <para>
301
 
      <literal>xk=getx(ns,offset)</literal> および
302
 
      <literal>yk=gety(ns,offset)</literal>. ただし,<literal>ns</literal> は
303
 
      セグメントの大きさ, <literal>offset</literal> は
304
 
      完全な信号のセグメントの最初の要素のインデックスです.
305
 
    </para>
306
 
  </refsection>
307
 
  <refsection>
308
 
    <title>参考文献</title>
309
 
    <para>Oppenheim, A.V., and R.W. Schafer. Discrete-Time Signal Processing,
310
 
      Upper Saddle River, NJ: Prentice-Hall, 1999
311
 
    </para>
312
 
  </refsection>
313
 
  <refsection>
314
 
    <title>例</title>
315
 
    <programlisting role="example"><![CDATA[ 
316
 
rand('normal');rand('seed',0);
317
 
x=rand(1:1024-33+1);
318
 
 
319
 
//make low-pass filter with eqfir
320
 
nf=33;bedge=[0 .1;.125 .5];des=[1 0];wate=[1 1];
321
 
h=eqfir(nf,bedge,des,wate);
322
 
 
323
 
//filter white data to obtain colored data 
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);y=real(fft(hf.*xf,1));
327
 
 
328
 
//plot magnitude of filter
329
 
h2=[h 0*ones(1:968)];hf2=fft(h2,-1);hf2=real(hf2.*conj(hf2));
330
 
hsize=max(size(hf2));fr=(1:hsize)/hsize;plot(fr,log(hf2));
331
 
 
332
 
//pspect example
333
 
sm=pspect(100,200,'tr',y);smsize=max(size(sm));fr=(1:smsize)/smsize;
334
 
plot(fr,log(sm));
335
 
rand('unif');
336
 
 ]]></programlisting>
337
 
  </refsection>
338
 
  <refsection role="see also">
339
 
    <title>参照</title>
340
 
    <simplelist type="inline">
341
 
      <member>
342
 
        <link linkend="cspect">cspect</link>
343
 
      </member>
344
 
      <member>
345
 
        <link linkend="pspect">pspect</link>
346
 
      </member>
347
 
      <member>
348
 
        <link linkend="mese">mese</link>
349
 
      </member>
350
 
      <member>
351
 
        <link linkend="window">window</link>
352
 
      </member>
353
 
    </simplelist>
354
 
  </refsection>
355
 
</refentry>