stftgydF4y2Ba

短時間フ，リエ変換gydF4y2Ba

ペ，ジ内をすべて折りたたむgydF4y2Ba

構文gydF4y2Ba

S = stft(x)gydF4y2Ba

S = stft(x,fs)gydF4y2Ba

S = stft(x,ts)gydF4y2Ba

S = stft(gydF4y2Ba___gydF4y2Ba、名称、值)gydF4y2Ba

[s,f] = stft(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

[s,f,t] = stft(gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

stft (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba

説明gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2BaxgydF4y2Ba)gydF4y2Baは,gydF4y2BaxgydF4y2Baの短時間フ，リエ変換(stft)を返します。gydF4y2Ba

例gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2BaxgydF4y2Ba，gydF4y2BafsgydF4y2Ba)gydF4y2Baは，サンプルレ，トgydF4y2BafsgydF4y2Baを使用してgydF4y2BaxgydF4y2BaのSTFTを返します。gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2BaxgydF4y2Ba，gydF4y2BatsgydF4y2Ba)gydF4y2Baは，サンプル時間gydF4y2BatsgydF4y2Baを使用してgydF4y2BaxgydF4y2BaのSTFTを返します。gydF4y2Ba

例gydF4y2Ba

年代gydF4y2Ba= stft (gydF4y2Ba___gydF4y2Ba，gydF4y2Ba名称,值gydF4y2Ba)gydF4y2Baは，名前と値のペアの引数を使用して追加オプションを指定します。オプションにはFFTのウィンドウと長さが含まれます。これらの引数を前の入力構文のいずれかに追加できます。gydF4y2Ba

例gydF4y2Ba

［gydF4y2Ba年代gydF4y2Ba，gydF4y2BafgydF4y2Ba= stft(gydF4y2Ba___gydF4y2Ba)gydF4y2Baは，stftが評価される周波数gydF4y2BafgydF4y2Baを返します。gydF4y2Ba

例gydF4y2Ba

［gydF4y2Ba年代gydF4y2Ba，gydF4y2BafgydF4y2Ba，gydF4y2BatgydF4y2Ba= stft(gydF4y2Ba___gydF4y2Ba)gydF4y2Baは，stftが評価される時間を返します。gydF4y2Ba

出力引数を設定せずにgydF4y2Bastft (gydF4y2Ba___gydF4y2Ba)gydF4y2Baを使用すると、現在の图ウィンドウにSTFTの振幅がプロットされます。gydF4y2Ba

例gydF4y2Ba

すべて折りたたむgydF4y2Ba

3d STFTの可視化gydF4y2Ba

ラ@ @ブスクリプトを開くgydF4y2Ba

10 kHzでサンプリングされる正弦波によって制御された,2秒間の電圧制御発振器出力を生成します。gydF4y2Ba

Fs = 10e3;T = 0:1/fs:2;X = vco(sin(2*pi*t)，[0.1 0.4]*fs,fs);gydF4y2Ba

信号のSTFTを計算してプロットします。長さ256で形状パラメタgydF4y2Ba $βgydF4y2Ba ＝gydF4y2Ba 5gydF4y2Ba$ のカ▪▪ザ▪▪ウィンドウを使用します。オ，バ，ラップの長さを220サンプル，dft長を512点として指定します。既定のカラ，マップとビュ，をも，STFTをプロットします。gydF4y2Ba

stft (x, fs,gydF4y2Ba“窗口”gydF4y2Ba凯瑟(256 5),gydF4y2Ba“OverlapLength”gydF4y2Ba, 220,gydF4y2Ba“FFTLength”gydF4y2Ba, 512);gydF4y2Ba

STFTをウォ，タ，フォ，ルプロットとして表示するようにビュ，を変更します。カラ，マップをgydF4y2Ba飞机gydF4y2Baに設定します。gydF4y2Ba

colormap视图(-45、65)gydF4y2Ba飞机gydF4y2Ba

2次チャ，プのSTFTgydF4y2Ba

ラ@ @ブスクリプトを開くgydF4y2Ba

1 kHzで2秒間サンプリングされた2次チャプを生成します。瞬時周波数は，gydF4y2Ba $tgydF4y2Ba ＝gydF4y2Ba 0gydF4y2Ba$ では100 Hzであり，gydF4y2Ba $tgydF4y2Ba ＝gydF4y2Ba 1gydF4y2Ba$ 秒で200hzになります。gydF4y2Ba

Ts = 0:1/1e3:2;F0 = 100;F1 = 200;X = chirp(ts,f0,1,f1，gydF4y2Ba“二次”gydF4y2Ba[],gydF4y2Ba“凹”gydF4y2Ba）;gydF4y2Ba

持続時間が1 msの2次チャ，プのSTFTを計算して表示します。gydF4y2Ba

D =秒(1e-3);赢=汉明(100，gydF4y2Ba“周期”gydF4y2Ba）;stft (x, d,gydF4y2Ba“窗口”gydF4y2Ba,赢了,gydF4y2Ba“OverlapLength”gydF4y2Ba, 98,gydF4y2Ba“FFTLength”gydF4y2Ba, 128);gydF4y2Ba

STFT周波数範囲gydF4y2Ba

ラ@ @ブスクリプトを開くgydF4y2Ba

5 kHzで4秒間サンプリングされた信号を生成します。信号は,振動する振幅の領域および増加傾向で変動する周波数の領域によって分離された、持続時間が減少していく一連のパルスで構成されています。信号をプロットします。

Fs = 5000;T = 0:1/fs:4-1/fs;X = besselj(0,600*(sin(2*pi*(t+1).^3/30).^5));情节(t, x)gydF4y2Ba

信号の片側，両側，および中心の短時間フ，リエ変換を計算します。いずれの場合も，形状係数gydF4y2Ba $βgydF4y2Ba ＝gydF4y2Ba 10gydF4y2Ba$ で202サンプルのカザウィンドウを使用して，信号セグメントにウィンドウを適用します。各変換の計算に使用される周波数範囲を表示します。gydF4y2Ba

RNGS = [gydF4y2Ba“单向的”gydF4y2Ba“双侧”gydF4y2Ba“中心”gydF4y2Ba];gydF4y2Ba为gydF4y2BaKj = 1:length(rngs) opts = {gydF4y2Ba“窗口”gydF4y2Ba凯瑟(202、10)gydF4y2Ba“FrequencyRange”gydF4y2Ba随机数生成器(kj)};[~，f] = stft(x,fs,opts{:});Subplot (length(rngs)，1,kj) stft(x,fs,opts{:}) title(sprintf(gydF4y2Ba" '%s': [%5.3f， %5.3f] kHz'gydF4y2Ba随机数生成器(kj), (f (1) f(结束)]/ 1000)gydF4y2Ba结束gydF4y2Ba

計算を繰り返しますが，今回はカe e e e e e e e e e e e e e e e e e e e e e e e e e e e egydF4y2Ba双侧的gydF4y2Baの周波数範囲は変わりません。他の2の周波数範囲は，上限でオプンになります。gydF4y2Ba

为gydF4y2BaKj = 1:length(rngs) opts = {gydF4y2Ba“窗口”gydF4y2Ba凯瑟(203、10)gydF4y2Ba“FrequencyRange”gydF4y2Ba随机数生成器(kj)};[~，f] = stft(x,fs,opts{:});Subplot (length(rngs)，1,kj) stft(x,fs,opts{:}) title(sprintf(gydF4y2Ba" '%s': [%5.3f， %5.3f] kHz'gydF4y2Ba随机数生成器(kj), (f (1) f(结束)]/ 1000)gydF4y2Ba结束gydF4y2Ba

マルチチャネル信号のSTFTgydF4y2Ba

ラ@ @ブスクリプトを開くgydF4y2Ba

1 kHzで1秒間サンプリングされた3つの異なるチャープで構成される3チャネルの信号を生成します。gydF4y2Ba

最初のチャネルは凹二次チャ，プで構成され，瞬時周波数はgydF4y2Bat =gydF4y2Ba0では100 Hzであり，gydF4y2Bat =gydF4y2Ba1秒では300hzになります。初期位相は45度です。gydF4y2Ba
2番目のチャネルは凸二次チャ，プで構成され，瞬時周波数はgydF4y2Bat =gydF4y2Ba0では200 Hzであり，gydF4y2Bat =gydF4y2Ba1秒では600hzになります。gydF4y2Ba
3番目のチャネルは対数チャプで構成され，瞬時周波数はgydF4y2Bat =gydF4y2Ba0では300 Hzであり，gydF4y2Bat =gydF4y2Ba1秒では500hzになります。gydF4y2Ba

128年長さがの周期的ハミングウィンドウと50サンプルのオーバーラップ長を使用して,マルチチャネル信号のSTFTを計算します。gydF4y2Ba

Fs = 1e3;T = 0:1/fs:1-1/fs;X = [chirp(t,100, 1300，gydF4y2Ba“二次”gydF4y2Ba45岁的gydF4y2Ba“凹”gydF4y2Ba）;唧唧声(t, 200, 1600,gydF4y2Ba“二次”gydF4y2Ba[],gydF4y2Ba“凸”gydF4y2Ba）;唧唧声(t, 300, 1500,gydF4y2Ba“对数”gydF4y2Ba)];[S,F,T] = stft(x,fs，gydF4y2Ba“窗口”gydF4y2Ba128年,汉明(gydF4y2Ba“周期”gydF4y2Ba)，gydF4y2Ba“OverlapLength”gydF4y2Ba, 50);gydF4y2Ba

各チャネルのSTFTをウォ，タフォ，ルプロットとして可視化します。補助関数gydF4y2BahelperGraphicsOptgydF4y2Baを使用して，座標軸の動作を制御します。gydF4y2Ba

瀑布(F T abs (S (:,: 1))) helperGraphicsOpt (1)gydF4y2Ba

瀑布(F T abs (S (:,: 2))) helperGraphicsOpt (2)gydF4y2Ba

瀑布(F T abs (S (:,:, 3))) helperGraphicsOpt (3)gydF4y2Ba

この補助関数は，現在の座標軸の外観や動作を設定します。gydF4y2Ba

函数gydF4y2Bahelpergraphicopt (ChannelId) ax = gca;斧子。XDir =“反向”gydF4y2Ba；斧子。ZLim = [0 30];ax.Title.String = [gydF4y2Ba输入通道:gydF4y2Banum2str (ChannelId)];ax.XLabel.String =gydF4y2Ba的频率(赫兹)gydF4y2Ba；ax.YLabel.String =gydF4y2Ba的时间(秒)gydF4y2Ba；斧子。View = [30 45];gydF4y2Ba结束gydF4y2Ba

入力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

`xgydF4y2Ba`- - - - - -gydF4y2Ba入力信号gydF4y2Ba
ベクトルgydF4y2Ba|gydF4y2Ba行列gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

入力信号。ベクトル、行列、またはMATLABgydF4y2Ba^®gydF4y2Ba时间表gydF4y2Baとして指定します。gydF4y2Ba

メモgydF4y2Ba

istftgydF4y2Baを使用してgydF4y2Ba年代gydF4y2Baを反転させて，その結果をgydF4y2BaxgydF4y2Baと同じ長さになるようにする場合，gydF4y2Ba长度(长度(x) -noverlap) /((窗口)-noverlap)gydF4y2Baの値は整数でなければなりません。gydF4y2Ba

入力に複数のチャネルがある場合は，各列がチャネルに対応する行列としてgydF4y2BaxgydF4y2Baを指定します。gydF4y2Ba
时间表入力の場合，gydF4y2BaxgydF4y2Baに等間隔に増加する有限の行時間を含めなければなりません。时间表が欠損している場合や時間点が重複している場合、欠損または重複する時間および非等間隔の時間をも时间表の整理gydF4y2Baのヒントを使用して修正できます。gydF4y2Ba
マルチチャネル时间表入力の場合は,行列を含む単一の変数を使用した时间表,またはそれぞれが列ベクトルを含む複数の変数を使用した时间表としてgydF4y2BaxgydF4y2Baを指定します。すべての変数は同じ精度でなければなりません。

xgydF4y2Baの各チャネルの長さは，ウィンドウの長さより大きくなければなりません。gydF4y2Ba

例:gydF4y2Ba唧唧声(0:1/4e3:2,250, 1500,“二次”)gydF4y2Baは，単一チャネルのチャ，プを指定します。gydF4y2Ba

例:gydF4y2Ba时间表(兰德(5,2)' SampleRate ', 1)gydF4y2Baは1 Hzで4秒間サンプリングされた2チャネルの確率変数を指定します。gydF4y2Ba

例:gydF4y2Ba时间表(兰德(5,1),兰德(5、1),“SampleRate”,1)gydF4y2Baは1 Hzで4秒間サンプリングされた2チャネルの確率変数を指定します。gydF4y2Ba

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

`fsgydF4y2Ba`- - - - - -gydF4y2Baサンプルレ，トgydF4y2Ba
`2πgydF4y2Ba`(既定値) |gydF4y2Ba正のスカラgydF4y2Ba

サンプルレ，ト。正のスカラで指定します。この引数は、xgydF4y2Baがベクトルまたは行列の場合にのみ適用されます。gydF4y2Ba

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

`tsgydF4y2Ba`- - - - - -gydF4y2Baサンプル時間gydF4y2Ba
持续时间スカラgydF4y2Ba

サンプル時間。gydF4y2Ba持续时间gydF4y2Baスカラ，で指定します。この引数は，gydF4y2BaxgydF4y2Baがベクトルまたは行列の場合にのみ適用されます。gydF4y2Ba

例:gydF4y2Ba秒(1)gydF4y2Baは，連続する信号サンプル間の1秒間の時間差を表すgydF4y2Ba持续时间gydF4y2Baスカラ，です。gydF4y2Ba

デ，タ型:gydF4y2Ba持续时间gydF4y2Ba

名前と値のペアの引数gydF4y2Ba

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

例:gydF4y2Bastft(“窗口”,赢得“OverlapLength”,50岁,FFTLength, 128)gydF4y2Baは,隣接するセグメント間で50サンプルがオーバーラップし128点のFFTをもつウィンドウgydF4y2Ba赢得gydF4y2Baを使用してデ，タをウィンドウ処理します。gydF4y2Ba

`窗口gydF4y2Ba`- - - - - -gydF4y2BaスペクトルウィンドウgydF4y2Ba
`损害(128年,“周期性”)gydF4y2Ba`(既定値) |gydF4y2BaベクトルgydF4y2Ba

スペクトルウィンドウ。gydF4y2Ba“窗口”gydF4y2Baとベクトルから構成されるコンマ区切りのペアとして指定します。ウィンドウを指定しない場合、またはウィンドウを空として指定する場合、関数は長さが 128 のハンウィンドウを使用します。“窗口”gydF4y2Baの長さは2以上でなければなりません。gydF4y2Ba

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

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

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

`OverlapLengthgydF4y2Ba`- - - - - -gydF4y2Baオ，バ，ラップするサンプル数gydF4y2Ba
ウィンドウの長さのgydF4y2Ba`75％gydF4y2Ba`(既定値) |gydF4y2Ba非負の整数gydF4y2Ba

オ，バ，ラップするサンプル数。gydF4y2Ba“OverlapLength”gydF4y2BaとgydF4y2Ba“窗口”gydF4y2Baの長さより小さな正の整数で構成されるコンマ区切りのペアとして指定します。gydF4y2Ba“OverlapLength”gydF4y2Baを省略するか，空として指定する場合，ウィンドウの長さが75%より少ない最大整数に設定されます。これは，既定のハンウィンドウの96サンプルです。gydF4y2Ba

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

`FFTLengthgydF4y2Ba`- - - - - -gydF4y2BaDFT点の数gydF4y2Ba
`128gydF4y2Ba`(既定値) |gydF4y2Ba正の整数gydF4y2Ba

DFT点の数。gydF4y2Ba“FFTLength”gydF4y2Baと正の整数で構成されるコンマ区切りのペアとして指定します。この値は，ウィンドウの長さ以上でなければなりません。入力信号の長さが DFT 長より小さい場合、データにはゼロがパディングされます。

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

`FrequencyRangegydF4y2Ba`- - - - - -gydF4y2BaSTFT周波数範囲gydF4y2Ba
`“中心”gydF4y2Ba`(既定値) |gydF4y2Ba`双侧的gydF4y2Ba`|gydF4y2Ba`“单向的”gydF4y2Ba`

STFT周波数範囲。gydF4y2Ba“FrequencyRange”gydF4y2Baと,gydF4y2Ba“中心”gydF4y2Ba、gydF4y2Ba双侧的gydF4y2BaまたはgydF4y2Ba“单向的”gydF4y2Baで構成されるコンマ区切りのペアとして指定します。gydF4y2Ba

“中心”gydF4y2Ba-中央揃えの両側STFTを計算します。gydF4y2BanfftgydF4y2Baが偶数の場合，gydF4y2Ba年代gydF4y2Baの計算区間はgydF4y2Ba(π-π)gydF4y2Baラジアン/サンプルです。gydF4y2BanfftgydF4y2Baが奇数の場合，gydF4y2Ba年代gydF4y2Baの計算区間はgydF4y2Ba(π-π)gydF4y2Baラジアン/サンプルです。時間情報を指定すると，計算区間はそれぞれgydF4y2Ba(- fgydF4y2Ba_{年代gydF4y2Ba}fgydF4y2Ba_{年代gydF4y2Ba}/ 2)gydF4y2Baサ电子邮箱クル/単位時間，gydF4y2Ba(- fgydF4y2Ba_{年代gydF4y2Ba}fgydF4y2Ba_{年代gydF4y2Ba}/ 2)gydF4y2Baサ▪▪クル/単位時間となります。ここで，fgydF4y2Ba_{年代gydF4y2Ba}はサンプルレ，トです。gydF4y2Ba
双侧的gydF4y2Ba-区間gydF4y2Ba[0, 2π)gydF4y2Baラジアン/サンプルで両側STFTを計算します。時間情報を指定した場合，計算区間はgydF4y2Ba[0 fgydF4y2Ba_{年代gydF4y2Ba})gydF4y2Baサ▪▪クル/単位時間となります。gydF4y2Ba
“单向的”gydF4y2Ba-片側STFTを計算します。gydF4y2BanfftgydF4y2Baが偶数の場合，gydF4y2Ba年代gydF4y2Baの計算区間はgydF4y2Ba[0,π]gydF4y2Baラジアン/サンプルです。gydF4y2BanfftgydF4y2Baが奇数の場合，gydF4y2Ba年代gydF4y2Baの計算区間はgydF4y2Ba[0,π)gydF4y2Baラジアン/サンプルです。時間情報を指定すると，計算区間はそれぞれgydF4y2Ba[0 fgydF4y2Ba_{年代gydF4y2Ba}/ 2)gydF4y2Baサ电子邮箱クル/単位時間，gydF4y2Ba[0 fgydF4y2Ba_{年代gydF4y2Ba}/ 2)gydF4y2Baサ▪▪クル/単位時間となります。ここで，fgydF4y2Ba_{年代gydF4y2Ba}はサンプルレ，トです。このオプションは，実信号でのみ有効です。gydF4y2Ba

メモgydF4y2Ba

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

例にいては，gydF4y2BaSTFT周波数範囲gydF4y2Baを参照してください。gydF4y2Ba

デ，タ型:gydF4y2Ba字符gydF4y2Ba|gydF4y2Ba字符串gydF4y2Ba

`OutputTimeDimensiongydF4y2Ba`- - - - - -gydF4y2Ba出力時間次元gydF4y2Ba
`“acrosscolumns”gydF4y2Ba`(既定値) |gydF4y2Ba`“downrows”gydF4y2Ba`

出力時間次元。gydF4y2Ba“OutputTimeDimension”gydF4y2Baと,gydF4y2Ba“acrosscolumns”gydF4y2BaまたはgydF4y2Ba“downrows”gydF4y2Baから構成されるコンマ区切りのペアとして指定します。行に沿ったgydF4y2Ba年代gydF4y2Baの時間次元と，列に沿った周波数次元が必要な場合は，この値をgydF4y2Ba“downrows”gydF4y2Baに設定します。列に沿ったgydF4y2Ba年代gydF4y2Baの時間次元と，行に沿った周波数次元が必要な場合は，この値をgydF4y2Ba“acrosscolumns”gydF4y2Baに設定します。関数が出力引数なしで呼び出された場合、この入力は無視されます。

出力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

`年代gydF4y2Ba`—短時間フリエ変換gydF4y2Ba
行列 | 3 次元配列gydF4y2Ba

短時間フ，リエ変換。行列または3.次元配列として返されます。時間は年代gydF4y2Baの列方向に，周波数は行方向に下がって増加します。3.番目の次元が存在する場合は、入力チャネルに対応します。

信号gydF4y2BaxgydF4y2BaにNgydF4y2Ba_xgydF4y2Ba個の時間サンプルがある場合，gydF4y2Ba年代gydF4y2Baにはk列があります。ここで，gydF4y2Bak =⌊(NgydF4y2Ba_xgydF4y2Ba- l) /(马丁)⌋gydF4y2Baです。MはgydF4y2Ba“窗口”gydF4y2Baの長さで，lはgydF4y2Ba“OverlapLength”gydF4y2Baであり，⌊⌋記号は床関数を表します。gydF4y2Ba
年代gydF4y2Baの行数はgydF4y2Ba“FFTLength”gydF4y2Baで指定された値と等しくなります。gydF4y2Ba

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

`fgydF4y2Ba`-周波数gydF4y2Ba
ベクトルgydF4y2Ba

STFTが評価される周波数。ベクトルとして返されます。gydF4y2Ba

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

`tgydF4y2Ba`-時点gydF4y2Ba
ベクトルgydF4y2Ba

時点。ベクトルとして返されます。gydF4y2BatgydF4y2Baには,短時間のパワースペクトルの推定の計算に使用されるデータセグメントの中心に対応する時間値が含まれます。gydF4y2Ba

サンプルレ，トgydF4y2BafsgydF4y2Baが指定された場合，ベクトルには時間値(秒単位)が含まれます。gydF4y2Ba
サンプル時間gydF4y2BatsgydF4y2Baが指定された場合、ベクトルは、入力と同じ時間形式をも运行时长配列です。gydF4y2Ba
時間情報が指定されない場合，ベクトルはサンプル数を含みます。gydF4y2Ba

デ，タ型:gydF4y2Ba双gydF4y2Ba|gydF4y2Ba单gydF4y2Ba

詳細gydF4y2Ba

すべて折りたたむgydF4y2Ba

短時間フ，リエ変換gydF4y2Ba

短時間フーリエ変換(STFT)を使用して,非定常信号の周波数成分が時間の経過と共に変化する様子を解析します。gydF4y2Ba

信号のSTFTは，信号上の長さgydF4y2Ba $米gydF4y2Ba$ のgydF4y2Ba“解析ウィンドウ”gydF4y2Baをスラドして，ウィンドウが適用されたデタの離散フリエ変換を計算することによって計算されます。ウィンドウは，gydF4y2Ba $RgydF4y2Ba$ サンプルの間隔で元の信号を飛び越えます。ほとんどのウィンドウ関数は，スペクトルリンギングを回避するためにエッジで小さくなります。非ゼロのオ，バ，ラップ長gydF4y2Ba $lgydF4y2Ba$ が指定されている場合,ウィンドウが適用されたセグメントのオーバーラップ加算がウィンドウエッジでの信号の減衰を補正します。ウィンドウが適用された各セグメントのDFTは,時間と周波数の各点の振幅と位相を含む行列に対して追加されます。STFT配列の行数はDFT点の数と同じで，列数はgydF4y2Ba

$kgydF4y2Ba ＝gydF4y2Ba ⌊gydF4y2Ba \frac{{NgydF4y2Ba}_{xgydF4y2Ba} -gydF4y2Ba lgydF4y2Ba}{米gydF4y2Ba -gydF4y2Ba lgydF4y2Ba} ⌋gydF4y2Ba ，gydF4y2Ba$

によって指定されます。gydF4y2Ba

ここで，gydF4y2Ba ${NgydF4y2Ba}_{xgydF4y2Ba}$ は元の信号gydF4y2Ba $xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba)gydF4y2Ba$ の長さです。gydF4y2Ba⌊⌋gydF4y2Ba記号は床関数を表します。gydF4y2Ba

STFT行列は，この行列のgydF4y2Ba $米gydF4y2Ba$ 番目の要素がgydF4y2Ba

${XgydF4y2Ba}_{米gydF4y2Ba} （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba ＝gydF4y2Ba {\sumgydF4y2Ba}_{ngydF4y2Ba ＝gydF4y2Ba -gydF4y2Ba ∞gydF4y2Ba}^{∞gydF4y2Ba} xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba)gydF4y2Ba ggydF4y2Ba （gydF4y2Ba ngydF4y2Ba -gydF4y2Ba 米gydF4y2Ba RgydF4y2Ba)gydF4y2Ba {egydF4y2Ba}^{-gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba ngydF4y2Ba} ，gydF4y2Ba$

であるgydF4y2Ba $XgydF4y2Ba （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba ＝gydF4y2Ba ［gydF4y2Ba \begin{matrix} {XgydF4y2Ba}_{1gydF4y2Ba} （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba & {XgydF4y2Ba}_{2gydF4y2Ba} （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba & {XgydF4y2Ba}_{3.gydF4y2Ba} （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba & \dotsgydF4y2Ba & {XgydF4y2Ba}_{kgydF4y2Ba} （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba \end{matrix}]gydF4y2Ba$ によって指定されます。gydF4y2Ba

ここで，gydF4y2Ba

$ggydF4y2Ba （gydF4y2Ba ngydF4y2Ba)gydF4y2Ba$ -長さgydF4y2Ba $米gydF4y2Ba$ のウィンドウ関数。gydF4y2Ba
${XgydF4y2Ba}_{米gydF4y2Ba} （gydF4y2Ba fgydF4y2Ba)gydF4y2Ba$ -時間gydF4y2Ba $米gydF4y2Ba RgydF4y2Ba$ 付近を中心としたウィンドウが適用されたデ，タのdft。gydF4y2Ba
$RgydF4y2Ba$ -連続するDFT間のホップサズ。ホップサesc escズは，ウィンドウの長さgydF4y2Ba $米gydF4y2Ba$ とオ，バ，ラップ長gydF4y2Ba $lgydF4y2Ba$ 間の差異です。gydF4y2Ba

STFTの振幅二乗では，関数のパワ，スペクトル密度のgydF4y2Ba光谱图gydF4y2Ba表現が得られます。gydF4y2Ba

完全再構成gydF4y2Ba

一般に，入力信号のSTFTを計算して逆変換しても，完全再構成は行われません。ISTFTの出力を元の入力信号と可能な限り一致させるには,信号とウィンドウが以下の条件を満たす必要があります。gydF4y2Ba

入力サesc escズ-gydF4y2BaistftgydF4y2Baを使用してgydF4y2BastftgydF4y2Baの出力を逆変換して，その結果を入力信号gydF4y2BaxgydF4y2Baと同じ長さにする場合，gydF4y2Ba $kgydF4y2Ba ＝gydF4y2Ba \frac{（gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba （gydF4y2Ba xgydF4y2Ba)gydF4y2Ba -gydF4y2Ba ngydF4y2Ba ogydF4y2Ba vgydF4y2Ba egydF4y2Ba rgydF4y2Ba lgydF4y2Ba 一个gydF4y2Ba pgydF4y2Ba)gydF4y2Ba}{（gydF4y2Ba lgydF4y2Ba egydF4y2Ba ngydF4y2Ba ggydF4y2Ba tgydF4y2Ba hgydF4y2Ba （gydF4y2Ba wgydF4y2Ba 我gydF4y2Ba ngydF4y2Ba dgydF4y2Ba ogydF4y2Ba wgydF4y2Ba)gydF4y2Ba -gydF4y2Ba ngydF4y2Ba ogydF4y2Ba vgydF4y2Ba egydF4y2Ba rgydF4y2Ba lgydF4y2Ba 一个gydF4y2Ba pgydF4y2Ba)gydF4y2Ba}$ の値は整数でなければなりません。gydF4y2Ba
可乐準拠——可乐準拠ウィンドウを使用して,信号の短時間フーリエ変換を変更していないと仮定します。gydF4y2Ba
パディング- kの値が整数にならないような入力信号の長さの場合,短時間フーリエ変換を計算する前に信号をゼロパディングします。信号を逆変換した後に，余分なゼロを削除します。gydF4y2Ba

参照gydF4y2Ba

Mitra, Sanjit K.数字信号处理:基于计算机的方法，第二版，纽约:McGraw-Hill, 2001。gydF4y2Ba

夏普，布鲁斯。叠加处理的可逆性。gydF4y2Bahttps://gauss256.github.io/blog/cola.htmlgydF4y2Ba，于2019年7月访问。gydF4y2Ba

史密斯，朱利叶斯·奥利翁。频谱音频信号处理。gydF4y2Bahttps://ccrma.stanford.edu/~jos/sasp/gydF4y2Ba，在线书籍，2011年版，2018年11月访问。gydF4y2Ba

拡張機能gydF4y2Ba

高配列gydF4y2Ba
メモリの許容量を超えるような多数の行を含む配列を計算します。gydF4y2Ba

使用上の注意および制限:gydF4y2Ba

“OutputTimeDimension”gydF4y2Baは常に指定しなければならず，gydF4y2Ba“downrows”gydF4y2Baに設定しなければなりません。gydF4y2Ba

詳細にいては，gydF4y2Ba高配列gydF4y2Baを参照してください。gydF4y2Ba

C/ c++コ，ド生成gydF4y2Ba
MATLAB®Coder™を使用してCおよびc++コドを生成します。gydF4y2Ba

使用上の注意および制限:gydF4y2Ba

コ，ド生成では，时间表はサポ，トされていません。gydF4y2Ba

Gpuコ，ド生成gydF4y2Ba
GPU编码器™を使用してNVIDIA GPU®のためのCUDA®コードを生成します。gydF4y2Ba

使用上の注意および制限:gydF4y2Ba

コ，ド生成では，时间表はサポ，トされていません。gydF4y2Ba

Gpu配列gydF4y2Ba
并行计算工具箱™を使用してグラフィックス処理装置(GPU)上で実行することにより,コードを高速化します。gydF4y2Ba

この関数は，gpu配列を完全にサポ，トします。詳細にいては，gydF4y2BaGpuでのmatlab関数の実行gydF4y2Ba(并行计算工具箱)gydF4y2Baを参照してください。gydF4y2Ba

参考gydF4y2Ba

関数gydF4y2Ba

pspectrumgydF4y2Ba|gydF4y2BaistftgydF4y2Ba|gydF4y2BaiscolagydF4y2Ba

トピックgydF4y2Ba

時間-周波数ギャラリgydF4y2Ba

R2019aで導入gydF4y2Ba

stftgydF4y2Ba

構文gydF4y2Ba

説明gydF4y2Ba

例gydF4y2Ba

3d STFTの可視化gydF4y2Ba

2次チャ，プのSTFTgydF4y2Ba

STFT周波数範囲gydF4y2Ba

マルチチャネル信号のSTFTgydF4y2Ba

入力引数gydF4y2Ba

`xgydF4y2Ba`- - - - - -gydF4y2Ba入力信号gydF4y2Ba
ベクトルgydF4y2Ba|gydF4y2Ba行列gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

`fsgydF4y2Ba`- - - - - -gydF4y2Baサンプルレ，トgydF4y2Ba
`2πgydF4y2Ba`(既定値) |gydF4y2Ba正のスカラgydF4y2Ba

`tsgydF4y2Ba`- - - - - -gydF4y2Baサンプル時間gydF4y2Ba
持续时间スカラgydF4y2Ba

名前と値のペアの引数gydF4y2Ba

`窗口gydF4y2Ba`- - - - - -gydF4y2BaスペクトルウィンドウgydF4y2Ba
`损害(128年,“周期性”)gydF4y2Ba`(既定値) |gydF4y2BaベクトルgydF4y2Ba

`OverlapLengthgydF4y2Ba`- - - - - -gydF4y2Baオ，バ，ラップするサンプル数gydF4y2Ba
ウィンドウの長さのgydF4y2Ba`75％gydF4y2Ba`(既定値) |gydF4y2Ba非負の整数gydF4y2Ba

`FFTLengthgydF4y2Ba`- - - - - -gydF4y2BaDFT点の数gydF4y2Ba
`128gydF4y2Ba`(既定値) |gydF4y2Ba正の整数gydF4y2Ba

`FrequencyRangegydF4y2Ba`- - - - - -gydF4y2BaSTFT周波数範囲gydF4y2Ba
`“中心”gydF4y2Ba`(既定値) |gydF4y2Ba`双侧的gydF4y2Ba`|gydF4y2Ba`“单向的”gydF4y2Ba`

`OutputTimeDimensiongydF4y2Ba`- - - - - -gydF4y2Ba出力時間次元gydF4y2Ba
`“acrosscolumns”gydF4y2Ba`(既定値) |gydF4y2Ba`“downrows”gydF4y2Ba`

出力引数gydF4y2Ba

`年代gydF4y2Ba`—短時間フリエ変換gydF4y2Ba
行列 | 3 次元配列gydF4y2Ba

`fgydF4y2Ba`-周波数gydF4y2Ba
ベクトルgydF4y2Ba

`tgydF4y2Ba`-時点gydF4y2Ba
ベクトルgydF4y2Ba

詳細gydF4y2Ba

短時間フ，リエ変換gydF4y2Ba

完全再構成gydF4y2Ba

参照gydF4y2Ba

拡張機能gydF4y2Ba

高配列gydF4y2Ba
メモリの許容量を超えるような多数の行を含む配列を計算します。gydF4y2Ba

C/ c++コ，ド生成gydF4y2Ba
MATLAB®Coder™を使用してCおよびc++コドを生成します。gydF4y2Ba

Gpuコ，ド生成gydF4y2Ba
GPU编码器™を使用してNVIDIA GPU®のためのCUDA®コードを生成します。gydF4y2Ba

Gpu配列gydF4y2Ba
并行计算工具箱™を使用してグラフィックス処理装置(GPU)上で実行することにより,コードを高速化します。gydF4y2Ba

参考gydF4y2Ba

関数gydF4y2Ba

トピックgydF4y2Ba

信号处理工具箱ドキュメンテ，ションgydF4y2Ba

サポトgydF4y2Ba

Matlabによる信号処理向けディ，プラ，ニングgydF4y2Ba

stftgydF4y2Ba

構文gydF4y2Ba

説明gydF4y2Ba

例gydF4y2Ba

3d STFTの可視化gydF4y2Ba

2次チャ，プのSTFTgydF4y2Ba

STFT周波数範囲gydF4y2Ba

マルチチャネル信号のSTFTgydF4y2Ba

入力引数gydF4y2Ba

xgydF4y2Ba- - - - - -gydF4y2Ba入力信号gydF4y2BaベクトルgydF4y2Ba|gydF4y2Ba行列gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

fsgydF4y2Ba- - - - - -gydF4y2Baサンプルレ，トgydF4y2Ba2πgydF4y2Ba(既定値) |gydF4y2Ba正のスカラgydF4y2Ba

tsgydF4y2Ba- - - - - -gydF4y2Baサンプル時間gydF4y2Ba持续时间スカラgydF4y2Ba

名前と値のペアの引数gydF4y2Ba

窗口gydF4y2Ba- - - - - -gydF4y2BaスペクトルウィンドウgydF4y2Ba损害(128年,“周期性”)gydF4y2Ba(既定値) |gydF4y2BaベクトルgydF4y2Ba

OverlapLengthgydF4y2Ba- - - - - -gydF4y2Baオ，バ，ラップするサンプル数gydF4y2Baウィンドウの長さのgydF4y2Ba75％gydF4y2Ba(既定値) |gydF4y2Ba非負の整数gydF4y2Ba

FFTLengthgydF4y2Ba- - - - - -gydF4y2BaDFT点の数gydF4y2Ba128gydF4y2Ba(既定値) |gydF4y2Ba正の整数gydF4y2Ba

FrequencyRangegydF4y2Ba- - - - - -gydF4y2BaSTFT周波数範囲gydF4y2Ba“中心”gydF4y2Ba(既定値) |gydF4y2Ba双侧的gydF4y2Ba|gydF4y2Ba“单向的”gydF4y2Ba

OutputTimeDimensiongydF4y2Ba- - - - - -gydF4y2Ba出力時間次元gydF4y2Ba“acrosscolumns”gydF4y2Ba(既定値) |gydF4y2Ba“downrows”gydF4y2Ba

出力引数gydF4y2Ba

年代gydF4y2Ba—短時間フリエ変換gydF4y2Ba行列 | 3 次元配列gydF4y2Ba

fgydF4y2Ba-周波数gydF4y2BaベクトルgydF4y2Ba

tgydF4y2Ba-時点gydF4y2BaベクトルgydF4y2Ba

詳細gydF4y2Ba

短時間フ，リエ変換gydF4y2Ba

完全再構成gydF4y2Ba

参照gydF4y2Ba

拡張機能gydF4y2Ba

高配列gydF4y2Baメモリの許容量を超えるような多数の行を含む配列を計算します。gydF4y2Ba

C/ c++コ，ド生成gydF4y2BaMATLAB®Coder™を使用してCおよびc++コドを生成します。gydF4y2Ba

Gpuコ，ド生成gydF4y2BaGPU编码器™を使用してNVIDIA GPU®のためのCUDA®コードを生成します。gydF4y2Ba

Gpu配列gydF4y2Ba并行计算工具箱™を使用してグラフィックス処理装置(GPU)上で実行することにより,コードを高速化します。gydF4y2Ba

参考gydF4y2Ba

関数gydF4y2Ba

トピックgydF4y2Ba

信号处理工具箱ドキュメンテ，ションgydF4y2Ba

サポトgydF4y2Ba

Matlabによる信号処理向けディ，プラ，ニングgydF4y2Ba

`xgydF4y2Ba`- - - - - -gydF4y2Ba入力信号gydF4y2Ba
ベクトルgydF4y2Ba|gydF4y2Ba行列gydF4y2Ba|gydF4y2Ba时间表gydF4y2Ba

`fsgydF4y2Ba`- - - - - -gydF4y2Baサンプルレ，トgydF4y2Ba
`2πgydF4y2Ba`(既定値) |gydF4y2Ba正のスカラgydF4y2Ba

`tsgydF4y2Ba`- - - - - -gydF4y2Baサンプル時間gydF4y2Ba
持续时间スカラgydF4y2Ba

`窗口gydF4y2Ba`- - - - - -gydF4y2BaスペクトルウィンドウgydF4y2Ba
`损害(128年,“周期性”)gydF4y2Ba`(既定値) |gydF4y2BaベクトルgydF4y2Ba

`OverlapLengthgydF4y2Ba`- - - - - -gydF4y2Baオ，バ，ラップするサンプル数gydF4y2Ba
ウィンドウの長さのgydF4y2Ba`75％gydF4y2Ba`(既定値) |gydF4y2Ba非負の整数gydF4y2Ba

`FFTLengthgydF4y2Ba`- - - - - -gydF4y2BaDFT点の数gydF4y2Ba
`128gydF4y2Ba`(既定値) |gydF4y2Ba正の整数gydF4y2Ba

`FrequencyRangegydF4y2Ba`- - - - - -gydF4y2BaSTFT周波数範囲gydF4y2Ba
`“中心”gydF4y2Ba`(既定値) |gydF4y2Ba`双侧的gydF4y2Ba`|gydF4y2Ba`“单向的”gydF4y2Ba`

`OutputTimeDimensiongydF4y2Ba`- - - - - -gydF4y2Ba出力時間次元gydF4y2Ba
`“acrosscolumns”gydF4y2Ba`(既定値) |gydF4y2Ba`“downrows”gydF4y2Ba`

`年代gydF4y2Ba`—短時間フリエ変換gydF4y2Ba
行列 | 3 次元配列gydF4y2Ba

`fgydF4y2Ba`-周波数gydF4y2Ba
ベクトルgydF4y2Ba

`tgydF4y2Ba`-時点gydF4y2Ba
ベクトルgydF4y2Ba

高配列gydF4y2Ba
メモリの許容量を超えるような多数の行を含む配列を計算します。gydF4y2Ba

C/ c++コ，ド生成gydF4y2Ba
MATLAB®Coder™を使用してCおよびc++コドを生成します。gydF4y2Ba

Gpuコ，ド生成gydF4y2Ba
GPU编码器™を使用してNVIDIA GPU®のためのCUDA®コードを生成します。gydF4y2Ba

Gpu配列gydF4y2Ba
并行计算工具箱™を使用してグラフィックス処理装置(GPU)上で実行することにより,コードを高速化します。gydF4y2Ba