seqperiod

マルチチャネル信号を生成し、各列の周期を求めます。

x = [4 0 1 6; 
     2 0 2 7; 
     4 0 1 5; 
     2 0 5 6];

p = seqperiod(x)

p = 1×4

     2     1     4     3

x の最初の列は、周期 2 をもっています。x の 2 列目は、周期 1 をもっています。x の 3 列目は、周期性を示していません。そのため、p(3) は、単に x の行数になります。x の 4 列目は、周期 3 をもっています。ただし、周期的なシーケンスの 2 番目の繰り返しは部分的なものです。

各周期的シーケンスが繰り返される回数を計算します。

[~,nr] = seqperiod(x)

nr = 1×4

    2.0000    4.0000    1.0000    1.3333

x の最初の列に、周期的シーケンスが 2 回現れます。2 列目には、1 サンプルのシーケンスがサンプル数だけ繰り返されます。3 列目には、繰り返しはありません。4 列目の繰り返しの数は、残りのサンプルによって表されるシーケンス長の小数部プラス 1 になります。

1 つのチャネルがサンプリング間隔に 4 つの周期をもち、他のチャネルが 2 つの周期をもつように、2 チャネルの正弦波を生成します。正弦波をプロットします。

n = 0:31;
x = cos(2*pi./[8;16].*n)';

plot(n,x,'.-')
axis tight

繰り返されたサブシーケンスの長さと繰り返し数を計算します。1e-5 の絶対許容誤差を指定します。

[p,nr] = seqperiod(x,1e-5)

p = 1×2

     8    16

nr = 1×2

     4     2

最初の 2 つの次元のサイズが 1 となる配列を作成します。3 番目の次元に沿って、配列は繰り返しシーケンスをもちます。

a = permute([5 4 3 5 4 3 5 4],[3 1 2])

a = 
a(:,:,1) =

     5


a(:,:,2) =

     4


a(:,:,3) =

     3


a(:,:,4) =

     5


a(:,:,5) =

     4


a(:,:,6) =

     3


a(:,:,7) =

     5


a(:,:,8) =

     4

繰り返しシーケンスの周期と、配列に含まれる繰り返し数を計算します。関数は、予想どおり、3 番目の次元に沿って機能します。

[p,nr] = seqperiod(a)

p = 3

nr = 2.6667