f6_5.m

来自「《Matlab时频分析技术及其应用》（人民邮电2006年出版）一书的MATALB」· M 代码 · 共 17 行

%产生线性调频信号
sig=fmlin(200,0.1,0.4);
%计算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);
%计算一阶和二阶矩
[fm,B2]=momttfr(tfr,'tfrwv');

%结果显示
figure(1);
plot(fm,'LineWidth',2);
xlabel('时间 t');
ylabel('一阶矩 fm');
figure(2);
plot(B2,'LineWidth',2);
xlabel('时间 t');
ylabel('二阶矩 B2');