%產(chǎn)生線性調(diào)頻信號
sig=fmlin(200,0.1,0.4);
%計算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);
%計算一階和二階矩
[fm,B2]=momttfr(tfr,'tfrwv');

%結(jié)果顯示
figure(1);
plot(fm,'LineWidth',2);
xlabel('時間 t');
ylabel('一階矩 fm');
figure(2);
plot(B2,'LineWidth',2);
xlabel('時間 t');
ylabel('二階矩 B2');