wigner

引数

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)));