wigner
'時間-周波数' ウィグナースペクトラム
呼び出し手順
[tab]=wigner(x,h,deltat,zp)
引数
- tab
ウィグナースペクトラム (時間変数と同じ行数)
- x
解析後の信号
- h
データウインドウ
- deltat
解析用の時間ステップ (単位:サンプル)
- zp
FFTの長さ.
%pi/zpが周波数増分となります.
説明
信号の '時間-周波数' ウィグナースペクトラムを計算する関数です.
例
例
a=[488^2 488 1;408^2 408 1;568^2 568 1]; b=[1.28;0;0]; x=a\b; t=408:568; p=x'*[t.*t;t;ones(t)]; // 単位ステップ関数 u=[0*ones(408:487) ones(488:568)]; // 有限時間サイン波 s=p.*sin(2*%pi/16*t+u*%pi); // 解析する信号 s=[0*ones(0:407) s 0*ones(569:951)]; // 64点矩形ウインドウ h=ones(1,64); // ウィグナースペクトラム w=wigner(s,h,12,128); scf(); plot3d(1:69,1:64,abs(w(1:69,1:64)));
| Report an issue | ||
| << wiener | Filters | window >> |