主要内容

bandpower

説明

p= bandpower(x)では、入力信号xの平均が返され。。xが行列の場合、bandpowerは各列の平均強度を個別に計算します。

p= bandpower(x,fs,freqrange)2要素要素ベクトル指定れる周波周波范囲范囲freqrangeの平均返さます。さた周波数内のパワーをを返すにには,fsを入力しなければなりません。bandpowerは修正を使用してfreqrange内の平均パワーを求めます。

p= bandpower(pxx,f,'psd')では、パワー スペクトル密度 (PSD) 推定を積分して計算された平均パワーpxxが返されます。積分値は長方形近似法によって近似されます。入力のfpxxの PSD 推定に対応する周波数のベクトルです。'psd'オプションは入力が PSD 推定であり、時系列データではないことを示します。

p= bandpower(pxx,f,freqrange,'psd')では、周波数範囲freqrangeに含まれる平均パワーが返されます。freqrangeの周波数がfの値と一致しない場合は、最も近い値が使用されます。平均パワーは、パワー スペクトル密度 (PSD) 推定pxxを積分することで計算されます。積分値は長方形近似法によって近似されます。'psd'オプションは入力が PSD 推定であり、時系列データではないことを示します。

すべて折りたたむ

N(0,1)加法性ホワイトガウスノイズを伴う100 Hzの正弦波で構成される信号を生成します。サンプリング周波数は 1 kHz です。平均パワーを求めて 2 ノルムと比較します。

t = 0:0.001:1-0.001;x = cos(2*pi*100*t)+randn(size(t));P = Bandpower(X)
p = 1.5264
l2norm = norm(x,2)^2/numel(x)
l2norm = 1.5264

指定さ周波数内のパワーのををますます。

N(0,1)加法性ホワイトガウスノイズを伴う100 Hzの正弦波で構成される信号を生成します。サンプリング周波数は 1 kHz です。50 Hz ~ 150 Hz の周波数範囲内のパワー全体のパーセント比を求めます。再現可能な結果が必要な場合は、乱数発生器をリセットします。

rng('默认') t = 0:0.001:1-0.001; x = cos(2*pi*100*t)+randn(size(t)); pband = bandpower(x,1000,[50 150]); ptot = bandpower(x,1000,[0 500]); per_power = 100*(pband/ptot)
per_power = 51.9591

ピリオドグラムで PSD 推定を最初に計算することによって平均パワーを求めます。bandpowerに PSD 推定を入力します。

N(0,1)加法性ホワイトガウスノイズを伴う100 Hzの正弦波で構成される信号を生成します。サンプリング周波数は 1 kHz です。ピリオドグラムを求めて'psd'psd推定推定推定で平均计算します结果を时间时间领域で计算しした平均平均パワーととと

t = 0:0.001:1-0.001;Fs = 1000; x = cos(2*pi*100*t)+randn(size(t)); [Pxx,F] = periodogram(x,rectwin(length(x)),length(x),Fs); p = bandpower(Pxx,F,'psd')
p = 1.5264
avpow = norm(x,2)^2/numel(x)
avpow = 1.5264

ピリオドグラムを入力として使い、指定された周波数範囲内のパワー全体のパーセント比を求めます。

N(0,1)100 Hzの伴う伴う伴う伴うので信号信号を生成生成ししますます。。はははははははははははですですですですですですピリオドグラムピリオドグラムピリオドグラムと50 Hz ~150 Hzのの数内パワーのパーセント比を求め。

Fs = 1000; t = 0:1/Fs:1-0.001; x = cos(2*pi*100*t)+randn(size(t)); [Pxx,F] = periodogram(x,rectwin(length(x)),length(x),Fs); pBand = bandpower(Pxx,F,[50 150],'psd');ptot = bandpower(pxx,f,'psd');per_power = 100*(pBand/pTot)
per_power = 49.1798

N(0,1) 加法性ホワイト ガウス ノイズを伴う 3 個の正弦波から構成されるマルチチャネル信号を作成します。正弦波の周波数は 100 Hz、200 Hz および 300 Hz です。サンプリング周波数は 1 kHz で、信号は 1 秒間持続します。

Fs = 1000; t = 0:1/Fs:1-1/Fs; f = [100;200;300]; x = cos(2*pi*f*t)'+randn(length(t),3);

信号の平均強度を求めて 2 ノルムと比較します。

P = Bandpower(X)
p =1×31.5264 1.5382 1.4717
l2norm = dot(x,x)/length(x)
l2norm =1×31.5264 1.5382 1.4717

入力引数

すべて折りたたむ

入力時系列データ、行ベクトルまたは列ベクトルもしくは行列として指定する。xが行列场合そのは独立チャネル扱われます。

例:cos(pi/4*(0:159))'+randn(160,1)は単一チャネルの列ベクトル信号です。

例:cos(pi./[4;2]*(0:159))'+randn(160,2)は 2 チャネルのノイズを含んだ正弦波です。

データ::double|single
複素数のサポート:あり

正のスカラーで指定する入力時系列データのサンプリング周波数

データ::double|single

2 要素の実数値の行ベクトルまたは列ベクトルで指定する、帯域パワーを計算する周波数範囲。入力信号xがn个个を含む,,,freqrangeは次区间でなけれなりません。

  • xが実数値で N が偶数の場合 [0,fs/2]

  • xが実数でnが奇数场合[0,(n - 1)fs/(2N)]

  • xが複素数値で N が偶数の場合 [–(N – 2)fs/(2n),fs/2]

  • xが複素数値で N が奇数の場合 [–(N – 1)fs/(2n),(N – 1)fs/(2N)]

データ::double|single

片側または両側 PSD 推定値。非負の要素をもつ実数値列ベクトルまたは行列として指定します。

パワー スペクトル密度は、デシベル単位ではなく、線形単位で表さなければなりません。デシベル値をパワー値に変換するには、db2powを使用します。

例:[pxx,f] = periodogram(cos(pi./[4;2]*(0:159))'+randn(160,2))2πhzでサンプリングさたノイズノイズノイズノイズノイズチャネル波のピリオドグラムピリオドグラムピリオドグラムピリオドグラム推定推定推定値とについてについて计算されるれるれる周波周波数

データ::double|single

列ベクトルとして指定する周波数ベクトル。周波数ベクトルfpxxの PSD 推定に対応する周波数を含みます。

データ::double|single

出力引数

すべて折りたたむ

非负のとして返す平均帯域パワー。

データ::double|single

参照

[1] Hayes, Monson H. Statistical Digital Signal Processing and Modeling. New York: John Wiley & Sons, 1996.

[2] Stoica, Petre, and Randolph Moses. Spectral Analysis of Signals. Upper Saddle River, NJ: Prentice Hall, 2005.

拡張機能

バージョン履歴

R2013a で導入

参考

|