主要内容gydF4y2Ba

このペ,ジの翻訳は最新ではありません。ここをクリックして,英語の最新版を参照してください。gydF4y2Ba

光谱图gydF4y2Ba

短時間フ,リエ変換を使用したスペクトログラムgydF4y2Ba

説明gydF4y2Ba

例gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba)gydF4y2Baは,入力信号gydF4y2BaxgydF4y2Baの短時間フ,リエ変換を返します。gydF4y2Ba年代gydF4y2Baの各列は,gydF4y2BaxgydF4y2Baの短時間の,時間が局所化された周波数成分の推定を含みます。gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba)gydF4y2Baは,gydF4y2Ba窗口gydF4y2Baを使用して信号をセグメントに分割し,ウィンドウ処理を実行します。gydF4y2Ba

例gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba)gydF4y2Baは,隣り合ったセグメント間でgydF4y2BanoverlapgydF4y2Ba個のサンプルのオ,バ,ラップを使用します。gydF4y2Ba

例gydF4y2Ba

年代gydF4y2Ba=光谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BanfftgydF4y2Ba)gydF4y2Baは,gydF4y2BanfftgydF4y2Baサンプリング点を使用して離散フ,リエ変換を計算します。gydF4y2Ba

[gydF4y2Ba年代gydF4y2Ba,gydF4y2BawgydF4y2Ba,gydF4y2BatgydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Baは,正規化周波数gydF4y2BawgydF4y2Baのベクトルとスペクトログラムを計算した時点gydF4y2BatgydF4y2Baのベクトルを返します。この構文には,前の構文の入力引数を任意に組み合わせて含めることができます。gydF4y2Ba

例gydF4y2Ba

[gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Baは,巡回周波数gydF4y2BafgydF4y2Baのベクトルをサンプルレ,トgydF4y2BafsgydF4y2Baで表して返します。gydF4y2BafsgydF4y2BaはgydF4y2Ba光谱图gydF4y2Baの5番目の入力でなければなりません。サンプルレ,トを入力した場合でも、前のオプション引数の既定値を使用するには、これらの引数を空[]gydF4y2Baとして指定します。gydF4y2Ba

例gydF4y2Ba

[gydF4y2Ba年代gydF4y2Ba,gydF4y2BawgydF4y2Ba,gydF4y2BatgydF4y2Ba=谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BawgydF4y2Ba)gydF4y2Baは,gydF4y2BawgydF4y2Baで指定した正規化周波数におけるスペクトログラムを返します。gydF4y2Ba

[gydF4y2Ba年代gydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BatgydF4y2Ba=谱图(gydF4y2BaxgydF4y2Ba,gydF4y2Ba窗口gydF4y2Ba,gydF4y2BanoverlapgydF4y2Ba,gydF4y2BafgydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Baは,gydF4y2BafgydF4y2Baで指定した巡回周波数におけるスペクトログラムを返します。gydF4y2Ba

[gydF4y2Ba___gydF4y2Ba,gydF4y2BapsgydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Baは,各セグメントのパワ,スペクトル密度(psd)推定またはパワ,スペクトル推定を含む行列gydF4y2BapsgydF4y2Baも返します。gydF4y2Ba

例gydF4y2Ba

[gydF4y2Ba___gydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba“重新分配”)gydF4y2Baは,各PSDまたはパワ,スペクトル推定を,そのエネルギ,の中心の位置に再代入します。十分に局所化された時相成分またはスペクトル成分が信号に含まれる場合,このオプションは,よりシャープなスペクトログラムを生成します。gydF4y2Ba

例gydF4y2Ba

[gydF4y2Ba___gydF4y2Ba,gydF4y2BapsgydF4y2Ba,gydF4y2Ba足球俱乐部gydF4y2Ba,gydF4y2BatcgydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Baは,各PSDまたはパワースペクトル推定のエネルギー中心における周波数および時間を含む2つの行列gydF4y2Ba足球俱乐部gydF4y2BaおよびgydF4y2BatcgydF4y2Baも返します。gydF4y2Ba

[gydF4y2Ba___gydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafreqrangegydF4y2Ba)gydF4y2Baは,gydF4y2BafreqrangegydF4y2Baで指定した周波数範囲でのPSDまたはパワ,スペクトル推定を返します。gydF4y2BafreqrangegydF4y2Baの有効なオプションは,gydF4y2Ba“单向的”gydF4y2Ba、gydF4y2Ba双侧的gydF4y2BaおよびgydF4y2Ba“中心”gydF4y2Baです。gydF4y2Ba

例gydF4y2Ba

[gydF4y2Ba___gydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2Ba名称,值gydF4y2Ba)gydF4y2Baは,名前と値の組の引数を使用して追加オプションを指定します。オプションには,最小しきい値と出力時間次元が含まれます。gydF4y2Ba

[gydF4y2Ba___gydF4y2Ba=谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BaspectrumtypegydF4y2Ba)gydF4y2Baでは,gydF4y2BaspectrumtypegydF4y2BaがgydF4y2Bapsd的gydF4y2Baに指定されている場合はPSD推定が,gydF4y2BaspectrumtypegydF4y2BaがgydF4y2Ba“权力”gydF4y2Baに指定されている場合はパワ,スペクトル推定が返されます。gydF4y2Ba

例gydF4y2Ba

出力引数なしでgydF4y2Ba光谱图(gydF4y2Ba___gydF4y2Ba)gydF4y2Baを使用すると、現在の图ウィンドウにスペクトルグラムがプロットされます。gydF4y2Ba

例gydF4y2Ba

光谱图(gydF4y2Ba___gydF4y2Ba,gydF4y2BafreqlocgydF4y2Ba)gydF4y2Baでは,周波数をプロットする軸を指定します。gydF4y2Ba

例gydF4y2Ba

すべて折りたたむgydF4y2Ba

正弦波の和から構成される信号のサンプルをgydF4y2Ba NgydF4y2Ba xgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 2gydF4y2Ba 4gydF4y2Ba 個生成します。正弦波の正規化周波数は,gydF4y2Ba 2gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 5gydF4y2Ba ラジアン/サンプルおよびgydF4y2Ba 4gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 5gydF4y2Ba ラジアン/サンプルです。高周波数の正弦波の振幅は,他の正弦波の振幅の10倍です。gydF4y2Ba

N = 1024;n = 0: n -1;W0 = 2*pi/5;X = sin(w0*n)+10*sin(2*w0*n);gydF4y2Ba

関数の既定の設定を使用して,短時間フ,リエ変換を計算します。スペクトログラムをプロットします。gydF4y2Ba

S =谱图(x);光谱图(x,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

計算を繰り返します。gydF4y2Ba

  • 長さgydF4y2Ba ngydF4y2Ba 年代gydF4y2Ba cgydF4y2Ba =gydF4y2Ba ⌊gydF4y2Ba NgydF4y2Ba xgydF4y2Ba /gydF4y2Ba 4gydF4y2Ba .gydF4y2Ba 5gydF4y2Ba ⌋gydF4y2Ba のセクションに信号を分割します。gydF4y2Ba

  • ハミングウィンドウを使用して,セクションにウィンドウを適用します。gydF4y2Ba

  • 隣接するセクション間で50%のオバラップを指定します。gydF4y2Ba

  • FFTを計算するには,gydF4y2Ba 马克斯gydF4y2Ba (gydF4y2Ba 2gydF4y2Ba 5gydF4y2Ba 6gydF4y2Ba ,gydF4y2Ba 2gydF4y2Ba pgydF4y2Ba )gydF4y2Ba 点を使用します。ここで,gydF4y2Ba pgydF4y2Ba =gydF4y2Ba ⌈gydF4y2Ba 日志gydF4y2Ba 2gydF4y2Ba ngydF4y2Ba 年代gydF4y2Ba cgydF4y2Ba ⌉gydF4y2Ba です。gydF4y2Ba

2の方法の結果が同じになることを確認します。gydF4y2Ba

Nx =长度(x);nsc =下限(Nx/4.5);Nov =楼层(nsc/2);NFF = max(256,2^nextpow2(nsc));T =谱图(x,hamming(nsc),nov,nff);Maxerr = max(abs(abs(t(:)))-abs(s(:))))gydF4y2Ba
Maxerr = 0gydF4y2Ba

長さが等しく,セクション間のオーバーラップが50%である8つのセクションに信号を分割します。前のステップと同じFFT長を指定します。短時間フ:リエ変換を計算し,前の2:の手順と結果が同じであることを確認します。gydF4y2Ba

Ns = 8;Ov = 0.5;lsc = floor(Nx/(ns-(ns-1)*ov));T =谱图(x,lsc,floor(ov*lsc),nff);Maxerr = max(abs(abs(t(:)))-abs(s(:))))gydF4y2Ba
Maxerr = 0gydF4y2Ba

1 kHzで2秒間サンプリングされた2次チャプgydF4y2BaxgydF4y2Baを生成します。チャ,プの初期周波数は100 Hzで,gydF4y2BatgydF4y2Ba= 1秒で200hzになります。gydF4y2Ba

T = 0:0.001:2;X = chirp(t,100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

xgydF4y2Baのスペクトログラムを計算して表示します。gydF4y2Ba

  • 信号にハミングウィンドウを適用し,長さ128のセクションに分割します。gydF4y2Ba

  • 隣り合ったセクション間のオ,バ,ラップのサンプルを120個に指定します。gydF4y2Ba

  • ⌊gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 周波数およびgydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 2gydF4y2Ba 3.gydF4y2Ba 5gydF4y2Ba 時間ビンでスペクトルを評価します。gydF4y2Ba

光谱图(x, 128120128年,1 e3)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

ハミングウィンドウをブラックマンウィンドウに置き換えます。オ,バ,ラップを60サンプルに減らします。値が上から下へ向かって増加するように時間軸をプロットします。gydF4y2Ba

谱图(x,blackman(128),60,128,1e3) ax = gca;斧子。YDir =gydF4y2Ba“反向”gydF4y2Ba;gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

100 Hzで開始し,gydF4y2BatgydF4y2Ba= 1秒で200hzになる2次チャ,プの各セグメントのPSDを計算して表示します。1kHz のサンプルレート、128 サンプルのセグメント長、および 120 サンプルのオーバーラップを指定します。128 個の DFT 点と既定のハミング ウィンドウを使用します。

Fs = 1000;T = 0:1/fs:2;X = chirp(t,100, 1200,gydF4y2Ba“二次”gydF4y2Ba);fs光谱图(x, 128120128年,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba二次唧唧喳喳的gydF4y2Ba)gydF4y2Ba

图中包含一个轴。标题为二次啁啾的轴包含一个图像类型的对象。gydF4y2Ba

Dcで開始し,gydF4y2BatgydF4y2Ba= 1秒で150 Hzになる,1 kHzでサンプリングされる線形チャープの各セグメントのPSDを計算して表示します。256サンプルのセグメント長,および250サンプルのオバラップを指定します。既定のハミングウィンドウと256個のDFT点を使用します。gydF4y2Ba

X = chirp(t,0,1,150);fs光谱图(x, 256250256年,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba“线性啁啾”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。标题为Linear Chirp的轴包含一个图像类型的对象。gydF4y2Ba

20hzから開始し,gydF4y2BatgydF4y2Ba= 1秒で60 Hzになる1 kHzでサンプリングされる対数チャープの各セグメントのPSDを計算して表示します。256サンプルのセグメント長,および250サンプルのオバラップを指定します。既定のハミングウィンドウと256個のDFT点を使用します。gydF4y2Ba

X = chirp(t,20,1,60,gydF4y2Ba“对数”gydF4y2Ba);光谱图(x, 256250年,[],fs,gydF4y2Ba“桠溪”gydF4y2Ba)标题(gydF4y2Ba对数唧唧喳喳的gydF4y2Ba)gydF4y2Ba

图中包含一个轴。标题为“对数啁啾”的轴包含一个图像类型的对象。gydF4y2Ba

周波数軸に対数スケ,ルを使用します。スペクトログラムは直線になります。gydF4y2Ba

Ax = gca;斧子。YScale =gydF4y2Ba“日志”gydF4y2Ba;gydF4y2Ba

图中包含一个轴。标题为“对数啁啾”的轴包含一个类型为surface的对象。gydF4y2Ba

関数gydF4y2Ba光谱图gydF4y2Baを使用して信号の瞬時周波数を測定し,追跡します。gydF4y2Ba

1 kHzで2秒間サンプリングされた2次チャプを生成します。初期周波数が100hzで,1秒後には200hzに増加するようにチャプを指定します。gydF4y2Ba

Fs = 1000;T = 0:1/fs:2-1/fs;Y = chirp(t,100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

関数gydF4y2Ba光谱图gydF4y2Baに実装された短時間フ,リエ変換を使用して,チャ,プのスペクトルを推定します。信号にハミングウィンドウを適用し,長さ100のセクションに分割します。隣り合ったセクション間のオ,バ,ラップのサンプル80個を指定し,周波数gydF4y2Ba ⌊gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 5gydF4y2Ba 1gydF4y2Ba におけるスペクトルを評価します。gydF4y2Ba

fs光谱图(y, 100年,80100年,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

⌊gydF4y2Ba (gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba -gydF4y2Ba 8gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 0gydF4y2Ba 0gydF4y2Ba -gydF4y2Ba 8gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 9gydF4y2Ba 6gydF4y2Ba の時間点にわたる最大エネルギーをもつ時間——周波数リッジを検出することにより,チャープの周波数を追跡します。スペクトログラムプロットに瞬時周波数を重ね合わせます。gydF4y2Ba

[~,f,t,p] =谱图(y,100,80,100,fs);[冰箱,~,lr] = tfridge(p,f);持有gydF4y2Ba在gydF4y2Baplot3 (t、冰箱、abs (p (lr)),gydF4y2Ba“线宽”gydF4y2Ba, 4)gydF4y2Ba从gydF4y2Ba

图中包含一个轴。坐标轴包含image、line类型的2个对象。gydF4y2Ba

正弦関数的に変化する周波数成分をも512。gydF4y2Ba

N = 512;n = 0: n -1;x = exp(1j*pi*sin(8*n/ n)*32);gydF4y2Ba

チャ,プの中央揃えの短時間フ,リエ変換を計算します。信号を16サンプルがオ,バ,ラップする32サンプルのセグメントに分割します。64個のDFTポントを指定します。スペクトログラムをプロットします。gydF4y2Ba

[标量,fs,ts] =谱图(x,32,16,64,gydF4y2Ba“中心”gydF4y2Ba);64年光谱图(x, 32岁,16日,gydF4y2Ba“中心”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

区間gydF4y2Ba (gydF4y2Ba -gydF4y2Ba πgydF4y2Ba ,gydF4y2Ba πgydF4y2Ba ]gydF4y2Ba に対して64等間隔周波数のスペクトログラムを計算すると,同じ結果が得られます。gydF4y2Ba“中心”gydF4y2Baオプションは必要ありません。gydF4y2Ba

Fintv = -pi+pi/32:pi/32:pi;[vector,fv,tv] =谱图(x,32,16,fintv);fintv光谱图(x, 32岁,16日,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

1 kHzで2秒間サンプリングされたチャプ信号を生成します。初期周波数が100hzで,1秒後には200hzに増加するようにチャプを指定します。gydF4y2Ba

Fs = 1000;t = 0:1/Fs:2;Y = chirp(t,100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

再割り当てされた信号のスペクトログラムを推定します。gydF4y2Ba

  • 信号に形状パラメタgydF4y2Ba βgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 8gydF4y2Ba のカ@ @ザ@ @ウィンドウを適用し,長さ@ @ 128 @ @のセクションに分割します。gydF4y2Ba

  • 隣り合ったセクション間のオ,バ,ラップのサンプルを120個に指定します。gydF4y2Ba

  • ⌊gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 周波数およびgydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 2gydF4y2Ba 3.gydF4y2Ba 5gydF4y2Ba 時間ビンでスペクトルを評価します。gydF4y2Ba

光谱图(y,凯瑟(128年,18),120128年,Fs,gydF4y2Ba“重新分配”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

1 kHzで2秒間サンプリングされたチャプ信号を生成します。初期周波数が100hzで,1秒後には200hzに増加するようにチャプを指定します。gydF4y2Ba

Fs = 1000;t = 0:1/Fs:2;Y = chirp(t,100, 1200,gydF4y2Ba“二次”gydF4y2Ba);gydF4y2Ba

信号の時間依存のパワ,スペクトル密度(psd)を推定します。gydF4y2Ba

  • 信号に形状パラメタgydF4y2Ba βgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 8gydF4y2Ba のカ@ @ザ@ @ウィンドウを適用し,長さ@ @ 128 @ @のセクションに分割します。gydF4y2Ba

  • 隣り合ったセクション間のオ,バ,ラップのサンプルを120個に指定します。gydF4y2Ba

  • ⌊gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 周波数およびgydF4y2Ba ⌊gydF4y2Ba (gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba (gydF4y2Ba xgydF4y2Ba )gydF4y2Ba -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba /gydF4y2Ba (gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 8gydF4y2Ba -gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 0gydF4y2Ba )gydF4y2Ba ⌋gydF4y2Ba =gydF4y2Ba 2gydF4y2Ba 3.gydF4y2Ba 5gydF4y2Ba 時間ビンでスペクトルを評価します。gydF4y2Ba

各PSD推定の重心の周波数および時間を出力します。gydF4y2Ba -gydF4y2Ba 3.gydF4y2Ba 0gydF4y2Ba dBより小さいPSDの要素をゼロに設定します。gydF4y2Ba

[~, ~, ~, pxx, fc, tc) =光谱图(y,凯瑟(128年,18),120128年,Fs,gydF4y2Ba...gydF4y2Ba“MinThreshold”gydF4y2Ba, -30);gydF4y2Ba

重心の周波数および時間の関数として非ゼロの要素をプロットします。gydF4y2Ba

情节(tc (pxx > 0), fc (pxx > 0),gydF4y2Ba“。”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个line类型的对象。gydF4y2Ba

1024 Hzで2秒間サンプリングされた信号を生成します。gydF4y2Ba

nSamp = 2048;Fs = 1024;t = (0:nSamp-1)'/Fs;gydF4y2Ba

最初の1秒間,信号は400hzの正弦波と凹の2次チャ,プで構成されています。開始および終了時の周波数が250 Hzで,最小でも150 Hzに到達するチャープを区間の中間点で対称となるように指定します。gydF4y2Ba

t1 = t(1:nSamp/2);X11 = sin(2*pi*400*t1);x12 = chirp(t1-t1(nSamp/4),150,nSamp/Fs,1750,gydF4y2Ba“二次”gydF4y2Ba);X1 = x11+x12;gydF4y2Ba

信号の他の部分は,減少する周波数の2の線形チャプで構成されています。一方のチャ,プの初期周波数は250hzで,100hzまで減少します。もう一方のチャ,プの初期周波数は400hzで,250hzまで減少します。gydF4y2Ba

t2 = t(nSamp/2+1:nSamp);x21 = chirp(t2400,nSamp/Fs,100);x22 = chirp(t2,550,nSamp/Fs,250);X2 = x21+x22;gydF4y2Ba

ホワ▪▪トガウスノ▪▪ズを信号に付加します。20 dBのS/N比を指定します。再現可能な結果が必要な場合は,乱数発生器をリセットします。gydF4y2Ba

信噪比= 20;rng (gydF4y2Ba“默认”gydF4y2Ba) sig = [x1;x2];sig = sig + randn(size(sig))*std(sig)/db2mag(SNR);gydF4y2Ba

信号のスペクトログラムを計算してプロットします。カイザーウィンドウを長さ63年形状パラメーターgydF4y2Ba βgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 7gydF4y2Ba に,隣接するセクション間のオーバーラップのサンプルを10個少ない数に,FFT長を256にそれぞれ指定します。gydF4y2Ba

Nwin = 63;风= kaiser(nwin,17);Nlap = nwin-10;NFFT = 256;光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

信噪比より小さい値をも要素がゼロに設定されるように,スペクトログラムをしきい値処理します。gydF4y2Ba

光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“MinThreshold”gydF4y2Ba信噪比,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

各PSD推定を,そのエネルギ,中心の位置に再代入します。gydF4y2Ba

光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“再分配”gydF4y2Ba,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

信噪比より小さい値をもつ要素がゼロに設定されるように,再割り当てされたスペクトログラムをしきい値処理します。gydF4y2Ba

光谱图(团体、风、nlap nfft, Fs,gydF4y2Ba“再分配”gydF4y2Ba,gydF4y2Ba“MinThreshold”gydF4y2Ba信噪比,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

減少する2のチャプと広帯域のスプラッタ音を含むオディオ信号を読み込みます。短時間フ,リエ変換を計算します。波形を300サンプルがオ,バ,ラップする400サンプルのセグメントに分割します。スペクトログラムをプロットします。gydF4y2Ba

负载gydF4y2Ba长条木板gydF4y2Ba输入soundsc(y,Fs)gydF4y2BaSg = 400;Ov = 300;光谱图(y, sg, ov [], Fs,gydF4y2Ba“桠溪”gydF4y2Ba) colormapgydF4y2Ba骨gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

関数gydF4y2Ba光谱图gydF4y2Baを使用し,信号のパワ,スペクトル密度(psd)を出力します。gydF4y2Ba

[s,f,t,p] =谱图(y,sg,ov,[],Fs);gydF4y2Ba

関数gydF4y2BamedfreqgydF4y2Baを使用して2のチャプを追跡します。より強い低周波のチャープを検出するには,探索を100 Hzより高い周波数と広帯域音が始まる前の時間に制限します。gydF4y2Ba

F1 = f > 100;T1 = t < 0.75;M1 = medfreq(p(f1,t1),f(f1));gydF4y2Ba

微弱な高周波のチャープを検出するには,探索を2500 Hzより高い周波数と0.3秒~ 0.65秒の時間に制限します。gydF4y2Ba

F2 = f > 2500;T2 = t > 0.3 & t < 0.65;M2 = medfreq(p(f2,t2),f(f2));gydF4y2Ba

結果をスペクトログラムに重ね合わせます。kHzの単位で表すために周波数の値を1000で除算します。gydF4y2Ba

持有gydF4y2Ba在gydF4y2Ba情节(t (t1)、m1/1000gydF4y2Ba“线宽”gydF4y2Ba4)图(t (t2)、m2/1000,gydF4y2Ba“线宽”gydF4y2Ba, 4)gydF4y2Ba从gydF4y2Ba

图中包含一个轴。坐标轴包含图像、直线类型的3个对象。gydF4y2Ba

10 kHzでサンプリングされた2秒間の信号を生成します。信号の瞬時周波数を時間の三角形関数として指定します。gydF4y2Ba

Fs = 10e3;T = 0:1/fs:2;X1 = vco(锯齿(2*pi*t,0.5),[0.1 0.4]*fs,fs);gydF4y2Ba

信号のスペクトログラムを計算してプロットします。長さ256で形状パラメタgydF4y2Ba βgydF4y2Ba =gydF4y2Ba 5gydF4y2Ba のカ▪▪ザ▪▪ウィンドウを使用します。セクション間のオ,バ,ラップのサンプル220個とDFT点512個を指定します。周波数をgydF4y2BaygydF4y2Ba軸にプロットします。既定のカラ,マップとビュ,を使用します。gydF4y2Ba

光谱图(x1,凯瑟(256 5),220512年,fs,gydF4y2Ba“桠溪”gydF4y2Ba)gydF4y2Ba

图中包含一个轴。坐标轴包含一个image类型的对象。gydF4y2Ba

スペクトログラムをウォ,タ,フォ,ルプロットとして表示するようにビュ,を変更します。カラ,マップをgydF4y2Ba骨gydF4y2Baに設定します。gydF4y2Ba

colormap视图(-45、65)gydF4y2Ba骨gydF4y2Ba

图中包含一个轴。坐标轴包含一个surface类型的对象。gydF4y2Ba

入力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

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

例:gydF4y2Ba因为(π/ 4 * (0:159))+ randn (1160)gydF4y2Baは,ホワaaplトガウスノaaplズに含まれる正弦波を指定します。gydF4y2Ba

デ,タ型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba
複素数のサポ,ト:gydF4y2BaありgydF4y2Ba

ウィンドウ。整数、あるいは行ベクトルまたは列ベクトルとして指定します。窗口gydF4y2Baは信号をセグメントに分割するために使用します。gydF4y2Ba

  • 窗口gydF4y2Baが整数の場合,gydF4y2Ba光谱图gydF4y2BaはgydF4y2BaxgydF4y2Baを長さgydF4y2Ba窗口gydF4y2Baのセグメントに分割し,各セグメントにその長さのハミングウィンドウを適用します。gydF4y2Ba

  • 窗口gydF4y2Baがベクトルの場合,gydF4y2Ba光谱图gydF4y2BaはgydF4y2BaxgydF4y2Baをベクトルと同じ長さのセグメントに分割し,gydF4y2Ba窗口gydF4y2Baを使用して各セグメントにウィンドウを適用します。gydF4y2Ba

xgydF4y2Baの長さがgydF4y2BanoverlapgydF4y2Ba個のオバラップサンプルをも整数個のセグメントに厳密に分割できない場合,gydF4y2BaxgydF4y2Baはそれに応じた長さで打切られます。gydF4y2Ba

窗口gydF4y2Baを空として指定した場合,gydF4y2Ba光谱图gydF4y2Baは,gydF4y2BaxgydF4y2BaがgydF4y2BanoverlapgydF4y2Ba個のオーバーラップサンプルをもつ8個のセグメントに分割されているハミングウィンドウを使用します。gydF4y2Ba

利用可能なウィンドウのリストにいては,gydF4y2BaウィンドウgydF4y2Baを参照してください。gydF4y2Ba

例:gydF4y2Ba损害(N + 1)gydF4y2BaとgydF4y2Ba(1-cos(2 *π* (0:N) / N)) / 2gydF4y2Baは,いずれも長さgydF4y2BaNgydF4y2Ba+ 1のハンウィンドウを指定します。gydF4y2Ba

オ,バ,ラップするサンプル数。正の整数で指定します。gydF4y2Ba

  • 窗口gydF4y2Baがスカラ,の場合,gydF4y2BanoverlapgydF4y2BaはgydF4y2Ba窗口gydF4y2Baより小さくなければなりません。gydF4y2Ba

  • 窗口gydF4y2Baがベクトルの場合,gydF4y2BanoverlapgydF4y2BaはgydF4y2Ba窗口gydF4y2Baの長さより小さくなければなりません。gydF4y2Ba

noverlapgydF4y2Baを空として指定した場合,gydF4y2Ba光谱图gydF4y2Baはセグメント間で50%のオバラップが発生する数を使用します。セグメントの長さを指定していない場合,関数によりgydF4y2BanoverlapgydF4y2BaがgydF4y2Ba⌊NgydF4y2BaxgydF4y2Ba/ 4.5⌋gydF4y2Baに設定されます。ここで,ngydF4y2BaxgydF4y2Baは入力信号の長さで,gydF4y2Ba⌊⌋gydF4y2Ba記号は床関数を表します。gydF4y2Ba

DFT点の数。正の整数スカラとして指定します。nfftgydF4y2Baを空として指定した場合,gydF4y2Ba光谱图gydF4y2Baによりパラメ,タ,がgydF4y2Ba马克斯(256 2gydF4y2BapgydF4y2Ba)gydF4y2Baに設定されます。ここでgydF4y2BaP =⌈loggydF4y2Ba2gydF4y2BaNgydF4y2BawgydF4y2Ba⌉gydF4y2Baで,gydF4y2Ba⌈⌉gydF4y2Ba記号は天井関数を表し,gydF4y2Ba

  • 窗口gydF4y2Baがスカラ,の場合は,ngydF4y2BawgydF4y2Ba=gydF4y2Ba窗口gydF4y2Baです。gydF4y2Ba

  • 窗口gydF4y2Baがベクトルの場合は,ngydF4y2BawgydF4y2Ba=gydF4y2Ba长度(gydF4y2Ba窗口gydF4y2Ba)gydF4y2Baです。gydF4y2Ba

正規化周波数。ベクトルとして指定します。wgydF4y2Baは少なくとも2の要素をもたなければなりません。そうでない場合は,関数がgydF4y2BanfftgydF4y2Baとして解釈するためです。正規化周波数の単位はラジアン/サンプルです。

例:gydF4y2Baπ。/[24]gydF4y2Ba

巡回周波数。ベクトルとして指定します。fgydF4y2Baは少なくとも2の要素をもたなければなりません。そうでない場合は,関数がgydF4y2BanfftgydF4y2Baとして解釈するためです。gydF4y2BafgydF4y2Baの単位はサンプルレ,トgydF4y2BafsgydF4y2Baにより指定されます。gydF4y2Ba

サンプルレ,ト。正のスカラで指定します。サンプルレ,トは単位時間あたりのサンプル数です。時間の単位が秒の場合、サンプルレートの単位は Hz です。

PSD推定の周波数範囲。gydF4y2Ba“单向的”gydF4y2Ba、gydF4y2Ba双侧的gydF4y2BaまたはgydF4y2Ba“中心”gydF4y2Baで指定します。実数値信号の場合,既定の設定はgydF4y2Ba“单向的”gydF4y2Baです。複素数値信号の場合,既定の設定はgydF4y2Ba双侧的gydF4y2Baで,gydF4y2Ba“单向的”gydF4y2Baを指定するとエラ,になります。gydF4y2Ba

  • “单向的”gydF4y2Ba-実数の入力信号の片側スペクトログラムを返します。gydF4y2BanfftgydF4y2Baが偶数の場合,gydF4y2BapsgydF4y2BaはgydF4y2BanfftgydF4y2Ba/2 + 1行をも,計算区間はgydF4y2Ba[0,π]gydF4y2Baラジアン/サンプルです。gydF4y2BanfftgydF4y2Baが奇数の場合,gydF4y2BapsgydF4y2Baは(gydF4y2BanfftgydF4y2Ba+ 1)/2行をも,計算区間はgydF4y2Ba[0,π)gydF4y2Baラジアン/サンプルです。gydF4y2BafsgydF4y2Baを指定すると,それぞれの場合の計算区間は[0,gydF4y2BafsgydF4y2Ba/2]サ▪▪クル/単位時間,[0,gydF4y2BafsgydF4y2Ba/2)サ▪▪クル/単位時間となります。gydF4y2Ba

    メモgydF4y2Ba

    この引数がgydF4y2Ba“单向的”gydF4y2Baに設定されている場合,gydF4y2Ba光谱图gydF4y2Baは,正のナ。gydF4y2Ba

  • 双侧的gydF4y2Ba-実信号または複素信号の両側スペクトログラムを返します。gydF4y2BapsgydF4y2BaはgydF4y2BanfftgydF4y2Ba行をも,計算区間はgydF4y2Ba[0, 2π)gydF4y2Baラジアン/サンプルです。gydF4y2BafsgydF4y2Baを指定した場合,計算区間は[0,gydF4y2BafsgydF4y2Ba)サ▪▪クル/単位時間となります。gydF4y2Ba

  • “中心”gydF4y2Ba-中央に揃えた、実信号または複素信号の両側スペクトログラムを返します。gydF4y2BapsgydF4y2BaはgydF4y2BanfftgydF4y2Ba行をもます。gydF4y2BanfftgydF4y2Baが偶数の場合,gydF4y2BapsgydF4y2Baの計算区間はgydF4y2Ba(π-π)gydF4y2Baラジアン/サンプルです。gydF4y2BanfftgydF4y2Baが奇数の場合,gydF4y2BapsgydF4y2Baの計算区間はgydF4y2Ba(π-π)gydF4y2Baラジアン/サンプルです。gydF4y2BafsgydF4y2Baを指定すると,それぞれの場合の計算区間は(-gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/2]サ▪▪クル/単位時間,(-gydF4y2BafsgydF4y2Ba/ 2,gydF4y2BafsgydF4y2Ba/2)サ▪▪クル/単位時間となります。gydF4y2Ba

パワ,スペクトルのスケ,リング。gydF4y2Bapsd的gydF4y2BaまたはgydF4y2Ba“权力”gydF4y2Baで指定します。gydF4y2Ba

  • spectrumtypegydF4y2Baを省略するか,gydF4y2Bapsd的gydF4y2Baを指定すると,パワ,スペクトル密度が返されます。gydF4y2Ba

  • “权力”gydF4y2Baを指定すると,ウィンドウの等価ノズ帯域幅ごとにPSD推定をスケリングします。結果は,各周波数のパワ,の推定です。gydF4y2Ba“重新分配”gydF4y2Baオプションがオンの場合,関数は,再代入の前にPSDを各周波数ビンの幅で積分します。gydF4y2Ba

周波数の表示軸。gydF4y2Ba“xaxis”gydF4y2BaまたはgydF4y2Ba“桠溪”gydF4y2Baで指定します。gydF4y2Ba

  • “xaxis”gydF4y2Ba-周波数がgydF4y2BaxgydF4y2Ba軸に,時間がgydF4y2BaygydF4y2Ba軸に表示されます。gydF4y2Ba

  • “桠溪”gydF4y2Ba-周波数がgydF4y2BaygydF4y2Ba軸に,時間がgydF4y2BaxgydF4y2Ba軸に表示されます。gydF4y2Ba

この引数は出力引数でgydF4y2Ba光谱图gydF4y2Baを呼び出している場合に無視されます。gydF4y2Ba

名前と値のペアの引数gydF4y2Ba

オプションの引数gydF4y2Ba名称,值gydF4y2Baのコンマ区切りペアを指定します。gydF4y2Ba的名字gydF4y2Baは引数名で,gydF4y2Ba价值gydF4y2Baは対応する値です。gydF4y2Ba的名字gydF4y2Baは引用符で囲まなければなりません。gydF4y2BaName1, Value1,…,的家gydF4y2Baのように,複数の名前と値のペアの引数を,任意の順番で指定できます。gydF4y2Ba

例:gydF4y2Ba光谱图(x, 100,“OutputTimeDimension”,“downrows”)gydF4y2BaはgydF4y2BaxgydF4y2Baを長さ100のセグメントに分割し,各セグメントにその長さのハミングウィンドウを適用します。スペクトログラムの出力には,行方向に沿った時間次元があります。gydF4y2Ba

しきい値。gydF4y2BaMinThresholdgydF4y2Baで構成されるコンマ区切りのペアと,デシベル単位で表される実数スカラ,として指定します。gydF4y2Ba光谱图gydF4y2Baは,10 loggydF4y2Ba10gydF4y2Ba(gydF4y2Ba年代gydF4y2Ba)≤gydF4y2Ba打gydF4y2BaとなるgydF4y2Ba年代gydF4y2Baの要素をゼロに設定します。gydF4y2Ba

出力時間次元。gydF4y2BaOutputTimeDimensiongydF4y2Baと,gydF4y2BaacrosscolumnsgydF4y2BaまたはgydF4y2BadownrowsgydF4y2Baから構成されるコンマ区切りのペアとして指定します。行に沿ったgydF4y2Ba年代gydF4y2Ba、gydF4y2BapsgydF4y2Ba、gydF4y2Ba足球俱乐部gydF4y2Ba,およびgydF4y2BatcgydF4y2Baの時間次元と,列に沿った周波数次元が必要な場合は,この値をgydF4y2BadownrowsgydF4y2Baに設定します。列に沿ったgydF4y2Ba年代gydF4y2Ba、gydF4y2BapsgydF4y2Ba、gydF4y2Ba足球俱乐部gydF4y2Ba,およびgydF4y2BatcgydF4y2Baの時間次元と,行に沿った周波数次元が必要な場合は,この値をgydF4y2BaacrosscolumnsgydF4y2Baに設定します。関数が出力引数なしで呼び出された場合、この入力は無視されます。

出力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

短時間フ,リエ変換。行列として返されます。時間は年代gydF4y2Baの列方向に,周波数は行方向に下がって0から増加します。gydF4y2Ba

  • xgydF4y2Baが長さngydF4y2BaxgydF4y2Baの信号の場合,gydF4y2Ba年代gydF4y2Baはk列になります。ここでgydF4y2Ba

    • 窗口gydF4y2Baがスカラ,の場合は,k =⌊(NgydF4y2BaxgydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)/(gydF4y2Ba窗口gydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)⌋。gydF4y2Ba

    • 窗口gydF4y2Baがベクトルの場合は,k =⌊(NgydF4y2BaxgydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)/(gydF4y2Ba长度(gydF4y2Ba窗口gydF4y2Ba)gydF4y2Ba- - - - - -gydF4y2BanoverlapgydF4y2Ba)⌋。gydF4y2Ba

  • xgydF4y2Baが実数でgydF4y2BanfftgydF4y2Baが偶数の場合,gydF4y2Ba年代gydF4y2Baは(gydF4y2BanfftgydF4y2Ba/2 + 1)行になります。gydF4y2Ba

  • xgydF4y2Baが実数でgydF4y2BanfftgydF4y2Baが奇数の場合,gydF4y2Ba年代gydF4y2Baは(gydF4y2BanfftgydF4y2Ba+ 1)/2行になります。gydF4y2Ba

  • xgydF4y2Baが複素数の場合,gydF4y2Ba年代gydF4y2BaはgydF4y2BanfftgydF4y2Ba行になります。gydF4y2Ba

年代gydF4y2BaはgydF4y2Ba“重新分配”gydF4y2Baオプションの影響は受けません。gydF4y2Ba

正規化周波数。ベクトルとして返されます。wgydF4y2Baの長さはgydF4y2Ba年代gydF4y2Baの行数と等しくなります。gydF4y2Ba

時点。ベクトルとして返されます。gydF4y2BatgydF4y2Baの時間値は,各セグメントの中間点に対応します。gydF4y2Ba

巡回周波数。ベクトルとして返されます。fgydF4y2Baの長さはgydF4y2Ba年代gydF4y2Baの行数と等しくなります。gydF4y2Ba

パワスペクトル密度(psd)またはパワスペクトル。行列として返されます。gydF4y2Ba

  • xgydF4y2Baが実数の場合,gydF4y2BapsgydF4y2Baは各セグメントのPSDまたはパワ,スペクトルの片側修正ピリオドグラム推定を含みます。gydF4y2Ba

  • xgydF4y2Baが複素数の場合または周波数のベクトルを指定した場合には,gydF4y2BapsgydF4y2Baは各セグメントのPSDまたはパワ,スペクトルの両側修正ピリオドグラム推定を含みます。gydF4y2Ba

エネルギ,中心の周波数および時間。短時間フリエ変換と同じサズの行列として返されます。サンプルレ,トを指定しない場合,gydF4y2Ba足球俱乐部gydF4y2Baの要素は正規化周波数として返されます。gydF4y2Ba

ヒントgydF4y2Ba

短時間フ,リエ変換にゼロがある場合,デシベルに変換すると,プロットできない負の無限大になります。この問題の発生を避けるため,gydF4y2Ba光谱图gydF4y2Baを出力引数なしで呼び出した場合には短時間フ,リエ変換にgydF4y2Ba每股收益gydF4y2Baが加えられます。gydF4y2Ba

参照gydF4y2Ba

奥本海姆、艾伦五世、罗纳德·w·谢弗和约翰·r·巴克。离散时间信号处理。第2版。上马鞍河,新泽西州:普伦蒂斯霍尔,1999。gydF4y2Ba

Rabiner, Lawrence R.和Ronald W. Schafer。语音信号的数字处理。恩格尔伍德悬崖,新泽西州:Prentice-Hall, 1978。gydF4y2Ba

[3] Chassande-Motin, Éric, François Auger和Patrick Flandrin。“重新分配”。时频分析:概念与方法。弗朗茨·赫拉瓦奇和François奥格编辑。伦敦:ISTE/John Wiley and Sons, 2008。gydF4y2Ba

[4]富洛普,肖恩·A.和凯利·菲兹。计算经时间校正的瞬时频率(重新分配)谱图的算法及其应用。美国声学学会杂志。第119卷,2006年1月,360-371页。gydF4y2Ba

拡張機能gydF4y2Ba

R2006aより前に導入gydF4y2Ba