主要内容gydF4y2Ba

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

stftmag2siggydF4y2Ba

STFT振幅からの信号の再構成gydF4y2Ba

説明gydF4y2Ba

例gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba年代gydF4y2Ba,gydF4y2BanfftgydF4y2Ba)gydF4y2Baは,再構成された時間領域の実信号gydF4y2BaxgydF4y2Baを返します。これは,格里芬林アルゴリズムに基づきgydF4y2Ba短時間フ,リエ変換gydF4y2Ba(stft)振幅gydF4y2Ba年代gydF4y2Baから推定されたものです。関数は,gydF4y2Ba年代gydF4y2Baが離散フ,リエ変換(dft)長gydF4y2BanfftgydF4y2Baを使用して計算されたと仮定します。gydF4y2Ba

例gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba年代gydF4y2Ba,gydF4y2BanfftgydF4y2Ba,gydF4y2BafsgydF4y2Ba)gydF4y2Baは,gydF4y2Ba年代gydF4y2BaがレトgydF4y2BafsgydF4y2Baでサンプリングされたと仮定して,再構成された信号を返します。gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba年代gydF4y2Ba,gydF4y2BanfftgydF4y2Ba,gydF4y2BatsgydF4y2Ba)gydF4y2Baは,gydF4y2Ba年代gydF4y2Baがサンプル時間gydF4y2BatsgydF4y2Baでサンプリングされたと仮定して,再構成された信号を返します。gydF4y2Ba

xgydF4y2Ba= stftmag2sig (gydF4y2Ba___gydF4y2Ba,gydF4y2Ba名称,值gydF4y2Ba)gydF4y2Baは,名前と値の組の引数を使用して追加オプションを指定します。オプションには,fftウィンドウや,初期位相を指定するメソッドなどが含まれます。これらの引数を前の入力構文のいずれかに追加できます。たとえば,gydF4y2Ba“FrequencyRange”、“单向的”,“InitializePhaseMethod”、“随机”gydF4y2Baは,ランダムな初期位相の片側STFTから信号が再構成されることを指定します。gydF4y2Ba

[gydF4y2BaxgydF4y2Ba,gydF4y2BatgydF4y2Ba,gydF4y2Ba信息gydF4y2Ba] = stftmag2sig(gydF4y2Ba___gydF4y2Ba)gydF4y2Baは,信号が再構成される時点や,再構成プロセスにいての情報を含む構造体も返します。gydF4y2Ba

例gydF4y2Ba

すべて折りたたむgydF4y2Ba

正規化周波数がgydF4y2Ba πgydF4y2Ba /gydF4y2Ba 60gydF4y2Ba ラジアン/サンプル,dc値が1の正弦波サンプル512個にいて考えます。信号のSTFTを計算します。gydF4y2Ba

N = 512;X = cos(/60*(0:n-1)')+1;S = stft(x);gydF4y2Ba

STFTの振幅から正弦波を再構成します。元の信号と再構成後の信号をプロットします。gydF4y2Ba

xr = stftmag2sig(abs(S),size(S,1));情节(x)gydF4y2Ba在gydF4y2Ba情节(xr,gydF4y2Ba“——”gydF4y2Ba,gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba传奇(gydF4y2Ba“原始”gydF4y2Ba,gydF4y2Ba“重建”gydF4y2Ba)gydF4y2Ba

图中包含一个轴对象。axis对象包含2个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

計算を繰り返しますが,今回は信号をゼロでパディングして,エッジの影響を減らします。gydF4y2Ba

Xz = circshift([x;0 (n, 1)), n / 2);Sz = stft(xz);xr = stftmag2sig(abs(Sz),size(Sz,1));Xz = Xz (n/2+(1:n));Xr = Xr (n/2+(1:n));情节(xz)gydF4y2Ba在gydF4y2Ba情节(xr,gydF4y2Ba“——”gydF4y2Ba,gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba传奇(gydF4y2Ba“原始”gydF4y2Ba,gydF4y2Ba“重建”gydF4y2Ba)gydF4y2Ba

图中包含一个轴对象。axis对象包含2个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

計算を繰り返しますが,今回はgydF4y2BaxgydF4y2Baを,長さが2倍の信号のセグメントだと仮定して,エッジの影響を減らします。gydF4y2Ba

Xx = cos(/60*(-n/2:n/2+n-1)')+1;Sx = stft(xx);xr = stftmag2sig(abs(Sx),size(Sx,1));Xx = Xx (n/2+(1:n));Xr = Xr (n/2+(1:n));情节(xx)gydF4y2Ba在gydF4y2Ba情节(xr,gydF4y2Ba“——”gydF4y2Ba,gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba传奇(gydF4y2Ba“原始”gydF4y2Ba,gydF4y2Ba“重建”gydF4y2Ba)gydF4y2Ba

图中包含一个轴对象。axis对象包含2个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

減少する2のチャプと広帯域のスプラッタ音を含むオディオ信号を読み込みます。信号は8192 Hzでサンプリングされています。信号のSTFTをプロットします。128年波形をサンプルのセグメントに分割し,ハミングウィンドウを使用してセグメントにウィンドウを適用します。隣接するセグメント間のオ,バ,ラップを64サンプル,fft点を1024に指定します。gydF4y2Ba

负载gydF4y2Ba长条木板gydF4y2Baty =(0:长度(y)-1)/Fs;gydF4y2Ba要听到,输入声音(y,Fs)gydF4y2Ba风=汉明(128);奥伦= 64;NFFT = 1024;stft (y, Fs,gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉,gydF4y2Ba“FFTLength”gydF4y2Banfft)gydF4y2Ba

图中包含一个轴对象。标题为短时傅里叶变换的坐标轴对象包含一个图像类型的对象。gydF4y2Ba

STFTの振幅と位相を計算します。gydF4y2Ba

s = stft(y,Fs,gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉,gydF4y2Ba“FFTLength”gydF4y2Ba, nfft);Smag = abs(s);SPHS =角度(s);gydF4y2Ba

STFTの振幅に基づいて信号を再構成します。STFTを計算したときと同じパラメーターを使用します。既定では、stftmag2siggydF4y2Baは,位相をゼロに初期化し,100回の最適化反復を使用します。gydF4y2Ba

[x,tx,info] = stftmag2sig(smag,nfft,Fs,gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉);gydF4y2Ba要收听,请输入声音(x,Fs)gydF4y2Ba

元の信号と再構成後の信号をプロットします。より良く比較するために,再構成された信号を右上にオフセットします。gydF4y2Ba

情节(泰,y, tx + 500 / Fs x + 1)传说(gydF4y2Ba“原始”gydF4y2Ba,gydF4y2Ba“重建”gydF4y2Ba,gydF4y2Ba“位置”gydF4y2Ba,gydF4y2Ba“最佳”gydF4y2Ba)gydF4y2Ba

图中包含一个轴对象。axis对象包含2个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

最後の2の繰り返し間の収束に向けた相対的改善を出力します。gydF4y2Ba

Impr = info。不一致gydF4y2Ba
Impr = 0.0424gydF4y2Ba

最適化反復の数を2倍にし,初期位相をSTFTから実際の位相に設定して,再構成を改善します。元の信号と再構成後の信号をプロットします。より良く比較するために,再構成された信号の負の値をプロットし,それを右上にオフセットします。gydF4y2Ba

[x,tx,info] = stftmag2sig(smag,nfft,Fs,gydF4y2Ba“窗口”gydF4y2Ba、风能、gydF4y2Ba“OverlapLength”gydF4y2Ba克拉,gydF4y2Ba...gydF4y2Ba“MaxIterations”gydF4y2Ba, 200,gydF4y2Ba“InitialPhase”gydF4y2Ba, sph);gydF4y2Ba要收听,请输入声音(x,Fs)gydF4y2Ba情节(泰,y, tx + 500 / Fs - x + 1)传说(gydF4y2Ba“原始”gydF4y2Ba,gydF4y2Ba“重建”gydF4y2Ba,gydF4y2Ba“位置”gydF4y2Ba,gydF4y2Ba“最佳”gydF4y2Ba)gydF4y2Ba

图中包含一个轴对象。axis对象包含2个line类型的对象。这些对象表示原始的、重建的。gydF4y2Ba

最後の2の繰り返し間の収束に向けた相対的改善を出力します。gydF4y2Ba

Impr = info。不一致gydF4y2Ba
Impr = 1.3874e-16gydF4y2Ba

入力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

STFT振幅。行列として指定します。年代gydF4y2Baは,単一チャネル,実数値の信号に対応しなければなりません。gydF4y2Ba

例:gydF4y2Ba腹肌gydF4y2Ba(gydF4y2BastftgydF4y2Ba(gydF4y2Ba罪gydF4y2Ba(π/ 2 *(0:255)),“FFTLength”,128年)gydF4y2Baは,正弦波のSTFT振幅を指定します。gydF4y2Ba

例:gydF4y2Ba腹肌gydF4y2Ba(gydF4y2BastftgydF4y2Ba(gydF4y2Ba尖声地说gydF4y2Ba1 (0:1/1e3:1, 25日,50)))gydF4y2Baは,1 kHzでサンプリングされたチャ,プのSTFT振幅を指定します。gydF4y2Ba

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

DFT点の数。正の整数スカラとして指定します。この引数が常に必要です。

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

サンプルレ,ト。正の数値スカラとして指定します。

サンプル時間。gydF4y2Ba持续时间gydF4y2Baスカラ,で指定します。gydF4y2BatsgydF4y2Baの指定は,サンプルレ,トgydF4y2BafgydF4y2Ba年代gydF4y2Ba= 1 /gydF4y2BatsgydF4y2Baの設定と等価です。gydF4y2Ba

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

デ,タ型:gydF4y2Ba持续时间gydF4y2Ba

名前と値の引数gydF4y2Ba

オプションの引数のペアをgydF4y2BaName1 = Value1,…,以=家gydF4y2Baとして指定します。ここで,gydF4y2Ba的名字gydF4y2Baは引数名で,gydF4y2Ba价值gydF4y2Baは対応する値です。名前と値の引数は他の引数の後に指定しなければなりませんが、ペアの順序は重要ではありません。

R2021aより前は,名前と値をそれぞれコンマを使って区切り,gydF4y2Ba的名字gydF4y2Baを引用符で囲みました。gydF4y2Ba

例:gydF4y2Ba“FrequencyRange”、“单向的”,“InitializePhaseMethod”、“随机”gydF4y2Baは,ランダムな初期位相の片側STFTから信号が再構成されることを指定します。gydF4y2Ba

不整合の表示オプション。gydF4y2Ba“显示”gydF4y2Baと逻辑値で構成されるコンマ区切りのペアとして指定します。このオプションをgydF4y2Ba真正的gydF4y2Baに設定した場合,gydF4y2Bastftmag2siggydF4y2Baは20回の最適化反復ごとに,正規化された不整合を表示します。また,実行の最後に停止情報も表示します。gydF4y2Ba

デ,タ型:gydF4y2Ba逻辑gydF4y2Ba

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

  • “中心”gydF4y2Ba- - - - - -gydF4y2Ba年代gydF4y2Baを中央揃えの両側STFTの振幅として扱います。gydF4y2BanfftgydF4y2Baが偶数の場合,gydF4y2Ba年代gydF4y2Baは区間gydF4y2Ba(π-π)gydF4y2Baラジアン/サンプルで計算されています。gydF4y2BanfftgydF4y2Baが奇数の場合,gydF4y2Ba年代gydF4y2Baは区間gydF4y2Ba(π-π)gydF4y2Baラジアン/サンプルで計算されています。時間情報を指定すると,計算区間はそれぞれgydF4y2Ba(- fgydF4y2Ba年代gydF4y2BafgydF4y2Ba年代gydF4y2Ba/ 2)gydF4y2Baサ电子邮箱クル/単位時間,gydF4y2Ba(- fgydF4y2Ba年代gydF4y2BafgydF4y2Ba年代gydF4y2Ba/ 2)gydF4y2Baサ▪▪クル/単位時間となります。ここで,fgydF4y2Ba年代gydF4y2Baはサンプルレ,トです。gydF4y2Ba

  • 双侧的gydF4y2Ba- - - - - -gydF4y2Ba年代gydF4y2Baを区間gydF4y2Ba[0, 2π)gydF4y2Baラジアン/サンプルで計算された両側STFTの振幅として扱います。時間情報を指定した場合,計算区間はgydF4y2Ba[0 fgydF4y2Ba年代gydF4y2Ba)gydF4y2Baサ▪▪クル/単位時間となります。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字符串gydF4y2Ba

再構成プロセスの不整合許容誤差。gydF4y2Ba“InconsistencyTolerance”gydF4y2Baと正のスカラ,とで構成されるコンマ区切りペアとして指定します。gydF4y2Ba正規化された不整合gydF4y2Baが許容誤差よりも低い場合,再構成プロセスは停止します。gydF4y2Ba

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

位相の初期化。gydF4y2Ba“InitializePhaseMethod”gydF4y2Baと,gydF4y2Ba“零”gydF4y2BaまたはgydF4y2Ba“随机”gydF4y2Baで構成されるコンマ区切りのペアとして指定します。gydF4y2Ba“InitializePhaseMethod”gydF4y2BaとgydF4y2Ba“InitialPhase”gydF4y2Baのいずれか1のみを指定します。gydF4y2Ba

  • “零”gydF4y2Ba-関数は、位相をゼロとして初期化します。gydF4y2Ba

  • “随机”gydF4y2Ba-関数は、区間gydF4y2Ba(-ππ)gydF4y2Baで一様分布した乱数として位相を初期化します。gydF4y2Ba

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

初期位相。gydF4y2Ba“InitialPhase”gydF4y2Baと,範囲gydF4y2Ba(-ππ)gydF4y2Baの実数値行列で構成されるコンマ区切りペアとして指定します。この行列は,gydF4y2Ba年代gydF4y2Baと同じサ@ @ズでなければなりません。gydF4y2Ba“InitializePhaseMethod”gydF4y2BaとgydF4y2Ba“InitialPhase”gydF4y2Baのいずれか1のみを指定します。gydF4y2Ba

例:gydF4y2Ba角(gydF4y2BastftgydF4y2Ba(gydF4y2BarandngydF4y2Ba(1000,1)))gydF4y2Baは,ランダムな信号の短時間フ,リエ変換の位相を指定します。gydF4y2Ba

例:gydF4y2Ba2 *π*(兰德(大小(stft (randn (1000,1)))) 1/2)gydF4y2Baは,区間gydF4y2Ba(-ππ)gydF4y2Baで一様分布したランダムな位相の行列を指定します。この行列のサ@ @ズは,ランダムな信号の短時間フ@ @リエ変換と同じです。gydF4y2Ba

.gydF4y2Ba

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

入力時間次元。gydF4y2Ba“InputTimeDimension”gydF4y2Baと,gydF4y2Ba“acrosscolumns”gydF4y2BaまたはgydF4y2Ba“downrows”gydF4y2Baから構成されるコンマ区切りのペアとして指定します。gydF4y2Ba

  • “acrosscolumns”gydF4y2Ba-関数はgydF4y2Ba年代gydF4y2Baの時間次元が列に沿い,周波数次元が行に沿うと見なします。gydF4y2Ba

  • “downrows”gydF4y2Ba-関数はgydF4y2Ba年代gydF4y2Baの時間次元が行に沿い,周波数次元が列に沿うと見なします。gydF4y2Ba

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

最大最適化反復回数。gydF4y2Ba“MaxIterations”gydF4y2Baと正の整数スカラ,で構成されるコンマ区切りのペアとして指定します。反復回数がgydF4y2Ba“MaxIterations”gydF4y2Baを超えると,再構成プロセスは停止します。gydF4y2Ba

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

信号の再構成アルゴリズム。gydF4y2Ba“方法”gydF4y2Baと次の値のいずれかで構成されるコンマ区切りのペアとして指定します。gydF4y2Ba

  • “杯子”gydF4y2Ba-格里芬と林によって提唱されたオリジナルの再構成アルゴリズム(gydF4y2Ba[1]gydF4y2Baを参照)。gydF4y2Ba

  • “fgla”gydF4y2Ba——罗丹,Balazs Søndergaardによって提唱された高速Griffin-Limアルゴリズム(gydF4y2Ba[2]gydF4y2Baを参照)。gydF4y2Ba

  • “legla”gydF4y2Ba——勒Roux Kameoka小野,Sagayamaによって提唱された高速アルゴリズム(gydF4y2Ba[3]gydF4y2Baを参照)。gydF4y2Ba

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

隣接するセグメント間でオ,バ,ラップするサンプルの数。gydF4y2Ba“OverlapLength”gydF4y2Baと,gydF4y2Ba“窗口”gydF4y2Baの長さより小さな正の整数で構成されるコンマ区切りのペアとして指定します。信号を正常に再構成させるには,gydF4y2Ba“OverlapLength”gydF4y2Baを,stft振幅の生成に使用されるオ,バ,ラップしたセグメントの数と一致させる必要があります。gydF4y2Ba“OverlapLength”gydF4y2Baは,省略するか空として指定すると,ウィンドウの長さの75%以下となる最大の整数に設定されます。これは,既定のハンウィンドウの96サンプルです。gydF4y2Ba

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

“legla”gydF4y2Baの更新規則の打切り次数。gydF4y2Ba“TruncationOrder”gydF4y2Baと正の整数で構成されるコンマ区切りのペアとして指定します。この引数は,gydF4y2Ba“方法”gydF4y2BaがgydF4y2Ba“legla”gydF4y2Baに設定されている場合にのみ適用され,そのメソッドの各反復で更新される位相値の数を制御します。指定しない場合,gydF4y2Ba“TruncationOrder”gydF4y2Baは適応アルゴリズムを使用して決定されます。gydF4y2Ba

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

高速Griffin-Limアルゴリズムの更新パラメ,タ,。gydF4y2Ba“UpdateParameter”gydF4y2Baと正のスカラ,とで構成されるコンマ区切りペアで指定します。この引数は,gydF4y2Ba“方法”gydF4y2BaがgydF4y2Ba“fgla”gydF4y2Baに設定されている場合にのみ適用され,そのメソッドの更新規則に対するパラメ,タ,を指定します。gydF4y2Ba

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

スペクトルウィンドウ。gydF4y2Ba“窗口”gydF4y2Baとベクトルから構成されるコンマ区切りのペアとして指定します。信号を正常に再構成させるには,gydF4y2Ba“窗口”gydF4y2Baを,stft振幅の生成に使用されるウィンドウと一致させる必要があります。ウィンドウを指定しない場合、またはウィンドウを空として指定する場合、関数は長さが 128 の周期的ハン ウィンドウを使用します。“窗口”gydF4y2Baの長さは2以上でなければなりません。gydF4y2Ba

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

例:gydF4y2Ba损害(128年,“周期性”)gydF4y2BaとgydF4y2Ba(1-cos(2 *π* (128:1:1)/ 128))/ 2gydF4y2Baは両方とも,gydF4y2Bastftmag2siggydF4y2Baで使用される既定のウィンドウを指定します。gydF4y2Ba

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

出力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

再構成された時間領域信号。ベクトルとして返されます。gydF4y2Ba

信号が再構成される時点。ベクトルとして返されます。

再構成プロセス情報。次のフィ,ルドを含む構造体として返されます。gydF4y2Ba

  • ExitFlaggydF4y2Ba-終了フラグ。gydF4y2Ba

    • 値がgydF4y2Ba0gydF4y2Baであれば,最大反復回数に到達してアルゴリズムが停止したことを示します。gydF4y2Ba

    • 値がgydF4y2Ba1gydF4y2Baであれば,相対許容誤差を満たしてアルゴリズムが停止したことを示します。gydF4y2Ba

  • NumIterationsgydF4y2Ba-合計反復回数。gydF4y2Ba

  • 不一致gydF4y2Ba—最後の2の繰り返し間の収束に向けた相対的改善の平均。gydF4y2Ba

  • ReconstructedPhasegydF4y2Ba-最後の反復における再構成された位相。gydF4y2Ba

  • ReconstructedSTFTgydF4y2Ba—最後の反復における再構成された短時間フリエ変換。gydF4y2Ba

詳細gydF4y2Ba

すべて折りたたむgydF4y2Ba

短時間フ,リエ変換gydF4y2Ba

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

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

kgydF4y2Ba =gydF4y2Ba ⌊gydF4y2Ba NgydF4y2Ba xgydF4y2Ba −gydF4y2Ba lgydF4y2Ba 米gydF4y2Ba −gydF4y2Ba lgydF4y2Ba ⌋gydF4y2Ba

ここで,gydF4y2Ba NgydF4y2Ba xgydF4y2Ba は元の信号gydF4y2Ba xgydF4y2Ba (gydF4y2Ba ngydF4y2Ba )gydF4y2Ba の長さです。gydF4y2Ba⌊⌋gydF4y2Ba記号は床関数を表します。行列内の行数は、中央変換および両側変換の場合は DFT 点の数であるNgydF4y2BaDFTgydF4y2Baと同じで,片側変換の場合はgydF4y2Ba⌊NgydF4y2BaDFTgydF4y2Ba/2⌋+ 1gydF4y2Baと同じです。gydF4y2Ba

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

XgydF4y2Ba 米gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba 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 XgydF4y2Ba 1gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba XgydF4y2Ba 2gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba XgydF4y2Ba 3.gydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba ⋯gydF4y2Ba XgydF4y2Ba kgydF4y2Ba (gydF4y2Ba fgydF4y2Ba )gydF4y2Ba ]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

“正規化された不整合”gydF4y2Baは,連続する最適化反復における再構成プロセスの収束に向けた改善を測定します。gydF4y2Ba

正規化された不整合は次で定義されます。gydF4y2Ba

不一致gydF4y2Ba =gydF4y2Ba 为gydF4y2Ba STFTgydF4y2Ba (gydF4y2Ba ISTFTgydF4y2Ba (gydF4y2Ba 年代gydF4y2Ba 美国东部时间gydF4y2Ba )gydF4y2Ba )gydF4y2Ba −gydF4y2Ba 年代gydF4y2Ba 美国东部时间gydF4y2Ba 为gydF4y2Ba 为gydF4y2Ba 年代gydF4y2Ba 美国东部时间gydF4y2Ba 为gydF4y2Ba ,gydF4y2Ba

ここで,gydF4y2Ba年代gydF4y2Ba美国东部时间gydF4y2Baは各反復で推定された複素短時間フ,リエ変換で,双柱は行列ノルム,gydF4y2BaSTFTgydF4y2Baは短時間フ,リエ変換,gydF4y2BaISTFTgydF4y2Baはその逆を表します。gydF4y2Bastftmag2siggydF4y2Baは,matlabgydF4y2Ba®gydF4y2Ba関数gydF4y2Ba规范gydF4y2Baを使用して行列ノルムを計算します。STFTとその逆の詳細については、短時間フ,リエ変換gydF4y2BaとgydF4y2Ba逆短時間フ,リエ変換gydF4y2Baを参照してください。gydF4y2Ba

参照gydF4y2Ba

[1]格里芬,丹尼尔W.和杰s林。修正短时间傅里叶变换下的信号估计IEEE声学、语音与信号处理汇刊。第32卷,第2期,1984年4月,第236-243页。https://doi.org/10.1109/TASSP.1984.1164317。gydF4y2Ba

[2] Perraudin, Nathanaël, Peter Balazs和Peter L. Søndergaard。“快速格里芬-利姆算法。”2013年IEEE信号处理在音频和声学中的应用研讨会,New Paltz, NY, 2013年10月20-23日。https://doi.org/10.1109/WASPAA.2013.6701851。gydF4y2Ba

[3]勒鲁,乔纳森,龟冈博和,小野信孝,坂山茂树。基于谱图一致性的星等STFT谱图快速信号重建。第13届国际数字音频效果会议(DAFx-10),格拉茨,奥地利,2010年9月6日至10日。gydF4y2Ba

拡張機能gydF4y2Ba

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

バ,ジョン履歴gydF4y2Ba

R2020bで導入gydF4y2Ba

参考gydF4y2Ba

関数gydF4y2Ba