%產(chǎn)生非平穩(wěn)信號
sig=amgauss(128).*fmlin(128);
%計算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);
%計算邊緣特性
[margt,margf,E]=margtfr(tfr);
%結(jié)果顯示
figure(1);
plot(t,margt,'LineWidth',2);
xlabel('頻率 f');
ylabel('時間邊緣 m(t)');
figure(2);
plot(f,margf,'LineWidth',2);
xlabel('頻率 f');
ylabel('頻率邊緣 m(f)');
E