複素チャープの瞬時周波数

この例では,フーリエシンクロスクイーズド変換を使用して,信号の瞬時周波数を計算する方法を示します。

正弦関数的に変化する周波数成分をもつチャープを生成します。信号はホワイトガウスノイズに組み込まれ,3千赫で1秒間サンプリングされます。

fs = 3000;t = 0:1 / fs: 1 - 1 / f;x = exp (2 j *π* 100 * cos(2 *π* 2 * t)) + randn(大小(t)) / 100;

信号のフーリエシンクロスクイーズド変換を計算し,プロットします。時間がx軸に,周波数がy軸に表示されます。

fsst (x, fs,“桠溪”）

信号の瞬時周波数を求めるために,フーリエシンクロスクイーズド変換の最大エネルギーの時間——周波数リッジを抽出します。

[海温,f, tfs] = fsst (x, fs);冰箱= tfridge (sst);

リッジを変換プロットに重ね合わせます。時間をミリ秒に,周波数をkHzに変換します。

持有在情节(t * 1000冰箱/ 1000,“r”)举行从

実信号の場合,関数instfreqを使用して瞬時周波数をより簡単に検出できます。たとえば,解析信号を計算し,その位相を微分することで複素チャープの実数部の瞬時周波数を表示します。

ax =真正的(x);instfreq (ax, fs,“方法”，希尔伯特的）

参考