convmtx

畳み込み行列

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

構文

一个= convmtx (h, n)

説明

例

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

例

すべて折りたたむ

畳み込みの効率的な計算

ライブスクリプトを開く

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

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

元= 30;Na = 1000;Nb = 1000;tcnv = 0;tmtx = 0;为kj = 1:Nt a = randn(Na,1);b = randn(注1);Tic n = conv(a,b);tcnv = tcnv + toc;tic c = convmtx(b,Na);d = c *;tmtx = tmtx + toc;结束t1col = [tcnv tmtx]/Nt

t1col =1×20.0007 - 0.0201

t1rat = tcnv \ tmtx

t1rat = 29.0764

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

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

Nchan = 1000;tcnv = 0;tmtx = 0;n = 0 (Na + Nb-1 Nchan);为kj = 1:Nt a = randn(Na,Nchan);b = randn(注1);抽搐为k = 1:Nchan n(:，k) = conv(a(:，k)，b);结束tcnv = tcnv + toc;tic c = convmtx(b,Na);d = c *;tmtx = tmtx + toc;结束tmcol = [tcnv tmtx]/Nt

tmcol =1×20.2539 - 0.0893

tmrat = tcnv / tmtx

tmrat = 2.8431

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

入力引数

すべて折りたたむ

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

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

データ型:单|双

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

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

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