wigner

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