convmtx

畳み込み行列

ページ内をすべて折りたたむ

構文

A = convmtx(h,n)

説明

例

A = convmtx(h,n) では、畳み込み行列 A が返されます。A と n 要素ベクトル x との積は h と x の畳み込みになります。

例

すべて折りたたむ

畳み込みの効率的な計算

ライブスクリプトを開く

信号がベクトルの場合、一般に、畳み込みの計算に conv を使用すると convmtx を使用するより効率的です。マルチチャネル信号では、convmtx の方が効率的な場合があります。

conv と convmtx の両方を使用して、2 つのランダムなベクトル a と b の畳み込みを計算します。信号のサンプルはそれぞれ 1000 個です。この 2 つの関数が費やした時間を比較します。計算を 30 回繰り返して平均化することで不規則な変動を排除します。

Nt = 30;
Na = 1000;
Nb = 1000;

tcnv = 0;
tmtx = 0;

for kj = 1:Nt
    a = randn(Na,1);
    b = randn(Nb,1);

    tic
    n = conv(a,b);
    tcnv = tcnv+toc;

    tic
    c = convmtx(b,Na);
    d = c*a;
    tmtx = tmtx+toc;
end

t1col = [tcnv tmtx]/Nt

t1col = 1×2

    0.0004    0.0210

t1rat = tcnv\tmtx

t1rat = 48.4760

conv の方がほぼ 2 桁効率的です。

a が 1000 チャネルのマルチチャネル信号である場合について、この実行を繰り返します。事前割り当てを行い、conv のパフォーマンスを最適化します。

Nchan = 1000;

tcnv = 0;
tmtx = 0;

n = zeros(Na+Nb-1,Nchan);

for kj = 1:Nt
    a = randn(Na,Nchan);
    b = randn(Nb,1);
    
    tic
    for k = 1:Nchan
        n(:,k) = conv(a(:,k),b);
    end
    tcnv = tcnv+toc;

    tic
    c = convmtx(b,Na);
    d = c*a;
    tmtx = tmtx+toc;
end

tmcol = [tcnv tmtx]/Nt

tmcol = 1×2

    0.2357    0.1020

tmrat = tcnv/tmtx

tmrat = 2.3104

convmtx は conv の約 3 倍効率的です。

入力引数

すべて折りたたむ

`h` — 入力ベクトル
ベクトル

入力ベクトル。行または列として指定します。

データ型: single | double

`n` — 畳み込みを行うベクトルの長さ
正の整数

畳み込みを行うベクトルの長さ。正の整数として指定します。

h が長さ m の列ベクトルの場合、A は (m+n-1) 行 n 列の行列であり、A と長さ n の列ベクトル x との積は、h と x の畳み込みになります。
h が長さ m の行ベクトルの場合、A は n 行 (m+n-1) 列の行列であり、A と長さ n の行ベクトル x との積は、h と x の畳み込みになります。