主要内容gydF4y2Ba

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

xwvdgydF4y2Ba

交差能量分布と交差平滑化疑似能量分布gydF4y2Ba

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

構文gydF4y2Ba

d = xwvd (x, y)gydF4y2Ba
d = xwvd (x, y, fs)gydF4y2Ba
d = xwvd (x, y, ts)gydF4y2Ba
d = xwvd (gydF4y2Ba＿＿＿gydF4y2Ba, ' smoothedPseudo ')gydF4y2Ba
d = xwvd (gydF4y2Ba＿＿＿gydF4y2Ba“smoothedPseudo”,双胞胎,fwin)gydF4y2Ba
d = xwvd (gydF4y2Ba＿＿＿gydF4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,nf)gydF4y2Ba
d = xwvd (gydF4y2Ba＿＿＿gydF4y2Ba、“MinThreshold”打)gydF4y2Ba
[d, f t] = xwvd (gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba
xwvd (gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba

説明gydF4y2Ba

例gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba，gydF4y2BaygydF4y2Ba）gydF4y2Baは,gydF4y2BaxgydF4y2BaとgydF4y2BaygydF4y2Baの交差能量分布を返します。gydF4y2Ba

例gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba，gydF4y2BaygydF4y2Ba，gydF4y2BafsgydF4y2Ba）gydF4y2Baは,gydF4y2BaxgydF4y2BaとgydF4y2BaygydF4y2BaがレートgydF4y2BafsgydF4y2Baでサンプリングされたときの交差能量分布を返します。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba，gydF4y2BaygydF4y2Ba，gydF4y2BatsgydF4y2Ba）gydF4y2Baは,gydF4y2BaxgydF4y2BaとgydF4y2BaygydF4y2Baがサンプル間の時間間隔gydF4y2BatsgydF4y2Baでサンプリングされたときの交差能量分布を返します。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba＿＿＿gydF4y2Ba, ' smoothedPseudo ')gydF4y2Baは,gydF4y2BaxgydF4y2BaとgydF4y2BaygydF4y2Baの交差平滑化疑似能量分布を返します。この関数は,入力信号の長さを使用して,時間と周波数の平滑化に使用されるウィンドウの長さを選択します。この構文には,前の構文の入力引数を任意に組み合わせて含めることができます。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba＿＿＿gydF4y2Ba“smoothedPseudo”,gydF4y2Ba双胞胎gydF4y2Ba，gydF4y2BafwingydF4y2Ba）gydF4y2Baは,平滑化に使用される時間ウィンドウ(gydF4y2Ba双胞胎gydF4y2Ba)と周波数ウィンドウ(gydF4y2BafwingydF4y2Ba)を指定します。時間または周波数の平滑化に既定のウィンドウを使用するには,対応する引数を空(gydF4y2Ba［］gydF4y2Ba)に指定します。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba＿＿＿gydF4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,gydF4y2BanfgydF4y2Ba）gydF4y2Baは,gydF4y2BanfgydF4y2Ba個の周波数点を使用して交差平滑化疑似能量分布を計算します。この構文では,gydF4y2Ba双胞胎gydF4y2BaとgydF4y2BafwingydF4y2Baを指定できます。また,省略することもできます。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba＿＿＿gydF4y2Ba“MinThreshold”,gydF4y2Ba打gydF4y2Ba）gydF4y2Baは,振幅がgydF4y2Ba打gydF4y2Baより小さいgydF4y2BadgydF4y2Baの要素をゼロに設定します。この構文は,交差能量分布と交差平滑化疑似能量分布の両方に適用されます。gydF4y2Ba

［gydF4y2BadgydF4y2Ba，gydF4y2BafgydF4y2Ba，gydF4y2BatgydF4y2Ba) = xwvd (gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Baも周波数のベクトル(gydF4y2BafgydF4y2Ba)と時間のベクトル(gydF4y2BatgydF4y2Ba)を返し,これによりgydF4y2BadgydF4y2Baが計算されます。gydF4y2Ba

例gydF4y2Ba

出力引数なしでgydF4y2Baxwvd (gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Baを使用すると,現在の图に交差能量分布または交差平滑化疑似能量分布の実数部がプロットされます。gydF4y2Ba

例gydF4y2Ba

すべて折りたたむgydF4y2Ba

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

1 kHzで1秒間サンプリングされ,ホワイトノイズに組み込まれた2つの信号を生成します。片方の信号は,周波数150 Hzの正弦波です。もう片方の信号は,200 Hzと400 Hzの間で周波数が正弦関数的に変化するチャープです。ノイズの分散はgydF4y2Ba 0gydF4y2Ba ．gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba です。gydF4y2Ba

fs = 1000;t = (0:1 / fs: 1) ';X = cos(2*pi*t*150) + 0.1*randn(size(t));y = vco (cos(3 *π* t) (200 400), fs) + 0.1 * randn(大小(t));gydF4y2Ba

信号の和の能量分布を計算します。gydF4y2Ba

项(x + y, fs)gydF4y2Ba

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

信号の交差能量分布を計算してプロットします。交差分布は,能量分布の交差項に対応します。gydF4y2Ba

fs xwvd (x, y)gydF4y2Ba

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

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

2つのチャープで構成された2チャネル信号を生成します。信号は3 kHzで1秒間サンプリングされます。最初のチャープの初期周波数は400 Hzで,サンプリングの最後には800 Hzに到達します。2番目のチャープは500 Hzから開始し,最後には1000 Hzに到達します。2番目のチャープの振幅は最初のチャープ2倍です。gydF4y2Ba

fs = 3000;t = (0:1 / fs: 1 - 1 / fs) ';x1 =唧唧喳喳(800 t, 1400 t(结束),);x2 = 2 *唧唧喳喳(1000 t, 200 t(结束),);gydF4y2Ba

信号を时间表として保存します。2つのチャネルの交差 Wigner-Ville 分布を計算してプロットします。

xt =时间表(秒(t), x1, x2);xwvd (xt (: 1), xt (:, 2))gydF4y2Ba

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

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

既知の基準信号と交差能量分布を使用して,信号の瞬時周波数を計算します。gydF4y2Ba

1 kHzで1秒間サンプリングされたガウス原子で構成される基準信号を作成します。ガウス原子は,ガウスにより変調された正弦波です。50赫兹の正弦波周波数を指定します。64年ガウスの中心はミリ秒で,gydF4y2Ba 0gydF4y2Ba ．gydF4y2Ba 01gydF4y2Ba 2gydF4y2Ba の分散があります。gydF4y2Ba

fs = 1 e3;t = (0:1 / fs: 1 - 1 / fs) ';μ= 0.064;σ= 0.01;傅氏国际= 50;(xr = exp () - tμ。^ 2 /(2 *σ^ 2))。*罪(2 *π*国际* t);gydF4y2Ba

チャープで構成される,解析用の”未知の“信号を作成します。この信号は0.4秒で突然開始し,0.5秒後に突然終了します。この間に,チャープの周波数は400 Hzから100 Hzに線形的に減少します。gydF4y2Ba

f0 = 400;f1 = 100;xa = 0(大小(t));xa (t > 0.4 & t < = 0.9) =唧唧声((0:1 / fs: 0.5 - 1 / fs)”,f0, 0.5, f1);gydF4y2Ba

未知の信号と基準信号の和で構成される2成分信号を作成します。結果の平滑化疑似能量分布は,“理想的な”時間——周波数表現を提供します。gydF4y2Ba

平滑化疑似能量分布を計算して表示します。gydF4y2Ba

w =项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba）;项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba）gydF4y2Ba

图中包含一个坐标轴。标题为平滑伪维格纳-维尔分布的轴包含一个类型为image的对象。gydF4y2Ba

未知の信号と基準信号の交差能量分布を計算します。分布の絶対値を取り,振幅が10より小さい要素をゼロに設定します。交差能量分布は,2成分信号の交差項に等しくなります。

交差能量分布の実数部をプロットします。gydF4y2Ba

[c、fc、tc] = xwvd (xa xr, fs);c = abs (c);c (c < 10) = 0;xwvd (xa xr, fs)gydF4y2Ba

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

理想的な時間——周波数表現を交差能量分布に追加することで,能量交差項を強調します。能量分布の交差項は,基準信号と未知の信号の中間に発生します。gydF4y2Ba

D = w + c;d = abs(真正的(d));显示亮度图像(tc、fc、d)轴gydF4y2BaxygydF4y2BacolorbargydF4y2Ba

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

交差項に対応する高エネルギーリッジを特定してプロットします。リッジを分離するには,交差分布に非ゼロのエネルギーがある時間値を検索します。gydF4y2Ba

ff = tfridge (c, fc);电视=和(c) > 0;ff = ff(电视);tc = tc(电视);持有gydF4y2Ba在gydF4y2Ba情节(tc、ff,gydF4y2Ba“r——”gydF4y2Ba，gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba

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

リッジと参照関数を使用して,未知の信号の瞬時周波数を再構成します。瞬時周波数を時間の関数としてプロットします。gydF4y2Ba

测试= 2*tc - mu;= 2*ff - fsin;情节(测试、节日)gydF4y2Ba

图中包含一个坐标轴。轴包含一个线型对象。gydF4y2Ba

入力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

入力信号。ベクトルまたはそれぞれが単一のベクトル変数を含む MATLAB®gydF4y2Baの时间表として指定します。gydF4y2BaxgydF4y2BaとgydF4y2BaygydF4y2Baは,両方ともベクトルであるか両方とも时间表である必要があり,同じ長さを持たなければなりません。gydF4y2Ba

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

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

例:gydF4y2Ba时间表(秒(0:5)”,兰德(6,1))gydF4y2Baは1 Hzで4秒間サンプリングされた確率変数を指定します。gydF4y2Ba

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

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

サンプル時間。gydF4y2Ba持续时间gydF4y2Baスカラーで指定します。gydF4y2Ba

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

  • 双胞胎gydF4y2Baの既定の長さは,gydF4y2Ba轮gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 10)gydF4y2Ba以上の最も小さい奇数の整数です。gydF4y2Ba

  • fwingydF4y2Baの既定の長さは,gydF4y2BanfgydF4y2Ba／4gydF4y2Ba以上の最も小さい奇数の整数です。gydF4y2Ba

各ウィンドウの長さは,gydF4y2Ba2 *gydF4y2Ba装天花板gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 2)gydF4y2Ba以下でなければなりません。gydF4y2Ba

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

周波数点の数。整数として指定します。この引数は、周波数のオーバーサンプリングの程度を制御します。周波数点の数は、（gydF4y2Ba长度gydF4y2Ba（gydF4y2BafwingydF4y2Ba) + 1) / 2gydF4y2Ba以上でなければならず,既定よりも大きくすることはできません。gydF4y2Ba

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

出力引数gydF4y2Ba

すべて折りたたむgydF4y2Ba

交差能量分布。行列として返されます。時間はgydF4y2BadgydF4y2Baの列方向に,周波数は行方向に下がって増加します。行列のサイズはgydF4y2BaNgydF4y2BafgydF4y2Ba×NgydF4y2BatgydF4y2Baです。ここで,gydF4y2BaNgydF4y2BafgydF4y2BaはgydF4y2BafgydF4y2Baの長さ,gydF4y2BaNgydF4y2BatgydF4y2BaはgydF4y2BatgydF4y2Baの長さです。gydF4y2Ba

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

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

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

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

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

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

時間点の数は,gydF4y2Ba4 *gydF4y2Ba装天花板gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 2)gydF4y2Baに固定されます。gydF4y2Ba

詳細gydF4y2Ba

すべて折りたたむgydF4y2Ba

交差能量分布gydF4y2Ba

連続信号gydF4y2Bax (t)gydF4y2BaおよびgydF4y2Bay (t)gydF4y2Baでは”,交差能量分布”gydF4y2Baは次のように定義されます。gydF4y2Ba

XWVDgydF4y2Ba xgydF4y2Ba ，gydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba ，gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∫gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba ygydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba tgydF4y2Ba −gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba τgydF4y2Ba dgydF4y2Ba τgydF4y2Ba ．gydF4y2Ba

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

XWVDgydF4y2Ba xgydF4y2Ba ，gydF4y2Ba ygydF4y2Ba （gydF4y2Ba ngydF4y2Ba ，gydF4y2Ba kgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba ＝gydF4y2Ba −gydF4y2Ba NgydF4y2Ba NgydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba +gydF4y2Ba 米gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba ygydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba kgydF4y2Ba 米gydF4y2Ba /gydF4y2Ba NgydF4y2Ba ．gydF4y2Ba

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

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

XSPWVDgydF4y2Ba xgydF4y2Ba ，gydF4y2Ba ygydF4y2Ba 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 2gydF4y2Ba ）gydF4y2Ba ygydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba tgydF4y2Ba −gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba τgydF4y2Ba dgydF4y2Ba τgydF4y2Ba ．gydF4y2Ba

参照gydF4y2Ba

[1]科恩,莱昂。时频分析:理论与应用。Englewood Cliffs, NJ: Prentice-Hall, 1995。gydF4y2Ba

[2] Mallat,史蒂芬。信号处理的小波变换。第二版。加州圣地亚哥:学术出版社，1999。gydF4y2Ba

Malnar, Damir, Victor Sucic, Boualem Boashash。“一种基于交叉项几何的方法，利用交叉维格纳-维尔分布进行组件瞬时频率估计。”第11届国际信息科学会议，信号处理及其应用(ISSPA)，第1217-1222页。蒙特利尔:IEEEgydF4y2Ba®gydF4y2Ba, 2012年。gydF4y2Ba

拡張機能gydF4y2Ba

参考gydF4y2Ba

関数gydF4y2Ba

トピックgydF4y2Ba

R2018bで導入gydF4y2Ba

小波工具箱ドキュメンテーションgydF4y2Ba

サポートgydF4y2Ba

仿真软件MAT金宝appLAB,その他の製品をお試しくださいgydF4y2Ba

無料評価版のダウンロードgydF4y2Ba