stFtmag2sig

stft振幅振幅信号の构成构成GydF4y2Ba

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

构文GydF4y2Ba

x = stftmag2sig（s，，，，nFFt）GydF4y2Ba

x = stftmag2sig（s，，，，nFFt，，，，Fs）GydF4y2Ba

x = stftmag2sig（s，，，，nFFt，，，，ts）GydF4y2Ba

x = stftmag2sig（GydF4y2Ba___GydF4y2Ba，，，，name,Value)

[x，t，info] = stftmag2sig（GydF4y2Ba___GydF4y2Ba）GydF4y2Ba

说明GydF4y2Ba

例

XGydF4y2Ba=stFtmag2sig(sGydF4y2Ba，，，，GydF4y2BanfftGydF4y2Ba）GydF4y2Baは，构成れた领域の実信号GydF4y2BaXGydF4y2Baを返しますこれ，，，，，，アルゴリズムGydF4y2Ba短時間フーリエ変換（（STFT) 振幅sGydF4y2Baから推定されたものです。関数は、sGydF4y2Baが離散フーリエ変換 (DFT) 長nfftGydF4y2Baを使用てされと仮定します。GydF4y2Ba

例

XGydF4y2Ba=stFtmag2sig(sGydF4y2Ba，，，，GydF4y2BanfftGydF4y2Ba，，，，GydF4y2BaFSGydF4y2Ba）GydF4y2Baは，GydF4y2BasGydF4y2BaがレートFSGydF4y2Baでサンプリングされたと仮定して、再構成された信号を返します。

XGydF4y2Ba=stFtmag2sig(sGydF4y2Ba，，，，GydF4y2BanfftGydF4y2Ba，，，，GydF4y2BatsGydF4y2Ba）GydF4y2Baは，GydF4y2BasGydF4y2Baがサンプル时间GydF4y2BatsGydF4y2Baでサンプリングされたと仮定して、再構成された信号を返します。

XGydF4y2Ba=stFtmag2sig(___GydF4y2Ba，，，，GydF4y2Ba名称，价值GydF4y2Ba）GydF4y2Baは名前値组の引数をして追加追加オプションししますます。。。オプションオプションオプションにに，，，，ははははウィンドウウィンドウウィンドウウィンドウウィンドウや，，初期位相をを指定するするするメソッドメソッドなど含ま含まれかに追加ます。たとえば，GydF4y2Ba“FrequencyRange”、“单向的”、“InitializePhaseMethod','random'は,ランダムな初期位相の片側 STFT から信号が再構成されることを指定します。

[[GydF4y2BaXGydF4y2Ba，，，，GydF4y2BatGydF4y2Ba，，，，GydF4y2Ba信息GydF4y2Ba] = stftmag2sig（GydF4y2Ba___GydF4y2Ba）GydF4y2Baは，信号が再構成される時点や、再構成プロセスについての情報を含む構造体も返します。

例

すべて折りたたむGydF4y2Ba

stft振幅からの正弦波の再構成

ライブスクリプトを開く

正規化周波数が $π /GydF4y2Ba 60$ ラジアン/サンプルサンプル，値値ががの正弦波サンプルサンプルサンプルサンプル个个についてます。。。信号信号ののGydF4y2Ba

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

stftのからを构成し。の信号再构成后の信号信号をプロットプロットますGydF4y2Ba

xr = stftmag2sig（abs（s），size（s，1））;情节（x）保留GydF4y2Ba上GydF4y2Ba图（XR，GydF4y2Ba' - 'GydF4y2Ba，，，，GydF4y2Ba'LineWidth'，2）保持GydF4y2Ba离开GydF4y2Balegend('Original'，，，，GydF4y2Ba“重建”GydF4y2Ba）GydF4y2Ba

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original, Reconstructed.

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

xz =庭院（[x; zeros（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'LineWidth'，2）保持GydF4y2Ba离开GydF4y2Balegend('Original'，，，，GydF4y2Ba“重建”GydF4y2Ba）GydF4y2Ba

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original, Reconstructed.

計算を繰り返しますが、今回はXGydF4y2Ba2倍がが倍ののセグメントと仮定し，エッジの影响影响を减ら减らます。GydF4y2Ba

xx = cos（pi/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'LineWidth'，2）保持GydF4y2Ba离开GydF4y2Balegend('Original'，，，，GydF4y2Ba“重建”GydF4y2Ba）GydF4y2Ba

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original, Reconstructed.

stft振幅からのオーディオ信号の再構成

ライブスクリプトを開く

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

加载GydF4y2BasplatGydF4y2Baty =（0：长度（y）-1）/fs;GydF4y2Ba百分比听到，输入声音（y，fs）GydF4y2Ba风=锤子（128）;Olen = 64;NFFT = 1024;stft（y，fs，GydF4y2Ba'Window'，风，GydF4y2Ba'OverlapLength'，奥伦，GydF4y2Ba“ fftlength”GydF4y2Ba，，，，nFFt）GydF4y2Ba

Figure contains an axes object. The axes object with title Short-Time Fourier Transform contains an object of type image.

stftの振幅位相计算し。。GydF4y2Ba

s=stft（y，fs，GydF4y2Ba'Window'，风，GydF4y2Ba'OverlapLength'，奥伦，GydF4y2Ba“ fftlength”GydF4y2Ba，nfft）;smag = abs（s）;sphs =角度（s）;GydF4y2Ba

stftのにて信号再。。ますををしたときと同じパラメーターを使用しますます。。GydF4y2BastFtmag2sigは，をに初期化しし，，，回の化をを使用しますGydF4y2Ba

[[X，，，，tX，，，，信息这是给予的=stFtmag2sig(smag,nfft,Fs,'Window'，风，GydF4y2Ba'OverlapLength'，奥伦）;GydF4y2Ba百分比听到，输入声音（x，fs）GydF4y2Ba

元の再后のをを。良く比较比较するためにに，再构成されれたた信号信号信号をGydF4y2Ba

plot(ty,y,tx+500/Fs,x+1) legend('Original'，，，，GydF4y2Ba“重建”GydF4y2Ba，，，，GydF4y2Ba'Location'，，，，GydF4y2Ba'best'）GydF4y2Ba

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original, Reconstructed.

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

impr = info.Inconsistency

impr = 0.0424

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

[[X，，，，tX，，，，信息这是给予的=stFtmag2sig(smag,nfft,Fs,'Window'，风，GydF4y2Ba'OverlapLength'，奥伦，GydF4y2Ba。。。GydF4y2Ba“最大”GydF4y2Ba，200，GydF4y2Ba'初始阶段'GydF4y2Ba，SPHS）;GydF4y2Ba百分比听到，输入声音（x，fs）GydF4y2Baplot(ty,y,tx+500/Fs,-x+1) legend('Original'，，，，GydF4y2Ba“重建”GydF4y2Ba，，，，GydF4y2Ba'Location'，，，，GydF4y2Ba'best'）GydF4y2Ba

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original, Reconstructed.

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

impr = info.Inconsistency

impr = 1.3874e-16

入力引数

すべて折りたたむGydF4y2Ba

`sGydF4y2Ba`-GydF4y2Bastft振幅GydF4y2Ba
行列

stft振幅。行列として指定します。sGydF4y2Baは，，値の信号になければなりませ。。GydF4y2Ba

例:腹肌GydF4y2Ba（（GydF4y2BastFtGydF4y2Ba（（GydF4y2Ba罪GydF4y2Ba（（pi/2*(0:255)),'FFTLength',128))は，正弦波の振幅振幅をしし。。GydF4y2Ba

例:腹肌GydF4y2Ba（（GydF4y2BastFtGydF4y2Ba（（GydF4y2BachirGydF4y2Ba（0：1/1E3：1,25,1,50））））））GydF4y2Baは，1kHz でサンプリングされたチャープの STFT 振幅を指定します。

データ型:罪Gle|GydF4y2Badouble

`nfftGydF4y2Ba`-GydF4y2BaDFT点の数
正の整数スカラーGydF4y2Ba

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

データ型:罪Gle|GydF4y2Badouble

`FSGydF4y2Ba`-GydF4y2BaサンプルレートGydF4y2Ba
2πGydF4y2Ba（既定）|GydF4y2Ba正の数値スカラー

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

`tsGydF4y2Ba`-GydF4y2Baサンプル時間
`期间GydF4y2Ba`スカラーGydF4y2Ba

サンプル時間。期间GydF4y2Baスカラーで指定ます。GydF4y2BatsGydF4y2Baの指定はサンプルレートGydF4y2BaFGydF4y2Ba_sGydF4y2Ba= 1/GydF4y2BatsGydF4y2Baの设定とです。GydF4y2Ba

例:秒（1）GydF4y2Baは，するサンプル间のの秒间秒间秒间时间を表す表すGydF4y2Ba期间GydF4y2Baスカラーです。GydF4y2Ba

データ型:期间GydF4y2Ba

名前と値引数GydF4y2Ba

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

オプションの引数名称，价值GydF4y2Baのコンマ切りペア指定します。GydF4y2Banameは引数名で、价值GydF4y2Baは対応するです。GydF4y2Banameは引用符で囲まなければなりません。name1,Value1,...,NameN,ValueNのに复数名前とのペア引数を任意の顺番で指定指定。。。GydF4y2Ba

`Display`-GydF4y2Ba不整合の表示オプション
`错误的GydF4y2Ba`（既定）|GydF4y2Ba`真的GydF4y2Ba`

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

データ型:逻辑GydF4y2Ba

`FrequencyRange`-GydF4y2Bastft振幅の周波范囲GydF4y2Ba
`'centered'`（既定）|GydF4y2Ba`“二GydF4y2Ba`|GydF4y2Ba`'片面'GydF4y2Ba`

stft振幅振幅周波范囲。GydF4y2Ba'FrequencyRange'と，GydF4y2Ba'centered'，GydF4y2Ba“二GydF4y2Ba，GydF4y2Ba'片面'GydF4y2Baのいずれ构成さコンマ区切りペア指定指定します。GydF4y2Ba

'centered'-GydF4y2BasGydF4y2Baを中央揃え両侧両侧ののとして扱い。。。GydF4y2BanfftGydF4y2Baが偶数の場合、sGydF4y2Baは区間（（–π, π]ラジアン/サンプルサンプルされてい。。GydF4y2BanfftGydF4y2Baが奇数の，GydF4y2BasGydF4y2Baは区間（（–π, π)ラジアン/サンプルサンプルされてい。。時間情報を指定すると、計算区間はそれぞれ（-FGydF4y2Ba_sGydF4y2Ba， FGydF4y2Ba_sGydF4y2Ba/2这是给予的GydF4y2Baサイクル/単位時間、（-FGydF4y2Ba_sGydF4y2Ba， FGydF4y2Ba_sGydF4y2Ba/2）GydF4y2Baサイクル/単位時間となります。ここで、f_sGydF4y2Baはサンプルレート。GydF4y2Ba
“二GydF4y2Ba-GydF4y2BasGydF4y2Baを区间GydF4y2Ba[[0，，，，2π）GydF4y2Baラジアン/サンプルで計算された両側 STFT の振幅として扱います。時間情報を指定した場合、計算区間は[[0， FGydF4y2Ba_sGydF4y2Ba）GydF4y2Baサイクル/単位時間となります。
'片面'GydF4y2Ba-GydF4y2BasGydF4y2Baを片侧stftののとして扱い。。。GydF4y2BanfftGydF4y2Baが偶数の場合、sGydF4y2Baは区間[0,π]GydF4y2Baラジアン/サンプルサンプルされてい。。GydF4y2BanfftGydF4y2Baが奇数の，GydF4y2BasGydF4y2Baは区間[[0，，，，π)ラジアン/サンプルサンプルされてい。。時間情報を指定すると、計算区間はそれぞれ[[0， FGydF4y2Ba_sGydF4y2Ba/2这是给予的GydF4y2Baサイクル/単位時間、[[0， FGydF4y2Ba_sGydF4y2Ba/2）GydF4y2Baサイクル/単位時間となります。ここで、f_sGydF4y2Baはサンプルレート。GydF4y2Ba

データ型:char|GydF4y2Bastring

`不一致Tolerance`-GydF4y2Ba再構成プロセスの不整合許容誤差
`1e-4`（既定）|GydF4y2Ba正のスカラー

再構成プロセスの不整合許容誤差。“不一致的验证”GydF4y2Baと正と构成されるコンマペアとして指定し。。GydF4y2Ba正規化された不整合が许容より场合，再构成は停止します。GydF4y2Ba

データ型:罪Gle|GydF4y2Badouble

`InitializePhaseMethod`-GydF4y2Ba位相の初期化
`“零”GydF4y2Ba`（既定）|GydF4y2Ba`'随机的'GydF4y2Ba`

位相の初期化。“初始izephasemethod”GydF4y2Baと，GydF4y2Ba“零”GydF4y2BaまたはGydF4y2Ba'随机的'GydF4y2Baで構成されるコンマ区切りのペアとして指定します。“初始izephasemethod”GydF4y2BaとGydF4y2Ba'初始阶段'GydF4y2Baののかつのみつのみ指定ます。。GydF4y2Ba

“零”GydF4y2Ba-关数は，位相をゼロとして初期化します。
'随机的'GydF4y2Ba-关数は，区間[[–π, π]で一分布たとして位相を化します。GydF4y2Ba

データ型:char|GydF4y2Bastring

`初始阶段GydF4y2Ba`-GydF4y2Ba初期位相
範囲 [–π, π] の実数値行列

初期位相。GydF4y2Ba'初始阶段'GydF4y2Baと，范囲GydF4y2Ba[[–π, π]の実数値行列で構成されるコンマ区切りペアとして指定します。この行列は、sGydF4y2Baと同じサイズでなければなりません。“初始izephasemethod”GydF4y2BaとGydF4y2Ba'初始阶段'GydF4y2Baののかつのみつのみ指定ます。。GydF4y2Ba

例:angle(stFtGydF4y2Ba（（GydF4y2Barandn（1000,1））））GydF4y2Baは，信号の时间フーリエ変换位相指定指定します。GydF4y2Ba

例:2*pi*(rand(size(stft(randn(1000,1))))-1/2)は，区间GydF4y2Ba[[–π, π]で一しランダムな位相のします。この行列のサイズサイズは，ランダムランダムな信号のの短短短时间GydF4y2Ba

。GydF4y2Ba

データ型:罪Gle|GydF4y2Badouble

`InputTimeDimension`-GydF4y2Ba入力时间次元GydF4y2Ba
`“跨列”GydF4y2Ba`（既定）|GydF4y2Ba`“下班”GydF4y2Ba`

入力时间次元。GydF4y2Ba“ inputtimedimension”GydF4y2Baと，GydF4y2Ba“跨列”GydF4y2BaまたはGydF4y2Ba“下班”GydF4y2Baから构成コンマ区のペアとしてし。。。GydF4y2Ba

“跨列”GydF4y2Ba-关数はsGydF4y2Baの时间列に，周波数がに沿うと见なします。GydF4y2Ba
“下班”GydF4y2Ba-关数はsGydF4y2Baの時間次元が行に沿い、周波数次元が列に沿うと見なします。

データ型:char|GydF4y2Bastring

`最大值GydF4y2Ba`-GydF4y2Ba最大最适化回数GydF4y2Ba
`100GydF4y2Ba`（既定）|GydF4y2Ba正の整数スカラーGydF4y2Ba

最大最适化回数。“最大”GydF4y2Baと正スカラー构成さコンマ区のペア指定します。反复回数がGydF4y2Ba“最大”GydF4y2Baを超える，构成プロセス停止します。GydF4y2Ba

データ型:罪Gle|GydF4y2Badouble

`方法GydF4y2Ba`-GydF4y2Ba信号の再アルゴリズムGydF4y2Ba
`'gla'`（既定）|GydF4y2Ba`'fgla'GydF4y2Ba`|GydF4y2Ba`'legla'`

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

'gla'-Griffin と Lim によって提唱されたオリジナルの再構成アルゴリズム ([[1这是给予的GydF4y2Baを参照）。GydF4y2Ba
'fgla'GydF4y2Ba-Perraudin、Balazs、Søndergaard によって提唱された高速 Griffin-Lim アルゴリズム ([2]GydF4y2Baを参照）。GydF4y2Ba
'legla'-leroux、Kameoka、Ono、Sagayama によって提唱された高速アルゴリズム ([3]GydF4y2Baを参照）。GydF4y2Ba

データ型:char|GydF4y2Bastring

`OverlapLength`-GydF4y2Baオーバーラップする数GydF4y2Ba
ウィンドウの长のGydF4y2Ba`75％GydF4y2Ba`（既定）|GydF4y2Ba非負の整数

邻接する间オーバーラップサンプルの数。GydF4y2Ba'OverlapLength'と，GydF4y2Ba'Window'の长正の整数构成构成区ののペアとして指定します。。信号を正常正常に再再再构成构成GydF4y2Ba'OverlapLength'を、STFT 振幅の生成に使用されるオーバーラップしたセグメントの数と一致させる必要があります。'OverlapLength'75％以下さささとするする指定指定指定すると，长ととなるなるなる最大最大のののの整数整数にににさされれれれますますますますます。。。GydF4y2Ba

データ型:double|GydF4y2Ba罪Gle

`截断者GydF4y2Ba`-GydF4y2Ba`'legla'`の更新规则打ち切り次数GydF4y2Ba
正の整数GydF4y2Ba

'legla'の更新の打ち切り次数。GydF4y2Ba'TruncationOrder'と正でされる区切りペアとしてします。この引数は，GydF4y2Ba'Method'がGydF4y2Ba'legla'に设定て场合にのみ适用，の各各でで更新さされる位相値のの数ををを制御制御GydF4y2Ba'TruncationOrder'は适応をしてされます。GydF4y2Ba

データ型:罪Gle|GydF4y2Badouble

`更新参数仪GydF4y2Ba`-GydF4y2Ba高速griffin-limアルゴリズムアルゴリズムのパラメーターパラメーターGydF4y2Ba
`0.99GydF4y2Ba`（既定）|GydF4y2Ba正のスカラー

高速griffin-limアルゴリズムアルゴリズムのパラメーターパラメーター。'UpdateParameter'と正と构成さコンマ区ペアでします。この引数は，GydF4y2Ba'Method'がGydF4y2Ba'fgla'GydF4y2Baに設定されている場合にのみ適用され、そのメソッドの更新規則に対するパラメーターを指定します。

データ型:罪Gle|GydF4y2Badouble
复素数サポート：GydF4y2BaありGydF4y2Ba

`Window`-GydF4y2Baスペクトルウィンドウ
`hann(128,'periodic')`（既定）|GydF4y2BaベクトルGydF4y2Ba

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

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

例:hann(128,'periodic')とGydF4y2Ba（1-COS（2*PI*（128：-1：1）'/128））/2GydF4y2Baは両方と，GydF4y2BastFtmag2sigで使用れるのウィンドウ指定します。GydF4y2Ba

データ型:double|GydF4y2Ba罪Gle

出力引数

すべて折りたたむGydF4y2Ba

`XGydF4y2Ba`- 再构成れた领域信号信号GydF4y2Ba
ベクトルGydF4y2Ba

再构成れ时间信号。ベクトル返されます。GydF4y2Ba

`tGydF4y2Ba`- 时点GydF4y2Ba
ベクトルGydF4y2Ba

信号がされる。ベクトルとしてれ。。。GydF4y2Ba

`信息GydF4y2Ba`-再構成プロセス情報
構造体

再构成。次フィールドを含む体返さ返されます。GydF4y2Ba

出口GydF4y2Ba-終了フラグ。
- 値がGydF4y2Ba0GydF4y2Baであれば、最大反復回数に到達してアルゴリズムが停止したことを示します。
- 値がGydF4y2Ba1GydF4y2Baであれば许容を満たしてアルゴリズムしたことを示し。GydF4y2Ba
numIterations-合計反復回数。
不一致GydF4y2Ba- 2つのつの间收束にた的改善の平均。GydF4y2Ba
reconstructedPhase- 最后の再构成れた。。GydF4y2Ba
重建GydF4y2Ba-最後の反復における再構成された短時間フーリエ変換。

详细GydF4y2Ba

すべて折りたたむGydF4y2Ba

短時間フーリエ変換

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

信号の STFT は、信号上の長さ $mGydF4y2Ba$ のGydF4y2Ba"解析ウィンドウ"をスライドして、ウィンドウが適用されたデータの離散フーリエ変換を計算することによって計算されます。ウィンドウは、 $rGydF4y2Ba$ サンプルの間隔で元の信号を飛び越えます。ほとんどのウィンドウ関数は、スペクトルリンギングを回避するためにエッジで小さくなります。非ゼロのオーバーラップ長 $lGydF4y2Ba$ 指定さいる，ウィンドウが适用されたセグメントのラップラップ加算加算がががウィンドウウィンドウウィンドウエッジエッジでの信号信号の减衰减衰减衰减衰减衰ををを补正补正补正しますますます。。。ウィンドウウィンドウウィンドウウィンドウ适用适用さささたたたたたた位相をに対して追加れます。。行列行列のは次のように求め求めます。。GydF4y2Ba

$kGydF4y2Ba =GydF4y2Ba ⌊ \frac{{nGydF4y2Ba}_{XGydF4y2Ba} - lGydF4y2Ba}{mGydF4y2Ba - lGydF4y2Ba} ⌋ ，，，，GydF4y2Ba$

ここで，GydF4y2Ba ${nGydF4y2Ba}_{XGydF4y2Ba}$ は元の信号GydF4y2Ba $XGydF4y2Ba （（GydF4y2Ba nGydF4y2Ba ）GydF4y2Ba$ の長さです。⌊⌋GydF4y2Ba记号はを表し。行列内行数は，中央変换および両侧変换変换ののの场合场合场合场合场合场合はGydF4y2BanGydF4y2Ba_DFTGydF4y2Baと同じ，片侧変换ははGydF4y2Ba⌊N_DFTGydF4y2Ba/2⌋ + 1と同じです。GydF4y2Ba

stft行列行列，行列のGydF4y2Ba $mGydF4y2Ba$ 番目の要素が

${XGydF4y2Ba}_{mGydF4y2Ba} （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba =GydF4y2Ba \sum_{nGydF4y2Ba =GydF4y2Ba - \infty}^{\infty} XGydF4y2Ba （（GydF4y2Ba nGydF4y2Ba ）GydF4y2Ba GGydF4y2Ba （（GydF4y2Ba nGydF4y2Ba - mGydF4y2Ba rGydF4y2Ba ）GydF4y2Ba {eGydF4y2Ba}^{- jGydF4y2Ba 2GydF4y2Ba π FGydF4y2Ba nGydF4y2Ba} ，，，，GydF4y2Ba$

である $XGydF4y2Ba （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba =GydF4y2Ba [[GydF4y2Ba \begin{matrix} {XGydF4y2Ba}_{1GydF4y2Ba} （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba & {XGydF4y2Ba}_{2GydF4y2Ba} （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba & {XGydF4y2Ba}_{3GydF4y2Ba} （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba & \dots & {XGydF4y2Ba}_{kGydF4y2Ba} （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba \end{matrix} 这是给予的GydF4y2Ba$ によって指定されます。

ここで，GydF4y2Ba

$GGydF4y2Ba （（GydF4y2Ba nGydF4y2Ba ）GydF4y2Ba$ -長さ $mGydF4y2Ba$ のウィンドウ关数。GydF4y2Ba
${XGydF4y2Ba}_{mGydF4y2Ba} （（GydF4y2Ba FGydF4y2Ba ）GydF4y2Ba$ - 时间GydF4y2Ba $mGydF4y2Ba rGydF4y2Ba$ 付近を中心としたウィンドウが適用されたデータの DFT。
$rGydF4y2Ba$ - 连続するdft间のサイズ。サイズ，ウィンドウのささGydF4y2Ba $mGydF4y2Ba$ とオーバーラップ长GydF4y2Ba $lGydF4y2Ba$ 間の差異です。

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

正規化された不整合

"正規化された不整合"は，最适反复における构成プロセス收束にた改善を测定します。GydF4y2Ba

正规化れ不は次でされます。GydF4y2Ba

$不一致GydF4y2Ba =GydF4y2Ba \frac{”GydF4y2Ba STFT （（GydF4y2Ba ISTFT （（GydF4y2Ba {sGydF4y2Ba}_{美东时间GydF4y2Ba} ）GydF4y2Ba ）GydF4y2Ba - {sGydF4y2Ba}_{美东时间GydF4y2Ba} ”GydF4y2Ba}{”GydF4y2Ba {sGydF4y2Ba}_{美东时间GydF4y2Ba} ”GydF4y2Ba} ，，，，GydF4y2Ba$

ここで，GydF4y2BasGydF4y2Ba_{美东时间GydF4y2Ba}は各推定さた复素短フーリエで，双柱は行列ノルム，GydF4y2BaSTFTは短フーリエ変换，GydF4y2BaISTFTはそのを表します。GydF4y2BastFtmag2sigは，mATLAB^®GydF4y2Ba关数GydF4y2Banormを使用して行列ノルムを計算します。STFTとその逆の詳細については,GydF4y2Ba短時間フーリエ変換とGydF4y2Ba逆短時間フーリエ変換を参照してください。

参照

[1] Griffin，Daniel W.和Jae S. Lim。“经过修改的短期傅立叶变换的信号估计。”IEEE关于声学，语音和信号处理的交易。卷。32，第2号，1984年4月，第236–243页。https://doi.org/10.1109/tassp.1984.1164317。GydF4y2Ba

[2]Perraudin, Nathanaël, Peter Balazs, and Peter L. Søndergaard. "A Fast Griffin-Lim Algorithm." In 2013 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, NY, October 20–23, 2013. https://doi.org/10.1109/WASPAA.2013.6701851.

[3]leroux, Jonathan, Hirokazu Kameoka, Nobutaka Ono, and Shigeki Sagayama. "Fast Signal Reconstruction from Magnitude STFT Spectrogram Based on Spectrogram Consistency." In Proceedings of the 13th International Conference on Digital Audio Effects (DAFx-10), Graz, Austria, September 6–10, 2010.

拡張機能

C/C++ コード生成
mATLAB® Coder™ を使用して C および C++ コードを生成します。

GPU 配列
Parallel Computing Toolbox™ を使用してグラフィックス処理装置 (GPU) 上で実行することにより、コードを高速化します。

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

'legla'メソッドはされません。GydF4y2Ba

详细については、GPU での MATLAB 関数の実行（（Parallel Computing Toolbox)を参照してください。

バージョン履歴

R2020Bで导入GydF4y2Ba

参考GydF4y2Ba

关数GydF4y2Ba

iscola|GydF4y2Baistft|GydF4y2Bapspectrum|GydF4y2Ba频谱图GydF4y2Ba|GydF4y2BastFtGydF4y2Ba

stFtmag2sig

构文GydF4y2Ba

说明GydF4y2Ba

例

stft振幅からの正弦波の再構成

stft振幅からのオーディオ信号の再構成

入力引数

sGydF4y2Ba-GydF4y2Bastft振幅GydF4y2Ba行列

nfftGydF4y2Ba-GydF4y2BaDFT点の数正の整数スカラーGydF4y2Ba

FSGydF4y2Ba-GydF4y2BaサンプルレートGydF4y2Ba2πGydF4y2Ba（既定）|GydF4y2Ba正の数値スカラー

tsGydF4y2Ba-GydF4y2Baサンプル時間期间GydF4y2BaスカラーGydF4y2Ba

名前と値引数GydF4y2Ba

Display-GydF4y2Ba不整合の表示オプション错误的GydF4y2Ba（既定）|GydF4y2Ba真的GydF4y2Ba

FrequencyRange-GydF4y2Bastft振幅の周波范囲GydF4y2Ba'centered'（既定）|GydF4y2Ba“二GydF4y2Ba|GydF4y2Ba'片面'GydF4y2Ba

不一致Tolerance-GydF4y2Ba再構成プロセスの不整合許容誤差1e-4（既定）|GydF4y2Ba正のスカラー

InitializePhaseMethod-GydF4y2Ba位相の初期化“零”GydF4y2Ba（既定）|GydF4y2Ba'随机的'GydF4y2Ba

初始阶段GydF4y2Ba-GydF4y2Ba初期位相範囲 [–π, π] の実数値行列

InputTimeDimension-GydF4y2Ba入力时间次元GydF4y2Ba“跨列”GydF4y2Ba（既定）|GydF4y2Ba“下班”GydF4y2Ba

最大值GydF4y2Ba-GydF4y2Ba最大最适化回数GydF4y2Ba100GydF4y2Ba（既定）|GydF4y2Ba正の整数スカラーGydF4y2Ba

方法GydF4y2Ba-GydF4y2Ba信号の再アルゴリズムGydF4y2Ba'gla'（既定）|GydF4y2Ba'fgla'GydF4y2Ba|GydF4y2Ba'legla'

OverlapLength-GydF4y2Baオーバーラップする数GydF4y2Baウィンドウの长のGydF4y2Ba75％GydF4y2Ba（既定）|GydF4y2Ba非負の整数

截断者GydF4y2Ba-GydF4y2Ba'legla'の更新规则打ち切り次数GydF4y2Ba正の整数GydF4y2Ba

更新参数仪GydF4y2Ba-GydF4y2Ba高速griffin-limアルゴリズムアルゴリズムのパラメーターパラメーターGydF4y2Ba0.99GydF4y2Ba（既定）|GydF4y2Ba正のスカラー

Window-GydF4y2Baスペクトル ウィンドウhann(128,'periodic')（既定）|GydF4y2BaベクトルGydF4y2Ba

出力引数

XGydF4y2Ba- 再构成れた领域信号信号GydF4y2BaベクトルGydF4y2Ba

tGydF4y2Ba- 时点GydF4y2BaベクトルGydF4y2Ba

信息GydF4y2Ba-再構成プロセス情報構造体

详细GydF4y2Ba

短時間フーリエ変換

正規化された不整合

参照

拡張機能

C/C++ コード生成mATLAB® Coder™ を使用して C および C++ コードを生成します。

GPU 配列Parallel Computing Toolbox™ を使用してグラフィックス処理装置 (GPU) 上で実行することにより、コードを高速化します。

バージョン履歴

参考GydF4y2Ba

关数GydF4y2Ba

`sGydF4y2Ba`-GydF4y2Bastft振幅GydF4y2Ba
行列

`nfftGydF4y2Ba`-GydF4y2BaDFT点の数
正の整数スカラーGydF4y2Ba

`FSGydF4y2Ba`-GydF4y2BaサンプルレートGydF4y2Ba
2πGydF4y2Ba（既定）|GydF4y2Ba正の数値スカラー

`tsGydF4y2Ba`-GydF4y2Baサンプル時間
`期间GydF4y2Ba`スカラーGydF4y2Ba

`Display`-GydF4y2Ba不整合の表示オプション
`错误的GydF4y2Ba`（既定）|GydF4y2Ba`真的GydF4y2Ba`

`FrequencyRange`-GydF4y2Bastft振幅の周波范囲GydF4y2Ba
`'centered'`（既定）|GydF4y2Ba`“二GydF4y2Ba`|GydF4y2Ba`'片面'GydF4y2Ba`

`不一致Tolerance`-GydF4y2Ba再構成プロセスの不整合許容誤差
`1e-4`（既定）|GydF4y2Ba正のスカラー

`InitializePhaseMethod`-GydF4y2Ba位相の初期化
`“零”GydF4y2Ba`（既定）|GydF4y2Ba`'随机的'GydF4y2Ba`

`初始阶段GydF4y2Ba`-GydF4y2Ba初期位相
範囲 [–π, π] の実数値行列

`InputTimeDimension`-GydF4y2Ba入力时间次元GydF4y2Ba
`“跨列”GydF4y2Ba`（既定）|GydF4y2Ba`“下班”GydF4y2Ba`

`最大值GydF4y2Ba`-GydF4y2Ba最大最适化回数GydF4y2Ba
`100GydF4y2Ba`（既定）|GydF4y2Ba正の整数スカラーGydF4y2Ba

`方法GydF4y2Ba`-GydF4y2Ba信号の再アルゴリズムGydF4y2Ba
`'gla'`（既定）|GydF4y2Ba`'fgla'GydF4y2Ba`|GydF4y2Ba`'legla'`

`OverlapLength`-GydF4y2Baオーバーラップする数GydF4y2Ba
ウィンドウの长のGydF4y2Ba`75％GydF4y2Ba`（既定）|GydF4y2Ba非負の整数

`截断者GydF4y2Ba`-GydF4y2Ba`'legla'`の更新规则打ち切り次数GydF4y2Ba
正の整数GydF4y2Ba

`更新参数仪GydF4y2Ba`-GydF4y2Ba高速griffin-limアルゴリズムアルゴリズムのパラメーターパラメーターGydF4y2Ba
`0.99GydF4y2Ba`（既定）|GydF4y2Ba正のスカラー

`Window`-GydF4y2Baスペクトルウィンドウ
`hann(128,'periodic')`（既定）|GydF4y2BaベクトルGydF4y2Ba

`XGydF4y2Ba`- 再构成れた领域信号信号GydF4y2Ba
ベクトルGydF4y2Ba

`tGydF4y2Ba`- 时点GydF4y2Ba
ベクトルGydF4y2Ba

`信息GydF4y2Ba`-再構成プロセス情報
構造体

C/C++ コード生成
mATLAB® Coder™ を使用して C および C++ コードを生成します。

GPU 配列
Parallel Computing Toolbox™ を使用してグラフィックス処理装置 (GPU) 上で実行することにより、コードを高速化します。