f6_6.m

机械测试信号处理 附书的matlab源码。

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

%结果显示
figure(1);
plot(f,tm,'LineWidth',2);
xlabel('频率 f');
ylabel('一阶矩 tm');
figure(2);
plot(f,T2,'LineWidth',2);
xlabel('频率 f');
ylabel('二阶矩 T2');