変換の説明

この関数は, パラメータ値optionで指定したシフト量を用いて 順方向または逆方向の1次元またはN次元離散サイン変換を 行います:

• "dst1"の場合, この関数は 正規化しないDST-I変換を計算します. 長さn のベクトルの1次元変換は以下のようになります:

  

• "dst2"の場合, この関数は 正規化しないDST-II変換を計算します. 長さn のベクトルの1次元変換は以下のようになります:

  

• "dst3"の場合, この関数は 正規化しないDST-III変換を計算します. 長さn のベクトルの1次元変換は以下のようになります:

  

• "dst4"の場合, この関数は 正規化しないDST-IV変換を計算します. 長さn のベクトルの1次元変換は以下のようになります:

  

• "dst"の場合, この関数は 正規化されたDST-I変換を計算します. 長さn のベクトルの1次元変換は以下のようになります:

  

• "idst"の場合, この関数は 正規化されたDST-I変換を計算します. 長さn のベクトルの1次元変換は以下のようになります:

  

多次元のDST変換は,一般に,配列の各次元方向の1次元変換の分離可能な積です. 正規化しない変換の場合,順方向変換の後に逆方向多次元変換を行うと,元の配列が各次元の大きさの積で拡大されたものとなります.

構文の説明

短縮構文

順方向

X=dst(A,-1 [,option]) または X=dst(A [,option]) により, オプション値を指定した順方向変換が得られます. デフォルトは,正規化された DST-I順方向変換です.

Aがベクトルの場合 (1より大きい次元が1つだけの場合), 1次元の変換が行われ, その他の場合にはn次元変換が行われます.

(引数-1は"inverse"ではなく, 指数の符号を意味します).

逆方向

X=dst(A,1 [,option])または X=idst(A [,option])は 逆変換を行います.

Aがベクトルの場合 (1より大きい次元が1つだけの場合), 1次元の変換が行われ, その他の場合にはn次元変換が行われます.

指定した次元方向のDSTの長い構文

• X=dst(A,sign,selection [,option]) により, 選択した次元方向のAの"スライス"の 順方向または逆方向dstを 効率的に計算することができます.

  例えば, A が3次元配列の場合, X=dst(A,-1,2) は以下と等価になります:

  for i1=1:size(A,1),
    for i3=1:size(A,3),
      X(i1,:,i3)=dst(A(i1,:,i3),-1);
    end
  end

  そして X=dst(A,-1,[1 3]) は 以下と等価になります:

  for i2=1:size(A,2),
    X(:,i2,:)=dst(A(:,i2,:),-1);
  end

• X=dst(A,sign,dims,incr) は 古い構文であり,この方法でも 指定した次元方向のAのスライスの 順方向または逆方向のdctを行うことができます.

  例えば,Aが n1*n2*n3 個の要素を有する配列の場合, X=dst(A,-1,n1,1) は X=dst(matrix(A,[n1,n2,n3]),-1,1) と等価で, X=dst(A,-1,[n1 n3],[1 n1*n2]) は X=dst(matrix(A,[n1,n2,n3]),-1,[1,3]) と等価です.

dstを最適化

注意: この関数は直近のパラメータをメモリに保持し,2回目に再利用します. これにより(同じパラメータで)連続的にコールした場合の 計算時間が著しく改善されます.

get_fftw_wisdom, set_fftw_wisdom 関数により dstを更に最適化することができます.