主要内容GydF4y2Ba

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

wvdGydF4y2Ba

维格纳维尔分布と平滑化疑似 维格纳维尔分布GydF4y2Ba

説明GydF4y2Ba

例GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2BaxGydF4y2Ba)GydF4y2Baは、GydF4y2BaxGydF4y2Baの 维格纳维尔分布を返します。GydF4y2Ba

例GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2BaxGydF4y2Ba,GydF4y2Ba财政司司长GydF4y2Ba)GydF4y2Baは、GydF4y2BaxGydF4y2BaがレートGydF4y2Ba财政司司长GydF4y2Baでサンプリングされたときの 维格纳维尔分布を返します。GydF4y2Ba

例GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2BaxGydF4y2Ba,GydF4y2BatsGydF4y2Ba)GydF4y2Baは、GydF4y2BaxGydF4y2Baがサンプル間の時間間隔GydF4y2BatsGydF4y2Baでサンプリングされたときの 维格纳维尔分布を返します。GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2Ba___GydF4y2Ba,'Smoothdpseudo')GydF4y2Baは、GydF4y2BaxGydF4y2Baの平滑化疑似 维格纳维尔分布を返します。この関数は、入力信号の長さを使用して、時間と周波数の平滑化に使用されるウィンドウの長さを選択します。この構文には、前の構文の入力引数を任意に組み合わせて含めることができます。GydF4y2Ba

例GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2Ba___GydF4y2Ba,'Smoothdpseudo',GydF4y2Ba成双的GydF4y2Ba,GydF4y2Ba弗温GydF4y2Ba)GydF4y2Baは、平滑化に使用される時間ウィンドウ (GydF4y2Ba成双的GydF4y2Ba) と周波数ウィンドウ (GydF4y2Ba弗温GydF4y2Ba) を指定します。時間または周波数の平滑化に既定のウィンドウを使用するには、対応する引数を空 (GydF4y2Ba[]GydF4y2Ba) に指定します。GydF4y2Ba

例GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2Ba___GydF4y2Ba,'Smoothdpseudo',GydF4y2Ba名称、值GydF4y2Ba)GydF4y2Baは、名前と値のペアの引数を使用して、平滑化疑似 维格纳维尔分布の追加オプションを指定します。この構文では、GydF4y2Ba成双的GydF4y2BaとGydF4y2Ba弗温GydF4y2Baを指定できます。また、省略することもできます。GydF4y2Ba

例GydF4y2Ba

DGydF4y2Ba=wvd(GydF4y2Ba___GydF4y2Ba“MinThreshold”,GydF4y2Ba脱粒GydF4y2Ba)GydF4y2Baは、振幅がGydF4y2Ba脱粒GydF4y2Baより小さいGydF4y2BaDGydF4y2Baの要素をゼロに設定します。この構文は、维格纳维尔分布と平滑化疑似 维格纳维尔分布の両方に適用されます。GydF4y2Ba

例GydF4y2Ba

[GydF4y2BaDGydF4y2Ba,GydF4y2BaFGydF4y2Ba,GydF4y2BaTGydF4y2Ba]=wvd(GydF4y2Ba___GydF4y2Ba)GydF4y2Baも周波数のベクトル (GydF4y2BaFGydF4y2Ba) と時間のベクトル (GydF4y2BaTGydF4y2Ba) を返し、これによりGydF4y2BaDGydF4y2Baが計算されます。GydF4y2Ba

出力引数なしでGydF4y2Bawvd(GydF4y2Ba___GydF4y2Ba)GydF4y2Baを使用すると、現在の 图形に 维格纳维尔分布または平滑化疑似 维格纳维尔分布がプロットされます。GydF4y2Ba

例GydF4y2Ba

すべて折りたたむGydF4y2Ba

πGydF4y2Ba /GydF4y2Ba 2.GydF4y2Ba の正規化周波数をもつ、1000サンプル インパルスと 1000サンプル トーンを生成します。2.つの信号の和の 维格纳维尔分布を計算します。GydF4y2Ba

x=零(1001,1);x(500)=10;y=sin(π*(0:1000)/2)';[d,f,t]=wvd(x+y);GydF4y2Ba

维格纳维尔分布をプロットします。GydF4y2Ba

图像SC(t、f、d)轴GydF4y2BaxyGydF4y2Ba色条GydF4y2Ba

图中包含一个轴。这些轴包含一个image类型的对象。GydF4y2Ba

出力引数なしでGydF4y2BawvdGydF4y2Baを呼び出して結果を再現します。GydF4y2Ba

wvd(x+y)GydF4y2Ba

图中包含一个轴。标题为Wigner Ville Distribution的轴包含一个image类型的对象。GydF4y2Ba

1千赫で 1.5秒間サンプリングされた 200赫兹の正弦波で構成される信号を生成します。GydF4y2Ba

fs=1000;t=(0:1/fs:1.5)';x=cos(2*pi*t*200);GydF4y2Ba

信号の 维格纳维尔分布を計算します。GydF4y2Ba

wvd(x,fs)GydF4y2Ba

图中包含一个轴。标题为Wigner Ville Distribution的轴包含一个image类型的对象。GydF4y2Ba

250赫兹と 450赫兹の間で周波数が正弦関数的に変化するチャープを信号に追加します。信号を MATLAB®の 时间表に変換します。维格纳维尔分布を計算します。GydF4y2Ba

x=x+vco(cos(2*pi*t),[250450],fs);xt=时刻表(秒(t),x);wvd(xt)GydF4y2Ba

图中包含一个轴。标题为Wigner Ville Distribution的轴包含一个image类型的对象。GydF4y2Ba

振幅が 0より小さい分布要素をゼロに設定します。GydF4y2Ba

wvd(xt,GydF4y2Ba“MinThreshold”GydF4y2Ba,0)GydF4y2Ba

图中包含一个轴。标题为Wigner Ville Distribution的轴包含一个image类型的对象。GydF4y2Ba

1千赫で 1.秒間サンプリングされた信号を生成します。信号の 1.つの成分は、測定中に周波数が 100赫兹から 400赫兹に 2.次的に増加するチャープです。信号のもう 1.つの成分は、同じ時間で周波数が 350赫兹から 50赫兹に線形的に減少するチャープです。GydF4y2Ba

信号を 时间表に保存します。GydF4y2Ba

fs=1000;t=0:1/fs:1;x=chirp(t,100,1400,GydF4y2Ba“二次”GydF4y2Ba)+啁啾(t,350,1,50);GydF4y2Ba

信号の 维格纳维尔分布を計算します。GydF4y2Ba

wvd(x,fs)GydF4y2Ba

图中包含一个轴。标题为Wigner Ville Distribution的轴包含一个image类型的对象。GydF4y2Ba

信号の平滑化疑似 维格纳维尔分布を計算します。501個の周波数点と 502個の時間点を指定します。GydF4y2Ba

wvd(x,fs,GydF4y2Ba“平滑假想”GydF4y2Ba,GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,501,GydF4y2Ba“NumTimePoints”GydF4y2Ba,502)GydF4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含一个image类型的对象。GydF4y2Ba

時間点の数を増やして、2.次チャープが表示されるようにします。GydF4y2Ba

wvd(x,fs,GydF4y2Ba“平滑假想”GydF4y2Ba,GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,501,GydF4y2Ba“NumTimePoints”GydF4y2Ba,522)GydF4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含一个image类型的对象。GydF4y2Ba

シャープなイメージを得るために、周波数点と時間点を増やします。GydF4y2Ba

wvd(x,fs,GydF4y2Ba“平滑假想”GydF4y2Ba,GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,1000,GydF4y2Ba“NumTimePoints”GydF4y2Ba,1502)GydF4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含一个image类型的对象。GydF4y2Ba

3千赫で 1.秒間サンプリングされた 2.成分信号を生成します。最初の成分は 2.次チャープで、測定中に周波数が 300赫兹から 1300赫兹に増加します。2.番目の成分は、正弦関数的に変化する周波数成分を持つチャープです。信号はホワイト ガウス ノイズに組み込まれます。連続するサンプル間の時間をGydF4y2Ba期间GydF4y2Baスカラーとして表します。GydF4y2Ba

fs=3000;t=0:1/fs:1-1/fs;dt=秒(t(2)-t(1));x1=啁啾(t,300,t(结束),1300,GydF4y2Ba“二次”GydF4y2Ba)x2=exp(2j*pi*100*cos(2*pi*2*t));x=x1+x2+randn(size(t))/10;GydF4y2Ba

信号の平滑化疑似 维格纳维尔を計算してプロットします。601サンプルのハミング ウィンドウを使用して時間領域で、305サンプルの箱型ウィンドウを使用して周波数領域で、それぞれ分布にウィンドウを適用します。表示用に 600個の周波数点を使用します。振幅がGydF4y2Ba -GydF4y2Ba 50GydF4y2Ba より小さい分布の成分をゼロに設定します。GydF4y2Ba

wvd(x,dt,GydF4y2Ba“平滑假想”GydF4y2Ba,哈明(601),雷克特温(305),GydF4y2Ba...GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,600,GydF4y2Ba“MinThreshold”GydF4y2Ba,-50)GydF4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含一个image类型的对象。GydF4y2Ba

4.つのガウス原子で構成された信号を生成します。各原子は、ガウスにより変調された正弦波で構成されています。正弦波には 100赫兹と 400赫兹の周波数があります。ガウスの中心は 150ミリ秒と 350ミリ秒で、GydF4y2Ba 0GydF4y2Ba .GydF4y2Ba 01GydF4y2Ba 2.GydF4y2Ba の分散があります。すべての原子に単位振幅があります。信号は、1千赫で 0.5秒間サンプリングされます。GydF4y2Ba

fs=1000;t=(0:1/fs:0.5);(f1=100;f2=400;mu1=0.15;mu2=0.35;gaussFun=@(A,x,mu,f)exp(-(x-mu)。^2/(2*0.01^2)).*sin(2*pi*f.*x)*A';s=gaussFun([1-11],t[mu1-mu2-mu2],[f1-f2]);GydF4y2Ba

信号の 维格纳维尔分布を計算して表示します。負の値を持つ可能性がある干渉項が、自動項の各ペアの中間に表示されます。GydF4y2Ba

西九龙文娱艺术区(s,fs)GydF4y2Ba

图中包含一个轴。标题为Wigner Ville Distribution的轴包含一个image类型的对象。GydF4y2Ba

信号の平滑化疑似 维格纳维尔分布を計算して表示します。時間と周波数の平滑化により干渉項が減衰されます。GydF4y2Ba

西九龙文娱艺术区(s,fs,GydF4y2Ba“平滑假想”GydF4y2Ba)GydF4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含一个image类型的对象。GydF4y2Ba

入力引数GydF4y2Ba

すべて折りたたむGydF4y2Ba

入力信号。ベクトルまたは単一のベクトル変数を含む MATLABGydF4y2Ba®GydF4y2Baの 时间表として指定します。GydF4y2Ba

入力信号の長さが奇数の場合は、関数によりゼロが追加され、長さが偶数になります。GydF4y2Ba

例:GydF4y2Bacos(pi/8*(0:159))'+randn(160,1)/10GydF4y2Baは、ホワイト ノイズに含まれる正弦波を指定します。GydF4y2Ba

例:GydF4y2Ba时间表(秒(0:5),兰特(6,1))GydF4y2Baは 1赫兹で 5.秒間サンプリングされた確率変数を指定します。GydF4y2Ba

データ型:GydF4y2Ba仅有一个的GydF4y2Ba|GydF4y2Ba双重的GydF4y2Ba
複素数のサポート:GydF4y2BaありGydF4y2Ba

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

サンプル時間。GydF4y2Ba期间GydF4y2Baスカラーで指定します。GydF4y2Ba

平滑化に使用される時間ウィンドウと周波数ウィンドウ。長さが奇数のベクトルとして指定します。既定では、GydF4y2BawvdGydF4y2Baは形状係数GydF4y2Baβ = 20GydF4y2Baのカイザー ウィンドウを使用します。GydF4y2Ba

  • 成双的GydF4y2Baの既定の長さは、GydF4y2Ba圆形的GydF4y2Ba(GydF4y2Ba长GydF4y2Ba(GydF4y2BaxGydF4y2Ba)/10)GydF4y2Ba以上の最も小さい奇数の整数です。GydF4y2Ba

  • 弗温GydF4y2Baの既定の長さは、GydF4y2Ba核因子GydF4y2Ba/4GydF4y2Ba以上の最も小さい奇数の整数です。GydF4y2Ba核因子GydF4y2BaはGydF4y2BaNumFrequencyPointsGydF4y2Baを使用して指定します。GydF4y2Ba

各ウィンドウの長さは、GydF4y2Ba2*GydF4y2Ba(用熟石膏、木板等)装天花板GydF4y2Ba(GydF4y2Ba长GydF4y2Ba(GydF4y2BaxGydF4y2Ba)/2)GydF4y2Ba以下でなければなりません。GydF4y2Ba

例:GydF4y2Ba凯撒GydF4y2Ba(65,0.5)GydF4y2Baは、0.5の形状係数を持つ 65サンプルのカイザー ウィンドウを指定します。GydF4y2Ba

最小非ゼロ値。実数スカラーとして指定します。関数は、振幅がGydF4y2Ba脱粒GydF4y2Baより小さいGydF4y2BaDGydF4y2Baの要素をゼロに設定します。GydF4y2Ba

名前と値のペアの引数GydF4y2Ba

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

例:GydF4y2Ba“NumFrequencyPoints”,201,“NumTimePoints”,300GydF4y2Baは、201個の周波数点と 300個の時間点における 维格纳维尔分布を計算します。GydF4y2Ba

周波数点の数。GydF4y2Ba“NumFrequencyPoints”GydF4y2Baと整数で構成されるコンマ区切りのペアとして指定します。この引数は、周波数のオーバーサンプリングの程度を制御します。周波数点の数は、GydF4y2Ba(GydF4y2Ba长GydF4y2Ba(GydF4y2Ba弗温GydF4y2Ba)+1)/2GydF4y2Ba以上でなければならず、既定よりも大きくすることはできません。GydF4y2Ba

時間点の数。GydF4y2Ba“NumTimePoints”GydF4y2Baと偶数の整数で構成されるコンマ区切りのペアとして指定します。この引数は、時間のオーバーサンプリングの程度を制御しますGydF4y2Ba[3]GydF4y2Ba。時間点の数は、GydF4y2Ba2*GydF4y2Ba长GydF4y2Ba(GydF4y2Ba成双的GydF4y2Ba)GydF4y2Ba以上でなければならず、既定よりも大きくすることはできません。GydF4y2Ba

ヒントGydF4y2Ba

入力信号が大きい場合、メモリ要求を下げて計算速度を上げるために時間点の数を減らします。GydF4y2Ba

出力引数GydF4y2Ba

すべて折りたたむGydF4y2Ba

维格纳维尔分布。行列として返されます。時間はGydF4y2BaDGydF4y2Baの列方向に、周波数は行方向に下がって増加します。行列のサイズはGydF4y2BaNGydF4y2BaFGydF4y2Ba×NGydF4y2BaTGydF4y2Baです。ここで、GydF4y2BaNGydF4y2BaFGydF4y2BaはGydF4y2BaFGydF4y2Baの長さ、GydF4y2BaNGydF4y2BaTGydF4y2BaはGydF4y2BaTGydF4y2Baの長さです。GydF4y2Ba

ベクトルとして返される周波数。GydF4y2Ba

  • 入力に時間情報がある場合、GydF4y2BaFGydF4y2Baは 赫兹単位で表される周波数を含みます。GydF4y2Ba

  • 入力に時間情報がない場合、GydF4y2BaFGydF4y2Baはラジアン/サンプル単位で表される正規化周波数を含みます。GydF4y2Ba

時点。ベクトルとして返されます。GydF4y2Ba

  • 入力に時間情報がある場合、GydF4y2BaTGydF4y2Baは秒単位で表される時間値を含みます。GydF4y2Ba

  • 入力に時間情報がない場合、GydF4y2BaTGydF4y2Baはサンプル数を含みます。GydF4y2Ba

詳細GydF4y2Ba

すべて折りたたむGydF4y2Ba

维格纳维尔分布GydF4y2Ba

“维格纳维尔分布"GydF4y2Baは、信号の高分解能な時間-周波数表現を提供します。この分布は、信号の可視化、検出、および推定に応用されます。GydF4y2Ba

連続信号GydF4y2Bax(t)GydF4y2Baでは、维格纳维尔分布は次のように定義されます。GydF4y2Ba

WVDGydF4y2Ba xGydF4y2Ba (GydF4y2Ba TGydF4y2Ba ,GydF4y2Ba FGydF4y2Ba )GydF4y2Ba =GydF4y2Ba ∫GydF4y2Ba −GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba xGydF4y2Ba (GydF4y2Ba TGydF4y2Ba +GydF4y2Ba τGydF4y2Ba 2.GydF4y2Ba )GydF4y2Ba xGydF4y2Ba *GydF4y2Ba (GydF4y2Ba TGydF4y2Ba −GydF4y2Ba τGydF4y2Ba 2.GydF4y2Ba )GydF4y2Ba EGydF4y2Ba −GydF4y2Ba JGydF4y2Ba 2.GydF4y2Ba πGydF4y2Ba FGydF4y2Ba τGydF4y2Ba DGydF4y2Ba τGydF4y2Ba .GydF4y2Ba

N個のサンプルを持つ離散信号では、分布は次のようになります。GydF4y2Ba

WVDGydF4y2Ba xGydF4y2Ba (GydF4y2Ba NGydF4y2Ba ,GydF4y2Ba KGydF4y2Ba )GydF4y2Ba =GydF4y2Ba ∑GydF4y2Ba MGydF4y2Ba =GydF4y2Ba −GydF4y2Ba NGydF4y2Ba NGydF4y2Ba xGydF4y2Ba (GydF4y2Ba NGydF4y2Ba +GydF4y2Ba MGydF4y2Ba /GydF4y2Ba 2.GydF4y2Ba )GydF4y2Ba xGydF4y2Ba *GydF4y2Ba (GydF4y2Ba NGydF4y2Ba −GydF4y2Ba MGydF4y2Ba /GydF4y2Ba 2.GydF4y2Ba )GydF4y2Ba EGydF4y2Ba −GydF4y2Ba JGydF4y2Ba 2.GydF4y2Ba πGydF4y2Ba KGydF4y2Ba MGydF4y2Ba /GydF4y2Ba NGydF4y2Ba .GydF4y2Ba

奇数値の Mでは、この定義は半整数サンプル値での信号の評価を必要とします。これにより内挿が必要となるため、エイリアシングを避けるために離散フーリエ変換をゼロ パディングする必要があります。GydF4y2Ba

维格纳维尔分布には、その解釈を複雑にしがちな干渉項が含まれています。分布をシャープにするために、ローパス ウィンドウで定義をフィルター処理できます。平滑化疑似 维格纳维尔分布は、時間と周波数の平滑化に個別のウィンドウを使用します。GydF4y2Ba

SPWVDGydF4y2Ba xGydF4y2Ba GGydF4y2Ba ,GydF4y2Ba HGydF4y2Ba (GydF4y2Ba TGydF4y2Ba ,GydF4y2Ba FGydF4y2Ba )GydF4y2Ba =GydF4y2Ba ∫GydF4y2Ba −GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba GGydF4y2Ba (GydF4y2Ba TGydF4y2Ba )GydF4y2Ba HGydF4y2Ba (GydF4y2Ba FGydF4y2Ba )GydF4y2Ba xGydF4y2Ba (GydF4y2Ba TGydF4y2Ba +GydF4y2Ba τGydF4y2Ba 2.GydF4y2Ba )GydF4y2Ba xGydF4y2Ba *GydF4y2Ba (GydF4y2Ba TGydF4y2Ba −GydF4y2Ba τGydF4y2Ba 2.GydF4y2Ba )GydF4y2Ba EGydF4y2Ba −GydF4y2Ba JGydF4y2Ba 2.GydF4y2Ba πGydF4y2Ba FGydF4y2Ba τGydF4y2Ba DGydF4y2Ba τGydF4y2Ba .GydF4y2Ba

参照GydF4y2Ba

[1] 《时频分析:理论与应用》。新泽西州恩格伍德悬崖:普伦蒂斯大厅,1995年。GydF4y2Ba

[2] 《信号处理的小波之旅》,第二版。加利福尼亚州圣地亚哥:学术出版社,1999年。GydF4y2Ba

[3] 计算二次时频分布的快速和内存有效算法〉《应用和计算谐波分析》,第35卷,第2期,2013年,第350-358页。GydF4y2Ba

拡張機能GydF4y2Ba

参考GydF4y2Ba

関数GydF4y2Ba

R2018bで導入GydF4y2Ba