alfa=0.3;bita=0.85;wo=12;T=10;fc=2*bita*T;
fs=2*fc;Ts=1/fs;N=T*fs;
%T=20秒,wo=KHz),s=exp(j*wo*t)

nt=0:N-1;
x=(nt*Ts).*exp((-1*alfa*(nt*Ts).^2)/2+(j*bita*(nt*Ts).^2)/2+j*wo*nt*Ts);
wigner(x,fs);