このページ最新ではありませ。をクリックし,英语のの最新版版を参照参照し。
ヒルベルトでれるの,単一成分の信号の瞬时周波数です。。単一単一成分成分成分成分ののの信号ははは,単一のリッジリッジリッジリッジリッジリッジリッジリッジによっては単一正弦波チャープようなが含まれます。
1 kHzで2秒间したチャープを生成ます。数数数がががががががががががががが100hzで100秒后200hzに秒后秒后秒后秒后秒后秒后秒后秒后
FS = 1000;t = 0:1/fs:2-1/fs;y = chirp(t,100,1,200);
关数pspectrum
に実装た时间フーリエ変换ををて,のスペクトスペクトログラムをを推定推定しししますます信号信号は,,各时点各时点各时点ののの単一単一の
pspectrum(y,fs,“频谱图”)
瞬时周波测定ために,解析计算位相ををしします。。导关数をスケーリングスケーリングすることことことでで
z =希尔伯特(y);instfrq = fs/(2*pi)*diff(unwrap(angle(z)));clf图(t(2:end),instfrq)ylim([0 fs/2])
关数Instfreq
1ステップステップ瞬时数を计算て表示し。。
Instfreq(y,fs,'方法',,,,“希尔伯特”)
信号がではない场合,メソッドがし。。。
1023 Hzで260 Hzと60Hzと90Hzの22つのつの正弦波の和をしますますます。。。スペクトスペクトをを计算计算しててててプロットますますます各时点2
FS = 1023;t = 0:1/fs:2-1/fs;x = sin(2*pi*60*t)+sin(2*pi*90*t);pspectrum(x,fs,“频谱图”)yticks([60 90])
解析を,その位相を微分し。正弦波周波をを囲む领域领域でででズームズームインししますます。解析解析信号により,,波周波数波周波数のの
z =希尔伯特(x);instfrq = fs/(2*pi)*diff(unwrap(angle(z)));图(t(2:end),instfrq)ylim([60 90])xlabel(“时间)”)ylabel(“频率(Hz)”)
关数Instfreq
は,の推定も。。
Instfreq(x,fs,'方法',,,,“希尔伯特”)
时间の両方の数を推定に,,,频谱图
を使用て,パワーととTfridge
をを探しつのつのを追迹します。Tfridge
で,数変更に対してペナルティペナルティををををにに设定し。。
[s,f,tt] = pspectrum(x,fs,“频谱图”);numcomp = 2;[冰箱,〜,lr] = tfridge(s,f,0.1,'numridges',numcomp);pspectrum(x,fs,“频谱图”) 抓住在Plot3(TT,冰箱,ABS(S(LR)),'行宽',4)保持离开Yticks([60 90])